Missing heritability and strategies for finding the underlying causes of complex disease. Nat Rev Genet. 2010;11:446-450. [PMID: 20479774 DOI: 10.1038/nrg2809]

Marginal interaction test for detecting interactions between genetic marker sets and environment in genome-wide studies

As human complex diseases are influenced by the interaction between genetics and the environment, identifying gene-environment interactions (G×E) is crucial for understanding disease mechanisms and predicting risk. Developing robust quantitative tools for G×E analysis can enhance the study of complex diseases. However, many existing methods that explore G×E focus on the interplay between an environmental factor and genetic variants, exclusively for common or rare variants. In this study, we developed MAGEIT_RAN and MAGEIT_FIX to identify interactions between an environmental factor and a set of genetic markers, including both rare and common variants, based on the MinQue for Summary statistics. The genetic main effects in MAGEIT_RAN and MAGEIT_FIX are modeled as random and fixed effects, respectively. Simulation studies showed that both tests had type I error under control, with MAGEIT_RAN being the most powerful test. Applying MAGEIT to a genome-wide analysis of gene-alcohol interactions on hypertension and seated systolic blood pressure in the Multiethnic Study of Atherosclerosis revealed genes like EIF2AK2, CCNDBP1, and EPB42 influencing blood pressure through alcohol interaction. Pathway analysis identified 1 apoptosis and survival pathway involving PKR and 2 signal transduction pathways associated with hypertension and alcohol intake, demonstrating MAGEIT_RAN's ability to detect biologically relevant gene-environment interactions.

Powerful Rare-Variant Association Analysis of Secondary Phenotypes

Most genome-wide association studies are based on case-control designs, which provide abundant resources for secondary phenotype analyses. However, such studies suffer from biased sampling of primary phenotypes, and the traditional statistical methods can lead to seriously distorted analysis results when they are applied to secondary phenotypes without accounting for the biased sampling mechanism. To our knowledge, there are no statistical methods specifically tailored for rare variant association analysis with secondary phenotypes. In this article, we proposed two novel joint test statistics for identifying secondary-phenotype-associated rare variants based on prospective likelihood and retrospective likelihood, respectively. We also exploit the assumption of gene-environment independence in retrospective likelihood to improve the statistical power and adopt a two-step strategy to balance statistical power and robustness. Simulations and a real-data application are conducted to demonstrate the superior performance of our proposed methods.

Immunological roads diverged: mapping tuberculosis outcomes in mice

The journey from phenotypic observation to causal genetic mechanism is a long and challenging road. For pathogens like Mycobacterium tuberculosis (Mtb), which causes tuberculosis (TB), host-pathogen coevolution has spanned millennia, costing millions of human lives. Mammalian models can systematically recapitulate host genetic variation, producing a spectrum of disease outcomes. Leveraging genome sequences and deep phenotyping data from infected mouse genetic reference populations (GRPs), quantitative trait locus (QTL) mapping approaches have successfully identified host genomic regions associated with TB phenotypes. Here, we review the ongoing optimization of QTL mapping study design alongside advances in mouse GRPs. These next-generation resources and approaches have enabled identification of novel host-pathogen interactions governing one of the most prevalent infectious diseases in the world today.

Scrutinizing neurodegenerative diseases: decoding the complex genetic architectures through a multi-omics lens

Neurodegenerative diseases present complex genetic architectures, reflecting a continuum from monogenic to oligogenic and polygenic models. Recent advances in multi-omics data, coupled with systems genetics, have significantly refined our understanding of how these data impact neurodegenerative disease mechanisms. To contextualize these genetic discoveries, we provide a comprehensive critical overview of genetic architecture concepts, from Mendelian inheritance to the latest insights from oligogenic and omnigenic models. We explore the roles of common and rare genetic variants, gene-gene and gene-environment interactions, and epigenetic influences in shaping disease phenotypes. Additionally, we emphasize the importance of multi-omics layers including genomic, transcriptomic, proteomic, epigenetic, and metabolomic data in elucidating the molecular mechanisms underlying neurodegeneration. Special attention is given to missing heritability and the contribution of rare variants, particularly in the context of pleiotropy and network pleiotropy. We examine the application of single-cell omics technologies, transcriptome-wide association studies, and epigenome-wide association studies as key approaches for dissecting disease mechanisms at tissue- and cell-type levels. Our review introduces the OmicPeak Disease Trajectory Model, a conceptual framework for understanding the genetic architecture of neurodegenerative disease progression, which integrates multi-omics data across biological layers and time points. This review highlights the critical importance of adopting a systems genetics approach to unravel the complex genetic architecture of neurodegenerative diseases. Finally, this emerging holistic understanding of multi-omics data and the exploration of the intricate genetic landscape aim to provide a foundation for establishing more refined genetic architectures of these diseases, enhancing diagnostic precision, predicting disease progression, elucidating pathogenic mechanisms, and refining therapeutic strategies for neurodegenerative conditions.

Research Review: A review of the past decade of family and genomic studies on adolescent mental health

BACKGROUND

Mental health problems and traits capturing psychopathology are common and often begin during adolescence. Decades of twin studies indicate that genetic factors explain around 50% of individual differences in adolescent psychopathology. In recent years, significant advances, particularly in genomics, have moved this work towards more translational findings.

METHODS

This review provides an overview of the past decade of genetically sensitive studies on adolescent development, covering both family and genomic studies in adolescents aged 10-24 years. We focus on five research themes: (1) co-occurrence or comorbidity between psychopathologies, (2) stability and change over time, (3) intergenerational transmission, (4) gene-environment interplay, and (5) psychological treatment outcomes.

RESULTS

First, research shows that much of the co-occurrence of psychopathologies in adolescence is explained by genetic factors, with widespread pleiotropic influences on many traits. Second, stability in psychopathology across adolescence is largely explained by persistent genetic influences, whereas change is explained by emerging genetic and environmental influences. Third, contemporary twin-family studies suggest that different co-occurring genetic and environmental mechanisms may account for the intergenerational transmission of psychopathology, with some differences across psychopathologies. Fourth, genetic influences on adolescent psychopathology are correlated with a wide range of environmental exposures. However, the extent to which genetic factors interact with the environment remains unclear, as findings from both twin and genomic studies are inconsistent. Finally, a few studies suggest that genetic factors may play a role in psychological treatment response, but these findings have not yet been replicated.

CONCLUSIONS

Genetically sensitive research on adolescent psychopathology has progressed significantly in the past decade, with family and twin findings starting to be replicated at the genomic level. However, important gaps remain in the literature, and we conclude by providing suggestions of research questions that still need to be addressed.

PsyRiskMR: a comprehensive resource for identifying psychiatric disorders risk factors through Mendelian randomization

BACKGROUND

Psychiatric disorders pose an enormous economic and emotional burden on individuals, their families and society. Given that the current analysis of the pathogenesis of psychiatric disorders remains challenging and time-consuming, elucidating the modifiable risk factors becomes crucial for the diagnosis and management of psychiatric disorders. However, inferring the causal risk factors in these disorders from disparate data sources is challenging due to constraints in data collection and analytical capabilities.

METHODS

By leveraging the largest available genome-wide association studies (GWAS) summary statistics for ten psychiatric disorders and compiling an extensive set of risk factor datasets, including 71 psychiatric disorders-specific phenotypes, 3,935 brain imaging traits, and over 30 brain tissue/cell-specific xQTL datasets (covering 6 types of QTLs), we performed comprehensive Mendelian randomization (MR) analyses to explore the potential causal links between various exposures and psychiatric outcomes using genetic variants as instrumental variables.

RESULTS

After Bonferroni correction for multiple testing, we identified multiple potential risk factors for psychiatric disorders (including phenotypic level and molecular level traits), and provided robust MR evidence supporting these associations utilizing rigorous sensitivity analyses and colocalization analyses. Furthermore. we have established the PsyRiskMR database (http://bioinfo.ahu.edu.cn/PsyRiskMR/), which serves as an interactive platform for showcasing and querying risk factors for psychiatric disorders.

CONCLUSIONS

Our study offered a user-friendly PsyRiskMR database for the research community to browse, search, and download all MR results, potentially revealing new insights into the biological etiology of psychiatric disorders.

Understanding genetic variants in context

Over the last three decades, human genetics has gone from dissecting high-penetrance Mendelian diseases to discovering the vast and complex genetic etiology of common human diseases. In tackling this complexity, scientists have discovered the importance of numerous genetic processes - most notably functional regulatory elements - in the development and progression of these diseases. Simultaneously, scientists have increasingly used multiplex assays of variant effect to systematically phenotype the cellular consequences of millions of genetic variants. In this article, we argue that the context of genetic variants - at all scales, from other genetic variants and gene regulation to cell biology to organismal environment - are critical components of how we can employ genomics to interpret these variants, and ultimately treat these diseases. We describe approaches to extend existing experimental assays and computational approaches to examine and quantify the importance of this context, including through causal analytic approaches. Having a unified understanding of the molecular, physiological, and environmental processes governing the interpretation of genetic variants is sorely needed for the field, and this perspective argues for feasible approaches by which the combined interpretation of cellular, animal, and epidemiological data can yield that knowledge.

A phenome-wide association study of methylated GC-rich repeats identifies a GCC repeat expansion in AFF3 associated with intellectual disability

GC-rich tandem repeat expansions (TREs) are often associated with DNA methylation, gene silencing and folate-sensitive fragile sites, and underlie several congenital and late-onset disorders. Through a combination of DNA-methylation profiling and tandem repeat genotyping, we identified 24 methylated TREs and investigated their effects on human traits using phenome-wide association studies in 168,641 individuals from the UK Biobank, identifying 156 significant TRE-trait associations involving 17 different TREs. Of these, a GCC expansion in the promoter of AFF3 was associated with a 2.4-fold reduced probability of completing secondary education, an effect size comparable to several recurrent pathogenic microdeletions. In a cohort of 6,371 probands with neurodevelopmental problems of suspected genetic etiology, we observed a significant enrichment of AFF3 expansions compared with controls. With a population prevalence that is at least fivefold higher than the TRE that causes fragile X syndrome, AFF3 expansions represent a major cause of neurodevelopmental delay.

Overcoming the "feast or famine" effect: improved interaction testing in genome-wide association studies

In genetic association analysis of complex traits, detection of interaction (either GxG or GxE) can help to elucidate the genetic architecture and biological mechanisms underlying the trait. Detection of interaction in a genome-wide interaction study (GWIS) can be methodologically challenging for various reasons, including a high burden of multiple comparisons when testing for epistasis between all possible pairs of a set of genomewide variants, as well as heteroscedasticity effects occurring in the presence of GxG or GxE interaction. In this paper, we address the problem of an even more striking phenomenon that we call the "feast or famine" effect that occurs when testing interaction in a genomewide context. We show that in any given GxE GWIS, the type 1 error of standard interaction tests performed genomewide can vary widely from the nominal level, where the actual type 1 error in any given GWIS varies as a predictable function of the observed trait and environmental values. Using standard methods, some GWISs will have systematically underinflated p-values ("feast"), and others will have systematically overinflated p-values ("famine"), which can lead to false detection of interaction, reduced power, inconsistent results across studies, and failure to replicate true signal. This startling phenomenon is specific to detection of interaction in a GWIS, and it may partly explain why such detection has often proved challenging and difficult to replicate. We show that the feast or famine effect occurs across a wide range of GxE analysis methods, including but not limited to (1) testing interaction in a linear or linear mixed model (LMM) using standard approaches such as t-tests/Wald tests, likelihood ratio tests, or score tests; (2) doing a combined interaction-association test in a linear model or LMM using standard approaches such as F-tests or likelihood ratio tests; (3) testing interaction with multiple environments or multiple SNPs, where these are modeled as random effects in a LMM using standard approaches; (4) performing tests of interaction in a GWIS where significance is assessed using permutation of the trait residuals. We show theoretically that the key cause of this phenomenon is which variables are conditioned on in the analysis, and this suggests an approach to correct the problem by changing the way the conditioning is done. Using this insight, we have developed the TINGA method to adjust the interaction test statistics to make their p-values closer to uniform under the null hypothesis. In simulations we show that TINGA both controls type 1 error and improves power. TINGA allows for covariates and population structure through use of a linear mixed model and accounts for heteroscedasticity. We apply TINGA to detection of epistasis in a study of flowering time in Arabidopsis thaliana.

[Single nucleotide polymorphism heritability of non-syndromic cleft lip with or without cleft palate in Chinese population]

OBJECTIVE

To delve into the intricate relationship between common genetic variations across the entire genome and the risk of non-syndromic cleft lip with or without cleft palate (NSCL/P).

METHODS

Utilizing summary statistics data from genome-wide association studies (GWAS), a thorough investigation to evaluate the impact of common variations on the genome were undertook. This involved assessing single nucleotide polymorphism (SNP) heritability across the entire genome, as well as within specific genomic regions. To ensure the robustness of our analysis, stringent quality control measures were applied to the GWAS summary statistics data. Criteria for inclusion encompassed the absence of missing values, a minor allele frequency ≥1%, P-values falling within the range of 0 to 1, and clear SNP strand orientation. SNP meeting these stringent criteria were then meticulously included in our analysis. The SNP heritability of NSCL/P was calculated using linkage disequilibrium score regression. Additionally, hierarchical linkage disequilibrium score regression to partition SNP heritability within coding regions, promoters, introns, enhancers, and super enhancers were employed, and the enrichment levels within different genomic regions using LDSC (v1.0.1) software were further elucidated.

RESULTS

Our study drew upon GWAS summary statistics data obtained from 806 NSCL/P trios, comprising a total of 2 418 individuals from the Chinese population. Following rigorous quality control procedures, 490 593 out of 492 993 SNP were deemed suitable for inclusion in SNP heritability calculations. The observed SNP heritability of NSCL/P was 0.55 (95%CI: 0.28-0.82). Adjusting for the elevated disease pre-valence within our sample, the SNP heritability scaled down to 0.37 (95%CI: 0.19-0.55) based on the prevalence observed in the general Chinese population. Notably, our enrichment analysis unveiled significant enrichment of SNP heritability within enhancer regions (15.70, P=0.04) and super enhancer regions (3.18, P=0.03).

CONCLUSION

Our study sheds light on the intricate interplay between common genetic variations and the risk of NSCL/P in the Chinese population. By elucidating the SNP heritability landscape across different genomic regions, we contribute valuable insights into the genetic basis of NSCL/P. The significant enrichment of SNP heritability within enhancer and super enhancer regions underscores the potential role of these regulatory elements in shaping the genetic susceptibility to NSCL/P. This paves the way for further research aimed at uncovering novel genetic pathogenic factors underlying NSCL/P pathogenesis.

A scalable adaptive quadratic kernel method for interpretable epistasis analysis in complex traits

Our knowledge of the contribution of genetic interactions (epistasis) to variation in human complex traits remains limited, partly due to the lack of efficient, powerful, and interpretable algorithms to detect interactions. Recently proposed approaches for set-based association tests show promise in improving the power to detect epistasis by examining the aggregated effects of multiple variants. Nevertheless, these methods either do not scale to large Biobank data sets or lack interpretability. We propose QuadKAST, a scalable algorithm focused on testing pairwise interaction effects (quadratic effects) within small to medium-sized sets of genetic variants (window size ≤100) on a trait and provide quantified interpretation of these effects. Comprehensive simulations show that QuadKAST is well-calibrated. Additionally, QuadKAST is highly sensitive in detecting loci with epistatic signals and accurate in its estimation of quadratic effects. We applied QuadKAST to 52 quantitative phenotypes measured in ≈300,000 unrelated white British individuals in the UK Biobank to test for quadratic effects within each of 9515 protein-coding genes. We detect 32 trait-gene pairs across 17 traits and 29 genes that demonstrate statistically significant signals of quadratic effects (accounting for the number of genes and traits tested). Across these trait-gene pairs, the proportion of trait variance explained by quadratic effects is comparable to additive effects, with five pairs having a ratio >1. Our method enables the detailed investigation of epistasis on a large scale, offering new insights into its role and importance.

BLUPmrMLM: A Fast mrMLM Algorithm in Genome-wide Association Studies

Multilocus genome-wide association study has become the state-of-the-art tool for dissecting the genetic architecture of complex and multiomic traits. However, most existing multilocus methods require relatively long computational time when analyzing large datasets. To address this issue, in this study, we proposed a fast mrMLM method, namely, best linear unbiased prediction multilocus random-SNP-effect mixed linear model (BLUPmrMLM). First, genome-wide single-marker scanning in mrMLM was replaced by vectorized Wald tests based on the best linear unbiased prediction (BLUP) values of marker effects and their variances in BLUPmrMLM. Then, adaptive best subset selection (ABESS) was used to identify potentially associated markers on each chromosome to reduce computational time when estimating marker effects via empirical Bayes. Finally, shared memory and parallel computing schemes were used to reduce the computational time. In simulation studies, BLUPmrMLM outperformed GEMMA, EMMAX, mrMLM, and FarmCPU as well as the control method (BLUPmrMLM with ABESS removed), in terms of computational time, power, accuracy for estimating quantitative trait nucleotide positions and effects, false positive rate, false discovery rate, false negative rate, and F1 score. In the reanalysis of two large rice datasets, BLUPmrMLM significantly reduced the computational time and identified more previously reported genes, compared with the aforementioned methods. This study provides an excellent multilocus model method for the analysis of large-scale and multiomic datasets. The software mrMLM v5.1 is available at BioCode (https://ngdc.cncb.ac.cn/biocode/tool/BT007388) or GitHub (https://github.com/YuanmingZhang65/mrMLM).

The effect of family structure on the still-missing heritability and genomic prediction accuracy of type 2 diabetes

This study aims to assess the effect of familial structures on the still-missing heritability estimate and prediction accuracy of Type 2 Diabetes (T2D) using pedigree estimated risk values (ERV) and genomic ERV. We used 11,818 individuals (T2D cases: 2,210) with genotype (649,932 SNPs) and pedigree information from the ongoing periodic cohort study of the Iranian population project. We considered three different familial structure scenarios, including (i) all families, (ii) all families with ≥ 1 generation, and (iii) families with ≥ 1 generation in which both case and control individuals are presented. Comprehensive simulation strategies were implemented to quantify the difference between estimates of [Formula: see text] and [Formula: see text]. A proportion of still-missing heritability in T2D could be explained by overestimation of pedigree-based heritability due to the presence of families with individuals having only one of the two disease statuses. Our research findings underscore the significance of including families with only case/control individuals in cohort studies. The presence of such family structures (as observed in scenarios i and ii) contributes to a more accurate estimation of disease heritability, addressing the underestimation that was previously overlooked in prior research. However, when predicting disease risk, the absence of these families (as seen in scenario iii) can yield the highest prediction accuracy and the strongest correlation with Polygenic Risk Scores. Our findings represent the first evidence of the important contribution of familial structure for heritability estimations and genomic prediction studies in T2D.

Short Tandem Repeats in the era of next-generation sequencing: from historical loci to population databases

In this study, we explore the landscape of short tandem repeats (STRs) within the human genome through the lens of evolving technologies to detect genomic variations. STRs, which encompass approximately 3% of our genomic DNA, are crucial for understanding human genetic diversity, disease mechanisms, and evolutionary biology. The advent of high-throughput sequencing methods has revolutionized our ability to accurately map and analyze STRs, highlighting their significance in genetic disorders, forensic science, and population genetics. We review the current available methodologies for STR analysis, the challenges in interpreting STR variations across different populations, and the implications of STRs in medical genetics. Our findings underscore the urgent need for comprehensive STR databases that reflect the genetic diversity of global populations, facilitating the interpretation of STR data in clinical diagnostics, genetic research, and forensic applications. This work sets the stage for future studies aimed at harnessing STR variations to elucidate complex genetic traits and diseases, reinforcing the importance of integrating STRs into genetic research and clinical practice.

Chemical Synthesis of Human Proteoforms and Application in Biomedicine

Limited understanding of human proteoforms with complex posttranslational modifications and the underlying mechanisms poses a major obstacle to research on human health and disease. This Outlook discusses opportunities and challenges of de novo chemical protein synthesis in human proteoform studies. Our analysis suggests that to develop a comprehensive, robust, and cost-effective methodology for chemical synthesis of various human proteoforms, new chemistries of the following types need to be developed: (1) easy-to-use peptide ligation chemistries allowing more efficient de novo synthesis of protein structural domains, (2) robust temporary structural support strategies for ligation and folding of challenging targets, and (3) efficient transpeptidative protein domain-domain ligation methods for multidomain proteins. Our analysis also indicates that accurate chemical synthesis of human proteoforms can be applied to the following aspects of biomedical research: (1) dissection and reconstitution of the proteoform interaction networks, (2) structural mechanism elucidation and functional analysis of human proteoform complexes, and (3) development and evaluation of drugs targeting human proteoforms. Overall, we suggest that through integrating chemical protein synthesis with in vivo functional analysis, mechanistic biochemistry, and drug development, synthetic chemistry would play a pivotal role in human proteoform research and facilitate the development of precision diagnostics and therapeutics.

Genetic Risk of Second Malignant Neoplasm after Childhood Cancer Treatment: A Systematic Review

Second malignant neoplasm (SMN) is one of the most severe long-term risks for childhood cancer survivors (CCS), significantly impacting long-term patient survival. While radiotherapy and chemotherapy are known risk factors, the observed inter-individual variability suggests a genetic component contributing to the risk of SMN. This article aims to conduct a systematic review of genetic factors implicated in the SMN risk among CCS. Searches were performed in PubMed, Scopus, and Web of Sciences. Eighteen studies were included (eleven candidate gene studies, three genome-wide association studies, and four whole exome/genome sequencing studies). The included studies were based on different types of first cancers, investigated any or specific types of SMN, and focused mainly on genes involved in drug metabolism and DNA repair pathways. These differences in study design and methods used to characterize genetic variants limit the scope of the results and highlight the need for further extensive and standardized investigations. However, this review provides a valuable compilation of SMN risk-associated variants and genes, facilitating efficient replication and advancing our understanding of the genetic basis for this major risk for CCS.

Combining Hyperspectral Techniques and Genome-Wide Association Studies to Predict Peanut Seed Vigor and Explore Associated Genetic Loci

Seed vigor significantly affects peanut breeding and agricultural yield by influencing seed germination and seedling growth and development. Traditional vigor testing methods are inadequate for modern high-throughput assays. Although hyperspectral technology shows potential for monitoring various crop traits, its application in predicting peanut seed vigor is still limited. This study developed and validated a method that combines hyperspectral technology with genome-wide association studies (GWAS) to achieve high-throughput detection of seed vigor and identify related functional genes. Hyperspectral phenotyping data and physiological indices from different peanut seed populations were used as input data to construct models using machine learning regression algorithms to accurately monitor changes in vigor. Model-predicted phenotypic data from 191 peanut varieties were used in GWAS, gene-based association studies, and haplotype analyses to screen for functional genes. Real-time fluorescence quantitative PCR (qPCR) was used to analyze the expression of functional genes in three high-vigor and three low-vigor germplasms. The results indicated that the random forest and support vector machine models provided effective phenotypic data. We identified Arahy.VMLN7L and Arahy.7XWF6F, with Arahy.VMLN7L negatively regulating seed vigor and Arahy.7XWF6F positively regulating it, suggesting distinct regulatory mechanisms. This study confirms that GWAS based on hyperspectral phenotyping reveals genetic relationships in seed vigor levels, offering novel insights and directions for future peanut breeding, accelerating genetic improvements, and boosting agricultural yields. This approach can be extended to monitor and explore germplasms and other key variables in various crops.

Genome-wide analyses reveal a potential role for the MAPT, MOBP, and APOE loci in sporadic frontotemporal dementia

Frontotemporal dementia (FTD) is the second most common cause of early-onset dementia after Alzheimer disease (AD). Efforts in the field mainly focus on familial forms of disease (fFTDs), while studies of the genetic etiology of sporadic FTD (sFTD) have been less common. In the current work, we analyzed 4,685 sFTD cases and 15,308 controls looking for common genetic determinants for sFTD. We found a cluster of variants at the MAPT (rs199443; p = 2.5 × 10-12, OR = 1.27) and APOE (rs6857; p = 1.31 × 10-12, OR = 1.27) loci and a candidate locus on chromosome 3 (rs1009966; p = 2.41 × 10-8, OR = 1.16) in the intergenic region between RPSA and MOBP, contributing to increased risk for sFTD through effects on expression and/or splicing in brain cortex of functionally relevant in-cis genes at the MAPT and RPSA-MOBP loci. The association with the MAPT (H1c clade) and RPSA-MOBP loci may suggest common genetic pleiotropy across FTD and progressive supranuclear palsy (PSP) (MAPT and RPSA-MOBP loci) and across FTD, AD, Parkinson disease (PD), and cortico-basal degeneration (CBD) (MAPT locus). Our data also suggest population specificity of the risk signals, with MAPT and APOE loci associations mainly driven by Central/Nordic and Mediterranean Europeans, respectively. This study lays the foundations for future work aimed at further characterizing population-specific features of potential FTD-discriminant APOE haplotype(s) and the functional involvement and contribution of the MAPT H1c haplotype and RPSA-MOBP loci to pathogenesis of sporadic forms of FTD in brain cortex.

Low-input PacBio sequencing generates high-quality individual fly genomes and characterizes mutational processes

Long-read sequencing, exemplified by PacBio, revolutionizes genomics, overcoming challenges like repetitive sequences. However, the high DNA requirement ( > 1 µg) is prohibitive for small organisms. We develop a low-input (100 ng), low-cost, and amplification-free library-generation method for PacBio sequencing (LILAP) using Tn5-based tagmentation and DNA circularization within one tube. We test LILAP with two Drosophila melanogaster individuals, and generate near-complete genomes, surpassing preexisting single-fly genomes. By analyzing variations in these two genomes, we characterize mutational processes: complex transpositions (transposon insertions together with extra duplications and/or deletions) prefer regions characterized by non-B DNA structures, and gene conversion of transposons occurs on both DNA and RNA levels. Concurrently, we generate two complete assemblies for the endosymbiotic bacterium Wolbachia in these flies and similarly detect transposon conversion. Thus, LILAP promises a broad PacBio sequencing adoption for not only mutational studies of flies and their symbionts but also explorations of other small organisms or precious samples.

Discovering non-additive heritability using additive GWAS summary statistics

LD score regression (LDSC) is a method to estimate narrow-sense heritability from genome-wide association study (GWAS) summary statistics alone, making it a fast and popular approach. In this work, we present interaction-LD score (i-LDSC) regression: an extension of the original LDSC framework that accounts for interactions between genetic variants. By studying a wide range of generative models in simulations, and by re-analyzing 25 well-studied quantitative phenotypes from 349,468 individuals in the UK Biobank and up to 159,095 individuals in BioBank Japan, we show that the inclusion of a cis-interaction score (i.e. interactions between a focal variant and proximal variants) recovers genetic variance that is not captured by LDSC. For each of the 25 traits analyzed in the UK Biobank and BioBank Japan, i-LDSC detects additional variation contributed by genetic interactions. The i-LDSC software and its application to these biobanks represent a step towards resolving further genetic contributions of sources of non-additive genetic effects to complex trait variation.

Genome-wide association studies (GWAS) and post-GWAS analyses of impulsivity: A systematic review

Impulsivity is related to a host of mental and behavioral problems. It is a complex construct with many different manifestations, most of which are heritable. The genetic compositions of these impulsivity manifestations, however, remain unclear. A number of genome-wide association studies (GWAS) and post-GWAS analyses have tried to address this issue. We conducted a systematic review of all GWAS and post-GWAS analyses of impulsivity published up to December 2023. Available data suggest that single nucleotide polymorphisms (SNPs) in more than a dozen of genes (e.g., CADM2, CTNNA2, GPM6B) are associated with different measures of impulsivity at genome-wide significant levels. Post-GWAS analyses further show that different measures of impulsivity are subject to different degrees of genetic influence, share few genetic variants, and have divergent genetic overlap with basic personality traits such as extroversion and neuroticism, cognitive ability, psychiatric disorders, substance use, and obesity. These findings shed light on controversies in the conceptualization and measurement of impulsivity, while providing new insights on the underlying mechanisms that yoke impulsivity to psychopathology.

Naturally segregating genetic variants contribute to thermal tolerance in a Drosophila melanogaster model system

Thermal tolerance is a fundamental physiological complex trait for survival in many species. For example, everyday tasks such as foraging, finding a mate, and avoiding predation are highly dependent on how well an organism can tolerate extreme temperatures. Understanding the general architecture of the natural variants within the genes that control this trait is of high importance if we want to better comprehend thermal physiology. Here, we take a multipronged approach to further dissect the genetic architecture that controls thermal tolerance in natural populations using the Drosophila Synthetic Population Resource as a model system. First, we used quantitative genetics and Quantitative Trait Loci mapping to identify major effect regions within the genome that influences thermal tolerance, then integrated RNA-sequencing to identify differences in gene expression, and lastly, we used the RNAi system to (1) alter tissue-specific gene expression and (2) functionally validate our findings. This powerful integration of approaches not only allows for the identification of the genetic basis of thermal tolerance but also the physiology of thermal tolerance in a natural population, which ultimately elucidates thermal tolerance through a fitness-associated lens.

BridGE: a pathway-based analysis tool for detecting genetic interactions from GWAS

Genetic interactions have the potential to modulate phenotypes, including human disease. In principle, genome-wide association studies (GWAS) provide a platform for detecting genetic interactions; however, traditional methods for identifying them, which tend to focus on testing individual variant pairs, lack statistical power. In this protocol, we describe a novel computational approach, called Bridging Gene sets with Epistasis (BridGE), for discovering genetic interactions between biological pathways from GWAS data. We present a Python-based implementation of BridGE along with instructions for its application to a typical human GWAS cohort. The major stages include initial data processing and quality control, construction of a variant-level genetic interaction network, measurement of pathway-level genetic interactions, evaluation of statistical significance using sample permutations and generation of results in a standardized output format. The BridGE software pipeline includes options for running the analysis on multiple cores and multiple nodes for users who have access to computing clusters or a cloud computing environment. In a cluster computing environment with 10 nodes and 100 GB of memory per node, the method can be run in less than 24 h for typical human GWAS cohorts. Using BridGE requires knowledge of running Python programs and basic shell script programming experience.

Gene-based association tests in family samples using GWAS summary statistics

Genome-wide association studies (GWAS) have led to rapid growth in detecting genetic variants associated with various phenotypes. Owing to a great number of publicly accessible GWAS summary statistics, and the difficulty in obtaining individual-level genotype data, many existing gene-based association tests have been adapted to require only GWAS summary statistics rather than individual-level data. However, these association tests are restricted to unrelated individuals and thus do not apply to family samples directly. Moreover, due to its flexibility and effectiveness, the linear mixed model has been increasingly utilized in GWAS to handle correlated data, such as family samples. However, it remains unknown how to perform gene-based association tests in family samples using the GWAS summary statistics estimated from the linear mixed model. In this study, we show that, when family size is negligible compared to the total sample size, the diagonal block structure of the kinship matrix makes it possible to approximate the correlation matrix of marginal Z scores by linkage disequilibrium matrix. Based on this result, current methods utilizing summary statistics for unrelated individuals can be directly applied to family data without any modifications. Our simulation results demonstrate that this proposed strategy controls the type 1 error rate well in various situations. Finally, we exemplify the usefulness of the proposed approach with a dental caries GWAS data set.

Identification of a strong genetic risk factor for major depressive disorder in the human virome

The heritability of major depressive disorder (MDD) is reportedly 30-50%. However, the genetic basis of its heritability remains unknown. Within SITH-1, a risk factor for MDD in human herpesvirus 6B (HHV-6B), we discovered a gene polymorphism with a large odds ratio for an association with MDD. It was a sequence whose number of repeats was inversely correlated with SITH-1 expression. This number was significantly lower in MDD patients. Rates for 17 or fewer repeats of the sequence were 67.9% for MDD and 28.6% for normal controls, with an odds ratio of 5.28. For patients with 17 or less repeats, the rate for presence of another MDD patient in their families was 47.4%, whereas there were no MDD patients in the families of patients with more than 17 repeats. Since HHV-6B is transmitted primarily mother to child and within families and persists for life, this gene polymorphism could potentially influence heritability of MDD.

Integration of whole-exome sequencing and structural neuroimaging analysis in major depressive disorder: a joint study

Major depressive disorder (MDD) is a common mental illness worldwide and is triggered by an intricate interplay between environmental and genetic factors. Although there are several studies on common variants in MDD, studies on rare variants are relatively limited. In addition, few studies have examined the genetic contributions to neurostructural alterations in MDD using whole-exome sequencing (WES). We performed WES in 367 patients with MDD and 161 healthy controls (HCs) to detect germline and copy number variations in the Korean population. Gene-based rare variants were analyzed to investigate the association between the genes and individuals, followed by neuroimaging-genetic analysis to explore the neural mechanisms underlying the genetic impact in 234 patients with MDD and 135 HCs using diffusion tensor imaging data. We identified 40 MDD-related genes and observed 95 recurrent regions of copy number variations. We also discovered a novel gene, FRMPD3, carrying rare variants that influence MDD. In addition, the single nucleotide polymorphism rs771995197 in the MUC6 gene was significantly associated with the integrity of widespread white matter tracts. Moreover, we identified 918 rare exonic missense variants in genes associated with MDD susceptibility. We postulate that rare variants of FRMPD3 may contribute significantly to MDD, with a mild penetration effect.

SNP-based and haplotype-based genome-wide association on drug dependence in Han Chinese

BACKGROUND

Drug addiction is a serious problem worldwide and is influenced by genetic factors. The present study aimed to investigate the association between genetics and drug addiction among Han Chinese.

METHODS

A total of 1000 Chinese users of illicit drugs and 9693 healthy controls were enrolled and underwent single nucleotide polymorphism (SNP)-based and haplotype-based association analyses via whole-genome genotyping.

RESULTS

Both single-SNP and haplotype tests revealed associations between illicit drug use and several immune-related genes in the major histocompatibility complex (MHC) region (SNP association: log10BF = 15.135, p = 1.054e-18; haplotype association: log10BF = 20.925, p = 2.065e-24). These genes may affect the risk of drug addiction via modulation of the neuroimmune system. The single-SNP test exclusively reported genome-wide significant associations between rs3782886 (SNP association: log10BF = 8.726, p = 4.842e-11) in BRAP and rs671 (SNP association: log10BF = 7.406, p = 9.333e-10) in ALDH2 and drug addiction. The haplotype test exclusively reported a genome-wide significant association (haplotype association: log10BF = 7.607, p = 3.342e-11) between a region with allelic heterogeneity on chromosome 22 and drug addiction, which may be involved in the pathway of vitamin B12 transport and metabolism, indicating a causal link between lower vitamin B12 levels and methamphetamine addiction.

CONCLUSIONS

These findings provide new insights into risk-modeling and the prevention and treatment of methamphetamine and heroin dependence, which may further contribute to potential novel therapeutic approaches.

Inter- and transgenerational heritability of preconception chronic stress or alcohol exposure: Translational outcomes in brain and behavior

Chronic stress and alcohol (ethanol) use are highly interrelated and can change an individual's behavior through molecular adaptations that do not change the DNA sequence, but instead change gene expression. A recent wealth of research has found that these nongenomic changes can be transmitted across generations, which could partially account for the "missing heritability" observed in genome-wide association studies of alcohol use disorder and other stress-related neuropsychiatric disorders. In this review, we summarize the molecular and behavioral outcomes of nongenomic inheritance of chronic stress and ethanol exposure and the germline mechanisms that could give rise to this heritability. In doing so, we outline the need for further research to: (1) Investigate individual germline mechanisms of paternal, maternal, and biparental nongenomic chronic stress- and ethanol-related inheritance; (2) Synthesize and dissect cross-generational chronic stress and ethanol exposure; (3) Determine cross-generational molecular outcomes of preconception ethanol exposure that contribute to alcohol-related disease risk, using cancer as an example. A detailed understanding of the cross-generational nongenomic effects of stress and/or ethanol will yield novel insight into the impact of ancestral perturbations on disease risk across generations and uncover actionable targets to improve human health.

Genetic and molecular architecture of complex traits

Human genetics has emerged as one of the most dynamic areas of biology, with a broadening societal impact. In this review, we discuss recent achievements, ongoing efforts, and future challenges in the field. Advances in technology, statistical methods, and the growing scale of research efforts have all provided many insights into the processes that have given rise to the current patterns of genetic variation. Vast maps of genetic associations with human traits and diseases have allowed characterization of their genetic architecture. Finally, studies of molecular and cellular effects of genetic variants have provided insights into biological processes underlying disease. Many outstanding questions remain, but the field is well poised for groundbreaking discoveries as it increases the use of genetic data to understand both the history of our species and its applications to improve human health.

Regulatory Variants on the Leukocyte Immunoglobulin-Like Receptor Gene Cluster are Associated with Crohn's Disease and Interact with Regulatory Variants for TAP2

BACKGROUND AND AIMS

Crohn's disease [CD] has a complex polygenic aetiology with high heritability. There is ongoing effort to identify novel variants associated with susceptibility to CD through a genome-wide association study [GWAS] in large Korean populations.

METHODS

Genome-wide variant data from 902 Korean patients with CD and 72 179 controls were used to assess the genetic associations in a meta-analysis with previous Korean GWAS results from 1621 patients with CD and 4419 controls. Epistatic interactions between CD-risk variants of interest were tested using a multivariate logistic regression model with an interaction term.

RESULTS

We identified two novel genetic associations with the risk of CD near ZBTB38 and within the leukocyte immunoglobulin-like receptor [LILR] gene cluster [p < 5 × 10-8], with highly consistent effect sizes between the two independent Korean cohorts. CD-risk variants in the LILR locus are known quantitative trait loci [QTL] for multiple LILR genes, of which LILRB2 directly interacts with various ligands including MHC class I molecules. The LILR lead variant exhibited a significant epistatic interaction with CD-associated regulatory variants for TAP2 involved in the antigen presentation of MHC class I molecules [p = 4.11 × 10-4], showing higher CD-risk effects of the TAP2 variant in individuals carrying more risk alleles of the LILR lead variant (odds ratio [OR] = 0.941, p = 0.686 in non-carriers; OR = 1.45, p = 2.51 × 10-4 in single-copy carriers; OR = 2.38, p = 2.76 × 10-6 in two-copy carriers).

CONCLUSIONS

This study demonstrated that genetic variants at two novel susceptibility loci and the epistatic interaction between variants in LILR and TAP2 loci confer a risk of CD.

Modeling genotype-protein interaction and correlation for Alzheimer's disease: a multi-omics imaging genetics study

Integrating and analyzing multiple omics data sets, including genomics, proteomics and radiomics, can significantly advance researchers' comprehensive understanding of Alzheimer's disease (AD). However, current methodologies primarily focus on the main effects of genetic variation and protein, overlooking non-additive effects such as genotype-protein interaction (GPI) and correlation patterns in brain imaging genetics studies. Importantly, these non-additive effects could contribute to intermediate imaging phenotypes, finally leading to disease occurrence. In general, the interaction between genetic variations and proteins, and their correlations are two distinct biological effects, and thus disentangling the two effects for heritable imaging phenotypes is of great interest and need. Unfortunately, this issue has been largely unexploited. In this paper, to fill this gap, we propose $\textbf{M}$ulti-$\textbf{T}$ask $\textbf{G}$enotype-$\textbf{P}$rotein $\textbf{I}$nteraction and $\textbf{C}$orrelation disentangling method ($\textbf{MT-GPIC}$) to identify GPI and extract correlation patterns between them. To ensure stability and interpretability, we use novel and off-the-shelf penalties to identify meaningful genetic risk factors, as well as exploit the interconnectedness of different brain regions. Additionally, since computing GPI poses a high computational burden, we develop a fast optimization strategy for solving MT-GPIC, which is guaranteed to converge. Experimental results on the Alzheimer's Disease Neuroimaging Initiative data set show that MT-GPIC achieves higher correlation coefficients and classification accuracy than state-of-the-art methods. Moreover, our approach could effectively identify interpretable phenotype-related GPI and correlation patterns in high-dimensional omics data sets. These findings not only enhance the diagnostic accuracy but also contribute valuable insights into the underlying pathogenic mechanisms of AD.

Genetic interactions of schizophrenia using gene-based statistical epistasis exclusively identify nervous system-related pathways and key hub genes

Background: The relationship between genotype and phenotype is governed by numerous genetic interactions (GIs), and the mapping of GI networks is of interest for two main reasons: 1) By modelling biological robustness, GIs provide a powerful opportunity to infer compensatory biological mechanisms via the identification of functional relationships between genes, which is of interest for biological discovery and translational research. Biological systems have evolved to compensate for genetic (i.e., variations and mutations) and environmental (i.e., drug efficacy) perturbations by exploiting compensatory relationships between genes, pathways and biological processes; 2) GI facilitates the identification of the direction (alleviating or aggravating interactions) and magnitude of epistatic interactions that influence the phenotypic outcome. The generation of GIs for human diseases is impossible using experimental biology approaches such as systematic deletion analysis. Moreover, the generation of disease-specific GIs has never been undertaken in humans. Methods: We used our Indian schizophrenia case-control (case-816, controls-900) genetic dataset to implement the workflow. Standard GWAS sample quality control procedure was followed. We used the imputed genetic data to increase the SNP coverage to analyse epistatic interactions across the genome comprehensively. Using the odds ratio (OR), we identified the GIs that increase or decrease the risk of a disease phenotype. The SNP-based epistatic results were transformed into gene-based epistatic results. Results: We have developed a novel approach by conducting gene-based statistical epistatic analysis using an Indian schizophrenia case-control genetic dataset and transforming these results to infer GIs that increase the risk of schizophrenia. There were ∼9.5 million GIs with a p-value ≤ 1 × 10-5. Approximately 4.8 million GIs showed an increased risk (OR > 1.0), while ∼4.75 million GIs had a decreased risk (OR <1.0) for schizophrenia. Conclusion: Unlike model organisms, this approach is specifically viable in humans due to the availability of abundant disease-specific genome-wide genotype datasets. The study exclusively identified brain/nervous system-related processes, affirming the findings. This computational approach fills a critical gap by generating practically non-existent heritable disease-specific human GIs from human genetic data. These novel datasets can train innovative deep-learning models, potentially surpassing the limitations of conventional GWAS.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LIPID TRANSFER PROTEIN 2 (LTP2) greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild-type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

Genome-wide exploration of genetic interactions for bladder cancer risk

Although GWASs have been conducted to investigate genetic variation of bladder tumorigenesis, little is known about genetic interactions that may influence bladder cancer (BC) risk. By leveraging large-scale participants from UK Biobank, we established a discovery database with 4000 Caucasian participants (2000 cases vs 2000 non-cases), a database with 1648 Caucasian participants (824 cases vs 824 non-cases) and 856 non-Caucasian participants (428 cases vs 428 non-cases) as validation. We then performed a genome-wide SNP-SNP interaction investigation related to BC risk based a machine learning approach (ie, GenEpi). Moreover, we used the selected interactions to build a BC screening model with an integrated interaction-empowered polygenic risk score (iPRS) based on Cox proportional hazard model. With Bonferroni correction, we identified 10 statistically significant pairs of SNPs, which located in 17 chromosomes. Of these, four SNP-SNP interactions were found to be positively associated with BC risk among Caucasian participants (ORs 1.57-2.03), while six SNP-SNP interactions showed negatively associated with BC risk (ORs 0.54-0.65). Only four of the SNP-SNP interactions were consistently identified in non-Caucasian participants located in ST7L-ADSS2, FHIT-CHDH, LARP4B-LHPP and RBFOX3-MPRIP. In addition, the iPRS showed a HR of 1.81 (95% CI: 1.46-2.09) compared the highest tertile to the lowest tertile, with an enhanced AUC (0.91; 95% CI:0.85-0.97) than PRS (AUC: 0.86; 95% CI:0.76-0.95; P-DeLong test  = 2.2 × 10-4 ). In summary, this study identified several important SNP-SNP interactions for BC risk, and developed an iPRS model for BC screening, which may help to identify the people at high-risk state of BC before early manifestation.

The Gene: An appraisal

The gene can be described as the foundational concept of modern biology. As such, it has spilled over into daily discourse, yet it is acknowledged among biologists to be ill-defined. Here, following a short history of the gene, I analyse critically its role in inheritance, evolution, development, and morphogenesis. Wilhelm Johannsen's genotype-conception, formulated in 1910, has been adopted as the foundation stone of genetics, giving the gene a higher degree of prominence than is justified by the evidence. An analysis of the results of the Long-Term Evolution Experiment (LTEE) with E. coli bacteria, grown over 60,000 generations, does not support spontaneous gene mutation as the source of variance for natural selection. From this it follows that the gene is not Mendel's unit of inheritance: that must be Johannsen's transmission-conception at the gamete phenotype level, a form of inheritance that Johannsen did not consider. Alternatively, I contend that biology viewed on the bases of thermodynamics, complex system dynamics, and self-organisation, provides a new framework for the foundations of biology. In this framework, the gene plays a passive role as a vital information store: it is the phenotype that plays the active role in inheritance, evolution, development, and morphogenesis.

A phenome-wide association study of methylated GC-rich repeats identifies a GCC repeat expansion in AFF3 as a significant cause of intellectual disability

GC-rich tandem repeat expansions (TREs) are often associated with DNA methylation, gene silencing and folate-sensitive fragile sites and underlie several congenital and late-onset disorders. Through a combination of DNA methylation profiling and tandem repeat genotyping, we identified 24 methylated TREs and investigated their effects on human traits using PheWAS in 168,641 individuals from the UK Biobank, identifying 156 significant TRE:trait associations involving 17 different TREs. Of these, a GCC expansion in the promoter of AFF3 was linked with a 2.4-fold reduced probability of completing secondary education, an effect size comparable to several recurrent pathogenic microdeletions. In a cohort of 6,371 probands with neurodevelopmental problems of suspected genetic etiology, we observed a significant enrichment of AFF3 expansions compared to controls. With a population prevalence that is at least 5-fold higher than the TRE that causes fragile X syndrome, AFF3 expansions represent a significant cause of neurodevelopmental delay.

A multidimensional social risk atlas of depression and anxiety: An observational and genome-wide environmental interaction study

Background

Mental disorders are largely socially determined, yet the combined impact of multidimensional social factors on the two most common mental disorders, depression and anxiety, remains unclear.

Methods

We constructed a polysocial risk score (PsRS), a multidimensional social risk indicator including components from three domains: socioeconomic status, neighborhood and living environment and psychosocial factors. Supported by the UK Biobank cohort, we randomly divided 110 332 participants into the discovery cohort (60%; n = 66 200) and the replication cohort (40%; n = 44 134). We tested the associations between 13 single social factors with Patient Health Questionnaire (PHQ) score, Generalized Anxiety Disorder Scale (GAD) score and self-reported depression and anxiety. The significant social factors were used to calculate PsRS for each mental disorder by considering weights from the multivariable linear model. Generalized linear models were applied to explore the association between PsRS and depression and anxiety. Genome-wide environmental interaction study (GWEIS) was further performed to test the effect of interactions between PsRS and SNPs on the risk of mental phenotypes.

Results

In the discovery cohort, PsRS was positively associated with PHQ score (β = 0.37; 95% CI = 0.35-0.38), GAD score (β = 0.27; 95% CI = 0.25-0.28), risk of self-reported depression (OR = 1.29; 95% CI = 1.28-1.31) and anxiety (OR = 1.19; 95% CI = 1.19-1.23). Similar results were observed in the replication cohort. Emotional stress, lack of social support and low household income were significantly associated with the development of depression and anxiety. GWEIS identified multiple candidate loci for PHQ score, such as rs149137169 (ST18) (Pdiscovery = 1.08 × 10-8, Preplication = 3.25 × 10-6) and rs3759812 (MYO9A) (Pdiscovery = 3.87 × 10-9, Preplication = 6.21 × 10-5). Additionally, seven loci were detected for GAD score, such as rs114006170 (TMPRSS11D) (Pdiscovery = 1.14 × 10-9, Preplication = 7.36 × 10-5) and rs77927903 (PIP4K2A) (Pdiscovery = 2.40 × 10-9, Preplication = 0.002).

Conclusions

Our findings reveal the positive effects of multidimensional social factors on the risk of depression and anxiety. It is important to address key social disadvantage in mental health promotion and treatment.

Genome-wide association studies for epistatic genetic effects on fertility and reproduction traits in Holstein cattle

Non-additive genetic effects are well known to play an important role in the phenotypic expression of complex traits, such as fertility and reproduction. In this study, a genome scan was performed using 41,640 single nucleotide polymorphism (SNP) markers to identify genomic regions associated with epistatic (additive-by-additive) effects in fertility and reproduction traits in Holstein cattle. Nine fertility and reproduction traits were analysed on 5825 and 6090 Holstein heifers and cows with phenotypes and genotypes, respectively. The Marginal Epistasis Test (MAPIT) was used to identify SNPs with significant marginal epistatic effects at a chromosome-wise 5% and 10% false discovery rate (FDR) level. The -log10 (p) values were adjusted by the genomic inflation factor (λ) to correct for the potential bias on the p-values and minimize the possible effects of population stratification. After adjustments, MAPIT enabled the identification of genomic regions with significant marginal epistatic effects for heifers on BTA5 for age at first insemination, BTA3 and BTA24 for non-return rate (NRR); BTA16 and BTA28 for gestation length (GL); BTA1, BTA4 and BTA17 for stillbirth (SB). For the cow traits, MAPIT enabled the identification of regions on BTA11 for GL, BTA11 and BTA16 for SB and BTA19 for calf size (CZ). An additional approach for mapping epistasis in a genome-wide association study was also proposed, in which the genome scan was performed using estimates of epistatic values as the input pseudo-phenotypes, computed using single-trait animal models. Significant SNPs were identified at the chromosome-wise 5% and 10% FDR levels for all traits. For the heifer traits, significant regions were found on BTA7 for AFS; BTA12 for NRR; BTA14 and BTA19 for GL; BTA19 for calving ease (CE); BTA5, BTA24, BTA25 and in the X chromosome for SB; BTA23 and in the X chromosome for CZ and in the X chromosome for the number of services (NS). For the cow traits, significant regions were found on BTA29 and in the X chromosome for NRR, BTA11, BTA16 and in the X chromosome for SB, BTA2 for GL, BTA28 for CZ, BTA19 for calving to first insemination, and in the X chromosome for NS and first insemination to conception. The results suggest that the epistatic genetic effects are likely due to many loci with a small effect rather than few loci with a large effect and/or a single SNP marker alone do not capture the epistatic effects well. The genomic architecture of fertility and reproduction traits is complex, and these results should be validated in independent dairy cattle populations and using alternative statistical models.

Genome wide copy number variations using Porcine 60K SNP Beadchip in Landlly pigs

In the present study, Porcine 60K SNP genotype data from 69 Landlly pigs were used to explore Copy Number Variations (CNVs) across the autosomes. A total of 386 CNVs were identified using Hidden Markov Model (HMM) in PennCNV software, which were subsequently aggregated to 115 CNV regions (CNVRs). Among the total detected CNVRs, 58 gain, 49 were loss type while remaining 8 events were both gain and loss types. Identified CNVRs covered 12.5 Mb (0.55%) of Sus scrofa reference 11.1 genome. Comparison of our results with previous investigations on pigs revealed that approximately 75% CNVRs were novel, which may be due to differences in genetic background, environment and implementation of artificial selection in Landlly pigs. Functional annotation and pathway analysis showed the significant enrichment of 267 well-annotated Sus scrofa genes in CNVRs. These genes were involved in different biological functions like sensory perception, meat quality traits, back fat thickness and immunity. Additionally, KIT and FUT1 were two major genes detected on CNVR in our population. This investigation provided a comprehensive overview of CNV distribution in the Indian porcine genome for the first time, which may be useful for further investigating the association of important quantitative traits in Landlly pigs.Highlights115 CNVRs were identified in 69 Landlly pig population.Approximately 75% detected CNVRs were novel for Landlly population.Significant enrichment of 267 well-annotated Sus scrofa genes observed in these CNVRs.These genes were involved in different biological functions like sensory perception, meat quality traits, back fat thickness and immunity.Comprehensive CNV map in the Indian porcine genome developed for the first time.

Precision medicine of obesity as an integral part of type 2 diabetes management - past, present, and future

Obesity is a complex and heterogeneous condition that leads to various metabolic complications, including type 2 diabetes. Unfortunately, for some, treatment options to date for obesity are insufficient, with many people not reaching sustained weight loss or having improvements in metabolic health. In this Review, we discuss advances in the genetics of obesity from the past decade-with emphasis on developments from the past 5 years-with a focus on metabolic consequences, and their potential implications for precision management of the disease. We also provide an overview of the potential role of genetics in guiding weight loss strategies. Finally, we propose a vision for the future of precision obesity management that includes developing an obesity-centred multidisease management algorithm that targets both obesity and its comorbidities. However, further collaborative efforts and research are necessary to fully realise its potential and improve metabolic health outcomes.

The pan-genome and local adaptation of Arabidopsis thaliana

Arabidopsis thaliana serves as a model species for investigating various aspects of plant biology. However, the contribution of genomic structural variations (SVs) and their associate genes to the local adaptation of this widely distribute species remains unclear. Here, we de novo assemble chromosome-level genomes of 32 A. thaliana ecotypes and determine that variable genes expand the gene pool in different ecotypes and thus assist local adaptation. We develop a graph-based pan-genome and identify 61,332 SVs that overlap with 18,883 genes, some of which are highly involved in ecological adaptation of this species. For instance, we observe a specific 332 bp insertion in the promoter region of the HPCA1 gene in the Tibet-0 ecotype that enhances gene expression, thereby promotes adaptation to alpine environments. These findings augment our understanding of the molecular mechanisms underlying the local adaptation of A. thaliana across diverse habitats.

Hematologic setpoints are a stable and patient-specific deep phenotype

The complete blood count is an important screening tool for healthy adults and is the most commonly ordered test at periodic physical exams. However, results are usually interpreted relative to one-size-fits-all reference intervals, undermining the goal of precision medicine to tailor medical care to the needs of individual patients based on their unique characteristics. Here we show that standard complete blood count indices in healthy adults have robust homeostatic setpoints that are patient-specific and stable, with the typical healthy adult's set of 9 blood count setpoints distinguishable from 98% of others, and with these differences persisting for decades. These setpoints reflect a deep physiologic phenotype, enabling improved detection of both acquired and genetic determinants of hematologic regulation, including discovery of multiple novel loci via GWAS analyses. Patient-specific reference intervals derived from setpoints enable more accurate personalized risk assessment, and the setpoints themselves are significantly correlated with mortality risk, providing new opportunities to enhance patient-specific screening and early intervention. This study shows complete blood count setpoints are sufficiently stable and patient-specific to help realize the promise of precision medicine for healthy adults.

Association detection between multiple traits and rare variants based on family data via a nonparametric method

Background

The rapid development of next-generation sequencing technologies allow people to analyze human complex diseases at the molecular level. It has been shown that rare variants play important roles for human diseases besides common variants. Thus, effective statistical methods need to be proposed to test for the associations between traits (e.g., diseases) and rare variants. Currently, more and more rare genetic variants are being detected throughout the human genome, which demonstrates the possibility to study rare variants. Yet complex diseases are usually measured as a variety of forms, such as binary, ordinal, quantitative, or some mixture of them. Therefore, the genetic mapping problem can be attributable to the association detection between multiple traits and multiple loci, with sufficiently considering the correlated structure among multiple traits.

Methods

In this article, we construct a new non-parametric statistic by the generalized Kendall's τ theory based on family data. The new test statistic has an asymptotic distribution, it can be used to study the associations between multiple traits and rare variants, which broadens the way to identify genetic factors of human complex diseases.

Results

We apply our method (called Nonp-FAM) to analyze simulated data and GAW17 data, and conduct comprehensive comparison with some existing methods. Experimental results show that the proposed family-based method is powerful and robust for testing associations between multiple traits and rare variants, even if the data has some population stratification effect.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LTP2 greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs, and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

Genetic risk assessment based on association and prediction studies

The genetic basis of phenotypic emergence provides valuable information for assessing individual risk. While association studies have been pivotal in identifying genetic risk factors within a population, complementing it with insights derived from predictions studies that assess individual-level risk offers a more comprehensive approach to understanding phenotypic expression. In this study, we established personalized risk assessment models using single-nucleotide polymorphism (SNP) data from 200 Korean patients, of which 100 experienced hepatitis B surface antigen (HBsAg) seroclearance and 100 patients demonstrated high levels of HBsAg. The risk assessment models determined the predictive power of the following: (1) genome-wide association study (GWAS)-identified candidate biomarkers considered significant in a reference study and (2) machine learning (ML)-identified candidate biomarkers with the highest feature importance scores obtained by using random forest (RF). While utilizing all features yielded 64% model accuracy, using relevant biomarkers achieved higher model accuracies: 82% for 52 GWAS-identified candidate biomarkers, 71% for three GWAS-identified biomarkers, and 80% for 150 ML-identified candidate biomarkers. Findings highlight that the joint contributions of relevant biomarkers significantly influence phenotypic emergence. On the other hand, combining ML-identified candidate biomarkers into the pool of GWAS-identified candidate biomarkers resulted in the improved predictive accuracy of 90%, demonstrating the capability of ML as an auxiliary analysis to GWAS. Furthermore, some of the ML-identified candidate biomarkers were found to be linked with hepatocellular carcinoma (HCC), reinforcing previous claims that HCC can still occur despite the absence of HBsAg.

A biobank-scale test of marginal epistasis reveals genome-wide signals of polygenic epistasis

The contribution of epistasis (interactions among genes or genetic variants) to human complex trait variation remains poorly understood. Methods that aim to explicitly identify pairs of genetic variants, usually single nucleotide polymorphisms (SNPs), associated with a trait suffer from low power due to the large number of hypotheses tested while also having to deal with the computational problem of searching over a potentially large number of candidate pairs. An alternate approach involves testing whether a single SNP modulates variation in a trait against a polygenic background. While overcoming the limitation of low power, such tests of polygenic or marginal epistasis (ME) are infeasible on Biobank-scale data where hundreds of thousands of individuals are genotyped over millions of SNPs. We present a method to test for ME of a SNP on a trait that is applicable to biobank-scale data. We performed extensive simulations to show that our method provides calibrated tests of ME. We applied our method to test for ME at SNPs that are associated with 53 quantitative traits across ≈ 300 K unrelated white British individuals in the UK Biobank (UKBB). Testing 15, 601 trait-loci associations that were significant in GWAS, we identified 16 trait-loci pairs across 12 traits that demonstrate strong evidence of ME signals (p-value p < 5 × 10 - 8 53 ). We further partitioned the significant ME signals across the genome to identify 6 trait-loci pairs with evidence of local (within-chromosome) ME while 15 show evidence of distal (cross-chromosome) ME. Across the 16 trait-loci pairs, we document that the proportion of trait variance explained by ME is about 12x as large as that explained by the GWAS effects on average (range: 0.59 to 43.89). Our results show, for the first time, evidence of interaction effects between individual genetic variants and overall polygenic background modulating complex trait variation.

Genome-wide association studies: an intuitive solution for SNP identification and gene mapping in trees

Analysis of natural diversity in wild/cultivated plants can be used to understand the genetic basis for plant breeding programs. Recent advancements in DNA sequencing have expanded the possibilities for genetically altering essential features. There have been several recently disclosed statistical genetic methods for discovering the genes impacting target qualities. One of these useful methods is the genome-wide association study (GWAS), which effectively identifies candidate genes for a variety of plant properties by examining the relationship between a molecular marker (such as SNP) and a target trait. Conventional QTL mapping with highly structured populations has major limitations. The limited number of recombination events results in poor resolution for quantitative traits. Only two alleles at any given locus can be studied simultaneously. Conventional mapping approach fails to work in perennial plants and vegetatively propagated crops. These limitations are sidestepped by association mapping or GWAS. The flexibility of GWAS comes from the fact that the individuals being examined need not be linked to one another, allowing for the use of all meiotic and recombination events to increase resolution. Phenotyping, genotyping, population structure analysis, kinship analysis, and marker-trait association analysis are the fundamental phases of GWAS. With the rapid development of sequencing technologies and computational methods, GWAS is becoming a potent tool for identifying the natural variations that underlie complex characteristics in crops. The use of high-throughput sequencing technologies along with genotyping approaches like genotyping-by-sequencing (GBS) and restriction site associated DNA (RAD) sequencing may be highly useful in fast-forward mapping approach like GWAS. Breeders may use GWAS to quickly unravel the genomes through QTL and association mapping by taking advantage of natural variances. The drawbacks of conventional linkage mapping can be successfully overcome with the use of high-resolution mapping and the inclusion of multiple alleles in GWAS.

Family studies: A useful tool to better understand spondyloarthritis

Spondyloarthritis (SpA) is an immune-mediated disease characterized by a high heritability, reflected by strong familial aggregation. Therefore, family studies are a powerful tool for elucidating the genetic basis of SpA. First, they helped to assess the relative importance of genetic and environmental factors and established the polygenic character of the disease. Family-based designs were also historically used to identify genetic factors of susceptibility through linkage analyses. In SpA, three whole-genome linkage studies were published in the 1990's, unfortunately with few consistent results. After having been put aside for several years in favour of case-control GWAS, there is a renewed interest in family-based designs in particular to detect rare variant associations. This review aims at summarizing what family studies have brought to the field of SpA genetics, from genetic epidemiology studies to the most recent rare variant analyses. It also highlights the potential interest of family history of SpA to help diagnosis and detection of patients at high risk to develop the disease.

Complex traits and candidate genes: estimation of genetic variance components across multiple genetic architectures

Large-effect loci-those statistically significant loci discovered by genome-wide association studies or linkage mapping-associated with key traits segregate amidst a background of minor, often undetectable, genetic effects in wild and domesticated plants and animals. Accurately attributing mean differences and variance explained to the correct components in the linear mixed model analysis is vital for selecting superior progeny and parents in plant and animal breeding, gene therapy, and medical genetics in humans. Marker-assisted prediction and its successor, genomic prediction, have many advantages for selecting superior individuals and understanding disease risk. However, these two approaches are less often integrated to study complex traits with different genetic architectures. This simulation study demonstrates that the average semivariance can be applied to models incorporating Mendelian, oligogenic, and polygenic terms simultaneously and yields accurate estimates of the variance explained for all relevant variables. Our previous research focused on large-effect loci and polygenic variance separately. This work aims to synthesize and expand the average semivariance framework to various genetic architectures and the corresponding mixed models. This framework independently accounts for the effects of large-effect loci and the polygenic genetic background and is universally applicable to genetics studies in humans, plants, animals, and microbes.

Gene Interactions in Human Disease Studies-Evidence Is Mounting

Despite monumental advances in molecular technology to generate genome sequence data at scale, there is still a considerable proportion of heritability in most complex diseases that remains unexplained. Because many of the discoveries have been single-nucleotide variants with small to moderate effects on disease, the functional implication of many of the variants is still unknown and, thus, we have limited new drug targets and therapeutics. We, and many others, posit that one primary factor that has limited our ability to identify novel drug targets from genome-wide association studies may be due to gene interactions (epistasis), gene-environment interactions, network/pathway effects, or multiomic relationships. We propose that many of these complex models explain much of the underlying genetic architecture of complex disease. In this review, we discuss the evidence from multiple research avenues, ranging from pairs of alleles to multiomic integration studies and pharmacogenomics, that supports the need for further investigation of gene interactions (or epistasis) in genetic and genomic studies of human disease. Our goal is to catalog the mounting evidence for epistasis in genetic studies and the connections between genetic interactions and human health and disease that could enable precision medicine of the future.