- Open Access
Perturbation of semaphorin and VEGF signaling in ACDMPV lungs due to FOXF1 deficiency
Respiratory Research volume 22, Article number: 212 (2021)
Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a rare lethal congenital lung disorder in neonates characterized by severe progressive respiratory failure and refractory pulmonary hypertension, resulting from underdevelopment of the peripheral pulmonary tree. Causative heterozygous single nucleotide variants (SNVs) or copy-number variant (CNV) deletions involving FOXF1 or its distant lung-specific enhancer on chromosome 16q24.1 have been identified in 80–90% of ACDMPV patients. FOXF1 maps closely to and regulates the oppositely oriented FENDRR, with which it also shares regulatory elements.
To better understand the transcriptional networks downstream of FOXF1 that are relevant for lung organogenesis, using RNA-seq, we have examined lung transcriptomes in 12 histopathologically verified ACDMPV patients with or without pathogenic variants in the FOXF1 locus and analyzed gene expression profile in FENDRR-depleted fetal lung fibroblasts, IMR-90.
RNA-seq analyses in ACDMPV neonates revealed changes in the expression of several genes, including semaphorins (SEMAs), neuropilin 1 (NRP1), and plexins (PLXNs), essential for both epithelial branching and vascular patterning. In addition, we have found deregulation of the vascular endothelial growth factor (VEGF) signaling that also controls pulmonary vasculogenesis and a lung-specific endothelial gene TMEM100 known to be essential in vascular morphogenesis. Interestingly, we have observed a substantial difference in gene expression profiles between the ACDMPV samples with different types of FOXF1 defect. Moreover, partial overlap between transcriptome profiles of ACDMPV lungs with FOXF1 SNVs and FENDRR-depleted IMR-90 cells suggests contribution of FENDRR to ACDMPV etiology.
Our transcriptomic data imply potential crosstalk between several lung developmental pathways, including interactions between FOXF1-SHH and SEMA-NRP or VEGF/VEGFR2 signaling, and provide further insight into complexity of lung organogenesis in humans.
Alveolar Capillary Dysplasia with Misalignment of Pulmonary Veins (ACDMPV, MIM #265380) is a rarely diagnosed lethal lung developmental disorder (LLDD) in neonates [1,2,3,4,5,6]. Affected newborns manifest severe progressive hypoxemic respiratory failure and pulmonary arterial hypertension (PAH), refractory to treatment. Most ACDMPV infants show disease symptoms within the first 24–48 h of life and die usually within the first month of life. Very rarely, the course of ACDMPV is milder and patients survive beyond the neonatal period [7,8,9,10,11,12,13,14,15]. Histopathological findings in ACDMPV include reduced number of pulmonary capillaries, majority of which do not make contact with the alveolar epithelium, thickening of intra-alveolar septa and reduced alveolarization, medial hypertrophy of small peripheral pulmonary arteries and arterioles, and the presence of vascular anastomoses. In addition, most patients with ACDMPV often manifest anomalies of the cardiovascular, gastrointestinal, and genitourinary systems [3,4,5, 16].
In approximately 80–90% of ACDMPV patients, heterozygous single nucleotide variants (SNVs) or copy-number variant (CNV) deletions involving FOXF1 (MIM #601089) and/or its non-coding regulatory regions on chromosome 16q24.1 have been found [16,17,18]. FOXF1 encodes a forkhead-box (FOX) family transcription factor [19, 20]. In developing lungs, FOXF1 is expressed in peripheral lung mesenchyme where it mediates sonic hedgehog (SHH) signaling from epithelial cells of branching tubular structures [21,22,23,24,25,26,27,28,29]. FOXF1 expression in the lungs is controlled by a distant lung-specific enhancer region located ~ 270 kb upstream to the gene (chr16:86,178,434–86,238,313; hg38) [16, 18, 30,31,32].
FOXF1 maps 1.7 kb from the oppositely oriented long noncoding RNA (lncRNA) gene, FENDRR, with which it shares regulatory regions that include the promoters and the enhancer. Most recently, FENDRR expression was found to be regulated both in cis by the FOXF1 distant lung-specific enhancer and in trans by FOXF1, implying involvement of FENDRR in FOXF1-linked diseases including ACDMPV .
Here, we have examined lung transcriptomes of ACDMPV patients with or without pathogenic variants in the FOXF1 locus on chromosome 16q24.1 to identify gene expression profiles specific for ACDMPV. Using comparative analyses and an over-representation approach, we have investigated deregulated pathways to elucidate potential crosstalk between FOXF1 and other genes required for normal lung development in humans. We have compared also transcriptome profiles of ACDMPV lungs with transcriptomes of FENDRR-depleted fetal lung fibroblasts, IMR-90, to further consider the possibility of the contribution of FENDRR deficiency to etiology of ACDMPV.
Patients and samples
Peripheral blood, saliva, and/or frozen or formalin-fixed paraffin-embedded (FFPE) lung biopsy or autopsy samples were obtained from 12 unrelated patients with histopathologically verified ACDMPV (Additional file 1, Additional file 2). Three normal lung tissues from age-matched individuals were used as controls. The study protocol was approved by the Institutional Review Board for Human Subject Research at Baylor College of Medicine (BCM; H-8712).
Human fetal lung fibroblasts IMR-90 (ATCC, Manassas, VA) were cultured at 37 °C in 5% CO2 in Eagle's Minimal Essential Medium (ATCC), supplemented with 2 mM l-glutamine, 1 mM sodium pyruvate, and 10% fetal bovine serum.
DNA and RNA extraction
DNA was isolated from blood and saliva using Gentra Purgene Blood Kit (Qiagen, Germantown, MD), or from frozen lung tissue using DNeasy Blood and Tissue Kit (Qiagen). RNA from frozen lung or FFPE lung tissues was extracted using miRNeasy Mini Kit (Qiagen) or Quick-RNA FFPE Kit (Zymo Research, Irvine, CA), respectively. DNA and RNA from cultured IMR-90 cells was isolated using DNaesy Blood and Tissue Kit (Qiagen) and miRNeasy Mini Kit (Qiagen), respectively.
SNVs were identified by Sanger sequencing of PCR fragments amplified from genomic DNA. Pathogenic CNV deletions were identified and mapped by array comparative genomic hybridization (array CGH), using high resolution custom-designed 4 × 180 K oligo arrays (Agilent Technologies, Santa Clara, CA), followed by Sanger sequencing of long-range-PCR-amplified deletion junctions. Parental origin of the deletion-bearing chromosome 16q24.1 was determined using informative microsatellites or single nucleotide polymorphisms (SNPs).
Based on genomic findings, the ACDMPV patients were divided into five groups, four of which included cases featuring FOXF1-related defects at chromosome 16q24.1: CNV deletions involving FOXF1 enhancer (n = 4; group 1), FOXF1 and its enhancer (n = 3; group 2) or FOXF1 (n = 2; group 3), as well as FOXF1 mutation (n = 2; group 4). The fifth group contains a single ACDMPV case with no known FOXF1 abnormalities. In particular, heterozygous CNV deletion on 16q24.1, involving FOXF1 and its enhancer, was identified in patient 3.2. The deletion arose de novo on maternal chromosome 16. The genomic findings in the remaining ACDMPV neonates were described previously [16,17,18, 34,35,36]. Detailed sample description is presented in Additional file 1.
FENDRR knockout in lung fibroblasts
FENDRR knockout in IMR-90 cells was performed using CRISPR/Cas9 system by deletion of a 0.6 kb region of chromosome 16q24.1 containing two thirds of the putative FENDRR promoter (candidate cis-Regulatory Element, cCRE, in ENCODE; http://www.encodeproject.org), the entire exon 1, and adjacent to it, a portion of the intron 1 (Additional file 3). Upstream and downstream single guide RNA (sgRNA) sequences (5’ TGTGAGTTCAGATCCGAGCG; chr16: 86,509,150–86,509,169, hg38 and 5’ GGTCATCGCGAAGAGCCTAG; chr16: 86,508,542–86,508,561, hg38, respectively), flanking the region to be deleted, were designed using the CRISPOR Design Tool (http://crispor.tefor.net/; ). Double stranded oligonucleotides (Integrated DNA Technologies, Coralville, IA), encoding guide RNAs, were cloned in the BbsI sites of the pSpCas9(BB)-2A-GFP (pX458) plasmids (Addgene, Watertown, MA), which constitutively express Cas9 and GFP. A culture of 2 × 106 IMR-90 lung fibroblasts was co-transfected with two pX458-derived vectors, expressing either upstream or downstream sgRNA, using Lipofectamine 3000 (Invitrogen, Carlsbad, CA). GFP-positive cells were collected 48 h later applying Fluorescence Assisted Cell Sorting (FACS) at the BCM Cytometry and Cell Sorting Core (on FACSAriaII sorter with nozzle size of 70 µm). The deletion was confirmed by PCR amplification and Sanger sequencing of the genomic region across the deletion junction. GFP-positive cells were subsequently processed to isolate total RNA for RNA-seq.
RNA-seq based transcriptome profiling
The RNA-seq-based transcriptomic profiling was performed in 12 ACDMPV-affected lung samples and three age-matched normal lung controls (Additional file 1). Transcriptomic analyses were also performed in IMR-90 fetal lung fibroblasts with FENDRR knockout obtained using CRISPR/Cas9 system, and in IMR-90 cells in which FENDRR expression was reduced by applying RNA interference (RNAi) targeting spliced and un-spliced FENDRR transcripts . RNA-seq experiments were performed using the Illumina NovaSeq 6000 System at BCM Genomic and RNA Profiling Core.
RNA-seq raw reads were mapped to the human genome (GRCh38) using STAR aligner (v 2.7.0). To generate the counts of reads uniquely mapped to the known and annotated genes, the featureCount method (R package Rsubread) and the Ensembl annotation file GRCh38.95.gtf were used. The count table of the uniquely mapped reads was generated and differential expression was tested using the DESeq2 (R Bioconductor package). Heatmaps were generated using the ClustVis online tool at https://biit.cs.ut.ee/clustvis/ . To perform the Principal Component Analysis (PCA), the count data were transformed using the VarianceStabilizingTransformation method and visualized using the plotPCA method (both implemented in DESeq2 package).
Pathway over-representation analysis across differentially expressed genes
The over-representation analysis of molecular pathways in each group was performed using the ConsensusPathDB server  accessible at http://ConsensusPathDB.org. Collection of annotated gene sets for further pathway assessment were obtained through the Molecular Signatures Database (MSigDB, v7.2, http://www.gsea-msigdb.org/gsea/msigdb/index.jsp) [40, 41]. Functional annotation of deregulated genes in Groups 1–5 was determined by Gene Ontology (GO) enrichment analyses using the clusterProfiler .
Validation of RNA-Seq data by targeted gene expression analysis using NanoString
Verification of RNA-seq data was done using nCounter platform (NanoString Technologies, Seattle, WA). A hundred and twenty selected genes, whose expression was found by RNA-seq to be mis-regulated in ACDMPV patients, were included in the custom code sets. GAPDH, ACTB, and PGK1 were used as the reference for gene expression normalization. To compare NanoString and RNA-seq data, the average signal obtained in NanoString experiment was separately calculated for three patients with ACDMPV (1.1, 4.1, and 4.2) and three controls (C1, C3, and C4). All RNA samples used for NanoString experiment have been assessed also in the prior RNA-Seq studies. The ratio between the average signal from NanoString has been compared to the ratio of the average normalized RNA-seq read count in patients and controls.
Comparison of RNA-seq and microarray data
To compare the results of this RNA-seq and of our previous expression microarrays analysis , we have selected the genes that were consistently either up-regulated or down-regulated in both datasets. We have compared protein coding genes deregulated in each ACDMPV group (adjusted p-value < 0.01; log2 fold change > 1.0 or < − 1.0) with the genes deregulated in ACDMPV patients in the microarray study (adjusted p-value < 0.01). For the RNA-seq data, we have required that both direction of deregulation are consistent among all five ACDMPV groups.
Perturbations of ACDMPV transcriptomes
The number of reads that were uniquely mapped to the reference human genome varied between 5.6 and 40.1 million per sample in Group 1, 6.5 and 13.8 million in Group 2, 4.8 and 16.8 million in Group 3, 30.9 and 36.2 million in Group 4, and were 120.5 million in Group 5. Statistical analyses of RNA-seq data revealed on average ~ 2,600 differentially expressed genes in each ACDMPV group compared to normal control samples (adjusted p-value < 0.01; log2 fold change > 1.0 or < − 1.0). The highest number of deregulated protein-coding genes was observed in group 2 (n = 7525), while the lowest number was detected in group 1 (n = 682) (Additional file 4). The visual summary of the differentially expressed genes in each group is presented in Fig. 1A.
PCA showed that all analyzed ACDMPV cases clustered away from the normal lung controls (Fig. 2A). While group 2 (CNV deletion of FOXF1 and its enhancer) and group 4 (FOXF1 mutations) form separate clusters, group 1 (CNV deletion involving FOXF1 enhancer) and group 3 (CNV deletion involving FOXF1) are scattered (Fig. 2A).
Ninety-four genes were commonly deregulated in all ACDMPV individuals (Additional file 5). Notably, a changed expression of genes involved in vascular patterning or epithelial branching, TMEM100, SEMA3B, and EPHB6, was found in each patient. Deregulation of several genes previously described as involved in lung extracellular matrix formation (COL16A1 and MEGF6), lung cancer (IL7R and YES1), pulmonary fibrosis (MMP23B, PTP4A1, and MUC4), or chronic obstructive pulmonary disease (LRRC45 and MTCL1) was also observed in each group (Fig. 2B).
Pathway over-representation analysis across differentially expressed genes in ACDMPV patients
Over-representation analyses performed on the combined lists of upregulated and downregulated transcripts revealed significant enrichment of differentially expressed protein-coding genes from several molecular pathways linked to lung development and function (Additional file 6).
In group 1, a significant over-representation includes the genes encoding axon guidance proteins—ephrins (EPHs), semaphorins (SEMAs), and plexins (PLXNs), hedgehog signaling pathway—SHH and KIF7, and proteins involved in O2/CO2 transport—HBA1, CA1, and HBB, whereas in group 2 those involved in TGF-beta and WNT signaling (TGFB1, TGFBR1-3, BMPR2, CTNNB1, ACVR2A, ACVR2B, SMAD, and MAPK genes, FGF6 and ANGPTL1), PDGF and vascular endothelial growth factor (VEGF) signaling (PDGFC, VEGFC, and VEGFD), SEMA-neuropilin (SEMA-NRP) signaling (SEMA3B and NRP1) as well as genes involved in NOTCH (NOTCH2) and SHH (SHH) signaling were detected. In group 3, we have found deregulated genes that encode proteins involved in NOTCH3 signaling (APH1A and SNW1), TGF-beta, and SMAD, and MAPK signaling (TGFBR1, SP1, and MAPK1) and in group 4 those, involved in axon guidance and SEMA and PLXN receptors, and IL-17 signaling pathway (MAPK and IL17 receptors). Among pathways affected in group 5, there were axon guidance and SEMA signaling (PLXN genes), SHH and NOTCH signaling (GLI2, KIF7, and NOTCH3), ephrin signaling (EPHA2, EPHA1, EPHB3, and YES1), VEGF signaling (MMP14, HDAC5, and JUN), and MAPK signaling (MAP2K2 and TLN1) (Additional file 6).
Differentially expressed genes were subjected to the GO enrichment analyses. We have found 20 enriched GO terms for Group 1, 619 for Group 2, 254 for Group 3, 82 for Group 4, and 282 for Group 5. The top overrepresented GO biological processes (BP) include regulation of GTPase activity (Groups 1 and 4), mRNA catabolic process (Group 2), translational initiation (Group 3), and extracellular structure organization (Group 5) (Additional file 7). Among other significant GO BP terms was semaphorin–plexin signaling pathway involved in axon guidance (Group 4). Within GO molecular function (MF), the most enriched ones were associated with small GTPase binding (Groups 1 and 4), cadherin binding (Group 2), structural constituent of ribosome (Group 3), and actin binding (Group5) (Additional file 7).
Hierarchical clustering of genes from commonly deregulated pathways: axon guidance (R-HSA-422475), SEMA interactions (R-HSA-373755), and VEGF signaling (R-HSA-194138) showed separation of ACDMPV and non-ACDMPV samples (Additional file 8, Additional file 9, Additional file 10).
Deregulation of genes in FENDRR-depleted lung fibroblasts
Analysis of the RNA-seq data for IMR-90 cells with FENDRR knockdown and knockout  revealed deregulation of 886 protein-coding genes with log2 fold change > 1.0 or < − 1.0. These genes were significantly over-represented in several molecular pathways, including signaling by NOTCH2 (HES5 and FCER2) or JAK-STAT signaling (CSF2RA, IL10, IL13, and IL23R) (Additional file 11). Comparison of gene expression patterns in FENDRR-depleted IMR-90 cells (knockout and knockdown) and lung tissues (ACDMPV patients with FENDRR deletion; Groups 2 and 3) revealed 65 commonly deregulated genes (Fig. 1B).
Validation of the RNA-Seq results by NanoString assay
The 120 values for analyzed genes obtained by both methods showed a strong correlation between the two sets of results (Additional file 12).
Comparison of RNA-seq data with microarray-based expression results
We have found a partial overlap between the RNA-seq results and our previous expression microarray data . In Group 1 we have identified 46 genes that were commonly deregulated in both datasets. In Groups 2, 3, 4, and 5, we have detected 379, 77, 76, and 144 differentially expressed genes, respectively, that were also found deregulated in the microarray study. (Additional file 13).
In this study, we have performed a comprehensive RNA-seq profiling in the lung samples from 12 ACDMPV patients, classified into five separate categories based on the type of the FOXF1 defect. While we have identified the differences between the transcriptomic profiles in ACDMPV lungs when compared to the control tissues, a substantial difference in gene expression levels among ACDMPV study groups was also observed. Nonetheless, 94 genes were significantly deregulated in all ACDMPV individuals, suggesting their contribution to the disease. Those genes can be used as a reference set in the diagnosis of ACDMPV cases in which neither causative SNV nor CNV deletion involving the FOXF1 locus is detected.
In all ACDMPV lung specimens, we have found an increased expression of SEMA3B. SEMA3B belongs to the semaphorin family which, besides being involved in axon guidance signaling, was also found to be responsible for murine lung bud formation and normal alveolar development . While the role of SEMA3B in lung growth is not well characterized, SEMA3A was found to inhibit branching morphogenesis of fetal mouse lung  whereas SEMA3C and SEMA3F were reported to promote airway branching and increase proliferation of terminal epithelial cells [45,46,47]. Of note, in this study, SEMA3C and other SEMA genes were also significantly deregulated in ACDMPV patients with deletion of FOXF1 and its distant lung-specific enhancer.
Semaphorins act mostly through binding with PLXNs and NRPs, their primary receptors and co-receptors, respectively [48, 49]. Notably, in all but one studied ACDMPV groups, PLXNs, including PLXNA3, PLXNB1-3, and PLXND1, were found to be significantly deregulated. Also, the expression level of the NRP1 receptor, recently shown to facilitate SARS-CoV-2 cell entry , was decreased in groups 2 and 3. In support of this notion, disruption of SEMA-NRP1 signaling in murine lungs was reported to cause instability of the alveolar-capillary interface and hypertensive remodeling reminiscent of ACD [51, 52]. Our data on SEMA-NRP1 deregulation further imply the role of this signaling pathway in lung development.
Deregulation of SEMAs, PLXNs, and NRP1 in ACDMPV patients with various FOXF1 abnormalities is supported by our previous ChIP-seq analyses performed in E18.5 lungs of mice overexpressing Foxf1 that revealed direct interactions of FOXF1 with their loci . In addition, in two unrelated ACDMPV neonates without any detectable variant in FOXF1, we previously reported the c.631C > G (p.Leu211Val) and c.3256G > A (p.Ala1086Thr) variants in PLXNB2, inherited from their healthy fathers . Overall, these data imply potential interaction between the SHH-FOXF1 and SEMA-NRP signaling pathways.
In addition to SEMAs, also VEGFs are ligands to PLXNs, and NRP1. VEGF controls pulmonary vasculogenesis and its inhibition in rats resulted in alveolar simplification and loss of lung capillaries . Decreased expression of the VEGF receptor gene, Flt1, was found in the Foxf1 conditional knockout mouse endothelium, implying that FOXF1 can regulate endothelial genes involved in VEGF signaling . Here, downregulation of FLT1 was observed in ACDMPV patients with the FOXF1 enhancer deletions, further suggesting that this gene can be regulated by the FOXF1 pathway. Of note, VEGF signaling is connected to the axon guidance pathway through NRPs, which can also bind members of the VEGF family .
EPHB6, a member of the ephrin family that is upregulated in all ACDMPV patients, similarly to SEMAs, is also involved in axon guidance signaling. A possible role of EPHs as molecular stimulators of lung morphogenesis was reported in rat lungs . In humans, EPHB6 is implicated in inhibition of metastasis in several types of cancer, including non-small cell lung cancer .
In all ACDMPV newborns, we have found a decreased level of lung-specific TMEM100 encoding Transmembrane protein 100. In humans and mice, this gene was found to be most abundantly expressed in the lungs (Additional file 14) . It was demonstrated that TMEM100 acts downstream of BMP9/BMP10-ALK1 (ACVRL1) signaling pathway, necessary for murine embryonic arterial endothelium differentiation, vascular development, and early lung morphogenesis [60, 61]. Interestingly, mutations in genes in the BMP pathway are associated with PAH , which is often part of the clinical manifestation of ACDMPV. Importantly, TMEM100 is also a potential downstream target of FOXF1 detected by ChIP-seq in mice overexpressing Foxf1 , suggesting that FOXF1 abnormalities may trigger expression changes of this gene in ACDMPV patients.
Among the most prominent pathways with the over-represented genes detected in one or more studied ACDMPV groups is SHH signaling, one of the key pathways regulating lung formation . During mouse pulmonary morphogenesis, molecular signaling from epithelial SHH to the mesenchymal FOXF1 is critical for normal pulmonary vasculature development [23, 63]. The increase of SHH pathway expression in ACDMPV could represent a compensatory response to the loss of FOXF1. In addition to FOXF1, SHH is thought to control expression of the TBX genes, which are also associated with LLDD and PAH in humans [64,65,66,67,68,69] and are known regulators of lung branching in mice [70, 71]. Of note, in our recent ChIP-seq study performed in IMR-90 fetal lung fibroblasts, we demonstrated that TBX2 and/or TBX4 bind to FOXF1, FENDRR, their promoter(s) and lung-specific enhancer, as well as to genes from the axon guidance/SEMA signaling (SEMAs, PLXNs, and NRP), that were found to be deregulated in ACDMPV patients in this study. This suggests that SHH-FOXF1, SEMA-NRP, and the TBX2-TBX4 signaling may be involved in a crosstalk responsible for normal human lung development, disruption of which leads to LLDD.
We have found also a decrease of WNT signaling, especially the canonical WNT signaling to β-catenin, one of the pathways mediating crosstalk between lung epithelia and mesenchyme during branching and important for pulmonary angiogenesis .
Noteworthy is also a decrease of S1PR1 expression. S1PR1 is essential for blood/air barrier functioning and its murine orthologue was found downregulated in the Foxf1+/- mice . The detected increase of NKX2-1 expression might account for both lung and cardiac abnormalities in ACDMPV patients. NKX2-1 has been associated with severe abnormalities in early development of lungs  and found mutated in patients with cardiac septum defects . NKX2-5, which has been associated with hypoplastic left heart syndrome , was upregulated in ACDMPV patients with deletion of FOXF1 and its enhancer.
FENDRR shares with FOXF1 the distant lung-specific enhancer and promoter-containing region and its transcription was found to be positively regulated by FOXF1 . To better characterize the role of FENDRR in lung development, we have analyzed transcriptome perturbations in the IMR-90 fetal lung fibroblast cells depleted of FENDRR using siRNA or CRISPR/Cas9 since we did not have a patient-derived material with an isolated FENDRR deletion. We have chosen IMR-90 cells because of their relatively high expression of FENDRR and FOXF1 and their origin from normal fetal lungs in pseudoglandular and canalicular stages of lung development at which ACDMPV abnormalities are first observed. One of the upregulated genes in cells with knockdown and knockout of FENDRR is HES5, a downstream target of Notch signaling . High level of the HES5 protein was observed in lungs of patients with PAH. Another gene deregulated in lung fibroblasts with the reduced level of FENDRR is CSF2RA that was previously involved in pulmonary alveolar proteinosis . Comparison between fibroblasts with FENDRR knockout and knockdown and lung tissues from ACDMPV patients with FENDRR deletion revealed a common dysregulation of 65 genes. Nearly 50% of these genes have the same direction of changes in both groups whereas the other half of genes were deregulated in the opposite direction. That could be caused by a specific expression profiles among different cell types.
Interestingly, many gene expression changes identified in lung transcriptome from the ACDMPV patient with normal levels of FOXF1 and FENDRR transcripts and without any pathogenic variant involving FOXF1, FENDRR, or their distant enhancer (Group 5) overlapped with those found in ACDMPV cases related to FOXF1 deficiency, suggesting a defect involving the FOXF1 target in patient 5.1.
Here, we have expanded our previously reported targeted expression array-based ACDMPV study  by global transcriptome RNA-seq profiling. A partial overlap between both datasets was observed. However, since we have used fewer samples in our previous experiment than in the current RNA-seq analysis, the differences in gene expression levels in RNA-seq and microarray are not unexpected. Overall, RNA-seq provided more detailed information on specific transcript expression patterns, confirmed by a validation experiment using NanoString, and showed a significant correlation between both methods.
This is the first study of ACDMPV patients providing the comprehensive transcriptome profiles of samples with various genetic variants involving FOXF1 and/or its regulatory region. However, our analyses are based on a limited number of samples and further studies involving more patients should be performed to draw more robust conclusions. While the observed differences in transcriptome patterns of patients may be partially explained by variability in type and/or size of genomic defects within the 16q24.1 locus, the gene expression profiles may be affected by the discrepancies in the number of sequencing reads in each sample due to the quality differences of obtained tissues. The transcriptome profiles of the samples with a low number of reads (< 10 M) are less accurate when compared to those with a higher read depth and their inclusion could reduce the statistical power by increasing the number of false negatives. Nevertheless, we expect that a low number of reads in the selected samples should not impact the false positive rate that may be much better controlled by a strict algorithm applied for the selection of genes with differential expression.
Our transcriptomic analyses provide further insight into the complexity of the network of signaling pathways affected by FOXF1 deficiency in ACDMPV patients. Overall, consistent with the absence of major phenotypic differences among patients with CNV deletions and point mutations in FOXF1 , our RNA-seq analyses further indicate that FOXF1 deficiency is pathogenic for ACDMPV. The partial overlap of transcriptome changes between ACDMPV lungs and FENDRR-depleted IMR-90 cells suggests the involvement of this lncRNA in the ACDMPV etiology.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Alveolar capillary dysplasia with misalignment of pulmonary veins
- Array CGH:
Array comparative genomic hybridization
Baylor College of Medicine
Forkhead-box family transcription factor
Long noncoding RNA
Lethal lung developmental disorder
Pulmonary arterial hypertension
Single guide RNA
Single nucleotide variant
Single nucleotide polymorphism
Vascular endothelial growth factor
Janney CG, Askin FB, Kuhn C. Congenital alveolar capillary dysplasia—an unusual cause of respiratory distress in the newborn. Am J Clin Pathol. 1981;76:722–7.
Langston C. Misalignment of pulmonary veins and alveolar capillary dysplasia. Pediatr Pathol. 1991;11:163–70.
Sen P, Thakur N, Stockton DW, Langston C, Bejjani BA. Expanding the phenotype of alveolar capillary dysplasia (ACD). J Pediatr. 2004;145:646–51.
Bishop NB, Stankiewicz P, Steinhorn RH. Alveolar capillary dysplasia. Am J Respir Crit Care Med. 2011;184:172–9.
Galambos C, Sims-Lucas S, Ali N, Gien J, Dishop MK, Abman SH. Intrapulmonary vascular shunt pathways in alveolar capillary dysplasia with misalignment of pulmonary veins. Thorax. 2015;70:84–5.
Slot E, Edel G, Cutz E, van Heijst A, Post M, Schnater M, et al. Alveolar capillary dysplasia with misalignment of the pulmonary veins: clinical, histological, and genetic aspects. Pulm Circ. 2018;8:2045894018795143.
Abdallah HI, Karmazin N, Marks LA. Late presentation of misalignment of lung vessels with alveolar capillary dysplasia. Crit Care Med. 1993;21:628–30.
Michalsky MP, Arca MJ, Groenman F, Hammond S, Tibboel D, Caniano DA. Alveolar capillary dysplasia: a logical approach to a fatal disease. J Pediatr Surg. 2005;40:1100–5.
Shankar V, Haque A, Johnson J, Pietsch J. Late presentation of alveolar capillary dysplasia in an infant. Pediatr Crit Care Med. 2006;7:177–9.
Ahmed S, Ackerman V, Faught P, Langston C. Profound hypoxemia and pulmonary hypertension in a 7-month-old infant: late presentation of alveolar capillary dysplasia. Pediatr Crit Care Med. 2008;9:e43-46.
Kodama Y, Tao K, Ishida F, Kawakami T, Tsuchiya K, Ishida K, et al. Long survival of congenital alveolar capillary dysplasia patient with NO inhalation and epoprostenol: effect of sildenafil, beraprost and bosentan. Pediatr Int. 2012;54:923–6.
Ito Y, Akimoto T, Cho K, Yamada M, Tanino M, Dobata T, et al. A late presenter and long-term survivor of alveolar capillary dysplasia with misalignment of the pulmonary veins. Eur J Pediatr. 2015;174:1123–6.
Towe CT, White FV, Grady RM, Sweet SC, Eghtesady P, Wegner DJ, et al. Infants with atypical presentations of alveolar capillary dysplasia with misalignment of the pulmonary veins who underwent bilateral lung transplantation. J Pediatr. 2018;194:158-164.e1.
Edwards JJ, Murali C, Pogoriler J, Frank DB, Handler SS, Deardorff MA, et al. Histopathologic and genetic features of alveolar capillary dysplasia with atypical late presentation and prolonged survival. J Pediatr. 2019;210:214–219.e2.
Szafranski P, Liu Q, Karolak JA, Song X, de Leeuw N, Faas B, et al. Association of rare non-coding SNVs in the lung-specific FOXF1 enhancer with a mitigation of the lethal ACDMPV phenotype. Hum Genet. 2019;138:1301–11.
Szafranski P, Gambin T, Dharmadhikari AV, Akdemir KC, Jhangiani SN, Schuette J, et al. Pathogenetics of alveolar capillary dysplasia with misalignment of pulmonary veins. Hum Genet. 2016;135:569–86.
Stankiewicz P, Sen P, Bhatt SS, Storer M, Xia Z, Bejjani BA, et al. Genomic and genic deletions of the FOX gene cluster on 16q24.1 and inactivating mutations of FOXF1 cause alveolar capillary dysplasia and other malformations. Am J Hum Genet. 2009;84:780–91.
Szafranski P, Dharmadhikari AV, Brosens E, Gurha P, Kolodziejska KE, Zhishuo O, et al. Small noncoding differentially methylated copy-number variants, including lncRNA genes, cause a lethal lung developmental disorder. Genome Res. 2013;23:23–33.
Murphy DB, Wiese S, Burfeind P, Schmundt D, Mattei MG, Schulz-Schaeffer W, et al. Human brain factor 1, a new member of the fork head gene family. Genomics. 1994;21:551–7.
Pierrou S, Hellqvist M, Samuelsson L, Enerbäck S, Carlsson P. Cloning and characterization of seven human forkhead proteins: binding site specificity and DNA bending. EMBO J. 1994;13:5002–12.
Mahlapuu M, Pelto-Huikko M, Aitola M, Enerbäck S, Carlsson P. FREAC-1 contains a cell-type-specific transcriptional activation domain and is expressed in epithelial-mesenchymal interfaces. Dev Biol. 1998;202:183–95.
Costa RH, Kalinichenko VV, Lim L. Transcription factors in mouse lung development and function. Am J Physiol Lung Cell Mol Physiol. 2001;280:L823-838.
Mahlapuu M, Enerbäck S, Carlsson P. Haploinsufficiency of the forkhead gene Foxf1, a target for sonic hedgehog signaling, causes lung and foregut malformations. Development. 2001;128:2397–406.
Kalinichenko VV, Lim L, Stolz DB, Shin B, Rausa FM, Clark J, et al. Defects in pulmonary vasculature and perinatal lung hemorrhage in mice heterozygous null for the Forkhead Box f1 transcription factor. Dev Biol. 2001;235:489–506.
Lim L, Kalinichenko VV, Whitsett JA, Costa RH. Fusion of lung lobes and vessels in mouse embryos heterozygous for the forkhead box f1 targeted allele. Am J Physiol Lung Cell Mol Physiol. 2002;282:L1012-1022.
Astorga J, Carlsson P. Hedgehog induction of murine vasculogenesis is mediated by Foxf1 and Bmp4. Development. 2007;134:3753–61.
Sen P, Dharmadhikari AV, Majewski T, Mohammad MA, Kalin TV, Zabielska J, et al. Comparative analyses of lung transcriptomes in patients with alveolar capillary dysplasia with misalignment of pulmonary veins and in Foxf1 heterozygous knockout mice. PLoS ONE. 2014;9: e94390.
Fernandes-Silva H, Correia-Pinto J, Moura RS. Canonical sonic hedgehog signaling in early lung development. J Dev Biol. 2017;5:3.
Ho UY, Wainwright BJ. Patched1 patterns fibroblast growth factor 10 and Forkhead box F1 expression during pulmonary branch formation. Mech Dev. 2017;147:37–48.
Szafranski P, Herrera C, Proe LA, Coffman B, Kearney DL, Popek E, et al. Narrowing the FOXF1 distant enhancer region on 16q241 critical for ACDMPV. Clin Epigenet. 2016;8:112.
Dello Russo P, Franzoni A, Baldan F, Puppin C, De Maglio G, Pittini C, et al. A 16q deletion involving FOXF1 enhancer is associated to pulmonary capillary hemangiomatosis. BMC Med Genet. 2015;16:94.
Seo H, Kim J, Park G-H, Kim Y, Cho S-W. Long-range enhancers modulate Foxf1 transcription in blood vessels of pulmonary vascular network. Histochem Cell Biol. 2016;146:289–300.
Szafranski P, Gambin T, Karolak JA, Popek E, Stankiewicz P. Lung-specific distant enhancer cis-regulates expression of FOXF1 and lncRNA FENDRR. Hum Mutat. 2021;42:694–8.
Sen P, Yang Y, Navarro C, Silva I, Szafranski P, Kolodziejska KE, et al. Novel FOXF1 mutations in sporadic and familial cases of alveolar capillary dysplasia with misaligned pulmonary veins imply a role for its DNA binding domain. Hum Mutat. 2013;34:801–11.
Prothro SL, Plosa E, Markham M, Szafranski P, Stankiewicz P, Killen SAS. Prenatal diagnosis of alveolar capillary dysplasia with misalignment of pulmonary veins. J Pediatr. 2016;170:317–8.
Szafranski P, Kośmider E, Liu Q, Karolak JA, Currie L, Parkash S, et al. LINE- and Alu-containing genomic instability hotspot at 16q24.1 associated with recurrent and nonrecurrent CNV deletions causative for ACDMPV. Hum Mutat. 2018;39:1916–25.
Concordet J-P, Haeussler M. CRISPOR: intuitive guide selection for CRISPR/Cas9 genome editing experiments and screens. Nucleic Acids Res. 2018;46:W242–5.
Metsalu T, Vilo J. ClustVis: a web tool for visualizing clustering of multivariate data using principal component analysis and heatmap. Nucleic Acids Res. 2015;43:W566-570.
Kamburov A, Stelzl U, Lehrach H, Herwig R. The ConsensusPathDB interaction database: 2013 update. Nucleic Acids Res. 2013;41:D793-800.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.
Liberzon A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov JP, Tamayo P. The molecular signatures database (MSigDB) hallmark gene set collection. Cell Syst. 2015;1:417–25.
Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7.
Hinck L. The versatile roles of “axon guidance” cues in tissue morphogenesis. Dev Cell. 2004;7:783–93.
Ito T, Kagoshima M, Sasaki Y, Li C, Udaka N, Kitsukawa T, et al. Repulsive axon guidance molecule Sema3A inhibits branching morphogenesis of fetal mouse lung. Mech Dev. 2000;97:35–45.
Kagoshima M, Ito T. Diverse gene expression and function of semaphorins in developing lung: positive and negative regulatory roles of semaphorins in lung branching morphogenesis. Genes Cells. 2001;6:559–71.
Goshima Y, Ito T, Sasaki Y, Nakamura F. Semaphorins as signals for cell repulsion and invasion. J Clin Invest. 2002;109:993–8.
Vadivel A, Alphonse RS, Collins JJP, van Haaften T, O’Reilly M, Eaton F, et al. The axonal guidance cue semaphorin 3C contributes to alveolar growth and repair. PLoS ONE. 2013;8: e67225.
Yazdani U, Terman JR. The semaphorins. Genome Biol. 2006;7:211.
Alto LT, Terman JR. Semaphorins and their signaling mechanisms. Methods Mol Biol. 2017;1493:1–25.
Cantuti-Castelvetri L, Ojha R, Pedro LD, Djannatian M, Franz J, Kuivanen S, et al. Neuropilin-1 facilitates SARS-CoV-2 cell entry and infectivity. Science. 2020;370:856–60.
Joza S, Wang J, Fox E, Hillman V, Ackerley C, Post M. Loss of semaphorin-neuropilin-1 signaling causes dysmorphic vascularization reminiscent of alveolar capillary dysplasia. Am J Pathol. 2012;181:2003–17.
Joza S, Wang J, Tseu I, Ackerley C, Post M. Fetal, but not postnatal, deletion of semaphorin-neuropilin-1 signaling affects murine alveolar development. Am J Respir Cell Mol Biol. 2013;49:627–36.
Dharmadhikari AV, Sun JJ, Gogolewski K, Carofino BL, Ustiyan V, Hill M, et al. Lethal lung hypoplasia and vascular defects in mice with conditional Foxf1 overexpression. Biol Open. 2016;5:1595–606.
Thébaud B, Ladha F, Michelakis ED, Sawicka M, Thurston G, Eaton F, et al. Vascular endothelial growth factor gene therapy increases survival, promotes lung angiogenesis, and prevents alveolar damage in hyperoxia-induced lung injury: evidence that angiogenesis participates in alveolarization. Circulation. 2005;112:2477–86.
Ren X, Ustiyan V, Pradhan A, Cai Y, Havrilak JA, Bolte CS, et al. FOXF1 transcription factor is required for formation of embryonic vasculature by regulating VEGF signaling in endothelial cells. Circ Res. 2014;115:709–20.
Raimondi C, Ruhrberg C. Neuropilin signalling in vessels, neurons and tumours. Semin Cell Dev Biol. 2013;24:172–8.
Peixoto FO, Pereira-Terra P, Moura RS, Carvalho-Dias E, Correia-Pinto J, Nogueira-Silva C. The role of ephrins-B1 and -B2 during fetal rat lung development. Cell Physiol Biochem. 2015;35:104–15.
Yu J, Bulk E, Ji P, Hascher A, Tang M, Metzger R, et al. The EPHB6 receptor tyrosine kinase is a metastasis suppressor that is frequently silenced by promoter DNA hypermethylation in non-small cell lung cancer. Clin Cancer Res. 2010;16:2275–83.
Moon E-H, Kim M-J, Ko KS, Kim YS, Seo J, Oh SP, et al. Generation of mice with a conditional and reporter allele for Tmem100. Genesis. 2010;48:673–8.
Sountoulidis A, Stavropoulos A, Giaglis S, Apostolou E, Monteiro R, Chuva de Sousa Lopes SM, et al. Activation of the canonical bone morphogenetic protein (BMP) pathway during lung morphogenesis and adult lung tissue repair. PLoS ONE. 2012;7: e41460.
Somekawa S, Imagawa K, Hayashi H, Sakabe M, Ioka T, Sato GE, et al. Tmem100, an ALK1 receptor signaling-dependent gene essential for arterial endothelium differentiation and vascular morphogenesis. Proc Natl Acad Sci U S A. 2012;109:12064–9.
Welch CL, Chung WK. Genetics and genomics of pediatric pulmonary arterial hypertension. Genes (Basel). 2020;11:E1213.
Bellusci S, Furuta Y, Rush MG, Henderson R, Winnier G, Hogan BL. Involvement of Sonic hedgehog (Shh) in mouse embryonic lung growth and morphogenesis. Development. 1997;124:53–63.
Nimmakayalu M, Major H, Sheffield V, Solomon DH, Smith RJ, Patil SR, et al. Microdeletion of 17q22q23.2 encompassing TBX2 and TBX4 in a patient with congenital microcephaly, thyroid duct cyst, sensorineural hearing loss, and pulmonary hypertension. Am J Med Genet A. 2011;155A:418–23.
Zhu N, Gonzaga-Jauregui C, Welch CL, Ma L, Qi H, King AK, et al. Exome sequencing in children with pulmonary arterial hypertension demonstrates differences compared with adults. Circ Genom Precis Med. 2018;11: e001887.
Karolak JA, Vincent M, Deutsch G, Gambin T, Cogné B, Pichon O, et al. Complex compound inheritance of lethal lung developmental disorders due to disruption of the TBX-FGF pathway. Am J Hum Genet. 2019;104:213–28.
Karolak JA, Gambin T, Honey EM, Slavik T, Popek E, Stankiewicz P. A de novo 2.2 Mb recurrent 17q23.1q23.2 deletion unmasks novel putative regulatory non-coding SNVs associated with lethal lung hypoplasia and pulmonary hypertension: a case report. BMC Med Genomics. 2020;13:34.
Suhrie K, Pajor NM, Ahlfeld SK, Dawson DB, Dufendach KR, Kitzmiller JA, et al. Neonatal lung disease associated with TBX4 mutations. J Pediatr. 2019;206:286-292.e1.
Galambos C, Mullen MP, Shieh JT, Schwerk N, Kielt MJ, Ullmann N, et al. Phenotype characterisation of TBX4 mutation and deletion carriers with neonatal and paediatric pulmonary hypertension. Eur Respir J. 2019;54:1801965.
Arora R, Metzger RJ, Papaioannou VE. Multiple roles and interactions of Tbx4 and Tbx5 in development of the respiratory system. PLoS Genet. 2012;8: e1002866.
Lüdtke TH-W, Farin HF, Rudat C, Schuster-Gossler K, Petry M, Barnett P, et al. Tbx2 controls lung growth by direct repression of the cell cycle inhibitor genes Cdkn1a and Cdkn1b. PLoS Genet. 2013;9: e1003189.
Whitsett JA, Kalin TV, Xu Y, Kalinichenko VV. Building and regenerating the lung cell by cell. Physiol Rev. 2019;99:513–54.
Cai Y, Bolte C, Le T, Goda C, Xu Y, Kalin TV, et al. FOXF1 maintains endothelial barrier function and prevents edema after lung injury. Sci Signal. 2016;9:ra40.
Serls AE, Doherty S, Parvatiyar P, Wells JM, Deutsch GH. Different thresholds of fibroblast growth factors pattern the ventral foregut into liver and lung. Development. 2005;132:35–47.
Thorwarth A, Schnittert-Hübener S, Schrumpf P, Müller I, Jyrch S, Dame C, et al. Comprehensive genotyping and clinical characterisation reveal 27 novel NKX2-1 mutations and expand the phenotypic spectrum. J Med Genet. 2014;51:375–87.
Elliott DA, Kirk EP, Yeoh T, Chandar S, McKenzie F, Taylor P, et al. Cardiac homeobox gene NKX2-5 mutations and congenital heart disease: associations with atrial septal defect and hypoplastic left heart syndrome. J Am Coll Cardiol. 2003;41:2072–6.
Li X, Zhang X, Leathers R, Makino A, Huang C, Parsa P, et al. Notch3 signaling promotes the development of pulmonary arterial hypertension. Nat Med. 2009;15:1289–97.
Suzuki T, Sakagami T, Young LR, Carey BC, Wood RE, Luisetti M, et al. Hereditary pulmonary alveolar proteinosis: pathogenesis, presentation, diagnosis, and therapy. Am J Respir Crit Care Med. 2010;182:1292–304.
Design of the study and collection, analysis, and interpretation of data and in writing the manuscript were supported by grants awarded by the US National Institutes of Health (NIH), National Heart Lung and Blood Institute (NHLBI) R01HL137203 (P.St.), National Science Centre in Poland 2019/35/D/ NZ5/02896 (J.A.K.), the Genomic and RNA Profiling Core at Baylor College of Medicine with funding from the NIH S10 grant (1S10OD023469), and from the NIH NCI grant (P30CA125123).
Ethics approval and consent to participate
The study protocol was approved by the Institutional Review Board for Human Subject Research at Baylor College of Medicine (H-8712). In compliance with the Declaration of Helsinki, informed written consent for genetic study was obtained for participated individuals.
Consent for publication
Written consent for research study was obtained from the parents.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
ACDMPV patients enrolled in transcriptomic studies.
Schematic representation of CNVs and point mutations detected within FOXF1 locus in ACDMPV patients enrolled in transcriptomic studies.
CRISPR/Cas9 knockout of FENDRR in IMR-90 cell line in IMR-90 cell line. A Schematic representation of the FENDRR locus and CRISPR/Cas9-mediated FENDRR knockout strategy. Red box represents FENDRR exon 1; white box represents maximal 5’ extension of the exon 1, caused by alternative splicing, that is present in GENCODE transcript FENDRR-206 (ENST00000599749.5). F and R indicate primers (F: 5’-GAAAATGAAGGCGTCCGAGC, R: 5’-GCACATGCAAACCCGGTAAA) used in PCR amplification from normal (WT) and edited (KO) lung fibroblasts IMR-90. B Agarose gel electrophoresis of PCR products that verified promoter-exon 1 deletion (pe1∆). There was a minor heterogeneity with respect to the exact position of each of the deletion breakpoints. Sequence through the deletion junction and deletion coordinates refer to the predominant deletion (a small fraction of the deletions was 2 to 4 nt shorter). cCRE abbreviates FENDRR putative promoter (ENCODE).
List of protein coding genes deregulated in ACDMPV lungs. Case grouping is described in Additional file 1.
List of 94 genes commonly deregulated in all ACDMPV individuals.
Molecular pathways affected in ACDMPV lungs, identified by RNA-seq (ConsensusPathDB analysis). Case grouping is described in Additional file 1.
Gene ontology analysis of differentially expressed genes in each group. A bar chart of the most significant gene ontology (GO) biological processes (left panel) and molecular function (right panel) terms for differentially expressed genes.
Heatmap indicating hierarchical clustering of ACDMPV and non-ACDMPV control samples based on the RNA-seq expression values of genes from axon guidance (R-HSA-422475).
Heatmap indicating hierarchical clustering of ACDMPV and non-ACDMPV control samples based on the RNA-seq expression values of genes from semaphorin interaction (R-HSA-373755).
Heatmap indicating hierarchical clustering of ACDMPV and non-ACDMPV control samples based on the RNA-seq expression values of genes from VEGF signaling (R-HSA-194138).
Molecular pathways affected in IMR-90 cells with FENDRR depletion, identified by RNA-seq (ConsensusPathDB analysis).
Mutual comparison of RNA-seq and NanoString nCounter-determined changes in ACDMPV transcriptomes. Both methods detected relatively similar changes in expression levels of selected 120 genes.
Comparison of RNA-seq data with microarray-based expression results (GSE54780). Case grouping is described in Additional file 1.
Bio-GPS graphs demonstrating similar fetal lung- and lung-specific expression pattern of TMEM100 and FOXF1.
About this article
Cite this article
Karolak, J.A., Gambin, T., Szafranski, P. et al. Perturbation of semaphorin and VEGF signaling in ACDMPV lungs due to FOXF1 deficiency. Respir Res 22, 212 (2021). https://doi.org/10.1186/s12931-021-01797-7
- Vascular development
- Lung branching
- Lethal lung developmental disorders