Exome-wide Association Study Research Articles

Whole exome sequencing identified six novel genes for depressive symptoms.

Previous genome-wide association studies of depression have primarily focused on common variants, limiting our comprehensive understanding of the genetic architecture. In contrast, whole-exome sequencing can capture rare coding variants, helping to explore the phenotypic consequences of altering protein-coding genes. Here, we conducted a large-scale exome-wide association study on 296,199 participants from the UK Biobank, assessing their depressive symptom scores through the Patient Health Questionnaire-4. We identified 22 genes associated with depressive symptoms, including 6 newly discovered genes (TRIM27, UBD, SVOP, ADGRB2, IRF2BPL, and ANKRD12). Both ontology enrichment analysis and plasma proteomics association analysis consistently revealed that the identified genes were associated with immune responses. Furthermore, we identified associations between these genes and brain regions related to depression, such as anterior cingulate cortex and orbitofrontal cortex. Additionally, phenome-wide association analysis demonstrated that TRIM27 and UBD were associated with neuropsychiatric, cognitive, biochemistry, and inflammatory traits. Our findings offer new insights into the potential mechanisms and genetic architecture of depressive symptoms.

Genome-Wide and Exome-Wide Association Study Identifies Genetic Underpinning of Comorbidity between Myocardial Infarction and Severe Mental Disorders.

Myocardial Infarction (MI) and severe mental disorders (SMDs) are two types of highly prevalent and complex disorders and seem to have a relatively high possibility of mortality. However, the contributions of common and rare genetic variants to their comorbidity arestill unclear. We conducted a combined genome-wide association study (GWAS) and exome-wide association study (EWAS) approach. Using gene-based and gene-set association analyses based on the results of GWAS, we found the common genetic underpinnings of nine genes (GIGYF2, KCNJ13, PCCB, STAG1, HLA-C, HLA-B, FURIN, FES, and SMG6) and nine pathways significantly shared between MI and SMDs. Through Mendelian randomization analysis, we found that twenty-seven genes were potential causal genes for SMDs and MI. Based on the exome sequencing data of MI and SMDs patients from the UK Biobank, we found that MUC2 was exome-wide significant in the two diseases. The gene-set analyses of the exome-wide association study indicated that pathways related to insulin processing androgen catabolic process and angiotensin receptor binding may be involved in the comorbidity between SMDs and MI. We also found that six candidate genes were reported to interact with known therapeutic drugs based on the drug-gene interaction information in DGIdb. Altogether, this study revealed the overlap of common and rare genetic underpinning between SMDs and MI and may provide useful insights for their mechanism study and therapeutic investigations.

Joint testing of rare variant burden scores using non-negative least squares

Gene-based burden tests are a popular and powerful approach for analysis of exome-wide association studies. These approaches combine sets of variants within a gene into a single burden score that is then tested for association. Typically, a range of burden scores are calculated and tested across a range of annotation classes and frequency bins. Correlation between these tests can complicate the multiple testing correction and hamper interpretation of the results. We introduce a method called the sparse burden association test (SBAT) that tests the joint set of burden scores under the assumption that causal burden scores act in the same effect direction. The method simultaneously assesses the significance of the model fit and selects the set of burden scores that best explain the association at the same time. Using simulated data, we show that the method is well calibrated and highlight scenarios where the test outperforms existing gene-based tests. We apply the method to 73 quantitative traits from the UK Biobank, showing that SBAT is a valuable additional gene-based test when combined with other existing approaches. This test is implemented in the REGENIE software.

Genetic Diversity in Schizophrenia: Developmental Implications of Ultra-Rare, Protein-Truncating Mutations.

The genetic basis of schizophrenia (SZ) remains elusive despite its characterization as a highly heritable disorder. This incomplete understanding has led to stagnation in therapeutics and treatment, leaving many suffering with insufficient relief from symptoms. However, recent large-cohort genome- and exome-wide association studies have provided insights into the underlying genetic machinery. The scale of these studies allows for the identification of ultra-rare mutations that confer substantial disease risk, guiding clinicians and researchers toward general classes of genes that are central to SZ etiology. One such large-scale collaboration effort by the Schizophrenia Exome Sequencing Meta-Analysis consortium identified ten, high-risk, ultra-rare, protein-truncating variants, providing the clearest picture to date of the dysfunctional gene products that substantially increase risk for SZ. While genetic studies of SZ provide valuable information regarding "what" genes are linked with the disorder, it is an open question as to "when" during brain development these genetic mutations impose deleterious effects. To shed light on this unresolved aspect of SZ etiology, we queried the BrainSpan developmental mRNA expression database for these ten high-risk genes and discovered three general expression trajectories throughout pre- and postnatal brain development. The elusiveness of SZ etiology, we infer, is not only borne out of the genetic heterogeneity across clinical cases, but also in our incomplete understanding of how genetic mutations perturb neurodevelopment during multiple critical periods. We contextualize this notion within the National Institute of Mental Health's Research Domain Criteria framework and emphasize the utility of considering both genetic variables and developmental context in future studies.

Lack of a genetic risk continuum between pubertal timing in the general population and idiopathic hypogonadotropic hypogonadism.

Pubertal timing is a highly heritable trait in the general population. Recently, a large-scale exome-wide association study has implicated rare variants in six genes (KDM4C, MC3R, MKRN3, PDE10A, TACR3, and ZNF483) as genetic determinants of pubertal timing within the general population. Two of the genes (TACR3, MKRN3) are already implicated in extreme disorders of pubertal timing. This observation suggests that there may be a pervasive "genetic risk continuum" wherein genes that govern pubertal timing in the general population, by extension, may also be causal for rare Mendelian disorders of pubertal timing. Hence, we hypothesized that the four novel genes linked to pubertal timing in the population will also contribute to idiopathic hypogonadotropic hypogonadism (IHH), a genetic disorder characterized by absent puberty. Exome sequencing data from 1322 unrelated IHH probands were reviewed for rare sequence variants (RSVs) (minor allele frequency bins: <1%; <0.1%; <0.01%) in the six genes linked to puberty in the general population. A gene-based rare variant association testing (RVAT) was performed between the IHH cohort and a reference public genomic sequences repository-the Genome Aggregation Database (gnomAD). As expected, RVAT analysis showed that RSVs in TACR3, a known IHH gene, were significantly enriched in the IHH cohort compared to gnomAD cohort across all three MAF bins. However, RVAT analysis of the remaining five genes failed to show any RSV enrichment in the IHH cohort across all MAF bins. Our findings argue strongly against a pervasive genetic risk continuum between pubertal timing in the general population and extreme pubertal phenotypes. The biologic basis of such distinct genetic architectures' merits further evaluation.

Genetic variation in CCDC93 is associated with elevated central systolic blood pressure, impaired arterial relaxation, and mitochondrial dysfunction.

Genetic studies of blood pressure (BP) traits to date have been performed on conventional measures by brachial cuff sphygmomanometer for systolic BP (SBP) and diastolic BP, integrating several physiologic occurrences. Genetic associations with central SBP (cSBP) have not been well-studied. Genetic discovery studies of BP have been most often performed in European-ancestry samples. Here, we investigated genetic associations with cSBP in a Chinese population and functionally validated the impact of a novel associated coiled-coil domain containing 93 (CCDC93) gene on BP regulation. An exome-wide association study (EWAS) was performed using a mixed linear model of non-invasive cSBP and peripheral BP traits in a Han Chinese population (N = 5,954) from Beijing, China genotyped with a customized Illumina ExomeChip array. We identified four SNP-trait associations with three SNPs, including two novel associations (rs2165468-SBP and rs33975708-cSBP). rs33975708 is a coding variant in the CCDC93 gene, c.535C>T, p.Arg179Cys (MAF = 0.15%), and was associated with increased cSBP (β = 29.3 mmHg, P = 1.23x10-7). CRISPR/Cas9 genome editing was used to model the effect of Ccdc93 loss in mice. Homozygous Ccdc93 deletion was lethal prior to day 10.5 of embryonic development. Ccdc93+/- heterozygous mice were viable and morphologically normal, with 1.3-fold lower aortic Ccdc93 protein expression (P = 0.0041) and elevated SBP as compared to littermate Ccdc93+/+ controls (110±8 mmHg vs 125±10 mmHg, P = 0.016). Wire myography of Ccdc93+/- aortae showed impaired acetylcholine-induced relaxation and enhanced phenylephrine-induced contraction. RNA-Seq transcriptome analysis of Ccdc93+/- mouse thoracic aortae identified significantly enriched pathways altered in fatty acid metabolism and mitochondrial metabolism. Plasma free fatty acid levels were elevated in Ccdc93+/- mice (96±7mM vs 124±13mM, P = 0.0031) and aortic mitochondrial dysfunction was observed through aberrant Parkin and Nix protein expression. Together, our genetic and functional studies support a novel role of CCDC93 in the regulation of BP through its effects on vascular mitochondrial function and endothelial function.

Whole exome sequencing analyses identified novel genes for Alzheimer's disease and related dementia.

The heritability of Alzheimer's disease (AD) is estimated to be 58%-79%. However, known genes can only partially explain the heritability. Here, we conducted gene-based exome-wide association study (ExWAS) of rare variants and single-variant ExWAS of common variants, utilizing data of 54,569 clinically diagnosed/proxy AD and related dementia (ADRD) and 295,421 controls from the UK Biobank. Gene-based ExWAS identified 11 genes predicting a higher ADRD risk, including five novel ones, namely FRMD8, DDX1, DNMT3L, MORC1, and TGM2, along with six previously reported ones, SORL1, GRN, PSEN1, ABCA7, GBA, and ADAM10. Single-variant ExWAS identified two ADRD-associated novel genes, SLCO1C1 and NDNF. The identified genes were predominantly enriched in amyloid-β process pathways, microglia, and brain regions like hippocampus. The druggability evidence suggests that DDX1, DNMT3L, TGM2, SLCO1C1, and NDNF could be effective drug targets. Our study contributes to the current body of evidence on the genetic etiology of ADRD. Gene-based analyses of rare variants identified five novel genes for Alzheimer's disease and related dementia (ADRD), including FRMD8, DDX1, DNMT3L, MORC1, and TGM2. Single-variant analyses of common variants identified two novel genes for ADRD, including SLCO1C1 and NDNF. The identified genes were predominantly enriched in amyloid-β process pathways, microglia, and brain regions like hippocampus. DDX1, DNMT3L, TGM2, SLCO1C1, and NDNF could be effective drug targets.

Exome-wide association study identifies KDELR3 mutations in extreme myopia

Extreme myopia (EM), defined as a spherical equivalent (SE) ≤ −10.00 diopters (D), is one of the leading causes of sight impairment. Known EM-associated variants only explain limited risk and are inadequate for clinical decision-making. To discover risk genes, we performed a whole-exome sequencing (WES) on 449 EM individuals and 9606 controls. We find a significant excess of rare protein-truncating variants (PTVs) in EM cases, enriched in the retrograde vesicle-mediated transport pathway. Employing single-cell RNA-sequencing (scRNA-seq) and a single-cell polygenic burden score (scPBS), we pinpointed PI16 + /SFRP4+ fibroblasts as the most relevant cell type. We observed that KDELR3 is highly expressed in scleral fibroblast and involved in scleral extracellular matrix (ECM) organization. The zebrafish model revealed that kdelr3 downregulation leads to elongated ocular axial length and increased lens diameter. Together, our study provides insight into the genetics of EM in humans and highlights KDELR3’s role in EM pathogenesis.

Exome sequencing identifies novel genetic variants associated with varicose veins.

Varicose veins (VV) are one of the common human diseases, but the role of genetics in its development is not fully understood. We conducted an exome-wide association study of VV using whole-exome sequencing data from the UK Biobank, and focused on common and rare variants using single-variant association analysis and gene-level collapsing analysis. A total of 13,823,269 autosomal genetic variants were obtained after quality control. We identified 36 VV-related independent common variants mapping to 34 genes by single-variant analysis and three rare variant genes (PIEZO1, ECE1, FBLN7) by collapsing analysis, and most associations between genes and VV were replicated in FinnGen. PIEZO1 was the closest gene associated with VV (P = 5.05 × 10-31), and it was found to reach exome-wide significance in both single-variant and collapsing analyses. Two novel rare variant genes (ECE1 and METTL21A) associated with VV were identified, of which METTL21A was associated only with females. The pleiotropic effects of VV-related genes suggested that body size, inflammation, and pulmonary function are strongly associated with the development of VV. Our findings highlight the importance of causal genes for VV and provide new directions for treatment.

Exome-wide association study for replication of rare variants affecting the severity of COVID-19 in the Russian population

BACKGROUND: Human genotype is a factor that determines the severity of COVID-19. Previously, a large-scale whole-genome association study of the COVID-19 Host Genetics Initiative (2021) investigated the association of genetic variants at multiple loci with COVID-19 severity. The genetic variants that have the greatest effect on COVID-19 severity are expected to have a low frequency in the population. Therefore, the study of rare variants may provide additional insights into the disease pathogenesis and thus help in the development of prevention and treatment options. AIM: To search for genes enriched for rare genetic variants associated with COVID-19 severity in the Russian population by replication analysis. METHODS: The clinical exome of a Russian cohort of patients was sequenced based on the St. Petersburg State Budgetary Institution “City Hospital No. 40” and St Petersburg University. The study used biomaterial from patients hospitalized at City Hospital No. 40 diagnosed with COVID-19 and healthy individuals (population control group). The severity of the course of COVID-19 was determined according to the results of lung computed tomography. The list of genes for subsequent replication was generated by a literature review. Burden test methods were used for the replication analysis of genes associated with COVID-19 severity. RESULTS: In total, 701 clinical exomes were sequenced from 263 individuals with severe COVID-19 and 438 healthy individuals. In the literature review, 18 genes associated with severe COVID-19 were included in the replication analysis. The replication analysis did not identify any genes whose association with severe COVID-19 was confirmed in the study cohort. CONCLUSION: The replication analysis did not identify any genes that showed a significant association between the functional variant enrichment and COVID-19 severity. However, the direction of the correlation was consistent with the findings of previous studies. Expanding the study cohort would increase the power of the tests and allow us to detect additional rare variants that influence the severity of COVID-19 progression.

Whole exome sequencing analyses reveal novel genes in telomere length and their biomedical implications.

Telomere length is a putative biomarker of aging and is associated with multiple age-related diseases. There are limited data on the landscape of rare genetic variations in telomere length. Here, we systematically characterize the rare variant associations with leukocyte telomere length (LTL) through exome-wide association study (ExWAS) among 390,231 individuals in the UK Biobank. We identified 18 robust rare-variant genes for LTL, most of which estimated effects on LTL were significant (> 0.2 standard deviation per allele). The biological functions of the rare-variant genes were associated with telomere maintenance and capping and several genes were specifically expressed in the testis. Three novel genes (ASXL1, CFAP58, and TET2) associated with LTL were identified. Phenotypic association analyses indicated significant associations of ASXL1 and TET2 with cancers, age-related diseases, blood assays, and cardiovascular traits. Survival analyses suggested that carriers of ASXL1 or TET2 variants were at increased risk for cancers; diseases of the circulatory, respiratory, and genitourinary systems; and all-cause and cause-specific deaths. The CFAP58 carriers were at elevated risk of deaths due to cancers. Collectively, the present whole exome sequencing study provides novel insights into the genetic landscape of LTL, identifying novel genes associated with LTL and their implications on human health and facilitating a better understanding of aging, thus pinpointing the genetic relevance of LTL with clonal hematopoiesis, biomedical traits, and health-related outcomes.

The association between HLA-B variants and amoxicillin-induced severe cutaneous adverse reactions in Chinese han population.

Amoxicillin (AMX) is among the most prescribed and the best tolerated antimicrobials worldwide. However, it can occasionally trigger severe cutaneous adverse reactions (SCAR) with a significant morbidity and mortality. The genetic factors that may be relevant to AMX-induced SCAR (AMX-SCAR) remain unclear. Identification of the genetic risk factor may prevent patients from the risk of AMX exposure and resume therapy with other falsely implicated drugs. Four patients with AMX-SCAR, 1,000 population control and 100 AMX-tolerant individuals were enrolled in this study. Both exome-wide and HLA-based association studies were conducted. Molecular docking analysis was employed to simulate the interactions between AMX and risk HLA proteins. Compared with AMX-tolerant controls, a significant association of HLA-B*15:01 with AMX-SCAR was validated [odds ratio (OR) = 22.9, 95% confidence interval (CI): 1.68-1275.67; p = 7.34 × 10-3]. Moreover, 75% carriers of HLA-B*15:01 in four patients with AMX-SCAR, and the carrier frequency of 10.7% in 1,000 control individuals and 11.0% in 100 AMX-tolerant controls, respectively. Within HLA-B protein, the S140 present in all cases and demonstrated the strongest association with AMX-SCAR [OR = 53.5, p = 5.18 × 10-4]. Molecular docking results also confirmed the interaction between AMX and S140 of the HLA-B protein, thus eliminating the false-positive results during in association analysis. Our findings suggest that genetic susceptibility may be involved in the development of AMX-SCAR in Han Chinese. However, whether the HLA-B variants observed in this study can be used as an effective genetic marker of AMX-induced SCAR still needs to be further explored in larger cohort studies and other ethnic populations.

Genetic associations of protein-coding variants in venous thromboembolism

Previous genetic studies of venous thromboembolism (VTE) have been largely limited to common variants, leaving the genetic determinants relatively incomplete. We performed an exome-wide association study of VTE among 14,723 cases and 334,315 controls. Fourteen known and four novel genes (SRSF6, PHPT1, CGN, and MAP3K2) were identified through protein-coding variants, with broad replication in the FinnGen cohort. Most genes we discovered exhibited the potential to predict future VTE events in longitudinal analysis. Notably, we provide evidence for the additive contribution of rare coding variants to known genome-wide polygenic risk in shaping VTE risk. The identified genes were enriched in pathways affecting coagulation and platelet activation, along with liver-specific expression. The pleiotropic effects of these genes indicated the potential involvement of coagulation factors, blood cell traits, liver function, and immunometabolic processes in VTE pathogenesis. In conclusion, our study unveils the valuable contribution of protein-coding variants in VTE etiology and sheds new light on its risk stratification.

Divergent role of Mitochondrial Amidoxime Reducing Component 1 (MARC1) in human and mouse.

Recent human genome-wide association studies have identified common missense variants in MARC1, p.Ala165Thr and p.Met187Lys, associated with lower hepatic fat, reduction in liver enzymes and protection from most causes of cirrhosis. Using an exome-wide association study we recapitulated earlier MARC1 p.Ala165Thr and p.Met187Lys findings in 540,000 individuals from five ancestry groups. We also discovered novel rare putative loss of function variants in MARC1 with a phenotype similar to MARC1 p.Ala165Thr/p.Met187Lys variants. In vitro studies of recombinant human MARC1 protein revealed Ala165Thr substitution causes protein instability and aberrant localization in hepatic cells, suggesting MARC1 inhibition or deletion may lead to hepatoprotection. Following this hypothesis, we generated Marc1 knockout mice and evaluated the effect of Marc1 deletion on liver phenotype. Unexpectedly, our study found that whole-body Marc1 deficiency in mouse is not protective against hepatic triglyceride accumulation, liver inflammation or fibrosis. In attempts to explain the lack of the observed phenotype, we discovered that Marc1 plays only a minor role in mouse liver while its paralogue Marc2 is the main Marc family enzyme in mice. Our findings highlight the major difference in MARC1 physiological function between human and mouse.

A cross-trait study of lung cancer and its related respiratory diseases based on large-scale exome sequencing population.

Genome-wide association studies (GWASs) explain the genetic susceptibility between diseases and common variants. Nevertheless, with the appearance of large-scale sequencing profiles, we could explore the rare coding variants in disease pathogenesis. We estimated the genetic correlation of nine respiratory diseases and lung cancer in the UK Biobank (UKB) by linkage disequilibrium score regression (LDSC). Then, we performed exome-wide association studies at single-variant level and gene-level for lung cancer and lung cancer-related respiratory diseases using the whole-exome sequencing (WES) data of 427,934 European participants. Cross-trait meta-analysis was conducted by association analysis based on subsets (ASSET) to identify the pleiotropic variants, while in-silico functional analysis was performed to explore their function. Causal mediation analysis was used to explore whether these pleiotropic variants lead to lung cancer is mediated by affecting the chronic respiratory diseases. Five respiratory diseases [emphysema, pneumonia, asthma, chronic obstructive pulmonary disease (COPD), and fibrosis] were genetically correlated with lung cancer. We identified 102 significant independent variants at single-variant levels for lung cancer and five lung cancer-related diseases. 15:78590583:G>A (missense variant in CHRNA5) was shared in lung cancer, emphysema, and COPD. Meanwhile, 14 significant genes and 87 suggestive genes were identified in gene-based association tests, including HSD3B7 (lung cancer), SRSF2 (pneumonia), TNXB (asthma), TERT (fibrosis), MOSPD3 (emphysema). Based on the cross-trait meta-analysis, we detected 145 independent pleiotropic variants. We further identified abundant pathways with significant enrichment effects, demonstrating that these pleiotropic genes were functional. Meanwhile, the proportion of mediation effects of these variants ranged from 6 to 23 (emphysema: 23%; COPD: 20%; pneumonia: 20%; fibrosis: 7%; asthma: 6%) through these five respiratory diseases to the incidence of lung cancer. The identified shared genetic variants, genes, biological pathways, and potential intermediate causal pathways provide a basis for further exploration of the relationship between lung cancer and respiratory diseases.

Exome sequencing identifies genes associated with sleep-related traits.

Sleep is vital for human health and has a moderate heritability. Previous genome-wide association studies have limitations in capturing the role of rare genetic variants in sleep-related traits. Here we conducted a large-scale exome-wide association study of eight sleep-related traits (sleep duration, insomnia symptoms, chronotype, daytime sleepiness, daytime napping, ease of getting up in the morning, snoring and sleep apnoea) among 450,000 participants from UK Biobank. We identified 22 new genes associated with chronotype (ADGRL4, COL6A3, CLK4 and KRTAP3-3), daytime sleepiness (ST3GAL1 and ANKRD12), daytime napping (PLEKHM1, ANKRD12 and ZBTB21), snoring (WDR59) and sleep apnoea (13 genes). Notably, 20 of these genes were confirmed to be significantly associated with sleep disorders in the FinnGen cohort. Enrichment analysis revealed that these discovered genes were enriched in circadian rhythm and central nervous system neurons. Phenotypic association analysis showed that ANKRD12 was associated with cognition and inflammatory traits. Our results demonstrate the value of large-scale whole-exome analysis in understanding the genetic architecture of sleep-related traits and potential biological mechanisms.

Human APOE4 Protects High-Fat and High-Sucrose Diet Fed Targeted Replacement Mice against Fatty Liver Disease Compared to APOE3.

Recent genome- and exome-wide association studies suggest that the human APOE ε4 allele protects against non-alcoholic fatty liver disease (NAFLD), while ε3 promotes hepatic steatosis and steatohepatitis. The present study aimed at examining the APOE genotype-dependent development of fatty liver disease and its underlying mechanisms in a targeted replacement mouse model. Male mice expressing the human APOE3 or APOE4 protein isoforms on a C57BL/6J background and unmodified C57BL/6J mice were chronically fed a high-fat and high-sucrose diet to induce obesity. After 7 months, body weight gain was more pronounced in human APOE than endogenous APOE expressing mice with elevated plasma biomarkers suggesting aggravated metabolic dysfunction. APOE3 mice exhibited the highest liver weights and, compared to APOE4, massive hepatic steatosis. An untargeted quantitative proteome analysis of the liver identified a high number of proteins differentially abundant in APOE3 versus APOE4 mice. The majority of the higher abundant proteins in APOE3 mice could be grouped to inflammation and damage-associated response, and lipid storage, amongst others. Results of the targeted qRT-PCR and Western blot analyses contribute to the overall finding that APOE3 as opposed to APOE4 promotes hepatic steatosis, inflammatory- and damage-associated response signaling and fibrosis in the liver of obese mice. Our experimental data substantiate the observation of an increased NAFLD-risk associated with the human APOEε3 allele, while APOEε4 appears protective. The underlying mechanisms of the protection possibly involve a higher capacity of nonectopic lipid deposition in subcutaneous adipose tissue and lower hepatic pathogen recognition in the APOE4 mice.

Identification of novel rare variants for anxiety: an exome-wide association study in the UK Biobank

BackgroundRare variants are believed to play a substantial role in the genetic architecture of mental disorders, particularly in coding regions. However, limited evidence supports the impact of rare variants on anxiety. MethodsUsing whole-exome sequencing data from 200,643 participants in the UK Biobank, we investigated the contribution of rare variants to anxiety. Firstly, we computed genetic risk score (GRS) of anxiety utilizing genotype data and summary data from a genome-wide association study (GWAS) on anxiety disorder. Subsequently, we identified individuals within the lowest 50% GRS, a subgroup more likely to carry pathogenic rare variants. Within this subgroup, we classified individuals with the highest 10% 7-item Generalized Anxiety Disorder scale (GAD-7) score as cases (N = 1869), and those with the lowest 10% GAD-7 score were designated as controls (N = 1869). Finally, we conducted gene-based burden tests and single-variant association analyses to assess the relationship between rare variants and anxiety. ResultsTotally, 47,800 variants with MAF ≤0.01 were annotated as non-benign coding variants, consisting of 42,698 nonsynonymous SNVs, 489 nonframeshift substitution, 236 frameshift substitution, 617 stop-gain and 40 stop-loss variants. After variation aggregation, 5066 genes were included in gene-based association analysis. Totally, 11 candidate genes were detected in burden test, such as RNF123 (PBonferroni adjusted = 3.40 × 10−6), MOAP1(PBonferroni adjusted = 4.35 × 10−4), CCDC110 (PBonferroni adjusted = 5.83 × 10−4). Single-variant test detected 9 rare variants, such as rs35726701(RNF123)(PBonferroni adjusted = 3.16 × 10−10) and rs16942615(CAMTA2) (PBonferroni adjusted = 4.04 × 10−4). Notably, RNF123, CCDC110, DNAH2, and CSKMT gene were identified in both tests. ConclusionsOur study identified novel candidate genes for anxiety in protein-coding regions, revealing the contribution of rare variants to anxiety.

The first exome wide association study in Tunisia: identification of candidate loci and pathways with biological relevance for type 2 diabetes.

Type 2 diabetes (T2D) is a multifactorial disease involving genetic and environmental components. Several genome-wide association studies (GWAS) have been conducted to decipher potential genetic aberrations promoting the onset of this metabolic disorder. These GWAS have identified over 400 associated variants, mostly in the intronic or intergenic regions. Recently, a growing number of exome genotyping or exome sequencing experiments have identified coding variants associated with T2D. Such studies were mainly conducted in European populations, and the few candidate-gene replication studies in North African populations revealed inconsistent results. In the present study, we aimed to discover the coding genetic etiology of T2D in the Tunisian population. We carried out a pilot Exome Wide Association Study (EWAS) on 50 Tunisian individuals. Single variant analysis was performed as implemented in PLINK on potentially deleterious coding variants. Subsequently, we applied gene-based and gene-set analyses using MAGMA software to identify genes and pathways associated with T2D. Potential signals were further replicated in an existing large in-silico dataset, involving up to 177116 European individuals. Our analysis revealed, for the first time, promising associations between T2D and variations in MYORG gene, implicated in the skeletal muscle fiber development. Gene-set analysis identified two candidate pathways having nominal associations with T2D in our study samples, namely the positive regulation of neuron apoptotic process and the regulation of mucus secretion. These two pathways are implicated in the neurogenerative alterations and in the inflammatory mechanisms of metabolic diseases. In addition, replication analysis revealed nominal associations of the regulation of beta-cell development and the regulation of peptidase activity pathways with T2D, both in the Tunisian subjects and in the European in-silico dataset. The present study is the first EWAS to investigate the impact of single genetic variants and their aggregate effects on T2D risk in Africa. The promising disease markers, revealed by our pilot EWAS, will promote the understanding of the T2D pathophysiology in North Africa as well as the discovery of potential treatments.

Genetic Factors and Long-term Treatment-Related Neurocognitive Deficits, Anxiety, and Depression in Childhood Leukemia Survivors: An Exome-Wide Association Study.

An increased risk of neurocognitive deficits, anxiety, and depression has been reported in childhood cancer survivors. We analyzed associations of neurocognitive deficits, as well as anxiety and depression, with common and rare genetic variants derived from whole-exome sequencing data of acute lymphoblastic leukemia (ALL) survivors from the PETALE cohort. In addition, significant associations were assessed using stratified and multivariable analyses. Next, top-ranking common associations were analyzed in an independent SJLIFE replication cohort of ALL survivors. Significant associations were identified in the entire discovery cohort (N = 229) between the AK8 gene and changes in neurocognitive function, whereas PTPRZ1, MUC16, TNRC6C-AS1 were associated with anxiety. Following stratification according to sex, the ZNF382 gene was linked to a neurocognitive deficit in males, whereas APOL2 and C6orf165 were associated with anxiety and EXO5 with depression. Following stratification according to prognostic risk groups, the modulatory effect of rare variants on depression was additionally found in the CYP2W1 and PCMTD1 genes. In the replication SJLIFE cohort (N = 688), the male-specific association in the ZNF382 gene was not significant; however, a P value<0.05 was observed when the entire SJLIFE cohort was analyzed. ZNF382 was significant in males in the combined cohorts as shown by meta-analyses as well as the depression-associated gene EXO5. Further research is needed to confirm whether the current findings, along with other known risk factors, may be valuable in identifying patients at increased risk of these long-term complications. Our results suggest that specific genes may be related to increased neuropsychological consequences.