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.

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.

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.

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.

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.

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.

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.

Cross-Site Evaluation of Commercial Sanger Sequencing Chemistries.

Sanger sequencing remains an essential tool utilized by researchers. Despite competition from commercial sequencing providers, many academic sequencing core facilities continue to offer these services based on a model of competitive pricing, knowledgeable technical support, and rapid turnaround time. In-house Sanger sequencing remains a viable core service and, until recently, Applied Biosystems BigDye Terminator chemistry was the only commercially available solution for Sanger DNA sequencing on Applied Biosystems (ABI) instruments; however, several new products employing novel dye chemistries and reaction configurations have entered the market. As a result, there is a need to benchmark the performance of these new chemistries on various DNA templates, including difficult-to-sequence templates, and their amenability to commonly employed cost-saving measures, such as dye dilution and reaction miniaturization. To evaluate these new reagents, a study was designed to compare the quality of Sanger sequencing data produced by ABI BigDye and commercially available kits from 2 other vendors using both control and difficult-to-sequence DNA templates under various reaction conditions. This study will serve as a valuable resource to core facilities conducting Sanger sequencing that wish to evaluate the use of an alternative chemistry in their sequencing core.

Multisite Evaluation of Next-Generation Methods for Small RNA Quantification.

Small RNAs (smRNAs) are important regulators of many biologic processes and are now most frequently characterized using Illumina sequencing. However, although standard RNA sequencing library preparation has become routine in most sequencing facilities, smRNA sequencing library preparation has historically been challenging because of high input requirements, laborious protocols involving gel purifications, inability to automate, and a lack of benchmarking standards. Additionally, studies have suggested that many of these methods are nonlinear and do not accurately reflect the amounts of smRNAs in vivo. Recently, a number of new kits have become available that permit lower input amounts and less laborious, gel-free protocol options. Several of these new kits claim to reduce RNA ligase-dependent sequence bias through novel adapter modifications and to lessen adapter-dimer contamination in the resulting libraries. With the increasing number of smRNA kits available, understanding the relative strengths of each method is crucial for appropriate experimental design. In this study, we systematically compared 9 commercially available smRNA library preparation kits as well as NanoString probe hybridization across multiple study sites. Although several of the new methodologies do reduce the amount of artificially over- and underrepresented microRNAs (miRNAs), we observed that none of the methods was able to remove all of the bias in the library preparation. Identical samples prepared with different methods show highly varied levels of different miRNAs. Even so, many methods excelled in ease of use, lower input requirement, fraction of usable reads, and reproducibility across sites. These differences may help users select the most appropriate methods for their specific question of interest.

Tobacco-Related Alterations in Airway Gene Expression are Rapidly Reversed Within Weeks Following Smoking-Cessation.

The physiologic response to tobacco smoke can be measured by gene-expression profiling of the airway epithelium. Temporal resolution of kinetics of gene-expression alterations upon smoking-cessation might delineate distinct biological processes that are activated during recovery from tobacco smoke exposure. Using whole genome gene-expression profiling of individuals initiating a smoking-cessation attempt, we sought to characterize the kinetics of gene-expression alterations in response to short-term smoking-cessation in the nasal epithelium. RNA was extracted from the nasal epithelial of active smokers at baseline and at 4, 8, 16, and 24-weeks after smoking-cessation and put onto Gene ST arrays. Gene-expression levels of 119 genes were associated with smoking-cessation (FDR < 0.05, FC ≥1.7) with a majority of the changes occurring by 8-weeks and a subset changing by 4-weeks. Genes down-regulated by 4- and 8-weeks post-smoking-cessation were involved in xenobiotic metabolism and anti-apoptotic functions respectively. These genes were enriched among genes previously found to be induced in smokers and following short-term in vitro exposure of airway epithelial cells to cigarette smoke (FDR < 0.05). Our findings suggest that the nasal epithelium can serve as a minimally-invasive tool to measure the reversible impact of smoking and broadly, may serve to assess the physiological impact of changes in smoking behavior.

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.

A Next-Generation Sequencing Primer-How Does It Work and What Can It Do?

Next-generation sequencing refers to a high-throughput technology that determines the nucleic acid sequences and identifies variants in a sample. The technology has been introduced into clinical laboratory testing and produces test results for precision medicine. Since next-generation sequencing is relatively new, graduate students, medical students, pathology residents, and other physicians may benefit from a primer to provide a foundation about basic next-generation sequencing methods and applications, as well as specific examples where it has had diagnostic and prognostic utility. Next-generation sequencing technology grew out of advances in multiple fields to produce a sophisticated laboratory test with tremendous potential. Next-generation sequencing may be used in the clinical setting to look for specific genetic alterations in patients with cancer, diagnose inherited conditions such as cystic fibrosis, and detect and profile microbial organisms. This primer will review DNA sequencing technology, the commercialization of next-generation sequencing, and clinical uses of next-generation sequencing. Specific applications where next-generation sequencing has demonstrated utility in oncology are provided.

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.

TLR3 Signaling Promotes the Induction of Unique Human BDCA-3 Dendritic Cell Populations.

Conventional and plasmacytoid dendritic cells (cDCs and pDCs) are the two populations of DCs that can be readily identified in human blood. Conventional DCs have been subdivided into CD1c(+), or blood dendritic cells antigen (BDCA) 1 and CD141(+), or BDCA-3, DCs, each having both unique gene expression profiles and functions. BDCA-3 DCs express high levels of toll-like receptor 3 and upon stimulation with Poly I:C secrete IFN-ß, CXCL10, and IL-12p70. In this article, we show that activation of human BDCA-3 DCs with Poly I:C induces the expression of activation markers (CD40, CD80, and CD86) and immunoglobulin-like transcript (ILT) 3 and 4. This Poly I:C stimulation results in four populations identifiable by flow cytometry based on their expression of ILT3 and ILT4. We focused our efforts on profiling the ILT4(-) and ILT4(+) DCs. These ILT-expressing BDCA-3 populations exhibit similar levels of activation as measured by CD40, CD80, and CD86; however, they exhibit differential cytokine secretion profiles, unique gene signatures, and vary in their ability to prime allogenic naïve T cells. Taken together, these data illustrate that within a pool of BDCA-3 DCs, there are cells poised to respond differently to a given input stimulus with unique output of immune functions.

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.

Defining the nasal transcriptome in granulomatosis with polyangiitis (Wegener's).

OBJECTIVE: To determine whether disease processes related to granulomatosis with polyangiitis (Wegener's) (GPA) are reflected in gene expression profiles of the nasal mucosa. METHODS: Nasal brushings of the inferior turbinate were obtained from 32 patients with GPA (10 with active nasal disease, 13 with prior nasal disease, and 9 with no history of nasal disease) and a composite comparator group with and without inflammatory nasal disease (12 healthy people, 15 with sarcoidosis, and 8 with allergic rhinitis). Differential gene expression was assessed between subgroups of GPA and comparators. RESULTS: A total of 339 genes were differentially expressed between the GPA and comparator groups (absolute fold change >1.5; false discovery rate <0.05). Top canonical pathways up-regulated in nasal brushings from patients with GPA included granulocyte adhesion and diapedesis (P = 8.6(-22) ), agranulocyte adhesion and diapedesis (P = 1.3(-14) ), IL10 signaling (P = 3.0(-11) ), LXR/RXR activation (P = 4.3(-11) ), and TREM1 signaling (P = 9.0(-11) ). A set of genes differentially expressed in GPA independently of nasal disease activity status included genes related to epithelial barrier integrity (fibronectin 1, desmosomal proteins) and several matricellular proteins (e.g., osteonectin, osteopontin). Significant overlap of differentially expressed genes was observed between active and prior nasal disease GPA subgroups. Peripheral blood neutrophil and mononuclear gene expression levels associated with GPA were similarly altered in the nasal gene expression profiles of patients with active or prior nasal disease. CONCLUSION: Profiling the nasal transcriptome in GPA reveals gene expression signatures related to innate immunity, inflammatory cell chemotaxis, extracellular matrix composition, and epithelial barrier integrity. Thus, airway-based expression profiling is feasible and informative in GPA.

Evaluation of commercially available RNA amplification kits for RNA sequencing using very low input amounts of total RNA.

This article includes supplemental data. Please visit http://www.fasebj.org to obtain this information.Multiple recent publications on RNA sequencing (RNA-seq) have demonstrated the power of next-generation sequencing technologies in whole-transcriptome analysis. Vendor-specific protocols used for RNA library construction often require at least 100 ng total RNA. However, under certain conditions, much less RNA is available for library construction. In these cases, effective transcriptome profiling requires amplification of subnanogram amounts of RNA. Several commercial RNA amplification kits are available for amplification prior to library construction for next-generation sequencing, but these kits have not been comprehensively field evaluated for accuracy and performance of RNA-seq for picogram amounts of RNA. To address this, 4 types of amplification kits were tested with 3 different concentrations, from 5 ng to 50 pg, of a commercially available RNA. Kits were tested at multiple sites to assess reproducibility and ease of use. The human total reference RNA used was spiked with a control pool of RNA molecules in order to further evaluate quantitative recovery of input material. Additional control data sets were generated from libraries constructed following polyA selection or ribosomal depletion using established kits and protocols. cDNA was collected from the different sites, and libraries were synthesized at a single site using established protocols. Sequencing runs were carried out on the Illumina platform. Numerous metrics were compared among the kits and dilutions used. Overall, no single kit appeared to meet all the challenges of small input material. However, it is encouraging that excellent data can be recovered with even the 50 pg input total RNA.

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.

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.