De novo mutations revealed by whole-exome sequencing are strongly associated with autism

Multiple studies have confirmed the contribution of rare de novo copy number variations to the risk for autism spectrum disorders. But whereas de novo single nucleotide variants have been identified in affected individuals, their contribution to risk has yet to be clarified. Specifically, the frequency and distribution of these mutations have not been well characterized in matched unaffected controls, and such data are vital to the interpretation of de novo coding mutations observed in probands. Here we show, using whole-exome sequencing of 928 individuals, including 200 phenotypically discordant sibling pairs, that highly disruptive (nonsense and splice-site) de novo mutations in brain-expressed genes are associated with autism spectrum disorders and carry large effects. On the basis of mutation rates in unaffected individuals, we demonstrate that multiple independent de novo single nucleotide variants in the same gene among unrelated probands reliably identifies risk alleles, providing a clear path forward for gene discovery. Among a total of 279 identified de novo coding mutations, there is a single instance in probands, and none in siblings, in which two independent nonsense variants disrupt the same gene, SCN2A (sodium channel, voltage-gated, type II, α subunit), a result that is highly unlikely by chance.

De novo mutations in histone-modifying genes in congenital heart disease

Congenital heart disease (CHD) is the most frequent birth defect, affecting 0.8% of live births. Many cases occur sporadically and impair reproductive fitness, suggesting a role for de novo mutations. Here we compare the incidence of de novo mutations in 362 severe CHD cases and 264 controls by analysing exome sequencing of parent-offspring trios. CHD cases show a significant excess of protein-altering de novo mutations in genes expressed in the developing heart, with an odds ratio of 7.5 for damaging (premature termination, frameshift, splice site) mutations. Similar odds ratios are seen across the main classes of severe CHD. We find a marked excess of de novo mutations in genes involved in the production, removal or reading of histone 3 lysine 4 (H3K4) methylation, or ubiquitination of H2BK120, which is required for H3K4 methylation. There are also two de novo mutations in SMAD2, which regulates H3K27 methylation in the embryonic left-right organizer. The combination of both activating (H3K4 methylation) and inactivating (H3K27 methylation) chromatin marks characterizes 'poised' promoters and enhancers, which regulate expression of key developmental genes. These findings implicate de novo point mutations in several hundreds of genes that collectively contribute to approximately 10% of severe CHD.

Genomic analysis of non-NF2 meningiomas reveals mutations in TRAF7, KLF4, AKT1, and SMO

We report genomic analysis of 300 meningiomas, the most common primary brain tumors, leading to the discovery of mutations in TRAF7, a proapoptotic E3 ubiquitin ligase, in nearly one-fourth of all meningiomas. Mutations in TRAF7 commonly occurred with a recurrent mutation (K409Q) in KLF4, a transcription factor known for its role in inducing pluripotency, or with AKT1(E17K), a mutation known to activate the PI3K pathway. SMO mutations, which activate Hedgehog signaling, were identified in ~5% of non-NF2 mutant meningiomas. These non-NF2 meningiomas were clinically distinctive-nearly always benign, with chromosomal stability, and originating from the medial skull base. In contrast, meningiomas with mutant NF2 and/or chromosome 22 loss were more likely to be atypical, showing genomic instability, and localizing to the cerebral and cerebellar hemispheres. Collectively, these findings identify distinct meningioma subtypes, suggesting avenues for targeted therapeutics.

Genetic and Pharmacologic Inactivation of ANGPTL3 and Cardiovascular Disease

BACKGROUND

Loss-of-function variants in the angiopoietin-like 3 gene (ANGPTL3) have been associated with decreased plasma levels of triglycerides, low-density lipoprotein (LDL) cholesterol, and high-density lipoprotein (HDL) cholesterol. It is not known whether such variants or therapeutic antagonism of ANGPTL3 are associated with a reduced risk of atherosclerotic cardiovascular disease.

METHODS

We sequenced the exons of ANGPTL3 in 58,335 participants in the DiscovEHR human genetics study. We performed tests of association for loss-of-function variants in ANGPTL3 with lipid levels and with coronary artery disease in 13,102 case patients and 40,430 controls from the DiscovEHR study, with follow-up studies involving 23,317 case patients and 107,166 controls from four population studies. We also tested the effects of a human monoclonal antibody, evinacumab, against Angptl3 in dyslipidemic mice and against ANGPTL3 in healthy human volunteers with elevated levels of triglycerides or LDL cholesterol.

RESULTS

In the DiscovEHR study, participants with heterozygous loss-of-function variants in ANGPTL3 had significantly lower serum levels of triglycerides, HDL cholesterol, and LDL cholesterol than participants without these variants. Loss-of-function variants were found in 0.33% of case patients with coronary artery disease and in 0.45% of controls (adjusted odds ratio, 0.59; 95% confidence interval, 0.41 to 0.85; P=0.004). These results were confirmed in the follow-up studies. In dyslipidemic mice, inhibition of Angptl3 with evinacumab resulted in a greater decrease in atherosclerotic lesion area and necrotic content than a control antibody. In humans, evinacumab caused a dose-dependent placebo-adjusted reduction in fasting triglyceride levels of up to 76% and LDL cholesterol levels of up to 23%.

CONCLUSIONS

Genetic and therapeutic antagonism of ANGPTL3 in humans and of Angptl3 in mice was associated with decreased levels of all three major lipid fractions and decreased odds of atherosclerotic cardiovascular disease. (Funded by Regeneron Pharmaceuticals and others; ClinicalTrials.gov number, NCT01749878 .).

A Protein-Truncating HSD17B13 Variant and Protection from Chronic Liver Disease

BACKGROUND

Elucidation of the genetic factors underlying chronic liver disease may reveal new therapeutic targets.

METHODS

We used exome sequence data and electronic health records from 46,544 participants in the DiscovEHR human genetics study to identify genetic variants associated with serum levels of alanine aminotransferase (ALT) and aspartate aminotransferase (AST). Variants that were replicated in three additional cohorts (12,527 persons) were evaluated for association with clinical diagnoses of chronic liver disease in DiscovEHR study participants and two independent cohorts (total of 37,173 persons) and with histopathological severity of liver disease in 2391 human liver samples.

RESULTS

A splice variant (rs72613567:TA) in HSD17B13, encoding the hepatic lipid droplet protein hydroxysteroid 17-beta dehydrogenase 13, was associated with reduced levels of ALT (P=4.2×10^-12) and AST (P=6.2×10^-10). Among DiscovEHR study participants, this variant was associated with a reduced risk of alcoholic liver disease (by 42% [95% confidence interval {CI}, 20 to 58] among heterozygotes and by 53% [95% CI, 3 to 77] among homozygotes), nonalcoholic liver disease (by 17% [95% CI, 8 to 25] among heterozygotes and by 30% [95% CI, 13 to 43] among homozygotes), alcoholic cirrhosis (by 42% [95% CI, 14 to 61] among heterozygotes and by 73% [95% CI, 15 to 91] among homozygotes), and nonalcoholic cirrhosis (by 26% [95% CI, 7 to 40] among heterozygotes and by 49% [95% CI, 15 to 69] among homozygotes). Associations were confirmed in two independent cohorts. The rs72613567:TA variant was associated with a reduced risk of nonalcoholic steatohepatitis, but not steatosis, in human liver samples. The rs72613567:TA variant mitigated liver injury associated with the risk-increasing PNPLA3 p.I148M allele and resulted in an unstable and truncated protein with reduced enzymatic activity.

CONCLUSIONS

A loss-of-function variant in HSD17B13 was associated with a reduced risk of chronic liver disease and of progression from steatosis to steatohepatitis. (Funded by Regeneron Pharmaceuticals and others.).

Exome sequencing and analysis of 454,787 UK Biobank participants

A major goal in human genetics is to use natural variation to understand the phenotypic consequences of altering each protein-coding gene in the genome. Here we used exome sequencing¹ to explore protein-altering variants and their consequences in 454,787 participants in the UK Biobank study². We identified 12 million coding variants, including around 1 million loss-of-function and around 1.8 million deleterious missense variants. When these were tested for association with 3,994 health-related traits, we found 564 genes with trait associations at P ≤ 2.18 × 10^-11. Rare variant associations were enriched in loci from genome-wide association studies (GWAS), but most (91%) were independent of common variant signals. We discovered several risk-increasing associations with traits related to liver disease, eye disease and cancer, among others, as well as risk-lowering associations for hypertension (SLC9A3R2), diabetes (MAP3K15, FAM234A) and asthma (SLC27A3). Six genes were associated with brain imaging phenotypes, including two involved in neural development (GBE1, PLD1). Of the signals available and powered for replication in an independent cohort, 81% were confirmed; furthermore, association signals were generally consistent across individuals of European, Asian and African ancestry. We illustrate the ability of exome sequencing to identify gene-trait associations, elucidate gene function and pinpoint effector genes that underlie GWAS signals at scale.

Somatic and germline CACNA1D calcium channel mutations in aldosterone-producing adenomas and primary aldosteronism

Adrenal aldosterone-producing adenomas (APAs) constitutively produce the salt-retaining hormone aldosterone and are a common cause of severe hypertension. Recurrent mutations in the potassium channel KCNJ5 that result in cell depolarization and Ca²⁺ influx cause ~40% of these tumors¹. We found five somatic mutations (four altering glycine 403, one altering isoleucine 770) in CACNA1D, encoding a voltage-gated calcium channel, among 43 non-KCNJ5-mutant APAs. These mutations lie in S6 segments that line the channel pore. Both result in channel activation at less depolarized potentials, and glycine 403 mutations also impair channel inactivation. These effects are inferred to cause increased Ca²⁺ influx, the sufficient stimulus for aldosterone production and cell proliferation in adrenal glomerulosa². Remarkably, we identified de novo mutations at the identical positions in two children with a previously undescribed syndrome featuring primary aldosteronism and neuromuscular abnormalities. These findings implicate gain of function Ca²⁺ channel mutations in aldosterone-producing adenomas and primary aldosteronism.

Distribution and clinical impact of functional variants in 50,726 whole-exome sequences from the DiscovEHR study

The DiscovEHR collaboration between the Regeneron Genetics Center and Geisinger Health System couples high-throughput sequencing to an integrated health care system using longitudinal electronic health records (EHRs). We sequenced the exomes of 50,726 adult participants in the DiscovEHR study to identify ~4.2 million rare single-nucleotide variants and insertion/deletion events, of which ~176,000 are predicted to result in a loss of gene function. Linking these data to EHR-derived clinical phenotypes, we find clinical associations supporting therapeutic targets, including genes encoding drug targets for lipid lowering, and identify previously unidentified rare alleles associated with lipid levels and other blood level traits. About 3.5% of individuals harbor deleterious variants in 76 clinically actionable genes. The DiscovEHR data set provides a blueprint for large-scale precision medicine initiatives and genomics-guided therapeutic discovery.

Inactivating Variants in ANGPTL4 and Risk of Coronary Artery Disease

BACKGROUND

Higher-than-normal levels of circulating triglycerides are a risk factor for ischemic cardiovascular disease. Activation of lipoprotein lipase, an enzyme that is inhibited by angiopoietin-like 4 (ANGPTL4), has been shown to reduce levels of circulating triglycerides.

METHODS

We sequenced the exons of ANGPTL4 in samples obtain from 42,930 participants of predominantly European ancestry in the DiscovEHR human genetics study. We performed tests of association between lipid levels and the missense E40K variant (which has been associated with reduced plasma triglyceride levels) and other inactivating mutations. We then tested for associations between coronary artery disease and the E40K variant and other inactivating mutations in 10,552 participants with coronary artery disease and 29,223 controls. We also tested the effect of a human monoclonal antibody against ANGPTL4 on lipid levels in mice and monkeys.

RESULTS

We identified 1661 heterozygotes and 17 homozygotes for the E40K variant and 75 participants who had 13 other monoallelic inactivating mutations in ANGPTL4. The levels of triglycerides were 13% lower and the levels of high-density lipoprotein (HDL) cholesterol were 7% higher among carriers of the E40K variant than among noncarriers. Carriers of the E40K variant were also significantly less likely than noncarriers to have coronary artery disease (odds ratio, 0.81; 95% confidence interval, 0.70 to 0.92; P=0.002). K40 homozygotes had markedly lower levels of triglycerides and higher levels of HDL cholesterol than did heterozygotes. Carriers of other inactivating mutations also had lower triglyceride levels and higher HDL cholesterol levels and were less likely to have coronary artery disease than were noncarriers. Monoclonal antibody inhibition of Angptl4 in mice and monkeys reduced triglyceride levels.

CONCLUSIONS

Carriers of E40K and other inactivating mutations in ANGPTL4 had lower levels of triglycerides and a lower risk of coronary artery disease than did noncarriers. The inhibition of Angptl4 in mice and monkeys also resulted in corresponding reductions in these values. (Funded by Regeneron Pharmaceuticals.).

Mutation of NLRC4 causes a syndrome of enterocolitis and autoinflammation

Upon detection of pathogen-associated molecular patterns, innate immune receptors initiate inflammatory responses. These receptors include cytoplasmic NOD-like receptors (NLRs), whose stimulation recruits and proteolytically activates caspase-1 within the inflammasome, a multi-protein complex. Caspase-1 mediates the production of interleukin-1 family cytokines (IL1FCs), leading to fever, and inflammatory cell death (pyroptosis)^1,2. Mutations that constitutively activate these pathways underlie several autoinflammatory diseases with diverse clinical features³. We describe a family with a previously unreported syndrome featuring neonatal-onset enterocolitis, periodic fever, and fatal/near-fatal episodes of autoinflammation caused by a de novo gain-of-function mutation (p.V341A) in the HD1 domain of NLRC4 that co-segregates with disease. Mutant NLRC4 causes constitutive Interleukin-1 family cytokine production and macrophage cell death. Infected patient macrophages are polarized toward pyroptosis and exhibit abnormal staining for inflammasome components. These findings describe and reveal the cause of a life-threatening but treatable autoinflammatory disease that underscores the divergent roles of the NLRC4 inflammasome.

Exome sequencing and characterization of 49,960 individuals in the UK Biobank

The UK Biobank is a prospective study of 502,543 individuals, combining extensive phenotypic and genotypic data with streamlined access for researchers around the world1. Here we describe the release of exome-sequence data for the first 49,960 study participants, revealing approximately 4 million coding variants (of which around 98.6% have a frequency of less than 1%). The data include 198,269 autosomal predicted loss-of-function (LOF) variants, a more than 14-fold increase compared to the imputed sequence. Nearly all genes (more than 97%) had at least one carrier with a LOF variant, and most genes (more than 69%) had at least ten carriers with a LOF variant. We illustrate the power of characterizing LOF variants in this population through association analyses across 1,730 phenotypes. In addition to replicating established associations, we found novel LOF variants with large effects on disease traits, including PIEZO1 on varicose veins, COL6A1 on corneal resistance, MEPE on bone density, and IQGAP2 and GMPR on blood cell traits. We further demonstrate the value of exome sequencing by surveying the prevalence of pathogenic variants of clinical importance, and show that 2% of this population has a medically actionable variant. Furthermore, we characterize the penetrance of cancer in carriers of pathogenic BRCA1 and BRCA2 variants. Exome sequences from the first 49,960 participants highlight the promise of genome sequencing in large population-based studies and are now accessible to the scientific community.

Genome-wide association study of alcohol consumption and use disorder in 274,424 individuals from multiple populations

Alcohol consumption level and alcohol use disorder (AUD) diagnosis are moderately heritable traits. We conduct genome-wide association studies of these traits using longitudinal Alcohol Use Disorder Identification Test-Consumption (AUDIT-C) scores and AUD diagnoses in a multi-ancestry Million Veteran Program sample (N = 274,424). We identify 18 genome-wide significant loci: 5 associated with both traits, 8 associated with AUDIT-C only, and 5 associated with AUD diagnosis only. Polygenic Risk Scores (PRS) for both traits are associated with alcohol-related disorders in two independent samples. Although a significant genetic correlation reflects the overlap between the traits, genetic correlations for 188 non-alcohol-related traits differ significantly for the two traits, as do the phenotypes associated with the traits' PRS. Cell type group partitioning heritability enrichment analyses also differentiate the two traits. We conclude that, although heavy drinking is a key risk factor for AUD, it is not a sufficient cause of the disorder.

Within-sibship genome-wide association analyses decrease bias in estimates of direct genetic effects

Estimates from genome-wide association studies (GWAS) of unrelated individuals capture effects of inherited variation (direct effects), demography (population stratification, assortative mating) and relatives (indirect genetic effects). Family-based GWAS designs can control for demographic and indirect genetic effects, but large-scale family datasets have been lacking. We combined data from 178,086 siblings from 19 cohorts to generate population (between-family) and within-sibship (within-family) GWAS estimates for 25 phenotypes. Within-sibship GWAS estimates were smaller than population estimates for height, educational attainment, age at first birth, number of children, cognitive ability, depressive symptoms and smoking. Some differences were observed in downstream SNP heritability, genetic correlations and Mendelian randomization analyses. For example, the within-sibship genetic correlation between educational attainment and body mass index attenuated towards zero. In contrast, analyses of most molecular phenotypes (for example, low-density lipoprotein-cholesterol) were generally consistent. We also found within-sibship evidence of polygenic adaptation on taller height. Here, we illustrate the importance of family-based GWAS data for phenotypes influenced by demographic and indirect genetic effects.

Sequencing of 640,000 exomes identifies GPR75 variants associated with protection from obesity

Large-scale human exome sequencing can identify rare protein-coding variants with a large impact on complex traits such as body adiposity. We sequenced the exomes of 645,626 individuals from the United Kingdom, the United States, and Mexico and estimated associations of rare coding variants with body mass index (BMI). We identified 16 genes with an exome-wide significant association with BMI, including those encoding five brain-expressed G protein-coupled receptors (CALCR, MC4R, GIPR, GPR151, and GPR75). Protein-truncating variants in GPR75 were observed in ~4/10,000 sequenced individuals and were associated with 1.8 kilograms per square meter lower BMI and 54% lower odds of obesity in the heterozygous state. Knock out of Gpr75 in mice resulted in resistance to weight gain and improved glycemic control in a high-fat diet model. Inhibition of GPR75 may provide a therapeutic strategy for obesity.

Common and rare variant associations with clonal haematopoiesis phenotypes

Clonal haematopoiesis involves the expansion of certain blood cell lineages and has been associated with ageing and adverse health outcomes1-5. Here we use exome sequence data on 628,388 individuals to identify 40,208 carriers of clonal haematopoiesis of indeterminate potential (CHIP). Using genome-wide and exome-wide association analyses, we identify 24 loci (21 of which are novel) where germline genetic variation influences predisposition to CHIP, including missense variants in the lymphocytic antigen coding gene LY75, which are associated with reduced incidence of CHIP. We also identify novel rare variant associations with clonal haematopoiesis and telomere length. Analysis of 5,041 health traits from the UK Biobank (UKB) found relationships between CHIP and severe COVID-19 outcomes, cardiovascular disease, haematologic traits, malignancy, smoking, obesity, infection and all-cause mortality. Longitudinal and Mendelian randomization analyses revealed that CHIP is associated with solid cancers, including non-melanoma skin cancer and lung cancer, and that CHIP linked to DNMT3A is associated with the subsequent development of myeloid but not lymphoid leukaemias. Additionally, contrary to previous findings from the initial 50,000 UKB exomes6, our results in the full sample do not support a role for IL-6 inhibition in reducing the risk of cardiovascular disease among CHIP carriers. Our findings demonstrate that CHIP represents a complex set of heterogeneous phenotypes with shared and unique germline genetic causes and varied clinical implications.

Patients With High Genome-Wide Polygenic Risk Scores for Coronary Artery Disease May Receive Greater Clinical Benefit From Alirocumab Treatment in the ODYSSEY OUTCOMES Trial

BACKGROUND

Alirocumab, an antibody that blocks PCSK9 (proprotein convertase subtilisin/kexin type 9), was associated with reduced major adverse cardiovascular events (MACE) and death in the ODYSSEY OUTCOMES trial (Evaluation of Cardiovascular Outcomes After an Acute Coronary Syndrome During Treatment With Alirocumab). In this study, higher baseline levels of low-density lipoprotein cholesterol (LDL-C) predicted greater benefit from alirocumab treatment. Recent studies indicate high polygenic risk scores (PRS) for coronary artery disease (CAD) identify individuals at higher risk who derive increased benefit from statins. We performed post hoc analyses to determine whether high PRS for CAD identifies higher-risk individuals, independent of baseline LDL-C and other known risk factors, who might derive greater benefit from alirocumab treatment.

METHODS

ODYSSEY OUTCOMES was a randomized, double-blind, placebo-controlled trial comparing alirocumab or placebo in 18 924 patients with acute coronary syndrome and elevated atherogenic lipoproteins despite optimized statin treatment. The primary endpoint (MACE) comprised death of CAD, nonfatal myocardial infarction, ischemic stroke, or unstable angina requiring hospitalization. A genome-wide PRS for CAD comprising 6 579 025 genetic variants was evaluated in 11 953 patients with available DNA samples. Analysis of MACE risk was performed in placebo-treated patients, whereas treatment benefit analysis was performed in all patients.

RESULTS

The incidence of MACE in the placebo group was related to PRS for CAD: 17.0% for high PRS patients (>90th percentile) and 11.4% for lower PRS patients (≤90th percentile; P<0.001); this PRS relationship was not explained by baseline LDL-C or other established risk factors. Both the absolute and relative reduction of MACE by alirocumab compared with placebo was greater in high versus low PRS patients. There was an absolute reduction by alirocumab in high versus low PRS groups of 6.0% and 1.5%, respectively, and a relative risk reduction by alirocumab of 37% in the high PRS group (hazard ratio, 0.63 [95% CI, 0.46-0.86]; P=0.004) versus a 13% reduction in the low PRS group (hazard ratio, 0.87 [95% CI, 0.78-0.98]; P=0.022; interaction P=0.04).

CONCLUSIONS

A high PRS for CAD is associated with elevated risk for recurrent MACE after acute coronary syndrome and a larger absolute and relative risk reduction with alirocumab treatment, providing an independent tool for risk stratification and precision medicine.

Exome Sequencing-Based Screening for BRCA1/2 Expected Pathogenic Variants Among Adult Biobank Participants

IMPORTANCE

Detection of disease-associated variants in the BRCA1 and BRCA2 (BRCA1/2) genes allows for cancer prevention and early diagnosis in high-risk individuals.

OBJECTIVES

To identify pathogenic and likely pathogenic (P/LP) BRCA1/2 variants in an unselected research cohort, and to characterize the features associated with P/LP variants.

DESIGN, SETTING, AND PARTICIPANTS

This is a cross-sectional study of adult volunteers (n = 50 726) who underwent exome sequencing at a single health care system (Geisinger Health System, Danville, Pennsylvania) from January 1, 2014, to March 1, 2016. Participants are part of the DiscovEHR cohort and were identified through the Geisinger MyCode Community Health Initiative. They consented to a research protocol that included sequencing and return of actionable test results. Clinical data from electronic health records and clinical visits were correlated with variants. Comparisons were made between those with (cases) and those without (controls) P/LP variants in BRCA1/2.

MAIN OUTCOMES

Prevalence of P/LP BRCA1/2 variants in cohort, proportion of variant carriers not previously ascertained through clinical testing, and personal and family history of relevant cancers among BRCA1/2 variant carriers and noncarriers.

RESULTS

Of the 50 726 health system patients who underwent exome sequencing, 50 459 (99.5%) had no expected pathogenic BRCA1/2 variants and 267 (0.5%) were BRCA1/2 carriers. Of the 267 cases (148 [55.4%] were women and 119 [44.6%] were men with a mean [range] age of 58.9 [23-90] years), 183 (68.5%) received clinically confirmed results in their electronic health record. Among the 267 participants with P/LP BRCA1/2 variants, 219 (82.0%) had no prior clinical testing, 95 (35.6%) had BRCA1 variants, and 172 (64.4%) had BRCA2 variants. Syndromic cancer diagnoses were present in 11 (47.8%) of the 23 deceased BRCA1/2 carriers and in 56 (20.9%) of all 267 BRCA1/2 carriers. Among women, 31 (20.9%) of 148 variant carriers had a personal history of breast cancer, compared with 1554 (5.2%) of 29 880 noncarriers (odds ratio [OR], 5.95; 95% CI, 3.88-9.13; P < .001). Ovarian cancer history was present in 15 (10.1%) of 148 variant carriers and in 195 (0.6%) of 29 880 variant noncarriers (OR, 18.30; 95% CI, 10.48-31.4; P < .001). Among 89 BRCA1/2 carriers without prior testing but with comprehensive personal and family history data, 44 (49.4%) did not meet published guidelines for clinical testing.

CONCLUSIONS AND RELEVANCE

This study found that compared with previous clinical care, exome sequencing-based screening identified 5 times as many individuals with P/LP BRCA1/2 variants. These findings suggest that genomic screening may identify BRCA1/2-associated cancer risk that might otherwise remain undetected within health care systems and may provide opportunities to reduce morbidity and mortality in patients.

Measuring Absolute RNA Copy Numbers at High Temporal Resolution Reveals Transcriptome Kinetics in Development

Transcript regulation is essential for cell function, and misregulation can lead to disease. Despite technologies to survey the transcriptome, we lack a comprehensive understanding of transcript kinetics, which limits quantitative biology. This is an acute challenge in embryonic development, where rapid changes in gene expression dictate cell fate decisions. By ultra-high-frequency sampling of Xenopus embryos and absolute normalization of sequence reads, we present smooth gene expression trajectories in absolute transcript numbers. During a developmental period approximating the first 8 weeks of human gestation, transcript kinetics vary by eight orders of magnitude. Ordering genes by expression dynamics, we find that "temporal synexpression" predicts common gene function. Remarkably, a single parameter, the characteristic timescale, can classify transcript kinetics globally and distinguish genes regulating development from those involved in cellular metabolism. Overall, our analysis provides unprecedented insight into the reorganization of maternal and embryonic transcripts and redefines our ability to perform quantitative biology.

De novo insertions and deletions of predominantly paternal origin are associated with autism spectrum disorder

Whole-exome sequencing (WES) studies have demonstrated the contribution of de novo loss-of-function single-nucleotide variants (SNVs) to autism spectrum disorder (ASD). However, challenges in the reliable detection of de novo insertions and deletions (indels) have limited inclusion of these variants in prior analyses. By applying a robust indel detection method to WES data from 787 ASD families (2,963 individuals), we demonstrate that de novo frameshift indels contribute to ASD risk (OR = 1.6; 95% CI = 1.0-2.7; p = 0.03), are more common in female probands (p = 0.02), are enriched among genes encoding FMRP targets (p = 6 × 10(-9)), and arise predominantly on the paternal chromosome (p < 0.001). On the basis of mutation rates in probands versus unaffected siblings, we conclude that de novo frameshift indels contribute to risk in approximately 3% of individuals with ASD. Finally, by observing clustering of mutations in unrelated probands, we uncover two ASD-associated genes: KMT2E (MLL5), a chromatin regulator, and RIMS1, a regulator of synaptic vesicle release.

Whole-exome sequencing characterizes the landscape of somatic mutations and copy number alterations in adrenocortical carcinoma

CONTEXT

Adrenocortical carcinoma (ACC) is a rare and lethal malignancy with a poorly defined etiology, and the molecular genetics of ACC are incompletely understood.

OBJECTIVE

To utilize whole-exome sequencing for genetic characterization of the underlying somatic mutations and copy number alterations present in ACC.

DESIGN

Screening for somatic mutation events and copy number alterations (CNAs) was performed by comparative analysis of tumors and matched normal samples from 41 patients with ACC.

RESULTS

In total, 966 nonsynonymous somatic mutations were detected, including 40 tumors with a mean of 16 mutations per sample and one tumor with 314 mutations. Somatic mutations in ACC-associated genes included TP53 (8/41 tumors, 19.5%) and CTNNB1 (4/41, 9.8%). Genes with potential disease-causing mutations included GNAS, NF2, and RB1, and recurrently mutated genes with unknown roles in tumorigenesis comprised CDC27, SCN7A, and SDK1. Recurrent CNAs included amplification at 5p15.33 including TERT (6/41, 14.6%) and homozygous deletion at 22q12.1 including the Wnt repressors ZNRF3 and KREMEN1 (4/41 9.8% and 3/41, 7.3%, respectively). Somatic mutations in ACC-established genes and recurrent ZNRF3 and TERT loci CNAs were mutually exclusive in the majority of cases. Moreover, gene ontology identified Wnt signaling as the most frequently mutated pathway in ACCs.

CONCLUSIONS

These findings highlight the importance of Wnt pathway dysregulation in ACC and corroborate the finding of homozygous deletion of Wnt repressors ZNRF3 and KREMEN1. Overall, mutations in either TP53 or CTNNB1 as well as focal CNAs at the ZNRF3 or TERT loci denote mutually exclusive events, suggesting separate mechanisms underlying the development of these tumors.

Genetic inactivation of ANGPTL4 improves glucose homeostasis and is associated with reduced risk of diabetes

Angiopoietin-like 4 (ANGPTL4) is an endogenous inhibitor of lipoprotein lipase that modulates lipid levels, coronary atherosclerosis risk, and nutrient partitioning. We hypothesize that loss of ANGPTL4 function might improve glucose homeostasis and decrease risk of type 2 diabetes (T2D). We investigate protein-altering variants in ANGPTL4 among 58,124 participants in the DiscovEHR human genetics study, with follow-up studies in 82,766 T2D cases and 498,761 controls. Carriers of p.E40K, a variant that abolishes ANGPTL4 ability to inhibit lipoprotein lipase, have lower odds of T2D (odds ratio 0.89, 95% confidence interval 0.85-0.92, p = 6.3 × 10-10), lower fasting glucose, and greater insulin sensitivity. Predicted loss-of-function variants are associated with lower odds of T2D among 32,015 cases and 84,006 controls (odds ratio 0.71, 95% confidence interval 0.49-0.99, p = 0.041). Functional studies in Angptl4-deficient mice confirm improved insulin sensitivity and glucose homeostasis. In conclusion, genetic inactivation of ANGPTL4 is associated with improved glucose homeostasis and reduced risk of T2D.

Association of the V122I Hereditary Transthyretin Amyloidosis Genetic Variant With Heart Failure Among Individuals of African or Hispanic/Latino Ancestry

Importance

Hereditary transthyretin (TTR) amyloid cardiomyopathy (hATTR-CM) due to the TTR V122I variant is an autosomal-dominant disorder that causes heart failure in elderly individuals of African ancestry. The clinical associations of carrying the variant, its effect in other African ancestry populations including Hispanic/Latino individuals, and the rates of achieving a clinical diagnosis in carriers are unknown.

Objective

To assess the association between the TTR V122I variant and heart failure and identify rates of hATTR-CM diagnosis among carriers with heart failure.

Design, Setting, and Participants

Cross-sectional analysis of carriers and noncarriers of TTR V122I of African ancestry aged 50 years or older enrolled in the Penn Medicine Biobank between 2008 and 2017 using electronic health record data from 1996 to 2017. Case-control study in participants of African and Hispanic/Latino ancestry with and without heart failure in the Mount Sinai BioMe Biobank enrolled between 2007 and 2015 using electronic health record data from 2007 to 2018.

Exposures

TTR V122I carrier status.

Main Outcomes and Measures

The primary outcome was prevalent heart failure. The rate of diagnosis with hATTR-CM among TTR V122I carriers with heart failure was measured.

Results

The cross-sectional cohort included 3724 individuals of African ancestry with a median age of 64 years (interquartile range, 57-71); 1755 (47%) were male, 2896 (78%) had a diagnosis of hypertension, and 753 (20%) had a history of myocardial infarction or coronary revascularization. There were 116 TTR V122I carriers (3.1%); 1121 participants (30%) had heart failure. The case-control study consisted of 2307 individuals of African ancestry and 3663 Hispanic/Latino individuals; the median age was 73 years (interquartile range, 68-80), 2271 (38%) were male, 4709 (79%) had a diagnosis of hypertension, and 1008 (17%) had a history of myocardial infarction or coronary revascularization. There were 1376 cases of heart failure. TTR V122I was associated with higher rates of heart failure (cross-sectional cohort: n = 51/116 TTR V122I carriers [44%], n = 1070/3608 noncarriers [30%], adjusted odds ratio, 1.7 [95% CI, 1.2-2.4], P = .006; case-control study: n = 36/1376 heart failure cases [2.6%], n = 82/4594 controls [1.8%], adjusted odds ratio, 1.8 [95% CI, 1.2-2.7], P = .008). Ten of 92 TTR V122I carriers with heart failure (11%) were diagnosed as having hATTR-CM; the median time from onset of symptoms to clinical diagnosis was 3 years.

Conclusions and Relevance

Among individuals of African or Hispanic/Latino ancestry enrolled in 2 academic medical center-based biobanks, the TTR V122I genetic variant was significantly associated with heart failure.

Rare variants in SOX17 are associated with pulmonary arterial hypertension with congenital heart disease

BACKGROUND

Pulmonary arterial hypertension (PAH) is a rare disease characterized by distinctive changes in pulmonary arterioles that lead to progressive pulmonary arterial pressures, right-sided heart failure, and a high mortality rate. Up to 30% of adult and 75% of pediatric PAH cases are associated with congenital heart disease (PAH-CHD), and the underlying etiology is largely unknown. There are no known major risk genes for PAH-CHD.

METHODS

To identify novel genetic causes of PAH-CHD, we performed whole exome sequencing in 256 PAH-CHD patients. We performed a case-control gene-based association test of rare deleterious variants using 7509 gnomAD whole genome sequencing population controls. We then screened a separate cohort of 413 idiopathic and familial PAH patients without CHD for rare deleterious variants in the top association gene.

RESULTS

We identified SOX17 as a novel candidate risk gene (p = 5.5e-7). SOX17 is highly constrained and encodes a transcription factor involved in Wnt/β-catenin and Notch signaling during development. We estimate that rare deleterious variants contribute to approximately 3.2% of PAH-CHD cases. The coding variants identified include likely gene-disrupting (LGD) and deleterious missense, with most of the missense variants occurring in a highly conserved HMG-box protein domain. We further observed an enrichment of rare deleterious variants in putative targets of SOX17, many of which are highly expressed in developing heart and pulmonary vasculature. In the cohort of PAH without CHD, rare deleterious variants of SOX17 were observed in 0.7% of cases.

CONCLUSIONS

These data strongly implicate SOX17 as a new risk gene contributing to PAH-CHD as well as idiopathic/familial PAH. Replication in other PAH cohorts and further characterization of the clinical phenotype will be important to confirm the precise role of SOX17 and better estimate the contribution of genes regulated by SOX17.