A common polymorphism in the Intelectin-1 gene influences mucus plugging in severe asthma.

By incompletely understood mechanisms, type 2 (T2) inflammation present in the airways of severe asthmatics drives the formation of pathologic mucus which leads to airway mucus plugging. Here we investigate the molecular role and clinical significance of intelectin-1 (ITLN-1) in the development of pathologic airway mucus in asthma. Through analyses of human airway epithelial cells we find that ITLN1 gene expression is highly induced by interleukin-13 (IL-13) in a subset of metaplastic MUC5AC+ mucus secretory cells, and that ITLN-1 protein is a secreted component of IL-13-induced mucus. Additionally, we find ITLN-1 protein binds the C-terminus of the MUC5AC mucin and that its deletion in airway epithelial cells partially reverses IL-13-induced mucostasis. Through analysis of nasal airway epithelial brushings, we find that ITLN1 is highly expressed in T2-high asthmatics, when compared to T2-low children. Furthermore, we demonstrate that both ITLN-1 gene expression and protein levels are significantly reduced by a common genetic variant that is associated with protection from the formation of mucus plugs in T2-high asthma. This work identifies an important biomarker and targetable pathways for the treatment of mucus obstruction in asthma.

IL-4/STAT6 axis observed to reverse proliferative defect in SCA3 patient-derived neural progenitor cells.

Spinocerebellar ataxia 3 (SCA3) is an incurable, neurodegenerative genetic disorder that leads to progressive cerebellar ataxia and other parkinsonian-like pathologies because of loss of cerebellar neurons. The role of an expanded polyglutamine aggregate on neural progenitor cells is unknown. Here, we show that SCA3 patient-specific induced neural progenitor cells (iNPCs) exhibit proliferative defects. Moreover, SCA3 iNPCs have reduced autophagic expression compared to control. Furthermore, although SCA3 iNPCs continue to proliferate, they do not survive subsequent passages compared to control iNPCs, indicating the likelihood that SCA3 iNPCs undergo rapid senescence. Exposure to interleukin-4 (IL-4), a type 2 cytokine produced by immune cells, resulted in an observed increase in expression of autophagic programs and a reduction in the proliferation defect observed in SCA3 iNPCs. Our results indicate a previously unobserved role of SCA3 disease ontology on the neural stem cell pool and a potential therapeutic strategy using IL-4 to ameliorate or delay disease pathology in the SCA3 neural progenitor cell population.

Gene expression in African Americans, Puerto Ricans and Mexican Americans reveals ancestry-specific patterns of genetic architecture.

We explored ancestry-related differences in the genetic architecture of whole-blood gene expression using whole-genome and RNA sequencing data from 2,733 African Americans, Puerto Ricans and Mexican Americans. We found that heritability of gene expression significantly increased with greater proportions of African genetic ancestry and decreased with higher proportions of Indigenous American ancestry, reflecting the relationship between heterozygosity and genetic variance. Among heritable protein-coding genes, the prevalence of ancestry-specific expression quantitative trait loci (anc-eQTLs) was 30% in African ancestry and 8% for Indigenous American ancestry segments. Most anc-eQTLs (89%) were driven by population differences in allele frequency. Transcriptome-wide association analyses of multi-ancestry summary statistics for 28 traits identified 79% more gene-trait associations using transcriptome prediction models trained in our admixed population than models trained using data from the Genotype-Tissue Expression project. Our study highlights the importance of measuring gene expression across large and ancestrally diverse populations for enabling new discoveries and reducing disparities.

Multi-omic approach associates blood methylome with bronchodilator drug response in pediatric asthma.

BACKGROUND: Albuterol is the drug most widely used as asthma treatment among African Americans despite having a lower bronchodilator drug response (BDR) than other populations. Although BDR is affected by gene and environmental factors, the influence of DNA methylation is unknown. OBJECTIVE: This study aimed to identify epigenetic markers in whole blood associated with BDR, study their functional consequences by multi-omic integration, and assess their clinical applicability in admixed populations with a high asthma burden. METHODS: We studied 414 children and young adults (8-21 years old) with asthma in a discovery and replication design. We performed an epigenome-wide association study on 221 African Americans and replicated the results on 193 Latinos. Functional consequences were assessed by integrating epigenomics with genomics, transcriptomics, and environmental exposure data. Machine learning was used to develop a panel of epigenetic markers to classify treatment response. RESULTS: We identified 5 differentially methylated regions and 2 CpGs genome-wide significantly associated with BDR in African Americans located in FGL2 (cg08241295, P = 6.8 × 10-9) and DNASE2 (cg15341340, P = 7.8 × 10-8), which were regulated by genetic variation and/or associated with gene expression of nearby genes (false discovery rate < 0.05). The CpG cg15341340 was replicated in Latinos (P = 3.5 × 10-3). Moreover, a panel of 70 CpGs showed good classification for those with response and nonresponse to albuterol therapy in African American and Latino children (area under the receiver operating characteristic curve for training, 0.99; for validation, 0.70-0.71). The DNA methylation model showed similar discrimination as clinical predictors (P > .05). CONCLUSIONS: We report novel associations of epigenetic markers with BDR in pediatric asthma and demonstrate for the first time the applicability of pharmacoepigenetics in precision medicine of respiratory diseases.

Admixture mapping of severe asthma exacerbations in Hispanic/Latino children and youth.

BACKGROUND: In the USA, genetically admixed populations have the highest asthma prevalence and severe asthma exacerbations rates. This could be explained not only by environmental factors but also by genetic variants that exert ethnic-specific effects. However, no admixture mapping has been performed for severe asthma exacerbations. OBJECTIVE: We sought to identify genetic variants associated with severe asthma exacerbations in Hispanic/Latino subgroups by means of admixture mapping analyses and fine mapping, and to assess their transferability to other populations and potential functional roles. METHODS: We performed an admixture mapping in 1124 Puerto Rican and 625 Mexican American children with asthma. Fine-mapping of the significant peaks was performed via allelic testing of common and rare variants. We performed replication across Hispanic/Latino subgroups, and the transferability to non-Hispanic/Latino populations was assessed in 1001 African Americans, 1250 Singaporeans and 941 Europeans with asthma. The effects of the variants on gene expression and DNA methylation from whole blood were also evaluated in participants with asthma and in silico with data obtained through public databases. RESULTS: Genomewide significant associations of Indigenous American ancestry with severe asthma exacerbations were found at 5q32 in Mexican Americans as well as at 13q13-q13.2 and 3p13 in Puerto Ricans. The single nucleotide polymorphism (SNP) rs1144986 (C5orf46) showed consistent effects for severe asthma exacerbations across Hispanic/Latino subgroups, but it was not validated in non-Hispanics/Latinos. This SNP was associated with DPYSL3 DNA methylation and SCGB3A2 gene expression levels. CONCLUSIONS: Admixture mapping study of asthma exacerbations revealed a novel locus that exhibited Hispanic/Latino-specific effects and regulated DPYSL3 and SCGB3A2.

SOX9 and SOX10 control fluid homeostasis in the inner ear for hearing through independent and cooperative mechanisms.

The in vivo mechanisms underlying dominant syndromes caused by mutations in SRY-Box Transcription Factor 9 (SOX9) and SOX10 (SOXE) transcription factors, when they either are expressed alone or are coexpressed, are ill-defined. We created a mouse model for the campomelic dysplasia SOX9Y440X mutation, which truncates the transactivation domain but leaves DNA binding and dimerization intact. Here, we find that SOX9Y440X causes deafness via distinct mechanisms in the endolymphatic sac (ES)/duct and cochlea. By contrast, conditional heterozygous Sox9-null mice are normal. During the ES development of Sox9Y440X/+ heterozygotes, Sox10 and genes important for ionic homeostasis are down-regulated, and there is developmental persistence of progenitors, resulting in fewer mature cells. Sox10 heterozygous null mutants also display persistence of ES/duct progenitors. By contrast, SOX10 retains its expression in the early Sox9Y440X/+ mutant cochlea. Later, in the postnatal stria vascularis, dominant interference by SOX9Y440X is implicated in impairing the normal cooperation of SOX9 and SOX10 in repressing the expression of the water channel Aquaporin 3, thereby contributing to endolymphatic hydrops. Our study shows that for a functioning endolymphatic system in the inner ear, SOX9 regulates Sox10, and depending on the cell type and target gene, it works either independently of or cooperatively with SOX10. SOX9Y440X can interfere with the activity of both SOXE factors, exerting effects that can be classified as haploinsufficient/hypomorphic or dominant negative depending on the cell/gene context. This model of disruption of transcription factor partnerships may be applicable to congenital deafness, which affects ∼0.3% of newborns, and other syndromic disorders.

Genome-wide association study in minority children with asthma implicates DNAH5 in bronchodilator responsiveness.

Variability in response to short-acting ß2-agonists (e.g., albuterol) among patients with asthma from diverse racial/ethnic groups may contribute to asthma disparities. We sought to identify genetic variants associated with bronchodilator response (BDR) to identify potential mechanisms of drug response and risk factors for worse asthma outcomes. Genome-wide association studies of bronchodilator response (BDR) were performed using TOPMed Whole Genome Sequencing data of the Asthma Translational Genomic Collaboration (ATGC), which corresponded to 1136 Puerto Rican, 656 Mexican and 4337 African American patients with asthma. With the population-specific GWAS results, a trans-ethnic meta-analysis was performed to identify BDR-associated variants shared across the three populations. Replication analysis was carried out in three pediatric asthma cohorts, including CAMP (Childhood Asthma Management Program; n = 560), GACRS (Genetics of Asthma in Costa Rica Study; n = 967) and HPR (Hartford-Puerto Rico; n = 417). A genome-wide significant locus (rs35661809; P = 3.61 × 10-8) in LINC02220, a non-coding RNA gene, was identified in Puerto Ricans. While this region was devoid of protein-coding genes, capture Hi-C data showed a distal interaction with the promoter of the DNAH5 gene in lung tissue. In replication analysis, the GACRS cohort yielded a nominal association (1-tailed P < 0.05). No genetic variant was associated with BDR at the genome-wide significant threshold in Mexicans and African Americans. Our findings help inform genetic underpinnings of BDR for understudied minority patients with asthma, but the limited availability of genetic data for racial/ethnic minority children with asthma remains a paramount challenge.

Nasal airway transcriptome-wide association study of asthma reveals genetically driven mucus pathobiology.

To identify genetic determinants of airway dysfunction, we performed a transcriptome-wide association study for asthma by combining RNA-seq data from the nasal airway epithelium of 681 children, with UK Biobank genetic association data. Our airway analysis identified 95 asthma genes, 58 of which were not identified by transcriptome-wide association analyses using other asthma-relevant tissues. Among these genes were MUC5AC, an airway mucin, and FOXA3, a transcriptional driver of mucus metaplasia. Muco-ciliary epithelial cultures from genotyped donors revealed that the MUC5AC risk variant increases MUC5AC protein secretion and mucus secretory cell frequency. Airway transcriptome-wide association analyses for mucus production and chronic cough also identified MUC5AC. These cis-expression variants were associated with trans effects on expression; the MUC5AC variant was associated with upregulation of non-inflammatory mucus secretory network genes, while the FOXA3 variant was associated with upregulation of type-2 inflammation-induced mucus-metaplasia pathway genes. Our results reveal genetic mechanisms of airway mucus pathobiology.

Epigenome-wide association study of lung function in Latino children and youth with asthma.

INTRODUCTION: DNA methylation studies have associated methylation levels at different CpG sites or genomic regions with lung function. Moreover, genetic ancestry has been associated with lung function in Latinos. However, no epigenome-wide association study (EWAS) of lung function has been performed in this population. Here, we aimed to identify DNA methylation patterns associated with lung function in pediatric asthma among Latinos. RESULTS: We conducted an EWAS in whole blood from 250 Puerto Rican and 148 Mexican American children and young adults with asthma. A total of five CpGs exceeded the genome-wide significance threshold of p = 1.17 × 10-7 in the combined analyses from Puerto Ricans and Mexican Americans: cg06035600 (MAP3K6, p = 6.13 × 10-8) showed significant association with pre-bronchodilator Tiffeneau-Pinelli index, the probes cg00914963 (TBC1D16, p = 1.04 × 10-7), cg16405908 (MRGPRE, p = 2.05 × 10-8), and cg07428101 (MUC2, p = 5.02 × 10-9) were associated with post-bronchodilator forced vital capacity (FVC), and cg20515679 (KCNJ6) with post-bronchodilator Tiffeneau-Pinelli index (p = 1.13 × 10-8). However, these markers did not show significant associations in publicly available data from Europeans (p > 0.05). A methylation quantitative trait loci analysis revealed that methylation levels at these CpG sites were regulated by genetic variation in Latinos and the Biobank-based Integrative Omics Studies (BIOS) consortium. Additionally, two differentially methylated regions in REXOC and AURKC were associated with pre-bronchodilator Tiffeneau-Pinelli index (adjusted p < 0.05) in Puerto Ricans and Mexican Americans. Moreover, we replicated some of the previous differentially methylated signals associated with lung function in non-Latino populations. CONCLUSIONS: We replicated previous associations of epigenetic markers with lung function in whole blood and identified novel population-specific associations shared among Latino subgroups.

Pharmacogenetic studies of long-acting beta agonist and inhaled corticosteroid responsiveness in randomised controlled trials of individuals of African descent with asthma.

BACKGROUND: Pharmacogenetic studies in asthma cohorts, primarily made up of White people of European descent, have identified loci associated with response to inhaled beta agonists and corticosteroids (ICSs). Differences exist in how individuals from different ancestral backgrounds respond to long-acting beta agonist (LABA) and ICSs. Therefore, we sought to understand the pharmacogenetic mechanisms regulating therapeutic responsiveness in individuals of African descent. METHODS: We did ancestry-based pharmacogenetic studies of children (aged 5-11 years) and adolescents and adults (aged 12-69 years) from the Best African Response to Drug (BARD) trials, in which participants with asthma uncontrolled with low-dose ICS (fluticasone propionate 50 µg in children, 100 µg in adolescents and adults) received different step-up combination therapies. The hierarchal composite outcome of pairwise superior responsiveness in BARD was based on asthma exacerbations, a 31-day difference in annualised asthma-control days, or a 5% difference in percentage predicted FEV1. We did whole-genome admixture mapping of 15 159 ancestral segments within 312 independent regions, stratified by the two age groups. The two co-primary outcome comparisons were the step up from low-dose ICS to the quintuple dose of ICS (5 × ICS: 250 µg twice daily in children and 500 µg twice daily in adolescents and adults) versus double dose (2-2·5 × ICS: 100 µg twice daily in children, 250 µg twice daily in adolescents and adults), and 5 × ICS versus 100 µg fluticasone plus a LABA (salmeterol 50 µg twice daily). We used a genome-wide significance threshold of p<1·6 × 10-4, and tested for replication using independent cohorts of individuals of African descent with asthma. FINDINGS: We included 249 unrelated children and 267 unrelated adolescents and adults in the BARD pharmacogenetic analysis. In children, we identified a significant admixture mapping peak for superior responsiveness to 5 × ICS versus 100 µg fluticasone plus salmeterol on chromosome 12 (odds ratio [ORlocal African] 3·95, 95% CI 2·02-7·72, p=6·1 × 10-5) fine mapped to a locus adjacent to RNFT2 and NOS1 (rs73399224, ORallele dose 0·17, 95% CI 0·07-0·42, p=8·4 × 10-5). In adolescents and adults, we identified a peak for superior responsiveness to 5 × ICS versus 2·5 × ICS on chromosome 22 (ORlocal African 3·35, 1·98-5·67, p=6·8 × 10-6) containing a locus adjacent to TPST2 (rs5752429, ORallele dose 0·21, 0·09-0·52, p=5·7 × 10-4). We replicated rs5752429 and nominally replicated rs73399224 in independent African American cohorts. INTERPRETATION: BARD is the first genome-wide pharmacogenetic study of LABA and ICS response in clinical trials of individuals of African descent to detect and replicate genome-wide significant loci. Admixture mapping of the composite BARD trial outcome enabled the identification of novel pharmacogenetic variation accounting for differential therapeutic responses in people of African descent with asthma. FUNDING: National Institutes of Health, National Heart, Lung, and Blood Institute.

Lymph node-resident dendritic cells drive TH2 cell development involving MARCH1.

Type 2 T helper (TH2) cells are protective against parasitic worm infections but also aggravate allergic inflammation. Although the role of dendritic cells (DCs) in TH2 cell differentiation is well established, the underlying mechanisms are largely unknown. Here, we show that DC induction of TH2 cells depends on membrane-associated RING-CH-1 (MARCH1) ubiquitin ligase. The pro-TH2 effect of MARCH1 relied on lymph node (LN)­resident DCs, which triggered T cell receptor (TCR) signaling and induced GATA-3 expression from naïve CD4+ T cells independent of tissue-driven migratory DCs. Mice with mutations in the ubiquitin acceptor sites of MHCII and CD86, the two substrates of MARCH1, failed to develop TH2 cells. These findings suggest that TH2 cell development depends on ubiquitin-mediated clearance of antigen-presenting and costimulatory molecules by LN-resident DCs and consequent control of TCR signaling.

Racial/ethnic differences in eligibility for asthma biologics among pediatric populations.

BACKGROUND: Asthma is a heterogeneous disease. Clinical blood parameters differ by race/ethnicity and are used to distinguish asthma subtypes and inform therapies. Differences in subtypes may explain population-specific trends in asthma outcomes. However, these differences in racial/ethnic minority pediatric populations are unclear. OBJECTIVE: We investigated the association of blood parameters and asthma subtypes with asthma outcomes and examined population-specific eligibility for biologic therapies in minority pediatric populations. METHODS: Using data from 2 asthma case-control studies of pediatric minority populations, we performed case-control (N = 3738) and case-only (N = 2743) logistic regressions to quantify the association of blood parameters and asthma subtypes with asthma outcomes. Heterogeneity of these associations was tested using an interaction term between race/ethnicity and each exposure. Differences in therapeutic eligibility were investigated using chi-square tests. RESULTS: Race/ethnicity modified the association between total IgE and asthma exacerbations. Elevated IgE level was associated with worse asthma outcomes in Puerto Ricans. Allergic asthma was associated with worse outcomes in Mexican Americans, whereas eosinophilic asthma was associated with worse outcomes in Puerto Ricans. A lower proportion of Puerto Ricans met dosing criteria for allergic asthma-directed biologic therapy than other groups. A higher proportion of Puerto Ricans qualified for eosinophilic asthma-directed biologic therapy than African Americans. CONCLUSIONS: We found population-specific associations between blood parameters and asthma subtypes with asthma outcomes. Our findings suggest that eligibility for asthma biologic therapies differs across pediatric racial/ethnic populations. These findings call for more studies in diverse populations for equitable treatment of minority patients with asthma.

Robust, flexible, and scalable tests for Hardy-Weinberg equilibrium across diverse ancestries.

Traditional Hardy-Weinberg equilibrium (HWE) tests (the χ2 test and the exact test) have long been used as a metric for evaluating genotype quality, as technical artifacts leading to incorrect genotype calls often can be identified as deviations from HWE. However, in data sets composed of individuals from diverse ancestries, HWE can be violated even without genotyping error, complicating the use of HWE testing to assess genotype data quality. In this manuscript, we present the Robust Unified Test for HWE (RUTH) to test for HWE while accounting for population structure and genotype uncertainty, and to evaluate the impact of population heterogeneity and genotype uncertainty on the standard HWE tests and alternative methods using simulated and real sequence data sets. Our results demonstrate that ignoring population structure or genotype uncertainty in HWE tests can inflate false-positive rates by many orders of magnitude. Our evaluations demonstrate different tradeoffs between false positives and statistical power across the methods, with RUTH consistently among the best across all evaluations. RUTH is implemented as a practical and scalable software tool to rapidly perform HWE tests across millions of markers and hundreds of thousands of individuals while supporting standard VCF/BCF formats. RUTH is publicly available at https://www.github.com/statgen/ruth.

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.

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.

Mapping the 17q12-21.1 Locus for Variants Associated with Early-Onset Asthma in African Americans.

Rationale: The 17q12-21.1 locus is one of the most highly replicated genetic associations with asthma. Individuals of African descent have lower linkage disequilibrium in this region, which could facilitate identifying causal variants.Objectives: To identify functional variants at 17q12-21.1 associated with early-onset asthma among African American individuals.Methods: We evaluated African American participants from SAPPHIRE (Study of Asthma Phenotypes and Pharmacogenomic Interactions by Race-Ethnicity) (n = 1,940), SAGE II (Study of African Americans, Asthma, Genes and Environment) (n = 885), and GCPD-A (Study of the Genetic Causes of Complex Pediatric Disorders-Asthma) (n = 2,805). Associations with asthma onset at ages under 5 years were meta-analyzed across cohorts. The lead signal was reevaluated considering haplotypes informed by genetic ancestry (i.e., African vs. European). Both an expression-quantitative trait locus analysis and a phenome-wide association study were performed on the lead variant.Measurements and Main Results: The meta-analyzed results from SAPPHIRE, SAGE II, and the GCPD-A identified rs11078928 as the top association for early-onset asthma. A haplotype analysis suggested that the asthma association partitioned most closely with the rs11078928 genotype. Genetic ancestry did not appear to influence the effect of this variant. In the expression-quantitative trait locus analysis, rs11078928 was related to alternative splicing of GSDMB (gasdermin-B) transcripts. The phenome-wide association study of rs11078928 suggested that this variant was predominantly associated with asthma and asthma-associated symptoms.Conclusions: A splice-acceptor polymorphism appears to be a causal variant for asthma at the 17q12-21.1 locus. This variant appears to have the same magnitude of effect in individuals of African and European descent.

Genome-wide association study reveals a novel locus for asthma with severe exacerbations in diverse populations.

BACKGROUND: Severe asthma exacerbations are a major cause of asthma morbidity and increased healthcare costs. Several studies have shown racial and ethnic differences in asthma exacerbation rates. We aimed to identify genetic variants associated with severe exacerbations in two high-risk populations for asthma. METHODS: A genome-wide association study of asthma in children and youth with severe exacerbations was performed in 1283 exacerbators and 2027 controls without asthma of Latino ancestry. Independent suggestive variants (P ≤ 5 × 10-6 ) were selected for replication in 448 African Americans exacerbators and 595 controls. Case-only analyses were performed comparing the exacerbators with additional 898 Latinos and 524 African Americans asthma patients without exacerbations, while adjusting by treatment category as a proxy of asthma severity. We analyzed the functionality of associated variants with in silico methods and by correlating genotypes with methylation levels in whole blood in a subset of 473 Latinos. RESULTS: We identified two genome-wide significant associations for susceptibility to asthma with severe exacerbations, including a novel locus located at chromosome 2p21 (rs4952375, odds ratio = 1.39, P = 3.8 × 10-8 ), which was also associated with asthma exacerbations in a case-only analysis (odds ratio = 1.25, P = 1.95 × 10-3 ). This polymorphism is an expression quantitative trait locus of the long intergenic non-protein coding RNA 1913 (LINC01913) in lung tissues (P = 1.3 × 10-7 ) and influences methylation levels of the protein kinase domain-containing cytoplasmic (PKDCC) gene in whole-blood cells (P = 9.8 × 10-5 ). CONCLUSION: We identified a novel susceptibility locus for severe asthma exacerbations in Hispanic/Latino and African American youths with functional effects in gene expression and methylation status of neighboring genes.

Asthma and its relationship to mitochondrial copy number: Results from the Asthma Translational Genomics Collaborative (ATGC) of the Trans-Omics for Precision Medicine (TOPMed) program.

BACKGROUND: Mitochondria support critical cellular functions, such as energy production through oxidative phosphorylation, regulation of reactive oxygen species, apoptosis, and calcium homeostasis. OBJECTIVE: Given the heightened level of cellular activity in patients with asthma, we sought to determine whether mitochondrial DNA (mtDNA) copy number measured in peripheral blood differed between individuals with and without asthma. METHODS: Whole genome sequence data was generated as part of the Trans-Omics for Precision Medicine (TOPMed) Program on participants from the Study of Asthma Phenotypes and Pharmacogenomic Interactions by Race-ethnicity (SAPPHIRE) and the Study of African Americans, Asthma, Genes, & Environment II (SAGE II). We restricted our analysis to individuals who self-identified as African American (3,651 asthma cases and 1,344 controls). Mitochondrial copy number was estimated using the sequencing read depth ratio for the mitochondrial and nuclear genomes. Respiratory complex expression was assessed using RNA-sequencing. RESULTS: Average mitochondrial copy number was significantly higher among individuals with asthma when compared with controls (SAPPHIRE: 218.60 vs. 200.47, P<0.001; SAGE II: 235.99 vs. 223.07, P<0.001). Asthma status was significantly associated with mitochondrial copy number after accounting for potential explanatory variables, such as participant age, sex, leukocyte counts, and mitochondrial haplogroup. Despite the consistent relationship between asthma status and mitochondrial copy number, the latter was not associated with time-to-exacerbation or patient-reported asthma control. Mitochondrial respiratory complex gene expression was disproportionately lower in individuals with asthma when compared with individuals without asthma and other protein-encoding genes. CONCLUSIONS: We observed a robust association between asthma and higher mitochondrial copy number. Asthma having an effect on mitochondria function was also supported by lower respiratory complex gene expression in this group.

Type 2 and interferon inflammation regulate SARS-CoV-2 entry factor expression in the airway epithelium.

Coronavirus disease 2019 (COVID-19) is caused by SARS-CoV-2, an emerging virus that utilizes host proteins ACE2 and TMPRSS2 as entry factors. Understanding the factors affecting the pattern and levels of expression of these genes is important for deeper understanding of SARS-CoV-2 tropism and pathogenesis. Here we explore the role of genetics and co-expression networks in regulating these genes in the airway, through the analysis of nasal airway transcriptome data from 695 children. We identify expression quantitative trait loci for both ACE2 and TMPRSS2, that vary in frequency across world populations. We find TMPRSS2 is part of a mucus secretory network, highly upregulated by type 2 (T2) inflammation through the action of interleukin-13, and that the interferon response to respiratory viruses highly upregulates ACE2 expression. IL-13 and virus infection mediated effects on ACE2 expression were also observed at the protein level in the airway epithelium. Finally, we define airway responses to common coronavirus infections in children, finding that these infections generate host responses similar to other viral species, including upregulation of IL6 and ACE2. Our results reveal possible mechanisms influencing SARS-CoV-2 infectivity and COVID-19 clinical outcomes.