Rare and common variants: twenty arguments

The Contribution of Mosaic Chromosomal Alterations to Schizophrenia

BACKGROUND

Mosaic chromosomal alterations are implicated in neuropsychiatric disorders, but the contribution to schizophrenia (SCZ) risk for somatic copy number variations (sCNVs) emerging in early developmental stages has not been fully established.

METHODS

We analyzed blood-derived genotype arrays from 9715 patients with SCZ and 28,822 control participants of Chinese descent using a computational tool (MoChA) based on long-range chromosomal information to detect mosaic chromosomal alterations. We focused on probable early developmental sCNVs through stringent filtering. We assessed the burden of sCNVs across varying cell fraction cutoffs, as well as the frequency with which genes were involved in sCNVs. We integrated this data with the PGC (Psychiatric Genomics Consortium) dataset, which comprises 12,834 SCZ cases and 11,648 controls of European descent, and complemented it with genotyping data from postmortem brain tissue of 936 participants (449 cases and 487 controls).

RESULTS

Patients with SCZ had a significantly higher somatic losses detection rate than control participants (1.00% vs. 0.52%; odds ratio = 1.91; 95% CI, 1.47-2.49; two-sided Fisher's exact test, p = 1.49 × 10-6). Further analysis indicated that the odds ratios escalated proportionately (from 1.91 to 2.78) with the increment in cell fraction cutoffs. Recurrent sCNVs associated with SCZ (odds ratio > 8; Fisher's exact test, p < .05) were identified, including notable regions at 10q21.1 (ZWINT), 3q26.1 (SLITRK3), 1q31.1 (BRINP3) and 12q21.31-21.32 (MGAT4C and NTS) in the Chinese cohort, and some regions were validated with PGC data. Cross-tissue validation pinpointed somatic losses at loci like 1p35.3-35.2 and 19p13.3-13.2.

CONCLUSIONS

The study highlights the significant impact of mosaic chromosomal alterations on SCZ, suggesting their pivotal role in the disorder's genetic etiology.

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.

SEAD reference panel with 22,134 haplotypes boosts rare variant imputation and genome-wide association analysis in Asian populations

Limited whole genome sequencing (WGS) studies in Asian populations result in a lack of representative reference panels, thus hindering the discovery of ancestry-specific variants. Here, we present the South and East Asian reference Database (SEAD) panel ( https://imputationserver.westlake.edu.cn/ ), which integrates WGS data for 11,067 individuals from various sources across 17 Asian countries. The SEAD panel, comprising 22,134 haplotypes and 88,294,957 variants, demonstrates improved imputation accuracy for South Asian populations compared to 1000 Genomes Project, TOPMed, and ChinaMAP panels, with a higher proportion of well-imputed rare variants. For East Asian populations, SEAD shows concordance comparable to ChinaMAP, but outperforming TOPMed. Additionally, we apply the SEAD panel to conduct a genome-wide association study for total hip (Hip) and femoral neck (FN) bone mineral density (BMD) traits in 5369 genotyped Chinese samples. The single-variant test suggests that rare variants near SNTG1 are associated with Hip BMD (rs60103302, MAF = 0.0092, P = 1.67 × 10-7), and variant-set analysis further supports the association (Pslide_window = 9.08 × 10-9, Pgene_centric = 5.27 × 10-8). This association was not reported previously and can only be detected by using Asian reference panels. Preliminary in vitro experiments for one of the rare variants identified provide evidence that it upregulates SNTG1 expression, which could in turn inhibit the proliferation and differentiation of preosteoblasts.

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.

Mendelian randomization in Alzheimer's disease and mild cognitive impairment: Hippocampal volume associations

This study investigates the association between cognitive dysfunction and hippocampal volumes in Alzheimer's Disease (AD) and Mild Cognitive Impairment (MCI) using Mendelian randomization. A meta-analysis of 503 healthy controls, 562 MCI patients, and 389 CE patients revealed significant reductions in hippocampal and subregion volumes in MCI and AD compared to controls. While various subregions showed volume reductions, no causal relationship between hippocampal volume and AD was established through Mendelian randomization analysis. In conclusion, significant volume reductions were observed in MCI and AD patients, highlighting the complexity of the relationship between hippocampal volume and cognitive impairment.

Polygenic risk of social isolation behavior and its influence on psychopathology and personality

Social isolation has been linked to a range of psychiatric issues, but the behavioral component that drives it is not well understood. Here, a genome-wide associations study (GWAS) was carried out to identify genetic variants that contribute specifically to social isolation behavior (SIB) in up to 449,609 participants from the UK Biobank. 17 loci were identified at genome-wide significance, contributing to a 4% SNP-based heritability estimate. Using the SIB GWAS, polygenic risk scores (PRS) were derived in ALSPAC, an independent, developmental cohort, and used to test for association with self-reported friendship scores, comprising items related to friendship quality and quantity, at age 12 and 18 to determine whether genetic predisposition manifests during childhood development. At age 18, friendship scores were associated with the SIB PRS, demonstrating that the genetic factors can predict related social traits in late adolescence. Linkage disequilibrium (LD) score correlation using the SIB GWAS demonstrated genetic correlations with autism spectrum disorder (ASD), schizophrenia, major depressive disorder (MDD), educational attainment, extraversion, and loneliness. However, no evidence of causality was found using a conservative Mendelian randomization approach between SIB and any of the traits in either direction. Genomic Structural Equation Modeling (SEM) revealed a common factor contributing to SIB, neuroticism, loneliness, MDD, and ASD, weakly correlated with a second common factor that contributes to psychiatric and psychotic traits. Our results show that SIB contributes a small heritable component, which is associated genetically with other social traits such as friendship as well as psychiatric disorders.

Modeling common and rare genetic risk factors of neuropsychiatric disorders in human induced pluripotent stem cells

Recent genome-wide association studies (GWAS) and whole-exome sequencing of neuropsychiatric disorders, especially schizophrenia, have identified a plethora of common and rare disease risk variants/genes. Translating the mounting human genetic discoveries into novel disease biology and more tailored clinical treatments is tied to our ability to causally connect genetic risk variants to molecular and cellular phenotypes. When combined with the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/CRISPR-associated (Cas) nuclease-mediated genome editing system, human induced pluripotent stem cell (hiPSC)-derived neural cultures (both 2D and 3D organoids) provide a promising tractable cellular model for bridging the gap between genetic findings and disease biology. In this review, we first conceptualize the advances in understanding the disease polygenicity and convergence from the past decade of iPSC modeling of different types of genetic risk factors of neuropsychiatric disorders. We then discuss the major cell types and cellular phenotypes that are most relevant to neuropsychiatric disorders in iPSC modeling. Finally, we critically review the limitations of iPSC modeling of neuropsychiatric disorders and outline the need for implementing and developing novel methods to scale up the number of iPSC lines and disease risk variants in a systematic manner. Sufficiently scaled-up iPSC modeling and a better functional interpretation of genetic risk variants, in combination with cutting-edge CRISPR/Cas9 gene editing and single-cell multi-omics methods, will enable the field to identify the specific and convergent molecular and cellular phenotypes in precision for neuropsychiatric disorders.

SNP Genotype Imputation in Forensics-A Performance Study

BACKGROUND/OBJECTIVES

Emerging forensic genetic applications, such as forensic investigative genetic genealogy (FIGG), advanced DNA phenotyping, and distant kinship inference, increasingly require dense SNP genotype datasets. However, forensic-grade DNA often contains missing genotypes due to its quality and quantity limitations, potentially hindering these applications. Genotype imputation, a method that predicts missing genotypes, is widely used in population and medical genetics, but its utility in forensic genetics has not been thoroughly explored. This study aims to assess the performance of genotype imputation in forensic contexts and determine the conditions under which it can be effectively applied.

METHODS

We employed a simulation-based approach to generate realistic forensic SNP genotype datasets with varying numbers, densities, and qualities of observed genotypes. Genotype imputation was performed using Beagle software, and the performance was evaluated based on the call rate and imputation accuracy across different datasets and imputation settings.

RESULTS

The results demonstrate that genotype imputation can significantly increase the number of SNP genotypes. However, imputation accuracy was dependent on factors such as the quality of the original genotype data and the characteristics of the reference population. Higher SNP density and fewer genotype errors generally resulted in improved imputation accuracy.

CONCLUSIONS

This study highlights the potential of genotype imputation to enhance forensic SNP datasets but underscores the importance of optimizing imputation parameters and understanding the limitations of the original data. These findings will inform the future application of imputation in forensic genetics, supporting its integration into forensic workflows.

Genetic Polymorphisms Associated with Fetal Hemoglobin (HbF) Levels and F-Cell Numbers: A Systematic Review of Genome-Wide Association Studies

Elevated fetal hemoglobin (HbF), which is partly controlled by genetic modifiers, ameliorates disease severity in β hemoglobinopathies. Understanding the genetic basis of this trait holds great promise for personalized therapeutic approaches. PubMed, MedRxiv, and the GWAS Catalog were searched up to May 2024 to identify eligible GWAS studies following PRISMA guidelines. Four independent reviewers screened, extracted, and synthesized data using narrative and descriptive methods. Study quality was assessed using a modified version of the Q-Genie tool. Pathway enrichment analysis was conducted on gene lists derived from the selected GWAS studies. Out of 113 initially screened studies, 62 underwent full-text review, and 16 met the inclusion criteria for quality assessment and data synthesis. A total of 939 significant SNP-trait associations (p-value < 1 × 10-5) were identified, mapping to 133 genes (23 with overlapping variant positions) and 103 intergenic sequences. Most SNP-trait associations converged around BCL11A (chr.2), HBS1L-MYB, (chr.6), olfactory receptor and beta globin (HBB) gene clusters (chr.11), with less frequent loci including FHIT (chr.3), ALDH8A1, BACH2, RPS6KA2, SGK1 (chr.6), JAZF1 (chr.7), MMP26 (chr.11), COCH (chr.14), ABCC1 (chr.16), CTC1, PFAS (chr.17), GCDH, KLF1, NFIX, and ZBTB7A (chr.19). Pathway analysis highlighted Gene Ontology (GO) terms and pathways related to olfaction, hemoglobin and haptoglobin binding, and oxygen carrier activity. This systematic review confirms established genetic modifiers of HbF level, while highlighting less frequently associated loci as promising areas for further research. Expanding research across ethnic populations is essential for advancing personalized therapies and enhancing outcomes for individuals with sickle cell disease or β-thalassemia.

PathVar: A Customisable NGS Variant Calling Algorithm Implicates Novel Candidate Genes and Pathways in Hemiplegic Migraine

The exponential growth of next-generation sequencing (NGS) data requires innovative bioinformatics approaches to unravel the genetic underpinnings of diseases. Hemiplegic migraine (HM), a debilitating neurological disorder with a genetic basis, is one such condition that warrants further investigation. Notably, the genetic heterogeneity of HM is underscored by the fact that approximately two-thirds of patients lack pathogenic variants in the known causal ion channel genes. In this context, we have developed PathVar, a novel bioinformatics algorithm that harnesses publicly available tools and software for pathogenic variant discovery in NGS data. PathVar integrates a suite of tools, including HaplotypeCaller from the Genome Analysis Toolkit (GATK) for variant calling, Variant Effect Predictor (VEP) and ANNOVAR for variant annotation, and TAPES for assigning the American College of Medical Genetics and Genomics (ACMG) pathogenicity labels. Applying PathVar to whole exome sequencing data from 184 HM patients, we detected 648 variants that are probably pathogenic in multiple patients. Moreover, we have identified several candidate genes for HM, many of which cluster around the Rho GTPases pathway. Future research can leverage PathVar to generate high quality, candidate pathogenic variants, which may enhance our understanding of HM and other complex diseases.

The performance of AlphaMissense to identify genes influencing disease

A novel algorithm, AlphaMissense, has been shown to have an improved ability to predict the pathogenicity of rare missense genetic variants. However, it is not known whether AlphaMissense improves the ability of gene-based testing to identify disease-influencing genes. Using whole-exome sequencing data from the UK Biobank, we compared gene-based association analysis strategies including sets of deleterious variants: predicted loss-of-function (pLoF) variants only, pLoF plus AlphaMissense pathogenic variants, pLoF with missense variants predicted to be deleterious by any of five commonly utilized annotation methods (Missense (1/5)) or only variants predicted to be deleterious by all five methods (Missense (5/5)). We measured performance to identify 519 previously identified positive control genes, which can lead to Mendelian diseases, or are the targets of successfully developed medicines. These strategies identified 0.85 million pLoF variants and 5 million deleterious missense variants, including 22,131 likely pathogenic missense variants identified exclusively by AlphaMissense. The gene-based association tests found 608 significant gene associations (at p < 1.25 × 10-7) across 24 common traits and diseases. Compared with pLoFs plus Missense (5/5), tests using pLoFs and AlphaMissense variants found slightly more significant gene-disease and gene-trait associations, albeit with a marginally lower proportion of positive control genes. Nevertheless, their overall performance was similar. Merging AlphaMissense with Missense (5/5), whether through their intersection or union, did not yield any further enhancement in performance. In summary, employing AlphaMissense to select deleterious variants for gene-based testing did not improve the ability to identify genes that are known to influence disease.

Functional Roles of Sensorimotor Alpha and Beta Oscillations in Overt Speech Production

Power decreases, or desynchronization, of sensorimotor alpha and beta oscillations (i.e., alpha and beta ERD) have long been considered as indices of sensorimotor control in overt speech production. However, their specific functional roles are not well understood. Hence, we first conducted a systematic review to investigate how these two oscillations are modulated by speech motor tasks in typically fluent speakers (TFS) and in persons who stutter (PWS). Eleven EEG/MEG papers with source localization were included in our systematic review. The results revealed consistent alpha and beta ERD in the sensorimotor cortex of TFS and PWS. Furthermore, the results suggested that sensorimotor alpha and beta ERD may be functionally dissociable, with alpha related to (somato-)sensory feedback processing during articulation and beta related to motor processes throughout planning and articulation. To (partly) test this hypothesis of a potential functional dissociation between alpha and beta ERD, we then analyzed existing intracranial electro-encephalography (iEEG) data from the primary somatosensory cortex (S1) of picture naming. We found moderate evidence for alpha, but not beta, ERD's sensitivity to speech movements in S1, lending supporting evidence for the functional dissociation hypothesis identified by the systematic review.

Divergent Flowering Time Responses to Increasing Temperatures Are Associated With Transcriptome Plasticity and Epigenetic Modification Differences at FLC Promoter Region of Arabidopsis thaliana

Understanding the genetic, and transcriptomic changes that drive the phenotypic plasticity of fitness traits is a central question in evolutionary biology. In this study, we utilised 152 natural Swedish Arabidopsis thaliana accessions with re-sequenced genomes, transcriptomes and methylomes and measured flowering times (FTs) under two temperature conditions (10°C and 16°C) to address this question. We revealed that the northern accessions exhibited advanced flowering in response to decreased temperature, whereas the southern accessions delayed their flowering, indicating a divergent flowering response. This contrast in flowering responses was associated with the isothermality of their native ranges, which potentially enables the northern accessions to complete their life cycle more rapidly in years with shorter growth seasons. At the transcriptome level, we observed extensive rewiring of gene co-expression networks, with the expression of 25 core genes being associated with the mean FT and its plastic variation. Notably, variations in FLC expression sensitivity between northern and southern accessions were found to be associated with the divergence FT response. Further analysis suggests that FLC expression sensitivity is associated with differences in CG, CHG and CHH methylation at the promoter region. Overall, our study revealed the association between transcriptome plasticity and flowering time plasticity among different accessions, providing evidence for its relevance in ecological adaptation. These findings offer deeper insights into the genetics of rapid responses to environmental changes and ecological adaptation.

Translation of genome-wide association study: from genomic signals to biological insights

Since the turn of the 21st century, genome-wide association study (GWAS) have successfully identified genetic signals associated with a myriad of common complex traits and diseases. As we transition from establishing robust genetic associations with diverse phenotypes, the central challenge is now focused on characterizing the underlying functional mechanisms driving these signals. Previous GWAS efforts have revealed multiple variants, each conferring relatively subtle susceptibility, collectively contributing to the pathogenesis of various common diseases. Such variants can further exhibit associations with multiple other traits and differ across ancestries, plus disentangling causal variants from non-causal due to linkage disequilibrium complexities can lead to challenges in drawing direct biological conclusions. Combined with cellular context considerations, such challenges can reduce the capacity to definitively elucidate the biological significance of GWAS signals, limiting the potential to define mechanistic insights. This review will detail current and anticipated approaches for functional interpretation of GWAS signals, both in terms of characterizing the underlying causal variants and the corresponding effector genes.

Network medicine-based epistasis detection in complex diseases: ready for quantum computing

Most heritable diseases are polygenic. To comprehend the underlying genetic architecture, it is crucial to discover the clinically relevant epistatic interactions (EIs) between genomic single nucleotide polymorphisms (SNPs) (1-3). Existing statistical computational methods for EI detection are mostly limited to pairs of SNPs due to the combinatorial explosion of higher-order EIs. With NeEDL (network-based epistasis detection via local search), we leverage network medicine to inform the selection of EIs that are an order of magnitude more statistically significant compared to existing tools and consist, on average, of five SNPs. We further show that this computationally demanding task can be substantially accelerated once quantum computing hardware becomes available. We apply NeEDL to eight different diseases and discover genes (affected by EIs of SNPs) that are partly known to affect the disease, additionally, these results are reproducible across independent cohorts. EIs for these eight diseases can be interactively explored in the Epistasis Disease Atlas (https://epistasis-disease-atlas.com). In summary, NeEDL demonstrates the potential of seamlessly integrated quantum computing techniques to accelerate biomedical research. Our network medicine approach detects higher-order EIs with unprecedented statistical and biological evidence, yielding unique insights into polygenic diseases and providing a basis for the development of improved risk scores and combination therapies.

Microevolutionary hypothesis of the obesity epidemic

The obesity epidemic represents potentially the largest phenotypic change in Homo sapiens since the origin of the species. Despite obesity's high heritability, it is generally presumed a change in the gene pool could not have caused the obesity epidemic. Here we advance the hypothesis that a rapid change in the obesogenic gene pool has occurred second to the introduction of modern obstetrics dramatically altering evolutionary pressures on obesity-the microevolutionary hypothesis of the obesity epidemic. Obesity is known to increase childbirth-related mortality several fold. Prior to modern obstetrics, childbirth-related mortality occurred in over 10% of women in their lifetime. After modern obstetrics, this mortality reduced to a fraction of a percent, thereby lifting a strong negative selection pressure. Regression analysis of data for ~ 190 countries was carried out to examine associations between 1990 lifetime maternal death rates (LMDR) and current obesity rates. Multivariate regression showed LMDR correlated more strongly with national obesity rates than GDP, calorie intake and physical inactivity. Analyses controlling for confounders via partial correlation show that LMDR explains approximately 11% of the variability of obesity rate between nations. For nations with LMDR above the median (>0.45%), LMDR explains 33% of obesity variance, while calorie intake, GDP and physical inactivity show no association with obesity in these nations. The microevolutionary hypothesis offers a parsimonious explanation of the global nature of the obesity epidemic.

Identification of SNP Markers and Candidate Genes Associated with Major Agronomic Traits in Coffea arabica

Genome-wide association studies (GWASs) allow for inferences about the relationships between genomic variants and phenotypic traits in natural or breeding populations. However, few have used this methodology in Coffea arabica. We aimed to identify chromosomal regions with significant associations between SNP markers and agronomic traits in C. arabica. We used a coffee panel consisting of 195 plants derived from 13 families in F2 generations and backcrosses of crosses between leaf rust-susceptible and -resistant genotypes. The plants were phenotyped for 18 agronomic markers and genotyped for 21,211 SNP markers. A GWAS enabled the identification of 110 SNPs with significant associations (p < 0.05) for several agronomic traits in C. arabica: plant height, plagiotropic branch length, number of vegetative nodes, canopy diameter, fruit size, cercosporiosis incidence, and rust incidence. The effects of each SNP marker associated with the traits were analyzed, such that they can be used for molecular marker-assisted selection. For the first time, a GWAS was used for these important agronomic traits in C. arabica, enabling applications in accelerated coffee breeding through marker-assisted selection and ensuring greater efficiency and time reduction. Furthermore, our findings provide preliminary knowledge to further confirm the genomic loci and potential candidate genes contributing to various structural and disease-related traits of C. arabica.

Somatic Variants Acquired Later in Life Associated with Thoracic Aortic Aneurysms: JAK2 V617F

The JAK2 V617F somatic variant is a well-known driver of myeloproliferative neoplasms (MPN) associated with an increased risk for athero-thrombotic cardiovascular disease. Recent studies have demonstrated its role in the development of thoracic aortic aneurysm (TAA). However, limited clinical information and level of JAK2 V617F burden have been provided for a comprehensive evaluation of potential confounders. A retrospective genotype-first study was conducted to identify carriers of the JAK2 V617F variant from an internal exome sequencing database in Yale DNA Diagnostics Lab. Additionally, the overall incidence of somatic variants in the JAK2 gene across various tissue types in the healthy population was carried out based on reanalysis of SomaMutDB and data from the UK Biobank (UKBB) cohort to compare our dataset to the population prevalence of the variant. In our database of 12,439 exomes, 594 (4.8%) were found to have a thoracic aortic aneurysm (TAA), and 12 (0.049%) were found to have a JAK2 V617F variant. Among the 12 JAK2 V617F variant carriers, five had a TAA (42%), among whom four had an ascending TAA and one had a descending TAA, with a variant allele fraction ranging from 11.2% to 20%. Among these five patients, 60% were female, and average age at diagnosis was 70 (49-79). The mean ascending aneurysm size was 5.05 cm (range 4.6-5.5 cm), and four patients had undergone surgical aortic replacement or repair. UKBB data revealed a positive correlation between the JAK2 V617F somatic variant and aortic valve disease (effect size 0.0086, p = 0.85) and TAA (effect size = 0.004, p = 0.92), although not statistically significant. An unexpectedly high prevalence of TAA in our dataset (5/594, 0.84%) is greater than the prevalence reported before for the general population, supporting its association with TAA. JAK2 V617F may contribute a meaningful proportion of otherwise unexplained aneurysm patients. Additionally, it may imply a potential JAK2-specific disease mechanism in the developmental of TAA, which suggests a possible target of therapy that warrants further investigation.

Spike-phase coupling of subthalamic neurons to posterior opercular cortex predicts speech sound accuracy

Speech provides a rich context for understanding how cortical interactions with the basal ganglia contribute to unique human behaviors, but opportunities for direct intracranial recordings across cortical-basal ganglia networks are rare. We recorded electrocorticographic signals in the cortex synchronously with single units in the basal ganglia during awake neurosurgeries where subjects spoke syllable repetitions. We discovered that individual STN neurons have transient (200ms) spike-phase coupling (SPC) events with multiple cortical regions. The spike timing of STN neurons was coordinated with the phase of theta-alpha oscillations in the posterior supramarginal and superior temporal gyrus during speech planning and production. Speech sound errors occurred when this STN-cortical interaction was delayed. Our results suggest that the STN supports mechanisms of speech planning and auditory-sensorimotor integration during speech production that are required to achieve high fidelity of the phonological and articulatory representation of the target phoneme. These findings establish a framework for understanding cortical-basal ganglia interaction in other human behaviors, and additionally indicate that firing-rate based models are insufficient for explaining basal ganglia circuit behavior.

Variant Effect Prediction in the Age of Machine Learning

Over the years, many computational methods have been created for the analysis of the impact of single amino acid substitutions resulting from single-nucleotide variants in genome coding regions. Historically, all methods have been supervised and thus limited by the inadequate sizes of experimentally curated data sets and by the lack of a standardized definition of variant effect. The emergence of unsupervised, deep learning (DL)-based methods raised an important question: Can machines learn the language of life from the unannotated protein sequence data well enough to identify significant errors in the protein "sentences"? Our analysis suggests that some unsupervised methods perform as well or better than existing supervised methods. Unsupervised methods are also faster and can, thus, be useful in large-scale variant evaluations. For all other methods, however, their performance varies by both evaluation metrics and by the type of variant effect being predicted. We also note that the evaluation of method performance is still lacking on less-studied, nonhuman proteins where unsupervised methods hold the most promise.

Plasminogen missense variants and their involvement in cardiovascular and inflammatory disease

Human plasminogen (PLG), the zymogen of the fibrinolytic protease, plasmin, is a polymorphic protein with two widely distributed codominant alleles, PLG/Asp453 and PLG/Asn453. About 15 other missense or non-synonymous single nucleotide polymorphisms (nsSNPs) of PLG show major, yet different, relative abundances in world populations. Although the existence of these relatively abundant allelic variants is generally acknowledged, they are often overlooked or assumed to be non-pathogenic. In fact, at least half of those major variants are classified as having conflicting pathogenicity, and it is unclear if they contribute to different molecular phenotypes. From those, PLG/K19E and PLG/A601T are examples of two relatively abundant PLG variants that have been associated with PLG deficiencies (PD), but their pathogenic mechanisms are unclear. On the other hand, approximately 50 rare and ultra-rare PLG missense variants have been reported to cause PD as homozygous or compound heterozygous variants, often leading to a debilitating disease known as ligneous conjunctivitis. The true abundance of PD-associated nsSNPs is unknown since they can remain undetected in heterozygous carriers. However, PD variants may also contribute to other diseases. Recently, the ultra-rare autosomal dominant PLG/K311E has been found to be causative of hereditary angioedema (HAE) with normal C1 inhibitor. Two other rare pathogenic PLG missense variants, PLG/R153G and PLG/V709E, appear to affect platelet function and lead to HAE, respectively. Herein, PLG missense variants that are abundant and/or clinically relevant due to association with disease are examined along with their world distribution. Proposed molecular mechanisms are discussed when known or can be reasonably assumed.

Speech-induced suppression and vocal feedback sensitivity in human cortex

Across the animal kingdom, neural responses in the auditory cortex are suppressed during vocalization, and humans are no exception. A common hypothesis is that suppression increases sensitivity to auditory feedback, enabling the detection of vocalization errors. This hypothesis has been previously confirmed in non-human primates, however a direct link between auditory suppression and sensitivity in human speech monitoring remains elusive. To address this issue, we obtained intracranial electroencephalography (iEEG) recordings from 35 neurosurgical participants during speech production. We first characterized the detailed topography of auditory suppression, which varied across superior temporal gyrus (STG). Next, we performed a delayed auditory feedback (DAF) task to determine whether the suppressed sites were also sensitive to auditory feedback alterations. Indeed, overlapping sites showed enhanced responses to feedback, indicating sensitivity. Importantly, there was a strong correlation between the degree of auditory suppression and feedback sensitivity, suggesting suppression might be a key mechanism that underlies speech monitoring. Further, we found that when participants produced speech with simultaneous auditory feedback, posterior STG was selectively activated if participants were engaged in a DAF paradigm, suggesting that increased attentional load can modulate auditory feedback sensitivity.

Recent Advances in Genomic Studies of Gestational Duration and Preterm Birth

Preterm birth (PTB) is the leading cause of infant mortality and morbidity. For several decades, extensive epidemiologic and genetic studies have highlighted the significant contribution of maternal and offspring genetic factors to PTB. This review discusses the challenges inherent in conventional genomic analyses of PTB and underscores the importance of adopting nonconventional approaches, such as analyzing the mother-child pair as a single analytical unit, to disentangle the intertwined maternal and fetal genetic influences. We elaborate on studies investigating PTB phenotypes through 3 levels of genetic analyses: single-variant, multi-variant, and genome-wide variants.

Structural Neuroplasticity Effects of Singing in Chronic Aphasia

Singing-based treatments of aphasia can improve language outcomes, but the neural benefits of group-based singing in aphasia are unknown. Here, we set out to determine the structural neuroplasticity changes underpinning group-based singing-induced treatment effects in chronic aphasia. Twenty-eight patients with at least mild nonfluent poststroke aphasia were randomized into two groups that received a 4-month multicomponent singing intervention (singing group) or standard care (control group). High-resolution T1 images and multishell diffusion-weighted MRI data were collected in two time points (baseline/5 months). Structural gray matter (GM) and white matter (WM) neuroplasticity changes were assessed using language network region of interest-based voxel-based morphometry (VBM) and quantitative anisotropy-based connectometry, and their associations to improved language outcomes (Western Aphasia Battery Naming and Repetition) were evaluated. Connectometry analyses showed that the singing group enhanced structural WM connectivity in the left arcuate fasciculus (AF) and corpus callosum as well as in the frontal aslant tract (FAT), superior longitudinal fasciculus, and corticostriatal tract bilaterally compared with the control group. Moreover, in VBM, the singing group showed GM volume increase in the left inferior frontal cortex (Brodmann area 44) compared with the control group. The neuroplasticity effects in the left BA44, AF, and FAT correlated with improved naming abilities after the intervention. These findings suggest that in the poststroke aphasia group, singing can bring about structural neuroplasticity changes in left frontal language areas and in bilateral language pathways, which underpin treatment-induced improvement in speech production.

Genetically transitional disease: conceptual understanding and applicability to rheumatic disease

In genomic medicine, the concept of genetically transitional disease (GTD) refers to cases in which gene mutation is necessary but not sufficient to cause disease. In this Perspective, we apply this novel concept to rheumatic diseases, which have been linked to hundreds of genetic variants via association studies. These variants are in the 'grey zone' between monogenic variants with large effect sizes and common susceptibility alleles with small effect sizes. Among genes associated with rare autoinflammatory diseases, many low-frequency and/or low-penetrance variants are known to increase susceptibility to systemic inflammation. In autoimmune diseases, hundreds of HLA and non-HLA genetic variants have been revealed to be modest- to moderate-risk alleles. These diseases can be reclassified as GTDs. The same concept could apply to many other human diseases. GTD could improve the reporting of genetic testing results, diagnostic yields, genetic counselling and selection of therapy, as well as facilitating research using a novel approach to human genetic diseases.

Association of serotonin receptor gene polymorphisms with anorexia nervosa: a systematic review and meta-analysis

PURPOSE

Several studies have investigated the association between anorexia nervosa and polymorphisms of genes regulating serotonin neurotransmission, with a focus on the rs6311 polymorphism of 5-HTR2A. However, inconsistent results of these studies and conflicting conclusions of existing meta-analyses complicate the understanding of a possible association. We have updated these results and evaluated the involvement of other serotonin receptor gene polymorphisms in anorexia nervosa.

METHODS

Adhering to PRISMA guidelines, we have searched studies on anorexia nervosa and serotonin-regulating genes published from 1997 to 2022, selected those concerning receptor genes and meta-analyzed the results from twenty candidate gene studies on the 5-HTR2A rs6311 polymorphism and the 5-HTR2C rs6318 polymorphism.

RESULTS

Present analyses reveal an association for the 5-HTR2A rs6311 polymorphism, with G and A alleles, across eighteen studies (2049 patients, 2877 controls; A vs. G allele, Odds Ratio = 1.24; 95% Confidence Interval = 1.06-1.47; p = 0.009). However, after geographic subgrouping, an association emerged only in a Southern European area, involving five studies (722 patients, 773 controls; A vs. G allele, Odds Ratio = 1.82; 95% Confidence Interval = 1.41-2.37; p < 0.00001). No association was observed for the 5-HTR2C rs6318 polymorphism across three studies.

CONCLUSIONS

To date, the involvement in the pathophysiology of anorexia nervosa of the 5-HTR2A rs6311 polymorphism appears limited to a specific genetic and/or environmental context, while that of the 5-HTR2C rs6318 polymorphism seems excluded. Genome-wide association studies and epigenetic studies will likely offer deeper insights of genetic and environmental factors possibly contributing to the disorder.

LEVEL OF EVIDENCE

III Evidence obtained from well-designed cohort or case-control analytic studies. Clinical trial registration PROSPERO registration number: CRD42021246122.

Identifying latent genetic interactions in genome-wide association studies using multiple traits

The "missing" heritability of complex traits may be partly explained by genetic variants interacting with other genes or environments that are difficult to specify, observe, and detect. We propose a new kernel-based method called Latent Interaction Testing (LIT) to screen for genetic interactions that leverages pleiotropy from multiple related traits without requiring the interacting variable to be specified or observed. Using simulated data, we demonstrate that LIT increases power to detect latent genetic interactions compared to univariate methods. We then apply LIT to obesity-related traits in the UK Biobank and detect variants with interactive effects near known obesity-related genes (URL: https://CRAN.R-project.org/package=lit ).

Diallel panel reveals a significant impact of low-frequency genetic variants on gene expression variation in yeast

Unraveling the genetic sources of gene expression variation is essential to better understand the origins of phenotypic diversity in natural populations. Genome-wide association studies identified thousands of variants involved in gene expression variation, however, variants detected only explain part of the heritability. In fact, variants such as low-frequency and structural variants (SVs) are poorly captured in association studies. To assess the impact of these variants on gene expression variation, we explored a half-diallel panel composed of 323 hybrids originated from pairwise crosses of 26 natural Saccharomyces cerevisiae isolates. Using short- and long-read sequencing strategies, we established an exhaustive catalog of single nucleotide polymorphisms (SNPs) and SVs for this panel. Combining this dataset with the transcriptomes of all hybrids, we comprehensively mapped SNPs and SVs associated with gene expression variation. While SVs impact gene expression variation, SNPs exhibit a higher effect size with an overrepresentation of low-frequency variants compared to common ones. These results reinforce the importance of dissecting the heritability of complex traits with a comprehensive catalog of genetic variants at the population level.

Unravelling the Genetic Landscape of Hemiplegic Migraine: Exploring Innovative Strategies and Emerging Approaches

Migraine is a severe, debilitating neurovascular disorder. Hemiplegic migraine (HM) is a rare and debilitating neurological condition with a strong genetic basis. Sequencing technologies have improved the diagnosis and our understanding of the molecular pathophysiology of HM. Linkage analysis and sequencing studies in HM families have identified pathogenic variants in ion channels and related genes, including CACNA1A, ATP1A2, and SCN1A, that cause HM. However, approximately 75% of HM patients are negative for these mutations, indicating there are other genes involved in disease causation. In this review, we explored our current understanding of the genetics of HM. The evidence presented herein summarises the current knowledge of the genetics of HM, which can be expanded further to explain the remaining heritability of this debilitating condition. Innovative bioinformatics and computational strategies to cover the entire genetic spectrum of HM are also discussed in this review.

py_ped_sim - A flexible forward genetic simulator for complex family pedigree analysis

Background

Large-scale family pedigrees are commonly used across medical, evolutionary, and forensic genetics. These pedigrees are tools for identifying genetic disorders, tracking evolutionary patterns, and establishing familial relationships via forensic genetic identification. However, there is a lack of software to accurately simulate different pedigree structures along with genomes corresponding to those individuals in a family pedigree. This limits simulation-based evaluations of methods that use pedigrees.

Results

We have developed a python command-line-based tool called py_ped_sim that facilitates the simulation of pedigree structures and the genomes of individuals in a pedigree. py_ped_sim represents pedigrees as directed acyclic graphs, enabling conversion between standard pedigree formats and integration with the forward population genetic simulator, SLiM. Notably, py_ped_sim allows the simulation of varying numbers of offspring for a set of parents, with the capacity to shift the distribution of sibship sizes over generations. We additionally add simulations for events of misattributed paternity, which offers a way to simulate half-sibling relationships. We validated the accuracy of our software by simulating genomes onto diverse family pedigree structures, showing that the estimated kinship coefficients closely approximated expected values.

Conclusions

py_ped_sim is a user-friendly and open-source solution for simulating pedigree structures and conducting pedigree genome simulations. It empowers medical, forensic, and evolutionary genetics researchers to gain deeper insights into the dynamics of genetic inheritance and relatedness within families.

Novel SNP markers for flowering and seed quality traits in faba bean (Vicia faba L.): characterization and GWAS of a diversity panel

Faba bean (Vicia faba L.) is a legume crop grown in diverse climates worldwide. It has a high potential for increased cultivation to meet the need for more plant-based proteins in human diets, a prerequisite for a more sustainable food production system. Characterization of diversity panels of crops can identify variation in and genetic markers for target traits of interest for plant breeding. In this work, we collected a diversity panel of 220 accessions of faba bean from around the world consisting of gene bank material and commercially available cultivars. The aims of this study were to quantify the phenotypic diversity in target traits to analyze the impact of breeding on these traits, and to identify genetic markers associated with traits through a genome-wide association study (GWAS). Characterization under field conditions at Nordic latitude across two years revealed a large genotypic variation and high broad-sense heritability for eleven agronomic and seed quality traits. Pairwise correlations showed that seed yield was positively correlated to plant height, number of seeds per plant, and days to maturity. Further, susceptibility to bean weevil damage was significantly higher for early flowering accessions and accessions with larger seeds. In this study, no yield penalty was found for higher seed protein content, but protein content was negatively correlated to starch content. Our results showed that while breeding advances in faba bean germplasm have resulted in increased yields and number of seeds per plant, they have also led to a selection pressure towards delayed onset of flowering and maturity. DArTseq genotyping identified 6,606 single nucleotide polymorphisms (SNPs) by alignment to the faba bean reference genome. These SNPs were used in a GWAS, revealing 51 novel SNP markers significantly associated with ten of the assessed traits. Three markers for days to flowering were found in predicted genes encoding proteins for which homologs in other plant species regulate flowering. Altogether, this work enriches the growing pool of phenotypic and genotypic data on faba bean as a valuable resource for developing efficient breeding strategies to expand crop cultivation.

Rare and Common Variants Associated with Alcohol Consumption Identify a Conserved Molecular Network

Genome-wide association studies (GWAS) have identified hundreds of common variants associated with alcohol consumption. In contrast, rare variants have only begun to be studied for their role in alcohol consumption. No studies have examined whether common and rare variants implicate the same genes and molecular networks. To address this knowledge gap, we used publicly available alcohol consumption GWAS summary statistics (GSCAN, N=666,978) and whole exome sequencing data (Genebass, N=393,099) to identify a set of common and rare variants for alcohol consumption. Gene-based analysis of each dataset have implicated 294 (common variants) and 35 (rare variants) genes, including ethanol metabolizing genes ADH1B and ADH1C, which were identified by both analyses, and ANKRD12, GIGYF1, KIF21B, and STK31, which were identified only by rare variant analysis, but have been associated with related psychiatric traits. We then used a network colocalization procedure to propagate the common and rare gene sets onto a shared molecular network, revealing significant overlap. The shared network identified gene families that function in alcohol metabolism, including ADH, ALDH, CYP, and UGT. 74 of the genes in the network were previously implicated in comorbid psychiatric or substance use disorders, but had not previously been identified for alcohol-related behaviors, including EXOC2, EPM2A, CACNB3, and CACNG4. Differential gene expression analysis showed enrichment in the liver and several brain regions supporting the role of network genes in alcohol consumption. Thus, genes implicated by common and rare variants identify shared functions relevant to alcohol consumption, which also underlie psychiatric traits and substance use disorders that are comorbid with alcohol use.

Genetic risk score for gestational weight gain

Gestational weight gain (GWG) involves health consequences for both mother and offspring. Genetic factors seem to play a role in the GWG trait. For small effect sizes of a single genetic polymorphism (SNP), a genetic risk score (GRS) summarizing risk-associated variation from multiple SNPs can serve as an effective approach to genetic association analysis. The aim of the study was to analyze the association between genetic risk score (GRS) and gestational weight gain (GWG). GWG was calculated for a total of 342 healthy Polish women of Caucasian origin, aged 19 to 45 years. The SNPs rs9939609 (FTO), rs6548238 (TMEM18), rs17782313 (MC4R), rs10938397 (GNPDA2), rs10913469 (SEC16B), rs1137101 (LEPR), rs7799039 (LEP), and rs5443 (GNB3) were genotyped using commercial TaqMan SNP assays. A simple genetic risk score was calculated into two ways: GRS1 based on the sum of risk alleles from each of the SNPs, while GRS2 based on the sum of risk alleles of FTO, LEPR, LEP, and GNB3. Positive association between GRS2 and GWG (β = 0.12, p = 0.029) was observed. Genetic risk variants of TMEM18 (p = 0.006, OR = 2.6) and GNB3 (p < 0.001, OR = 3.3) are more frequent in women with increased GWG, but a risk variant of GNPDA2 (p < 0.001, OR = 2.7) is more frequent in women with adequate GWG, and a risk variant of LEPR (p = 0.011, OR = 3.1) in women with decreased GWG. GRS2 and genetic variants of TMEM18, GNB3, GNPDA2, and LEPR are associated with weight gain during pregnancy.

Integration of risk factor polygenic risk score with disease polygenic risk score for disease prediction

Polygenic risk score (PRS) is useful for capturing an individual's genetic susceptibility. However, previous studies have not fully exploited the potential of the risk factor PRS (RFPRS) for disease prediction. We explored the potential of integrating disease-related RFPRSs with disease PRS to enhance disease prediction performance. We constructed 112 RFPRSs and analyzed the association of RFPRSs with diseases to identify disease-related RFPRSs in 700 diseases, using the UK Biobank dataset. We uncovered 6157 statistically significant associations between 247 diseases and 109 RFPRSs. We estimated the disease PRSs of 70 diseases that exhibited statistically significant heritability, to generate RFDiseasemetaPRS-a combined PRS integrating RFPRSs and disease PRS-and compare the prediction performance metrics between RFDiseasemetaPRS and disease PRS. RFDiseasemetaPRS showed better performance for Nagelkerke's pseudo-R2, odds ratio (OR) per 1 SD, net reclassification improvement (NRI) values and difference of R2 considered by variance of R2 in 31 out of 70 diseases. Additionally, we assessed risk classification between two models by examining OR between the top 10% and remaining 90% individuals for the 31 diseases; RFDiseasemetaPRS exhibited better R2, NRI and OR than disease PRS. These findings highlight the importance of utilizing RFDiseasemetaPRS, which can provide personalized healthcare and tailored prevention strategies.

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.

Comprehensive assessment of the genetic characteristics of small for gestational age newborns in NICU: from diagnosis of genetic disorders to prediction of prognosis

BACKGROUND

In China, ~1,072,100 small for gestational age (SGA) births occur annually. These SGA newborns are a high-risk population of developmental delay. Our study aimed to evaluate the genetic profile of SGA newborns in the newborn intensive care unit (NICU) and establish a prognosis prediction model by combining clinical and genetic factors.

METHODS

A cohort of 723 SGA and 1317 appropriate for gestational age (AGA) newborns were recruited between June 2018 and June 2020. Clinical exome sequencing was performed for each newborn. The gene-based rare-variant collapsing analyses and the gene burden test were applied to identify the risk genes for SGA and SGA with poor prognosis. The Gradient Boosting Machine framework was used to generate two models to predict the prognosis of SGA. The performance of two models were validated with an independent cohort of 115 SGA newborns without genetic diagnosis from July 2020 to April 2022. All newborns in this study were recruited through the China Neonatal Genomes Project (CNGP) and were hospitalized in NICU, Children's Hospital of Fudan University, Shanghai, China.

RESULTS

Among the 723 SGA newborns, 88(12.2%) received genetic diagnosis, including 42(47.7%) with monogenic diseases and 46(52.3%) with chromosomal abnormalities. SGA with genetic diagnosis showed higher rates in severe SGA(54.5% vs. 41.9%, P=0.0025) than SGA without genetic diagnosis. SGA with chromosomal abnormalities showed higher incidences of physical and neurodevelopmental delay compared to those with monogenic diseases (45.7% vs. 19.0%, P=0.012). We filtered out 3 genes (ITGB4, TXNRD2, RRM2B) as potential causative genes for SGA and 1 gene (ADIPOQ) as potential causative gene for SGA with poor prognosis. The model integrating clinical and genetic factors demonstrated a higher area under the receiver operating characteristic curve (AUC) over the model based solely on clinical factors in both the SGA-model generation dataset (AUC=0.9[95% confidence interval 0.84-0.96] vs. AUC=0.74 [0.64-0.84]; P=0.00196) and the independent SGA-validation dataset (AUC=0.76 [0.6-0.93] vs. AUC=0.53[0.29-0.76]; P=0.0117).

CONCLUSION

SGA newborns in NICU presented with roughly equal proportions of monogenic and chromosomal abnormalities. Chromosomal disorders were associated with poorer prognosis. The rare-variant collapsing analyses studies have the ability to identify potential causative factors associated with growth and development. The SGA prognosis prediction model integrating genetic and clinical factors outperformed that relying solely on clinical factors. The application of genetic sequencing in hospitalized SGA newborns may improve early genetic diagnosis and prognosis prediction.

A phenome-wide scan reveals convergence of common and rare variant associations

BACKGROUND

Common and rare variants contribute to the etiology of complex traits. However, the extent to which the phenotypic effects of common and rare variants involve shared molecular mediators remains poorly understood. The question is essential to the basic and translational goals of the science of genomics, with critical basic-science, methodological, and clinical consequences.

METHODS

Leveraging the latest release of whole-exome sequencing (WES, for rare variants) and genome-wide association study (GWAS, for common variants) data from the UK Biobank, we developed a metric, the COmmon variant and RAre variant Convergence (CORAC) signature, to quantify the convergence for a broad range of complex traits. We characterized the relationship between CORAC and effective sample size across phenome-wide association studies.

RESULTS

We found that the signature is positively correlated with effective sample size (Spearman ρ = 0.594, P < 2.2e - 16), indicating increased functional convergence of trait-associated genetic variation, across the allele frequency spectrum, with increased power. Sensitivity analyses, including accounting for heteroskedasticity and varying the number of detected association signals, further strengthened the validity of the finding. In addition, consistent with empirical data, extensive simulations showed that negative selection, in line with enhancing polygenicity, has a dampening effect on the convergence signature. Methodologically, leveraging the convergence leads to enhanced association analysis.

CONCLUSIONS

The presented framework for the convergence signature has important implications for fine-mapping strategies and drug discovery efforts. In addition, our study provides a blueprint for the expectation from future large-scale whole-genome sequencing (WGS)/WES and sheds methodological light on post-GWAS studies.

Leveraging Single-Cell Populations to Uncover the Genetic Basis of Complex Traits

The ease and throughput of single-cell genomics have steadily improved, and its current trajectory suggests that surveying single-cell populations will become routine. We discuss the merger of quantitative genetics with single-cell genomics and emphasize how this synergizes with advantages intrinsic to plants. Single-cell population genomics provides increased detection resolution when mapping variants that control molecular traits, including gene expression or chromatin accessibility. Additionally, single-cell population genomics reveals the cell types in which variants act and, when combined with organism-level phenotype measurements, unveils which cellular contexts impact higher-order traits. Emerging technologies, notably multiomics, can facilitate the measurement of both genetic changes and genomic traits in single cells, enabling single-cell genetic experiments. The implementation of single-cell genetics will advance the investigation of the genetic architecture of complex molecular traits and provide new experimental paradigms to study eukaryotic genetics.

The role of admixture in the rare variant contribution to inflammatory bowel disease

BACKGROUND

Identification of rare variants involved in complex, polygenic diseases like Crohn's disease (CD) has accelerated with the introduction of whole exome/genome sequencing association studies. Rare variants can be used in both diagnostic and therapeutic assessments; however, since they are likely to be restricted to specific ancestry groups, their contributions to risk assessment need to be evaluated outside the discovery population. Prior studies implied that the three known rare variants in NOD2 are absent in West African and Asian populations and only contribute in African Americans via admixture.

METHODS

Whole genome sequencing (WGS) data from 3418 African American individuals, 1774 inflammatory bowel disease (IBD) cases, and 1644 controls were used to assess odds ratios and allele frequencies (AF), as well as haplotype-specific ancestral origins of European-derived CD variants discovered in a large exome-wide association study. Local and global ancestry was performed to assess the contribution of admixture to IBD contrasting European and African American cohorts.

RESULTS

Twenty-five rare variants associated with CD in European discovery cohorts are typically five-fold lower frequency in African Americans. Correspondingly, where comparisons could be made, the rare variants were found to have a predicted four-fold reduced burden for IBD in African Americans, when compared to European individuals. Almost all of the rare CD European variants were found on European haplotypes in the African American cohort, implying that they contribute to disease risk in African Americans primarily due to recent admixture. In addition, proportion of European ancestry correlates the number of rare CD European variants each African American individual carry, as well as their polygenic risk of disease. Similar findings were observed for 23 mutations affecting 10 other common complex diseases for which the rare variants were discovered in European cohorts.

CONCLUSIONS

European-derived Crohn's disease rare variants are even more rare in African Americans and contribute to disease risk mainly due to admixture, which needs to be accounted for when performing cross-ancestry genetic assessments.

Genome interpretation in a federated learning context allows the multi-center exome-based risk prediction of Crohn's disease patients

High-throughput sequencing allowed the discovery of many disease variants, but nowadays it is becoming clear that the abundance of genomics data mostly just moved the bottleneck in Genetics and Precision Medicine from a data availability issue to a data interpretation issue. To solve this empasse it would be beneficial to apply the latest Deep Learning (DL) methods to the Genome Interpretation (GI) problem, similarly to what AlphaFold did for Structural Biology. Unfortunately DL requires large datasets to be viable, and aggregating genomics datasets poses several legal, ethical and infrastructural complications. Federated Learning (FL) is a Machine Learning (ML) paradigm designed to tackle these issues. It allows ML methods to be collaboratively trained and tested on collections of physically separate datasets, without requiring the actual centralization of sensitive data. FL could thus be key to enable DL applications to GI on sufficiently large genomics data. We propose FedCrohn, a FL GI Neural Network model for the exome-based Crohn's Disease risk prediction, providing a proof-of-concept that FL is a viable paradigm to build novel ML GI approaches. We benchmark it in several realistic scenarios, showing that FL can indeed provide performances similar to conventional ML on centralized data, and that collaborating in FL initiatives is likely beneficial for most of the medical centers participating in them.

Network medicine-based epistasis detection in complex diseases: ready for quantum computing

Most heritable diseases are polygenic. To comprehend the underlying genetic architecture, it is crucial to discover the clinically relevant epistatic interactions (EIs) between genomic single nucleotide polymorphisms (SNPs)1-3. Existing statistical computational methods for EI detection are mostly limited to pairs of SNPs due to the combinatorial explosion of higher-order EIs. With NeEDL (network-based epistasis detection via local search), we leverage network medicine to inform the selection of EIs that are an order of magnitude more statistically significant compared to existing tools and consist, on average, of five SNPs. We further show that this computationally demanding task can be substantially accelerated once quantum computing hardware becomes available. We apply NeEDL to eight different diseases and discover genes (affected by EIs of SNPs) that are partly known to affect the disease, additionally, these results are reproducible across independent cohorts. EIs for these eight diseases can be interactively explored in the Epistasis Disease Atlas (https://epistasis-disease-atlas.com). In summary, NeEDL is the first application that demonstrates the potential of seamlessly integrated quantum computing techniques to accelerate biomedical research. Our network medicine approach detects higher-order EIs with unprecedented statistical and biological evidence, yielding unique insights into polygenic diseases and providing a basis for the development of improved risk scores and combination therapies.

Genome-Wide Association Analysis of Protein-Coding Variants in IgA Nephropathy

SIGNIFICANCE STATEMENT

Genome-wide association studies have identified nearly 20 IgA nephropathy susceptibility loci. However, most nonsynonymous coding variants, particularly ones that occur rarely or at a low frequency, have not been well investigated. The authors performed a chip-based association study of IgA nephropathy in 8529 patients with the disorder and 23,224 controls. They identified a rare variant in the gene encoding vascular endothelial growth factor A (VEGFA) that was significantly associated with a two-fold increased risk of IgA nephropathy, which was further confirmed by sequencing analysis. They also identified a novel common variant in PKD1L3 that was significantly associated with lower haptoglobin protein levels. This study, which was well-powered to detect low-frequency variants with moderate to large effect sizes, helps expand our understanding of the genetic basis of IgA nephropathy susceptibility.

BACKGROUND

Genome-wide association studies have identified nearly 20 susceptibility loci for IgA nephropathy. However, most nonsynonymous coding variants, particularly those occurring rarely or at a low frequency, have not been well investigated.

METHODS

We performed a three-stage exome chip-based association study of coding variants in 8529 patients with IgA nephropathy and 23,224 controls, all of Han Chinese ancestry. Sequencing analysis was conducted to investigate rare coding variants that were not covered by the exome chip. We used molecular dynamic simulation to characterize the effects of mutations of VEGFA on the protein's structure and function. We also explored the relationship between the identified variants and the risk of disease progression.

RESULTS

We discovered a novel rare nonsynonymous risk variant in VEGFA (odds ratio, 1.97; 95% confidence interval [95% CI], 1.61 to 2.41; P = 3.61×10 -11 ). Further sequencing of VEGFA revealed twice as many carriers of other rare variants in 2148 cases compared with 2732 controls. We also identified a common nonsynonymous risk variant in PKD1L3 (odds ratio, 1.16; 95% CI, 1.11 to 1.21; P = 1.43×10 -11 ), which was associated with lower haptoglobin protein levels. The rare VEGFA mutation could cause a conformational change and increase the binding affinity of VEGFA to its receptors. Furthermore, this variant was associated with the increased risk of kidney disease progression in IgA nephropathy (hazard ratio, 2.99; 95% CI, 1.09 to 8.21; P = 0.03).

CONCLUSIONS

Our study identified two novel risk variants for IgA nephropathy in VEGFA and PKD1L3 and helps expand our understanding of the genetic basis of IgA nephropathy susceptibility.

LDLR and PCSK9 3´UTR variants and their putative effects on microRNA molecular interactions in familial hypercholesterolemia: a computational approach

BACKGROUND

Familial hypercholesterolemia (FH) is caused by pathogenic variants in low-density lipoprotein (LDL) receptor (LDLR) or its associated genes, including apolipoprotein B (APOB), proprotein convertase subtilisin/kexin type 9 (PCSK9), and LDLR adaptor protein 1 (LDLRAP1). However, approximately 40% of the FH patients clinically diagnosed (based on FH phenotypes) may not carry a causal variant in a FH-related gene. Variants located at 3' untranslated region (UTR) of FH-related genes could elucidate mechanisms involved in FH pathogenesis. This study used a computational approach to assess the effects of 3'UTR variants in FH-related genes on miRNAs molecular interactions and to explore the association of these variants with molecular diagnosis of FH.

METHODS AND RESULTS

Exons and regulatory regions of FH-related genes were sequenced in 83 FH patients using an exon-target gene sequencing strategy. In silico prediction tools were used to study the effects of 3´UTR variants on interactions between miRNAs and target mRNAs. Pathogenic variants in FH-related genes (molecular diagnosis) were detected in 44.6% FH patients. Among 59 3'UTR variants identified, LDLR rs5742911 and PCSK9 rs17111557 were associated with molecular diagnosis of FH, whereas LDLR rs7258146 and rs7254521 and LDLRAP1 rs397860393 had an opposite effect (p < 0.05). 3´UTR variants in LDLR (rs5742911, rs7258146, rs7254521) and PCSK9 (rs17111557) disrupt interactions with several miRNAs, and more stable bindings were found with LDLR (miR-4435, miR-509-3 and miR-502) and PCSK9 (miR-4796).

CONCLUSION

LDLR and PCSK9 3´UTR variants disturb miRNA:mRNA interactions that could affect gene expression and are potentially associated with molecular diagnosis of FH.

Implications of combined NOD2 and other gene mutations in autoinflammatory diseases

NOD-like receptors (NLRs) are intracellular sensors associated with systemic autoinflammatory diseases (SAIDs). We investigated the largest monocentric cohort of patients with adult-onset SAIDs for coinheritance of low frequency and rare mutations in NOD2 and other autoinflammatory genes. Sixty-three patients underwent molecular testing for SAID gene panels after extensive clinical workups. Whole exome sequencing data from the large Atherosclerosis Risk in Communities (ARIC) study of individuals of European-American ancestry were used as control. Of 63 patients, 44 (69.8%) were found to carry combined gene variants in NOD2 and another gene (Group 1), and 19 (30.2%) were carriers only for NOD2 variants (Group 2). The genetic variant combinations in SAID patients were digenic in 66% (NOD2/MEFV, NOD2/NLRP12, NOD2/NLRP3, and NOD2/TNFRSF1A) and oligogenic in 34% of cases. These variant combinations were either absent or significantly less frequent in the control population. By phenotype-genotype correlation, approximately 40% of patients met diagnostic criteria for a specific SAID, and 60% had mixed diagnoses. There were no statistically significant differences in clinical manifestations between the two patient groups except for chest pain. Due to overlapping phenotypes and mixed genotypes, we have suggested a new term, "Mixed NLR-associated Autoinflammatory Disease ", to describe this disease scenario. Gene variant combinations are significant in patients with SAIDs primarily presenting with mixed clinical phenotypes. Our data support the proposition that immunological disease expression is modified by genetic background and environmental exposure. We provide a preliminary framework in diagnosis, management, and interpretation of the clinical scenario.

The functional impact of rare variation across the regulatory cascade

Each human genome has tens of thousands of rare genetic variants; however, identifying impactful rare variants remains a major challenge. We demonstrate how use of personal multi-omics can enable identification of impactful rare variants by using the Multi-Ethnic Study of Atherosclerosis, which included several hundred individuals, with whole-genome sequencing, transcriptomes, methylomes, and proteomes collected across two time points, 10 years apart. We evaluated each multi-omics phenotype's ability to separately and jointly inform functional rare variation. By combining expression and protein data, we observed rare stop variants 62 times and rare frameshift variants 216 times as frequently as controls, compared to 13-27 times as frequently for expression or protein effects alone. We extended a Bayesian hierarchical model, "Watershed," to prioritize specific rare variants underlying multi-omics signals across the regulatory cascade. With this approach, we identified rare variants that exhibited large effect sizes on multiple complex traits including height, schizophrenia, and Alzheimer's disease.

Genomic approaches to explore susceptibility and pathogenesis of alcohol use disorder and alcohol-associated liver disease

Excessive alcohol use is a major risk factor for the development of an alcohol use disorder (AUD) and contributes to a wide variety of other medical illnesses, including alcohol-associated liver disease (ALD). Both AUD and ALD are complex and causally interrelated diseases, and multiple factors other than alcohol consumption are implicated in the disease pathogenesis. While the underlying pathophysiology of AUD and ALD is complex, there is substantial evidence for a genetic susceptibility of both diseases. Current genome-wide association studies indicate that the genes associated with clinical AUD only poorly overlap with the genes identified for heavy drinking and, in turn, neither overlap with the genes identified for ALD. Uncovering the main genetic factors will enable us to identify molecular drivers underlying the pathogenesis, discover potential targets for therapy, and implement patient care early in disease progression. In this review, we described multiple genomic approaches and their implications to investigate the susceptibility and pathogenesis of both AUD and ALD. We concluded our review with a discussion of the knowledge gaps and future research on genomic studies in these 2 diseases.

Genomic and transcriptomic analyses reveal polygenic architecture for ecologically important traits in aspen (Populus tremuloides Michx.)

Intraspecific genetic variation in foundation species such as aspen (Populus tremuloides Michx.) shapes their impact on forest structure and function. Identifying genes underlying ecologically important traits is key to understanding that impact. Previous studies, using single-locus genome-wide association (GWA) analyses to identify candidate genes, have identified fewer genes than anticipated for highly heritable quantitative traits. Mounting evidence suggests that polygenic control of quantitative traits is largely responsible for this "missing heritability" phenomenon. Our research characterized the genetic architecture of 30 ecologically important traits using a common garden of aspen through genomic and transcriptomic analyses. A multilocus association model revealed that most traits displayed a highly polygenic architecture, with most variation explained by loci with small effects (likely below the detection levels of single-locus GWA methods). Consistent with a polygenic architecture, our single-locus GWA analyses found only 38 significant SNPs in 22 genes across 15 traits. Next, we used differential expression analysis on a subset of aspen genets with divergent concentrations of salicinoid phenolic glycosides (key defense traits). This complementary method to traditional GWA discovered 1243 differentially expressed genes for a polygenic trait. Soft clustering analysis revealed three gene clusters (241 candidate genes) involved in secondary metabolite biosynthesis and regulation. Our work reveals that ecologically important traits governing higher-order community- and ecosystem-level attributes of a foundation forest tree species have complex underlying genetic structures and will require methods beyond traditional GWA analyses to unravel.

Weighted Selection Probability to Prioritize Susceptible Rare Variants in Multi-Phenotype Association Studies with Application to a Soybean Genetic Data Set

Rare variant association studies with multiple traits or diseases have drawn a lot of attention since association signals of rare variants can be boosted if more than one phenotype outcome is associated with the same rare variants. Most of the existing statistical methods to identify rare variants associated with multiple phenotypes are based on a group test, where a pre-specified genetic region is tested one at a time. However, these methods are not designed to locate susceptible rare variants within the genetic region. In this article, we propose new statistical methods to prioritize rare variants within a genetic region when a group test for the genetic region identifies a statistical association with multiple phenotypes. It computes the weighted selection probability (WSP) of individual rare variants and ranks them from largest to smallest according to their WSP. In simulation studies, we demonstrated that the proposed method outperforms other statistical methods in terms of true positive selection, when multiple phenotypes are correlated with each other. We also applied it to our soybean single nucleotide polymorphism (SNP) data with 13 highly correlated amino acids, where we identified some potentially susceptible rare variants in chromosome 19.