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, log₂ 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.

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.