Non-linear Mendelian randomization: detection of biases using negative controls with a focus on BMI, Vitamin D and LDL cholesterol.

Mendelian randomisation (MR) is an established technique in epidemiological investigation, using the principle of random allocation of genetic variants at conception to estimate the causal linear effect of an exposure on an outcome. Extensions to this technique include non-linear approaches that allow for differential effects of the exposure on the outcome depending on the level of the exposure. A widely used non-linear method is the residual approach, which estimates the causal effect within different strata of the non-genetically predicted exposure (i.e. the "residual" exposure). These "local" causal estimates are then used to make inferences about non-linear effects. Recent work has identified that this method can lead to estimates that are seriously biased, and a new method-the doubly-ranked method-has been introduced as a possibly more robust approach. In this paper, we perform negative control outcome analyses in the MR context. These are analyses with outcomes onto which the exposure should have no predicted causal effect. Using both methods we find clearly biased estimates in certain situations. We additionally examined a situation for which there are robust randomised controlled trial estimates of effects-that of low-density lipoprotein cholesterol (LDL-C) reduction onto myocardial infarction, where randomised trials have provided strong evidence of the shape of the relationship. The doubly-ranked method did not identify the same shape as the trial data, and for LDL-C and other lipids they generated some highly implausible findings. Therefore, we suggest there should be extensive simulation and empirical methodological examination of performance of both methods for NLMR under different conditions before further use of these methods. In the interim, use of NLMR methods needs justification, and a number of sanity checks (such as analysis of negative and positive control outcomes, sensitivity analyses excluding removal of strata at the extremes of the distribution, examination of biological plausibility and triangulation of results) should be performed.

The association between body mass index and metabolite response to a liquid mixed meal challenge: a Mendelian randomization study.

BACKGROUND: Metabolite abundance is a dynamic trait that varies in response to environmental stimuli and phenotypic traits, such as food consumption and body mass index (BMI, kg/m2). OBJECTIVES: In this study, we used the Netherlands Epidemiology of Obesity (NEO) study data to identify observational and causal associations between BMI and metabolite response to a liquid meal. METHODS: A liquid meal challenge was performed, and Nightingale Health metabolite profiles were collected in 5744 NEO participants. Observational and one-sample Mendelian randomization (MR) analysis were conducted to estimate the effect of BMI on metabolites (n = 229) in the fasting, postprandial, and response (or change in abundance) states. RESULTS: We observed 473 associations with BMI (175 fasting, 188 postprandial, and 110 response) in observational analyses. In MR analyses, we observed 20 metabolite traits (5 fasting, 12 postprandial, and 3 response) to be associated with BMI. MR associations included the glucogenic amino acid alanine, which was inversely associated with BMI in the response state (ß: -0.081; SE: 0.023; P = 5.91 × 10-4), suggesting that as alanine increased in postprandial abundance, that increase was attenuated with increasing BMI. CONCLUSIONS: Overall, this study showed that MR estimates were strongly correlated with observational effect estimates, suggesting that the broad associations seen between BMI and metabolite variation has a causal underpinning. Specific effects in previously unassessed postprandial and response states are detected, and these may likely mark novel life course risk exposures driven by regular nutrition.

Genome-wide characterization of circulating metabolic biomarkers.

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8-11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.

A genome-wide association study of neutrophil count in individuals associated to an African continental ancestry group facilitates studies of malaria pathogenesis.

BACKGROUND: 'Benign ethnic neutropenia' (BEN) is a heritable condition characterized by lower neutrophil counts, predominantly observed in individuals of African ancestry, and the genetic basis of BEN remains a subject of extensive research. In this study, we aimed to dissect the genetic architecture underlying neutrophil count variation through a linear-mixed model genome-wide association study (GWAS) in a population of African ancestry (N = 5976). Malaria caused by P. falciparum imposes a tremendous public health burden on people living in sub-Saharan Africa. Individuals living in malaria endemic regions often have a reduced circulating neutrophil count due to BEN, raising the possibility that reduced neutrophil counts modulate severity of malaria in susceptible populations. As a follow-up, we tested this hypothesis by conducting a Mendelian randomization (MR) analysis of neutrophil counts on severe malaria (MalariaGEN, N = 17,056). RESULTS: We carried out a GWAS of neutrophil count in individuals associated to an African continental ancestry group within UK Biobank, identifying 73 loci (r2 = 0.1) and 10 index SNPs (GCTA-COJO loci) associated with neutrophil count, including previously unknown rare loci regulating neutrophil count in a non-European population. BOLT-LMM was reliable when conducted in a non-European population, and additional covariates added to the model did not largely alter the results of the top loci or index SNPs. The two-sample bi-directional MR analysis between neutrophil count and severe malaria showed the greatest evidence for an effect between neutrophil count and severe anaemia, although the confidence intervals crossed the null. CONCLUSION: Our GWAS of neutrophil count revealed unique loci present in individuals of African ancestry. We note that a small sample-size reduced our power to identify variants with low allele frequencies and/or low effect sizes in our GWAS. Our work highlights the need for conducting large-scale biobank studies in Africa and for further exploring the link between neutrophils and severe malaria.

Impact of weight loss on cancer-related proteins in serum: results from a cluster randomised controlled trial of individuals with type 2 diabetes.

BACKGROUND: Type 2 diabetes is associated with higher risk of several cancer types. However, the biological intermediates driving this relationship are not fully understood. As novel interventions for treating and managing type 2 diabetes become increasingly available, whether they also disrupt the pathways leading to increased cancer risk is currently unknown. We investigated the effect of a type 2 diabetes intervention, in the form of intentional weight loss, on circulating proteins associated with cancer risk to gain insight into potential mechanisms linking type 2 diabetes and adiposity with cancer development. METHODS: Fasting serum samples from participants with diabetes enrolled in the Diabetes Remission Clinical Trial (DiRECT) receiving the Counterweight-Plus weight-loss programme (intervention, N = 117, mean weight-loss 10 kg, 46% diabetes remission) or best-practice care by guidelines (control, N = 143, mean weight-loss 1 kg, 4% diabetes remission) were subject to proteomic analysis using the Olink Oncology-II platform (48% of participants were female; 52% male). To identify proteins which may be altered by the weight-loss intervention, the difference in protein levels between groups at baseline and 1 year was examined using linear regression. Mendelian randomization (MR) was performed to extend these results to evaluate cancer risk and elucidate possible biological mechanisms linking type 2 diabetes and cancer development. MR analyses were conducted using independent datasets, including large cancer meta-analyses, UK Biobank, and FinnGen, to estimate potential causal relationships between proteins modified during intentional weight loss and the risk of colorectal, breast, endometrial, gallbladder, liver, and pancreatic cancers. FINDINGS: Nine proteins were modified by the intervention: glycoprotein Nmb; furin; Wnt inhibitory factor 1; toll-like receptor 3; pancreatic prohormone; erb-b2 receptor tyrosine kinase 2; hepatocyte growth factor; endothelial cell specific molecule 1 and Ret proto-oncogene (Holm corrected P-value <0.05). Mendelian randomization analyses indicated a causal relationship between predicted circulating furin and glycoprotein Nmb on breast cancer risk (odds ratio (OR) = 0.81, 95% confidence interval (CI) = 0.67-0.99, P-value = 0.03; and OR = 0.88, 95% CI = 0.78-0.99, P-value = 0.04 respectively), though these results were not supported in sensitivity analyses examining violations of MR assumptions. INTERPRETATION: Intentional weight loss among individuals with recently diagnosed diabetes may modify levels of cancer-related proteins in serum. Further evaluation of the proteins identified in this analysis could reveal molecular pathways that mediate the effect of adiposity and type 2 diabetes on cancer risk. FUNDING: The main sources of funding for this work were Diabetes UK, Cancer Research UK, World Cancer Research Fund, and Wellcome.

The metabolomic signature of weight loss and remission in the Diabetes Remission Clinical Trial (DiRECT).

AIMS/HYPOTHESIS: High-throughput metabolomics technologies in a variety of study designs have demonstrated a consistent metabolomic signature of overweight and type 2 diabetes. However, the extent to which these metabolomic patterns can be reversed with weight loss and diabetes remission has been weakly investigated. We aimed to characterise the metabolomic consequences of a weight-loss intervention in individuals with type 2 diabetes. METHODS: We analysed 574 fasted serum samples collected within an existing RCT (the Diabetes Remission Clinical Trial [DiRECT]) (N=298). In the trial, participating primary care practices were randomly assigned (1:1) to provide either a weight management programme (intervention) or best-practice care by guidelines (control) treatment to individuals with type 2 diabetes. Here, metabolomics analysis was performed on samples collected at baseline and 12 months using both untargeted MS and targeted 1H-NMR spectroscopy. Multivariable regression models were fitted to evaluate the effect of the intervention on metabolite levels. RESULTS: Decreases in branched-chain amino acids, sugars and LDL triglycerides, and increases in sphingolipids, plasmalogens and metabolites related to fatty acid metabolism were associated with the intervention (Holm-corrected p<0.05). In individuals who lost more than 9 kg between baseline and 12 months, those who achieved diabetes remission saw greater reductions in glucose, fructose and mannose, compared with those who did not achieve remission. CONCLUSIONS/INTERPRETATION: We have characterised the metabolomic effects of an integrated weight management programme previously shown to deliver weight loss and diabetes remission. A large proportion of the metabolome appears to be modifiable. Patterns of change were largely and strikingly opposite to perturbances previously documented with the development of type 2 diabetes. DATA AVAILABILITY: The data used for analysis are available on a research data repository ( https://researchdata.gla.ac.uk/ ) with access given to researchers subject to appropriate data sharing agreements. Metabolite data preparation, data pre-processing, statistical analyses and figure generation were performed in R Studio v.1.0.143 using R v.4.0.2. The R code for this study has been made publicly available on GitHub at: https://github.com/lauracorbin/metabolomics_of_direct .

Using trials of caloric restriction and bariatric surgery to explore the effects of body mass index on the circulating proteome.

Thousands of proteins circulate in the bloodstream; identifying those which associate with weight and intervention-induced weight loss may help explain mechanisms of diseases associated with adiposity. We aimed to identify consistent protein signatures of weight loss across independent studies capturing changes in body mass index (BMI). We analysed proteomic data from studies implementing caloric restriction (Diabetes Remission Clinical trial) and bariatric surgery (By-Band-Sleeve), using SomaLogic and Olink Explore1536 technologies, respectively. Linear mixed models were used to estimate the effect of the interventions on circulating proteins. Twenty-three proteins were altered in a consistent direction after both bariatric surgery and caloric restriction, suggesting that these proteins are modulated by weight change, independent of intervention type. We also integrated Mendelian randomisation (MR) estimates of the effect of BMI on proteins measured by SomaLogic from a UK blood donor cohort as a third line of causal evidence. These MR estimates provided further corroborative evidence for a role of BMI in regulating the levels of six proteins including alcohol dehydrogenase-4, nogo receptor and interleukin-1 receptor antagonist protein. These results indicate the importance of triangulation in interrogating causal relationships; further study into the role of proteins modulated by weight in disease is now warranted.

Application of Mendelian randomization to explore the causal role of the human gut microbiome in colorectal cancer.

The role of the human gut microbiome in colorectal cancer (CRC) is unclear as most studies on the topic are unable to discern correlation from causation. We apply two-sample Mendelian randomization (MR) to estimate the causal relationship between the gut microbiome and CRC. We used summary-level data from independent genome-wide association studies to estimate the causal effect of 14 microbial traits (n = 3890 individuals) on overall CRC (55,168 cases, 65,160 controls) and site-specific CRC risk, conducting several sensitivity analyses to understand the nature of results. Initial MR analysis suggested that a higher abundance of Bifidobacterium and presence of an unclassified group of bacteria within the Bacteroidales order in the gut increased overall and site-specific CRC risk. However, sensitivity analyses suggested that instruments used to estimate relationships were likely complex and involved in many potential horizontal pleiotropic pathways, demonstrating that caution is needed when interpreting MR analyses with gut microbiome exposures. In assessing reverse causality, we did not find strong evidence that CRC causally affected these microbial traits. Whilst our study initially identified potential causal roles for two microbial traits in CRC, importantly, further exploration of these relationships highlighted that these were unlikely to reflect causality.

Loci for insulin processing and secretion provide insight into type 2 diabetes risk.

Insulin secretion is critical for glucose homeostasis, and increased levels of the precursor proinsulin relative to insulin indicate pancreatic islet beta-cell stress and insufficient insulin secretory capacity in the setting of insulin resistance. We conducted meta-analyses of genome-wide association results for fasting proinsulin from 16 European-ancestry studies in 45,861 individuals. We found 36 independent signals at 30 loci (p value < 5 × 10-8), which validated 12 previously reported loci for proinsulin and ten additional loci previously identified for another glycemic trait. Half of the alleles associated with higher proinsulin showed higher rather than lower effects on glucose levels, corresponding to different mechanisms. Proinsulin loci included genes that affect prohormone convertases, beta-cell dysfunction, vesicle trafficking, beta-cell transcriptional regulation, and lysosomes/autophagy processes. We colocalized 11 proinsulin signals with islet expression quantitative trait locus (eQTL) data, suggesting candidate genes, including ARSG, WIPI1, SLC7A14, and SIX3. The NKX6-3/ANK1 proinsulin signal colocalized with a T2D signal and an adipose ANK1 eQTL signal but not the islet NKX6-3 eQTL. Signals were enriched for islet enhancers, and we showed a plausible islet regulatory mechanism for the lead signal in the MADD locus. These results show how detailed genetic studies of an intermediate phenotype can elucidate mechanisms that may predispose one to disease.

A multivariant recall-by-genotype study of the metabolomic signature of BMI.

OBJECTIVE: This study estimated the effect of BMI on circulating metabolites in young adults using a recall-by-genotype study design. METHODS: A recall-by-genotype study was implemented in the Avon Longitudinal Study of Parents and Children. Samples from 756 participants were selected for untargeted metabolomics analysis based on low versus high genetic liability for higher BMI defined by a genetic risk score (GRS). Regression analyses were performed to investigate associations between BMI GRS group and relative abundance of 973 metabolites. RESULTS: After correction for multiple testing, 29 metabolites were associated with BMI GRS group. Bilirubin was among the most strongly associated metabolites, with reduced levels measured in individuals in the high-BMI GRS group (ß = -0.32, 95% CI: -0.46 to -0.18, Benjamini-Hochberg adjusted p = 0.005). This study observed associations between BMI GRS group and the levels of several potentially diet-related metabolites, including hippurate, which had lower mean abundance in individuals in the high-BMI GRS group (ß = -0.29, 95% CI: -0.44 to -0.15, Benjamini-Hochberg adjusted p = 0.008). CONCLUSIONS: Together with existing literature, these results suggest that a genetic predisposition to higher BMI captures differences in metabolism leading to adiposity gain. In the absence of prospective data, separating these effects from the downstream consequences of weight gain is challenging.

metaboprep: an R package for preanalysis data description and processing.

MOTIVATION: Metabolomics is an increasingly common part of health research and there is need for preanalytical data processing. Researchers typically need to characterize the data and to exclude errors within the context of the intended analysis. Whilst some preprocessing steps are common, there is currently a lack of standardization and reporting transparency for these procedures. RESULTS: Here, we introduce metaboprep, a standardized data processing workflow to extract and characterize high quality metabolomics datasets. The package extracts data from preformed worksheets, provides summary statistics and enables the user to select samples and metabolites for their analysis based on a set of quality metrics. A report summarizing quality metrics and the influence of available batch variables on the data are generated for the purpose of open disclosure. Where possible, we provide users flexibility in defining their own selection thresholds. AVAILABILITY AND IMPLEMENTATION: metaboprep is an open-source R package available at https://github.com/MRCIEU/metaboprep. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A framework for research into continental ancestry groups of the UK Biobank.

BACKGROUND: The UK Biobank is a large prospective cohort, based in the UK, that has deep phenotypic and genomic data on roughly a half a million individuals. Included in this resource are data on approximately 78,000 individuals with "non-white British ancestry." While most epidemiology studies have focused predominantly on populations of European ancestry, there is an opportunity to contribute to the study of health and disease for a broader segment of the population by making use of the UK Biobank's "non-white British ancestry" samples. Here, we present an empirical description of the continental ancestry and population structure among the individuals in this UK Biobank subset. RESULTS: Reference populations from the 1000 Genomes Project for Africa, Europe, East Asia, and South Asia were used to estimate ancestry for each individual. Those with at least 80% ancestry in one of these four continental ancestry groups were taken forward (N = 62,484). Principal component and K-means clustering analyses were used to identify and characterize population structure within each ancestry group. Of the approximately 78,000 individuals in the UK Biobank that are of "non-white British" ancestry, 50,685, 6653, 2782, and 2364 individuals were associated to the European, African, South Asian, and East Asian continental ancestry groups, respectively. Each continental ancestry group exhibits prominent population structure that is consistent with self-reported country of birth data and geography. CONCLUSIONS: Methods outlined here provide an avenue to leverage UK Biobank's deeply phenotyped data allowing researchers to maximize its potential in the study of health and disease in individuals of non-white British ancestry.

GWAS of stool frequency provides insights into gastrointestinal motility and irritable bowel syndrome.

Gut dysmotility is associated with constipation, diarrhea, and functional gastrointestinal disorders like irritable bowel syndrome (IBS), although its molecular underpinnings are poorly characterized. We studied stool frequency (defined by the number of bowel movements per day, based on questionnaire data) as a proxy for gut motility in a GWAS meta-analysis including 167,875 individuals from UK Biobank and four smaller population-based cohorts. We identify 14 loci associated with stool frequency (p ≤ 5.0 × 10-8). Gene set and pathway analyses detected enrichment for genes involved in neurotransmitter/neuropeptide signaling and preferentially expressed in enteric motor neurons controlling peristalsis. PheWAS identified pleiotropic associations with dysmotility syndromes and the response to their pharmacological treatment. The genetic architecture of stool frequency correlates with that of IBS, and UK Biobank participants from the top 1% of stool frequency polygenic score distribution were associated with 5× higher risk of IBS with diarrhea. These findings pave the way for the identification of actionable pathological mechanisms in IBS and the dysmotility syndromes.

Correction: Telford et al. Expanding the Geographic Characterisation of Epstein-Barr Virus Variation through Gene-Based Approaches. Microorganisms 2020, 8, 1686.

The authors wish to make the following correction to this paper [...].

Genetic Studies of Metabolomics Change After a Liquid Meal Illuminate Novel Pathways for Glucose and Lipid Metabolism.

Humans spend the greater part of the day in a postprandial state. However, the genetic basis of postprandial blood measures is relatively uncharted territory. We examined the genetics of variation in concentrations of postprandial metabolites (t = 150 min) in response to a liquid mixed meal through genome-wide association studies (GWAS) performed in the Netherlands Epidemiology of Obesity (NEO) study (n = 5,705). The metabolite response GWAS identified an association between glucose change and rs10830963:G in the melatonin receptor 1B (ß [SE] -0.23 [0.03], P = 2.15 × 10-19). In addition, the ANKRD55 locus led by rs458741:C showed strong associations with extremely large VLDL (XXLVLDL) particle response (XXLVLDL total cholesterol: ß [SE] 0.17 [0.03], P = 5.76 × 10-10; XXLVLDL cholesterol ester: ß [SE] 0.17 [0.03], P = 9.74 × 10-10), which also revealed strong associations with body composition and diabetes in the UK Biobank (P < 5 × 10-8). Furthermore, the associations between XXLVLDL response and insulinogenic index, HOMA-ß, Matsuda insulin sensitivity index, and HbA1c in the NEO study implied the role of chylomicron synthesis in diabetes (with false discovery rate-corrected q <0.05). To conclude, genetic studies of metabolomics change after a liquid meal illuminate novel pathways for glucose and lipid metabolism. Further studies are warranted to corroborate biological pathways of the ANKRD55 locus underlying diabetes.

Effects of adiposity on the human plasma proteome: observational and Mendelian randomisation estimates.

BACKGROUND: Variation in adiposity is associated with cardiometabolic disease outcomes, but mechanisms leading from this exposure to disease are unclear. This study aimed to estimate effects of body mass index (BMI) on an extensive set of circulating proteins. METHODS: We used SomaLogic proteomic data from up to 2737 healthy participants from the INTERVAL study. Associations between self-reported BMI and 3622 unique plasma proteins were explored using linear regression. These were complemented by Mendelian randomisation (MR) analyses using a genetic risk score (GRS) comprised of 654 BMI-associated polymorphisms from a recent genome-wide association study (GWAS) of adult BMI. A disease enrichment analysis was performed using DAVID Bioinformatics 6.8 for proteins which were altered by BMI. RESULTS: Observationally, BMI was associated with 1576 proteins (P < 1.4 × 10-5), with particularly strong evidence for a positive association with leptin and fatty acid-binding protein-4 (FABP4), and a negative association with sex hormone-binding globulin (SHBG). Observational estimates were likely confounded, but the GRS for BMI did not associate with measured confounders. MR analyses provided evidence for a causal relationship between BMI and eight proteins including leptin (0.63 standard deviation (SD) per SD BMI, 95% CI 0.48-0.79, P = 1.6 × 10-15), FABP4 (0.64 SD per SD BMI, 95% CI 0.46-0.83, P = 6.7 × 10-12) and SHBG (-0.45 SD per SD BMI, 95% CI -0.65 to -0.25, P = 1.4 × 10-5). There was agreement in the magnitude of observational and MR estimates (R2 = 0.33) and evidence that proteins most strongly altered by BMI were enriched for genes involved in cardiovascular disease. CONCLUSIONS: This study provides evidence for a broad impact of adiposity on the human proteome. Proteins strongly altered by BMI include those involved in regulating appetite, sex hormones and inflammation; such proteins are also enriched for cardiovascular disease-related genes. Altogether, results help focus attention onto new proteomic signatures of obesity-related disease.

Loss-of-function mutations in the melanocortin 4 receptor in a UK birth cohort.

Mutations in the melanocortin 4 receptor gene (MC4R) are associated with obesity but little is known about the prevalence and impact of such mutations throughout human growth and development. We examined the MC4R coding sequence in 5,724 participants from the Avon Longitudinal Study of Parents and Children, functionally characterized all nonsynonymous MC4R variants and examined their association with anthropometric phenotypes from childhood to early adulthood. The frequency of heterozygous loss-of-function (LoF) mutations in MC4R was ~1 in 337 (0.30%), considerably higher than previous estimates. At age 18 years, mean differences in body weight, body mass index and fat mass between carriers and noncarriers of LoF mutations were 17.76 kg (95% CI 9.41, 26.10), 4.84 kg m-2 (95% CI 2.19, 7.49) and 14.78 kg (95% CI 8.56, 20.99), respectively. MC4R LoF mutations may be more common than previously reported and carriers of such variants may enter adult life with a substantial burden of excess adiposity.

Dietary intake of vitamin A, lung function and incident asthma in childhood.

BACKGROUND: Longitudinal epidemiological data are scarce on the relationship between dietary intake of vitamin A and respiratory outcomes in childhood. We investigated whether a higher intake of preformed vitamin A or pro-vitamin ß-carotene in mid-childhood is associated with higher lung function and with asthma risk in adolescence. METHODS: In the Avon Longitudinal Study of Parents and Children, dietary intakes of preformed vitamin A and ß-carotene equivalents were estimated by food frequency questionnaire at 7 years of age. Post-bronchodilator forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and forced expiratory flow at 25-75% of FVC (FEF25-75%) were measured at 15.5 years and transformed to z-scores. Incident asthma was defined by new cases of doctor-diagnosed asthma at age 11 or 14 years. RESULTS: In multivariable adjusted models, a higher intake of preformed vitamin A was associated with higher lung function and a lower risk of incident asthma: comparing top versus bottom quartiles of intake, regression coefficients for FEV1 and FEF25-75% were 0.21 (95% CI 0.05-0.38; ptrend=0.008) and 0.18 (95% CI 0.03-0.32; ptrend=0.02), respectively; odds ratios for FEV1/FVC below the lower limit of normal and incident asthma were 0.49 (95% CI 0.27-0.90; ptrend=0.04) and 0.68 (95% CI 0.47-0.99; ptrend=0.07), respectively. In contrast, there was no evidence for association with ß-carotene. We also found some evidence for modification of the associations between preformed vitamin A intake and lung function by BCMO1, NCOR2 and SCGB1A1 gene polymorphisms. CONCLUSION: A higher intake of preformed vitamin A, but not ß-carotene, in mid-childhood is associated with higher subsequent lung function and lower risk of fixed airflow limitation and incident asthma.

Variation of microRNA expression in the human placenta driven by population identity and sex of the newborn.

BACKGROUND: Analysis of lymphocyte cell lines revealed substantial differences in the expression of mRNA and microRNA (miRNA) among human populations. The extent of such population-associated differences in actual human tissues remains largely unexplored. The placenta is one of the few solid human tissues that can be collected in substantial numbers in a controlled manner, enabling quantitative analysis of transient biomolecules such as RNA transcripts. Here, we analyzed microRNA (miRNA) expression in human placental samples derived from 36 individuals representing four genetically distinct human populations: African Americans, European Americans, South Asians, and East Asians. All samples were collected at the same hospital following a unified protocol, thus minimizing potential biases that might influence the results. RESULTS: Sequence analysis of the miRNA fraction yielded 938 annotated and 70 novel miRNA transcripts expressed in the placenta. Of them, 82 (9%) of annotated and 11 (16%) of novel miRNAs displayed quantitative expression differences among populations, generally reflecting reported genetic and mRNA-expression-based distances. Several co-expressed miRNA clusters stood out from the rest of the population-associated differences in terms of miRNA evolutionary age, tissue-specificity, and disease-association characteristics. Among three non-environmental influenced demographic parameters, the second largest contributor to miRNA expression variation after population was the sex of the newborn, with 32 miRNAs (3% of detected) exhibiting significant expression differences depending on whether the newborn was male or female. Male-associated miRNAs were evolutionarily younger and correlated inversely with the expression of target mRNA involved in neuron-related functions. In contrast, both male and female-associated miRNAs appeared to mediate different types of hormonal responses. Demographic factors further affected reported imprinted expression of 66 placental miRNAs: the imprinting strength correlated with the mother's weight, but not height. CONCLUSIONS: Our results showed that among 12 assessed demographic variables, population affiliation and fetal sex had a substantial influence on miRNA expression variation among human placental samples. The effect of newborn-sex-associated miRNA differences further led to expression inhibition of the target genes clustering in specific functional pathways. By contrast, population-driven miRNA differences might mainly represent neutral changes with minimal functional impacts.