Impact of genome build on RNA-seq interpretation and diagnostics.

Transcriptomics is a powerful tool for unraveling the molecular effects of genetic variants and disease diagnosis. Prior studies have demonstrated that choice of genome build impacts variant interpretation and diagnostic yield for genomic analyses. To identify the extent genome build also impacts transcriptomics analyses, we studied the effect of the hg19, hg38, and CHM13 genome builds on expression quantification and outlier detection in 386 rare disease and familial control samples from both the Undiagnosed Diseases Network and Genomics Research to Elucidate the Genetics of Rare Disease Consortium. Across six routinely collected biospecimens, 61% of quantified genes were not influenced by genome build. However, we identified 1,492 genes with build-dependent quantification, 3,377 genes with build-exclusive expression, and 9,077 genes with annotation-specific expression across six routinely collected biospecimens, including 566 clinically relevant and 512 known OMIM genes. Further, we demonstrate that between builds for a given gene, a larger difference in quantification is well correlated with a larger change in expression outlier calling. Combined, we provide a database of genes impacted by build choice and recommend that transcriptomics-guided analyses and diagnoses are cross referenced with these data for robustness.

Loss of function of FAM177A1, a Golgi complex localized protein, causes a novel neurodevelopmental disorder.

PURPOSE: The function of FAM177A1 and its relationship to human disease is largely unknown. Recent studies have demonstrated FAM177A1 to be a critical immune-associated gene. One previous case study has linked FAM177A1 to a neurodevelopmental disorder in 4 siblings. METHODS: We identified 5 individuals from 3 unrelated families with biallelic variants in FAM177A1. The physiological function of FAM177A1 was studied in a zebrafish model organism and human cell lines with loss-of-function variants similar to the affected cohort. RESULTS: These individuals share a characteristic phenotype defined by macrocephaly, global developmental delay, intellectual disability, seizures, behavioral abnormalities, hypotonia, and gait disturbance. We show that FAM177A1 localizes to the Golgi complex in mammalian and zebrafish cells. Intersection of the RNA sequencing and metabolomic data sets from FAM177A1-deficient human fibroblasts and whole zebrafish larvae demonstrated dysregulation of pathways associated with apoptosis, inflammation, and negative regulation of cell proliferation. CONCLUSION: Our data shed light on the emerging function of FAM177A1 and defines FAM177A1-related neurodevelopmental disorder as a new clinical entity.

Transcriptional Responses to Priority Effects in Nectar Yeast.

Priority effects, where the order and timing of species arrival influence the assembly of ecological communities, have been observed in a variety of taxa and habitats. However, the genetic and molecular basis of priority effects remains unclear, hindering a better understanding of when priority effects will be strong. We sought to gain such an understanding for the nectar yeast Metschnikowia reukaufii commonly found in the nectar of our study plant, the hummingbird-pollinated Diplacus (Mimulus) aurantiacus. In this plant, M. reukaufii can experience strong priority effects when it reaches flowers after other nectar yeasts, such as M. rancensis. After inoculation into two contrasting types of synthetic nectar simulating early arrival of M. rancensis, we conducted whole-transcriptome sequencing of 108 strains of M. reukaufii. We found that several genes were differentially expressed in M. reukaufii strains when the nectar had been conditioned by growth of M. rancensis. Many of these genes were associated with amino acid metabolism, suggesting that M. reukaufii strains responded molecularly to the reduction in amino acid availability caused by M. rancensis. Furthermore, investigation of expression quantitative trait loci (eQTLs) revealed that genes involved in amino acid transport and resistance to antifungal compounds were enriched in some genetic variants of M. reukaufii. We also found that gene expression was associated with population growth rate, particularly when amino acids were limited. These results suggest that intraspecific genetic variation in the ability of nectar yeasts to respond to nutrient limitation and direct fungal competition underpins priority effects in this microbial system.

Integrative genomic analysis in African American children with asthma finds three novel loci associated with lung function.

Bronchodilator (BD) drugs are commonly prescribed for treatment and management of obstructive lung function present with diseases such as asthma. Administration of BD medication can partially or fully restore lung function as measured by pulmonary function tests. The genetics of baseline lung function measures taken before BD medication have been extensively studied, and the genetics of the BD response itself have received some attention. However, few studies have focused on the genetics of post-BD lung function. To address this gap, we analyzed lung function phenotypes in 1103 subjects from the Study of African Americans, Asthma, Genes, and Environment, a pediatric asthma case-control cohort, using an integrative genomic analysis approach that combined genotype, locus-specific genetic ancestry, and functional annotation information. We integrated genome-wide association study (GWAS) results with an admixture mapping scan of three pulmonary function tests (forced expiratory volume in 1 s [FEV1 ], forced vital capacity [FVC], and FEV1 /FVC) taken before and after albuterol BD administration on the same subjects, yielding six traits. We identified 18 GWAS loci, and five additional loci from admixture mapping, spanning several known and novel lung function candidate genes. Most loci identified via admixture mapping exhibited wide variation in minor allele frequency across genotyped global populations. Functional fine-mapping revealed an enrichment of epigenetic annotations from peripheral blood mononuclear cells, fetal lung tissue, and lung fibroblasts. Our results point to three novel potential genetic drivers of pre- and post-BD lung function: ADAMTS1, RAD54B, and EGLN3.

Integrating Gene Expression with Summary Association Statistics to Identify Genes Associated with 30 Complex Traits.

Although genome-wide association studies (GWASs) have identified thousands of risk loci for many complex traits and diseases, the causal variants and genes at these loci remain largely unknown. Here, we introduce a method for estimating the local genetic correlation between gene expression and a complex trait and utilize it to estimate the genetic correlation due to predicted expression between pairs of traits. We integrated gene expression measurements from 45 expression panels with summary GWAS data to perform 30 multi-tissue transcriptome-wide association studies (TWASs). We identified 1,196 genes whose expression is associated with these traits; of these, 168 reside more than 0.5 Mb away from any previously reported GWAS significant variant. We then used our approach to find 43 pairs of traits with significant genetic correlation at the level of predicted expression; of these, eight were not found through genetic correlation at the SNP level. Finally, we used bi-directional regression to find evidence that BMI causally influences triglyceride levels and that triglyceride levels causally influence low-density lipoprotein. Together, our results provide insight into the role of gene expression in the susceptibility of complex traits and diseases.

Ambient air pollution, asthma drug response, and telomere length in African American youth.

BACKGROUND: Telomere length (TL) can serve as a potential biomarker for conditions associated with chronic oxidative stress and inflammation, such as asthma. Air pollution can induce oxidative stress. Understanding the relationship between TL, asthma, and air pollution is important for identifying risk factors contributing to unhealthy aging in children. OBJECTIVES: We sought to investigate associations between exposures to ambient air pollutants and TL in African American children and adolescents and to examine whether African ancestry, asthma status, and steroid medication use alter the association. METHODS: Linear regression was used to examine associations between absolute telomere length (aTL) and estimated annual average residential ozone (O3) and fine particulate matter with a diameter of 2.5 µm or less (PM2.5) exposures in a cross-sectional analysis of 1072 children in an existing asthma case-control study. African ancestry, asthma status, and use of steroid medications were examined as effect modifiers. RESULTS: Participants' aTLs were measured by using quantitative PCR. A 1-ppb and 1 µg/m3 increase in annual average exposure to O3 and PM2.5 were associated with a decrease in aTL of 37.1 kilo-base pair (kb; 95% CI, -66.7 to -7.4 kb) and 57.1 kb (95% CI, -118.1 to 3.9 kb), respectively. African ancestry and asthma were not effect modifiers; however, exposure to steroid medications modified the relationships between TL and pollutants. Past-year exposure to O3 and PM2.5 was associated with shorter TLs in patients without steroid use. CONCLUSION: Exposure to air pollution was associated with shorter TLs in nonasthmatic children and adolescents. This was not the case for asthmatic children as a group, but those receiving steroid medication had less shortening than those not using steroids. Reduced exposure to air pollution in childhood might help to preserve TL.

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 < 3.53 × 10-7) and suggestive (P < 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.

Impact of genome build on RNA-seq interpretation and diagnostics.

Transcriptomics is a powerful tool for unraveling the molecular effects of genetic variants and disease diagnosis. Prior studies have demonstrated that choice of genome build impacts variant interpretation and diagnostic yield for genomic analyses. To identify the extent genome build also impacts transcriptomics analyses, we studied the effect of the hg19, hg38, and CHM13 genome builds on expression quantification and outlier detection in 386 rare disease and familial control samples from both the Undiagnosed Diseases Network (UDN) and Genomics Research to Elucidate the Genetics of Rare Disease (GREGoR) Consortium. We identified 2,800 genes with build-dependent quantification across six routinely-collected biospecimens, including 1,391 protein-coding genes and 341 known rare disease genes. We further observed multiple genes that only have detectable expression in a subset of genome builds. Finally, we characterized how genome build impacts the detection of outlier transcriptomic events. Combined, we provide a database of genes impacted by build choice, and recommend that transcriptomics-guided analyses and diagnoses are cross-referenced with these data for robustness.

Genetic architecture of cardiac dynamic flow volumes.

Cardiac blood flow is a critical determinant of human health. However, the definition of its genetic architecture is limited by the technical challenge of capturing dynamic flow volumes from cardiac imaging at scale. We present DeepFlow, a deep-learning system to extract cardiac flow and volumes from phase-contrast cardiac magnetic resonance imaging. A mixed-linear model applied to 37,653 individuals from the UK Biobank reveals genome-wide significant associations across cardiac dynamic flow volumes spanning from aortic forward velocity to aortic regurgitation fraction. Mendelian randomization reveals a causal role for aortic root size in aortic valve regurgitation. Among the most significant contributing variants, localizing genes (near ELN, PRDM6 and ADAMTS7) are implicated in connective tissue and blood pressure pathways. Here we show that DeepFlow cardiac flow phenotyping at scale, combined with genotyping data, reinforces the contribution of connective tissue genes, blood pressure and root size to aortic valve function.

Increasing equity in science requires better ethics training: A course by trainees, for trainees.

Despite the profound impacts of scientific research, few scientists have received the necessary training to productively discuss the ethical and societal implications of their work. To address this critical gap, we-a group of predominantly human genetics trainees-developed a course on genetics, ethics, and society. We intend for this course to serve as a template for other institutions and scientific disciplines. Our curriculum positions human genetics within its historical and societal context and encourages students to evaluate how societal norms and structures impact the conduct of scientific research. We demonstrate the utility of this course via surveys of enrolled students and provide resources and strategies for others hoping to teach a similar course. We conclude by arguing that if we are to work toward rectifying the inequities and injustices produced by our field, we must first learn to view our own research as impacting and being impacted by society.

Integration of transcriptomics and long-read genomics prioritizes structural variants in rare disease.

Rare structural variants (SVs) - insertions, deletions, and complex rearrangements - can cause Mendelian disease, yet they remain difficult to accurately detect and interpret. We sequenced and analyzed Oxford Nanopore long-read genomes of 68 individuals from the Undiagnosed Disease Network (UDN) with no previously identified diagnostic mutations from short-read sequencing. Using our optimized SV detection pipelines and 571 control long-read genomes, we detected 716 long-read rare (MAF < 0.01) SV alleles per genome on average, achieving a 2.4x increase from short-reads. To characterize the functional effects of rare SVs, we assessed their relationship with gene expression from blood or fibroblasts from the same individuals, and found that rare SVs overlapping enhancers were enriched (LOR = 0.46) near expression outliers. We also evaluated tandem repeat expansions (TREs) and found 14 rare TREs per genome; notably these TREs were also enriched near overexpression outliers. To prioritize candidate functional SVs, we developed Watershed-SV, a probabilistic model that integrates expression data with SV-specific genomic annotations, which significantly outperforms baseline models that don't incorporate expression data. Watershed-SV identified a median of eight high-confidence functional SVs per UDN genome. Notably, this included compound heterozygous deletions in FAM177A1 shared by two siblings, which were likely causal for a rare neurodevelopmental disorder. Our observations demonstrate the promise of integrating long-read sequencing with gene expression towards improving the prioritization of functional SVs and TREs in rare disease patients.

Transcriptomics and chromatin accessibility in multiple African population samples.

Mapping the functional human genome and impact of genetic variants is often limited to European-descendent population samples. To aid in overcoming this limitation, we measured gene expression using RNA sequencing in lymphoblastoid cell lines (LCLs) from 599 individuals from six African populations to identify novel transcripts including those not represented in the hg38 reference genome. We used whole genomes from the 1000 Genomes Project and 164 Maasai individuals to identify 8,881 expression and 6,949 splicing quantitative trait loci (eQTLs/sQTLs), and 2,611 structural variants associated with gene expression (SV-eQTLs). We further profiled chromatin accessibility using ATAC-Seq in a subset of 100 representative individuals, to identity chromatin accessibility quantitative trait loci (caQTLs) and allele-specific chromatin accessibility, and provide predictions for the functional effect of 78.9 million variants on chromatin accessibility. Using this map of eQTLs and caQTLs we fine-mapped GWAS signals for a range of complex diseases. Combined, this work expands global functional genomic data to identify novel transcripts, functional elements and variants, understand population genetic history of molecular quantitative trait loci, and further resolve the genetic basis of multiple human traits and disease.

Native American Ancestry and Air Pollution Interact to Impact Bronchodilator Response in Puerto Rican Children with Asthma.

Objective: Asthma is the most common chronic disease in children. Short-acting bronchodilator medications are the most commonly prescribed asthma treatment worldwide, regardless of disease severity. Puerto Rican children display the highest asthma morbidity and mortality of any US population. Alarmingly, Puerto Rican children with asthma display poor bronchodilator drug response (BDR). Reduced BDR may explain, in part, the increased asthma morbidity and mortality observed in Puerto Rican children with asthma. Gene-environment interactions may explain a portion of the heritability of BDR. We aimed to identify gene-environment interactions associated with BDR in Puerto Rican children with asthma. Setting: Genetic, environmental, and psycho-social data from the Genes-environments and Admixture in Latino Americans (GALA II) case-control study. Participants: Our discovery dataset consisted of 658 Puerto Rican children with asthma; our replication dataset consisted of 514 Mexican American children with asthma. Main Outcome Measures: We assessed the association of pairwise interaction models with BDR using ViSEN (Visualization of Statistical Epistasis Networks). Results: We identified a non-linear interaction between Native American genetic ancestry and air pollution significantly associated with BDR in Puerto Rican children with asthma. This interaction was robust to adjustment for age and sex but was not significantly associated with BDR in our replication population. Conclusions: Decreased Native American ancestry coupled with increased air pollution exposure was associated with increased BDR in Puerto Rican children with asthma. Our study acknowledges BDR's phenotypic complexity, and emphasizes the importance of integrating social, environmental, and biological data to further our understanding of complex disease.

Strategies for Pathway Analysis Using GWAS and WGS Data.

Single-allele study designs, commonly used in genome-wide association studies (GWAS) as well as the more recently developed whole genome sequencing (WGS) studies, are a standard approach for investigating the relationship of common variation within the human genome to a given phenotype of interest. However, single-allele association results published for many GWAS studies represent only the tip of the iceberg for the information that can be extracted from these datasets. The primary analysis strategy for GWAS entails association analysis in which only the single nucleotide polymorphisms (SNPs) with the strongest p-values are declared statistically significant due to issues arising from multiple testing and type I errors. Factors such as locus heterogeneity, epistasis, and multiple genes conferring small effects contribute to the complexity of the genetic models underlying phenotype expression. Thus, many biologically meaningful associations having lower effect sizes at individual genes are overlooked, making it difficult to separate true associations from a sea of false-positive associations. Organizing these individual SNPs into biologically meaningful groups to look at the overall effects of minor perturbations to genes and pathways is desirable. This pathway-based approach provides researchers with insight into the functional foundations of the phenotype being studied and allows testing of various genetic scenarios. © 2018 by John Wiley & Sons, Inc.

Dysregulated invertebrate tropomyosin-dectin-1 interaction confers susceptibility to allergic diseases.

The key factors underlying the development of allergic diseases-the propensity for a minority of individuals to develop dysfunctional responses to harmless environmental molecules-remain undefined. We report a pathway of immune counter-regulation that suppresses the development of aeroallergy and shrimp-induced anaphylaxis. In mice, signaling through epithelially expressed dectin-1 suppresses the development of type 2 immune responses through inhibition of interleukin-33 (IL-33) secretion and the subsequent recruitment of IL-13-producing innate lymphoid cells. Although this homeostatic pathway is functional in respiratory epithelial cells from healthy humans, it is dramatically impaired in epithelial cells from asthmatic and chronic rhinosinusitis patients, resulting in elevated IL-33 production. Moreover, we identify an association between a single-nucleotide polymorphism (SNP) in the dectin-1 gene loci and reduced pulmonary function in two cohorts of asthmatics. This intronic SNP is a predicted eQTL (expression quantitative trait locus) that is associated with reduced dectin-1 expression in human tissue. We identify invertebrate tropomyosin, a ubiquitous arthropod-derived molecule, as an immunobiologically relevant dectin-1 ligand that normally serves to restrain IL-33 release and dampen type 2 immunity in healthy individuals. However, invertebrate tropomyosin presented in the context of impaired dectin-1 function, as observed in allergic individuals, leads to unrestrained IL-33 secretion and skewing of immune responses toward type 2 immunity. Collectively, we uncover a previously unrecognized mechanism of protection against allergy to a conserved recognition element omnipresent in our environment.

Genetic Determinants of Telomere Length in African American Youth.

Telomere length (TL) is associated with numerous disease states and is affected by genetic and environmental factors. However, TL has been mostly studied in adult populations of European or Asian ancestry. These studies have identified 34 TL-associated genetic variants recently used as genetic proxies for TL. The generalizability of these associations to pediatric populations and racially diverse populations, specifically of African ancestry, remains unclear. Furthermore, six novel variants associated with TL in a population of European children have been identified but not validated. We measured TL from whole blood samples of 492 healthy African American youth (children and adolescents between 8 and 20 years old) and performed the first genome-wide association study of TL in this population. We were unable to replicate neither the 34 reported genetic associations found in adults nor the six genetic associations found in European children. However, we discovered a novel genome-wide significant association between TL and rs1483898 on chromosome 14. Our results underscore the importance of examining genetic associations with TL in diverse pediatric populations such as African Americans.