Large-Scale Exome Sequencing Study Implicates Both Developmental and Functional Changes in the Neurobiology of Autism.

We present the largest exome sequencing study of autism spectrum disorder (ASD) to date (n = 35,584 total samples, 11,986 with ASD). Using an enhanced analytical framework to integrate de novo and case-control rare variation, we identify 102 risk genes at a false discovery rate of 0.1 or less. Of these genes, 49 show higher frequencies of disruptive de novo variants in individuals ascertained to have severe neurodevelopmental delay, whereas 53 show higher frequencies in individuals ascertained to have ASD; comparing ASD cases with mutations in these groups reveals phenotypic differences. Expressed early in brain development, most risk genes have roles in regulation of gene expression or neuronal communication (i.e., mutations effect neurodevelopmental and neurophysiological changes), and 13 fall within loci recurrently hit by copy number variants. In cells from the human cortex, expression of risk genes is enriched in excitatory and inhibitory neuronal lineages, consistent with multiple paths to an excitatory-inhibitory imbalance underlying ASD.

Oral and non-oral lichen planus show genetic heterogeneity and differential risk for autoimmune disease and oral cancer.

Lichen planus (LP) is a T-cell-mediated inflammatory disease affecting squamous epithelia in many parts of the body, most often the skin and oral mucosa. Cutaneous LP is usually transient and oral LP (OLP) is most often chronic, so we performed a large-scale genetic and epidemiological study of LP to address whether the oral and non-oral subgroups have shared or distinct underlying pathologies and their overlap with autoimmune disease. Using lifelong records covering diagnoses, procedures, and clinic identity from 473,580 individuals in the FinnGen study, genome-wide association analyses were conducted on carefully constructed subcategories of OLP (n = 3,323) and non-oral LP (n = 4,356) and on the combined group. We identified 15 genome-wide significant associations in FinnGen and an additional 12 when meta-analyzed with UKBB (27 independent associations at 25 distinct genomic locations), most of which are shared between oral and non-oral LP. Many associations coincide with known autoimmune disease loci, consistent with the epidemiologic enrichment of LP with hypothyroidism and other autoimmune diseases. Notably, a third of the FinnGen associations demonstrate significant differences between OLP and non-OLP. We also observed a 13.6-fold risk for tongue cancer and an elevated risk for other oral cancers in OLP, in agreement with earlier reports that connect LP with higher cancer incidence. In addition to a large-scale dissection of LP genetics and comorbidities, our study demonstrates the use of comprehensive, multidimensional health registry data to address outstanding clinical questions and reveal underlying biological mechanisms in common but understudied diseases.

Discordant calls across genotype discovery approaches elucidate variants with systematic errors.

Large-scale high-throughput sequencing data sets have been transformative for informing clinical variant interpretation and for use as reference panels for statistical and population genetic efforts. Although such resources are often treated as ground truth, we find that in widely used reference data sets such as the Genome Aggregation Database (gnomAD), some variants pass gold-standard filters, yet are systematically different in their genotype calls across genotype discovery approaches. The inclusion of such discordant sites in study designs involving multiple genotype discovery strategies could bias results and lead to false-positive hits in association studies owing to technological artifacts rather than a true relationship to the phenotype. Here, we describe this phenomenon of discordant genotype calls across genotype discovery approaches, characterize the error mode of wrong calls, provide a list of discordant sites identified in gnomAD that should be treated with caution in analyses, and present a metric and machine learning classifier trained on gnomAD data to identify likely discordant variants in other data sets. We find that different genotype discovery approaches have different sets of variants at which this problem occurs, but there are characteristic variant features that can be used to predict discordant behavior. Discordant sites are largely shared across ancestry groups, although different populations are powered for the discovery of different variants. We find that the most common error mode is that of a variant being heterozygous for one approach and homozygous for the other, with heterozygous in the genomes and homozygous reference in the exomes making up the majority of miscalls.

Germline cancer susceptibility in individuals with melanoma.

BACKGROUND: Prior studies have estimated a small number of individuals with melanoma (2%-2.5%) have germline cancer predisposition, yet a recent twin study suggested melanoma has the highest hereditability among cancers. OBJECTIVE: To determine the incidence of hereditary melanoma and characterize the spectrum of cancer predisposition genes that may increase the risk of melanoma. METHODS: Four hundred individuals with melanoma and personal or family history of cancers underwent germline testing of >80 cancer predisposition genes. Comparative analysis of germline data was performed on 3 additional oncologic and dermatologic data sets. RESULTS: Germline pathogenic/likely pathogenic (P/LP) variants were identified in 15.3% (61) individuals with melanoma. Most variants (41, 67%) involved genes considered unrelated to melanoma (BLM, BRIP1, CHEK2, MLH1, MSH2, PMS2, RAD51C). A third (20, 33%) were in genes previously associated with familial melanoma (BAP1, BRCA2, CDKN2A, MITF, TP53). Nearly half (30, 46.9%) of P/LP variants were in homologous repair deficiency genes. Validation cohorts demonstrated P/LP rates of 10.6% from an unselected oncologic cohort, 15.8% from a selected commercial testing cohort, and 14.5% from a highly selected dermatologic study. LIMITATIONS: Cohorts with varying degrees of selection, some retrospective. CONCLUSION: Germline predisposition in individuals with melanoma is common, with clinically actionable findings diagnosed in 10.6% to 15.8%.

Haplotype Sharing Provides Insights into Fine-Scale Population History and Disease in Finland.

Finland provides unique opportunities to investigate population and medical genomics because of its adoption of unified national electronic health records, detailed historical and birth records, and serial population bottlenecks. We assembled a comprehensive view of recent population history (≤100 generations), the timespan during which most rare-disease-causing alleles arose, by comparing pairwise haplotype sharing from 43,254 Finns to that of 16,060 Swedes, Estonians, Russians, and Hungarians from geographically and linguistically adjacent countries with different population histories. We find much more extensive sharing in Finns, with at least one ≥ 5 cM tract on average between pairs of unrelated individuals. By coupling haplotype sharing with fine-scale birth records from more than 25,000 individuals, we find that although haplotype sharing broadly decays with geographical distance, there are pockets of excess haplotype sharing; individuals from northeast Finland typically share several-fold more of their genome in identity-by-descent segments than individuals from southwest regions. We estimate recent effective population-size changes through time across regions of Finland, and we find that there was more continuous gene flow as Finns migrated from southwest to northeast between the early- and late-settlement regions than was dichotomously described previously. Lastly, we show that haplotype sharing is locally enriched by an order of magnitude among pairs of individuals sharing rare alleles and especially among pairs sharing rare disease-causing variants. Our work provides a general framework for using haplotype sharing to reconstruct an integrative view of recent population history and gain insight into the evolutionary origins of rare variants contributing to disease.

Targeted gene sequencing in 6994 individuals with neurodevelopmental disorder with epilepsy.

PURPOSE: We aimed to gain insight into frequencies of genetic variants in genes implicated in neurodevelopmental disorder with epilepsy (NDD+E) by investigating large cohorts of patients in a diagnostic setting. METHODS: We analyzed variants in NDD+E using epilepsy gene panel sequencing performed between 2013 and 2017 by two large diagnostic companies. We compared variant frequencies in 6994 panels with another 8588 recently published panels as well as exome-wide de novo variants in 1942 individuals with NDD+E and 10,937 controls. RESULTS: Genes with highest frequencies of ultrarare variants in NDD+E comprised SCN1A, KCNQ2, SCN2A, CDKL5, SCN8A, and STXBP1, concordant with the two other epilepsy cohorts we investigated. In only 46% of the analyzed 262 dominant and X-linked panel genes ultrarare variants in patients were reported. Among genes with contradictory evidence of association with epilepsy, CACNB4, CLCN2, EFHC1, GABRD, MAGI2, and SRPX2 showed equal frequencies in cases and controls. CONCLUSION: We show that improvement of panel design increased diagnostic yield over time, but panels still display genes with low or no diagnostic yield. With our data, we hope to improve current diagnostic NDD+E panel design and provide a resource of ultrarare variants in individuals with NDD+E to the community.

Public platform with 39,472 exome control samples enables association studies without genotype sharing.

Acquiring a sufficiently powered cohort of control samples matched to a case sample can be time-consuming or, in some cases, impossible. Accordingly, an ability to leverage genetic data from control samples that were already collected elsewhere could dramatically improve power in genetic association studies. Sharing of control samples can pose significant challenges, since most human genetic data are subject to strict sharing regulations. Here, using the properties of singular value decomposition and subsampling algorithm, we developed a method allowing selection of the best-matching controls in an external pool of samples compliant with personal data protection and eliminating the need for genotype sharing. We provide access to a library of 39,472 exome sequencing controls at http://dnascore.net enabling association studies for case cohorts lacking control subjects. Using this approach, control sets can be selected from this online library with a prespecified matching accuracy, ensuring well-calibrated association analysis for both rare and common variants.

BRAF V600E in a preclinical model of pleomorphic Xanthoastrocytoma: Analysis of the tumor microenvironment and immune cell infiltration dynamics in vivo.

Low-grade glioma (LGG) is the most common brain tumor affecting pediatric patients (pLGG) and BRAF mutations constitute the most frequent genetic alterations. Within the spectrum of pLGGs, approximately 70%-80% of pediatric patients diagnosed with transforming pleomorphic xanthoastrocytoma (PXA) harbor the BRAF V600E mutation. However, the impact of glioma BRAF V600E cell regulation of tumor-infiltrating immune cells and their contribution to tumor progression remains unclear. Moreover, the efficacy of BRAF inhibitors in treating pLGGs is limited compared with their impact on BRAF-mutated melanoma. Here we report a novel immunocompetent RCAS-BRAF V600E murine glioma model. Pathological assessment indicates this model seems to be consistent with diffuse gliomas and morphological features of PXA. Our investigations revealed distinct immune cell signatures associated with increased trafficking and activation within the tumor microenvironment (TME). Intriguingly, immune system activation within the TME also generated a pronounced inflammatory response associated with dysfunctional CD8+ T cells, increased presence of immunosuppressive myeloid cells and regulatory T cells. Further, our data suggests tumor-induced inflammatory processes, such as cytokine storm. These findings suggest a complex interplay between tumor progression and the robust inflammatory response within the TME in preclinical BRAF V600E LGGs, which may significantly influence animal survival.

Whole-exome sequencing reveals an association of rs112065068 in TGOLN2 gene with distant metastasis of non-small cell lung cancer.

Early prediction and prevention of recurring illness is critical for improving the survival rates of patients with non-small cell lung cancer (NSCLC). Previously, we demonstrated that the presence of premalignant epithelial changes in the small bronchi distant to the primary tumor is associated with NSCLC progression: isolated basal cell hyperplasia (iBCH) indicates a high risk of distant metastasis, BCH combined with squamous metaplasia (BCHSM) - a high risk of locoregional recurrence. Here, we aimed to identify germline single nucleotide variants (SNVs) and insertions and deletions (InDels) associated with distant metastasis and locoregional recurrence in cases with iBCH and BCHSM using whole-exome sequencing of 172 NSCLC patients. The rs112065068 of the TGOLN2 gene was identified only in iBCH patients and was associated with a high risk of distant metastasis (P < .001) and worse metastasis-free survival (HR = 4.19 (95 %CI 1.97-8.93); P < .001). This variant was validated in a group of 109 NSCLC patients using real-time PCR and Sanger sequencing analyses. To our knowledge, this study is the first to identify a germline variant associated with NSCLC distant metastasis.

Autoimmune hypothyroidism GWAS reveals independent autoimmune and thyroid-specific contributions and an inverse relation with cancer risk.

The high prevalence of autoimmune hypothyroidism (AIHT) - more than 5% in human populations - provides a unique opportunity to unlock the most complete picture to date of genetic loci that underlie systemic and organ-specific autoimmunity. Using a meta-analysis of 81,718 AIHT cases in FinnGen and the UK Biobank, we dissect associations along axes of thyroid dysfunction and autoimmunity. This largest-to-date scan of hypothyroidism identifies 418 independent associations (p < 5x10 - 8 ), more than half of which have not previously been documented in thyroid disease. In 48 of these, a protein-coding variant is the lead SNP or is highly correlated (r 2 > 0.95) with the lead SNP at the locus, including low-frequency coding variants at LAG3 , ZAP70 , TG, TNFSF11, IRF3, S1PR4, HABP2, ZNF429 as well as established variants at ADCY7, IFIH1 and TYK2 . The variants at LAG3 (P67T), ZAP70 (T155M), and TG (Q655X) are highly enriched in Finland and functional experiments in T-cells demonstrate that the ZAP70 :T155M allele reduces T-cell activation. By employing a large-scale scan of non-thyroid autoimmunity and a published meta-analysis of TSH levels, we use a Bayesian classifier to dissect the associated loci into distinct groupings and from this estimate, a significant proportion are involved in systemic (i.e., general to multiple autoimmune conditions) autoimmunity (34%) and another subset in thyroid-specific dysfunction (17%). By comparing these association results further to other common disease endpoints, we identify a noteworthy overlap with skin cancer, with 10% of AIHT loci showing a consistent but opposite pattern of association where alleles that increase the risk of hypothyroidism have protective effects for skin cancer. The association results, including genes encoding checkpoint inhibitors and other genes affecting protein levels of PD1, bolster the causal role of natural variation in autoimmunity influencing cancer outcomes.

Homeostatic iron regulatory protein drives glioblastoma growth via tumor cell-intrinsic and sex-specific responses.

Background: Glioblastoma (GBM) displays alterations in iron that drive proliferation and tumor growth. Iron regulation is complex and involves many regulatory mechanisms, including the homeostatic iron regulator (HFE) gene, which encodes the homeostatic iron regulatory protein. While HFE is upregulated in GBM and correlates with poor survival outcomes, the function of HFE in GBM remains unclear. Methods: We interrogated the impact of cell-intrinsic Hfe expression on proliferation and survival of intracranially implanted animals through genetic gain- and loss-of-function approaches in syngeneic mouse glioma models, along with in vivo immune assessments. We also determined the expression of iron-associated genes and their relationship to survival in GBM using public data sets and used transcriptional profiling to identify differentially expressed pathways in control compared to Hfe-knockdown cells. Results: Overexpression of Hfe accelerated GBM proliferation and reduced animal survival, whereas suppression of Hfe induced apoptotic cell death and extended survival, which was more pronounced in females and associated with attenuation of natural killer cells and CD8+ T cell activity. Analysis of iron gene signatures in Hfe-knockdown cells revealed alterations in the expression of several iron-associated genes, suggesting global disruption of intracellular iron homeostasis. Further analysis of differentially expressed pathways revealed oxidative stress as the top pathway upregulated following Hfe loss. Hfe knockdown indeed resulted in enhanced 55Fe uptake and generation of reactive oxygen species. Conclusions: These findings reveal an essential function for HFE in GBM cell growth and survival, as well as a sex-specific interaction with the immune response.

Detecting biased validation of predictive models in the positive-unlabeled setting: disease gene prioritization case study.

Motivation: Positive-unlabeled data consists of points with either positive or unknown labels. It is widespread in medical, genetic, and biological settings, creating a high demand for predictive positive-unlabeled models. The performance of such models is usually estimated using validation sets, assumed to be selected completely at random (SCAR) from known positive examples. For certain metrics, this assumption enables unbiased performance estimation when treating positive-unlabeled data as positive/negative. However, the SCAR assumption is often adopted without proper justifications, simply for the sake of convenience. Results: We provide an algorithm that under the weak assumptions of a lower bound on the number of positive examples can test for the violation of the SCAR assumption. Applying it to the problem of gene prioritization for complex genetic traits, we illustrate that the SCAR assumption is often violated there, causing the inflation of performance estimates, which we refer to as validation bias. We estimate the potential impact of validation bias on performance estimation. Our analysis reveals that validation bias is widespread in gene prioritization data and can significantly overestimate the performance of models. This finding elucidates the discrepancy between the reported good performance of models and their limited practical applications. Availability and implementation: Python code with examples of application of the validation bias detection algorithm is available at github.com/ArtomovLab/ValidationBias.

Gene, cell type, and drug prioritization analysis suggest genetic basis for the utility of diuretics in treating Alzheimer disease.

We introduce a user-friendly tool for risk gene, cell type, and drug prioritization for complex traits: GCDPipe. It uses gene-level GWAS-derived data and gene expression data to train a model for the identification of disease risk genes and relevant cell types. Gene prioritization information is then coupled with known drug target data to search for applicable drug agents based on their estimated functional effects on the identified risk genes. We illustrate the utility of our approach in different settings: identification of the cell types, implicated in disease pathogenesis, was tested in inflammatory bowel disease (IBD) and Alzheimer disease (AD); gene target and drug prioritization was tested in IBD and schizophrenia. The analysis of phenotypes with known disease-affected cell types and/or existing drug candidates shows that GCDPipe is an effective tool to unify genetic risk factors with cellular context and known drug targets. Next, analysis of the AD data with GCDPipe suggested that gene targets of diuretics, as an Anatomical Therapeutic Chemical drug subgroup, are significantly enriched among the genes prioritized by GCDPipe, indicating their possible effect on the course of the disease.

Mendelian Randomization Analysis reveals Inverse Genetic Risks between Skin Cancers and Vitiligo.

Several observational studies have demonstrated a consistent pattern of decreased melanoma risk among patients with vitiligo. More recently, this finding has been supported by a suggested genetic relationship between the two entities, with certain variants significantly associated with an increased risk of melanoma, basal cell carcinoma, and squamous cell carcinoma but a decreased risk of vitiligo. We compared 48 associated variants from a recently published GWAS and identified three variants-located in the TYR, MC1R-DEF8, and RALY-EIF2S2-ASIP-AHCY-ITCH loci- that correlated with an increased risk for melanoma, basal cell carcinoma, and squamous cell carcinoma and a decreased risk for vitiligo. We then used results of skin cancers and vitiligo GWAS to compare the shared genetic properties between these two traits through an unbiased Mendelian randomization analysis. Our results suggest that the inverse genetic relationship between common skin cancers and vitiligo is broader than previously reported owing to the influence of shared genome-wide significant associations.

Transgenerational and intergenerational effects of early childhood famine exposure in the cohort of offspring of Leningrad Siege survivors.

Famine exposure during early life development can affect disease risk in late-life period, yet, transmission of phenotypic features from famine-exposed individuals to the next generations has not been well characterized. The purpose of our case-control study was to investigate the association of parental starvation in the perinatal period and the period of early childhood with the phenotypic features observed in two generations of descendants of Leningrad siege survivors. We examined 54 children and 30 grandchildren of 58 besieged Leningrad residents who suffered from starvation in early childhood and prenatal age during the Second World War. Controls from the population-based national epidemiological ESSE-RF study (n = 175) were matched on sex, age and body mass index (BMI). Phenotypes of controls and descendants (both generations, children and grandchildren separately) were compared, taking into account multiple testing. Comparison of two generations descendants with corresponding control groups revealed significantly higher creatinine and lower glomerular filtration rate (GFR), both in meta-analysis and in independent analyses. The mean values of GFR for all groups were within the normal range (GFR less than 60 mL/min/1.73 m2 was recorded in 2 controls and no one in DLSS). Additionally, independent of the creatinine level, differences in the eating pattern were detected: insufficient fish and excessive red meat consumption were significantly more frequent in the children of the Leningrad siege survivors compared with controls. Blood pressure, blood lipids and glucose did not differ between the groups. Parental famine exposure in early childhood may contribute to a decrease in kidney filtration capacity and altered eating pattern in the offspring of famine-exposed individuals.

CR1 variants contribute to FSGS susceptibility across multiple populations.

Focal segmental glomerulosclerosis (FSGS) is a common cause of nephrotic syndrome with an annual incidence in the United States in African-Americans compared to European-Americans of 24 cases and 5 cases per million, respectively. Among glomerular diseases in Europe and Latin-America, FSGS was the second most frequent diagnosis, and in Asia the fifth. We expand previous efforts in understanding genetics of FSGS by performing a case-control study involving ethnically-diverse groups FSGS cases (726) and a pool of controls (13,994), using panel sequencing of approximately 2,500 podocyte-expressed genes. Through rare variant association tests, we replicated known risk genes - KANK1, COL4A4, and APOL1. A novel significant association was observed for the gene encoding complement receptor 1 (CR1). High-risk rare variants in CR1 in the European-American cohort were commonly observed in Latin- and African-Americans. Therefore, a combined rare and common variant analysis was used to replicate the CR1 association in non-European populations. The CR1 risk variant, rs17047661, gives rise to the Sl1/Sl2 (R1601G) allele that was previously associated with protection against cerebral malaria. Pleiotropic effects of rs17047661 may explain the difference in allele frequencies across continental ancestries and suggest a possible role for genetically-driven alterations of adaptive immunity in the pathogenesis of FSGS.

Case report: Somatic mutations in microtubule dynamics-associated genes in patients with WNT-medulloblastoma tumors.

Medulloblastoma (MB) is the most common pediatric brain tumor which accounts for about 20% of all pediatric brain tumors and 63% of intracranial embryonal tumors. MB is considered to arise from precursor cell populations present during an early brain development. Most cases (~70%) of MB occur at the age of 1-4 and 5-9, but are also infrequently found in adults. Total annual frequency of pediatric tumors is about 5 cases per 1 million children. WNT-subtype of MB is characterized by a high probability of remission, with a long-term survival rate of about 90%. However, in some rare cases there may be increased metastatic activity, which dramatically reduces the likelihood of a favorable outcome. Here we report two cases of MB with a histological pattern consistent with desmoplastic/nodular (DP) and classic MB, and genetically classified as WNT-MB. Both cases showed putative causal somatic protein truncating mutations identified in microtubule-associated genes: ARID2, TUBB4A, and ANK3.

Genotype imputation and polygenic score estimation in northwestern Russian population.

Numerous studies demonstrated the lack of transferability of polygenic score (PGS) models across populations and the problem arising from unequal presentation of ancestries across genetic studies. However, even within European ancestry there are ethnic groups that are rarely presented in genetic studies. For instance, Russians, being one of the largest, diverse, and yet understudied group in Europe. In this study, we evaluated the reliability of genotype imputation for the Russian cohort by testing several commonly used imputation reference panels (e.g. HRC, 1000G, HGDP). HRC, in comparison with two other panels, showed the most accurate results based on both imputation accuracy and allele frequency concordance between masked and imputed genotypes. We built polygenic score models based on GWAS results from the UK biobank, measured the explained phenotypic variance in the Russian cohort attributed to polygenic scores for 11 phenotypes, collected in the clinic for each participant, and finally explored the role of allele frequency discordance between the UK biobank and the study cohort in the resulting PGS performance.

Case Report: Supernormal Vascular Aging in Leningrad Siege Survivors.

Age-related changes in the vascular system play an important role in the biological age and lifespan of a person and maybe affected from an early age onward. One of the indicators of changes in the vascular system is arterial wall stiffness and its main measure, i.e., carotid-femoral pulse wave velocity (cfPWV). We examined arterial wall stiffness in a sample of 305 Leningrad Siege survivors to assess how hunger and stressful conditions during fetal development and early childhood affected the state of the cardiovascular system at a later age and what factors may neutralize the negative impact sustained in early childhood. Here, we presented an evaluation of two unique patients with supernormal vascular aging (SUPERNOVA) phenotype from this cohort and described the details of congruence between hereditary resistance and practiced lifestyle yielding slower biological aging rate.

GWAS of depression in 4,520 individuals from the Russian population highlights the role of MAGI2 (S-SCAM) in the gut-brain axis.

We present the results of the depression Genome-wide association studies study performed on a cohort of Russian-descent individuals, which identified a novel association at chromosome 7q21 locus. Gene prioritization analysis based on already known depression risk genes indicated MAGI2 (S-SCAM) as the most probable gene from the locus and potential susceptibility gene for the disease. Brain and gut expression patterns were the main features highlighting functional relatedness of MAGI2 to the previously known depression risk genes. Local genetic covariance analysis, analysis of gene expression, provided initial suggestive evidence of hospital anxiety and depression scale and diagnostic and statistical manual of mental disorders scales having a different relationship with gut-brain axis disturbance. It should be noted, that while several independent methods successfully in silico validate the role of MAGI2, we were unable to replicate genetic association for the leading variant in the MAGI2 locus, therefore the role of rs521851 in depression should be interpreted with caution.