TBX2-positive cells represent a multi-potent mesenchymal progenitor pool in the developing lung

Spatiotemporal transcriptome atlas of developing mouse lung

The functional development of the mammalian lung is a complex process that relies on the spatial and temporal organization of multiple cell types and their states. However, a comprehensive spatiotemporal transcriptome atlas of the developing lung has not yet been reported. Here we apply high-throughput spatial transcriptomics to allow for a comprehensive assessment of mouse lung development comprised of two critical developmental events: branching morphogenesis and alveologenesis. We firstly generate a spatial molecular atlas of mouse lung development spanning from E12.5 to P0 based on the integration of published single cell RNA-sequencing data and identify 10 spatial domains critical for functional lung organization. Furthermore, we create a lineage trajectory connecting spatial clusters from adjacent time points in E12.5-P0 lungs and explore TF (transcription factor) regulatory networks for each lineage specification. We observe the establishment of pulmonary airways within the developing lung, accompanied by the proximal-distal patterning with distinct characteristics of gene expression, signaling landscape and transcription factors enrichment. We characterize the alveolar niche heterogeneity with maturation state differences during the later developmental stage around birth and demonstrate differentially expressed genes, such as Angpt2 and Epha3, which may perform a critical role during alveologenesis. In addition, multiple signaling pathways, including ANGPT, VEGF and EPHA, exhibit increased levels in more maturing alveolar niche. Collectively, by integrating the spatial transcriptome with corresponding single-cell transcriptome data, we provide a comprehensive molecular atlas of mouse lung development with detailed molecular domain annotation and communication, which would pave the way for understanding human lung development and respiratory regeneration medicine.

Single-Cell RNA-Sequencing Identifies Bone Marrow-Derived Progenitor Cells as a Main Source of Extracellular Matrix-Producing Cells Across Multiple Organ-Based Fibrotic Diseases

Fibrosis is characterized by the aberrant deposition of extracellular matrix (ECM) due to dysregulated tissue repair responses, imposing a significant global burden on fibrosis-related diseases. Although alpha-smooth muscle actin (α-SMA/ACTA2)-expressing myofibroblasts are considered as key player in fibrogenesis, the origin of ECM-producing cells remains controversial. To address this issue, we integrated and analyzed large-scale single-cell transcriptomic datasets from patients with distinct fibrotic diseases involving the heart, lung, liver, or kidney. Unexpectedly, not all ACTA2-expressing cells were ECM-producing cells identified by expressing collagen genes; instead, the majority of ECM-producing cells were myofibroblasts and fibroblasts derived from circulating bone marrow precursor, and to a lesser extent from local pericytes and vascular smooth cells in all fibrotic diseases. This was confirmed in sex-mismatched kidney transplants by the discovery that ECM-producing cells originated from recipient, not donor, bone marrow-derived progenitor cells (BMPCs). Moreover, these BMPCs-derived ECM-producing cells exhibited a proinflammatory phenotype. Thus, bone marrow-derived proinflammatory and profibrotic fibroblasts/myofibroblasts with stem cell properties serve as a major source of ECM-producing cells and may play a driving role in tissue fibrosis across a wide range of human fibrotic diseases. Targeting these ECM-producing cells may provide a novel therapy for diseases with fibrosis.

A spatial sequencing atlas of age-induced changes in the lung during influenza infection

Influenza virus infection causes increased morbidity and mortality in the elderly. Aging impairs the immune response to influenza, both intrinsically and because of altered interactions with endothelial and pulmonary epithelial cells. To characterize these changes, we performed single-cell RNA sequencing (scRNA-seq), spatial transcriptomics, and bulk RNA sequencing (bulk RNA-seq) on lung tissue from young and aged female mice at days 0, 3, and 9 post-influenza infection. Our analyses identified dozens of key genes differentially expressed in kinetic, age-dependent, and cell type-specific manners. Aged immune cells exhibited altered inflammatory, memory, and chemotactic profiles. Aged endothelial cells demonstrated characteristics of reduced vascular wound healing and a prothrombotic state. Spatial transcriptomics identified novel profibrotic and antifibrotic markers expressed by epithelial and non-epithelial cells, highlighting the complex networks that promote fibrosis in aged lungs. Bulk RNA-seq generated a timeline of global transcriptional activity, showing increased expression of genes involved in inflammation and coagulation in aged lungs. Our work provides an atlas of high-throughput sequencing methodologies that can be used to investigate age-related changes in the response to influenza virus, identify novel cell-cell interactions for further study, and ultimately uncover potential therapeutic targets to improve health outcomes in the elderly following influenza infection.

Single-cell transcriptomic analysis uncovers the origin and intratumoral heterogeneity of parotid pleomorphic adenoma

Pleomorphic adenoma (PA) is the most common benign tumour in the salivary gland and has high morphological complexity. However, the origin and intratumoral heterogeneity of PA are largely unknown. Here, we constructed a comprehensive atlas of PA at single-cell resolution and showed that PA exhibited five tumour subpopulations, three recapitulating the epithelial states of the normal parotid gland, and two PA-specific epithelial cell (PASE) populations unique to tumours. Then, six subgroups of PASE cells were identified, which varied in epithelium, bone, immune, metabolism, stemness and cell cycle signatures. Moreover, we revealed that CD36+ myoepithelial cells were the tumour-initiating cells (TICs) in PA, and were dominated by the PI3K-AKT pathway. Targeting the PI3K-AKT pathway significantly inhibited CD36+ myoepithelial cell-derived tumour spheres and the growth of PA organoids. Our results provide new insights into the diversity and origin of PA, offering an important clinical implication for targeting the PI3K-AKT signalling pathway in PA treatment.

Maternal blood pressure associates with placental DNA methylation both directly and through alterations in cell-type composition

BACKGROUND

Maternal blood pressure levels reflect cardiovascular adaptation to pregnancy and proper maternal-fetal exchanges through the placenta and are very sensitive to numerous environmental stressors. Maternal hypertension during pregnancy has been associated with impaired placental functions and with an increased risk for children to suffer from cardiovascular and respiratory diseases later on. Investigating changes in placental DNA methylation levels and cell-type composition in association with maternal blood pressure could help elucidate its relationships with placental and fetal development.

METHODS

Taking advantage of a large cohort of 666 participants, we investigated the association between epigenome-wide DNA methylation patterns in the placenta, measured using the Infinium HumanMethylation450 BeadChip, placental cell-type composition, estimated in silico, and repeated measurements of maternal steady and pulsatile blood pressure indicators during pregnancy.

RESULTS

At the site-specific level, no significant association was found between maternal blood pressure and DNA methylation levels after correction for multiple testing (false discovery rate < 0.05), but 5 out of 24 previously found CpG associations were replicated (p-value < 0.05). At the regional level, our analyses highlighted 64 differentially methylated regions significantly associated with at least one blood pressure component, including 35 regions associated with mean arterial pressure levels during late pregnancy. These regions were found enriched for genes implicated in lung development and diseases. Further mediation analyses show that a significant part of the association between steady blood pressure-but not pulsatile pressure-and placental methylation can be explained by alterations in placental cell-type composition. In particular, elevated blood pressure levels are associated with a decrease in the ratio between mesenchymal stromal cells and syncytiotrophoblasts, even in the absence of preeclampsia.

CONCLUSIONS

This study provides the first evidence that the association between maternal steady blood pressure during pregnancy and placental DNA methylation is both direct and partly explained by changes in cell-type composition. These results could hint at molecular mechanisms linking maternal hypertension to lung development and early origins of childhood respiratory problems and at the importance of controlling maternal blood pressure during pregnancy.

Tbx2 is a master regulator of inner versus outer hair cell differentiation

The cochlea uses two types of mechanosensory cell to detect sounds. A single row of inner hair cells (IHCs) synapse onto neurons to transmit sensory information to the brain, and three rows of outer hair cells (OHCs) selectively amplify auditory inputs¹. So far, two transcription factors have been implicated in the specific differentiation of OHCs, whereas, to our knowledge, none has been identified in the differentiation of IHCs^2-4. One such transcription factor for OHCs, INSM1, acts during a crucial embryonic period to consolidate the OHC fate, preventing OHCs from transdifferentiating into IHCs². In the absence of INSM1, embryonic OHCs misexpress a core set of IHC-specific genes, which we predict are involved in IHC differentiation. Here we find that one of these genes, Tbx2, is a master regulator of IHC versus OHC differentiation in mice. Ablation of Tbx2 in embryonic IHCs results in their development as OHCs, expressing early OHC markers such as Insm1 and eventually becoming completely mature OHCs in the position of IHCs. Furthermore, Tbx2 is epistatic to Insm1: in the absence of both genes, cochleae generate only OHCs, which suggests that TBX2 is necessary for the abnormal transdifferentiation of INSM1-deficient OHCs into IHCs, as well as for normal IHC differentiation. Ablation of Tbx2 in postnatal, largely differentiated IHCs makes them transdifferentiate directly into OHCs, replacing IHC features with those of mature and not embryonic OHCs. Finally, ectopic expression of Tbx2 in OHCs results in their transdifferentiation into IHCs. Hence, Tbx2 is both necessary and sufficient to make IHCs distinct from OHCs and maintain this difference throughout development.

A Competitive Endogenous RNA Network Based on Differentially Expressed lncRNA in Lipopolysaccharide-Induced Acute Lung Injury in Mice

Non-coding RNAs have remarkable roles in acute lung injury (ALI) initiation. Nevertheless, the significance of long non-coding RNAs (lncRNAs) in ALI is still unknown. Herein, we purposed to identify potential key genes in ALI and create a competitive endogenous RNA (ceRNA) modulatory network to uncover possible molecular mechanisms that affect lung injury. We generated a lipopolysaccharide-triggered ALI mouse model, whose lung tissue was subjected to RNA sequencing, and then we conducted bioinformatics analysis to select genes showing differential expression (DE) and to build a lncRNA-miRNA (microRNA)- mRNA (messenger RNA) modulatory network. Besides, GO along with KEGG assessments were conducted to identify major biological processes and pathways, respectively, involved in ALI. Then, RT-qPCR assay was employed to verify levels of major RNAs. A protein-protein interaction (PPI) network was created using the Search Tool for the Retrieval of Interacting Genes (STRING) database, and the hub genes were obtained with the Molecular Complex Detection plugin. Finally, a key ceRNA subnetwork was built from major genes and their docking sites. Overall, a total of 8,610 lncRNAs were identified in the normal and LPS groups. Based on the 308 DE lncRNAs [p-value < 0.05, |log2 (fold change) | > 1] and 3,357 DE mRNAs [p-value < 0.05, |log2 (fold change) | > 1], lncRNA-miRNA and miRNA-mRNA pairs were predicted using miRanda. The lncRNA-miRNA-mRNA network was created from 175 lncRNAs, 22 miRNAs, and 209 mRNAs in ALI. The RT-qPCR data keep in step with the RNA sequencing data. GO along with KEGG analyses illustrated that DE mRNAs in this network were mainly bound up with the inflammatory response, developmental process, cell differentiation, cell proliferation, apoptosis, and the NF-kappa B, PI3K-Akt, HIF-1, MAPK, Jak-STAT, and Notch signaling pathways. A PPI network on the basis of the 209 genes was established, and three hub genes (Nkx2-1, Tbx2, and Atf5) were obtained from the network. Additionally, a lncRNA-miRNA-hub gene subnetwork was built from 15 lncRNAs, 3 miRNAs, and 3 mRNAs. Herein, novel ideas are presented to expand our knowledge on the regulation mechanisms of lncRNA-related ceRNAs in the pathogenesis of ALI.

Single-Cell Epigenomics and Functional Fine-Mapping of Atherosclerosis GWAS Loci

Supplemental Digital Content is available in the text.

Genome-wide association studies have identified hundreds of loci associated with coronary artery disease (CAD). Many of these loci are enriched in cisregulatory elements but not linked to cardiometabolic risk factors nor to candidate causal genes, complicating their functional interpretation.

Combined genomic and proteomic approaches reveal DNA binding sites and interaction partners of TBX2 in the developing lung

BACKGROUND

Tbx2 encodes a transcriptional repressor implicated in the development of numerous organs in mouse. During lung development TBX2 maintains the proliferation of mesenchymal progenitors, and hence, epithelial proliferation and branching morphogenesis. The pro-proliferative function was traced to direct repression of the cell-cycle inhibitor genes Cdkn1a and Cdkn1b, as well as of genes encoding WNT antagonists, Frzb and Shisa3, to increase pro-proliferative WNT signaling. Despite these important molecular insights, we still lack knowledge of the DNA occupancy of TBX2 in the genome, and of the protein interaction partners involved in transcriptional repression of target genes.

METHODS

We used chromatin immunoprecipitation (ChIP)-sequencing and expression analyses to identify genomic DNA-binding sites and transcription units directly regulated by TBX2 in the developing lung. Moreover, we purified TBX2 containing protein complexes from embryonic lung tissue and identified potential interaction partners by subsequent liquid chromatography/mass spectrometry. The interaction with candidate proteins was validated by immunofluorescence, proximity ligation and individual co-immunoprecipitation analyses.

RESULTS

We identified Il33 and Ccn4 as additional direct target genes of TBX2 in the pulmonary mesenchyme. Analyzing TBX2 occupancy data unveiled the enrichment of five consensus sequences, three of which match T-box binding elements. The remaining two correspond to a high mobility group (HMG)-box and a homeobox consensus sequence motif. We found and validated binding of TBX2 to the HMG-box transcription factor HMGB2 and the homeobox transcription factor PBX1, to the heterochromatin protein CBX3, and to various members of the nucleosome remodeling and deacetylase (NuRD) chromatin remodeling complex including HDAC1, HDAC2 and CHD4.

CONCLUSION

Our data suggest that TBX2 interacts with homeobox and HMG-box transcription factors as well as with the NuRD chromatin remodeling complex to repress transcription of anti-proliferative genes in the pulmonary mesenchyme.

Transcriptome analysis reveals a protective role of liver X receptor alpha against silica particle-induced experimental silicosis

Silicosis, a severe and irreversible form of pulmonary fibrosis (PF) caused by long-term exposure to dust particles in production environments, is the biggest occupational health concern in China and most low-income countries. The transdifferentiation of pulmonary fibroblasts is the terminal event in silicosis, and specific transcription factors (TFs) play a crucial role in this condition. However, the relationship between TF-mediated regulation and silicosis remains unknown. We performed a transcriptomic analysis to elucidate this relationship, and our results revealed that two TFs, EGR2 and BHLHE40, were upregulated and five, i.e., TBX2, NR1H3 (LXRα), NR2F1, PPARG (PPARγ), and EPAS1, were downregulated in activated fibroblasts. Notably, PPARγ and LXRα expression was also decreased in an experimental mouse model of silicosis. The mechanism underlying these changes may involve TGF-β1 secretion from silica-exposed alveolar macrophages, causing PPARγ and LXRα downregulation, which in turn would result in aberrant α-SMA transcription. Our results suggest that LXRα is a potential target for the prevention of silicosis and PF.