Rare and common variants: twenty arguments

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.

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.

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.

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.

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.

The complete and fully-phased diploid genome of a male Han Chinese

Since the release of the complete human genome, the priority of human genomic study has now been shifting towards closing gaps in ethnic diversity. Here, we present a fully phased and well-annotated diploid human genome from a Han Chinese male individual (CN1), in which the assemblies of both haploids achieve the telomere-to-telomere (T2T) level. Comparison of this diploid genome with the CHM13 haploid T2T genome revealed significant variations in the centromere. Outside the centromere, we discovered 11,413 structural variations, including numerous novel ones. We also detected thousands of CN1 alleles that have accumulated high substitution rates and a few that have been under positive selection in the East Asian population. Further, we found that CN1 outperforms CHM13 as a reference genome in mapping and variant calling for the East Asian population owing to the distinct structural variants of the two references. Comparison of SNP calling for a large cohort of 8869 Chinese genomes using CN1 and CHM13 as reference respectively showed that the reference bias profoundly impacts rare SNP calling, with nearly 2 million rare SNPs miss-called with different reference genomes. Finally, applying the CN1 as a reference, we discovered 5.80 Mb and 4.21 Mb putative introgression sequences from Neanderthal and Denisovan, respectively, including many East Asian specific ones undetected using CHM13 as the reference. Our analyses reveal the advances of using CN1 as a reference for population genomic studies and paleo-genomic studies. This complete genome will serve as an alternative reference for future genomic studies on the East Asian population.

Translating pharmacogenomic sequencing data into drug response predictions-How to interpret variants of unknown significance

The rapid development of sequencing technologies during the past 20 years has provided a variety of methods and tools to interrogate human genomic variations at the population level. Pharmacogenes are well known to be highly polymorphic and a plethora of pharmacogenomic variants has been identified in population sequencing data. However, so far only a small number of these variants have been functionally characterized regarding their impact on drug efficacy and toxicity and the significance of the vast majority remains unknown. It is therefore of high importance to develop tools and frameworks to accurately infer the effects of pharmacogenomic variants and, eventually, aggregate the effect of individual variations into personalized drug response predictions. To address this challenge, we here first describe the technological advances, including sequencing methods and accompanying bioinformatic processing pipelines that have enabled reliable variant identification. Subsequently, we highlight advances in computational algorithms for pharmacogenomic variant interpretation and discuss the added value of emerging strategies, such as machine learning and the integrative use of omics techniques that have the potential to further contribute to the refinement of personalized pharmacological response predictions. Lastly, we provide an overview of experimental and clinical approaches to validate in silico predictions. We conclude that the iterative feedback between computational predictions and experimental validations is likely to rapidly improve the accuracy of pharmacogenomic prediction models, which might soon allow for an incorporation of the entire pharmacogenetic profile into personalized response predictions.

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.

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

Genome-wide association studies of complex traits frequently find that SNP-based estimates of heritability are considerably smaller than estimates from classic family-based studies. This 'missing' heritability may be partly explained by genetic variants interacting with other genes or environments that are difficult to specify, observe, and detect. To circumvent these challenges, we propose a new method to detect genetic interactions that leverages pleiotropy from multiple related traits without requiring the interacting variable to be specified or observed. Our approach, Latent Interaction Testing (LIT), uses the observation that correlated traits with shared latent genetic interactions have trait variance and covariance patterns that differ by genotype. LIT examines the relationship between trait variance/covariance patterns and genotype using a flexible kernel-based framework that is computationally scalable for biobank-sized datasets with a large number of traits. We first use simulated data to demonstrate that LIT substantially increases power to detect latent genetic interactions compared to a trait-by-trait univariate method. We then apply LIT to four obesity-related traits in the UK Biobank and detect genetic variants with interactive effects near known obesity-related genes. Overall, we show that LIT, implemented in the R package lit, uses shared information across traits to improve detection of latent genetic interactions compared to standard approaches.

Polygenic scores, and the genome-wide association studies they derive from, will have difficulty identifying genes that predispose one to develop a social behavioral trait

Polygenic scores (PGSs) have several limitations. They are confounded with environmental effects on behavior and cannot be used to study how mutations affect brain function and behavior. For this, mutations with large effects, which often arise in only one geographical population are needed. Genome-wide association studies (GWASs), commonly used for identifying mutations, have difficulty detecting these mutations. A strategy that overcomes this challenge is discussed.

Novelties in Brugada Syndrome: Complex Genetics, Risk Stratification, and Catheter Ablation

Brugada syndrome (BrS) is an inherited arrhythmia syndrome with distinctive electrocardiographic abnormalities in the right precordial leads and predisposes to ventricular arrhythmias and sudden cardiac death in otherwise healthy patients. Its complex genetic architecture and pathophysiological mechanism are not yet completely understood, and risk stratification remains challenging, particularly in patients at intermediate risk of arrhythmic events. Further understanding of its complex genetic architecture may help improving future risk stratification, and advances in management may contribute to alternatives to implantable cardioverter-defibrillators. Here, the authors review the latest insights and developments in BrS.

Human Exome Sequencing and Prospects for Predictive Medicine: Analysis of International Data and Own Experience

Today, whole-exome sequencing (WES) is used to conduct the massive screening of structural and regulatory genes in order to identify the allele frequencies of disease-associated polymorphisms in various populations and thus detect pathogenic genetic changes (mutations or polymorphisms) conducive to malfunctional protein sequences. With its extensive capabilities, exome sequencing today allows both the diagnosis of monogenic diseases (MDs) and the examination of seemingly healthy populations to reveal a wide range of potential risks prior to disease manifestation (in the future, exome sequencing may outpace costly and less informative genome sequencing to become the first-line examination technique). This review establishes the human genetic passport as a new WES-based clinical concept for the identification of new candidate genes, gene variants, and molecular mechanisms in the diagnosis, prediction, and treatment of monogenic, oligogenic, and multifactorial diseases. Various diseases are addressed to demonstrate the extensive potential of WES and consider its advantages as well as disadvantages. Thus, WES can become a general test with a broad spectrum pf applications, including opportunistic screening.

Sex differences in deleterious genetic variants in a haplodiploid social insect

Deleterious variants are selected against but can linger in populations at low frequencies for long periods of time, decreasing fitness and contributing to disease burden in humans and other species. Deleterious variants occur at low frequency but distinguishing deleterious variants from low-frequency neutral variation is challenging based on population genomics data alone. As a result, we have little sense of the number and identity of deleterious variants in wild populations. For haplodiploid species, it has been hypothesised that deleterious alleles will be directly exposed to selection in haploid males, but selection can be masked in diploid females when deleterious variants are recessive, resulting in more efficient purging of deleterious mutations in males. Therefore, comparisons of the differences between haploid and diploid genomes from the same population may be a useful method for inferring rare deleterious variants. This study provides the first formal test of this hypothesis. Using wild populations of Northern paper wasps (Polistes fuscatus), we find that males have fewer missense and nonsense variants per generation than females from the same population. Allele frequency differences are especially pronounced for rare missense and nonsense variants and these differences lead to a lower mutational load in males than females. Based on these data we infer that many highly deleterious mutations are segregating in the paper wasp population. Stronger selection against deleterious alleles in haploid males may have implications for adaptation in other haplodiploid insects and provides evidence that wild populations harbour abundant deleterious variants.

Predicting the effect of mutations to investigate recent events of selection across 60,472 Escherichia coli strains

Microbial genomics studies focusing on the dynamics of selection have often used a small number of distant genomes. As a result, they could only analyze mutations that had become fixed during the divergence between species. However, thousands of genomes of some species are now available in public databases, thanks to high-throughput sequencing. These data provide a more complete picture of the polymorphisms segregating within a species, offering a unique insight into the processes that shape the recent evolution of a species. In this study, we present GLASS (Gene-Level Amino-acid Score Shift), a selection test that is based on the predicted effects of amino acid changes. By comparing the distribution of effects of mutations observed in a gene to the expectation in the absence of selection, GLASS can quantify the intensity of selection. We applied GLASS to a dataset of 60,472 Escherichia coli strains and used this to reexamine the longstanding debate about the role of essentiality versus expression level in the rate of protein evolution. We found that selection has contrasting short-term and long-term dynamics, with essential genes being subject to strong purifying selection in the short term, while expression level determines the rate of gene evolution in the long term. GLASS also found an overrepresentation of inactivating mutations in specific transcription factors, such as efflux pump repressors, which is consistent with selection for antibiotic resistance. These gene-inactivating polymorphisms do not reach fixation, suggesting another contrast between short-term fitness gains and long-term counterselection.

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.

Genomics of founders for conservation breeding: the Jasper caribou case

Conservation breeding programs are increasingly used as recovery actions for wild animals; bringing founders into captivity to rear captive populations for future reintroduction into the wild. The International Union for the Conservation of Nature recommends that founders should come from genetically close populations and should have sufficient genetic diversity to avoid mating among relatives. Genomic data are highly informative for evaluating founders due to their high resolution and ability to capture adaptive divergence, yet, their application in that context remains limited. Woodland caribou are federally listed as a Species at Risk in Canada, with several populations facing extirpation, such as those in the Rocky Mountains of Alberta and British Columbia (BC). To prevent local extirpation, Jasper National Park (JNP) is proposing a conservation breeding program. We examined single nucleotide polymorphisms for 144 caribou from 11 populations encompassing a 200,0002 km area surrounding JNP to provide information useful for identifying appropriate founders for this program. We found that this area likely hosts a caribou metapopulation historically characterized by high levels of gene flow, which indicates that multiple sources of founders would be appropriate for initiating a breeding program. However, population structure and adaptive divergence analyses indicate that JNP caribou are closest to populations in the BC Columbia range, which also have suitable genetic diversity for conservation breeding. We suggest that collaboration among jurisdictions would be beneficial to implement the program to promote recovery of JNP caribou and possibly other caribou populations in the surrounding area, which is strategically at the periphery of the distribution of this endangered species.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10592-023-01540-3.

Invited review: Good practices in genome-wide association studies to identify candidate sequence variants in dairy cattle

Genotype data from dairy cattle selection programs have greatly facilitated GWAS to identify variants related to economic traits. Results can enhance the accuracy of genomic prediction, analyze more complex models that go beyond additive effects, elucidate the genetic architecture of a trait, and finally, decipher the underlying biology of traits. The entire process, comprising data generation, quality control, statistical analyses, interpretation of association results, and linking results to biology should be designed and executed to minimize the generation of false-positive and false-negative associations and misleading links to biological processes. This review aims to provide general guidelines for data analysis that address data quality control, association tests, adjustment for population stratification, and significance evaluation to improve the reliability of conclusions. We also provide guidance on post-GWAS strategy and the interpretation of results. These guidelines are tailored to dairy cattle, which are characterized by long-range linkage disequilibrium, large half-sib families, and routinely collected phenotypes, requiring different approaches than those applied in human GWAS. We discuss common limitations and challenges that have been overlooked in the analysis and interpretation of GWAS to identify candidate sequence variants in dairy cattle.

Functional characterization of human genomic variation linked to polygenic diseases

The burden of human disease lies predominantly in polygenic diseases. Since the early 2000s, genome-wide association studies (GWAS) have identified genetic variants and loci associated with complex traits. These have ranged from variants in coding sequences to mutations in regulatory regions, such as promoters and enhancers, as well as mutations affecting mediators of mRNA stability and other downstream regulators, such as 5' and 3'-untranslated regions (UTRs), long noncoding RNA (lncRNA), and miRNA. Recent research advances in genetics have utilized a combination of computational techniques, high-throughput in vitro and in vivo screening modalities, and precise genome editing to impute the function of diverse classes of genetic variants identified through GWAS. In this review, we highlight the vastness of genomic variants associated with polygenic disease risk and address recent advances in how genetic tools can be used to functionally characterize them.