Variants in the FTO and CDKAL1 loci have recessive effects on risk of obesity and type 2 diabetes, respectively

Causal relationship between affect disorders and endometrial cancer: a Mendelian randomisation study

BACKGROUND

The aim was to assess the causal relationship between depression and anxiety disorders and endometrial cancer.

METHOD

We performed two-sample Mendelian randomisation analysis using summary statistics from genome-wide association studies to assess associations of major depressive disorder, anxiety and stress-related disorders with endometrial cancer. The genome-wide association studies(GWASs) data were derived from participants of predominantly European ancestry included in the Genome-wide Association Research Collaboration. Inverse variance-weighted, MR-Egger and weighted median MR analyses were performed, together with a range of sensitivity analyses.

RESULTS

Mendelian randomisation analysis showed no statistically significant genetic responsibility effect of anxiety and stress-related disorders on any pathological type of endometrial cancer. Only the effect of major depressive disorder under the inverse variance weighting method increasing the risk of endometrial endometrial cancer (effect 0.004 p = 0.047) and the effect of major depressive disorder under the MR-Egger method decreasing endometrial cancer of all pathology types (effect -0.691 p = 0.015) were statistically significant. Other Mendelian randomisation analyses did not show a statistically significant effect.

CONCLUSION

Major depressive disorder(MDD), anxiety and stress-related disorders(ASRD) are not genetically responsible for endometrial cancer. We consider that emotional disorders may affect endometrial cancer indirectly by affecting body mass index. This study provides us with new insights to better understand the aetiology of endometrial cancer and inform prevention strategies.

Unsupervised clustering identified clinically relevant metabolic syndrome endotypes in UK and Taiwan Biobanks

Metabolic syndrome (MetS) is a collection of cardiovascular risk factors; however, the high prevalence and heterogeneity impede effective clinical management. We conducted unsupervised clustering on individuals from UK Biobank to reveal endotypes. Five MetS subgroups were identified: Cluster 1 (C1): non-descriptive, Cluster 2 (C2): hypertensive, Cluster 3 (C3): obese, Cluster 4 (C4): lipodystrophy-like, and Cluster 5 (C5): hyperglycemic. For all of the endotypes, we identified the corresponding cardiometabolic traits and their associations with clinical outcomes. Genome-wide association studies (GWASs) were conducted to identify associated genotypic traits. We then determined endotype-specific genotypic traits and constructed polygenic risk score (PRS) models specific to each endotype. GWAS of each MetS clusters revealed different genotypic traits. C1 GWAS revealed novel findings of TRIM63, MYBPC3, MYLPF, and RAPSN. Intriguingly, C1, C3, and C4 were associated with genes highly expressed in brain tissues. MetS clusters with comparable phenotypic and genotypic traits were identified in Taiwan Biobank.

Exploring the causal relationship between body mass index and keratoconus: a Mendelian randomization study

Background

Despite reports suggesting a link between obesity and keratoconus, the causal relationship is not fully understood.

Methods

We used genome-wide association study (GWAS) data from public databases for a two-sample Mendelian randomization analysis to investigate the causal link between body mass index (BMI) and keratoconus. The primary method was inverse variance weighted (IVW), complemented by different analytical techniques and sensitivity analyses to ensure result robustness. A meta-analysis was also performed to bolster the findings' reliability.

Results

Our study identified a significant causal relationship between BMI and keratoconus. Out of 20 Mendelian randomization (MR) analyses conducted, 9 showed heterogeneity or pleiotropy. Among the 11 analyses that met all three MR assumptions, 4 demonstrated a significant causal difference between BMI and keratoconus, while the remaining 7 showed a positive trend but were not statistically significant. Meta-analysis confirmed a significant causal relationship between BMI and keratoconus.

Conclusion

There is a significant causal relationship between BMI and keratoconus, suggesting that obesity may be a risk factor for keratoconus.

Genome sequence analyses identify novel risk loci for multiple system atrophy

Multiple system atrophy (MSA) is an adult-onset, sporadic synucleinopathy characterized by parkinsonism, cerebellar ataxia, and dysautonomia. The genetic architecture of MSA is poorly understood, and treatments are limited to supportive measures. Here, we performed a comprehensive analysis of whole genome sequence data from 888 European-ancestry MSA cases and 7,128 controls to systematically investigate the genetic underpinnings of this understudied neurodegenerative disease. We identified four significantly associated risk loci using a genome-wide association study approach. Transcriptome-wide association analyses prioritized USP38-DT, KCTD7, and lnc-KCTD7-2 as novel susceptibility genes for MSA within these loci, and single-nucleus RNA sequence analysis found that the associated variants acted as cis-expression quantitative trait loci for multiple genes across neuronal and glial cell types. In conclusion, this study highlights the role of genetic determinants in the pathogenesis of MSA, and the publicly available data from this study represent a valuable resource for investigating synucleinopathies.

Integrating Genetic and Transcriptomic Data to Identify Genes Underlying Obesity Risk Loci

Genome-wide association studies (GWAS) have identified numerous body mass index (BMI) loci. However, most underlying mechanisms from risk locus to BMI remain unknown. Leveraging omics data through integrative analyses could provide more comprehensive views of biological pathways on BMI. We analyzed genotype and blood gene expression data in up to 5,619 samples from the Framingham Heart Study (FHS). Using 3,992 single nucleotide polymorphisms (SNPs) at 97 BMI loci and 20,692 transcripts within 1 Mb, we performed separate association analyses of transcript with BMI and SNP with transcript (PBMI and PSNP, respectively) and then a correlated meta-analysis between the full summary data sets (PMETA). We identified transcripts that met Bonferroni-corrected significance for each omic, were more significant in the correlated meta-analysis than each omic, and were at least nominally associated with BMI in FHS data. Among 308 significant SNP-transcript-BMI associations, we identified seven genes (NT5C2, GSTM3, SNAPC3, SPNS1, TMEM245, YPEL3, and ZNF646) in five association regions. Using an independent sample of blood gene expression data, we validated results for SNAPC3 and YPEL3. We tested for generalization of these associations in hypothalamus, nucleus accumbens, and liver and observed significant (PMETA<0.05 & PMETA Collapse

Key Words Collapse

MESH Headings Collapse

Grants

• R01 HL143885 NHLBI NIH HHS
• R01 HL142302 NHLBI NIH HHS
• R01 HL151152 NHLBI NIH HHS
• HHSN268201500001I NHLBI NIH HHS
• N01HC25195 NHLBI NIH HHS
• R01 HG010297 NHGRI NIH HHS
• HHSN268201500001C NHLBI NIH HHS
• R01 HL163262 NHLBI NIH HHS
• R01 HL092577 NHLBI NIH HHS
• P20 MD000170 NIMHD NIH HHS
• 75N92019D00031 NHLBI NIH HHS
• R01 HD057194 NICHD NIH HHS
• UL1 TR000371 NCATS NIH HHS
• R01 DK122503 NIDDK NIH HHS

Collapse

Affiliation(s)

Hanfei Xu Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA
Shreyash Gupta Department of Population Health Sciences, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Ian Dinsmore Department of Genomic Health, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Abbey Kollu Department of Psychology and Neuroscience, University of North Carolina, 235 E. Cameron Avenue, Chapel Hill, NC, 27599, USA
Anne Marie Cawley Marsico Lung Institute, University of North Carolina, 125 Mason Farm Rd, Chapel Hill, NC, 27599, USA
Mohammad Y. Anwar Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina, 135 Dauer Drive, Chapel Hill, NC, 27599, USA
Hung-Hsin Chen Institute of Biomedical Sciences, Academia Sinica, No. 128, Section 2, Academia Rd., Taipei, Nangang District, 115201, Taiwan Vanderbilt Genetics Institute, Division of Genetic Medicine, Vanderbilt University Medical Center, 1161 21st Ave S, Nashville, TN, 37232, USA
Lauren E. Petty Vanderbilt Genetics Institute, Division of Genetic Medicine, Vanderbilt University Medical Center, 1161 21st Ave S, Nashville, TN, 37232, USA
Sudha Seshadri Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA Department of Neurology, School of Medicine, Boston University, 85 East Concord Street, Boston, MA, 02118, USA Glenn Biggs Institute for Alzheimer’s & Neurodegenerative Diseases, UT Health San Antonio, 8300 Floyd Curl Drive, San Antonio, TX, 78229, USA
Misa Graff Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina, 135 Dauer Drive, Chapel Hill, NC, 27599, USA
Piper Below Vanderbilt Genetics Institute, Division of Genetic Medicine, Vanderbilt University Medical Center, 1161 21st Ave S, Nashville, TN, 37232, USA
Jennifer A. Brody Department of Medicine, Cardiovascular Health Research Unit, University of Washington, 1730 Minor Ave, Seattle, WA, 98101, USA
Geetha Chittoor Department of Population Health Sciences, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Susan P. Fisher-Hoch Department of Epidemiology, School of Public Health, UT Health Houston, Regional Academic Health Center, One West University Blvd, Brownsville, TX, 78520, USA
Nancy L. Heard-Costa Framingham Heart Study, 73 Mt Wayte Ave, Framingham, MA, 01702, USA Department of Neurology, Chobanian & Avedisian School of Medicine, Boston University, 72 E Concord St, Boston, MA, 02118, USA
Daniel Levy Population Sciences Branch, National Heart, Lung, and Blood Institute of the National Institutes of Health, 6701 Rockledge Drive, Bethesda, MD, 20892, USA
Honghuang Lin Department of Medicine, University of Massachusetts Chan Medical School, 55 N Lake Ave, Worcester, MA, 01655, USA
Ruth JF. Loos Charles Bronfman Institute for Personalized Medicine at Mount Sinai, Icahn School of Medicine at Mount Sinai, 1 Gustave L. Levy Pl, New York, NY, 10029, USA Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and Medical Sciences, University of Copenhagen, Blegdamsvej 3A, 2200, Copenhagen, Denmark
Joseph B. Mccormick Department of Epidemiology, School of Public Health, UT Health Houston, Regional Academic Health Center, One West University Blvd, Brownsville, TX, 78520, USA
Jerome I. Rotter Department of Pediatrics, The Institute for Translational Genomics and Population Sciences, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, 1124 West Carson Street, Torrance, CA, 90502, USA
Tooraj Mirshahi Department of Genomic Health, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Christopher D. Still Center for Obesity and Metabolic Health, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Anita Destefano Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA Department of Neurology, School of Medicine, Boston University, 85 East Concord Street, Boston, MA, 02118, USA
L. Adrienne Cupples Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA
Karen L Mohlke Department of Genetics, School of Medicine, University of North Carolina, Chapel Hill, NC, 27599, USA
Kari E. North Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina, 135 Dauer Drive, Chapel Hill, NC, 27599, USA
Anne E. Justice Department of Population Health Sciences, Geisinger, 100 N. Academy Ave., Danville, PA, 17822, USA
Ching-Ti Liu Department of Biostatistics, School of Public Health, Boston University, 801 Massachusettes Ave, Boston, MA, 02118, USA

Collapse

Causal relationship between type 2 diabetes mellitus and aortic dissection: insights from two-sample Mendelian randomization and mediation analysis

Objective

Some evidence suggests a reduced prevalence of type 2 diabetes mellitus (T2DM) in patients with aortic dissection (AD), a catastrophic cardiovascular illness, compared to general population. However, the conclusions were inconsistent, and the causal relationship between T2DM and AD remains unclear.

Methods

In this study, we aimed to explore the causal relationship between T2DM and AD using bidirectional Mendelian randomization (MR) analysis. Mediation MR analysis was conducted to explore and quantify the possible mediation effects of 1400 metabolites in T2DM and AD.

Results

The results of 26 datasets showed no causal relationship between T2DM and AD (P>0.05). Only one dataset (ebi-a-GCST90006934) showed that T2DM was a protective factor for AD (I9-AORTDIS) (OR=0.815, 95%CI: 0.692-0.960, P=0.014), and did not show horizontal pleiotropy (P=0.808) and heterogeneity (P=0.525). Vanillic acid glycine plays a mediator in the causal relationship between T2DM and AD. The mediator effect for vanillic acid glycine levels was -0.023 (95%CI: -0.066-0.021).

Conclusion

From the perspective of MR analysis, there might not be a causal relationship between T2DM and AD, and T2DM might not be a protective factor for AD. If a causal relationship does exist between T2DM and AD, with T2DM serving as a protective factor, vanillic acid glycine may act as a mediator and enhance such a protective effect.

Inflammatory bowel disease and allergic diseases: A Mendelian randomization study

BACKGROUND

Inflammatory bowel disease (IBD) and allergic diseases possess similar genetic backgrounds and pathogenesis. Observational studies have shown a correlation, but the exact direction of cause and effect remains unclear. The aim of this Mendelian randomization (MR) study is to assess bidirectional causality between inflammatory bowel disease and allergic diseases.

METHOD

We comprehensively analyzed the causal relationship between inflammatory bowel disease (IBD), Crohn's disease (CD), ulcerative colitis (UC) and allergic disease (asthma, Hay fever, and eczema) as a whole, allergic conjunctivitis (AC), atopic dermatitis (AD), allergic asthma (AAS), and allergic rhinitis (AR) by performing a bidirectional Mendelian randomization study using summary-level data from genome-wide association studies. The analysis results mainly came from the random-effects model of inverse variance weighted (IVW-RE). In addition, multivariate Mendelian randomization (MVMR) analysis was conducted to adjust the effect of body mass index (BMI) on the instrumental variables.

RESULTS

The IVW-RE method revealed that IBD genetically increased the risk of allergic disease as a whole (OR = 1.03, 95% CI = 1.01-1.04, fdr.p = .015), AC (OR = 1.04, 95% CI = 1.01-1.06, fdr.p = .011), and AD (OR = 1.06, 95% CI = 1.02-1.09, fdr.p = .004). Subgroup analysis further confirmed that CD increased the risk of allergic disease as a whole (OR = 1.02, 95% CI = 1.00-1.03, fdr.p = .031), AC (OR = 1.03, 95% CI = 1.01-1.05, fdr.p = .012), AD (OR = 1.06, 95% CI = 1.02-1.09, fdr.p = 2E-05), AAS (OR = 1.05, 95% CI = 1.02-1.08, fdr.p = .002) and AR (OR = 1.03, 95% CI = 1.00-1.07, fdr.p = .025), UC increased the risk of AAS (OR = 1.02, 95% CI = 0.98-1.07, fdr.p = .038). MVMR results showed that after taking BMI as secondary exposure, the causal effects of IBD on AC, IBD on AD, CD on allergic disease as a whole, CD on AC, CD on AD, CD on AAS, and CD on AR were still statistically significant. No significant association was observed in the reverse MR analysis.

CONCLUSION

This Mendelian randomized study demonstrated that IBD is a risk factor for allergic diseases, which is largely attributed to its subtype CD increasing the risk of AC, AD, ASS, and AR. Further investigations are needed to explore the causal relationship between allergic diseases and IBD.

Genetic effect of metformin use on risk of cancers: evidence from Mendelian randomization analysis

BACKGROUND

Increasing number of studies reported the positive effect of metformin on the prevention and treatment of cancers. However, the genetic causal effect of metformin utilization on the risk of common cancers was not completely demonstrated.

METHODS

Two-sample Mendelian Randomization (two-sample MR) analysis was conducted to uncover the genetically predicted causal association between metformin use and 26 kinds of cancers. Besides, two-step Mendelian Randomization (two-step MR) assessment was applied to clarify the mediators which mediated the causal effect of metformin on certain cancer. We utilized five robust analytical methods, in which the inverse variance weighting (IVW) method served as the major one. Sensitivity, pleiotropy, and heterogeneity were assessed. The genetic statistics of exposure, outcomes, and mediators were downloaded from publicly available datasets, including the Open Genome-Wide Association Study (GWAS), FinnGen consortium (FinnGen), and UK Biobank (UKB).

RESULTS

Among 26 kinds of common cancers, HER-positive breast cancer was presented with a significant causal relationship with metformin use [Beta: - 4.0982; OR: 0.0166 (95% CI: 0.0008, 0.3376); P value: 0.0077], which indicated metformin could prevent people from HER-positive breast cancer. Other cancers only showed modest associations with metformin use. Potential mediators were included in two-step MR, among which total testosterone levels (mediating effect: 24.52%) displayed significant mediating roles. Leave-one-out, MR-Egger, and MR-PRESSO analyses produced consistent outcomes.

CONCLUSION

Metformin use exhibited a genetically protective effect on HER-positive breast cancer, which was partially mediated by total testosterone levels.

Type 1 diabetes, glycemic traits, and risk of dental caries: a Mendelian randomization study

Background: Regarding past epidemiological studies, there has been disagreement over whether type 1 diabetes (T1DM) is one of the risk factors for dental caries. The purpose of this study was to determine the causative links between genetic susceptibility to T1DM, glycemic traits, and the risk of dental caries using Mendelian randomization (MR) approaches. Methods: Summary-level data were collected on genome-wide association studies (GWAS) of T1DM, fasting glucose (FG), glycated hemoglobin (HbA1c), fasting insulin (FI), and dental caries. MR was performed using the inverse-variance weighting (IVW) method, and sensitivity analyses were conducted using the MR-Egger method, weighted median, weighted mode, replication cohort, and multivariable MR conditioning on potential mediators. Results: The risk of dental caries increased as a result of genetic susceptibility to T1DM [odds ratio (OR) = 1.044; 95% confidence interval (CI) = 1.015-1.074; p = 0.003], with consistent findings in the replication cohort. The relationship between T1DM and dental caries was stable when adjusted for BMI, smoking, alcohol intake, and type 2 diabetes (T2DM) in multivariable MR. However, no significant correlations between the risk of dental caries and FG, HbA1c, or FI were found. Conclusion: These results indicate that T1DM has causal involvement in the genesis of dental caries. Therefore, periodic reinforcement of oral hygiene instructions must be added to the management and early multidisciplinary intervention of T1DM patients, especially among adolescents and teenagers, who are more susceptible to T1DM.

Causality of particulate matter on cardiovascular diseases and cardiovascular biomarkers

Background

Previous observational studies have shown that the prevalence of cardiovascular diseases (CVDs) is related to particulate matter (PM). However, given the methodological limitations of conventional observational research, it is difficult to identify causality conclusively. To explore the causality of PM on CVDs and cardiovascular biomarkers, we conducted a Mendelian randomization (MR) analysis.

Method

In this study, we obtained summary-level data for CVDs and cardiovascular biomarkers including atrial fibrillation (AF), heart failure (HF), myocardial infarction (MI), ischemic stroke (IS), stroke subtypes, body mass index (BMI), lipid traits, fasting glucose, fasting insulin, and blood pressure from several large genome-wide association studies (GWASs). Then we used two-sample MR to assess the causality of PM on CVDs and cardiovascular biomarkers, 16 single nucleotide polymorphisms (SNPs) for PM2.5 and 6 SNPs for PM10 were obtained from UK Biobank participants. Inverse variance weighting (IVW) analyses under the fixed effects model were used as the main analytical method to calculate MR Estimates, followed by multiple sensitivity analyses to confirm the robustness of the results.

Results

Our study revealed increases in PM2.5 concentration were significantly related to a higher risk of MI (odds ratio (OR), 2.578; 95% confidence interval (CI), 1.611-4.127; p = 7.920 × 10-5). Suggestive evidence was found between PM10 concentration and HF (OR, 2.015; 95% CI, 1.082-3.753; p = 0.027) and IS (OR, 2.279; 95% CI,1.099-4.723; p = 0.027). There was no evidence for an effect of PM concentration on other CVDs. Furthermore, PM2.5 concentration increases were significantly associated with increases in triglyceride (TG) (OR, 1.426; 95% CI, 1.133-1.795; p = 2.469 × 10-3) and decreases in high-density lipoprotein cholesterol (HDL-C) (OR, 0.779; 95% CI, 0.615-0.986; p = 0.038). The PM10 concentration increases were also closely related to the decreases in HDL-C (OR, 0.563; 95% CI, 0.366-0.865; p = 8.756 × 10-3). We observed no causal effect of PM on other cardiovascular biomarkers.

Conclusion

At the genetic level, our study suggested the causality of PM2.5 on MI, TG, as well HDL-C, and revealed the causality of PM10 on HF, IS, and HDL-C. Our findings indicated the need for continued improvements in air pollution abatement for CVDs prevention.

Investigating causality in the association between DNA methylation and type 2 diabetes using bidirectional two-sample Mendelian randomisation

AIMS/HYPOTHESIS

Several studies have identified associations between type 2 diabetes and DNA methylation (DNAm). However, the causal role of these associations remains unclear. This study aimed to provide evidence for a causal relationship between DNAm and type 2 diabetes.

METHODS

We used bidirectional two-sample Mendelian randomisation (2SMR) to evaluate causality at 58 CpG sites previously detected in a meta-analysis of epigenome-wide association studies (meta-EWAS) of prevalent type 2 diabetes in European populations. We retrieved genetic proxies for type 2 diabetes and DNAm from the largest genome-wide association study (GWAS) available. We also used data from the Avon Longitudinal Study of Parents and Children (ALSPAC, UK) when associations of interest were not available in the larger datasets. We identified 62 independent SNPs as proxies for type 2 diabetes, and 39 methylation quantitative trait loci as proxies for 30 of the 58 type 2 diabetes-related CpGs. We applied the Bonferroni correction for multiple testing and inferred causality based on p<0.001 for the type 2 diabetes to DNAm direction and p<0.002 for the opposing DNAm to type 2 diabetes direction in the 2SMR analysis.

RESULTS

We found strong evidence of a causal effect of DNAm at cg25536676 (DHCR24) on type 2 diabetes. An increase in transformed residuals of DNAm at this site was associated with a 43% (OR 1.43, 95% CI 1.15, 1.78, p=0.001) higher risk of type 2 diabetes. We inferred a likely causal direction for the remaining CpG sites assessed. In silico analyses showed that the CpGs analysed were enriched for expression quantitative trait methylation sites (eQTMs) and for specific traits, dependent on the direction of causality predicted by the 2SMR analysis.

CONCLUSIONS/INTERPRETATION

We identified one CpG mapping to a gene related to the metabolism of lipids (DHCR24) as a novel causal biomarker for risk of type 2 diabetes. CpGs within the same gene region have previously been associated with type 2 diabetes-related traits in observational studies (BMI, waist circumference, HDL-cholesterol, insulin) and in Mendelian randomisation analyses (LDL-cholesterol). Thus, we hypothesise that our candidate CpG in DHCR24 may be a causal mediator of the association between known modifiable risk factors and type 2 diabetes. Formal causal mediation analysis should be implemented to further validate this assumption.

Causal associations between serum amino acid levels and osteoarthritis: A Mendelian randomization study

OBJECTIVE

The association between serum amino acid (AA) levels and osteoarthritis (OA) risk remains unclear.

METHOD

We performed a two-sample Mendelian randomization (MR) analysis to analyze the causal effects of. serum AA levels on the OA risk by using summary-level genome-wide association study (GWAS) data. Inverse variance weighted (IVW) and Wald ratio were used as the main analysis.We also applied MR-Egger, Weighted median and Robust Adjusted Profile Score (MR.RAPS) methods. Heterogeneity and horizontally pleiotropic outliers were checked. The causal effects of AAs on early-onset all OA were explored. We also performed multivariable MR (MVMR) and conducted the bidirectional MR.

RESULTS

The results suggested that genetically predicted alanine (Ala), tyrosine (Tyr) and isoleucine (Ile) levels were significantly associated with OA risk (e.g., association between Ala and hip/knee OA risk: OR = 0.82, 95% CI = 0.75 ∼ 0.90, p = 1.54E-05). The study yielded little evidence of associations between genetically predicted AA levels with early-onset all OA risk. When adjusting the BMI in the MVMR model, suggestive causal effects of Ala and Tyr were also identified, while the effects of Ile substantially attenuated with OA risk. No significant associations between OA and AA levels were observed after testing for bidirectionality.

CONCLUSIONS

Some AAs, such as Ala, Tyr and Ile likely affects the OA risk especially at hip or knee joints. The findings highlight the important role that serum AAs might play in the development of OA and provided new treatment approaches to OA.

Integration of non-additive genome-wide association study with a multi-tissue transcriptome analysis of growth and carcass traits in Duroc pigs

Growth and carcass traits are of economic importance in the pig production, which affect pork quality and profitability of finishing pig production. This study used whole-genome and transcriptome sequencing technologies to identify potential candidate genes affecting growth and carcass traits in Duroc pigs. The medium (50-60 k) single nucleotide polymorphism (SNP) arrays of 4 154 Duroc pigs from three populations were imputed to whole-genome sequence data, yielding 10 463 227 markers on 18 autosomes. The dominance heritabilities estimated for growth and carcass traits ranged from 0.000 ± 0.041 to 0.161 ± 0.054. Using non-additive genome-wide association study (GWAS), we identified 80 dominance quantitative trait loci for growth and carcass traits at genome-wide significance (false discovery rate < 5%), 15 of which were also detected in our additive GWAS. After fine mapping, 31 candidate genes for dominance GWAS were annotated, and 8 of them were highlighted that have been previously reported to be associated with growth and development (e.g. SNX14, RELN and ENPP2), autosomal recessive diseases (e.g. AMPH, SNX14, RELN and CACNB4) and immune response (e.g. UNC93B1 and PPM1D). By integrating the lead SNPs with RNA-seq data of 34 pig tissues from the Pig Genotype-Tissue Expression project (https://piggtex.farmgtex.org/), we found that the rs691128548, rs333063869, and rs1110730611 have significantly dominant effects for the expression of SNX14, AMPH and UNC93B1 genes in tissues related to growth and development for pig, respectively. Finally, the identified candidate genes were significantly enriched for biological processes involved in the cell and organ development, lipids catabolic process and phosphatidylinositol 3-kinase signalling (P < 0.05). These results provide new molecular markers for meat production and quality selection of pig as well as basis for deciphering the genetic mechanisms of growth and carcass traits.

GWAS-Identified Variants for Obesity Do Not Influence the Risk of Developing Multiple Myeloma: A Population-Based Study and Meta-Analysis

Multiple myeloma (MM) is an incurable disease characterized by the presence of malignant plasma cells in the bone marrow that secrete specific monoclonal immunoglobulins into the blood. Obesity has been associated with the risk of developing solid and hematological cancers, but its role as a risk factor for MM needs to be further explored. Here, we evaluated whether 32 genome-wide association study (GWAS)-identified variants for obesity were associated with the risk of MM in 4189 German subjects from the German Multiple Myeloma Group (GMMG) cohort (2121 MM cases and 2068 controls) and 1293 Spanish subjects (206 MM cases and 1087 controls). Results were then validated through meta-analysis with data from the UKBiobank (554 MM cases and 402,714 controls) and FinnGen cohorts (914 MM cases and 248,695 controls). Finally, we evaluated the correlation of these single nucleotide polymorphisms (SNPs) with cQTL data, serum inflammatory proteins, steroid hormones, and absolute numbers of blood-derived cell populations (n = 520). The meta-analysis of the four European cohorts showed no effect of obesity-related variants on the risk of developing MM. We only found a very modest association of the POC5_rs2112347G and ADCY3_rs11676272G alleles with MM risk that did not remain significant after correction for multiple testing (per-allele OR = 1.08, p = 0.0083 and per-allele OR = 1.06, p = 0.046). No correlation between these SNPs and functional data was found, which confirms that obesity-related variants do not influence MM risk.

The association between plasma chemokines and breast cancer risk and prognosis: A mendelian randomization study

Background: Despite the potential role of several chemokines in the migration of cytotoxic immune cells to prohibit breast cancer cell proliferation, a comprehensive view of chemokines and the risk and prognosis of breast cancer is scarce, and little is known about their causal associations. Methods: With a two-sample Mendelian randomization (MR) approach, genetic instruments associated with 30 plasma chemokines were created. Their genetic associations with breast cancer and its survival by molecular subtypes were extracted from the recent genome-wide association study of 133,384 breast cancer cases and 113,789 controls, with available survival information for 96,661 patients. We further tested the associations between the polygenic risk score (PRS) for chemokines and breast cancer in the UK Biobank cohort using logistic regression models, while the association with breast cancer survival was tested using Cox regression models. In addition, the association between chemokine expression in tumors and breast cancer survival was also analyzed in the TCGA cohort using Cox regression models. Results: Plasma CCL5 was causally associated with breast cancer in the MR analysis, which was significant in the luminal and HER-2 enriched subtypes and further confirmed using PRS analysis (OR = 0.94, 95% CI = 0.89-1.00). A potential causal association with breast cancer survival was only found for plasma CCL19, especially for ER-positive patients. Although not replicated in the UK Biobank, we still found an inverse association between CCL19 expression in tumors and breast cancer overall and relapse-free survival in the TCGA cohort (HR = 0.58, 95% CI = 0.35-0.95). Conclusion: We observed an inverse association between genetic predisposition to CCL5 and breast cancer, while CCL19 was associated with breast cancer survival. These associations suggested the potential of these chemokines as tools for breast cancer prevention and treatment.

Identification of shared genetic architecture between non-alcoholic fatty liver disease and type 2 diabetes: A genome-wide analysis

BACKGROUND

The incidence of complications of non-alcoholic fatty liver disease (NAFLD) and type 2 diabetes (T2D) has been increasing.

METHOD

In order to identify the shared genetic architecture of the two disease phenotypes of NAFLD and T2D, a European population-based GWAS summary and a cross-trait meta-analysis was used to identify significant shared genes for NAFLD and T2D. The enrichment of shared genes was then determined through the use of functional enrichment analysis to investigate the relationship between genes and phenotypes. Additionally, differential gene expression analysis was performed, significant differentially expressed genes in NAFLD and T2D were identified, genes that overlapped between those that were differentially expressed and cross-trait results were reported, and enrichment analysis was performed on the core genes that had been obtained in this way. Finally, the application of a bidirectional Mendelian randomization (MR) approach determined the causal link between NAFLD and T2D.

RESULT

A total of 115 genes were discovered to be shared between NAFLD and T2D in the GWAS analysis. The enrichment analysis of these genes showed that some were involved in the processes such as the decomposition and metabolism of lipids, phospholipids, and glycerophospholipids. Additionally, through the use of differential gene expression analysis, 15 core genes were confirmed to be linked to both T2D and NAFLD. They were correlated with carcinoma cells and inflammation. Furthermore, the bidirectional MR identified a positive causal relationship between NAFLD and T2D.

CONCLUSION

Our study determined the genetic structure shared between NAFLD and T2D, offering a new reference for the genetic pathogenesis and mechanism of NAFLD and T2D comorbidities.

The CRTh2 polymorphism rs533116 G > A associates with asthma severity in older females

Background

CRTh2 is G protein coupled receptor for prostaglandin D2 (PGD)₂ expressed by immune cells that drive type 2 inflammation such as CD4⁺ T cells (Th2), eosinophils and group 2 innate lymphoid cells (ILC2) as well as structural cells including smooth muscle and epithelium. CRTh2-expressing cells are increased in the blood and airways of asthmatics and severe asthma is characterized by increased activity of the PGD₂-CRTh2 pathway. The CRTh2 single nucleotide polymorphism (SNP) rs533116 G > A is associated with development of asthma and increased Th2 cell differentiation.

Objective

To examine whether CRTh2 rs533116G > A associates with asthma severity. Since severe asthma is more common in females than males, we performed a sex-stratified analysis.

Methods

Clinical data from asthmatics (n = 170) were obtained from clinic visits and chart review. Asthma severity was assessed according to ERS/ATS guidelines. Peripheral blood cells were characterized by flow cytometry and qRT-PCR. Genotyping was performed by TaqMan assay.

Results

Older females (≥45 years) homozygous for minor A allele of rs533116 were more likely to have severe asthma, lower FEV₁, a higher prescribed dose of inhaled corticosteroid and more type 2 inflammation than females carrying GA or GG genotypes. Comparing females and males with the AA genotype also revealed that women had more type 2 inflammation.

Conclusions and significance

The polymorphism CRTh2 rs533116 G > A associates with severe asthma and type 2 inflammation in older females. This study reveals a gene-sex-aging interaction influencing the effect of CRTh2 on asthma severity.

Celebrities in the heart, strangers in the pancreatic beta cell: Voltage-gated potassium channels Kv 7.1 and Kv 11.1 bridge long QT syndrome with hyperinsulinaemia as well as type 2 diabetes

Voltage‐gated potassium (K_v) channels play an important role in the repolarization of a variety of excitable tissues, including in the cardiomyocyte and the pancreatic beta cell. Recently, individuals carrying loss‐of‐function (LoF) mutations in KCNQ1, encoding K_v7.1, and KCNH2 (hERG), encoding K_v11.1, were found to exhibit post‐prandial hyperinsulinaemia and episodes of hypoglycaemia. These LoF mutations also cause the cardiac disorder long QT syndrome (LQTS), which can be aggravated by hypoglycaemia. Interestingly, patients with LQTS also have a higher burden of diabetes compared to the background population, an apparent paradox in relation to the hyperinsulinaemic phenotype, and KCNQ1 has been identified as a type 2 diabetes risk gene. This review article summarizes the involvement of delayed rectifier K⁺ channels in pancreatic beta cell function, with emphasis on K_v7.1 and K_v11.1, using the cardiomyocyte for context. The functional and clinical consequences of LoF mutations and polymorphisms in these channels on blood glucose homeostasis are explored using evidence from pre‐clinical, clinical and genome‐wide association studies, thereby evaluating the link between LQTS, hyperinsulinaemia and type 2 diabetes.

Recessive Genome-Wide Meta-analysis Illuminates Genetic Architecture of Type 2 Diabetes

Most genome-wide association studies (GWAS) of complex traits are performed using models with additive allelic effects. Hundreds of loci associated with type 2 diabetes have been identified using this approach. Additive models, however, can miss loci with recessive effects, thereby leaving potentially important genes undiscovered. We conducted the largest GWAS meta-analysis using a recessive model for type 2 diabetes. Our discovery sample included 33,139 case subjects and 279,507 control subjects from 7 European-ancestry cohorts, including the UK Biobank. We identified 51 loci associated with type 2 diabetes, including five variants undetected by prior additive analyses. Two of the five variants had minor allele frequency of <5% and were each associated with more than a doubled risk in homozygous carriers. Using two additional cohorts, FinnGen and a Danish cohort, we replicated three of the variants, including one of the low-frequency variants, rs115018790, which had an odds ratio in homozygous carriers of 2.56 (95% CI 2.05-3.19; P = 1 × 10-16) and a stronger effect in men than in women (for interaction, P = 7 × 10-7). The signal was associated with multiple diabetes-related traits, with homozygous carriers showing a 10% decrease in LDL cholesterol and a 20% increase in triglycerides; colocalization analysis linked this signal to reduced expression of the nearby PELO gene. These results demonstrate that recessive models, when compared with GWAS using the additive approach, can identify novel loci, including large-effect variants with pathophysiological consequences relevant to type 2 diabetes.

Simultaneous test and estimation of total genetic effect in eQTL integrative analysis through mixed models

Integration of expression quantitative trait loci (eQTL) into genome-wide association studies (GWASs) is a promising manner to reveal functional roles of associated single-nucleotide polymorphisms (SNPs) in complex phenotypes and has become an active research field in post-GWAS era. However, how to efficiently incorporate eQTL mapping study into GWAS for prioritization of causal genes remains elusive. We herein proposed a novel method termed as Mixed transcriptome-wide association studies (TWAS) and mediated Variance estimation (MTV) by modeling the effects of cis-SNPs of a gene as a function of eQTL. MTV formulates the integrative method and TWAS within a unified framework via mixed models and therefore includes many prior methods/tests as special cases. We further justified MTV from another two statistical perspectives of mediation analysis and two-stage Mendelian randomization. Relative to existing methods, MTV is superior for pronounced features including the processing of direct effects of cis-SNPs on phenotypes, the powerful likelihood ratio test for assessment of joint effects of cis-SNPs and genetically regulated gene expression (GReX), two useful quantities to measure relative genetic contributions of GReX and cis-SNPs to phenotypic variance, and the computationally efferent parameter expansion expectation maximum algorithm. With extensive simulations, we identified that MTV correctly controlled the type I error in joint evaluation of the total genetic effect and proved more powerful to discover true association signals across various scenarios compared to existing methods. We finally applied MTV to 41 complex traits/diseases available from three GWASs and discovered many new associated genes that had otherwise been missed by existing methods. We also revealed that a small but substantial fraction of phenotypic variation was mediated by GReX. Overall, MTV constructs a robust and realistic modeling foundation for integrative omics analysis and has the advantage of offering more attractive biological interpretations of GWAS results.

Posttranscriptional Regulation of Insulin Resistance: Implications for Metabolic Diseases

Insulin resistance defines an impairment in the biologic response to insulin action in target tissues, primarily the liver, muscle, adipose tissue, and brain. Insulin resistance affects physiology in many ways, causing hyperglycemia, hypertension, dyslipidemia, visceral adiposity, hyperinsulinemia, elevated inflammatory markers, and endothelial dysfunction, and its persistence leads to the development metabolic disease, including diabetes, obesity, cardiovascular disease, or nonalcoholic fatty liver disease (NAFLD), as well as neurological disorders such as Alzheimer’s disease. In addition to classical transcriptional factors, posttranscriptional control of gene expression exerted by microRNAs and RNA-binding proteins constitutes a new level of regulation with important implications in metabolic homeostasis. In this review, we describe miRNAs and RBPs that control key genes involved in the insulin signaling pathway and related regulatory networks, and their impact on human metabolic diseases at the molecular level, as well as their potential use for diagnosis and future therapeutics.

Genetic insights into biological mechanisms governing human ovarian ageing

Reproductive longevity is essential for fertility and influences healthy ageing in women1,2, but insights into its underlying biological mechanisms and treatments to preserve it are limited. Here we identify 290 genetic determinants of ovarian ageing, assessed using normal variation in age at natural menopause (ANM) in about 200,000 women of European ancestry. These common alleles were associated with clinical extremes of ANM; women in the top 1% of genetic susceptibility have an equivalent risk of premature ovarian insufficiency to those carrying monogenic FMR1 premutations3. The identified loci implicate a broad range of DNA damage response (DDR) processes and include loss-of-function variants in key DDR-associated genes. Integration with experimental models demonstrates that these DDR processes act across the life-course to shape the ovarian reserve and its rate of depletion. Furthermore, we demonstrate that experimental manipulation of DDR pathways highlighted by human genetics increases fertility and extends reproductive life in mice. Causal inference analyses using the identified genetic variants indicate that extending reproductive life in women improves bone health and reduces risk of type 2 diabetes, but increases the risk of hormone-sensitive cancers. These findings provide insight into the mechanisms that govern ovarian ageing, when they act, and how they might be targeted by therapeutic approaches to extend fertility and prevent disease.

Recessive/dominant model: Alternative choice in case-control-based genome-wide association studies

An additive genetic model is usually employed in case-control-based genome-wide association studies. The model usually encodes "AA", "Aa" and "aa" ("a" represents the minor allele) as three different numbers, implying the contribution of genotype "Aa" to the phenotype is different from "AA" and "aa". From the perspective of biological phenomena, the coding is reasonable since the phenotypes of lives are not "black and white". A case-control based study, however, has only two phenotypes, case and control, which means that the phenotypes are "black and white". It suggests that a recessive/dominant model may be an alternative to the additive model. In order to investigate whether the alternative is feasible, we conducted comparative experiments on several models used in those studies through chi-square test and logistic regression. Our simulation experiments demonstrate that a recessive model is better than the additive model. The area under the curve of the former has increased by 5% compared with the latter, the discrimination of identifying risk single nucleotide polymorphisms has been improved by 61%, and the precision has also reached 1.10 times that of the latter. Furthermore, the real data experiments show that the precision and area under the curve of the former are 16% and 20% higher than the latter respectively, and the area under the curve of dominant model of the former is 13% higher than the latter. The results indicate a recessive/dominant model may be an alternative to the additive model and suggest a new route for case-control-based studies.

The impact of non-additive genetic associations on age-related complex diseases

Genome-wide association studies (GWAS) are not fully comprehensive, as current strategies typically test only the additive model, exclude the X chromosome, and use only one reference panel for genotype imputation. We implement an extensive GWAS strategy, GUIDANCE, which improves genotype imputation by using multiple reference panels and includes the analysis of the X chromosome and non-additive models to test for association. We apply this methodology to 62,281 subjects across 22 age-related diseases and identify 94 genome-wide associated loci, including 26 previously unreported. Moreover, we observe that 27.7% of the 94 loci are missed if we use standard imputation strategies with a single reference panel, such as HRC, and only test the additive model. Among the new findings, we identify three novel low-frequency recessive variants with odds ratios larger than 4, which need at least a three-fold larger sample size to be detected under the additive model. This study highlights the benefits of applying innovative strategies to better uncover the genetic architecture of complex diseases.

Variable expression quantitative trait loci analysis of breast cancer risk variants

Genome wide association studies (GWAS) have identified more than 180 variants associated with breast cancer risk, however the underlying functional mechanisms and biological pathways which confer disease susceptibility remain largely unknown. As gene expression traits are under genetic regulation we hypothesise that differences in gene expression variability may identify causal breast cancer susceptibility genes. We performed variable expression quantitative trait loci (veQTL) analysis using tissue-specific expression data from the Genotype-Tissue Expression (GTEx) Common Fund Project. veQTL analysis identified 70 associations (p < 5 × 10^–8) consisting of 60 genes and 27 breast cancer risk variants, including 55 veQTL that were observed in breast tissue only. Pathway analysis of genes associated with breast-specific veQTL revealed an enrichment of four genes (CYP11B1, CYP17A1 HSD3B2 and STAR) involved in the C21-steroidal biosynthesis pathway that converts cholesterol to breast-related hormones (e.g. oestrogen). Each of these four genes were significantly more variable in individuals homozygous for rs11075995 (A/A) breast cancer risk allele located in the FTO gene, which encodes an RNA demethylase. The A/A allele was also found associated with reduced expression of FTO, suggesting an epi-transcriptomic mechanism may underlie the dysregulation of genes involved in hormonal biosynthesis leading to an increased risk of breast cancer. These findings provide evidence that genetic variants govern high levels of expression variance in breast tissue, thus building a more comprehensive insight into the underlying biology of breast cancer risk loci.

pgainsim: an R-package to assess the mode of inheritance for quantitative trait loci in GWAS

MOTIVATION

When performing genome-wide association studies conventionally the additive genetic model is used to explore whether a single nucleotide polymorphism (SNP) is associated with a quantitative trait. But for variants, which do not follow an intermediate mode of inheritance (MOI), the recessive or the dominant genetic model can have more power to detect associations and furthermore the MOI is important for downstream analyses and clinical interpretation. When multiple MOIs are modelled the question arises, which describes the true underlying MOI best.

RESULTS

We developed an R-package allowing for the first time to determine study specific critical values when one of the three models is more informative than the other ones for a quantitative trait locus. The package allows for user-friendly simulations to determine these critical values with predefined minor allele frequencies and study sizes. For application scenarios with extensive multiple testing we integrated an interpolation functionality to determine critical values already based on a moderate number of random draws.

AVAILABILITY AND IMPLEMENTATION

The R-package pgainsim is freely available for download on Github at https://github.com/genepi-freiburg/pgainsim.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

tRNA Biology in the Pathogenesis of Diabetes: Role of Genetic and Environmental Factors

The global rise in type 2 diabetes results from a combination of genetic predisposition with environmental assaults that negatively affect insulin action in peripheral tissues and impair pancreatic β-cell function and survival. Nongenetic heritability of metabolic traits may be an important contributor to the diabetes epidemic. Transfer RNAs (tRNAs) are noncoding RNA molecules that play a crucial role in protein synthesis. tRNAs also have noncanonical functions through which they control a variety of biological processes. Genetic and environmental effects on tRNAs have emerged as novel contributors to the pathogenesis of diabetes. Indeed, altered tRNA aminoacylation, modification, and fragmentation are associated with β-cell failure, obesity, and insulin resistance. Moreover, diet-induced tRNA fragments have been linked with intergenerational inheritance of metabolic traits. Here, we provide a comprehensive review of how perturbations in tRNA biology play a role in the pathogenesis of monogenic and type 2 diabetes.

The association of clinical phenotypes to known AD/FTD genetic risk loci and their inter-relationship

To elucidate how variants in genetic risk loci previously implicated in Alzheimer’s Disease (AD) and/or frontotemporal dementia (FTD) contribute to expression of disease phenotypes, a phenome-wide association study was performed in two waves. In the first wave, we explored clinical traits associated with thirteen genetic variants previously reported to be linked to disease risk using both the 23andMe and UKB cohorts. We tested 30 additional AD variants in UKB cohort only in the second wave. APOE variants defining ε2/ε3/ε4 alleles and rs646776 were identified to be significantly associated with metabolic/cardiovascular and longevity traits. APOE variants were also significantly associated with neurological traits. ABI3 variant rs28394864 was significantly associated with cardiovascular (e.g. (hypertension, ischemic heart disease, coronary atherosclerosis, angina) and immune-related trait asthma. Both APOE variants and CLU variant were significantly associated with nearsightedness. HLA- DRB1 variant was associated with diseases with immune-related traits. Additionally, variants from 10+ AD genes (BZRAP1-AS1, ADAMTS4, ADAM10, APH1B, SCIMP, ABI3, SPPL2A, ZNF232, GRN, CD2AP, and CD33) were associated with hematological measurements such as white blood cell (leukocyte) count, monocyte count, neutrophill count, platelet count, and/or mean platelet (thrombocyte) volume (an autoimmune disease biomarker). Many of these genes are expressed specifically in microglia. The associations of ABI3 variant with cardiovascular and immune-related traits are one of the novel findings from this study. Taken together, it is evidenced that at least some AD and FTD variants are associated with multiple clinical phenotypes and not just dementia. These findings were discussed in the context of causal relationship versus pleiotropy via Mendelian randomization analysis.

Genetics of Obesity in East Asians

Obesity has become a public health problem worldwide. Compared with Europe, people in Asia tend to suffer from type 2 diabetes with a lower body mass index (BMI). Genome-wide association studies (GWASs) have identified over 750 loci associated with obesity. Although the majority of GWAS results were conducted in individuals of European ancestry, a recent GWAS in individuals of Asian ancestry has made a significant contribution to the identification of obesity susceptibility loci. Indeed, owing to the multifactorial character of obesity with a strong environmental component, the revealed loci may have distinct contributions in different ancestral genetic backgrounds and in different environments as presented through diet and exercise among other factors. Uncovering novel, yet unrevealed genes in non-European ancestries may further contribute to explaining the missing heritability for BMI. In this review, we aimed to summarize recent advances in obesity genetics in individuals of Asian ancestry. We therefore compared proposed mechanisms underlying susceptibility loci for obesity associated with individuals of European and Asian ancestries and discussed whether known genetic variants might explain ethnic differences in obesity risk. We further acknowledged that GWAS implemented in individuals of Asian ancestries have not only validated the potential role of previously specified obesity susceptibility loci but also exposed novel ones, which have been missed in the initial genetic studies in individuals of European ancestries. Thus, multi-ethnic studies have a great potential not only to contribute to a better understanding of the complex etiology of human obesity but also potentially of ethnic differences in the prevalence of obesity, which may ultimately pave new avenues in more targeted and personalized obesity treatments.

Automated feature extraction from population wearable device data identified novel loci associated with sleep and circadian rhythms

Wearable devices have been increasingly used in research to provide continuous physical activity monitoring, but how to effectively extract features remains challenging for researchers. To analyze the generated actigraphy data in large-scale population studies, we developed computationally efficient methods to derive sleep and activity features through a Hidden Markov Model-based sleep/wake identification algorithm, and circadian rhythm features through a Penalized Multi-band Learning approach adapted from machine learning. Unsupervised feature extraction is useful when labeled data are unavailable, especially in large-scale population studies. We applied these two methods to the UK Biobank wearable device data and used the derived sleep and circadian features as phenotypes in genome-wide association studies. We identified 53 genetic loci with p<5×10-8 including genes known to be associated with sleep disorders and circadian rhythms as well as novel loci associated with Body Mass Index, mental diseases and neurological disorders, which suggest shared genetic factors of sleep and circadian rhythms with physical and mental health. Further cross-tissue enrichment analysis highlights the important role of the central nervous system and the shared genetic architecture with metabolism-related traits and the metabolic system. Our study demonstrates the effectiveness of our unsupervised methods for wearable device data when additional training data cannot be easily acquired, and our study further expands the application of wearable devices in population studies and genetic studies to provide novel biological insights.

Genetic architecture of Alzheimer's disease

Advances in genetic and genomic technologies over the last thirty years have greatly enhanced our knowledge concerning the genetic architecture of Alzheimer's disease (AD). Several genes including APP, PSEN1, PSEN2, and APOE have been shown to exhibit large effects on disease susceptibility, with the remaining risk loci having much smaller effects on AD risk. Notably, common genetic variants impacting AD are not randomly distributed across the genome. Instead, these variants are enriched within regulatory elements active in human myeloid cells, and to a lesser extent liver cells, implicating these cell and tissue types as critical to disease etiology. Integrative approaches are emerging as highly effective for identifying the specific target genes through which AD risk variants act and will likely yield important insights related to potential therapeutic targets in the coming years. In the future, additional consideration of sex- and ethnicity-specific contributions to risk as well as the contribution of complex gene-gene and gene-environment interactions will likely be necessary to further improve our understanding of AD genetic architecture.

mRNA Processing: An Emerging Frontier in the Regulation of Pancreatic β Cell Function

Robust endocrine cell function, particularly β cell function, is required to maintain blood glucose homeostasis. Diabetes can result from the loss or dysfunction of β cells. Despite decades of clinical and basic research, the precise regulation of β cell function and pathogenesis in diabetes remains incompletely understood. In this review, we highlight RNA processing of mRNAs as a rapidly emerging mechanism regulating β cell function and survival. RNA-binding proteins (RBPs) and RNA modifications are primed to be the next frontier to explain many of the poorly understood molecular processes that regulate β cell formation and function, and provide an exciting potential for the development of novel therapeutics. Here we outline the current understanding of β cell specific functions of several characterized RBPs, alternative splicing events, and transcriptome wide changes in RNA methylation. We also highlight several RBPs that are dysregulated in both Type 1 and Type 2 diabetes, and discuss remaining knowledge gaps in the field.

Non-additive (dominance) effects of genetic variants associated with refractive error and myopia

Genome-wide association studies (GWAS) have revealed that the genetic contribution to certain complex diseases is well-described by Fisher's infinitesimal model in which a vast number of polymorphisms each confer a small effect. Under Fisher's model, variants have additive effects both across loci and within loci. However, the latter assumption is at odds with the common observation of dominant or recessive rare alleles responsible for monogenic disorders. Here, we searched for evidence of non-additive (dominant or recessive) effects for GWAS variants known to confer susceptibility to the highly heritable quantitative trait, refractive error. Of 146 GWAS variants examined in a discovery sample of 228,423 individuals whose refractive error phenotype was inferred from their age-of-onset of spectacle wear, only 8 had even nominal evidence (p < 0.05) of non-additive effects. In a replication sample of 73,577 individuals who underwent direct assessment of refractive error, 1 of these 8 variants had robust independent evidence of non-additive effects (rs7829127 within ZMAT4, p = 4.76E-05) while a further 2 had suggestive evidence (rs35337422 in RD3L, p = 7.21E-03 and rs12193446 in LAMA2, p = 2.57E-02). Accounting for non-additive effects had minimal impact on the accuracy of a polygenic risk score for refractive error (R2 = 6.04% vs. 6.01%). Our findings demonstrate that very few GWAS variants for refractive error show evidence of a departure from an additive mode of action and that accounting for non-additive risk variants offers little scope to improve the accuracy of polygenic risk scores for myopia.

Bivariate traits association analysis using generalized estimating equations in family data

Genome wide association study (GWAS) is becoming fundamental in the arduous task of deciphering the etiology of complex diseases. The majority of the statistical models used to address the genes-disease association consider a single response variable. However, it is common for certain diseases to have correlated phenotypes such as in cardiovascular diseases. Usually, GWAS typically sample unrelated individuals from a population and the shared familial risk factors are not investigated. In this paper, we propose to apply a bivariate model using family data that associates two phenotypes with a genetic region. Using generalized estimation equations (GEE), we model two phenotypes, either discrete, continuous or a mixture of them, as a function of genetic variables and other important covariates. We incorporate the kinship relationships into the working matrix extended to a bivariate analysis. The estimation method and the joint gene-set effect in both phenotypes are developed in this work. We also evaluate the proposed methodology with a simulation study and an application to real data.

Benefits and limitations of genome-wide association studies

Genome-wide association studies (GWAS) involve testing genetic variants across the genomes of many individuals to identify genotype-phenotype associations. GWAS have revolutionized the field of complex disease genetics over the past decade, providing numerous compelling associations for human complex traits and diseases. Despite clear successes in identifying novel disease susceptibility genes and biological pathways and in translating these findings into clinical care, GWAS have not been without controversy. Prominent criticisms include concerns that GWAS will eventually implicate the entire genome in disease predisposition and that most association signals reflect variants and genes with no direct biological relevance to disease. In this Review, we comprehensively assess the benefits and limitations of GWAS in human populations and discuss the relevance of performing more GWAS.

Genome-wide analysis indicates association between heterozygote advantage and healthy aging in humans

Background

Genetic diversity is known to confer survival advantage in many species across the tree of life. Here, we hypothesize that such pattern applies to humans as well and could be a result of higher fitness in individuals with higher genomic heterozygosity.

Results

We use healthy aging as a proxy for better health and fitness, and observe greater heterozygosity in healthy-aged individuals. Specifically, we find that only common genetic variants show significantly higher excess of heterozygosity in the healthy-aged cohort. Lack of difference in heterozygosity for low-frequency variants or disease-associated variants excludes the possibility of compensation for deleterious recessive alleles as a mechanism. In addition, coding SNPs with the highest excess of heterozygosity in the healthy-aged cohort are enriched in genes involved in extracellular matrix and glycoproteins, a group of genes known to be under long-term balancing selection. We also find that individual heterozygosity rate is a significant predictor of electronic health record (EHR)-based estimates of 10-year survival probability in men but not in women, accounting for several factors including age and ethnicity.

Conclusions

Our results demonstrate that the genomic heterozygosity is associated with human healthspan, and that the relationship between higher heterozygosity and healthy aging could be explained by heterozygote advantage. Further characterization of this relationship will have important implications in aging-associated disease risk prediction.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0758-4) contains supplementary material, which is available to authorized users.

Genome-Wide Interaction and Pathway Association Studies for Body Mass Index

Objective: We investigated gene interactions (epistasis) for body mass index (BMI) in a European-American adult female cohort via genome-wide interaction analyses (GWIA) and pathway association analyses. Methods: Genome-wide pairwise interaction analyses were carried out for BMI in 493 extremely obese cases (BMI > 35 kg/m2) and 537 never-overweight controls (BMI < 25 kg/m2). To further validate the results, specific SNPs were selected based on the GWIA results for haplotype-based association studies. Pathway-based association analyses were performed using a modified Gene Set Enrichment Algorithm (GSEA) (GenGen program) to further explore BMI-related pathways using our genome wide association study (GWAS) data set, GIANT, ENGAGE, and DIAGRAM Consortia. Results: The EXOC4-1q23.1 interaction was associated with BMI, with the most significant epistasis between rs7800006 and rs10797020 (P = 2.63 × 10-11). In the pathway-based association analysis, Tob1 pathway showed the most significant association with BMI (empirical P < 0.001, FDR = 0.044, FWER = 0.040). These findings were further validated in different populations. Conclusion: Genome-wide pairwise SNP-SNP interaction and pathway analyses suggest that EXOC4 and TOB1-related pathways may contribute to the development of obesity.

Established and emerging strategies to crack the genetic code of obesity

Tremendous progress has been made in the genetic elucidation of obesity over the past two decades, driven largely by technological, methodological and organizational innovations. Current strategies for identifying obesity-predisposing loci/genes, including cytogenetics, linkage analysis, homozygosity mapping, admixture mapping, candidate gene studies, genome-wide association studies, custom genotyping arrays, whole-exome sequencing and targeted exome sequencing, have achieved differing levels of success, and the identified loci in aggregate explain only a modest fraction of the estimated heritability of obesity. This review outlines the successes and limitations of these approaches and proposes novel strategies, including the use of exceptionally large sample sizes, the study of diverse ethnic groups and deep phenotypes and the application of innovative methods and study designs, to identify the remaining obesity-predisposing genes. The use of both established and emerging strategies has the potential to crack the genetic code of obesity in the not-too-distant future. The resulting knowledge is likely to yield improvements in obesity prediction, prevention and care.

Genetic architecture of human thinness compared to severe obesity

The variation in weight within a shared environment is largely attributable to genetic factors. Whilst many genes/loci confer susceptibility to obesity, little is known about the genetic architecture of healthy thinness. Here, we characterise the heritability of thinness which we found was comparable to that of severe obesity (h2 = 28.07 vs 32.33% respectively), although with incomplete genetic overlap (r = -0.49, 95% CI [-0.17, -0.82], p = 0.003). In a genome-wide association analysis of thinness (n = 1,471) vs severe obesity (n = 1,456), we identified 10 loci previously associated with obesity, and demonstrate enrichment for established BMI-associated loci (pbinomial = 3.05x10-5). Simulation analyses showed that different association results between the extremes were likely in agreement with additive effects across the BMI distribution, suggesting different effects on thinness and obesity could be due to their different degrees of extremeness. In further analyses, we detected a novel obesity and BMI-associated locus at PKHD1 (rs2784243, obese vs. thin p = 5.99x10-6, obese vs. controls p = 2.13x10-6 pBMI = 2.3x10-13), associations at loci recently discovered with much larger sample sizes (e.g. FAM150B and PRDM6-CEP120), and novel variants driving associations at previously established signals (e.g. rs205262 at the SNRPC/C6orf106 locus and rs112446794 at the PRDM6-CEP120 locus). Our ability to replicate loci found with much larger sample sizes demonstrates the value of clinical extremes and suggest that characterisation of the genetics of thinness may provide a more nuanced understanding of the genetic architecture of body weight regulation and may inform the identification of potential anti-obesity targets.

Mendelian randomization evaluation of causal effects of fibrinogen on incident coronary heart disease

BACKGROUND

Fibrinogen is an essential hemostatic factor and cardiovascular disease risk factor. Early attempts at evaluating the causal effect of fibrinogen on coronary heart disease (CHD) and myocardial infraction (MI) using Mendelian randomization (MR) used single variant approaches, and did not take advantage of recent genome-wide association studies (GWAS) or multi-variant, pleiotropy robust MR methodologies.

METHODS AND FINDINGS

We evaluated evidence for a causal effect of fibrinogen on both CHD and MI using MR. We used both an allele score approach and pleiotropy robust MR models. The allele score was composed of 38 fibrinogen-associated variants from recent GWAS. Initial analyses using the allele score used a meta-analysis of 11 European-ancestry prospective cohorts, free of CHD and MI at baseline, to examine incidence CHD and MI. We also applied 2 sample MR methods with data from a prevalent CHD and MI GWAS. Results are given in terms of the hazard ratio (HR) or odds ratio (OR), depending on the study design, and associated 95% confidence interval (CI). In single variant analyses no causal effect of fibrinogen on CHD or MI was observed. In multi-variant analyses using incidence CHD cases and the allele score approach, the estimated causal effect (HR) of a 1 g/L higher fibrinogen concentration was 1.62 (CI = 1.12, 2.36) when using incident cases and the allele score approach. In 2 sample MR analyses that accounted for pleiotropy, the causal estimate (OR) was reduced to 1.18 (CI = 0.98, 1.42) and 1.09 (CI = 0.89, 1.33) in the 2 most precise (smallest CI) models, out of 4 models evaluated. In the 2 sample MR analyses for MI, there was only very weak evidence of a causal effect in only 1 out of 4 models.

CONCLUSIONS

A small causal effect of fibrinogen on CHD is observed using multi-variant MR approaches which account for pleiotropy, but not single variant MR approaches. Taken together, results indicate that even with large sample sizes and multi-variant approaches MR analyses still cannot exclude the null when estimating the causal effect of fibrinogen on CHD, but that any potential causal effect is likely to be much smaller than observed in epidemiological studies.

Influence of FTO rs9939609 and Mediterranean diet on body composition and weight loss: a randomized clinical trial

Background

The Mediterranean diet (MeD) plays a key role in the prevention of obesity. Among the genes involved in obesity, the Fat mass and obesity-associated gene (FTO) is one of the most known, but its interaction with MeD remained uncertain so far.

Methods

We carried out a study on a sample of 188 Italian subjects, analyzing their FTO rs9939609 alleles, and the difference in body composition between the baseline and a 4-weeks nutritional intervention. The sample was divided into two groups: the control group of 49 subjects, and the MeD group of 139 subjects.

Results

We found significant relations between MeD and both variation of total body fat (ΔTBFat) (p = 0.00) and gynoid body fat (p = 0.04). ∆TBFat (kg) demonstrated to have a significant relation with the interaction diet-gene (p = 0.04), whereas FTO was associated with the variation of total body water (p = 0.02).

Conclusions

MeD demonstrated to be a good nutritional treatment to reduce the body fat mass, whereas data about FTO remain uncertain. Confirming or rejecting the hypothesis of FTO and its influence on body tissues during nutritional treatments is fundamental to decide whether its effect has to be taken into consideration during both development of dietetic plans and patients monitoring.

Trial Registration ClinicalTrials.gov Id: NCT01890070. Registered 01 July 2013, https://clinicaltrials.gov/ct2/show/NCT01890070

Genetics of metabolic traits in Greenlanders: lessons from an isolated population

In this review, we describe the extraordinary population of Greenland, which differs from large outbred populations of Europe and Asia, both in terms of population history and living conditions. Many years in isolation, small population size and an extreme environment have shaped the genetic composition of the Greenlandic population. The unique genetic background combined with the transition from a traditional Inuit lifestyle and diet, to a more Westernized lifestyle, has led to an increase in the prevalence of metabolic conditions like obesity, where the prevalence from 1993 to 2010 has increased from 16.4% to 19.4% among men, and from 13.0% to 25.4% among women, type 2 diabetes and cardiovascular diseases. The genetic susceptibility to metabolic conditions has been explored in Greenlanders, as well as other isolated populations, taking advantage of population-genetic properties of these populations. During the last 10 years, these studies have provided examples of loci showing evidence of positive selection, due to adaption to Arctic climate and Inuit diet, including TBC1D4 and FADS/CPT1A, and have facilitated the discovery of several loci associated with metabolic phenotypes. Most recently, the c.2433-1G>A loss-of-function variant in ADCY3 associated with obesity and type 2 diabetes was described. This locus has provided novel biological insights, as it has been shown that reduced ADCY3 function causes obesity through disrupted function in primary cilia. Future studies of isolated populations will likely provide further genetic as well as biological insights.

Identification of novel high-impact recessively inherited type 2 diabetes risk variants in the Greenlandic population

AIMS/HYPOTHESIS

In a recent study using a standard additive genetic model, we identified a TBC1D4 loss-of-function variant with a large recessive impact on risk of type 2 diabetes in Greenlanders. The aim of the current study was to identify additional genetic variation underlying type 2 diabetes using a recessive genetic model, thereby increasing the power to detect variants with recessive effects.

METHODS

We investigated three cohorts of Greenlanders (B99, n = 1401; IHIT, n = 3115; and BBH, n = 547), which were genotyped using Illumina MetaboChip. Of the 4674 genotyped individuals passing quality control, 4648 had phenotype data available, and type 2 diabetes association analyses were performed for 317 individuals with type 2 diabetes and 2631 participants with normal glucose tolerance. Statistical association analyses were performed using a linear mixed model.

RESULTS

Using a recessive genetic model, we identified two novel loci associated with type 2 diabetes in Greenlanders, namely rs870992 in ITGA1 on chromosome 5 (OR 2.79, p = 1.8 × 10-8), and rs16993330 upstream of LARGE1 on chromosome 22 (OR 3.52, p = 1.3 × 10-7). The LARGE1 variant did not reach the conventional threshold for genome-wide significance (p < 5 × 10-8) but did withstand a study-wide Bonferroni-corrected significance threshold. Both variants were common in Greenlanders, with minor allele frequencies of 23% and 16%, respectively, and were estimated to have large recessive effects on risk of type 2 diabetes in Greenlanders, compared with additively inherited variants previously observed in European populations.

CONCLUSIONS/INTERPRETATION

We demonstrate the value of using a recessive genetic model in a historically small and isolated population to identify genetic risk variants. Our findings give new insights into the genetic architecture of type 2 diabetes, and further support the existence of high-effect genetic risk factors of potential clinical relevance, particularly in isolated populations.

DATA AVAILABILITY

The Greenlandic MetaboChip-genotype data are available at European Genome-Phenome Archive (EGA; https://ega-archive.org/ ) under the accession EGAS00001002641.

The MR-Base platform supports systematic causal inference across the human phenome

Results from genome-wide association studies (GWAS) can be used to infer causal relationships between phenotypes, using a strategy known as 2-sample Mendelian randomization (2SMR) and bypassing the need for individual-level data. However, 2SMR methods are evolving rapidly and GWAS results are often insufficiently curated, undermining efficient implementation of the approach. We therefore developed MR-Base (<ext-link ext-link-type="uri" xlink:href="http://www.mrbase.org">http://www.mrbase.org</ext-link>): a platform that integrates a curated database of complete GWAS results (no restrictions according to statistical significance) with an application programming interface, web app and R packages that automate 2SMR. The software includes several sensitivity analyses for assessing the impact of horizontal pleiotropy and other violations of assumptions. The database currently comprises 11 billion single nucleotide polymorphism-trait associations from 1673 GWAS and is updated on a regular basis. Integrating data with software ensures more rigorous application of hypothesis-driven analyses and allows millions of potential causal relationships to be efficiently evaluated in phenome-wide association studies.

The MR-Base platform supports systematic causal inference across the human phenome

Results from genome-wide association studies (GWAS) can be used to infer causal relationships between phenotypes, using a strategy known as 2-sample Mendelian randomization (2SMR) and bypassing the need for individual-level data. However, 2SMR methods are evolving rapidly and GWAS results are often insufficiently curated, undermining efficient implementation of the approach. We therefore developed MR-Base (http://www.mrbase.org): a platform that integrates a curated database of complete GWAS results (no restrictions according to statistical significance) with an application programming interface, web app and R packages that automate 2SMR. The software includes several sensitivity analyses for assessing the impact of horizontal pleiotropy and other violations of assumptions. The database currently comprises 11 billion single nucleotide polymorphism-trait associations from 1673 GWAS and is updated on a regular basis. Integrating data with software ensures more rigorous application of hypothesis-driven analyses and allows millions of potential causal relationships to be efficiently evaluated in phenome-wide association studies.

Fine-mapping of 98 obesity loci in Mexican children

BACKGROUND/OBJECTIVES

Mexico has one of the highest prevalence of childhood obesity in the world. Genome-wide association studies (GWAS) for obesity have identified multiple single-nucleotide polymorphisms (SNPs) in populations of European, East Asian, and African descent. The contribution of these loci to obesity in Mexican children is unclear. We assessed the transferability of 98 obesity loci in Mexican children and fine-mapped the association signals.

SUBJECTS/METHODS

The study included 405 and 390 Mexican children with normal weight and obesity. Participants were genotyped with a genome-wide dense SNP array designed for Latino populations, allowing for the analysis of GWAS index SNPs as well as fine-mapping SNPs, totaling 750 SNPs covering 98 loci. Two genetic risk scores (GRS) were constructed: a "discovery GRS" and a "best-associated GRS", representing the number of effect alleles at the GWAS index SNPs and at the best-associated SNPs after fine-mapping for each subject.

RESULTS

Seventeen obesity loci were significantly associated with obesity, and five had fine-mapping SNPs significantly better associated with obesity than their corresponding GWAS index SNPs in Mexican children. Six obesity-associated SNPs significantly departed from additive to dominant (N = 5) or recessive (N = 1) models, and a significant interaction was found between rs274609 (TNNI3K) and rs1010553 (ITIH4) on childhood obesity risk. The best-associated GRS was significantly more associated with childhood obesity (OR = 1.21 per additional risk allele [95%CI:1.17-1.25], P = 4.8 × 10^-25) than the discovery GRS (OR = 1.05 per additional risk allele [95%CI:1.02-1.08], P = 8.0 × 10^-4), and was also associated with waist-to-hip ratio, fasting glucose, fasting insulin and triglyceride levels, the association being mediated by obesity. An overall depletion of obesity risk alleles was observed in Mexican children with normal weight when compared to GWAS discovery populations.

CONCLUSIONS

Our study indicates a partial transferability of GWAS obesity loci in Mexican children, and supports the pertinence of post-GWAS fine-mapping experiments in the admixed Mexican population.

Obesity genetics: insights from the Pakistani population

The Pakistani population is extensively diverse, indicating a genetic admixture of European and Central/West Asian migrants with indigenous South Asian gene pools. Pakistanis are organized in different ethnicities/castes based on cultural, linguistic and geographical origin. While Pakistan is facing a rapid nutritional transition, the rising prevalence of obesity is driving a growing burden of health complications and mortality. This represents a unique opportunity for the research community to study the interplay between obesogenic environmental changes and obesity predisposing genes in the time frame of one generation. This review recapitulates the ancestral origins of Pakistani population, the societal determinants of the rise in obesity and its governmental management. We describe the contribution of syndromic, monogenic non-syndromic and polygenic obesity genes identified in the Pakistani population. We then discuss the utility of gene identification approaches based on large consanguineous families and original gene × environment interaction study designs in discovering new obesity genes and causal pathways. Elucidation of the genetic basis of obesity in the Pakistani population may result in improved methods of obesity prevention and treatment globally.

Emerging roles of RNA-binding proteins in diabetes and their therapeutic potential in diabetic complications

Diabetes is a debilitating health care problem affecting 422 million people around the world. Diabetic patients suffer from multisystemic complications that can cause mortality and morbidity. Recent advancements in high-throughput next-generation RNA-sequencing and computational algorithms led to the discovery of aberrant posttranscriptional gene regulatory programs in diabetes. However, very little is known about how these regulatory programs are mis-regulated in diabetes. RNA-binding proteins (RBPs) are important regulators of posttranscriptional RNA networks, which are also dysregulated in diabetes. Human genetic studies provide new evidence that polymorphisms and mutations in RBPs are linked to diabetes. Therefore, we will discuss the emerging roles of RBPs in abnormal posttranscriptional gene expression in diabetes. Questions that will be addressed are: Which posttranscriptional mechanisms are disrupted in diabetes? Which RBPs are responsible for such changes under diabetic conditions? How are RBPs altered in diabetes? How does dysregulation of RBPs contribute to diabetes? Can we target RBPs using RNA-based methods to restore gene expression profiles in diabetic patients? Studying the evolving roles of RBPs in diabetes is critical not only for a comprehensive understanding of diabetes pathogenesis but also to design RNA-based therapeutic approaches for diabetic complications. WIREs RNA 2018, 9:e1459. doi: 10.1002/wrna.1459 This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Splicing Regulation/Alternative Splicing Translation > Translation Regulation.

An Efficient Test for Gene-Environment Interaction in Generalized Linear Mixed Models with Family Data

Gene-environment (GE) interaction has important implications in the etiology of complex diseases that are caused by a combination of genetic factors and environment variables. Several authors have developed GE analysis in the context of independent subjects or longitudinal data using a gene-set. In this paper, we propose to analyze GE interaction for discrete and continuous phenotypes in family studies by incorporating the relatedness among the relatives for each family into a generalized linear mixed model (GLMM) and by using a gene-based variance component test. In addition, we deal with collinearity problems arising from linkage disequilibrium among single nucleotide polymorphisms (SNPs) by considering their coefficients as random effects under the null model estimation. We show that the best linear unbiased predictor (BLUP) of such random effects in the GLMM is equivalent to the ridge regression estimator. This equivalence provides a simple method to estimate the ridge penalty parameter in comparison to other computationally-demanding estimation approaches based on cross-validation schemes. We evaluated the proposed test using simulation studies and applied it to real data from the Baependi Heart Study consisting of 76 families. Using our approach, we identified an interaction between BMI and the Peroxisome Proliferator Activated Receptor Gamma (PPARG) gene associated with diabetes.