Alexander H. Li

TGFB2 mutations cause familial thoracic aortic aneurysms and dissections associated with mild systemic features of Marfan syndrome

A predisposition for thoracic aortic aneurysms leading to acute aortic dissections can be inherited in families in an autosomal dominant manner. Genome-wide linkage analysis of two large unrelated families with thoracic aortic disease followed by whole-exome sequencing of affected relatives identified causative mutations in TGFB2. These mutations—a frameshift mutation in exon 6 and a nonsense mutation in exon 4—segregated with disease with a combined logarithm of odds (LOD) score of 7.7. Sanger sequencing of 276 probands from families with inherited thoracic aortic disease identified 2 additional TGFB2 mutations. TGFB2 encodes transforming growth factor (TGF)-β2, and the mutations are predicted to cause haploinsufficiency for TGFB2; however, aortic tissue from cases paradoxically shows increased TGF-β2 expression and immunostaining. Thus, haploinsufficiency for TGFB2 predisposes to thoracic aortic disease, suggesting that the initial pathway driving disease is decreased cellular TGF-β2 levels leading to a secondary increase in TGF-β2 production in the diseased aorta.

Whole-genome sequence-based analysis of high-density lipoprotein cholesterol.

We describe initial steps for interrogating whole-genome sequence data to characterize the genetic architecture of a complex trait, levels of high-density lipoprotein cholesterol (HDL-C). We report whole-genome sequencing and analysis of 962 individuals from the Cohorts for Heart and Aging Research in Genetic Epidemiology (CHARGE) studies. From this analysis, we estimate that common variation contributes more to heritability of HDL-C levels than rare variation, and screening for mendelian variants for dyslipidemia identified individuals with extreme HDL-C levels. Whole-genome sequencing analyses highlight the value of regulatory and non-protein-coding regions of the genome in addition to protein-coding regions.

Analysis of loss-of-function variants and 20 risk factor phenotypes in 8,554 individuals identifies loci influencing chronic disease

A typical human exome harbors dozens of loss-of-function (LOF) variants, which can lower disease risk factor levels and affect drug efficacy. We hypothesized that LOF variants are enriched in genes influencing risk factor levels and the onset of common chronic diseases, such as cardiovascular disease and diabetes. To test this hypothesis, we sequenced the exomes of 8,554 individuals and analyzed the effects of predicted LOF variants on 20 chronic disease risk factor phenotypes. Analysis of this sample as discovery and replication strata of equal size verified two relationships in well-studied genes (PCSK9 and APOC3) and identified eight new loci. Previously unknown relationships included elevated fasting glucose in carriers of heterozygous LOF variation in TXNDC5, which encodes a biomarker for type 1 diabetes progression, and apparent recessive effects of C1QTNF8 on serum magnesium levels. These data demonstrate the utility of functional-variant annotation within a large sample of deeply phenotyped individuals for gene discovery.

Association of Rare Loss-Of-Function Alleles in HAL, Serum Histidine: Levels and Incident Coronary Heart Disease

Background— Histidine is a semiessential amino acid with antioxidant and anti-inflammatory properties. Few data are available on the associations between genetic variants, histidine levels, and incident coronary heart disease (CHD) in a population-based sample. Methods and Results— By conducting whole exome sequencing on 1152 African Americans in the Atherosclerosis Risk in Communities (ARIC) study and focusing on loss-of-function (LoF) variants, we identified 3 novel rare LoF variants in HAL , a gene that encodes histidine ammonia-lyase in the first step of histidine catabolism. These LoF variants had large effects on blood histidine levels ( β =0.26; P =1.2×10−13). The positive association with histidine levels was replicated by genotyping an independent sample of 718 ARIC African Americans (minor allele frequency=1%; P =1.2×10−4). In addition, high blood histidine levels were associated with reduced risk of developing incident CHD with an average of 21.5 years of follow-up among African Americans (hazard ratio=0.18; P =1.9×10−4). This finding was validated in an independent sample of European Americans from the Framingham Heart Study (FHS) Offspring Cohort. However, LoF variants in HAL were not directly significantly associated with incident CHD after meta-analyzing results from the CHARGE Consortium. Conclusions— Three LoF mutations in HAL were associated with increased histidine levels, which in turn were shown to be inversely related to the risk of CHD among both African Americans and European Americans. Future investigations on the association between HAL gene variation and CHD are warranted.Background—Histidine is a semiessential amino acid with antioxidant and anti-inflammatory properties. Few data are available on the associations between genetic variants, histidine levels, and incident coronary heart disease (CHD) in a population-based sample. Methods and Results—By conducting whole exome sequencing on 1152 African Americans in the Atherosclerosis Risk in Communities (ARIC) study and focusing on loss-of-function (LoF) variants, we identified 3 novel rare LoF variants in HAL, a gene that encodes histidine ammonia-lyase in the first step of histidine catabolism. These LoF variants had large effects on blood histidine levels (&bgr;=0.26; P=1.2×10−13). The positive association with histidine levels was replicated by genotyping an independent sample of 718 ARIC African Americans (minor allele frequency=1%; P=1.2×10−4). In addition, high blood histidine levels were associated with reduced risk of developing incident CHD with an average of 21.5 years of follow-up among African Americans (hazard ratio=0.18; P=1.9×10−4). This finding was validated in an independent sample of European Americans from the Framingham Heart Study (FHS) Offspring Cohort. However, LoF variants in HAL were not directly significantly associated with incident CHD after meta-analyzing results from the CHARGE Consortium. Conclusions—Three LoF mutations in HAL were associated with increased histidine levels, which in turn were shown to be inversely related to the risk of CHD among both African Americans and European Americans. Future investigations on the association between HAL gene variation and CHD are warranted.

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.)

WGSA: an annotation pipeline for human genome sequencing studies.

DNA sequencing technologies continue to make progress in increased throughput and quality, and decreased cost. As we transition from whole exome capture sequencing to whole genome sequencing (WGS), our ability to convert machine-generated variant calls, including single nucleotide variant (SNV) and insertion-deletion variants (indels), into human-interpretable knowledge has lagged far behind the ability to obtain enormous amounts of variants. To help narrow this gap, here we present WGSA (WGS annotator), a functional annotation pipeline for human genome sequencing studies, which is runnable out of the box on the Amazon Compute Cloud and is freely downloadable at (https://sites.google.com/site/jpopgen/wgsa/). Functional annotation is a key step in WGS analysis. In one way, annotation helps the analyst filter to a subset of elements of particular interest (eg, cell type specific enhancers), in another way annotation helps the investigators to increase the power of identifying phenotype-associated loci (eg, association test using functional prediction score as a weight) and interpret potentially interesting findings. Currently, there are several popular gene model based annotation tools, including ANNOVAR,1 SnpEff2 and the Ensembl Variant Effect Predictor (VEP).3 These can annotate a variety of protein coding and non-coding gene models from a range of species. It is well known among practitioners that different databases (eg, RefSeq4 and Ensembl5) use different models for …

FBN1 contributing to familial congenital diaphragmatic hernia

Congenital diaphragmatic hernia (CDH) is a relatively common, life‐threatening birth defect. We present a family with recurrent CDH—paraesophageal and central—for whom exome sequencing (ES) revealed a frameshift mutation (c.4969_4970insA, p.Ile1657Asnfs*30) in the fibrillin 1 gene (FBN1) that causes Marfan syndrome. A diagnosis of Marfan syndrome had not been considered previously in this family. However, a review of the literature demonstrated that FBN1 mutations have an unusual pattern of CDH in which paraesophageal hernias are particularly common. Subsequent clinical evaluations revealed evidence for ectopia lentis in affected family members supporting a clinical diagnosis of Marfan syndrome. Since only two other cases of familial CDH have been described in association with FBN1 mutations, we investigated an oligogenic hypothesis by examining ES data for deleterious sequence changes in other CDH‐related genes. This search revealed putatively deleterious sequence changes in four other genes that have been shown to cause diaphragm defects in humans and/or mice—FREM1, DES, PAX3 and MET. It is unclear whether these changes, alone or in aggregate, are contributing to the development of CDH in this family. However, their individual contribution is likely to be small compared to that of the frameshift mutation in FBN1. We conclude that ES can be used to identify both major and minor genetic factors that may contribute to CDH. These results also suggest that ES should be considered in the diagnostic evaluation of individuals and families with CDH, particularly when other diagnostic modalities have failed to reveal a molecular etiology.

Loss-of-function variants influence the human serum metabolome

Loss-of-function variants, which often lead to greatly truncated protein product, influence human metabolite levels. The metabolome is a collection of small molecules resulting from multiple cellular and biological processes that can act as biomarkers of disease, and African-Americans exhibit high levels of genetic diversity. Exome sequencing of a sample of deeply phenotyped African-Americans allowed us to analyze the effects of annotated loss-of-function (LoF) mutations on 308 serum metabolites measured by untargeted liquid and gas chromatography coupled with mass spectrometry. In an independent sample, we identified and replicated four genes harboring six LoF mutations that significantly affected five metabolites. These sites were related to a 19 to 45% difference in geometric mean metabolite levels, with an average effect size of 25%. We show that some of the affected metabolites are risk predictors or diagnostic biomarkers of disease and, using the principle of Mendelian randomization, are in the causal pathway of disease. For example, LoF mutations in SLCO1B1 elevate the levels of hexadecanedioate, a fatty acid significantly associated with increased blood pressure levels and risk of incident heart failure in both African-Americans and an independent sample of European-Americans. We show that SLCO1B1 LoF mutations significantly increase the risk of incident heart failure, thus implicating the metabolite in the causal pathway of disease. These results reveal new avenues into gene function and the understanding of disease etiology by integrating -omic technologies into a deeply phenotyped population study.

Whole exome sequencing in 342 congenital cardiac left sided lesion cases reveals extensive genetic heterogeneity and complex inheritance patterns

BackgroundLeft-sided lesions (LSLs) account for an important fraction of severe congenital cardiovascular malformations (CVMs). The genetic contributions to LSLs are complex, and the mutations that cause these malformations span several diverse biological signaling pathways: TGFB, NOTCH, SHH, and more. Here, we use whole exome sequence data generated in 342 LSL cases to identify likely damaging variants in putative candidate CVM genes.MethodsUsing a series of bioinformatics filters, we focused on genes harboring population-rare, putative loss-of-function (LOF), and predicted damaging variants in 1760 CVM candidate genes constructed a priori from the literature and model organism databases. Gene variants that were not observed in a comparably sequenced control dataset of 5492 samples without severe CVM were then subjected to targeted validation in cases and parents. Whole exome sequencing data from 4593 individuals referred for clinical sequencing were used to bolster evidence for the role of candidate genes in CVMs and LSLs.ResultsOur analyses revealed 28 candidate variants in 27 genes, including 17 genes not previously associated with a human CVM disorder, and revealed diverse patterns of inheritance among LOF carriers, including 9 confirmed de novo variants in both novel and newly described human CVM candidate genes (ACVR1, JARID2, NR2F2, PLRG1, SMURF1) as well as established syndromic CVM genes (KMT2D, NF1, TBX20, ZEB2). We also identified two genes (DNAH5, OFD1) with evidence of recessive and hemizygous inheritance patterns, respectively. Within our clinical cohort, we also observed heterozygous LOF variants in JARID2 and SMAD1 in individuals with cardiac phenotypes, and collectively, carriers of LOF variants in our candidate genes had a four times higher odds of having CVM (odds ratio = 4.0, 95% confidence interval 2.5–6.5).ConclusionsOur analytical strategy highlights the utility of bioinformatic resources, including human disease records and model organism phenotyping, in novel gene discovery for rare human disease. The results underscore the extensive genetic heterogeneity underlying non-syndromic LSLs, and posit potential novel candidate genes and complex modes of inheritance in this important group of birth defects.

Genetics of blood lipids among ~300,000 multi-ethnic participants of the Million Veteran Program

The Million Veteran Program (MVP) was established in 2011 as a national research initiative to determine how genetic variation influences the health of US military veterans. Here we genotyped 312,571 MVP participants using a custom biobank array and linked the genetic data to laboratory and clinical phenotypes extracted from electronic health records covering a median of 10.0 years of follow-up. Among 297,626 veterans with at least one blood lipid measurement, including 57,332 black and 24,743 Hispanic participants, we tested up to around 32 million variants for association with lipid levels and identified 118 novel genome-wide significant loci after meta-analysis with data from the Global Lipids Genetics Consortium (total n > 600,000). Through a focus on mutations predicted to result in a loss of gene function and a phenome-wide association study, we propose novel indications for pharmaceutical inhibitors targeting PCSK9 (abdominal aortic aneurysm), ANGPTL4 (type 2 diabetes) and PDE3B (triglycerides and coronary disease).Analysis of genetic data and blood lipid measurements from over 300,000 participants in the Million Veteran Program identifies new associations for blood lipid traits.