Transmission Dynamics and Rare Clustered Transmission Within an Urban University Population Before Widespread Vaccination.

BACKGROUND: Universities returned to in-person learning in 2021 while SARS-CoV-2 spread remained high. At the time, it was not clear whether in-person learning would be a source of disease spread. METHODS: We combined surveillance testing, universal contact tracing, and viral genome sequencing to quantify introductions and identify likely on-campus spread. RESULTS: Ninety-one percent of viral genotypes occurred once, indicating no follow-on transmission. Less than 5% of introductions resulted in >3 cases, with 2 notable exceptions of 40 and 47 cases. Both partially overlapped with outbreaks defined by contact tracing. In both cases, viral genomics eliminated over half the epidemiologically linked cases but added an equivalent or greater number of individuals to the transmission cluster. CONCLUSIONS: Public health interventions prevented within-university transmission for most SARS-CoV-2 introductions, with only 2 major outbreaks being identified January to May 2021. The genetically linked cases overlap with outbreaks identified by contact tracing; however, they persisted in the university population for fewer days and rounds of transmission than estimated via contact tracing. This underscores the effectiveness of test-trace-isolate strategies in controlling undetected spread of emerging respiratory infectious diseases. These approaches limit follow-on transmission in both outside-in and internal transmission conditions.

Bronchial gene expression alterations associated with radiological bronchiectasis.

OBJECTIVES: Discovering airway gene expression alterations associated with radiological bronchiectasis may improve the understanding of the pathobiology of early-stage bronchiectasis. METHODS: Presence of radiological bronchiectasis in 173 individuals without a clinical diagnosis of bronchiectasis was evaluated. Bronchial brushings from these individuals were transcriptomically profiled and analysed. Single-cell deconvolution was performed to estimate changes in cellular landscape that may be associated with early disease progression. RESULTS: 20 participants have widespread radiological bronchiectasis (three or more lobes). Transcriptomic analysis reflects biological processes associated with bronchiectasis including decreased expression of genes involved in cell adhesion and increased expression of genes involved in inflammatory pathways (655 genes, false discovery rate <0.1, log2 fold-change >0.25). Deconvolution analysis suggests that radiological bronchiectasis is associated with an increased proportion of ciliated and deuterosomal cells, and a decreased proportion of basal cells. Gene expression patterns separated participants into three clusters: normal, intermediate and bronchiectatic. The bronchiectatic cluster was enriched by participants with more lobes of radiological bronchiectasis (p<0.0001), more symptoms (p=0.002), higher SERPINA1 mutation rates (p=0.03) and higher computed tomography derived bronchiectasis scores (p<0.0001). CONCLUSIONS: Genes involved in cell adhesion, Wnt signalling, ciliogenesis and interferon-γ pathways had altered expression in the bronchus of participants with widespread radiological bronchiectasis, possibly associated with decreased basal and increased ciliated cells. This gene expression pattern is not only highly enriched among individuals with radiological bronchiectasis, but also associated with airway-related symptoms in those without discernible radiological bronchiectasis, suggesting that it reflects a bronchiectasis-associated, but non-bronchiectasis-specific lung pathophysiological process.

Bronchial gene expression signature associated with rate of subsequent FEV1 decline in individuals with and at risk of COPD.

BACKGROUND: COPD is characterised by progressive lung function decline. Leveraging prior work demonstrating bronchial airway COPD-associated gene expression alterations, we sought to determine if there are alterations associated with differences in the rate of FEV1 decline. METHODS: We examined gene expression among ever smokers with and without COPD who at baseline had bronchial brushings profiled by Affymetrix microarrays and had longitudinal lung function measurements (n=134; mean follow-up=6.38±2.48 years). Gene expression profiles associated with the rate of FEV1 decline were identified by linear modelling. RESULTS: Expression differences in 171 genes were associated with rate of FEV1 decline (false discovery rate <0.05). The FEV1 decline signature was replicated in an independent dataset of bronchial biopsies from patients with COPD (n=46; p=0.018; mean follow-up=6.76±1.32 years). Genes elevated in individuals with more rapid FEV1 decline are significantly enriched among the genes altered by modulation of XBP1 in two independent datasets (Gene Set Enrichment Analysis (GSEA) p<0.05) and are enriched in mucin-related genes (GSEA p<0.05). CONCLUSION: We have identified and replicated an airway gene expression signature associated with the rate of FEV1 decline. Aspects of this signature are related to increased expression of XBP1-regulated genes, a transcription factor involved in the unfolded protein response, and genes related to mucin production. Collectively, these data suggest that molecular processes related to the rate of FEV1 decline can be detected in airway epithelium, identify a possible indicator of FEV1 decline and make it possible to detect, in an early phase, ever smokers with and without COPD most at risk of rapid FEV1 decline.

Effect of long-term corticosteroid treatment on microRNA and gene-expression profiles in COPD.

The aim was to investigate whether microRNA (miRNA) expression is modulated by inhaled corticosteroid (ICS) treatmentWe performed genome-wide miRNA analysis on bronchial biopsies of 69 moderate/severe chronic obstructive pulmonary disease (COPD) patients at baseline and after 6- and 30-month treatment with the ICS fluticasone propionate or placebo. The effect of ICS on miRNA expression was validated in differentiated primary bronchial epithelial cultures, and functional studies were conducted in BEAS-2B cells. MiRNAs affected by ICS and their predicted targets were compared to an independent miRNA dataset of bronchial brushings from COPD patients and healthy controls.Treatment with ICS for both 6 and 30 months significantly altered the expression of four miRNAs, including miR-320d, which was increased during ICS treatment compared with placebo. The ICS-induced increase of miR-320d was confirmed in primary airway epithelial cells. MiR-320d negatively correlated targets were enriched for pro-inflammatory genes and were increased in the bronchial brushes of patients with lower lung function in the independent dataset. Overexpression of miR-320d in BEAS-2B cells dampened cigarette smoke extract-induced pro-inflammatory activity via inhibition of nuclear factor-κB.Collectively, we identified miR-320d as a novel mediator of ICS, regulating the pro-inflammatory response of the airway epithelium.

Cross-site comparison of ribosomal depletion kits for Illumina RNAseq library construction.

BACKGROUND: Ribosomal RNA (rRNA) comprises at least 90% of total RNA extracted from mammalian tissue or cell line samples. Informative transcriptional profiling using massively parallel sequencing technologies requires either enrichment of mature poly-adenylated transcripts or targeted depletion of the rRNA fraction. The latter method is of particular interest because it is compatible with degraded samples such as those extracted from FFPE and also captures transcripts that are not poly-adenylated such as some non-coding RNAs. Here we provide a cross-site study that evaluates the performance of ribosomal RNA removal kits from Illumina, Takara/Clontech, Kapa Biosystems, Lexogen, New England Biolabs and Qiagen on intact and degraded RNA samples. RESULTS: We find that all of the kits are capable of performing significant ribosomal depletion, though there are differences in their ease of use. All kits were able to remove ribosomal RNA to below 20% with intact RNA and identify ~ 14,000 protein coding genes from the Universal Human Reference RNA sample at >1FPKM. Analysis of differentially detected genes between kits suggests that transcript length may be a key factor in library production efficiency. CONCLUSIONS: These results provide a roadmap for labs on the strengths of each of these methods and how best to utilize them.

Assessment of microRNA differential expression and detection in multiplexed small RNA sequencing data.

Small RNA sequencing can be used to gain an unprecedented amount of detail into the microRNA transcriptome. The relatively high cost and low throughput of sequencing bases technologies can potentially be offset by the use of multiplexing. However, multiplexing involves a trade-off between increased number of sequenced samples and reduced number of reads per sample (i.e., lower depth of coverage). To assess the effect of different sequencing depths owing to multiplexing on microRNA differential expression and detection, we sequenced the small RNA of lung tissue samples collected in a clinical setting by multiplexing one, three, six, nine, or 12 samples per lane using the Illumina HiSeq 2000. As expected, the numbers of reads obtained per sample decreased as the number of samples in a multiplex increased. Furthermore, after normalization, replicate samples included in distinct multiplexes were highly correlated (R > 0.97). When detecting differential microRNA expression between groups of samples, microRNAs with average expression >1 reads per million (RPM) had reproducible fold change estimates (signal to noise) independent of the degree of multiplexing. The number of microRNAs detected was strongly correlated with the log2 number of reads aligning to microRNA loci (R = 0.96). However, most additional microRNAs detected in samples with greater sequencing depth were in the range of expression which had lower fold change reproducibility. These findings elucidate the trade-off between increasing the number of samples in a multiplex with decreasing sequencing depth and will aid in the design of large-scale clinical studies exploring microRNA expression and its role in disease.

Gene-expression profiling of buccal epithelium among non-smoking women exposed to household air pollution from smoky coal.

In China's rural counties of Xuanwei and Fuyuan, lung cancer rates are among the highest in the world. While the elevated disease risk in this population has been linked to the usage of smoky (bituminous) coal as compared to smokeless (anthracite) coal, the underlying molecular changes associated with this exposure remains unclear. To understand the physiologic effects of smoky coal exposure, we analyzed the genome-wide gene-expression profiles in buccal epithelial cells collected from healthy, non-smoking female residents of Xuanwei and Fuyuan who burn smoky (n = 26) and smokeless (n = 9) coal. Gene-expression was profiled via microarrays, and changes associated with coal type were correlated to household levels of fine particulate matter (PM2.5) and polycyclic aromatic hydrocarbons (PAHs). Expression levels of 282 genes were altered with smoky versus smokeless coal exposure (P < 0.005), including the 2-fold increase of proinflammatory IL8 and decrease of proapoptotic CASP3. This signature was more correlated with carcinogenic PAHs (e.g. Benzo[a]pyrene; r = 0.41) than with non-carcinogenic PAHs (e.g. Fluorene; r = 0.08) or PM2.5 (r = 0.05). Genes altered with smoky coal exposure were concordantly enriched with tobacco exposure in previously profiled buccal biopsies of smokers and non-smokers (GSEA, q < 0.05). This is the first study to identify a signature of buccal epithelial gene-expression that is associated with smoky coal exposure, which in part is similar to the molecular response to tobacco smoke, thereby lending biologic plausibility to prior epidemiological studies that have linked this exposure to lung cancer risk.

Low-dose radiation affects cardiac physiology: gene networks and molecular signaling in cardiomyocytes.

There are 160,000 cancer patients worldwide treated with particle radiotherapy (RT). With the advent of proton, and high (H) charge (Z) and energy (E) HZE ionizing particle RT, the cardiovascular diseases risk estimates are uncertain. In addition, future deep space exploratory-type missions will expose humans to unknown but low doses of particle irradiation (IR). We examined molecular responses using transcriptome profiling in left ventricular murine cardiomyocytes isolated from mice that were exposed to 90 cGy, 1 GeV proton ((1)H) and 15 cGy, 1 GeV/nucleon iron ((56)Fe) over 28 days after exposure. Unsupervised clustering analysis of gene expression segregated samples according to the IR response and time after exposure, with (56)Fe-IR showing the greatest level of gene modulation. (1)H-IR showed little differential transcript modulation. Network analysis categorized the major differentially expressed genes into cell cycle, oxidative responses, and transcriptional regulation functional groups. Transcriptional networks identified key nodes regulating expression. Validation of the signal transduction network by protein analysis and gel shift assay showed that particle IR clearly regulates a long-lived signaling mechanism for ERK1/2, p38 MAPK signaling and identified NFATc4, GATA4, STAT3, and NF-κB as regulators of the response at specific time points. These data suggest that the molecular responses and gene expression to (56)Fe-IR in cardiomyocytes are unique and long-lasting. Our study may have significant implications for the efforts of National Aeronautics and Space Administration to develop heart disease risk estimates for astronauts and for patients receiving conventional and particle RT via identification of specific HZE-IR molecular markers.

Distinct gene signatures in aortic tissue from ApoE-/- mice exposed to pathogens or Western diet.

BACKGROUND: Atherosclerosis is a progressive disease characterized by inflammation and accumulation of lipids in vascular tissue. Porphyromonas gingivalis (Pg) and Chlamydia pneumoniae (Cp) are associated with inflammatory atherosclerosis in humans. Similar to endogenous mediators arising from excessive dietary lipids, these Gram-negative pathogens are pro-atherogenic in animal models, although the specific inflammatory/atherogenic pathways induced by these stimuli are not well defined. In this study, we identified gene expression profiles that characterize P. gingivalis, C. pneumoniae, and Western diet (WD) at acute and chronic time points in aortas of Apolipoprotein E (ApoE-/-) mice. RESULTS: At the chronic time point, we observed that P. gingivalis was associated with a high number of unique differentially expressed genes compared to C. pneumoniae or WD. For the top 500 differentially expressed genes unique to each group, we observed a high percentage (76%) that exhibited decreased expression in P. gingivalis-treated mice in contrast to a high percentage (96%) that exhibited increased expression in WD mice. C. pneumoniae treatment resulted in approximately equal numbers of genes that exhibited increased and decreased expression. Gene Set Enrichment Analysis (GSEA) revealed distinct stimuli-associated phenotypes, including decreased expression of mitochondrion, glucose metabolism, and PPAR pathways in response to P. gingivalis but increased expression of mitochondrion, lipid metabolism, carbohydrate and amino acid metabolism, and PPAR pathways in response to C. pneumoniae; WD was associated with increased expression of immune and inflammatory pathways. DAVID analysis of gene clusters identified by two-way ANOVA at acute and chronic time points revealed a set of core genes that exhibited altered expression during the natural progression of atherosclerosis in ApoE-/- mice; these changes were enhanced in P. gingivalis-treated mice but attenuated in C. pneumoniae-treated mice. Notable differences in the expression of genes associated with unstable plaques were also observed among the three pro-atherogenic stimuli. CONCLUSIONS: Despite the common outcome of P. gingivalis, C. pneumoniae, and WD on the induction of vascular inflammation and atherosclerosis, distinct gene signatures and pathways unique to each pro-atherogenic stimulus were identified. Our results suggest that pathogen exposure results in dysregulated cellular responses that may impact plaque progression and regression pathways.

Airway gene expression in COPD is dynamic with inhaled corticosteroid treatment and reflects biological pathways associated with disease activity.

BACKGROUND: A core feature of chronic obstructive pulmonary disease (COPD) is the accelerated decline in forced expiratory volume in one second (FEV1). The recent Groningen and Leiden Universities study of Corticosteroids in Obstructive Lung Disease (GLUCOLD) study suggested that particular phenotypes of COPD benefit from fluticasone±salmeterol by reducing the rate of FEV1 decline, yet the underlying mechanisms are unknown. METHODS: Whole-genome gene expression profiling using the Affymetrix Gene ST array (V.1.0) was performed on 221 bronchial biopsies available from 89 COPD patients at baseline and after 6 and 30 months of fluticasone±salmeterol and placebo treatment in GLUCOLD. RESULTS: Linear mixed effects modelling revealed that the expression of 138 genes decreased, whereas the expression of 140 genes significantly upregulated after both 6 and 30 months of treatment with fluticasone±salmeterol versus placebo. A more pronounced treatment-induced change in the expression of 50 and 55 of these 278 genes was associated with a lower rate of decline in FEV1 and Saint George Respiratory Questionnaire, respectively. Genes decreasing with treatment were involved in pathways related to cell cycle, oxidative phosphorylation, epithelial cell signalling, p53 signalling and T cell signalling. Genes increasing with treatment were involved in pathways related to focal adhesion, gap junction and extracellular matrix deposition. Finally, the fluticasone-induced gene expression changes were enriched among genes that change in the airway epithelium in smokers with versus without COPD in an independent data set. CONCLUSIONS: The present study suggests that gene expression in biological pathways of COPD is dynamic with treatment and reflects disease activity. This study opens the gate to targeted and molecular phenotype-driven therapy of COPD.

Transcriptional responses of neonatal mouse lung to hyperoxia by Nrf2 status.

UNLABELLED: Hyperoxia exposure can inhibit alveolar growth in the neonatal lung through induction of p21/p53 pathways and is a risk factor for the development of bronchopulmonary dysplasia (BPD) in preterm infants. We previously found that activation of nuclear factor erythroid 2 p45-related factor (Nrf2) improved survival in neonatal mice exposed to hyperoxia likely due to increased expression of anti-oxidant response genes. It is not known however, whether hyperoxic induced Nrf2 activation attenuates the growth impairment caused by hyperoxia in neonatal lung. To determine if Nrf2 activation modulates cell cycle regulatory pathway genes associated with growth arrest we examined the gene expression in the lungs of Nrf2(-/-) and Nrf2(+/+) neonatal mice at one and 3days of hyperoxia exposure. METHODS: Microarray analysis was performed in neonatal Nrf2(+/+) and Nrf2(-/-) lungs exposed to one and 3days of hyperoxia. Sulforaphane, an inducer of Nrf2 was given to timed pregnant mice to determine if in utero exposure attenuated p21 and IL-6 gene expression in wildtype neonatal mice exposed to hyperoxia. RESULTS: Cell cycle regulatory genes were induced in Nrf2(-/-) lung at 1day of hyperoxia. At 3days of hyperoxia, induction of cell cycle regulatory genes was similar in Nrf2(+/+) and Nrf2(-/-) lungs, despite higher inflammatory gene expression in Nrf2(-/-) lung. CONCLUSION: p21/p53 pathways gene expression was not attenuated by Nrf2 activation in neonatal lung. In utero SUL did not attenuate p21 expression in wildtype neonatal lung exposed to hyperoxia. These findings suggest that although Nrf2 activation induces expression of anti-oxidant genes, it does not attenuate alveolar growth arrest caused by exposure to hyperoxia.

A dynamic bronchial airway gene expression signature of chronic obstructive pulmonary disease and lung function impairment.

RATIONALE: Molecular phenotyping of chronic obstructive pulmonary disease (COPD) has been impeded in part by the difficulty in obtaining lung tissue samples from individuals with impaired lung function. OBJECTIVES: We sought to determine whether COPD-associated processes are reflected in gene expression profiles of bronchial airway epithelial cells obtained by bronchoscopy. METHODS: Gene expression profiling of bronchial brushings obtained from 238 current and former smokers with and without COPD was performed using Affymetrix Human Gene 1.0 ST Arrays. MEASUREMENTS AND MAIN RESULTS: We identified 98 genes whose expression levels were associated with COPD status, FEV1% predicted, and FEV1/FVC. In silico analysis identified activating transcription factor 4 (ATF4) as a potential transcriptional regulator of genes with COPD-associated airway expression, and ATF4 overexpression in airway epithelial cells in vitro recapitulates COPD-associated gene expression changes. Genes with COPD-associated expression in the bronchial airway epithelium had similarly altered expression profiles in prior studies performed on small-airway epithelium and lung parenchyma, suggesting that transcriptomic alterations in the bronchial airway epithelium reflect molecular events found at more distal sites of disease activity. Many of the airway COPD-associated gene expression changes revert toward baseline after therapy with the inhaled corticosteroid fluticasone in independent cohorts. CONCLUSIONS: Our findings demonstrate a molecular field of injury throughout the bronchial airway of active and former smokers with COPD that may be driven in part by ATF4 and is modifiable with therapy. Bronchial airway epithelium may ultimately serve as a relatively accessible tissue in which to measure biomarkers of disease activity for guiding clinical management of COPD.

Endothelial Erg Regulates Expression of Pulmonary Lymphatic Junctional and Inflammation Genes in Mouse Lungs Impacting Lymphatic Transport.

The ETS transcription factor ERG is a master regulator of endothelial gene specificity and highly enriched in the capillary, vein, and arterial endothelial cells. ERG expression is critical for endothelial barrier function, permeability, and vascular inflammation. A dysfunctional vascular endothelial ERG has been shown to impair lung capillary homeostasis, contributing to pulmonary fibrosis as previously observed in IPF lungs. Our preliminary observations indicate that lymphatic endothelial cells (LEC) in the human IPF lung also lack ERG. To understand the role of ERG in pulmonary LECs, we developed LEC-specific inducible Erg-CKO and Erg-GFP-CKO conditional knockout (CKO) mice under Prox1 promoter. Whole lung microarray analysis, flow cytometry, and qPCR confirmed an inflammatory and pro-lymphvasculogenic predisposition in Erg-CKO lung. FITC-Dextran tracing analysis showed an increased pulmonary interstitial lymphatic fluid transport from the lung to the axial lymph node. Single-cell transcriptomics confirmed that genes associated with cell junction integrity were downregulated in Erg-CKO pre-collector and collector LECs. Integrating Single-cell transcriptomics and CellChatDB helped identify LEC specific communication pathways contributing to pulmonary inflammation, trans-endothelial migration, inflammation, and Endo-MT in Erg-CKO lung. Our findings suggest that downregulation of lymphatic Erg crucially affects LEC function, LEC permeability, pulmonary LEC communication pathways and lymphatic transcriptomics.

Identification of a gene expression signature of vascular invasion and recurrence in stage I lung adenocarcinoma via bulk and spatial transcriptomics.

Microscopic vascular invasion (VI) is predictive of recurrence and benefit from lobectomy in stage I lung adenocarcinoma (LUAD) but is difficult to assess in resection specimens and cannot be accurately predicted prior to surgery. Thus, new biomarkers are needed to identify this aggressive subset of stage I LUAD tumors. To assess molecular and microenvironment features associated with angioinvasive LUAD we profiled 162 resected stage I tumors with and without VI by RNA-seq and explored spatial patterns of gene expression in a subset of 15 samples by high-resolution spatial transcriptomics (stRNA-seq). Despite the small size of invaded blood vessels, we identified a gene expression signature of VI from the bulk RNA-seq discovery cohort (n=103) and found that it was associated with VI foci, desmoplastic stroma, and high-grade patterns in our stRNA-seq data. We observed a stronger association with high-grade patterns from VI+ compared with VI- tumors. Using the discovery cohort, we developed a transcriptomic predictor of VI, that in an independent validation cohort (n=60) was associated with VI (AUROC=0.86; p=5.42×10-6) and predictive of recurrence-free survival (HR=1.98; p=0.024), even in VI- LUAD (HR=2.76; p=0.003). To determine our VI predictor's robustness to intra-tumor heterogeneity we used RNA-seq data from multi-region sampling of stage I LUAD cases in TRACERx, where the predictor scores showed high correlation (R=0.87, p<2.2×10-16) between two randomly sampled regions of the same tumor. Our study suggests that VI-associated gene expression changes are detectable beyond the site of intravasation and can be used to predict the presence of VI. This may enable the prediction of angioinvasive LUAD from biopsy specimens, allowing for more tailored medical and surgical management of stage I LUAD.

Heat Inactivation of Nipah Virus for Downstream Single-Cell RNA Sequencing Does Not Interfere with Sample Quality.

Single-cell RNA sequencing (scRNA-seq) technologies are instrumental to improving our understanding of virus-host interactions in cell culture infection studies and complex biological systems because they allow separating the transcriptional signatures of infected versus non-infected bystander cells. A drawback of using biosafety level (BSL) 4 pathogens is that protocols are typically developed without consideration of virus inactivation during the procedure. To ensure complete inactivation of virus-containing samples for downstream analyses, an adaptation of the workflow is needed. Focusing on a commercially available microfluidic partitioning scRNA-seq platform to prepare samples for scRNA-seq, we tested various chemical and physical components of the platform for their ability to inactivate Nipah virus (NiV), a BSL-4 pathogen that belongs to the group of nonsegmented negative-sense RNA viruses. The only step of the standard protocol that led to NiV inactivation was a 5 min incubation at 85 °C. To comply with the more stringent biosafety requirements for BSL-4-derived samples, we included an additional heat step after cDNA synthesis. This step alone was sufficient to inactivate NiV-containing samples, adding to the necessary inactivation redundancy. Importantly, the additional heat step did not affect sample quality or downstream scRNA-seq results.

Single nuclei transcriptomics in human and non-human primate striatum in opioid use disorder.

In brain, the striatum is a heterogenous region involved in reward and goal-directed behaviors. Striatal dysfunction is linked to psychiatric disorders, including opioid use disorder (OUD). Striatal subregions are divided based on neuroanatomy, each with unique roles in OUD. In OUD, the dorsal striatum is involved in altered reward processing, formation of habits, and development of negative affect during withdrawal. Using single nuclei RNA-sequencing, we identified both canonical (e.g., dopamine receptor subtype) and less abundant cell populations (e.g., interneurons) in human dorsal striatum. Pathways related to neurodegeneration, interferon response, and DNA damage were significantly enriched in striatal neurons of individuals with OUD. DNA damage markers were also elevated in striatal neurons of opioid-exposed rhesus macaques. Sex-specific molecular differences in glial cell subtypes associated with chronic stress were found in OUD, particularly female individuals. Together, we describe different cell types in human dorsal striatum and identify cell type-specific alterations in OUD.

Integrative genetic and genomic networks identify microRNA associated with COPD and ILD.

Chronic obstructive pulmonary disease (COPD) and interstitial lung disease (ILD) are clinically and molecularly heterogeneous diseases. We utilized clustering and integrative network analyses to elucidate roles for microRNAs (miRNAs) and miRNA isoforms (isomiRs) in COPD and ILD pathogenesis. Short RNA sequencing was performed on 351 lung tissue samples of COPD (n = 145), ILD (n = 144) and controls (n = 64). Five distinct subclusters of samples were identified including 1 COPD-predominant cluster and 2 ILD-predominant clusters which associated with different clinical measurements of disease severity. Utilizing 262 samples with gene expression and SNP microarrays, we built disease-specific genetic and expression networks to predict key miRNA regulators of gene expression. Members of miR-449/34 family, known to promote airway differentiation by repressing the Notch pathway, were among the top connected miRNAs in both COPD and ILD networks. Genes associated with miR-449/34 members in the disease networks were enriched among genes that increase in expression with airway differentiation at an air-liquid interface. A highly expressed isomiR containing a novel seed sequence was identified at the miR-34c-5p locus. 47% of the anticorrelated predicted targets for this isomiR were distinct from the canonical seed sequence for miR-34c-5p. Overexpression of the canonical miR-34c-5p and the miR-34c-5p isomiR with an alternative seed sequence down-regulated NOTCH1 and NOTCH4. However, only overexpression of the isomiR down-regulated genes involved in Ras signaling such as CRKL and GRB2. Overall, these findings elucidate molecular heterogeneity inherent across COPD and ILD patients and further suggest roles for miR-34c in regulating disease-associated gene-expression.

Systematic dissection, preservation, and multiomics in whole human and bovine hearts.

OBJECTIVES: We sought to develop a rigorous, systematic protocol for the dissection and preservation of human hearts for biobanking that expands previous success in postmortem transcriptomics to multiomics from paired tissue. BACKGROUND: Existing cardiac biobanks consist largely of biopsy tissue or explanted hearts in select diseases and are insufficient for correlating whole organ phenotype with clinical data. METHODS: We demonstrate optimal conditions for multiomics interrogation (ribonucleic acid (RNA) sequencing, untargeted metabolomics) in hearts by evaluating the effect of technical variables (storage solution, temperature) and simulated postmortem interval (PMI) on RNA and metabolite stability. We used bovine (n=3) and human (n=2) hearts fixed in PAXgene or snap-frozen with liquid nitrogen. RESULTS: Using a paired Wald test, only two of the genes assessed were differentially expressed between left ventricular samples from bovine hearts stored in PAXgene at 0 and 12 hours PMI (FDR q<0.05). We obtained similar findings in human left ventricular samples, suggesting stability of RNA transcripts at PMIs up to 12 hours. Different library preparation methods (mRNA poly-A capture vs. rRNA depletion) resulted in similar quality metrics with both library preparations achieving >95% of reads properly aligning to the reference genomes across all PMIs for bovine and human hearts. PMI had no effect on RNA Integrity Number or quantity of RNA recovered at the time points evaluated. Of the metabolites identified (855 total) using untargeted metabolomics of human left ventricular tissue, 503 metabolites remained stable across PMIs (0, 4, 8, 12 hours). Most metabolic pathways retained several stable metabolites. CONCLUSIONS: Our data demonstrate a technically rigorous, reproducible protocol that will enhance cardiac biobanking practices and facilitate novel insights into human CVD. CONDENSED ABSTRACT: Cardiovascular disease (CVD) is the leading cause of mortality worldwide. Current biobanking practices insufficiently capture both the diverse array of phenotypes present in CVDs and the spatial heterogeneity across cardiac tissue sites. We have developed a rigorous and systematic protocol for the dissection and preservation of human cardiac biospecimens to enhance the availability of whole organ tissue for multiple applications. When combined with longitudinal clinical phenotyping, our protocol will enable multiomics in hearts to deepen our understanding of CVDs.

Clinical Study of Aspirin and Zileuton on Biomarkers of Tobacco-Related Carcinogenesis in Current Smokers.

The chemopreventive effect of aspirin and other non-steroidal anti-inflammatory drugs (NSAIDs) on lung cancer risk is supported by epidemiologic and preclinical studies. Zileuton, a 5-lipoxygenase inhibitor, has additive activity with NSAIDs against tobacco carcinogenesis in preclinical models. We hypothesized that cyclooxygenase plus 5-lipoxygenase inhibition would be more effective than a placebo in modulating the nasal epithelium gene signatures of tobacco exposure and lung cancer. We conducted a randomized, double-blinded study of low-dose aspirin plus zileuton vs. double placebo in current smokers to compare the modulating effects on nasal gene expression and arachidonic acid metabolism. In total, 63 participants took aspirin 81 mg daily plus zileuton (Zyflo CR) 600 mg BID or the placebo for 12 weeks. Nasal brushes from the baseline, end-of-intervention, and one-week post intervention were profiled via microarray. Aspirin plus zilueton had minimal effects on the modulation of the nasal or bronchial gene expression signatures of smoking, lung cancer, and COPD but favorably modulated a bronchial gene expression signature of squamous dysplasia. Aspirin plus zileuton suppressed urinary leukotriene but not prostaglandin E2, suggesting shunting through the cyclooxygenase pathway when combined with 5-lipoxygenase inhibition. Continued investigation of leukotriene inhibitors is needed to confirm these findings, understand the long-term effects on the airway epithelium, and identify the safest, optimally dosed agents.

Temporal and Quantitative Transcriptomic Differences Define Sexual Dimorphism in Murine Postnatal Bone Aging.

Time is a central element of the sexual dimorphic patterns of development, pathology, and aging of the skeleton. Because the transcriptome is a representation of the phenome, we hypothesized that both sex and sex-specific temporal, transcriptomic differences in bone tissues over an 18-month period would be informative to the underlying molecular processes that lead to postnatal sexual dimorphism. Regardless of age, sex-associated changes of the whole bone transcriptomes were primarily associated not only with bone but also vascular and connective tissue ontologies. A pattern-based approach used to screen the entire Gene Expression Omnibus (GEO) database against those that were sex-specific in bone identified two coordinately regulated gene sets: one related to high phosphate-induced aortic calcification and one induced by mechanical stimulation in bone. Temporal clustering of the transcriptome identified two skeletal tissue-associated, sex-specific patterns of gene expression. One set of genes, associated with skeletal patterning and morphology, showed peak expression earlier in females. The second set of genes, associated with coupled remodeling, had quantitatively higher expression in females and exhibited a broad peak between 3 to 12 months, concurrent with the animals' reproductive period. Results of phenome-level structural assessments of the tibia and vertebrae, and in vivo and in vitro analysis of cells having osteogenic potential, were consistent with the existence of functionally unique, skeletogenic cell populations that are separately responsible for appositional growth and intramedullary functions. These data suggest that skeletal sexual dimorphism arises through sex-specific, temporally different processes controlling morphometric growth and later coupled remodeling of the skeleton during the reproductive period of the animal. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.