Machine learning in time-lapse imaging to differentiate embryos from young vs old mice.

Time-lapse microscopy for embryos is a non-invasive technology used to characterize early embryo development. This study employs time-lapse microscopy and machine learning to elucidate changes in embryonic growth kinetics with maternal aging. We analyzed morphokinetic parameters of embryos from young and aged C57BL6/NJ mice via continuous imaging. Our findings show that aged embryos accelerated through cleavage stages (from 5-cells) to morula compared to younger counterparts, with no significant differences observed in later stages of blastulation. Unsupervised machine learning identified two distinct clusters comprising of embryos from aged or young donors. Moreover, in supervised learning, the XGBoost (extreme gradient boosting) algorithm successfully predicted the age-related phenotype with 0.78 accuracy, 0.81 precision, and 0.83 recall following hyperparameter tuning. These results highlight two main scientific insights: maternal aging affects embryonic development pace, and that AI can differentiate between embryos from aged and young maternal mice by a non-invasive approach. Thus, machine learning can be used to identify morphokinetics phenotypes for further studies. This study has potential for future applications in selecting human embryos for embryo transfer, without or in complement with preimplantation genetic testing.

A GREB1-steroid receptor feedforward mechanism governs differential GREB1 action in endometrial function and endometriosis.

Cellular responses to the steroid hormones, estrogen (E2), and progesterone (P4) are governed by their cognate receptor's transcriptional output. However, the feed-forward mechanisms that shape cell-type-specific transcriptional fulcrums for steroid receptors are unidentified. Herein, we found that a common feed-forward mechanism between GREB1 and steroid receptors regulates the differential effect of GREB1 on steroid hormones in a physiological or pathological context. In physiological (receptive) endometrium, GREB1 controls P4-responses in uterine stroma, affecting endometrial receptivity and decidualization, while not affecting E2-mediated epithelial proliferation. Of mechanism, progesterone-induced GREB1 physically interacts with the progesterone receptor, acting as a cofactor in a positive feedback mechanism to regulate P4-responsive genes. Conversely, in endometrial pathology (endometriosis), E2-induced GREB1 modulates E2-dependent gene expression to promote the growth of endometriotic lesions in mice. This differential action of GREB1 exerted by a common feed-forward mechanism with steroid receptors advances our understanding of mechanisms that underlie cell- and tissue-specific steroid hormone actions.

Generation of a humanized mAce2 and a conditional hACE2 mouse models permissive to SARS-COV-2 infection.

The Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) remains a public health concern and a subject of active research effort. Development of pre-clinical animal models is critical to study viral-host interaction, tissue tropism, disease mechanisms, therapeutic approaches, and long-term sequelae of infection. Here, we report two mouse models for studying SARS-CoV-2: A knock-in mAce2F83Y,H353K mouse that expresses a mouse-human hybrid form of the angiotensin-converting enzyme 2 (ACE2) receptor under the endogenous mouse Ace2 promoter, and a Rosa26 conditional knock-in mouse carrying the human ACE2 allele (Rosa26hACE2). Although the mAce2F83Y,H353K mice were susceptible to intranasal inoculation with SARS-CoV-2, they did not show gross phenotypic abnormalities. Next, we generated a Rosa26hACE2;CMV-Cre mouse line that ubiquitously expresses the human ACE2 receptor. By day 3 post infection with SARS-CoV-2, Rosa26hACE2;CMV-Cre mice showed significant weight loss, a variable degree of alveolar wall thickening and reduced survival rates. Viral load measurements confirmed inoculation in lung and brain tissues of infected Rosa26hACE2;CMV-Cre mice. The phenotypic spectrum displayed by our different mouse models translates to the broad range of clinical symptoms seen in the human patients and can serve as a resource for the community to model and explore both treatment strategies and long-term consequences of SARS-CoV-2 infection.

An oocyte-specific Cas9-expressing mouse for germline CRISPR/Cas9-mediated genome editing.

Cas9 transgenes can be employed for genome editing in mouse zygotes. However, using transgenic instead of exogenous Cas9 to produce gene-edited animals creates unique issues including ill-defined transgene integration sites, the potential for prolonged Cas9 expression in transgenic embryos, and increased genotyping burden. To overcome these issues, we generated mice harboring an oocyte-specific, Gdf9 promoter driven, Cas9 transgene (Gdf9-Cas9) targeted as a single copy into the Hprt1 locus. The X-linked Hprt1 locus was selected because it is a defined integration site that does not influence transgene expression, and breeding of transgenic males generates obligate transgenic females to serve as embryo donors. Using microinjections and electroporation to introduce sgRNAs into zygotes derived from transgenic dams, we demonstrate that Gdf9-Cas9 mediates genome editing as efficiently as exogenous Cas9 at several loci. We show that genome editing efficiency is independent of transgene inheritance, verifying that maternally derived Cas9 facilitates genome editing. We also show that paternal inheritance of Gdf9-Cas9 does not mediate genome editing, confirming that Gdf9-Cas9 is not expressed in embryos. Finally, we demonstrate that off-target mutagenesis is equally rare when using transgenic or exogenous Cas9. Together, these results show that the Gdf9-Cas9 transgene is a viable alternative to exogenous Cas9.

Delayed skeletal development and IGF-1 deficiency in a mouse model of lysinuric protein intolerance.

SLC7A7 deficiency, or lysinuric protein intolerance (LPI), causes loss of function of the y+LAT1 transporter critical for efflux of arginine, lysine and ornithine in certain cells. LPI is characterized by urea cycle dysfunction, renal disease, immune dysregulation, growth failure, delayed bone age and osteoporosis. We previously reported that Slc7a7 knockout mice (C57BL/6×129/SvEv F2) recapitulate LPI phenotypes, including growth failure. Our main objective in this study was to characterize the skeletal phenotype in these mice. Compared to wild-type littermates, juvenile Slc7a7 knockout mice demonstrated 70% lower body weights, 87% lower plasma IGF-1 concentrations and delayed skeletal development. Because poor survival prevents evaluation of mature knockout mice, we generated a conditional Slc7a7 deletion in mature osteoblasts or mesenchymal cells of the osteo-chondroprogenitor lineage, but no differences in bone architecture were observed. Overall, global Slc7a7 deficiency caused growth failure with low plasma IGF-1 concentrations and delayed skeletal development, but Slc7a7 deficiency in the osteoblastic lineage was not a major contributor to these phenotypes. Future studies utilizing additional tissue-specific Slc7a7 knockout models may help dissect cell-autonomous and non-cell-autonomous mechanisms underlying phenotypes in LPI.

Targeting AAV vectors to the central nervous system by engineering capsid-receptor interactions that enable crossing of the blood-brain barrier.

Viruses have evolved the ability to bind and enter cells through interactions with a wide variety of cell macromolecules. We engineered peptide-modified adeno-associated virus (AAV) capsids that transduce the brain through the introduction of de novo interactions with 2 proteins expressed on the mouse blood-brain barrier (BBB), LY6A or LY6C1. The in vivo tropisms of these capsids are predictable as they are dependent on the cell- and strain-specific expression of their target protein. This approach generated hundreds of capsids with dramatically enhanced central nervous system (CNS) tropisms within a single round of screening in vitro and secondary validation in vivo thereby reducing the use of animals in comparison to conventional multi-round in vivo selections. The reproducible and quantitative data derived via this method enabled both saturation mutagenesis and machine learning (ML)-guided exploration of the capsid sequence space. Notably, during our validation process, we determined that nearly all published AAV capsids that were selected for their ability to cross the BBB in mice leverage either the LY6A or LY6C1 protein, which are not present in primates. This work demonstrates that AAV capsids can be directly targeted to specific proteins to generate potent gene delivery vectors with known mechanisms of action and predictable tropisms.

Comprehensive ECG reference intervals in C57BL/6N substrains provide a generalizable guide for cardiac electrophysiology studies in mice.

Reference ranges provide a powerful tool for diagnostic decision-making in clinical medicine and are enormously valuable for understanding normality in pre-clinical scientific research that uses in vivo models. As yet, there are no published reference ranges for electrocardiography (ECG) in the laboratory mouse. The first mouse-specific reference ranges for the assessment of electrical conduction are reported herein generated from an ECG dataset of unprecedented scale. International Mouse Phenotyping Consortium data from over 26,000 conscious or anesthetized C57BL/6N wildtype control mice were stratified by sex and age to develop robust ECG reference ranges. Interesting findings include that heart rate and key elements from the ECG waveform (RR-, PR-, ST-, QT-interval, QT corrected, and QRS complex) demonstrate minimal sexual dimorphism. As expected, anesthesia induces a decrease in heart rate and was shown for both inhalation (isoflurane) and injectable (tribromoethanol) anesthesia. In the absence of pharmacological, environmental, or genetic challenges, we did not observe major age-related ECG changes in C57BL/6N-inbred mice as the differences in the reference ranges of 12-week-old compared to 62-week-old mice were negligible. The generalizability of the C57BL/6N substrain reference ranges was demonstrated by comparison with ECG data from a wide range of non-IMPC studies. The close overlap in data from a wide range of mouse strains suggests that the C57BL/6N-based reference ranges can be used as a robust and comprehensive indicator of normality. We report a unique ECG reference resource of fundamental importance for any experimental study of cardiac function in mice.

Generation of a novel Stra8-driven Cre recombinase strain for use in pre-meiotic germ cells in mice†.

The development of oocytes occurs over a broad time frame, starting at the earliest stages of embryogenesis and continuing into adulthood. Conditional knockout technologies such as the Cre/loxP recombination system are useful for analyzing oocyte development at specific stages, but not every time frame has appropriate Cre drivers, for instance, during oocyte meiotic initiation through early prophase I in the embryo. Here, we generated a novel knockin mouse line that produces a bicistronic transcript from the endogenous Stra8 locus that includes a "self-cleaving" 2A peptide upstream of cre. This allows for high efficiency cleavage and production of both proteins individually and results in expression of cre in both male and female gonads at the biologically relevant stage. Fluorescent reporter analysis confirms that this line recapitulates endogenous Stra8 expression in both sexes and does not affect fertility of heterozygous nor homozygous mice. This line, named Stra8P2Acre, adds to the repertoire of germ-cell specific cre driver lines and, importantly, allows for deletion of target genes during key embryonic oocyte developmental stages, including early events in meiosis. Summary Sentence Generation of a novel cre recombinase knockin to the Stra8 locus allows production of Stra8 and cre without affecting fertility.

Whole genome analysis for 163 gRNAs in Cas9-edited mice reveals minimal off-target activity.

Genome editing with CRISPR-associated (Cas) proteins holds exceptional promise for "correcting" variants causing genetic disease. To realize this promise, off-target genomic changes cannot occur during the editing process. Here, we use whole genome sequencing to compare the genomes of 50 Cas9-edited founder mice to 28 untreated control mice to assess the occurrence of S. pyogenes Cas9-induced off-target mutagenesis. Computational analysis of whole-genome sequencing data detects 26 unique sequence variants at 23 predicted off-target sites for 18/163 guides used. While computationally detected variants are identified in 30% (15/50) of Cas9 gene-edited founder animals, only 38% (10/26) of the variants in 8/15 founders validate by Sanger sequencing. In vitro assays for Cas9 off-target activity identify only two unpredicted off-target sites present in genome sequencing data. In total, only 4.9% (8/163) of guides tested have detectable off-target activity, a rate of 0.2 Cas9 off-target mutations per founder analyzed. In comparison, we observe ~1,100 unique variants in each mouse regardless of genome exposure to Cas9 indicating off-target variants comprise a small fraction of genetic heterogeneity in Cas9-edited mice. These findings will inform future design and use of Cas9-edited animal models as well as provide context for evaluating off-target potential in genetically diverse patient populations.

A CRISPR/Cas9-engineered mouse carrying a conditional knockout allele for the early growth response-1 transcription factor.

Early growth response 1 (EGR1) mediates transcriptional programs that are indispensable for cell division, differentiation, and apoptosis in numerous physiologies and pathophysiologies. Whole-body EGR1 knockouts in mice (Egr1KO ) have advanced our understanding of EGR1 function in an in vivo context. To extend the utility of the mouse to investigate EGR1 responses in a tissue- and/or cell-type-specific manner, we generated a mouse model in which exon 2 of the mouse Egr1 gene is floxed by CRISPR/Cas9 engineering. The floxed Egr1 alleles (Egr1f/f ) are designed to enable spatiotemporal control of Cre-mediated EGR1 ablation in the mouse. To confirm that the Egr1f/f alleles can be abrogated using a Cre driver, we crossed the Egr1f/f mouse with a global Cre driver to generate the Egr1 conditional knockout (Egr1d/d ) mouse in which EGR1 expression is ablated in all tissues. Genetic and protein analysis confirmed the absence of exon 2 and loss of EGR1 expression in the Egr1d/d mouse, respectively. Moreover, the Egr1d/d female exhibits overt reproductive phenotypes previously reported for the Egr1KO mouse. Therefore, studies described in this short technical report underscore the potential utility of the murine Egr1 floxed allele to further resolve EGR1 function at a tissue- and/or cell-type-specific level.

Genome-wide screening reveals the genetic basis of mammalian embryonic eye development.

BACKGROUND: Microphthalmia, anophthalmia, and coloboma (MAC) spectrum disease encompasses a group of eye malformations which play a role in childhood visual impairment. Although the predominant cause of eye malformations is known to be heritable in nature, with 80% of cases displaying loss-of-function mutations in the ocular developmental genes OTX2 or SOX2, the genetic abnormalities underlying the remaining cases of MAC are incompletely understood. This study intended to identify the novel genes and pathways required for early eye development. Additionally, pathways involved in eye formation during embryogenesis are also incompletely understood. This study aims to identify the novel genes and pathways required for early eye development through systematic forward screening of the mammalian genome. RESULTS: Query of the International Mouse Phenotyping Consortium (IMPC) database (data release 17.0, August 01, 2022) identified 74 unique knockout lines (genes) with genetically associated eye defects in mouse embryos. The vast majority of eye abnormalities were small or absent eyes, findings most relevant to MAC spectrum disease in humans. A literature search showed that 27 of the 74 lines had previously published knockout mouse models, of which only 15 had ocular defects identified in the original publications. These 12 previously published gene knockouts with no reported ocular abnormalities and the 47 unpublished knockouts with ocular abnormalities identified by the IMPC represent 59 genes not previously associated with early eye development in mice. Of these 59, we identified 19 genes with a reported human eye phenotype. Overall, mining of the IMPC data yielded 40 previously unimplicated genes linked to mammalian eye development. Bioinformatic analysis showed that several of the IMPC genes colocalized to several protein anabolic and pluripotency pathways in early eye development. Of note, our analysis suggests that the serine-glycine pathway producing glycine, a mitochondrial one-carbon donator to folate one-carbon metabolism (FOCM), is essential for eye formation. CONCLUSIONS: Using genome-wide phenotype screening of single-gene knockout mouse lines, STRING analysis, and bioinformatic methods, this study identified genes heretofore unassociated with MAC phenotypes providing models to research novel molecular and cellular mechanisms involved in eye development. These findings have the potential to hasten the diagnosis and treatment of this congenital blinding disease.

In vivo editing of the pan-endothelium by immunity evading simian adenoviral vector.

Biological applications deriving from the clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 site-specific nuclease system continue to impact and accelerate gene therapy strategies. Safe and effective in vivo co-delivery of the CRISPR/Cas9 system to target somatic cells is essential in the clinical therapeutic context. Both non-viral and viral vector systems have been applied for this delivery matter. Despite elegant proof-of-principle studies, available vector technologies still face challenges that restrict the application of CRISPR/Cas9-facilitated gene therapy. Of note, the mandated co-delivery of the gene-editing components must be accomplished in the potential presence of pre-formed anti-vector immunity. Additionally, methods must be sought to limit the potential of off-target editing. To this end, we have exploited the molecular promiscuities of adenovirus (Ad) to address the key requirements of CRISPR/Cas9-facilitated gene therapy. In this regard, we have endeavored capsid engineering of a simian (chimpanzee) adenovirus isolate 36 (SAd36) to achieve targeted modifications of vector tropism. The SAd36 vector with the myeloid cell-binding peptide (MBP) incorporated in the capsid has allowed selective in vivo modifications of the vascular endothelium. Importantly, vascular endothelium can serve as an effective non-hepatic cellular source of deficient serum factors relevant to several inherited genetic disorders. In addition to allowing for re-directed tropism, capsid engineering of nonhuman primate Ads provide the means to circumvent pre-formed vector immunity. Herein we have generated a SAd36. MBP vector that can serve as a single intravenously administered agent allowing effective and selective in vivo editing for endothelial target cells of the mouse spleen, brain and kidney. DATA AVAILABILITY: The data that support the findings of this study are available from the corresponding author upon reasonable request.

Analysis of genome-wide knockout mouse database identifies candidate ciliopathy genes.

We searched a database of single-gene knockout (KO) mice produced by the International Mouse Phenotyping Consortium (IMPC) to identify candidate ciliopathy genes. We first screened for phenotypes in mouse lines with both ocular and renal or reproductive trait abnormalities. The STRING protein interaction tool was used to identify interactions between known cilia gene products and those encoded by the genes in individual knockout mouse strains in order to generate a list of "candidate ciliopathy genes." From this list, 32 genes encoded proteins predicted to interact with known ciliopathy proteins. Of these, 25 had no previously described roles in ciliary pathobiology. Histological and morphological evidence of phenotypes found in ciliopathies in knockout mouse lines are presented as examples (genes Abi2, Wdr62, Ap4e1, Dync1li1, and Prkab1). Phenotyping data and descriptions generated on IMPC mouse line are useful for mechanistic studies, target discovery, rare disease diagnosis, and preclinical therapeutic development trials. Here we demonstrate the effective use of the IMPC phenotype data to uncover genes with no previous role in ciliary biology, which may be clinically relevant for identification of novel disease genes implicated in ciliopathies.

Mendelian gene identification through mouse embryo viability screening.

BACKGROUND: The diagnostic rate of Mendelian disorders in sequencing studies continues to increase, along with the pace of novel disease gene discovery. However, variant interpretation in novel genes not currently associated with disease is particularly challenging and strategies combining gene functional evidence with approaches that evaluate the phenotypic similarities between patients and model organisms have proven successful. A full spectrum of intolerance to loss-of-function variation has been previously described, providing evidence that gene essentiality should not be considered as a simple and fixed binary property. METHODS: Here we further dissected this spectrum by assessing the embryonic stage at which homozygous loss-of-function results in lethality in mice from the International Mouse Phenotyping Consortium, classifying the set of lethal genes into one of three windows of lethality: early, mid, or late gestation lethal. We studied the correlation between these windows of lethality and various gene features including expression across development, paralogy and constraint metrics together with human disease phenotypes. We explored a gene similarity approach for novel gene discovery and investigated unsolved cases from the 100,000 Genomes Project. RESULTS: We found that genes in the early gestation lethal category have distinct characteristics and are enriched for genes linked with recessive forms of inherited metabolic disease. We identified several genes sharing multiple features with known biallelic forms of inborn errors of the metabolism and found signs of enrichment of biallelic predicted pathogenic variants among early gestation lethal genes in patients recruited under this disease category. We highlight two novel gene candidates with phenotypic overlap between the patients and the mouse knockouts. CONCLUSIONS: Information on the developmental period at which embryonic lethality occurs in the knockout mouse may be used for novel disease gene discovery that helps to prioritise variants in unsolved rare disease cases.

Identifying genetic determinants of inflammatory pain in mice using a large-scale gene-targeted screen.

ABSTRACT: Identifying the genetic determinants of pain is a scientific imperative given the magnitude of the global health burden that pain causes. Here, we report a genetic screen for nociception, performed under the auspices of the International Mouse Phenotyping Consortium. A biased set of 110 single-gene knockout mouse strains was screened for 1 or more nociception and hypersensitivity assays, including chemical nociception (formalin) and mechanical and thermal nociception (von Frey filaments and Hargreaves tests, respectively), with or without an inflammatory agent (complete Freund's adjuvant). We identified 13 single-gene knockout strains with altered nocifensive behavior in 1 or more assays. All these novel mouse models are openly available to the scientific community to study gene function. Two of the 13 genes (Gria1 and Htr3a) have been previously reported with nociception-related phenotypes in genetically engineered mouse strains and represent useful benchmarking standards. One of the 13 genes (Cnrip1) is known from human studies to play a role in pain modulation and the knockout mouse reported herein can be used to explore this function further. The remaining 10 genes (Abhd13, Alg6, BC048562, Cgnl1, Cp, Mmp16, Oxa1l, Tecpr2, Trim14, and Trim2) reveal novel pathways involved in nociception and may provide new knowledge to better understand genetic mechanisms of inflammatory pain and to serve as models for therapeutic target validation and drug development.

GSDMB is increased in IBD and regulates epithelial restitution/repair independent of pyroptosis.

Gasdermins are a family of structurally related proteins originally described for their role in pyroptosis. Gasdermin B (GSDMB) is currently the least studied, and while its association with genetic susceptibility to chronic mucosal inflammatory disorders is well established, little is known about its functional relevance during active disease states. Herein, we report increased GSDMB in inflammatory bowel disease, with single-cell analysis identifying epithelial specificity to inflamed colonocytes/crypt top colonocytes. Surprisingly, mechanistic experiments and transcriptome profiling reveal lack of inherent GSDMB-dependent pyroptosis in activated epithelial cells and organoids but instead point to increased proliferation and migration during in vitro wound closure, which arrests in GSDMB-deficient cells that display hyper-adhesiveness and enhanced formation of vinculin-based focal adhesions dependent on PDGF-A-mediated FAK phosphorylation. Importantly, carriage of disease-associated GSDMB SNPs confers functional defects, disrupting epithelial restitution/repair, which, altogether, establishes GSDMB as a critical factor for restoration of epithelial barrier function and the resolution of inflammation.

AAV5 delivery of CRISPR-Cas9 supports effective genome editing in mouse lung airway.

Genome editing in the lung has the potential to provide long-term expression of therapeutic protein to treat lung genetic diseases. Yet efficient delivery of CRISPR to the lung remains a challenge. The NIH Somatic Cell Genome Editing (SCGE) Consortium is developing safe and effective methods for genome editing in disease tissues. Methods developed by consortium members are independently validated by the SCGE small animal testing center to establish rigor and reproducibility. We have developed and validated a dual adeno-associated virus (AAV) CRISPR platform that supports effective editing of a lox-stop-lox-Tomato reporter in mouse lung airway. After intratracheal injection of the AAV serotype 5 (AAV5)-packaged S. pyogenes Cas9 (SpCas9) and single guide RNAs (sgRNAs), we observed ∼19%-26% Tomato-positive cells in both large and small airways, including club and ciliated epithelial cell types. This highly effective AAV delivery platform will facilitate the study of therapeutic genome editing in the lung and other tissue types.

COPB2 loss of function causes a coatopathy with osteoporosis and developmental delay.

Coatomer complexes function in the sorting and trafficking of proteins between subcellular organelles. Pathogenic variants in coatomer subunits or associated factors have been reported in multi-systemic disorders, i.e., coatopathies, that can affect the skeletal and central nervous systems. We have identified loss-of-function variants in COPB2, a component of the coatomer complex I (COPI), in individuals presenting with osteoporosis, fractures, and developmental delay of variable severity. Electron microscopy of COPB2-deficient subjects' fibroblasts showed dilated endoplasmic reticulum (ER) with granular material, prominent rough ER, and vacuoles, consistent with an intracellular trafficking defect. We studied the effect of COPB2 deficiency on collagen trafficking because of the critical role of collagen secretion in bone biology. COPB2 siRNA-treated fibroblasts showed delayed collagen secretion with retention of type I collagen in the ER and Golgi and altered distribution of Golgi markers. copb2-null zebrafish embryos showed retention of type II collagen, disorganization of the ER and Golgi, and early larval lethality. Copb2+/- mice exhibited low bone mass, and consistent with the findings in human cells and zebrafish, studies in Copb2+/- mouse fibroblasts suggest ER stress and a Golgi defect. Interestingly, ascorbic acid treatment partially rescued the zebrafish developmental phenotype and the cellular phenotype in Copb2+/- mouse fibroblasts. This work identifies a form of coatopathy due to COPB2 haploinsufficiency, explores a potential therapeutic approach for this disorder, and highlights the role of the COPI complex as a regulator of skeletal homeostasis.

Perturbation of semaphorin and VEGF signaling in ACDMPV lungs due to FOXF1 deficiency.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.