National Institute of Health Dr. Ricardo Jorge
otherLisbon, Portugal
Research output, citation impact, and the most-cited recent papers from National Institute of Health Dr. Ricardo Jorge (Portugal). Aggregated across the NobleBlocks index of 300M+ scholarly works.
Top-cited papers from National Institute of Health Dr. Ricardo Jorge
Simulated gastro-intestinal digestion is widely employed in many fields of food and nutritional sciences, as conducting human trials are often costly, resource intensive, and ethically disputable. As a consequence, in vitro alternatives that determine endpoints such as the bioaccessibility of nutrients and non-nutrients or the digestibility of macronutrients (e.g. lipids, proteins and carbohydrates) are used for screening and building new hypotheses. Various digestion models have been proposed, often impeding the possibility to compare results across research teams. For example, a large variety of enzymes from different sources such as of porcine, rabbit or human origin have been used, differing in their activity and characterization. Differences in pH, mineral type, ionic strength and digestion time, which alter enzyme activity and other phenomena, may also considerably alter results. Other parameters such as the presence of phospholipids, individual enzymes such as gastric lipase and digestive emulsifiers vs. their mixtures (e.g. pancreatin and bile salts), and the ratio of food bolus to digestive fluids, have also been discussed at length. In the present consensus paper, within the COST Infogest network, we propose a general standardised and practical static digestion method based on physiologically relevant conditions that can be applied for various endpoints, which may be amended to accommodate further specific requirements. A frameset of parameters including the oral, gastric and small intestinal digestion are outlined and their relevance discussed in relation to available in vivo data and enzymes. This consensus paper will give a detailed protocol and a line-by-line, guidance, recommendations and justifications but also limitation of the proposed model. This harmonised static, in vitro digestion method for food should aid the production of more comparable data in the future.
Developing a mechanistic understanding of the impact of food structure and composition on human health has increasingly involved simulating digestion in the upper gastrointestinal tract. These simulations have used a wide range of different conditions that often have very little physiological relevance, and this impedes the meaningful comparison of results. The standardized protocol presented here is based on an international consensus developed by the COST INFOGEST network. The method is designed to be used with standard laboratory equipment and requires limited experience to encourage a wide range of researchers to adopt it. It is a static digestion method that uses constant ratios of meal to digestive fluids and a constant pH for each step of digestion. This makes the method simple to use but not suitable for simulating digestion kinetics. Using this method, food samples are subjected to sequential oral, gastric and intestinal digestion while parameters such as electrolytes, enzymes, bile, dilution, pH and time of digestion are based on available physiological data. This amended and improved digestion method (INFOGEST 2.0) avoids challenges associated with the original method, such as the inclusion of the oral phase and the use of gastric lipase. The method can be used to assess the endpoints resulting from digestion of foods by analyzing the digestion products (e.g., peptides/amino acids, fatty acids, simple sugars) and evaluating the release of micronutrients from the food matrix. The whole protocol can be completed in ~7 d, including ~5 d required for the determination of enzyme activities.
Disorders of the brain can exhibit considerable epidemiological comorbidity and often share symptoms, provoking debate about their etiologic overlap. We quantified the genetic sharing of 25 brain disorders from genome-wide association studies of 265,218 patients and 784,643 control participants and assessed their relationship to 17 phenotypes from 1,191,588 individuals. Psychiatric disorders share common variant risk, whereas neurological disorders appear more distinct from one another and from the psychiatric disorders. We also identified significant sharing between disorders and a number of brain phenotypes, including cognitive measures. Further, we conducted simulations to explore how statistical power, diagnostic misclassification, and phenotypic heterogeneity affect genetic correlations. These results highlight the importance of common genetic variation as a risk factor for brain disorders and the value of heritability-based methods in understanding their etiology.
BACKGROUND: The systematic evaluation of the results of time-series studies of air pollution is challenged by differences in model specification and publication bias. METHODS: ) with daily all-cause, cardiovascular, and respiratory mortality across multiple countries or regions. Daily data on mortality and air pollution were collected from 652 cities in 24 countries or regions. We used overdispersed generalized additive models with random-effects meta-analysis to investigate the associations. Two-pollutant models were fitted to test the robustness of the associations. Concentration-response curves from each city were pooled to allow global estimates to be derived. RESULTS: concentration were 0.68% (95% CI, 0.59 to 0.77), 0.55% (95% CI, 0.45 to 0.66), and 0.74% (95% CI, 0.53 to 0.95). These associations remained significant after adjustment for gaseous pollutants. Associations were stronger in locations with lower annual mean PM concentrations and higher annual mean temperatures. The pooled concentration-response curves showed a consistent increase in daily mortality with increasing PM concentration, with steeper slopes at lower PM concentrations. CONCLUSIONS: and daily all-cause, cardiovascular, and respiratory mortality in more than 600 cities across the globe. These data reinforce the evidence of a link between mortality and PM concentration established in regional and local studies. (Funded by the National Natural Science Foundation of China and others.).
Climate change affects human health; however, there have been no large-scale, systematic efforts to quantify the heat-related human health impacts that have already occurred due to climate change. Here, we use empirical data from 732 locations in 43 countries to estimate the mortality burdens associated with the additional heat exposure that has resulted from recent human-induced warming, during the period 1991–2018. Across all study countries, we find that 37.0% (range 20.5–76.3%) of warm-season heat-related deaths can be attributed to anthropogenic climate change and that increased mortality is evident on every continent. Burdens varied geographically but were of the order of dozens to hundreds of deaths per year in many locations. Our findings support the urgent need for more ambitious mitigation and adaptation strategies to minimize the public health impacts of climate change. Current and future climate change is expected to impact human health, both indirectly and directly, through increasing temperatures. Climate change has already had an impact and is responsible for 37% of warm-season heat-related deaths between 1991 and 2018, with increases in mortality observed globally.
Many factors are involved in determining the latitudinal and altitudinal spread of the important tick vector Ixodes ricinus (Acari: Ixodidae) in Europe, as well as in changes in the distribution within its prior endemic zones. This paper builds on published literature and unpublished expert opinion from the VBORNET network with the aim of reviewing the evidence for these changes in Europe and discusses the many climatic, ecological, landscape and anthropogenic drivers. These can be divided into those directly related to climatic change, contributing to an expansion in the tick's geographic range at extremes of altitude in central Europe, and at extremes of latitude in Scandinavia; those related to changes in the distribution of tick hosts, particularly roe deer and other cervids; other ecological changes such as habitat connectivity and changes in land management; and finally, anthropogenically induced changes. These factors are strongly interlinked and often not well quantified. Although a change in climate plays an important role in certain geographic regions, for much of Europe it is non-climatic factors that are becoming increasingly important. How we manage habitats on a landscape scale, and the changes in the distribution and abundance of tick hosts are important considerations during our assessment and management of the public health risks associated with ticks and tick-borne disease issues in 21(st) century Europe. Better understanding and mapping of the spread of I. ricinus (and changes in its abundance) is, however, essential to assess the risk of the spread of infections transmitted by this vector species. Enhanced tick surveillance with harmonized approaches for comparison of data enabling the follow-up of trends at EU level will improve the messages on risk related to tick-borne diseases to policy makers, other stake holders and to the general public.
Abstract The Tiled ECMWF Scheme for Surface Exchanges over Land (TESSEL) is used operationally in the Integrated Forecast System (IFS) for describing the evolution of soil, vegetation, and snow over the continents at diverse spatial resolutions. A revised land surface hydrology (H-TESSEL) is introduced in the ECMWF operational model to address shortcomings of the land surface scheme, specifically the lack of surface runoff and the choice of a global uniform soil texture. New infiltration and runoff schemes are introduced with a dependency on the soil texture and standard deviation of orography. A set of experiments in stand-alone mode is used to assess the improved prediction of soil moisture at the local scale against field site observations. Comparison with basin-scale water balance (BSWB) and Global Runoff Data Centre (GRDC) datasets indicates a consistently larger dynamical range of land water mass over large continental areas and an improved prediction of river runoff, while the effect on atmospheric fluxes is fairly small. Finally, the ECMWF data assimilation and prediction systems are used to verify the effect on surface and near-surface quantities in the atmospheric-coupled mode. A midlatitude error reduction is seen both in soil moisture and in 2-m temperature.
BACKGROUND: Exposure to cold or hot temperatures is associated with premature deaths. We aimed to evaluate the global, regional, and national mortality burden associated with non-optimal ambient temperatures. METHODS: In this modelling study, we collected time-series data on mortality and ambient temperatures from 750 locations in 43 countries and five meta-predictors at a grid size of 0·5° × 0·5° across the globe. A three-stage analysis strategy was used. First, the temperature-mortality association was fitted for each location by use of a time-series regression. Second, a multivariate meta-regression model was built between location-specific estimates and meta-predictors. Finally, the grid-specific temperature-mortality association between 2000 and 2019 was predicted by use of the fitted meta-regression and the grid-specific meta-predictors. Excess deaths due to non-optimal temperatures, the ratio between annual excess deaths and all deaths of a year (the excess death ratio), and the death rate per 100 000 residents were then calculated for each grid across the world. Grids were divided according to regional groupings of the UN Statistics Division. FINDINGS: Globally, 5 083 173 deaths (95% empirical CI [eCI] 4 087 967-5 965 520) were associated with non-optimal temperatures per year, accounting for 9·43% (95% eCI 7·58-11·07) of all deaths (8·52% [6·19-10·47] were cold-related and 0·91% [0·56-1·36] were heat-related). There were 74 temperature-related excess deaths per 100 000 residents (95% eCI 60-87). The mortality burden varied geographically. Of all excess deaths, 2 617 322 (51·49%) occurred in Asia. Eastern Europe had the highest heat-related excess death rate and Sub-Saharan Africa had the highest cold-related excess death rate. From 2000-03 to 2016-19, the global cold-related excess death ratio changed by -0·51 percentage points (95% eCI -0·61 to -0·42) and the global heat-related excess death ratio increased by 0·21 percentage points (0·13-0·31), leading to a net reduction in the overall ratio. The largest decline in overall excess death ratio occurred in South-eastern Asia, whereas excess death ratio fluctuated in Southern Asia and Europe. INTERPRETATION: Non-optimal temperatures are associated with a substantial mortality burden, which varies spatiotemporally. Our findings will benefit international, national, and local communities in developing preparedness and prevention strategies to reduce weather-related impacts immediately and under climate change scenarios. FUNDING: Australian Research Council and the Australian National Health and Medical Research Council.
). In this study, shotgun metagenomics allowed the rapid reconstruction and phylogenomic characterization of the first MPXV outbreak genome sequences, showing that this MPXV belongs to clade 3 and that the outbreak most likely has a single origin. Although 2022 MPXV (lineage B.1) clustered with 2018-2019 cases linked to an endemic country, it segregates in a divergent phylogenetic branch, likely reflecting continuous accelerated evolution. An in-depth mutational analysis suggests the action of host APOBEC3 in viral evolution as well as signs of potential MPXV human adaptation in ongoing microevolution. Our findings also indicate that genome sequencing may provide resolution to track the spread and transmission of this presumably slow-evolving double-stranded DNA virus.
BACKGROUND: Concomitant with the recent emergence of CTX-M-type extended-spectrum beta-lactamases (ESBLs), Escherichia coli has become the enterobacterial species most affected by ESBLs. Multiple locales are encountering CTX-M-positive E. coli, including specifically CTX-M-15. To gain insights into the mechanism underlying this phenomenon, we assessed clonality and diversity of virulence profiles within an international collection of CTX-M-15-positive E. coli. METHODS: Forty-one ESBL-positive E. coli isolates from eight countries and three continents (Europe, Asia and North America) were selected for study based on suspected clonality. Phylogenetic group, ERIC2 PCR profile, O H serotype, AmpC variant and antibiotic susceptibility were determined. Multilocus sequence typing (MLST) and PFGE provided additional discrimination. Virulence potential was inferred by detection of 46 virulence factor (VF) genes. RESULTS: Thirty-six (88%) of the 41 E. coli isolates exhibited the same set of core characteristics: phylogenetic group B2, ERIC2 PCR profile 1, serotype O25:H4, AmpC EC6, ciprofloxacin resistance and MLST profile ST131. By PFGE, the 36 isolates constituted one large cluster at the 68% similarity level; this comprised 17 PFGE groups (defined at 85% similarity), some of which included strains from different countries. The 36 isolates exhibited highly (91% to 100%) similar VF profiles. CONCLUSIONS: We describe a broadly disseminated, CTX-M-15-positive and virulent E. coli clonal group with highly homogeneous virulence genotypes and subgroups exhibiting highly similar PFGE profiles, suggesting recent emergence. Understanding how this clone has emerged and successfully disseminated within the hospital and community, including across national boundaries, should be a public health priority.
Abstract Body-mass index (BMI) has increased steadily in most countries in parallel with a rise in the proportion of the population who live in cities 1,2 . This has led to a widely reported view that urbanization is one of the most important drivers of the global rise in obesity 3–6 . Here we use 2,009 population-based studies, with measurements of height and weight in more than 112 million adults, to report national, regional and global trends in mean BMI segregated by place of residence (a rural or urban area) from 1985 to 2017. We show that, contrary to the dominant paradigm, more than 55% of the global rise in mean BMI from 1985 to 2017—and more than 80% in some low- and middle-income regions—was due to increases in BMI in rural areas. This large contribution stems from the fact that, with the exception of women in sub-Saharan Africa, BMI is increasing at the same rate or faster in rural areas than in cities in low- and middle-income regions. These trends have in turn resulted in a closing—and in some countries reversal—of the gap in BMI between urban and rural areas in low- and middle-income countries, especially for women. In high-income and industrialized countries, we noted a persistently higher rural BMI, especially for women. There is an urgent need for an integrated approach to rural nutrition that enhances financial and physical access to healthy foods, to avoid replacing the rural undernutrition disadvantage in poor countries with a more general malnutrition disadvantage that entails excessive consumption of low-quality calories.
OBJECTIVE: To use the rs1229984 variant in the alcohol dehydrogenase 1B gene (ADH1B) as an instrument to investigate the causal role of alcohol in cardiovascular disease. DESIGN: Mendelian randomisation meta-analysis of 56 epidemiological studies. PARTICIPANTS: 261 991 individuals of European descent, including 20 259 coronary heart disease cases and 10 164 stroke events. Data were available on ADH1B rs1229984 variant, alcohol phenotypes, and cardiovascular biomarkers. MAIN OUTCOME MEASURES: Odds ratio for coronary heart disease and stroke associated with the ADH1B variant in all individuals and by categories of alcohol consumption. RESULTS: Carriers of the A-allele of ADH1B rs1229984 consumed 17.2% fewer units of alcohol per week (95% confidence interval 15.6% to 18.9%), had a lower prevalence of binge drinking (odds ratio 0.78 (95% CI 0.73 to 0.84)), and had higher abstention (odds ratio 1.27 (1.21 to 1.34)) than non-carriers. Rs1229984 A-allele carriers had lower systolic blood pressure (-0.88 (-1.19 to -0.56) mm Hg), interleukin-6 levels (-5.2% (-7.8 to -2.4%)), waist circumference (-0.3 (-0.6 to -0.1) cm), and body mass index (-0.17 (-0.24 to -0.10) kg/m(2)). Rs1229984 A-allele carriers had lower odds of coronary heart disease (odds ratio 0.90 (0.84 to 0.96)). The protective association of the ADH1B rs1229984 A-allele variant remained the same across all categories of alcohol consumption (P=0.83 for heterogeneity). Although no association of rs1229984 was identified with the combined subtypes of stroke, carriers of the A-allele had lower odds of ischaemic stroke (odds ratio 0.83 (0.72 to 0.95)). CONCLUSIONS: Individuals with a genetic variant associated with non-drinking and lower alcohol consumption had a more favourable cardiovascular profile and a reduced risk of coronary heart disease than those without the genetic variant. This suggests that reduction of alcohol consumption, even for light to moderate drinkers, is beneficial for cardiovascular health.
Although autism spectrum disorders (ASDs) have a substantial genetic basis, most of the known genetic risk has been traced to rare variants, principally copy number variants (CNVs). To identify common risk variation, the Autism Genome Project (AGP) Consortium genotyped 1558 rigorously defined ASD families for 1 million single-nucleotide polymorphisms (SNPs) and analyzed these SNP genotypes for association with ASD. In one of four primary association analyses, the association signal for marker rs4141463, located within MACROD2, crossed the genome-wide association significance threshold of P < 5 × 10(-8). When a smaller replication sample was analyzed, the risk allele at rs4141463 was again over-transmitted; yet, consistent with the winner's curse, its effect size in the replication sample was much smaller; and, for the combined samples, the association signal barely fell below the P < 5 × 10(-8) threshold. Exploratory analyses of phenotypic subtypes yielded no significant associations after correction for multiple testing. They did, however, yield strong signals within several genes, KIAA0564, PLD5, POU6F2, ST8SIA2 and TAF1C.
Mycobacterium tuberculosis strains that are resistant to an increasing number of second-line drugs used to treat multidrug-resistant tuberculosis (MDR TB) are becoming a threat to public health worldwide. We surveyed the Network of Supranational Reference Laboratories for M. tuberculosis isolates that were resistant to second-line anti-TB drugs during 2000-2004. We defined extensively drug-resistant TB (XDR TB) as MDR TB with further resistance to > or = 3 of the 6 classes of second-line drugs. Of 23 eligible laboratories, 14 (61%) contributed data on 17,690 isolates, which reflected drug susceptibility results from 48 countries. Of 3,520 (19.9%) MDR TB isolates, 347 (9.9%) met criteria for XDR TB. Further investigation of population-based trends and expanded efforts to prevent drug resistance and effectively treat patients with MDR TB are crucial for protection of public health and control of TB.
Upstream open reading frames (uORFs) are major gene expression regulatory elements. In many eukaryotic mRNAs, one or more uORFs precede the initiation codon of the main coding region. Indeed, several studies have revealed that almost half of human transcripts present uORFs. Very interesting examples have shown that these uORFs can impact gene expression of the downstream main ORF by triggering mRNA decay or by regulating translation. Also, evidence from recent genetic and bioinformatic studies implicates disturbed uORF-mediated translational control in the etiology of many human diseases, including malignancies, metabolic or neurologic disorders, and inherited syndromes. In this review, we will briefly present the mechanisms through which uORFs regulate gene expression and how they can impact on the organism's response to different cell stress conditions. Then, we will emphasize the importance of these structures by illustrating, with specific examples, how disturbed uORF-mediated translational control can be involved in the etiology of human diseases, giving special importance to genotype-phenotype correlations. Identifying and studying more cases of uORF-altering mutations will help us to understand and establish genotype-phenotype associations, leading to advancements in diagnosis, prognosis, and treatment of many human disorders.
BACKGROUND: Various genome-wide association studies (GWAS) have been done in ischaemic stroke, identifying a few loci associated with the disease, but sample sizes have been 3500 cases or less. We established the METASTROKE collaboration with the aim of validating associations from previous GWAS and identifying novel genetic associations through meta-analysis of GWAS datasets for ischaemic stroke and its subtypes. METHODS: We meta-analysed data from 15 ischaemic stroke cohorts with a total of 12 389 individuals with ischaemic stroke and 62 004 controls, all of European ancestry. For the associations reaching genome-wide significance in METASTROKE, we did a further analysis, conditioning on the lead single nucleotide polymorphism in every associated region. Replication of novel suggestive signals was done in 13 347 cases and 29 083 controls. FINDINGS: We verified previous associations for cardioembolic stroke near PITX2 (p=2·8×10(-16)) and ZFHX3 (p=2·28×10(-8)), and for large-vessel stroke at a 9p21 locus (p=3·32×10(-5)) and HDAC9 (p=2·03×10(-12)). Additionally, we verified that all associations were subtype specific. Conditional analysis in the three regions for which the associations reached genome-wide significance (PITX2, ZFHX3, and HDAC9) indicated that all the signal in each region could be attributed to one risk haplotype. We also identified 12 potentially novel loci at p<5×10(-6). However, we were unable to replicate any of these novel associations in the replication cohort. INTERPRETATION: Our results show that, although genetic variants can be detected in patients with ischaemic stroke when compared with controls, all associations we were able to confirm are specific to a stroke subtype. This finding has two implications. First, to maximise success of genetic studies in ischaemic stroke, detailed stroke subtyping is required. Second, different genetic pathophysiological mechanisms seem to be associated with different stroke subtypes. FUNDING: Wellcome Trust, UK Medical Research Council (MRC), Australian National and Medical Health Research Council, National Institutes of Health (NIH) including National Heart, Lung and Blood Institute (NHLBI), the National Institute on Aging (NIA), the National Human Genome Research Institute (NHGRI), and the National Institute of Neurological Disorders and Stroke (NINDS).
Proteins are now widely produced in diverse microbial cell factories. The Escherichia coli is still the dominant host for recombinant protein production but, as a bacterial cell, it also has its issues: the aggregation of foreign proteins into insoluble inclusion bodies is perhaps the main limiting factor of the E. coli expression system. Conversely, E. coli benefits of cost, ease of use and scale make it essential to design new approaches directed for improved recombinant protein production in this host cell. With the aid of genetic and protein engineering novel tailored-made strategies can be designed to suit user or process requirements. Gene fusion technology has been widely used for the improvement of soluble protein production and/or purification in E. coli, and for increasing peptide's immunogenicity as well. New fusion partners are constantly emerging and complementing the traditional solutions, as for instance, the Fh8 fusion tag that has been recently studied and ranked among the best solubility enhancer partners. In this review, we provide an overview of current strategies to improve recombinant protein production in E. coli, including the key factors for successful protein production, highlighting soluble protein production, and a comprehensive summary of the latest available and traditionally used gene fusion technologies. A special emphasis is given to the recently discovered Fh8 fusion system that can be used for soluble protein production, purification, and immunogenicity in E. coli. The number of existing fusion tags will probably increase in the next few years, and efforts should be taken to better understand how fusion tags act in E. coli. This knowledge will undoubtedly drive the development of new tailored-made tools for protein production in this bacterial system.
Autism spectrum disorders (ASD) are a heterogeneous group of neurodevelopmental disorders with a complex inheritance pattern. While many rare variants in synaptic proteins have been identified in patients with ASD, little is known about their effects at the synapse and their interactions with other genetic variations. Here, following the discovery of two de novo SHANK2 deletions by the Autism Genome Project, we identified a novel 421 kb de novo SHANK2 deletion in a patient with autism. We then sequenced SHANK2 in 455 patients with ASD and 431 controls and integrated these results with those reported by Berkel et al. 2010 (n = 396 patients and n = 659 controls). We observed a significant enrichment of variants affecting conserved amino acids in 29 of 851 (3.4%) patients and in 16 of 1,090 (1.5%) controls (P = 0.004, OR = 2.37, 95% CI = 1.23-4.70). In neuronal cell cultures, the variants identified in patients were associated with a reduced synaptic density at dendrites compared to the variants only detected in controls (P = 0.0013). Interestingly, the three patients with de novo SHANK2 deletions also carried inherited CNVs at 15q11-q13 previously associated with neuropsychiatric disorders. In two cases, the nicotinic receptor CHRNA7 was duplicated and in one case the synaptic translation repressor CYFIP1 was deleted. These results strengthen the role of synaptic gene dysfunction in ASD but also highlight the presence of putative modifier genes, which is in keeping with the "multiple hit model" for ASD. A better knowledge of these genetic interactions will be necessary to understand the complex inheritance pattern of ASD.
Mycobacterium tuberculosis drug resistance (DR) challenges effective tuberculosis disease control. Current molecular tests examine limited numbers of mutations, and although whole genome sequencing approaches could fully characterise DR, data complexity has restricted their clinical application. A library (1,325 mutations) predictive of DR for 15 anti-tuberculosis drugs was compiled and validated for 11 of them using genomic-phenotypic data from 792 strains. A rapid online 'TB-Profiler' tool was developed to report DR and strain-type profiles directly from raw sequences. Using our DR mutation library, in silico diagnostic accuracy was superior to some commercial diagnostics and alternative databases. The library will facilitate sequence-based drug-susceptibility testing.
BACKGROUND: models. OBJECTIVES: NMs are available with a range of physicochemical characteristics, which allows a more systematic toxicological analysis. Therefore, the study of ultrafine particles (UFP, <100 nm in diameter) provides an opportunity to identify plausible health effects for NMs, and the study of NMs provides an opportunity to facilitate the understanding of the mechanism of toxicity of UFP. METHODS: A workshop of experts systematically analyzed the available information and identified 19 key lessons that can facilitate knowledge exchange between these discipline areas. DISCUSSION: Key lessons range from the availability of specific techniques and standard protocols for physicochemical characterization and toxicology assessment to understanding and defining dose and the molecular mechanisms of toxicity. This review identifies a number of key areas in which additional research prioritization would facilitate both research fields simultaneously. CONCLUSION: There is now an opportunity to apply knowledge from NM toxicology and use it to better inform PM health risk research and vice versa. https://doi.org/10.1289/EHP424.