Genetic Control of mRNA Splicing as a Potential Mechanism for Incomplete Penetrance of Rare Coding Variants

Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-sequencing data in GTEx V8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased whole genome sequencing data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants

Whole-Genome Sequencing of Pharmacogenetic Drug Response in Racially Diverse Children with Asthma

RATIONALE: Albuterol, a bronchodilator medication, is the first-line therapy for asthma worldwide. There are significant racial/ethnic differences in albuterol drug response. OBJECTIVES: To identify genetic variants important for bronchodilator drug response (BDR) in racially diverse children. METHODS: We performed the first whole-genome sequencing pharmacogenetics study from 1,441 children with asthma from the tails of the BDR distribution to identify genetic association with BDR. MEASUREMENTS AND MAIN RESULTS: We identified population-specific and shared genetic variants associated with BDR, including genome-wide significant (P \u3c 3.53 × 10 CONCLUSIONS: The lack of minority data, despite a collaboration of eight universities and 13 individual laboratories, highlights the urgent need for a dedicated national effort to prioritize diversity in research. Our study expands the understanding of pharmacogenetic analyses in racially/ethnically diverse populations and advances the foundation for precision medicine in at-risk and understudied minority populations

Whole-genome sequencing of pharmacogenetic drug response in racially diverse children with asthma

RATIONALE: Albuterol, a bronchodilator medication, is the first-line therapy for asthma worldwide. There are significant racial/ethnic differences in albuterol drug response. OBJECTIVES: To identify genetic variants important for bronchodilator drug response (BDR) in racially diverse children. METHODS: We performed the first whole-genome sequencing pharmacogenetics study from 1,441 children with asthma from the tails of the BDR distribution to identify genetic association with BDR. MEASUREMENTS AND MAIN RESULTS: We identified population-specific and shared genetic variants associated with BDR, including genome-wide significant (P \u3c 3.53 × 10-7) and suggestive (P \u3c 7.06 × 10-6) loci near genes previously associated with lung capacity (DNAH5), immunity (NFKB1 and PLCB1), and β-adrenergic signaling (ADAMTS3 and COX18). Functional analyses of the BDR-associated SNP in NFKB1 revealed potential regulatory function in bronchial smooth muscle cells. The SNP is also an expression quantitative trait locus for a neighboring gene, SLC39A8. The lack of other asthma study populations with BDR and whole-genome sequencing data on minority children makes it impossible to perform replication of our rare variant associations. Minority underrepresentation also poses significant challenges to identify age-matched and population-matched cohorts of sufficient sample size for replication of our common variant findings. CONCLUSIONS: The lack of minority data, despite a collaboration of eight universities and 13 individual laboratories, highlights the urgent need for a dedicated national effort to prioritize diversity in research. Our study expands the understanding of pharmacogenetic analyses in racially/ethnically diverse populations and advances the foundation for precision medicine in at-risk and understudied minority populations

Common Polymorphisms in MTNR1B, G6PC2 and GCK Are Associated with Increased Fasting Plasma Glucose and Impaired Beta-Cell Function in Chinese Subjects

BACKGROUND: Previous studies identified melatonin receptor 1B (MTNR1B), islet-specific glucose 6 phosphatase catalytic subunit-related protein (G6PC2), glucokinase (GCK) and glucokinase regulatory protein (GCKR) as candidate genes for type 2 diabetes (T2D) acting through elevated fasting plasma glucose (FPG). We examined the associations of the reported common variants of these genes with T2D and glucose homeostasis in three independent Chinese cohorts. METHODOLOGY/PRINCIPAL FINDINGS: Five single nucleotide polymorphisms (SNPs), MTNR1B rs10830963, G6PC2 rs16856187 and rs478333, GCK rs1799884 and GCKR rs780094, were genotyped in 1644 controls (583 adults and 1061 adolescents) and 1342 T2D patients. The G-allele of MTNR1B rs10830963 and the C-alleles of both G6PC2 rs16856187 and rs478333 were associated with higher FPG (0.0034<P<6.6x10(-5)) in healthy controls. In addition to our previous report for association with FPG, the A-allele of GCK rs1799884 was also associated with reduced homeostasis model assessment of beta-cell function (HOMA-B) (P=0.0015). Together with GCKR rs780094, the risk alleles of these SNPs exhibited dosage effect in their associations with increased FPG (P=2.9x10(-9)) and reduced HOMA-B (P=1.1x10(-3)). Meta-analyses strongly supported additive effects of MTNR1B rs10830963 and G6PC2 rs16856187 on FPG. CONCLUSIONS/SIGNIFICANCE: Common variants of MTNR1B, G6PC2 and GCK are associated with elevated FPG and impaired insulin secretion, both individually and jointly, suggesting that these risk alleles may precipitate or perpetuate hyperglycemia in predisposed individuals

Rare Variants in Long Non-Coding RNAs Are Associated With Blood Lipid Levels in the TOPMed Whole-Genome Sequencing Study

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs

System-wide transcriptome damage and tissue identity loss in COVID-19 patients

The molecular mechanisms underlying the clinical manifestations of coronavirus disease 2019 (COVID-19), and what distinguishes them from common seasonal influenza virus and other lung injury states such as acute respiratory distress syndrome, remain poorly understood. To address these challenges, we combine transcriptional profiling of 646 clinical nasopharyngeal swabs and 39 patient autopsy tissues to define body-wide transcriptome changes in response to COVID-19. We then match these data with spatial protein and expression profiling across 357 tissue sections from 16 representative patient lung samples and identify tissue-compartment-specific damage wrought by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, evident as a function of varying viral loads during the clinical course of infection and tissue-type-specific expression states. Overall, our findings reveal a systemic disruption of canonical cellular and transcriptional pathways across all tissues, which can inform subsequent studies to combat the mortality of COVID-19 and to better understand the molecular dynamics of lethal SARS-CoV-2 and other respiratory infections., • Across all organs, fibroblast, and immune cell populations increase in COVID-19 patients • Organ-specific cell types and functional markers are lost in all COVID-19 tissue types • Lung compartment identity loss correlates with SARS-CoV-2 viral loads • COVID-19 uniquely disrupts co-occurrence cell type clusters (different from IAV/ARDS) , Park et al. report system-wide transcriptome damage and tissue identity loss wrought by SARS-CoV-2, influenza, and bacterial infection across multiple organs (heart, liver, lung, kidney, and lymph nodes) and provide a spatiotemporal landscape of COVID-19 in the lung