Single-cell resolution of human airway epithelial cells exposed to bronchiolitis obliterans-associated chemicals.

Bronchiolitis obliterans (BO) is a fibrotic lung disease characterized by progressive luminal narrowing and obliteration of the small airways. In the nontransplant population, inhalation exposure to certain chemicals is associated with BO; however, the mechanisms contributing to disease induction remain poorly understood. This study's objective was to use single-cell RNA sequencing for the identification of transcriptomic signatures common to primary human airway epithelial cells after chemical exposure to BO-associated chemicals-diacetyl or nitrogen mustard-to help explain BO induction. Primary airway epithelial cells were cultured at air-liquid interface and exposed to diacetyl, nitrogen mustard, or control vapors. Cultures were dissociated and sequenced for single-cell RNA. Differential gene expression and functional pathway analyses were compared across exposures. In total, 75,663 single cells were captured and sequenced from all exposure conditions. Unbiased clustering identified 11 discrete phenotypes, including 5 basal, 2 ciliated, and 2 secretory cell clusters. With chemical exposure, the proportion of cells assigned to keratin 5+ basal cells decreased, whereas the proportion of cells aligned to secretory cell clusters increased compared with control exposures. Functional pathway analysis identified interferon signaling and antigen processing/presentation as pathways commonly upregulated after diacetyl or nitrogen mustard exposure in a ciliated cell cluster. Conversely, the response of airway basal cells differed significantly with upregulation of the unfolded protein response in diacetyl-exposed basal cells, not seen in nitrogen mustard-exposed cultures. These new insights provide early identification of airway epithelial signatures common to BO-associated chemical exposures.NEW & NOTEWORTHY Bronchiolitis obliterans (BO) is a devastating fibrotic lung disease of the small airways, or bronchioles. This original manuscript uses single-cell RNA sequencing for identifying common signatures of chemically exposed airway epithelial cells in BO induction. Chemical exposure reduced the proportion of keratin 5+ basal cells while increasing the proportion of keratin 4+ suprabasal cells. Functional pathways contributory to these shifts differed significantly across exposures. These new results highlight similarities and differences in BO induction across exposures.

A systems genomics approach uncovers molecular associates of RSV severity.

Respiratory syncytial virus (RSV) infection results in millions of hospitalizations and thousands of deaths each year. Variations in the adaptive and innate immune response appear to be associated with RSV severity. To investigate the host response to RSV infection in infants, we performed a systems-level study of RSV pathophysiology, incorporating high-throughput measurements of the peripheral innate and adaptive immune systems and the airway epithelium and microbiota. We implemented a novel multi-omic data integration method based on multilayered principal component analysis, penalized regression, and feature weight back-propagation, which enabled us to identify cellular pathways associated with RSV severity. In both airway and immune cells, we found an association between RSV severity and activation of pathways controlling Th17 and acute phase response signaling, as well as inhibition of B cell receptor signaling. Dysregulation of both the humoral and mucosal response to RSV may play a critical role in determining illness severity.

Airway gene-expression classifiers for respiratory syncytial virus (RSV) disease severity in infants.

BACKGROUND: A substantial number of infants infected with RSV develop severe symptoms requiring hospitalization. We currently lack accurate biomarkers that are associated with severe illness. METHOD: We defined airway gene expression profiles based on RNA sequencing from nasal brush samples from 106 full-tem previously healthy RSV infected subjects during acute infection (day 1-10 of illness) and convalescence stage (day 28 of illness). All subjects were assigned a clinical illness severity score (GRSS). Using AIC-based model selection, we built a sparse linear correlate of GRSS based on 41 genes (NGSS1). We also built an alternate model based upon 13 genes associated with severe infection acutely but displaying stable expression over time (NGSS2). RESULTS: NGSS1 is strongly correlated with the disease severity, demonstrating a naïve correlation (ρ) of ρ = 0.935 and cross-validated correlation of 0.813. As a binary classifier (mild versus severe), NGSS1 correctly classifies disease severity in 89.6% of the subjects following cross-validation. NGSS2 has slightly less, but comparable, accuracy with a cross-validated correlation of 0.741 and classification accuracy of 84.0%. CONCLUSION: Airway gene expression patterns, obtained following a minimally-invasive procedure, have potential utility for development of clinically useful biomarkers that correlate with disease severity in primary RSV infection.

Airway Gene Expression Correlates of Respiratory Syncytial Virus Disease Severity and Microbiome Composition in Infants.

BACKGROUND: Respiratory syncytial virus (RSV) is the leading cause of severe respiratory disease in infants. The causes and correlates of severe illness in the majority of infants are poorly defined. METHODS: We recruited a cohort of RSV-infected infants and simultaneously assayed the molecular status of their airways and the presence of airway microbiota. We used rigorous statistical approaches to identify gene expression patterns associated with disease severity and microbiota composition, separately and in combination. RESULTS: We measured comprehensive airway gene expression patterns in 106 infants with primary RSV infection. We identified an airway gene expression signature of severe illness dominated by excessive chemokine expression. We also found an association between Haemophilus influenzae, disease severity, and airway lymphocyte accumulation. Exploring the time of onset of clinical symptoms revealed acute activation of interferon signaling following RSV infection in infants with mild or moderate illness, which was absent in subjects with severe illness. CONCLUSIONS: Our data reveal that airway gene expression patterns distinguish mild/moderate from severe illness. Furthermore, our data identify biomarkers that may be therapeutic targets or useful for measuring efficacy of intervention responses.

CX3CR1 as a respiratory syncytial virus receptor in pediatric human lung.

BACKGROUND: Data on the host factors that contribute to infection of young children by respiratory syncytial virus (RSV) are limited. The human chemokine receptor, CX3CR1, has recently been implicated as an RSV receptor. Here we evaluate a role for CX3CR1 in pediatric lung RSV infections. METHODS: CX3CR1 transcript levels in the upper and lower pediatric airways were assessed. Tissue localization and cell-specific expression was confirmed using in situ hybridization and immunohistochemistry. The role of CX3CR1 in RSV infection was also investigated using a novel physiological model of pediatric epithelial cells. RESULTS: Low levels of CX3CR1 transcript were often, but not always, expressed in both upper (62%) and lower airways (36%) of pediatric subjects. CX3CR1 transcript and protein expression was detected in epithelial cells of normal human pediatric lung tissues. CX3CR1 expression was readily detected on primary cultures of differentiated pediatric/infant human lung epithelial cells. RSV demonstrated preferential infection of CX3CR1-positive cells, and blocking CX3CR1/RSV interaction significantly decreased viral load. CONCLUSION: CX3CR1 is present in the airways of pediatric subjects where it may serve as a receptor for RSV infection. Furthermore, CX3CR1 appears to play a mechanistic role in mediating viral infection of pediatric airway epithelial cells in vitro.

Microbiome-Transcriptome Interactions Related to Severity of Respiratory Syncytial Virus Infection.

Respiratory syncytial virus (RSV) is a major cause of lower respiratory tract infections and hospital visits during infancy and childhood. Although risk factors for RSV infection have been identified, the role of microbial species in the respiratory tract is only partially known. We aimed to understand the impact of interactions between the nasal microbiome and host transcriptome on the severity and clinical outcomes of RSV infection. We used 16 S rRNA sequencing to characterize the nasal microbiome of infants with RSV infection. We used RNA sequencing to interrogate the transcriptome of CD4+ T cells obtained from the same set of infants. After dimension reduction through principal component (PC) analysis, we performed an integrative analysis to identify significant co-variation between microbial clade and gene expression PCs. We then employed LIONESS (Linear Interpolation to Obtain Network Estimates for Single Samples) to estimate the clade-gene association patterns for each infant. Our network-based integrative analysis identified several clade-gene associations significantly related to the severity of RSV infection. The microbial taxa with the highest loadings in the implicated clade PCs included Moraxella, Corynebacterium, Streptococcus, Haemophilus influenzae, and Staphylococcus. Interestingly, many of the genes with the highest loadings in the implicated gene PCs are encoded in mitochondrial DNA, while others are involved in the host immune response. This study on microbiome-transcriptome interactions provides insights into how the host immune system mounts a response against RSV and specific infectious agents in nasal microbiota.

The Healthy Infant Nasal Transcriptome: A Benchmark Study.

Responses by resident cells are likely to play a key role in determining the severity of respiratory disease. However, sampling of the airways poses a significant challenge, particularly in infants and children. Here, we report a reliable method for obtaining nasal epithelial cell RNA from infants for genome-wide transcriptomic analysis, and describe baseline expression characteristics in an asymptomatic cohort. Nasal epithelial cells were collected by brushing of the inferior turbinates, and gene expression was interrogated by RNA-seq analysis. Reliable recovery of RNA occurred in the absence of adverse events. We observed high expression of epithelial cell markers and similarity to the transcriptome for intrapulmonary airway epithelial cells. We identified genes displaying low and high expression variability, both inherently, and in response to environmental exposures. The greatest gene expression differences in this asymptomatic cohort were associated with the presence of known pathogenic viruses and/or bacteria. Robust bacteria-associated gene expression patterns were significantly associated with the presence of Moraxella. In summary, we have developed a reliable method for interrogating the infant airway transcriptome by sampling the nasal epithelium. Our data demonstrates both the fidelity and feasibility of our methodology, and describes normal gene expression and variation within a healthy infant cohort.

MicroRNA-494 is a master epigenetic regulator of multiple invasion-suppressor microRNAs by targeting ten eleven translocation 1 in invasive human hepatocellular carcinoma tumors.

UNLABELLED: Vascular invasion provides a direct route for tumor metastasis. The degree to which microRNA (miRNA) expression plays a role in tumor vascular invasion is unclear. Here, we report that miR-494 is up-regulated in human hepatocellular carcinoma (HCC) tumors with vascular invasion and can promote HCC cell invasiveness by gene inactivation of multiple invasion-suppressor miRNAs. Our results show that ten eleven translocation (TET) methylcytosine dioxygenase, predominantly TET1 in HCC cells, is a direct target of miR-494. The reduced 5'-hydroxymethylcytosine levels observed in the proximal cytosine-phosphate-guanine (CpG) regions of multiple invasion-suppressor miRNA genes are strongly associated with their transcriptional repression upon miR-494 overexpression, whereas enforced DNA demethylation can abolish the repression. Furthermore, TET1 knockdown shows a similar effect as miR-494 overexpression. Conversely, miR-494 inhibition or enforced TET1 expression is able to restore invasion-suppressor miRNAs and inhibit miR-494-mediated HCC cell invasion. CONCLUSIONS: miR-494 can trigger gene silencing of multiple invasion-suppressor miRNAs by inhibiting genomic DNA demethylation by direct targeting of TET1, thereby leading to tumor vascular invasion.

The genome-wide transcriptional response to neonatal hyperoxia identifies Ahr as a key regulator.

Premature infants requiring supplemental oxygen are at increased risk for developing bronchopulmonary dysplasia (BPD). Rodent models involving neonatal exposure to excessive oxygen concentrations (hyperoxia) have helped to identify mechanisms of BPD-associated pathology. Genome-wide assessments of the effects of hyperoxia in neonatal mouse lungs could identify novel BPD-related genes and pathways. Newborn C57BL/6 mice were exposed to 100% oxygen for 10 days, and whole lung tissue RNA was used for high-throughput, sequencing-based transcriptomic analysis (RNA-Seq). Significance Analysis of Microarrays and Ingenuity Pathway Analysis were used to identify genes and pathways affected. Expression patterns for selected genes were validated by qPCR. Mechanistic relationships between genes were further tested in cultured mouse lung epithelial cells. We identified 300 genes significantly and substantially affected following acute neonatal hyperoxia. Canonical pathways dysregulated in hyperoxia lungs included nuclear factor (erythryoid-derived-2)-like 2-mediated oxidative stress signaling, p53 signaling, eNOS signaling, and aryl hydrocarbon receptor (Ahr) pathways. Cluster analysis identified Ccnd1, Cdkn1a, and Ahr as critical regulatory nodes in the response to hyperoxia, with Ahr serving as the major effector node. A mechanistic role for Ahr was assessed in lung epithelial cells, and we confirmed its ability to regulate the expression of multiple hyperoxia markers, including Cdkn1a, Pdgfrb, and A2m. We conclude that a global assessment of gene regulation in the acute neonatal hyperoxia model of BPD-like pathology has identified Ahr as one driver of gene dysregulation.

Neutropenia with impaired host defense against microbial infection in mice lacking androgen receptor.

Neutrophils, the major phagocytes that form the first line of cell-mediated defense against microbial infection, are produced in the bone marrow and released into the circulation in response to granulocyte-colony stimulating factor (G-CSF). Here, we report that androgen receptor knockout (ARKO) mice are neutropenic and susceptible to acute bacterial infection, whereas castration only results in moderate neutrophil reduction in mice and humans. Androgen supplement can restore neutrophil counts via stabilizing AR in castrated mice, but not in ARKO and testicular feminization mutant (Tfm) mice. Our results show that deletion of the AR gene does not influence myeloid lineage commitment, but significantly reduces the proliferative activity of neutrophil precursors and retards neutrophil maturation. CXCR2-dependent migration is also decreased in ARKO neutrophils as compared with wild-type controls. G-CSF is unable to delay apoptosis in ARKO neutrophils, and ARKO mice show a poor granulopoietic response to exogenous G-CSF injection. In addition, AR can restore G-CSF-dependent granulocytic differentiation upon transduction into ARKO progenitors. We further found that AR augments G-CSF signaling by activating extracellular signal-regulated kinase 1/2 and also by sustaining Stat3 activity via diminishing the inhibitory binding of PIAS3 to Stat3. Collectively, our findings demonstrate an essential role for AR in granulopoiesis and host defense against microbial infection.

Big results from small samples: evaluation of amplification protocols for gene expression profiling.

Microarrays have revolutionized many areas of biology due to our technical ability to quantify tens of thousands of transcripts within a single experiment. However, there are still many areas that cannot benefit from this technology due to the amount of biological material needed for microarray analysis. In response to this demand, chemistries have been developed that boast the capability of generating targets from nanogram amounts of total RnA, reflecting minimal amounts of biological material, on the order of several hundred or thousand cells. Herein, we describe the evaluation of four chemistries for RnA amplification in terms of reproducibility, sensitivity, accuracy, and comparability to results from a single round of T7 amplification. No evidence for false-positive measurements of differential expression was observed. In contrast, clear differences between chemistries in sensitivity and accuracy were detected. PCR validation showed an interaction of probe sequence on the array and target labeling chemistry, resulting in a chemistry-dependent probe set sensitivity varying over an order of magnitude.

9-cis-retinoic acid inhibits androgen receptor activity through activation of retinoid X receptor.

Although the retinoic X receptor (RXR) forms heterodimers with many members of the estrogen receptor subfamily, the interaction between RXR and the members of the glucocorticoid receptor subfamily remains unclear. Here we show that the RXR can form a heterodimer with the androgen receptor (AR) under in vitro and in vivo conditions. Functional analyses further demonstrated that the AR, in the presence or absence of androgen, can function as a repressor to suppress RXR target genes, thereby preventing the RXR binding to the RXR DNA response element. In contrast, RXR can function as a repressor to suppress AR target genes in the presence of 9-cis-retinoic acid, but unliganded RXR can function as a weak coactivator to moderately enhance AR transactivation. Together, these results not only reveal a unique interaction between members of the two nuclear receptor subfamilies, but also represent the first evidence showing a nuclear receptor (RXR) may function as either a repressor or a coactivator based on the ligand binding status.