Validated WGS and WES protocols proved saliva-derived gDNA as an equivalent to blood-derived gDNA for clinical and population genomic analyses.

BACKGROUND: Whole exome sequencing (WES) and whole genome sequencing (WGS) have become standard methods in human clinical diagnostics as well as in population genomics (POPGEN). Blood-derived genomic DNA (gDNA) is routinely used in the clinical environment. Conversely, many POPGEN studies and commercial tests benefit from easy saliva sampling. Here, we evaluated the quality of variant call sets and the level of genotype concordance of single nucleotide variants (SNVs) and small insertions and deletions (indels) for WES and WGS using paired blood- and saliva-derived gDNA isolates employing genomic reference-based validated protocols. METHODS: The genomic reference standard Coriell NA12878 was repeatedly analyzed using optimized WES and WGS protocols, and data calls were compared with the truth dataset published by the Genome in a Bottle Consortium. gDNA was extracted from the paired blood and saliva samples of 10 participants and processed using the same protocols. A comparison of paired blood-saliva call sets was performed in the context of WGS and WES genomic reference-based technical validation results. RESULTS: The quality pattern of called variants obtained from genomic-reference-based technical replicates correlates with data calls of paired blood-saliva-derived samples in all levels of tested examinations despite a higher rate of non-human contamination found in the saliva samples. The F1 score of 10 blood-to-saliva-derived comparisons ranged between 0.8030-0.9998 for SNVs and between 0.8883-0.9991 for small-indels in the case of the WGS protocol, and between 0.8643-0.999 for SNVs and between 0.7781-1.000 for small-indels in the case of the WES protocol. CONCLUSION: Saliva may be considered an equivalent material to blood for genetic analysis for both WGS and WES under strict protocol conditions. The accuracy of sequencing metrics and variant-detection accuracy is not affected by choosing saliva as the gDNA source instead of blood but much more significantly by the genomic context, variant types, and the sequencing technology used.

Ovol2 promoter mutations in mice and human illuminate species-specific phenotypic divergence.

Pathogenic variants in the highly conserved OVOL2 promoter region cause posterior polymorphous corneal dystrophy (PPCD) 1 by inducing an ectopic expression of the endothelial OVOL2 mRNA. Here we produced an allelic series of Ovol2 promoter mutations in the mouse model including the heterozygous c.-307T>C variant (RefSeq NM_021220.4) causing PPCD1 in humans. Despite the high evolutionary conservation of the Ovol2 promoter, only some alterations of its sequence had phenotypic consequences in mice. Four independent sequence variants in the distal part of the Ovol2 promoter had no significant effect on endothelial Ovol2 mRNA level or caused any ocular phenotype. In contrast, the mutation c.-307T>C resulted in increased Ovol2 expression in the corneal endothelium. However, only a small fraction of adult mice c.-307T>C heterozygotes developed ocular phenotypes such as irido-corneal adhesions, and corneal opacity. Interestingly, phenotypic penetrance was increased at embryonic stages. Notably, c.-307T>C mutation is located next to the Ovol1/Ovol2 transcription factor binding site. Mice carrying an allele with a deletion encompassing the Ovol2 binding site c.-307_-320del showed significant Ovol2 gene upregulation in the cornea endothelium and exhibited phenotypes similar to the c.-307T>C mutation. In conclusion, although the mutations c.-307T>C and -307_-320del lead to a comparably strong increase in endothelial Ovol2 expression as seen in PPCD1 patients, endothelial dystrophy was not observed in the mouse model, implicating species-specific differences in endothelial cell biology. Nonetheless, the emergence of dominant ocular phenotypes associated with Ovol2 promoter variants in mice implies a potential role of this gene in eye development and disease.

Conserved enhancers control notochord expression of vertebrate Brachyury.

The cell type-specific expression of key transcription factors is central to development and disease. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three conserved Brachyury-controlling notochord enhancers, T3, C, and I, in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, in cis deletion of all three enhancers in mouse abolishes Brachyury/T/Tbxt expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. The three Brachyury-driving notochord enhancers are conserved beyond mammals in the brachyury/tbxtb loci of fishes, dating their origin to the last common ancestor of jawed vertebrates. Our data define the vertebrate enhancers for Brachyury/T/TBXTB notochord expression through an auto-regulatory mechanism that conveys robustness and adaptability as ancient basis for axis development.

Conserved enhancer logic controls the notochord expression of vertebrate Brachyury.

The cell type-specific expression of key transcription factors is central to development. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three Brachyury-controlling notochord enhancers T3, C, and I in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, deletion of all three enhancers in mouse abolishes Brachyury/T expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. Sequence and functional conservation of Brachyury-driving notochord enhancers with the brachyury/tbxtb loci from diverse lineages of fishes dates their origin to the last common ancestor of jawed vertebrates. Our data define the enhancers for Brachyury/T/TBXTB notochord expression as ancient mechanism in axis development.

Chromatin Remodeling Enzyme Snf2h Is Essential for Retinal Cell Proliferation and Photoreceptor Maintenance.

Chromatin remodeling complexes are required for many distinct nuclear processes such as transcription, DNA replication, and DNA repair. However, the contribution of these complexes to the development of complex tissues within an organism is poorly characterized. Imitation switch (ISWI) proteins are among the most evolutionarily conserved ATP-dependent chromatin remodeling factors and are represented by yeast Isw1/Isw2, and their vertebrate counterparts Snf2h (Smarca5) and Snf2l (Smarca1). In this study, we focused on the role of the Snf2h gene during the development of the mammalian retina. We show that Snf2h is expressed in both retinal progenitors and post-mitotic retinal cells. Using Snf2h conditional knockout mice (Snf2h cKO), we found that when Snf2h is deleted, the laminar structure of the adult retina is not retained, the overall thickness of the retina is significantly reduced compared with controls, and the outer nuclear layer (ONL) is completely missing. The depletion of Snf2h did not influence the ability of retinal progenitors to generate all the differentiated retinal cell types. Instead, the Snf2h function is critical for the proliferation of retinal progenitor cells. Cells lacking Snf2h have a defective S-phase, leading to the entire cell division process impairments. Although all retinal cell types appear to be specified in the absence of the Snf2h function, cell-cycle defects and concomitantly increased apoptosis in Snf2h cKO result in abnormal retina lamination, complete destruction of the photoreceptor layer, and consequently, a physiologically non-functional retina.

Chromatin remodeling enzyme Snf2h is essential for retinal cell proliferation and photoreceptor maintenance.

Chromatin remodeling complexes are required for many distinct nuclear processes such as transcription, DNA replication and DNA repair. However, how these complexes contribute to the development of complex tissues within an organism is poorly characterized. Imitation switch (ISWI) proteins are among the most evolutionarily conserved ATP-dependent chromatin remodeling factors and are represented by yeast Isw1/Isw2, and their vertebrate counterparts Snf2h (Smarca5) and Snf2l (Smarca1). In this study, we focused on the role of the Snf2h gene during development of the mammalian retina. We show that Snf2h is expressed in both retinal progenitors and post-mitotic retinal cells. Using Snf2h conditional knockout mice ( Snf2h cKO), we found that when Snf2h is deleted the laminar structure of the adult retina is not retained, the overall thickness of the retina is significantly reduced compared with controls, and the outer nuclear layer (ONL) is completely missing. Depletion of Snf2h did not influence the ability of retinal progenitors to generate all of the differentiated retinal cell types. Instead, Snf2h function is critical for proliferation of retinal progenitor cells. Cells lacking Snf2h have a defective S-phase, leading to the entire cell division process impairments. Although, all retinal cell types appear to be specified in the absence of Snf2h function, cell cycle defects and concomitantly increased apoptosis in Snf2h cKO result in abnormal retina lamination, complete destruction of the photoreceptor layer and, consequently, in a physiologically non-functional retina.

Identification of Nodal-dependent enhancer of amphioxus Chordin sufficient to drive gene expression into the chordate dorsal organizer.

The core molecular mechanisms of dorsal organizer formation during gastrulation are highly conserved within the chordate lineage. One of the key characteristics is that Nodal signaling is required for the organizer-specific gene expression. This feature appears to be ancestral, as evidenced by the presence in the most basally divergent chordate amphioxus. To provide a better understanding of the evolution of organizer-specific gene regulation in chordates, we analyzed the cis-regulatory sequence of amphioxus Chordin in the context of the vertebrate embryo. First, we generated stable zebrafish transgenic lines, and by using light-sheet fluorescent microscopy, characterized in detail the expression pattern of GFP driven by the cis-regulatory sequences of amphioxus Chordin. Next, we performed a 5'deletion analysis and identified an enhancer sufficient to drive the expression of the reporter gene into a chordate dorsal organizer. Finally, we found that the identified enhancer element strongly depends on Nodal signaling, which is consistent with the well-established role of this pathway in the regulation of the expression of dorsal organizer-specific genes across chordates. The enhancer identified in our study may represent a suitable simple system to study the interplay of the evolutionarily conserved regulatory mechanisms operating during early chordate development.

Multiple roles of Pax6 in postnatal cornea development.

Mammalian corneal development is a multistep process, including formation of the corneal epithelium (CE), endothelium and stroma during embryogenesis, followed by postnatal stratification of the epithelial layers and continuous renewal of the epithelium to replace the outermost corneal cells. Here, we employed the Cre-loxP system to conditionally deplete Pax6 proteins in two domains of ocular cells, i.e., the ocular surface epithelium (cornea, limbus and conjunctiva) (OSE) or postnatal CE via K14-cre or Aldh3-cre, respectively. Earlier and broader inactivation of Pax6 in the OSE resulted in thickened OSE with CE and limbal cells adopting the conjunctival keratin expression pattern. More restricted depletion of Pax6 in postnatal CE resulted in an abnormal cornea marked by reduced epithelial thickness despite increased epithelial cell proliferation. Immunofluorescence studies revealed loss of intermediate filament Cytokeratin 12 and diffused expression of adherens junction components, together with reduced tight junction protein, Zonula occludens-1. Furthermore, the expression of Cytokeratin 14, a basal cell marker in apical layers, indicates impaired differentiation of CE cells. Collectively, our data demonstrate that Pax6 is essential for maintaining proper differentiation and strong intercellular adhesion in postnatal CE cells, whereas limbal Pax6 is required to prevent the outgrowth of conjunctival cells to the cornea.

Conditional knockout of hephaestin in the neural retina disrupts retinal iron homeostasis.

Iron accumulation has been implicated in degenerative retinal diseases. It can catalyze the production of damaging reactive oxygen species. Previous work has demonstrated iron accumulation in multiple retinal diseases, including age-related macular degeneration and diabetic retinopathy. In mice, systemic knockout of the ferroxidases ceruloplasmin (Cp) and hephaestin (Heph), which oxidize iron, results in retinal iron accumulation and iron-induced degeneration. To determine the role of Heph in the retina, we generated a neural retina-specific Heph knockout on a background of systemic Cp knockout. This resulted in elevated neural retina iron. Conversely, retinal ganglion cells had elevated transferrin receptor and decreased ferritin, suggesting diminished iron levels. The retinal degeneration observed in systemic Cp-/-, Heph-/- mice did not occur. These findings indicate that Heph has a local role in regulating neural retina iron homeostasis, but also suggest that preserved Heph function in either the RPE or systemically mitigates the degeneration phenotype observed in the systemic Cp-/-, Heph-/- mice.

Minimal effect of conditional ferroportin KO in the neural retina implicates ferrous iron in retinal iron overload and degeneration.

Iron-induced oxidative stress can cause or exacerbate retinal degenerative diseases. Retinal iron overload has been reported in several mouse disease models with systemic or neural retina-specific knockout (KO) of homologous ferroxidases ceruloplasmin (Cp) and hephaestin (Heph). Cp and Heph can potentiate ferroportin (Fpn) mediated cellular iron export. Here, we used retina-specific Fpn KO mice to test the hypothesis that retinal iron overload in Cp/Heph DKO mice is caused by impaired iron export from neurons and glia. Surprisingly, there was no indication of retinal iron overload in retina-specific Fpn KO mice: the mRNA levels of transferrin receptor in the retina were not altered at 7-10-months age. Consistent with this, levels and localization of ferritin light chain were unchanged. To "stress the system", we injected iron intraperitoneally into Fpn KO mice with or without Cp KO. Only mice with both retina-specific Fpn KO and Cp KO had modestly elevated retinal iron levels. These results suggest that impaired iron export through Fpn is not sufficient to explain the retinal iron overload in Cp/Heph DKO mice. An increase in the levels of retinal ferrous iron caused by the absence of these ferroxidases, followed by uptake into cells by ferrous iron importers, is most likely necessary.

mTORC1-induced retinal progenitor cell overproliferation leads to accelerated mitotic aging and degeneration of descendent Müller glia.

Retinal progenitor cells (RPCs) divide in limited numbers to generate the cells comprising vertebrate retina. The molecular mechanism that leads RPC to the division limit, however, remains elusive. Here, we find that the hyperactivation of mechanistic target of rapamycin complex 1 (mTORC1) in an RPC subset by deletion of tuberous sclerosis complex 1 (Tsc1) makes the RPCs arrive at the division limit precociously and produce Müller glia (MG) that degenerate from senescence-associated cell death. We further show the hyperproliferation of Tsc1-deficient RPCs and the degeneration of MG in the mouse retina disappear by concomitant deletion of hypoxia-induced factor 1-alpha (Hif1a), which induces glycolytic gene expression to support mTORC1-induced RPC proliferation. Collectively, our results suggest that, by having mTORC1 constitutively active, an RPC divides and exhausts mitotic capacity faster than neighboring RPCs, and thus produces retinal cells that degenerate with aging-related changes.

Optical Clearing and Light Sheet Microscopy Imaging of Amphioxus.

Cephalochordates (amphioxi or lancelets) are representatives of the most basally divergent group of the chordate phylum. Studies of amphioxus development and anatomy hence provide a key insight into vertebrate evolution. More widespread use of amphioxus in the evo-devo field would be greatly facilitated by expanding the methodological toolbox available in this model system. For example, evo-devo research on amphioxus requires deep understanding of animal anatomy. Although conventional confocal microscopy can visualize transparent amphioxus embryos and early larvae, the imaging of later developmental stages is problematic because of the size and opaqueness of the animal. Here, we show that light sheet microscopy combined with tissue clearing methods enables exploration of large amphioxus specimens while keeping the surface and the internal structures intact. We took advantage of the phenomenon of autofluorescence of amphioxus larva to highlight anatomical details. In order to investigate molecular markers at the single-cell level, we performed antibody-based immunodetection of melanopsin and acetylated-α-tubulin to label rhabdomeric photoreceptors and the neuronal scaffold. Our approach that combines light sheet microscopy with the clearing protocol, autofluorescence properties of amphioxus, and antibody immunodetection allows visualizing anatomical structures and even individual cells in the 3D space of the entire animal body.

MEIS-WNT5A axis regulates development of fourth ventricle choroid plexus.

The choroid plexus (ChP) produces cerebrospinal fluid and forms an essential brain barrier. ChP tissues form in each brain ventricle, each one adopting a distinct shape, but remarkably little is known about the mechanisms underlying ChP development. Here, we show that epithelial WNT5A is crucial for determining fourth ventricle (4V) ChP morphogenesis and size in mouse. Systemic Wnt5a knockout, or forced Wnt5a overexpression beginning at embryonic day 10.5, profoundly reduced ChP size and development. However, Wnt5a expression was enriched in Foxj1-positive epithelial cells of 4V ChP plexus, and its conditional deletion in these cells affected the branched, villous morphology of the 4V ChP. We found that WNT5A was enriched in epithelial cells localized to the distal tips of 4V ChP villi, where WNT5A acted locally to activate non-canonical WNT signaling via ROR1 and ROR2 receptors. During 4V ChP development, MEIS1 bound to the proximal Wnt5a promoter, and gain- and loss-of-function approaches demonstrated that MEIS1 regulated Wnt5a expression. Collectively, our findings demonstrate a dual function of WNT5A in ChP development and identify MEIS transcription factors as upstream regulators of Wnt5a in the 4V ChP epithelium.

Wnt/ß-Catenin Signaling Promotes Differentiation of Ischemia-Activated Adult Neural Stem/Progenitor Cells to Neuronal Precursors.

Modulating endogenous regenerative processes may represent a suitable treatment for central nervous system (CNS) injuries, such as stroke or trauma. Neural stem/progenitor cells (NS/PCs), which naturally reside in the subventricular zone (SVZ) of the adult brain, proliferate and differentiate to other cell types, and therefore may compensate the negative consequences of ischemic injury. The fate of NS/PCs in the developing brain is largely influenced by Wingless/Integrated (Wnt) signaling; however, its role in the differentiation of adult NS/PCs under ischemic conditions is still enigmatic. In our previous study, we identified the Wnt/ß-catenin signaling pathway as a factor promoting neurogenesis at the expense of gliogenesis in neonatal mice. In this study, we used adult transgenic mice in order to assess the impact of the canonical Wnt pathway modulation (inhibition or hyper-activation) on NS/PCs derived from the SVZ, and combined it with the middle cerebral artery occlusion (MCAO) to disclose the effect of focal cerebral ischemia (FCI). Based on the electrophysiological properties of cultured cells, we first identified three cell types that represented in vitro differentiated NS/PCs - astrocytes, neuron-like cells, and precursor cells. Following FCI, we detected fewer neuron-like cells after Wnt signaling inhibition. Furthermore, the immunohistochemical analysis revealed an overall higher expression of cell-type-specific proteins after FCI, indicating increased proliferation and differentiation rates of NS/PCs in the SVZ. Remarkably, Wnt signaling hyper-activation increased the abundance of proliferating and neuron-like cells, while Wnt pathway inhibition had the opposite effect. Finally, the expression profiling at the single cell level revealed an increased proportion of neural stem cells and neuroblasts after FCI. These observations indicate that Wnt signaling enhances NS/PCs-based regeneration in the adult mouse brain following FCI, and supports neuronal differentiation in the SVZ.

Meis homeobox genes control progenitor competence in the retina.

The vertebrate eye is derived from the neuroepithelium, surface ectoderm, and extracellular mesenchyme. The neuroepithelium forms an optic cup in which the spatial separation of three domains is established, namely, the region of multipotent retinal progenitor cells (RPCs), the ciliary margin zone (CMZ)-which possesses both a neurogenic and nonneurogenic potential-and the optic disk (OD), the interface between the optic stalk and the neuroretina. Here, we show by genetic ablation in the developing optic cup that Meis1 and Meis2 homeobox genes function redundantly to maintain the retinal progenitor pool while they simultaneously suppress the expression of genes characteristic of CMZ and OD fates. Furthermore, we demonstrate that Meis transcription factors bind regulatory regions of RPC-, CMZ-, and OD-specific genes, thus providing a mechanistic insight into the Meis-dependent gene regulatory network. Our work uncovers the essential role of Meis1 and Meis2 as regulators of cell fate competence, which organize spatial territories in the vertebrate eye.

Identification of retinal ganglion cell types and brain nuclei expressing the transcription factor Brn3c/Pou4f3 using a Cre recombinase knock-in allele.

Members of the POU4F/Brn3 transcription factor family have an established role in the development of retinal ganglion cell (RGCs) types, the main transducers of visual information from the mammalian eye to the brain. Our previous work using sparse random recombination of a conditional knock-in reporter allele expressing alkaline phosphatase (AP) and intersectional genetics had identified three types of Brn3c positive (Brn3c+ ) RGCs. Here, we describe a novel Brn3cCre mouse allele generated by serial Dre to Cre recombination and use it to explore the expression overlap of Brn3c with Brn3a and Brn3b and the dendritic arbor morphologies and visual stimulus response properties of Brn3c+ RGC types. Furthermore, we explore brain nuclei that express Brn3c or receive input from Brn3c+ neurons. Our analysis reveals a much larger number of Brn3c+ RGCs and more diverse set of RGC types than previously reported. Most RGCs expressing Brn3c during development are still Brn3c positive in the adult, and all express Brn3a while only about half express Brn3b. Genetic Brn3c-Brn3b intersection reveals an area of increased RGC density, extending from dorsotemporal to ventrolateral across the retina and overlapping with the mouse binocular field of view. In addition, we report a Brn3c+ RGC projection to the thalamic reticular nucleus, a visual nucleus that was not previously shown to receive retinal input. Furthermore, Brn3c+ neurons highlight a previously unknown subdivision of the deep mesencephalic nucleus. Thus, our newly generated allele provides novel biological insights into RGC type classification, brain connectivity, and cytoarchitectonic.

Molecular Fingerprint of Amphioxus Frontal Eye Illuminates the Evolution of Homologous Cell Types in the Chordate Retina.

The evolution of the vertebrate eye remains so far unresolved. Amphioxus frontal eye pigment cells and photoreceptors were proposed to be homologous to vertebrate photoreceptors and retinal pigmented epithelium, based on ultrastructural morphology and gene expression analysis in B. floridae. Here, we present comparative molecular data using two additional amphioxus species, a closely related B. lanceolatum, and the most divergent A. lucayanum. Taking advantage of a unique set of specific antibodies we characterized photoreceptors and putative interneurons of the frontal eye and investigated its neuronal circuitry. Our results corroborate generally conserved molecular fingerprint among cephalochordate species. Furthermore, we performed pharmacological perturbations and found that the Notch signaling pathway, a key regulator of retina development in vertebrates, is required for correct ratios among frontal eye cell types. In summary, our study provides a valuable insight into cell-type relationships in chordate visual organs and strengthens the previously proposed homology between amphioxus frontal eye and vertebrate eyes.

Dose-dependent regulation of horizontal cell fate by Onecut family of transcription factors.

Genome duplication leads to an emergence of gene paralogs that are essentially free to undergo the process of neofunctionalization, subfunctionalization or degeneration (gene loss). Onecut1 (Oc1) and Onecut2 (Oc2) transcription factors, encoded by paralogous genes in mammals, are expressed in precursors of horizontal cells (HCs), retinal ganglion cells and cone photoreceptors. Previous studies have shown that ablation of either Oc1 or Oc2 gene in the mouse retina results in a decreased number of HCs, while simultaneous deletion of Oc1 and Oc2 leads to a complete loss of HCs. Here we study the genetic redundancy between Oc1 and Oc2 paralogs and focus on how the dose of Onecut transcription factors influences abundance of individual retinal cell types and overall retina physiology. Our data show that reducing the number of functional Oc alleles in the developing retina leads to a gradual decrease in the number of HCs, progressive thinning of the outer plexiform layer and diminished electrophysiology responses. Taken together, these observations indicate that in the context of HC population, the alleles of Oc1/Oc2 paralogous genes are mutually interchangeable, function additively to support proper retinal function and their molecular evolution does not follow one of the typical routes after gene duplication.

Neural crest cells require Meis2 for patterning the mandibular arch via the Sonic hedgehog pathway.

Cranial neural crest cells (cNCCs) originate in the anterior neural tube and populate pharyngeal arches in which they contribute to formation of bone and cartilage. This cell population also provides molecular signals for the development of tissues of non-neural crest origin, such as the tongue muscles, teeth enamel or gland epithelium. Here we show that the transcription factor Meis2 is expressed in the oral region of the first pharyngeal arch (PA1) and later in the tongue primordium. Conditional inactivation of Meis2 in cNCCs resulted in loss of Sonic hedgehog signalling in the oropharyngeal epithelium and impaired patterning of PA1 along the lateral-medial and oral-aboral axis. Failure of molecular specification of PA1, illustrated by altered expression of Hand1/2, Dlx5, Barx1, Gsc and other markers, led to hypoplastic tongue and ectopic ossification of the mandible. Meis2-mutant mice thus display craniofacial defects that are reminiscent of several human syndromes and patients with mutations in the Meis2 gene.