Subcutaneous adipose tissue splice quantitative trait loci reveal differences in isoform usage associated with cardiometabolic traits.

Alternate splicing events can create isoforms that alter gene function, and genetic variants associated with alternate gene isoforms may reveal molecular mechanisms of disease. We used subcutaneous adipose tissue of 426 Finnish men from the METSIM study and identified splice junction quantitative trait loci (sQTLs) for 6,077 splice junctions (FDR < 1%). In the same individuals, we detected expression QTLs (eQTLs) for 59,443 exons and 15,397 genes (FDR < 1%). We identified 595 genes with an sQTL and exon eQTL but no gene eQTL, which could indicate potential isoform differences. Of the significant sQTL signals, 2,114 (39.8%) included at least one proxy variant (linkage disequilibrium r2 > 0.8) located within an intron spanned by the splice junction. We identified 203 sQTLs that colocalized with 141 genome-wide association study (GWAS) signals for cardiometabolic traits, including 25 signals for lipid traits, 24 signals for body mass index (BMI), and 12 signals for waist-hip ratio adjusted for BMI. Among all 141 GWAS signals colocalized with an sQTL, we detected 26 that also colocalized with an exon eQTL for an exon skipped by the sQTL splice junction. At a GWAS signal for high-density lipoprotein cholesterol colocalized with an NR1H3 sQTL splice junction, we show that the alternative splice product encodes an NR1H3 transcription factor that lacks a DNA binding domain and fails to activate transcription. Together, these results detect splicing events and candidate mechanisms that may contribute to gene function at GWAS loci.

Evaluating the clinical utility of early exome sequencing in diverse pediatric outpatient populations in the North Carolina Clinical Genomic Evaluation of Next-generation Exome Sequencing (NCGENES) 2 study: a randomized controlled trial.

BACKGROUND: Exome sequencing (ES) has probable utility for shortening the diagnostic odyssey of children with suspected genetic disorders. This report describes the design and methods of a study evaluating the potential of ES as a routine clinical tool for pediatric patients who have suspected genetic conditions and who are in the early stages of the diagnostic odyssey. METHODS: The North Carolina Clinical Genomic Evaluation by Next-generation Exome Sequencing (NCGENES) 2 study is an interdisciplinary, multi-site Phase III randomized controlled trial of two interventions: educational pre-visit preparation (PVP) and offer of first-line ES. In this full-factorial design, parent-child dyads are randomly assigned to one of four study arms (PVP + usual care, ES + usual care, PVP + ES + usual care, or usual care alone) in equal proportions. Participants are recruited from Pediatric Genetics or Neurology outpatient clinics in three North Carolina healthcare facilities. Eligible pediatric participants are < 16 years old and have a first visit to a participating clinic, a suspected genetic condition, and an eligible parent/guardian to attend the clinic visit and complete study measures. The study oversamples participants from underserved and under-represented populations. Participants assigned to the PVP arms receive an educational booklet and question prompt list before clinical interactions. Randomization to offer of first-line ES is revealed after a child's clinic visit. Parents complete measures at baseline, pre-clinic, post-clinic, and two follow-up timepoints. Study clinicians provide phenotypic data and complete measures after the clinic visit and after returning results. Reportable study-related research ES results are confirmed in a CLIA-certified clinical laboratory. Results are disclosed to the parent by the clinical team. A community consultation team contributed to the development of study materials and study implementation methods and remains engaged in the project. DISCUSSION: NCGENES 2 will contribute valuable knowledge concerning technical, clinical, psychosocial, and health economic issues associated with using early diagnostic ES to shorten the diagnostic odyssey of pediatric patients with likely genetic conditions. Results will inform efforts to engage diverse populations in genomic medicine research and generate evidence that can inform policy, practice, and future research related to the utility of first-line diagnostic ES in health care. TRIAL REGISTRATION: ClinicalTrials.gov NCT03548779 . Registered on June 07, 2018.

Variant Classification Concordance using the ACMG-AMP Variant Interpretation Guidelines across Nine Genomic Implementation Research Studies.

Harmonization of variant pathogenicity classification across laboratories is important for advancing clinical genomics. The two CLIA-accredited Electronic Medical Record and Genomics Network sequencing centers and the six CLIA-accredited laboratories and one research laboratory performing genome or exome sequencing in the Clinical Sequencing Evidence-Generating Research Consortium collaborated to explore current sources of discordance in classification. Eight laboratories each submitted 20 classified variants in the ACMG secondary finding v.2.0 genes. After removing duplicates, each of the 158 variants was annotated and independently classified by two additional laboratories using the ACMG-AMP guidelines. Overall concordance across three laboratories was assessed and discordant variants were reviewed via teleconference and email. The submitted variant set included 28 P/LP variants, 96 VUS, and 34 LB/B variants, mostly in cancer (40%) and cardiac (27%) risk genes. Eighty-six (54%) variants reached complete five-category (i.e., P, LP, VUS, LB, B) concordance, and 17 (11%) had a discordance that could affect clinical recommendations (P/LP versus VUS/LB/B). 21% and 63% of variants submitted as P and LP, respectively, were discordant with VUS. Of the 54 originally discordant variants that underwent further review, 32 reached agreement, for a post-review concordance rate of 84% (118/140 variants). This project provides an updated estimate of variant concordance, identifies considerations for LP classified variants, and highlights ongoing sources of discordance. Continued and increased sharing of variant classifications and evidence across laboratories, and the ongoing work of ClinGen to provide general as well as gene- and disease-specific guidance, will lead to continued increases in concordance.

Genomic Sequencing for Newborn Screening: Results of the NC NEXUS Project.

Newborn screening (NBS) was established as a public health program in the 1960s and is crucial for facilitating detection of certain medical conditions in which early intervention can prevent serious, life-threatening health problems. Genomic sequencing can potentially expand the screening for rare hereditary disorders, but many questions surround its possible use for this purpose. We examined the use of exome sequencing (ES) for NBS in the North Carolina Newborn Exome Sequencing for Universal Screening (NC NEXUS) project, comparing the yield from ES used in a screening versus a diagnostic context. We enrolled healthy newborns and children with metabolic diseases or hearing loss (106 participants total). ES confirmed the participant's underlying diagnosis in 15 out of 17 (88%) children with metabolic disorders and in 5 out of 28 (∼18%) children with hearing loss. We discovered actionable findings in four participants that would not have been detected by standard NBS. A subset of parents was eligible to receive additional information for their child about childhood-onset conditions with low or no clinical actionability, clinically actionable adult-onset conditions, and carrier status for autosomal-recessive conditions. We found pathogenic variants associated with hereditary breast and/or ovarian cancer in two children, a likely pathogenic variant in the gene associated with Lowe syndrome in one child, and an average of 1.8 reportable variants per child for carrier results. These results highlight the benefits and limitations of using genomic sequencing for NBS and the challenges of using such technology in future precision medicine approaches.

Enhancer deletion and allelic effects define a regulatory molecular mechanism at the VLDLR cholesterol GWAS locus.

Total cholesterol (TC) and low-density lipoprotein cholesterol (LDL-C) are heritable risk factors for cardiovascular disease, yet the molecular mechanisms underlying the majority of blood lipid-associated genome-wide association studies signals remain elusive. One association signal is located in intron 3 of VLDLR; rs3780181-A is a risk allele associated (P ≤ 2 × 10-9) with increased TC and LDL-C. We investigated variants, genes and mechanisms underlying this association signal. We used a functional genetic approach to show that the intronic region spanning rs3780181 exhibited 1.6-7.6-fold enhancer activity in human HepG2 hepatocyte, THP-1 monocyte and Simpson-Golabi-Behmel Syndrome (SGBS) preadipocyte cells and that the rs3780181-A risk allele showed significantly less enhancer activity compared with the G allele, consistent with the direction of an expression quantitative trait locus in liver. In addition, rs3780181 alleles showed differential binding to multiple nuclear proteins, including stronger IRF2 binding to the rs3780181 G allele. We used a CRISPR-cas9 approach to delete 475 and 663 bp of the putative enhancer element in HEK293T kidney cells; compared to expression of mock-edited cell lines, the homozygous enhancer deletion cell lines showed 1.2-fold significantly (P < 0.04) decreased expression of VLDLR, as well as 1.5-fold decreased expression of SMARCA2, located 388 kb away. Together, these results identify an enhancer of VLDLR expression and suggest that altered binding of one or more factors bound to rs3780181 alleles decreases enhancer activity and reduces at least VLDLR expression, leading to increased TC and LDL-C.

Navigating the nuances of clinical sequence variant interpretation in Mendelian disease.

Understanding clinical genetic test results in the era of next-generation sequencing has become increasingly complex, necessitating clear and thorough guidelines for sequence variant interpretation. To meet this need the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) published guidelines for a systematic approach for sequence variant interpretation in 2015. This framework is intended to be adaptable to any Mendelian condition, promoting transparency and consistency in variant interpretation, yet its comprehensive nature yields important challenges and caveats that end users must understand. In this review, we address some of these nuances and discuss the evolving efforts to refine and adapt this framework. We also consider the added complexity of distinguishing between variant-level interpretations and case-level conclusions, particularly in the context of the large gene panel approach to clinical diagnostics.

Functional genomics and assays of regulatory activity detect mechanisms at loci for lipid traits and coronary artery disease.

Many genome-wide association studies (GWAS) have identified signals located in non-coding regions, and an increasing number of functional genomics annotations of regulatory elements and assays of regulatory activity have been used to investigate mechanisms. Genome-wide datasets that characterize chromatin structure help detect potential regulatory elements. Assays to experimentally assess candidate variants include transcriptional reporter assays, and recently, massively parallel reporter assays (MPRAs). Additionally, the effect of candidate regulatory elements and variants on gene expression and function can be evaluated using genomic editing with the CRISPR-Cas9 technology. We highlight some recent studies that employed these strategies to identify variant effects and elucidate molecular and/or biological mechanisms at GWAS loci for lipid traits and coronary artery disease.

A Type 2 Diabetes-Associated Functional Regulatory Variant in a Pancreatic Islet Enhancer at the ADCY5 Locus.

Molecular mechanisms remain unknown for most type 2 diabetes genome-wide association study identified loci. Variants associated with type 2 diabetes and fasting glucose levels reside in introns of ADCY5, a gene that encodes adenylate cyclase 5. Adenylate cyclase 5 catalyzes the production of cyclic AMP, which is a second messenger molecule involved in cell signaling and pancreatic ß-cell insulin secretion. We demonstrated that type 2 diabetes risk alleles are associated with decreased ADCY5 expression in human islets and examined candidate variants for regulatory function. rs11708067 overlaps a predicted enhancer region in pancreatic islets. The type 2 diabetes risk rs11708067-A allele showed fewer H3K27ac ChIP-seq reads in human islets, lower transcriptional activity in reporter assays in rodent ß-cells (rat 832/13 and mouse MIN6), and increased nuclear protein binding compared with the rs11708067-G allele. Homozygous deletion of the orthologous enhancer region in 832/13 cells resulted in a 64% reduction in expression level of Adcy5, but not adjacent gene Sec22a, and a 39% reduction in insulin secretion. Together, these data suggest that rs11708067-A risk allele contributes to type 2 diabetes by disrupting an islet enhancer, which results in reduced ADCY5 expression and impaired insulin secretion.

Association analyses of East Asian individuals and trans-ancestry analyses with European individuals reveal new loci associated with cholesterol and triglyceride levels.

Large-scale meta-analyses of genome-wide association studies (GWAS) have identified >175 loci associated with fasting cholesterol levels, including total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG). With differences in linkage disequilibrium (LD) structure and allele frequencies between ancestry groups, studies in additional large samples may detect new associations. We conducted staged GWAS meta-analyses in up to 69,414 East Asian individuals from 24 studies with participants from Japan, the Philippines, Korea, China, Singapore, and Taiwan. These meta-analyses identified (P < 5 × 10-8) three novel loci associated with HDL-C near CD163-APOBEC1 (P = 7.4 × 10-9), NCOA2 (P = 1.6 × 10-8), and NID2-PTGDR (P = 4.2 × 10-8), and one novel locus associated with TG near WDR11-FGFR2 (P = 2.7 × 10-10). Conditional analyses identified a second signal near CD163-APOBEC1. We then combined results from the East Asian meta-analysis with association results from up to 187,365 European individuals from the Global Lipids Genetics Consortium in a trans-ancestry meta-analysis. This analysis identified (log10Bayes Factor ≥6.1) eight additional novel lipid loci. Among the twelve total loci identified, the index variants at eight loci have demonstrated at least nominal significance with other metabolic traits in prior studies, and two loci exhibited coincident eQTLs (P < 1 × 10-5) in subcutaneous adipose tissue for BPTF and PDGFC. Taken together, these analyses identified multiple novel lipid loci, providing new potential therapeutic targets.

Multiple Hepatic Regulatory Variants at the GALNT2 GWAS Locus Associated with High-Density Lipoprotein Cholesterol.

Genome-wide association studies (GWASs) have identified more than 150 loci associated with blood lipid and cholesterol levels; however, the functional and molecular mechanisms for many associations are unknown. We examined the functional regulatory effects of candidate variants at the GALNT2 locus associated with high-density lipoprotein cholesterol (HDL-C). Fine-mapping and conditional analyses in the METSIM study identified a single locus harboring 25 noncoding variants (r(2) > 0.7 with the lead GWAS variants) strongly associated with total cholesterol in medium-sized HDL (e.g., rs17315646, p = 3.5 × 10(-12)). We used luciferase reporter assays in HepG2 cells to test all 25 variants for allelic differences in regulatory enhancer activity. rs2281721 showed allelic differences in transcriptional activity (75-fold [T] versus 27-fold [C] more than the empty-vector control), as did a separate 780-bp segment containing rs4846913, rs2144300, and rs6143660 (49-fold [AT(-) haplotype] versus 16-fold [CC(+) haplotype] more). Using electrophoretic mobility shift assays, we observed differential CEBPB binding to rs4846913, and we confirmed this binding in a native chromatin context by performing chromatin-immunoprecipitation (ChIP) assays in HepG2 and Huh-7 cell lines of differing genotypes. Additionally, sequence reads in HepG2 DNase-I-hypersensitivity and CEBPB ChIP-seq signals spanning rs4846913 showed significant allelic imbalance. Allelic-expression-imbalance assays performed with RNA from primary human hepatocyte samples and expression-quantitative-trait-locus (eQTL) data in human subcutaneous adipose tissue samples confirmed that alleles associated with increased HDL-C are associated with a modest increase in GALNT2 expression. Together, these data suggest that at least rs4846913 and rs2281721 play key roles in influencing GALNT2 expression at this HDL-C locus.

New genetic loci link adipose and insulin biology to body fat distribution.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.