A powerful approach to identify replicable variants in genome-wide association studies.

Replicability is the cornerstone of modern scientific research. Reliable identifications of genotype-phenotype associations that are significant in multiple genome-wide association studies (GWASs) provide stronger evidence for the findings. Current replicability analysis relies on the independence assumption among single-nucleotide polymorphisms (SNPs) and ignores the linkage disequilibrium (LD) structure. We show that such a strategy may produce either overly liberal or overly conservative results in practice. We develop an efficient method, ReAD, to detect replicable SNPs associated with the phenotype from two GWASs accounting for the LD structure. The local dependence structure of SNPs across two heterogeneous studies is captured by a four-state hidden Markov model (HMM) built on two sequences of p values. By incorporating information from adjacent locations via the HMM, our approach provides more accurate SNP significance rankings. ReAD is scalable, platform independent, and more powerful than existing replicability analysis methods with effective false discovery rate control. Through analysis of datasets from two asthma GWASs and two ulcerative colitis GWASs, we show that ReAD can identify replicable genetic loci that existing methods might otherwise miss.

Longitudinal analysis of epigenome-wide DNA methylation reveals novel loci associated with BMI change in East Asians.

BACKGROUND: Obesity is a global public health concern linked to chronic diseases such as cardiovascular disease and type 2 diabetes (T2D). Emerging evidence suggests that epigenetic modifications, particularly DNA methylation, may contribute to obesity. However, the molecular mechanism underlying the longitudinal change of BMI has not been well-explored, especially in East Asian populations. METHODS: This study performed a longitudinal epigenome-wide association analysis of DNA methylation to uncover novel loci associated with BMI change in 533 individuals across two Chinese cohorts with repeated DNA methylation and BMI measurements over four years. RESULTS: We identified three novel CpG sites (cg14671384, cg25540824, and cg10848724) significantly associated with BMI change. Two of the identified CpG sites were located in regions previously associated with body shape and basal metabolic rate. Annotation of the top 20 BMI change-associated CpGs revealed strong connections to obesity and T2D. Notably, these CpGs exhibited active regulatory roles and located in genes with high expression in the liver and digestive tract, suggesting a potential regulatory pathway from genome to phenotypes of energy metabolism and absorption via DNA methylation. Cross-sectional and longitudinal EWAS comparisons indicated different mechanisms between CpGs related to BMI and BMI change. CONCLUSION: This study enhances our understanding of the epigenetic dynamics underlying BMI change and emphasizes the value of longitudinal analyses in deciphering the complex interplay between epigenetics and obesity.

A perspective on genetic and polygenic risk scores-advances and limitations and overview of associated tools.

Polygenetic Risk Scores are used to evaluate an individual's vulnerability to developing specific diseases or conditions based on their genetic composition, by taking into account numerous genetic variations. This article provides an overview of the concept of Polygenic Risk Scores (PRS). We elucidate the historical advancements of PRS, their advantages and shortcomings in comparison with other predictive methods, and discuss their conceptual limitations in light of the complexity of biological systems. Furthermore, we provide a survey of published tools for computing PRS and associated resources. The various tools and software packages are categorized based on their technical utility for users or prospective developers. Understanding the array of available tools and their limitations is crucial for accurately assessing and predicting disease risks, facilitating early interventions, and guiding personalized healthcare decisions. Additionally, we also identify potential new avenues for future bioinformatic analyzes and advancements related to PRS.

Stacked neural network for predicting polygenic risk score.

In recent years, the utility of polygenic risk scores (PRS) in forecasting disease susceptibility from genome-wide association studies (GWAS) results has been widely recognised. Yet, these models face limitations due to overfitting and the potential overestimation of effect sizes in correlated variants. To surmount these obstacles, we devised the Stacked Neural Network Polygenic Risk Score (SNPRS). This novel approach synthesises outputs from multiple neural network models, each calibrated using genetic variants chosen based on diverse p-value thresholds. By doing so, SNPRS captures a broader array of genetic variants, enabling a more nuanced interpretation of the combined effects of these variants. We assessed the efficacy of SNPRS using the UK Biobank data, focusing on the genetic risks associated with breast and prostate cancers, as well as quantitative traits like height and BMI. We also extended our analysis to the Korea Genome and Epidemiology Study (KoGES) dataset. Impressively, our results indicate that SNPRS surpasses traditional PRS models and an isolated deep neural network in terms of accuracy, highlighting its promise in refining the efficacy and relevance of PRS in genetic studies.

Investigating neonatal health risk variables through cell-type specific methylome-wide association studies.

Adverse neonatal outcomes are a prevailing risk factor for both short- and long-term mortality and morbidity in infants. Given the importance of these outcomes, refining their assessment is paramount for improving prevention and care. Here we aim to enhance the assessment of these often correlated and multifaceted neonatal outcomes. To achieve this, we employ factor analysis to identify common and unique effects and further confirm these effects using criterion-related validity testing. This validation leverages methylome-wide profiles from neonatal blood. Specifically, we investigate nine neonatal health risk variables, including gestational age, Apgar score, three indicators of body size, jaundice, birth diagnosis, maternal preeclampsia, and maternal age. The methylomic profiles used for this research capture data from nearly all 28 million methylation sites in human blood, derived from the blood spot collected from 333 neonates, within 72 h post-birth. Our factor analysis revealed two common factors, size factor, that captured the shared effects of weight, head size, height, and gestational age and disease factor capturing the orthogonal shared effects of gestational age, combined with jaundice and birth diagnosis. To minimize false positives in the validation studies, validation was limited to variables with significant cumulative association as estimated through an in-sample replication procedure. This screening resulted in that the two common factors and the unique effects for gestational age, jaundice and Apgar were further investigated with full-scale cell-type specific methylome-wide association analyses. Highly significant, cell-type specific, associations were detected for both common effect factors and for Apgar. Gene Ontology analyses revealed multiple significant biologically relevant terms for the five fully investigated neonatal health risk variables. Given the established links between adverse neonatal outcomes and both immediate and long-term health, the distinct factor effects (representing the common and unique effects of the risk variables) and their biological profiles confirmed in our work, suggest their potential role as clinical biomarkers for assessing health risks and enhancing personalized care.

iGWAS: Image-based genome-wide association of self-supervised deep phenotyping of retina fundus images.

Existing imaging genetics studies have been mostly limited in scope by using imaging-derived phenotypes defined by human experts. Here, leveraging new breakthroughs in self-supervised deep representation learning, we propose a new approach, image-based genome-wide association study (iGWAS), for identifying genetic factors associated with phenotypes discovered from medical images using contrastive learning. Using retinal fundus photos, our model extracts a 128-dimensional vector representing features of the retina as phenotypes. After training the model on 40,000 images from the EyePACS dataset, we generated phenotypes from 130,329 images of 65,629 British White participants in the UK Biobank. We conducted GWAS on these phenotypes and identified 14 loci with genome-wide significance (p<5×10-8 and intersection of hits from left and right eyes). We also did GWAS on the retina color, the average color of the center region of the retinal fundus photos. The GWAS of retina colors identified 34 loci, 7 are overlapping with GWAS of raw image phenotype. Our results establish the feasibility of this new framework of genomic study based on self-supervised phenotyping of medical images.

A novel method for multiple phenotype association studies based on genotype and phenotype network.

Joint analysis of multiple correlated phenotypes for genome-wide association studies (GWAS) can identify and interpret pleiotropic loci which are essential to understand pleiotropy in diseases and complex traits. Meanwhile, constructing a network based on associations between phenotypes and genotypes provides a new insight to analyze multiple phenotypes, which can explore whether phenotypes and genotypes might be related to each other at a higher level of cellular and organismal organization. In this paper, we first develop a bipartite signed network by linking phenotypes and genotypes into a Genotype and Phenotype Network (GPN). The GPN can be constructed by a mixture of quantitative and qualitative phenotypes and is applicable to binary phenotypes with extremely unbalanced case-control ratios in large-scale biobank datasets. We then apply a powerful community detection method to partition phenotypes into disjoint network modules based on GPN. Finally, we jointly test the association between multiple phenotypes in a network module and a single nucleotide polymorphism (SNP). Simulations and analyses of 72 complex traits in the UK Biobank show that multiple phenotype association tests based on network modules detected by GPN are much more powerful than those without considering network modules. The newly proposed GPN provides a new insight to investigate the genetic architecture among different types of phenotypes. Multiple phenotypes association studies based on GPN are improved by incorporating the genetic information into the phenotype clustering. Notably, it might broaden the understanding of genetic architecture that exists between diagnoses, genes, and pleiotropy.

MR-GGI: accurate inference of gene-gene interactions using Mendelian randomization.

BACKGROUND: Researchers have long studied the regulatory processes of genes to uncover their functions. Gene regulatory network analysis is one of the popular approaches for understanding these processes, requiring accurate identification of interactions among the genes to establish the gene regulatory network. Advances in genome-wide association studies and expression quantitative trait loci studies have led to a wealth of genomic data, facilitating more accurate inference of gene-gene interactions. However, unknown confounding factors may influence these interactions, making their interpretation complicated. Mendelian randomization (MR) has emerged as a valuable tool for causal inference in genetics, addressing confounding effects by estimating causal relationships using instrumental variables. In this paper, we propose a new statistical method, MR-GGI, for accurately inferring gene-gene interactions using Mendelian randomization. RESULTS: MR-GGI applies one gene as the exposure and another as the outcome, using causal cis-single-nucleotide polymorphisms as instrumental variables in the inverse-variance weighted MR model. Through simulations, we have demonstrated MR-GGI's ability to control type 1 error and maintain statistical power despite confounding effects. MR-GGI performed the best when compared to other methods using the F1 score on the DREAM5 dataset. Additionally, when applied to yeast genomic data, MR-GGI successfully identified six clusters. Through gene ontology analysis, we have confirmed that each cluster in our study performs distinct functional roles by gathering genes with specific functions. CONCLUSION: These findings demonstrate that MR-GGI accurately inferences gene-gene interactions despite the confounding effects in real biological environments.

Causal relationship of interleukin-6 and its receptor on sarcopenia traits using mendelian randomization.

BACKGROUND: Previous research has extensively examined the role of interleukin 6 (IL-6) in sarcopenia. However, the presence of a causal relationship between IL-6, its receptor (IL-6R), and sarcopenia remains unclear. METHOD: In this study, we utilized summary-level data from genome-wide association studies (GWAS) focused on appendicular lean mass (ALM), hand grip strength, and walking pace. Single nucleotide polymorphisms (SNPs) were employed as genetic instruments for IL-6 and IL-6R to estimate the causal effect of sarcopenia traits. We adopted the Mendelian randomization (MR) approach to investigate these associations using the inverse variance weighted (IVW) method as the primary analytical approach. Additionally, we performed sensitivity analyses to validate the reliability of the MR results. RESULT: This study revealed a significant negative association between main IL-6R and eQTL IL-6R on the left grip strength were - 0.013 (SE = 0.004, p < 0.001) and -0.029 (SE = 0.007, p < 0.001), respectively. While for the right grip strength, the estimates were - 0.011 (SE = 0.001, p < 0.001) and - 0.021 (SE = 0.008, p = 0.005). However, no evidence of an association for IL-6R with ALM and walking pace. In addition, IL-6 did not affect sarcopenia traits. CONCLUSION: Our study findings suggest a negative association between IL-6R and hand grip strength. Additionally, targeting IL-6R may hold potential value as a therapeutic approach for the treatment of hand grip-related issues.

Integrated analysis of single cell-RNA sequencing and Mendelian randomization identifies lactate dehydrogenase B as a target of melatonin in ischemic stroke.

AIMS: Despite the success of single-cell RNA sequencing in identifying cellular heterogeneity in ischemic stroke, clarifying the mechanisms underlying these associations of differently expressed genes remains challenging. Several studies that integrate gene expression and gene expression quantitative trait loci (eQTLs) with genome wide-association study (GWAS) data to determine their causal role have been proposed. METHODS: Here, we combined Mendelian randomization (MR) framework and single cell (sc) RNA sequencing to study how differently expressed genes (DEGs) mediating the effect of gene expression on ischemic stroke. The hub gene was further validated in the in vitro model. RESULTS: We identified 2339 DEGs in 10 cell clusters. Among these DEGs, 58 genes were associated with the risk of ischemic stroke. After external validation with eQTL dataset, lactate dehydrogenase B (LDHB) is identified to be positively associated with ischemic stroke. The expression of LDHB has also been validated in sc RNA-seq with dominant expression in microglia and astrocytes, and melatonin is able to reduce the LDHB expression and activity in vitro ischemic models. CONCLUSION: Our study identifies LDHB as a novel biomarker for ischemic stroke via combining the sc RNA-seq and MR analysis.

Interpretation of 10 years of Alzheimer's disease genetic findings in the perspective of statistical heterogeneity.

Common genetic variants and susceptibility loci associated with Alzheimer's disease (AD) have been discovered through large-scale genome-wide association studies (GWAS), GWAS by proxy (GWAX) and meta-analysis of GWAS and GWAX (GWAS+GWAX). However, due to the very low repeatability of AD susceptibility loci and the low heritability of AD, these AD genetic findings have been questioned. We summarize AD genetic findings from the past 10 years and provide a new interpretation of these findings in the context of statistical heterogeneity. We discovered that only 17% of AD risk loci demonstrated reproducibility with a genome-wide significance of P < 5.00E-08 across all AD GWAS and GWAS+GWAX datasets. We highlighted that the AD GWAS+GWAX with the largest sample size failed to identify the most significant signals, the maximum number of genome-wide significant genetic variants or maximum heritability. Additionally, we identified widespread statistical heterogeneity in AD GWAS+GWAX datasets, but not in AD GWAS datasets. We consider that statistical heterogeneity may have attenuated the statistical power in AD GWAS+GWAX and may contribute to explaining the low repeatability (17%) of genome-wide significant AD susceptibility loci and the decreased AD heritability (40-2%) as the sample size increased. Importantly, evidence supports the idea that a decrease in statistical heterogeneity facilitates the identification of genome-wide significant genetic loci and contributes to an increase in AD heritability. Collectively, current AD GWAX and GWAS+GWAX findings should be meticulously assessed and warrant additional investigation, and AD GWAS+GWAX should employ multiple meta-analysis methods, such as random-effects inverse variance-weighted meta-analysis, which is designed specifically for statistical heterogeneity.

New Insights into the Genetic Basis of Lysine Accumulation in Rice Revealed by Multi-Model GWAS.

Lysine is an essential amino acid that cannot be synthesized in humans. Rice is a global staple food for humans but has a rather low lysine content. Identification of the quantitative trait nucleotides (QTNs) and genes underlying lysine content is crucial to increase lysine accumulation. In this study, five grain and three leaf lysine content datasets and 4,630,367 single nucleotide polymorphisms (SNPs) of 387 rice accessions were used to perform a genome-wide association study (GWAS) by ten statistical models. A total of 248 and 71 common QTNs associated with grain/leaf lysine content were identified. The accuracy of genomic selection/prediction RR-BLUP models was up to 0.85, and the significant correlation between the number of favorable alleles per accession and lysine content was up to 0.71, which validated the reliability and additive effects of these QTNs. Several key genes were uncovered for fine-tuning lysine accumulation. Additionally, 20 and 30 QTN-by-environment interactions (QEIs) were detected in grains/leaves. The QEI-sf0111954416 candidate gene LOC_Os01g21380 putatively accounted for gene-by-environment interaction was identified in grains. These findings suggested the application of multi-model GWAS facilitates a better understanding of lysine accumulation in rice. The identified QTNs and genes hold the potential for lysine-rich rice with a normal phenotype.

Genome-wide association analyses of breast cancer in women of African ancestry identify new susceptibility loci and improve risk prediction.

We performed genome-wide association studies of breast cancer including 18,034 cases and 22,104 controls of African ancestry. Genetic variants at 12 loci were associated with breast cancer risk (P < 5 × 10-8), including associations of a low-frequency missense variant rs61751053 in ARHGEF38 with overall breast cancer (odds ratio (OR) = 1.48) and a common variant rs76664032 at chromosome 2q14.2 with triple-negative breast cancer (TNBC) (OR = 1.30). Approximately 15.4% of cases with TNBC carried six risk alleles in three genome-wide association study-identified TNBC risk variants, with an OR of 4.21 (95% confidence interval = 2.66-7.03) compared with those carrying fewer than two risk alleles. A polygenic risk score (PRS) showed an area under the receiver operating characteristic curve of 0.60 for the prediction of breast cancer risk, which outperformed PRS derived using data from females of European ancestry. Our study markedly increases the population diversity in genetic studies for breast cancer and demonstrates the utility of PRS for risk prediction in females of African ancestry.

Dietary factors and Alzheimer's disease risk: a Mendelian randomization study.

BACKGROUND: Prior observational research has investigated the association between dietary patterns and Alzheimer's disease (AD) risk. Nevertheless, due to constraints in past observational studies, establishing a causal link between dietary habits and AD remains challenging. METHODS: Methodology involved the utilization of extensive cohorts sourced from publicly accessible genome-wide association study (GWAS) datasets of European descent for conducting Mendelian randomization (MR) analyses. The principal analytical technique utilized was the inverse-variance weighted (IVW) method. RESULTS: The MR analysis conducted in this study found no statistically significant causal association between 20 dietary habits and the risk of AD (All p > 0.05). These results were consistent across various MR methods employed, including MR-Egger, weighted median, simple mode, and weighted mode approaches. Moreover, there was no evidence of horizontal pleiotropy detected (All p > 0.05). CONCLUSION: In this MR analysis, our finding did not provide evidence to support the causal genetic relationships between dietary habits and AD risk.

Identifying potential drug targets for idiopathic pulmonary fibrosis: a mendelian randomization study based on the druggable genes.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) is a chronic fibrotic interstitial lung disease characterized by progressive dyspnea and decreased lung function, yet its exact etiology remains unclear. It is of great significance to discover new drug targets for IPF. METHODS: We obtained the cis-expression quantitative trait locus (cis-eQTL) of druggable genes from eQTLGen Consortium as exposure and the genome wide association study (GWAS) of IPF from the International IPF Genetics Consortium as outcomes to simulate the effects of drugs on IPF by employing mendelian randomization analysis. Then colocalization analysis was performed to calculate the probability of both cis-eQTL of druggable genes and IPF sharing a causal variant. For further validation, we conducted protein quantitative trait locus (pQTL) analysis to reaffirm our findings. RESULTS: The expression of 45 druggable genes was significantly associated with IPF susceptibility at FDR < 0.05. The expression of 23 and 15 druggable genes was significantly associated with decreased forced vital capacity (FVC) and diffusing capacity of the lungs for carbon monoxide (DLco) in IPF patients, respectively. IPF susceptibility and two significant genes (IL-7 and ABCB2) were likely to share a causal variant. The results of the pQTL analysis demonstrated that high levels of IL-7 in plasma are associated with a reduced risk of IPF (OR = 0.67, 95%CI: 0.47-0.97). CONCLUSION: IL-7 stands out as the most promising potential drug target to mitigate the risk of IPF. Our study not only sheds light on potential drug targets but also provides a direction for future drug development in IPF.

Genetic Variants Underlying Plasticity in Natural Populations of Spadefoot Toads: Environmental Assessment versus Phenotypic Response.

Many organisms facultatively produce different phenotypes depending on their environment, yet relatively little is known about the genetic bases of such plasticity in natural populations. In this study, we describe the genetic variation underlying an extreme form of plasticity--resource polyphenism--in Mexican spadefoot toad tadpoles, Spea multiplicata. Depending on their environment, these tadpoles develop into one of two drastically different forms: a carnivore morph or an omnivore morph. We collected both morphs from two ponds that differed in which morph had an adaptive advantage and performed genome-wide association studies of phenotype (carnivore vs. omnivore) and adaptive plasticity (adaptive vs. maladaptive environmental assessment). We identified four quantitative trait loci associated with phenotype and nine with adaptive plasticity, two of which exhibited signatures of minor allele dominance and two of which (one phenotype locus and one adaptive plasticity locus) did not occur as minor allele homozygotes. Investigations into the genetics of plastic traits in natural populations promise to provide novel insights into how such complex, adaptive traits arise and evolve.

Mitochondrial DNA: Inherent Complexities Relevant to Genetic Analyses.

Mitochondrial DNA (mtDNA) exhibits distinct characteristics distinguishing it from the nuclear genome, necessitating specific analytical methods in genetic studies. This comprehensive review explores the complex role of mtDNA in a variety of genetic studies, including genome-wide, epigenome-wide, and phenome-wide association studies, with a focus on its implications for human traits and diseases. Here, we discuss the structure and gene-encoding properties of mtDNA, along with the influence of environmental factors and epigenetic modifications on its function and variability. Particularly significant are the challenges posed by mtDNA's high mutation rate, heteroplasmy, and copy number variations, and their impact on disease susceptibility and population genetic analyses. The review also highlights recent advances in methodological approaches that enhance our understanding of mtDNA associations, advocating for refined genetic research techniques that accommodate its complexities. By providing a comprehensive overview of the intricacies of mtDNA, this paper underscores the need for an integrated approach to genetic studies that considers the unique properties of mitochondrial genetics. Our findings aim to inform future research and encourage the development of innovative methodologies to better interpret the broad implications of mtDNA in human health and disease.

The Landscape of Presence/Absence Variations during the Improvement of Rice.

Rice is one of the most important staple crops in the world; therefore, the improvement of rice holds great significance for enhancing agricultural production and addressing food security challenges. Although there have been numerous studies on the role of single-nucleotide polymorphisms (SNPs) in rice improvement with the development of next-generation sequencing technologies, research on the role of presence/absence variations (PAVs) in the improvement of rice is limited. In particular, there is a scarcity of studies exploring the traits and genes that may be affected by PAVs in rice. Here, we extracted PAVs utilizing resequencing data from 148 improved rice varieties distributed in Asia. We detected a total of 33,220 PAVs and found that the number of variations decreased gradually as the length of the PAVs increased. The number of PAVs was the highest on chromosome 1. Furthermore, we identified a 6 Mb hotspot region on chromosome 11 containing 1091 PAVs in which there were 29 genes related to defense responses. By conducting a genome-wide association study (GWAS) using PAV variation data and phenotypic data for five traits (flowering time, plant height, flag leaf length, flag leaf width, and panicle number) across all materials, we identified 186 significantly associated PAVs involving 20 cloned genes. A haplotype analysis and expression analysis of candidate genes revealed that important genes might be affected by PAVs, such as the flowering time gene OsSFL1 and the flag leaf width gene NAL1. Our work investigated the pattern in PAVs and explored important PAV key functional genes associated with agronomic traits. Consequently, these results provide potential and exploitable genetic resources for rice breeding.

Genetic Diversity and Genome-Wide Association Analysis of the Hulled/Naked Trait in a Barley Collection from Shanghai Agricultural Gene Bank.

Barley is one of the most important cereal crops in the world, and its value as a food is constantly being revealed, so the research into and the use of barley germplasm are very important for global food security. Although a large number of barley germplasm samples have been collected globally, their specific genetic compositions are not well understood, and in many cases their origins are even disputed. In this study, 183 barley germplasm samples from the Shanghai Agricultural Gene Bank were genotyped using genotyping-by-sequencing (GBS) technology, SNPs were identified and their genetic parameters were estimated, principal component analysis (PCA) was preformed, and the phylogenetic tree and population structure of the samples were also analyzed. In addition, a genome-wide association study (GWAS) was carried out for the hulled/naked grain trait, and a KASP marker was developed using an associated SNP. The results showed that a total of 181,906 SNPs were identified, and these barley germplasm samples could be roughly divided into three categories according to the phylogenetic analysis, which was generally consistent with the classification of the traits of row type and hulled/naked grain. Population structure analysis showed that the whole barley population could be divided into four sub-populations (SPs), the main difference from previous classifications being that the two-rowed and the hulled genotypes were sub-divided into two SPs. The GWAS analysis of the hulled/naked trait showed that many associated loci were unrelated to the Nud/nud locus, indicating that there might be new loci controlling the trait. A KASP marker was developed for one exon-type SNP on chromosome 7. Genotyping based on the KASP assay was consistent with that based on SNPs, indicating that the gene of this locus might be associated with the hulled/naked trait. The above work not only lays a good foundation for the future utilization of this barley germplasm population but it provides new loci and candidate genes for the hulled/naked trait.