Discovery of genes required for body axis and limb formation by global identification of retinoic acid-regulated epigenetic marks.

Identification of target genes that mediate required functions downstream of transcription factors is hampered by the large number of genes whose expression changes when the factor is removed from a specific tissue and the numerous binding sites for the factor in the genome. Retinoic acid (RA) regulates transcription via RA receptors bound to RA response elements (RAREs) of which there are thousands in vertebrate genomes. Here, we combined chromatin immunoprecipitation sequencing (ChIP-seq) for epigenetic marks and RNA-seq on trunk tissue from wild-type and Aldh1a2-/- embryos lacking RA synthesis that exhibit body axis and forelimb defects. We identified a relatively small number of genes with altered expression when RA is missing that also have nearby RA-regulated deposition of histone H3 K27 acetylation (H3K27ac) (gene activation mark) or histone H3 K27 trimethylation (H3K27me3) (gene repression mark) associated with conserved RAREs, suggesting these genes function downstream of RA. RA-regulated epigenetic marks were identified near RA target genes already known to be required for body axis and limb formation, thus validating our approach; plus, many other candidate RA target genes were found. Nuclear receptor 2f1 (Nr2f1) and nuclear receptor 2f2 (Nr2f2) in addition to Meis homeobox 1 (Meis1) and Meis homeobox 2 (Meis2) gene family members were identified by our approach, and double knockouts of each family demonstrated previously unknown requirements for body axis and/or limb formation. A similar epigenetic approach can be used to determine the target genes for any transcriptional regulator for which a knockout is available.

Mouse but not zebrafish requires retinoic acid for control of neuromesodermal progenitors and body axis extension.

In mouse, retinoic acid (RA) is required for the early phase of body axis extension controlled by a population of neuromesodermal progenitors (NMPs) in the trunk called expanding-NMPs, but not for the later phase of body axis extension controlled by a population of NMPs in the tail called depleting-NMPs. Recent observations suggest that zebrafish utilize depleting-NMPs but not expanding-NMPs for body axis extension. In zebrafish, a role for RA in body axis extension was not supported by previous studies on aldh1a2 (raldh2) mutants lacking RA synthesis. Here, by treating zebrafish embryos with an RA synthesis inhibitor, we also found that body axis extension and somitogenesis was not perturbed, although loss of pectoral fin and cardiac edema were observed consistent with previous studies. The conclusion that zebrafish diverges from mouse in not requiring RA for body axis extension is consistent with zebrafish lacking early expanding-NMPs to generate the trunk. We suggest that RA control of body axis extension was added to higher vertebrates during evolution of expanding-NMPs.

Mutations in MYT1, encoding the myelin transcription factor 1, are a rare cause of OAVS.

BACKGROUND: Oculo-auriculo-vertebral spectrum (OAVS) is a developmental disorder involving first and second branchial arches derivatives, mainly characterised by asymmetric ear anomalies, hemifacial microsomia, ocular defects and vertebral malformations. Although numerous chromosomal abnormalities have been associated with OAVS, no causative gene has been identified so far. OBJECTIVES: We aimed to identify the first causative gene for OAVS. METHODS: As sporadic cases are mostly described in Goldenhar syndrome, we have performed whole exome sequencing (WES) on selected affected individuals and their unaffected parents, looking for de novo mutations. Candidate gene was tested through transient knockdown experiment in zebrafish using a morpholino-based approach. A functional test was developed in cell culture in order to assess deleterious consequences of mutations. RESULTS: By WES, we identified a heterozygous nonsense mutation in one patient in the myelin transcription factor 1 (MYT1) gene. Further, we detected one heterozygous missense mutation in another patient among a cohort of 169 patients with OAVS. This gene encodes the MYT1. Functional studies by transient knockdown of myt1a, homologue of MYT1 in zebrafish, led to specific craniofacial cartilage alterations. Treatment with all-trans retinoic acid (RA), a known teratogenic agent causing OAVS, led to an upregulation of cellular endogenous MYT1 expression. Additionally, cellular wild-type MYT1 overexpression induced a downregulation of RA receptor ß (RARB), whereas mutated MYT1 did not. CONCLUSION: We report MYT1 as the first gene implicated in OAVS, within the RA signalling pathway.

Retinoic acid, RARs and early development.

Vitamin A (retinol) is an important nutrient for embryonic development and adult health. Early studies identified retinoic acid (RA) as a metabolite of retinol, however, its importance was not apparent. Later, it was observed that RA treatment of vertebrate embryos had teratogenic effects on limb development. Subsequently, the discovery of nuclear RA receptors (RARs) revealed that RA controls gene expression directly at the transcriptional level through a process referred to as RA signaling. This important discovery led to further studies demonstrating that RA and RARs are required for normal embryonic development. The determination of RA function during normal development has been challenging as RA gain-of-function studies often lead to conclusions about normal development that conflict with RAR or RA loss-of-function studies. However, genetic loss-of-function studies have identified direct target genes of endogenous RA/RAR that are required for normal development of specific tissues. Thus, genetic loss-of-function studies that eliminate RARs or RA-generating enzymes have been instrumental in revealing that RA signaling is required for normal early development of many organs and tissues, including the hindbrain, posterior body axis, somites, spinal cord, forelimbs, heart, and eye.

Role of Retinoic Acid Signaling, FGF Signaling and Meis Genes in Control of Limb Development.

The function of retinoic acid (RA) during limb development is still debated, as loss and gain of function studies led to opposite conclusions. With regard to limb initiation, genetic studies demonstrated that activation of FGF10 signaling is required for the emergence of limb buds from the trunk, with Tbx5 and RA signaling acting upstream in the forelimb field, whereas Tbx4 and Pitx1 act upstream in the hindlimb field. Early studies in chick embryos suggested that RA as well as Meis1 and Meis2 (Meis1/2) are required for subsequent proximodistal patterning of both forelimbs and hindlimbs, with RA diffusing from the trunk, functioning to activate Meis1/2 specifically in the proximal limb bud mesoderm. However, genetic loss of RA signaling does not result in loss of limb Meis1/2 expression and limb patterning is normal, although Meis1/2 expression is reduced in trunk somitic mesoderm. More recent studies demonstrated that global genetic loss of Meis1/2 results in a somite defect and failure of limb bud initiation. Other new studies reported that conditional genetic loss of Meis1/2 in the limb results in proximodistal patterning defects, and distal FGF8 signaling represses Meis1/2 to constrain its expression to the proximal limb. In this review, we hypothesize that RA and Meis1/2 both function in the trunk to initiate forelimb bud initiation, but that limb Meis1/2 expression is activated proximally by a factor other than RA and repressed distally by FGF8 to generate proximodistal patterning.

Functional and genetic analyses of ZYG11B provide evidences for its involvement in OAVS.

BACKGROUND: The Oculo-Auriculo-Vertebral Spectrum (OAVS) or Goldenhar Syndrome is an embryonic developmental disorder characterized by hemifacial microsomia associated with auricular, ocular and vertebral malformations. The clinical heterogeneity of this spectrum and its incomplete penetrance limited the molecular diagnosis. In this study, we describe a novel causative gene, ZYG11B. METHODS: A sporadic case of OAVS was analyzed by whole exome sequencing in trio strategy. The identified candidate gene, ZYG11B, was screened in 143 patients by next generation sequencing. Overexpression and immunofluorescence of wild-type and mutated ZYG11B forms were performed in Hela cells. Moreover, morpholinos were used for transient knockdown of its homologue in zebrafish embryo. RESULTS: A nonsense de novo heterozygous variant in ZYG11B, (NM_024646, c.1609G>T, p.Glu537*) was identified in a single OAVS patient. This variant leads in vitro to a truncated protein whose subcellular localization is altered. Transient knockdown of the zebrafish homologue gene confirmed its role in craniofacial cartilages architecture and in notochord development. Moreover, ZYG11B expression regulates a cartilage master regulator, SOX6, and is regulated by Retinoic Acid, a known developmental toxic molecule leading to clinical features of OAVS. CONCLUSION: Based on genetic, cellular and animal model data, we proposed ZYG11B as a novel rare causative gene for OAVS.

Prenatal retinoic acid exposure reveals candidate genes for craniofacial disorders.

Syndromes that display craniofacial anomalies comprise a major class of birth defects. Both genetic and environmental factors, including prenatal retinoic acid (RA) exposure, have been associated with these syndromes. While next generation sequencing has allowed the discovery of new genes implicated in these syndromes, some are still poorly characterized such as Oculo-Auriculo-Vertebral Spectrum (OAVS). Due to the lack of clear diagnosis for patients, developing new strategies to identify novel genes involved in these syndromes is warranted. Thus, our study aimed to explore the link between genetic and environmental factors. Owing to a similar phenotype of OAVS reported after gestational RA exposures in humans and animals, we explored RA targets in a craniofacial developmental context to reveal new candidate genes for these related disorders. Using a proteomics approach, we detected 553 dysregulated proteins in the head region of mouse embryos following their exposure to prenatal RA treatment. This novel proteomic approach implicates changes in proteins that are critical for cell survival/apoptosis and cellular metabolism which could ultimately lead to the observed phenotype. We also identified potential molecular links between three major environmental factors known to contribute to craniofacial defects including maternal diabetes, prenatal hypoxia and RA exposure. Understanding these links could help reveal common key pathogenic mechanisms leading to craniofacial disorders. Using both in vitro and in vivo approaches, this work identified two new RA targets, Gnai3 and Eftud2, proteins known to be involved in craniofacial disorders, highlighting the power of this proteomic approach to uncover new genes whose dysregulation leads to craniofacial defects.

Genomic Knockout of Two Presumed Forelimb Tbx5 Enhancers Reveals They Are Nonessential for Limb Development.

A standard approach in the identification of transcriptional enhancers is the use of transgenic animals carrying DNA elements joined to reporter genes inserted randomly in the genome. We examined elements near Tbx5, a gene required for forelimb development in humans and other vertebrates. Previous transgenic studies reported a mammalian Tbx5 forelimb enhancer located in intron 2 containing a putative retinoic acid response element and a zebrafish tbx5a forelimb (pectoral fin) enhancer located downstream that is conserved from fish to mammals. We used CRISPR/Cas9 gene editing to knockout the endogenous elements and unexpectedly found that deletion of the intron 2 and downstream elements, either singly or together in double knockouts, resulted in no effect on forelimb development. Our findings show that reporter transgenes may not identify endogenous enhancers and that in vivo genetic loss-of-function studies are required, such as CRISPR/Cas9, which is similar in effort to production of animals carrying reporter transgenes.

A novel de novo mutation in MYT1, the unique OAVS gene identified so far.

Oculo-auriculo-vertebral spectrum (OAVS) is a developmental disorder characterized by hemifacial microsomia associated with ear, eyes and vertebrae malformations showing highly variable expressivity. Recently, MYT1, encoding the myelin transcription factor 1, was reported as the first gene involved in OAVS, within the retinoic acid (RA) pathway. Fifty-seven OAVS patients originating from Brazil were screened for MYT1 variants. A novel de novo missense variant affecting function, c.323C>T (p.(Ser108Leu)), was identified in MYT1, in a patient presenting with a severe form of OAVS. Functional studies showed that MYT1 overexpression downregulated all RA receptors genes (RARA, RARB, RARG), involved in RA-mediated transcription, whereas no effect was observed on CYP26A1 expression, the major enzyme involved in RA degradation, Moreover, MYT1 variants impacted significantly the expression of these genes, further supporting their pathogenicity. In conclusion, a third variant affecting function in MYT1 was identified as a cause of OAVS. Furthermore, we confirmed MYT1 connection to RA signaling pathway.