Fgf10-CRISPR mosaic mutants demonstrate the gene dose-related loss of the accessory lobe and decrease in the number of alveolar type 2 epithelial cells in mouse lung

CRISPR/Cas9-mediated gene editing often generates founder generation (F0) mice that exhibit somatic mosaicism in the targeted gene(s). It has been known that Fibroblast growth factor 10 (Fgf10)-null mice exhibit limbless and lungless phenotypes, while intermediate limb phenotypes (variable defective limbs) are observed in the Fgf10-CRISPR F0 mice. However, how the lung phenotype in the Fgf10-mosaic mutants is related to the limb phenotype and genotype has not been investigated. In this study, we examined variable lung phenotypes in the Fgf10-targeted F0 mice to determine if the lung phenotype was correlated with percentage of functional Fgf10 genotypes. Firstly, according to a previous report, Fgf10-CRISPR F0 embryos on embryonic day 16.5 (E16.5) were classified into three types: type I, no limb; type II, limb defect; and type III, normal limbs. Cartilage and bone staining showed that limb truncations were observed in the girdle, (type I), stylopodial, or zeugopodial region (type II). Deep sequencing of the Fgf10-mutant genomes revealed that the mean proportion of codons that encode putative functional FGF10 was 8.3 ± 6.2% in type I, 25.3 ± 2.7% in type II, and 54.3 ± 9.5% in type III (mean ± standard error of the mean) mutants at E16.5. Histological studies showed that almost all lung lobes were absent in type I embryos. The accessory lung lobe was often absent in type II embryos with other lobes dysplastic. All lung lobes formed in type III embryos. The number of terminal tubules was significantly lower in type I and II embryos, but unchanged in type III embryos. To identify alveolar type 2 epithelial (AECII) cells, known to be reduced in the Fgf10-heterozygous mutant, immunostaining using anti-surfactant protein C (SPC) antibody was performed: In the E18.5 lungs, the number of AECII was correlated to the percentage of functional Fgf10 genotypes. These data suggest the Fgf10 gene dose-related loss of the accessory lobe and decrease in the number of alveolar type 2 epithelial cells in mouse lung. Since dysfunction of AECII cells has been implicated in the pathogenesis of parenchymal lung diseases, the Fgf10-CRISPR F0 mouse would present an ideal experimental system to explore it.


Introduction
The recently developed CRISPR/Cas9 system provides a highly efficient means for editing the genomes of model and non-model organisms. This powerful tool can help elucidate pathophysiological mechanisms underlying various genetic diseases (reviewed in [1,2]). CRISPR/ Cas9 promises the possibility of an ultimate cure for genetic diseases by enabling replacement of mutated genes with normal alleles (reviewed in [3]); however, genome-edited founder mice often exhibit somatic mosaicism in the targeted gene, meaning more than two mutated alleles for that gene are mixed in the same mouse [4][5][6]. Such mosaicism may be undesirable when it complicates phenotypic analysis [7]. However, mosaic mice may help the study of genes whose constitutive mutations are lethal, or the study of interactions between mutant and normal cells in the same individual [8]. Furthermore, recent DNA sequencing advances have enabled us to identify genetic mosaicism even in phenotypically normal individuals, the pathological significance of which is still unclear (reviewed in [9,10]).
Previously, we reported that somatic mosaicism of a homeobox gene Pax6 mutation causes variable eye phenotypes correlated with its gene dosage [11]. In this study, we focused on the variable phenotypes of genome-edited Fgf10 founder generation (F0) mice. The Fgf10 gene is required for limb and lung formation; Fgf10 null mutants suffer embryonic lethality at birth due to lung agenesis [12,13]. Recent studies have shown that a fraction of Fgf10 genome-edited F0 mice exhibits typical limbless and lungless phenotypes, and that the severity of the limbless phenotype depends on the Fgf10 mutation rate [7,14,15]. However, the lung phenotype has not been described. In this study, we sought to determine the relationship between lung and limb phenotypes, and correlation of these phenotypes to putative functional FGF10 dosage.

Materials and methods
Reagents and equipment used in this study are listed in S1 Table.

Mice and ethical statement
The animal experimental design was approved by the Committee of Animal Experiments of Tokushima University, Tokushima, Japan, and by the Animal Care and Use Committee, Okayama University, Okayama, Japan (Permit numbers: T27-16, T30-8, OKU-2017404, OKU2018605). All surgeries were performed under reversible anesthesia with medetomidinemidazolam-butorphanol and atipamezole, and all efforts were made to minimize suffering. In vitro fertilized (IVF) eggs were prepared from male and female BDF1 (C57BL/6 x DBA2 F1) mice as previously reported [14]. CRISPR/Cas9 genome editing was performed on one-cell stage zygotes by microinjection (analysis on E16.5 embryos) with Cas9 mRNA and sgRNA targeted at exon 3 of the Fgf10 gene according to Yasue et al. (2014) [14] (S1 Fig). After culturing genome-edited zygotes until the two-cell stage, embryos were transferred to the uterine tube of foster mice (Jcl: MCH (ICR) strain developed by CLEA Japan, Inc. (Tokyo, Japan); Jcl for Japan clea, MCH for Multi-Cross Hybrid, and ICR for Institute of Cancer Research) and allowed to develop until E16. 5. To obtain E18.5 F0 embryos, Cas9 mRNA (100 ng/μL) and synthesized guide RNA (50 ng/μL) were electroporated into one-cell stage IVF eggs and further manipulation of embryos was performed according to Hashimoto et al. (2016) [7]. CRISPR-control embryos were obtained by electroporation of medium only into IVF eggs before processing them for sequencing and histological examination. To know developmental variations, the weight of E18.5 embryos was measured with a standard electronic scale.

Enzyme mismatch cleavage assay
Approximately 25 μg genomic DNA were extracted from neck skin tissue (epidermis and dermis) of collected embryos using a Qiagen DNeasy Blood & Tissue Kit. The genomic region encompassing the CRISPR/Cas9 target site was amplified by polymerase chain reaction (PCR) using a primer set specific for the Fgf10 gene (S3A Fig, S1 Table) [14]. The amplicons were processed for a mismatch cleavage assay using the Guide-it Mutation Detection Kit and the products were analyzed by 2% agarose gel electrophoresis (S3B Fig).

DNA sequencing
To verify the presence of the target site mutation, single bands of PCR amplicons of approximately 500 bp derived from the collected embryos at E16.5 (neck skin tissues) were subjected to Sanger sequencing. The resulting DNA sequence waves were dissociated using a web tool for TIDER data analysis available at https://tide.nki.nl [16].
Deep sequencing analysis using DNA extracted from neck skin, lung, and/or limb dermal tissues (Table 1 for E18.5; neck skin only for E16.5 embryos in Table 3; neck skin only for additional E18.5 embryos in S3 Table) for the Fgf10 target site was performed according to the previous report on Illumina MiSeq system [17]. Briefly, the Fgf10 on-target site from each individual embryo was amplified using custom barcode primers shown in S1 Table. Library construction and sequencing were performed at the National Institute for Basic Biology (NIBB, Japan) and Bioengineering Lab (Kanagawa, Japan). The pooled sequence data was demultiplexed into each sample and analyzed subsequently using a web-based tool, CLiCKAR [17].

Correction and calculation of putative somatic mutation rate in the case of large insertions
PCR amplicons for deep sequencing were analyzed using TAPEStation 4200 system. When larger DNA bands were detected above the main bands, their approximate molecular size was recorded. "Peak Molarity" value for each DNA band was documented and the percentage of large insertions was estimated by calculating the total value. The putative somatic mutation rate for large insertions was calculated proportionally. In embryos #3 (type II, E16.5), #23 (type III, E18.5), #11_2 (type II, E18.5), and #22_2 (type III, E18.5), large insertions were detected and corrected accordingly.

Cartilage and bone staining
After Cesarean dissection, E18.5 embryos, removed their skin and internal organs, were fixed in 95% ethanol and processed for whole mount skeletal staining with Alcian blue 8GX and Alizarin red S according to standard procedures [18].    embedded embryos and mounted on microscope slides. After deparaffinization and rehydration, sections were stained with Mayer's hematoxylin solution and 0.125% eosin. Following staining, sections were washed, dehydrated, and mounted with DPX. After HE staining, photos were taken with a 10x objective lens for each embryonic lung (n = 3 for each type, at E16.5) (Fig 5E-5H) and tiled using ImageJ (https://imagej.nih.gov/ij/index.html)) (RRID: SCR_003070). Terminal tubules were scored in the tiled full-field image (n = 1 for each embryo). The area for the full-field was calculated by ImageJ and the number of terminal tubules per mm 2 were obtained and further analyzed statistically.

Immunohistochemistry and analysis of SPC-positive cells
Immunohistochemistry was performed on the 5 μm-thick deparaffinized sections using Imm-PRESS Polymer Detection Kit for IHC. A rabbit polyclonal antibody against a synthetic peptide fragment of human prosurfactant protein C (SPC) (within residues 1-100) was used at 1:5,000 dilution. The specificity of the antibody for AECIIs in the mouse lung has been validated [19]. The public identifier from the Antibody Registry is RRID: AB_10674024. After immunohistochemistry, SPC-positive cells were scored in random portions of a section in eight photomicrographs using a 40x objective lens. The total number of 6627 cells on average was counted per sample. The results were presented as a ratio of the number of positive cells per total number of cells [20]. This method has a limitation, in which as Fgf10 mutations have been known to affect epithelial tube formation and branching [21] and therefore the number of distal airway epithelial cells in total may be impacted, the accuracy of the data would be ensured by normalizing AECII cell number to total epithelial cell number.

Microscopy and image processing
Images for sections were collected with a Nikon DS-Fi1 camera on a Leica DM5000B microscope. Embryos, skeletons, and intestines were observed and imaged using a Leica DFC310FX camera on a Leica M165FC stereomicroscope. Image manipulation such as levels and color balance adjustments, made to some images, and assembly of figures were performed with Adobe Photoshop CS6 Extended (RRID:SCR_014199).

Reverse transcription quantitative polymerase chain reaction (RT-qPCR)
Dissected embryonic lungs (E18.5) were immersed in RNAlater at 4˚C overnight and stored at -80˚C until further use. Total RNA was extracted using NucleoZOL. Five micrograms of RNA were reverse-transcribed to cDNA with FastGene cDNA Synthesis 5x ReadyMix. RT-qPCR was performed in duplicate wells with BrightGreen 2X qPCR MasterMix-No Dye and LightCycler Nano System using primers as shown in S4

Morphometric analysis of the cecum
Wild type (n = 3), type I (n = 4), and type II (n = 2) dissected intestines containing the cecum, colons, and part of small intestines at E18.5 were photographed at the same magnification (S7 The length of each cecum was measured using ImageJ. Briefly, multi-points were set from the most proximal point abutting the small intestine and colon to the most distal cecum, approximately 16 to 18 points in the case of wild type cecum. The sum of distance between each point was regarded as the length of the cecum.

Statistical analysis
Significance was determined by one-way analysis of variance (ANOVA) (Fig 4B), Dunnett's test (Fig 5I), or two-tailed unpaired Student's t-test (S5L and S7J Figs). Correlation analysis was done by calculating correlation coefficient (Fig 6J-6L, S5M Fig). Chi-square test was used for the comparison of deep sequencing data on DNA from two or three different tissues (Table 4). Data were presented as mean ± standard error of the mean (SEM). Values of p<0.05 were considered significant otherwise stated. Data analysis was essentially performed using Microsoft Excel except for Dunnett's test, which was done using R (https://www.r-project.org) (RRID:SCR_001905).

Generation of Fgf10-CRISPR founder embryos and their limb skeletal structures
To generate Fgf10-CRISPR mice, we used a single guide RNA (sgRNA) to target exon 3, as previously described (S1 Fig) [14]. According to Hashimoto and Takemoto (2015) [15], the resultant Fgf10-CRISPR F0 embryos were classified by the limb phenotype into three types: type I, no limb; type II, limb defect; and type III, normal limbs. Since skeletal pattern of these three limb phenotypes has not been revealed, cartilage and bone staining was performed on E18.5 Fgf10-crispants (n = 4 for wild type; n = 11 for mutant embryos) ( Table 2, Fig 1, S2 Fig). Normally, the appendicular skeleton consists of girdle ("pectoral" for forelimb and "pelvic" for hindlimb) and limb elements. The skeletons of the fore-and hindlimbs consist of stylopod for humerus/femur, zeugopod for the radius/ulna and tibia/fibula, and autopod for hand/foot  [22]). In four out of six type I embryos, the size of scapula was reduced and its spine was not formed, indicating that the posterior blade was lost ( Fig 1A). As for the pelvic girdle, the ischial and pubic bones were lost, while iliac bones reduced bilaterally (Table 2 and Fig 1B). In the remaining two embryos, the scapula was intact with a spine, while the pubis and ischium were lost or reduced with the ilium intact bilaterally ( Table 2 and Fig 1C and 1D; S2C and S2D Fig). Thus, we regarded the latter phenotype as a type between I and II (type I/ II). In type II embryos (n = 2), all limb elements were lost in one or two appendages and in the remaining appendages, limb truncations were observed at the stylopod or zeugopod level (  (Fig 1P-1R). These limb skeletal patterns indicated that limb truncations were observed in the stylopodial, zeugopodial (type II), or pelvic region (type I). Formation of the clavicle, and the superior part of the scapula and pelvic girdle were not affected by lack of Fgf10 gene products.

Genomic analysis of Fgf10-CRISPR founder embryos
To prescreen mutations at the Fgf10 locus, we performed an enzyme mismatch cleavage assay (S3A and S3B Fig) and Sanger sequencing with the use of DNA from neck skin tissues of the E16.5 embryos. We estimated mutation frequencies based on band intensities [23], and found an embryo (type II, #4) with obvious limb defects but showed faint band intensities for mutated DNA after endonuclease treatment (   To get more accurate mutation frequencies, we performed deep sequencing of PCR amplicons containing the Fgf10 target site by use of embryonic (E16.5) neck DNA. By multiple amplicon analysis, we could distinguish between wild type, in-frame and frameshift mutations (Fig 2). We also deduced amino acid sequences encoded by the targeted Fgf10 locus (Fig 3). Among in-frame mutations, we found that codons that encode Lys196 and His201 were preserved or those that encode His201 alone were preserved as in #3, #5, and #9 embryos; both residues are known to be indispensable for FGF receptor dimerization [24]. We calculated the percentage of the gene sequence that still encodes putative functional FGF10 proteins with both residues in each type of mutant at E16.5, finding 8.3 ± 6.2% for type I, 25.3 ± 2.7% for type II, and 54.3 ± 9.5% for type III, n = 3 for each type, mean ± SEM) ( Table 3) (Fig 4A and  4B). In the case of E18.5 neck DNA, it was 3.4 ± 2.9% for type I (n = 3), 33% for type II (n = 1), and 66.5 ± 8.8% for type III, respectively (n = 8) (Table 1).
To know whether genotyping results using DNA from different tissues were similar, we performed deep sequencing on DNA extracted from lung and limb dermal tissues as well as neck skin tissues of the E18.5 embryos (Table 1; S2 Table). We evaluated deep sequencing results on type III embryos (n = 10) with Chi-square test (Table 4). We found that in 6  crispants. wt, wild type Fgf10 target nucleotide sequence; in, in-frame mutations by small insertion or deletion; out, frameshift mutations. Among in-frame mutations, the percentage of those that preserve the codons for Lys196 and His201 is shown in green, while the percentage that eliminate either of them or both is shown in yellow. B, percentage of wild type and in-frame mutations in which nucleotides for both Lys196 and His201 are retained: 8.3 ± 6.2% for type I, 25.3 ± 2.7% for type II, and 54.3 ± 9.5% for type III. Data are presented as means ± SEM. � p<0.01 (p = 0.008). https://doi.org/10.1371/journal.pone.0240333.g004

PLOS ONE
Lung of Fgf10-CRISPR mosaic mouse embryos there was a significant difference in the frequencies of the functional Fgf10 genotypes from the different tissues and in 4 embryos there was no significant difference in the frequencies of the functional Fgf10 genotypes from the different tissues (Table 4). In type I and type II embryos, there is a significant decrease in the number compared with that of wild type. Data are presented as means ± SEM. Source data are available in Table 3. Ac, accessory lobe; Ca, caudal lobe; Cr, cranial lobe; L, left lobe; Mi, middle lobe. Scale bars: 500 μm in (A-D), and 100 μm in (E-H). https://doi.org/10.1371/journal.pone.0240333.g005

Analysis of lung phenotypes
We next examined the lung histology of Fgf10-mosaic F0 embryos on E16.5 (Table 3, Fig 5, S4  Fig). In the normal mouse lung, there are five lung lobes, four on the right (cranial, middle, caudal, and accessory) and one on the left (Fig 5A). The lung histology of type I embryos demonstrated either the absence of all lobes (two of three embryos) or only the presence of a residual lobe (one of three) (Fig 5B). We found that the accessory lobe was absent in two out of three type II embryos, and that either middle or caudal lobe was additionally absent or hypoplastic (Fig 5C). In type III embryos, all five lung lobes were present (Fig 5D). In E18.5 type II embryos, the accessory lung lobe was absent from all three embryos examined (embryo #11_2, 24_2, 25_2; S5I-S5K Fig).
Mouse lung development is divided into five stages, namely embryonic, pseudoglandular (E9.5 to 16.5), canalicular (E16.5 to 17.5), saccular (E17.5 to P5), and alveolar (P5 to P28) [25]. At E16.5, after branching morphogenesis the lung is in the transition of pseudoglandular to canalicular stage. We found that the number of terminal tubules was significantly decreased in type I and II residual lungs (Fig 5F, 5G and 5I), but not significantly altered in type III lungs (Fig 5H and 5I), compared with wild type (Fig 5E and 5I). By E18.5, the lung enters the saccular stage and cuboidal cells that produce surfactant protein C (SPC), i.e. alveolar type 2 epithelial (AECII) cells, differentiate on the wall of the developing airway space. Since AECII cell differentiation is perturbed in the Fgf10-heterozygous lung [26,27], we performed Fgf10 genotypes in limb (J) or neck DNA (K). L, the number of SPC-positive cells in type III embryos is not to correlated to their weight (K). Source data for (J-L) are available in Table 1. M, Schematic drawings to show percentage of total reads for the type III Fgf10 crispants. wt, wild type Fgf10 target nucleotide sequence; in, in-frame mutations by small insertion or deletion; out, out-of-frame mutations. Among in-frame mutations, the percentage of those that preserve the codons for Lys196 and His201 is shown in green, while the percentage that eliminate either of them is shown in yellow.
https://doi.org/10.1371/journal.pone.0240333.g006  immunohistochemistry on E18.5 lung sections to detect AECII. We counted the number of SPC-positive cells in lungs of type III embryos and examined whether it was correlated to the percentage of functional Fgf10 genotype. Representative data of immunostaining for type III lungs are shown (Fig 6B and 6E) along with those for type II lungs (Fig 6A and 6D) and positive/negative controls (Fig 6C, 6F and 6G-6I). We found that the number of SPC-positive cells in lungs of type III embryos was correlated to the percentage of functional Fgf10 genotype ( Fig  6J, correlation coefficient is 0.909 for limb DNA; Fig 6K, 0.924 for neck DNA). Since the delay in development can be caused by variation between embryos even from a same litter due to crowding in the uterus, we also examined whether the number of SPC-positive cells was correlated to the weight of each embryo, as an index for developmental stage (Table 1). It seemed that there was no correlation between them (Fig 6L; correlation coefficient is 0.301).
On the other hand, since it was reported that the mesenchymal FGF10 is also important for proper mesenchymal lineage formation during lung development as revealed by Fgf10-hypomorphic lungs [20,28], we examined expression of mesenchymal and its related marker genes [29] in type III lungs at E18.5 by quantitative PCR (Table 1, S4 Table). However, we could not detect a significant decrease in expression of these genes (S6 Fig), possibly because the percentage for functional Fgf10 genotypes was too high to mimic a Fgf10-hypomorphic phenotype in lung mesenchymal lineage formation (61.19% in #11, 65.98% in #12, and 70.11% in #17 embryos). This expression analysis should be further clarified on earlier stage lungs, and in those which have more loss-of-function mutations in the Fgf10 gene.

Discussion
Fgf10 genome edited F0 mice are classified into three limb phenotypes: no limb (type I), limb defect (type II), and normal limbs (type III). In this study, we first examined appendicular skeletal patterns of the Fgf10 mosaic mutant mice by cartilage and bone staining. We found that varying degrees of limb and girdle bone truncations were observed in type I and type II embryos. We next estimated the rate of putative functional Fgf10 genotypes in each mutant type and the mean percentage for type II embryos is 25.3 ± 2.7 by deep sequencing of neck DNA. Comparison of deep sequencing results on DNA from different tissues (neck, lung, and limb dermis) revealed that in half of the embryos the rate of putative functional Fgf10 genotypes were significantly varied between the three DNA. Lung phenotypes were examined on two embryonic stages, E16.5 and E18.5. In type II embryos, the accessory lobe was lost and the number of terminal tubules was significantly reduced. In type II and "normal limb (type III)" embryos, the number of alveolar type 2 cells, immunostained by anti-surfactant C protein, was correlated with the rate of putative functional Fgf10 genotypes.

Truncation patterns of girdle and limb bones in Fgf10-mosaic mutants
Skeletal staining in this study revealed that in type II embryos, limb bone structures are truncated in more distal structures rather than forming a miniature of whole limb and girdle structures. Also, in the mutants, the zeugopodial element consists of one ulna-or tibia-like bone: This pattern is different from those after extirpation of apical ectodermal ridge in chicken limb buds, where both zeugopod bones are truncated [30]. Although the number of type II embryos examined in this study is only two and therefore further studies are needed, a reduced number of Fgf10-expressing cells in the mosaic mutants does not seem to know the morphological pattern to be formed. Looking the expression of zeugopod patterning genes could help to elucidate the underlying mechanisms.
Regarding girdle bones, it has been known that the posterior (inferior) half of the scapula is absent and most of the pelvic structures is lacking in Fgf10-null embryos [12,13]. This study showed that in type I/II Fgf10-mosaic mutants the posterior scapula with the spine retained. In the mouse, while the pelvic girdle and most of the scapular blade derives solely from the lateral plate mesoderm (LPM), the medial edge of the blade arises from dermomyotome cells and the neural crest cells contribute to the scapular spine [31][32][33]. Thus, it is likely that a reduced number of Fgf10-expressing cells in the LPM and neural crest cells [34] contributed to the formation of the posterior (inferior) scapula in type I/II embryos. In the case of pelvic girdle, a reduced number of Fgf10-expressing cells in the LPM contributed to the formation of the inferior pelvic structures in type I/II embryos.

Fgf10 mutation rate may not be identical in different tissues
Here we showed that there is a correlation between limb phenotypes, lung phenotypes, and overall defective mutation rate in Fgf10-CRISPR knockout F0 mice. It has been shown that even if CRISPR-genome editing materials are introduced into fertilized eggs at the one-cell stage, mosaicism in the target genotype is caused under certain experimental conditions until the target sequences are mutated [6,7]. During the early embryonic stage, mutated cells and non-mutated cells are intermingled and distributed to lung and limb primordia proportionally [35], which we believe contributes to the correlation between limb and lung phenotypes. However, statistical analysis on deep sequencing data of the genotypes shows in more than half cases there is a significant difference of the defective mutation rate in DNA from different tissues. In this study, we could not identify the reason why there were differences in the rate of the functional Fgf10 genotypes. In some developmental settings, Fgf10-expressing normal cells might have a growth advantage than Fgf10-deficient cells, which should be further explored in the next studies [10].

Fgf10-CRISPR F0 mice possess degrees of lung dysgenesis that correlate with putative functional FGF10 dosage
Accessory lobe formation was firstly impaired due to reduction in putative functional FGF10 dosage. Ramasamy et al. (2007) reported that Fgf10-hypomorphic mice, in which Fgf10 expression is reduced by 27%, compared with Fgf10-heterozygous mutants, lack accessory lobe formation [20]. Since Fgf10 expression is highest in the mesenchyme of the accessory lobe on E11.5 and persists as such until at least E18.5 [21,29], high Fgf10 dosage was believed to be required for accessory lobe formation. Our result shows that more than 25.3 ± 2.7% (from the analysis of E16.5 embryos) of functional Fgf10 gene product is required for accessory lobe formation. Since reduction or loss of accessory lobe formation has been observed in loss-of-function mutants of Gli3, Fog2, and Gata4 [36][37][38], FGF10 signaling may be related to the expression of these transcription factors during accessory lobe formation.
Lineage-tracing analysis of Fgf10-expressing cells during mouse lung development has revealed that there are two waves of Fgf10 expression; the first begins from E11.5 and the second from E15.5 [29]. We think that a reduction in the latter cells would cause the decrease in the number of terminal tubules in type I and type II Fgf10-CRISPR F0 embryos. A recent study has shown that AECII as well as AECI independent progenitors are present at E13.5 as revealed by single cell RNA sequencing analysis [39]. By the saccular stage, the lung mesenchyme surrounding the epithelium becomes thinner and cuboidal AECII cells, characterized by the production of surfactant proteins, begin to differentiate (summarized in [40]). The number of AECII cells is found to be lower in the Fgf10-heterozygous lung [26,27], which is similar to our result in type III Fgf10-mosaic mice. Taken together, Fgf10 expression levels regulate the number of alveolar type 2 epithelial cells in mouse lung in somewhat dose-dependent manner.
Regarding FGF10 dosage, the Fgf10-hypomorphic mouse aforementioned [20,28] has also clarified that higher FGF10 dosage is required for full development of colonic crypts [41]. Our preliminary observation (S7 Fig) showed that in type I embryos the length of the cecum was significantly shorter than the wild type and there was an atresia of the colon as reported in the Fgf10-null mutants [41,42]. The length of the blind colon varied depending upon the embryos (S7D-S7G Fig), suggesting a correlation with their FGF10 dosage.
Fgf10-CRISPR F0 mice can serve as a series of lung disease models Volckaert et al. (2013) reported that lung branching morphogenesis does not require localized Fgf10 expression in the distal mesenchyme, because ubiquitous Fgf10 overexpression can induce lung formation in Fgf10 knockout mice [43]. Furthermore, Fgf10 expression after lung initiation is required for branching and proximal-distal differentiation by regulating Sox2/9 expression in the epithelium [43]. Our result also supports this property of Fgf10 as the number of terminal tubules decreased in the Fgf10-mosaic mutants where functional FGF10 dosage is reduced.
Recent papers have reported that mutations and single nucleotide polymorphisms (SNPs) in the Fgf10 gene are correlated with human lung disease (reviewed in [44]). For example, the absence of the right medial-basal airway is associated with a type of chronic obstructive pulmonary diseases and two types of SNPs are found within the same intron of Fgf10 in those cohorts [45]. Although some molecular and cellular differences have been identified between mouse and human lungs (reviewed in [46]), the F0 mice generated by CRISPR/Cas9-mediated Fgf10 gene editing can become a model animal to study the pathophysiology of human pulmonary hypoplasia and related chronic lung diseases that may be rooted in the developmental stage as recently postulated [47]).
Supporting information S1  Table 4) are shown for reference. , and pelvic girdle (K). au, autopod; f, femur; fi, fibula; h, humerus; il, ilium; is, ischium; p, pubis; r, radius; sc, scapula; sp, spine; st, stylopod; t, tibia; u, ulna; ze, zeugopod. Scale bars: 1 mm. (TIF) S3 Fig. Mismatch cleavage assay. A, PCR primers were designed in the upstream region of the exon 3 and in the exon 3, giving rise to a PCR amplicon size of 501 bp. DNA fragments of 309 bp and 192 bp are generated by the Resolvase when the Fgf10 genome has been cleaved by Cas9 and non-homologous end joining has been achieved. B, Electrophoresis of the enzymetreated mouse genomic DNA from the Fgf10-CRISPR F0 embryonic necks. The DNA ladder for DNA size reference and a result of DNA from a wild type (WT) mouse are shown on the left. Three DNA fragments of approximately 500 bp (▼), 300 bp (5 in gray), and 200 bp (5) are seen in all the lanes except for the wild type and #4 lanes. In embryo #3, an extra band for large insertion (328 base) is shown (see Fig 2). C, Genomic analysis of the #4 embryo as  (c), colon, and small intestine (si) are shown from three embryos examined. Ileum and colon were cut at dissection. D-G, type I cecum is reduced compared with the wild type. Whether the cecum epithelium is absent or not cannot be identified from these photos. Type I embryos show an atresia of the colon, but the length varies depending upon the embryos. H-I, type II embryos examined (n = 2) do not exhibit a reduced cecum or an atresia of the colon. The colons presented here were cut as distally as possible. J, the approximate length of the cecum. The length of type I cecum is significantly decreased compared with the wild type. The length of two type II embryos examined is also shown for reference.