Functional dissection of ultraconserved non-coding regions present near the Hoxd complex in vertebrates

Technical Summary: Functional Dissection of Ultraconserved Non-coding Regions near the Hoxd Complex

Problem Statement
Ultraconserved non-coding elements (UCEs) exhibit near-perfect sequence conservation across vertebrate genomes, suggesting severe evolutionary constraints and essential biological functions. Near the vertebrate Hoxd complex, a master regulator of axial patterning and limb morphogenesis, three specific Conserved Regions (CR1, CR2, and CR3) have been identified. These regions are longer than 200 bp and show exceptional sequence conservation across human, mouse, and rat genomes. While transient reporter assays have demonstrated that these regions possess robust regulatory activity (acting as enhancers or repressors), their precise in vivo functions, target genes, and mechanisms of action remain unresolved. A central question in the field is whether these elements are universally indispensable or if they exhibit functional redundancy, particularly given that previous studies in mice showed that deleting some UCEs resulted in viable, fertile animals without overt phenotypes.

Methodology
The study employed a dual bioinformatic and genetic approach using the zebrafish (Danio rerio) model system:

1. Bioinformatic Analysis: The authors utilized public datasets (UCSC Genome Browser, ChIP-seq, ChromHMM) to analyze the epigenetic and transcriptional signatures of CR1, CR2, and CR3 across multiple vertebrate species. They performed multiple sequence alignments to assess conservation patterns, analyzed SNP density, and mapped transcription factor binding sites and histone modification profiles (e.g., H3K27me3, H3K4me1/3).
2. CRISPR-Cas9 Genome Editing: To investigate in vivo function, the researchers designed multiple gRNAs targeting CR2 and CR3 individually, in combination, and spanning the intergenic regions. They generated mosaic founder embryos and screened thousands of individuals to identify germline-transmissible deletions.
3. Phenotypic and Molecular Characterization: For the successfully established mutant line (CR2∆626), the authors conducted:
   • Morphological Screening: Assessment of body axis formation, somite segmentation (via xirp2a staining), craniofacial development (Alizarin red/Alcian blue staining), and pectoral fin development.
   • Genetic Screening: Analysis of Mendelian ratios and breeding viability.
   • Expression Profiling: Spatial-temporal expression of Hoxd genes (hoxd3, hoxd4, hoxd9-13, and evx2) was analyzed using RNA in situ hybridization (RNA-ISH) at 16, 24, and 48 hours post-fertilization (hpf).
   • Quantitative Analysis: RT-qPCR was performed to detect subtle quantitative changes in Hoxd gene expression levels that might be missed by RNA-ISH.
   • Transcriptional Status: RNA-ISH with sense and antisense probes was used to determine if CR2 and CR3 are transcribed into non-coding RNAs.

Key Results

• Germline Transmission and Viability: A clear size- and conservation-dependent relationship was observed. Complete deletions of core CR2, CR3, or the combined CR2/3 loci could not be transmitted through the germline, despite extensive screening of founder embryos. In contrast, a smaller, less-conserved 626-bp subregion deletion within CR2 (termed CR2∆626) yielded a stable, homozygous-viable line.
• Phenotypic Analysis: The CR2∆626 homozygous zebrafish line displayed no overt developmental defects. There were no significant differences in:
  • Primary and secondary body axis formation.
  • Somite number and morphology.
  • Craniofacial skeletal and cartilage development.
  • Pectoral fin morphogenesis.
  • Breeding rates, clutch size, or survival ratios (Mendelian ratios were maintained).
• Gene Expression:
  • Spatial: RNA-ISH revealed that the spatial expression domains of Hoxd genes (anterior to posterior) and evx2 remained unchanged in CR2∆626 embryos compared to wild-type (WT) at 24 and 48 hpf, including in the pectoral fin buds.
  • Quantitative: While RNA-ISH showed no significant changes, RT-qPCR revealed a consistent, modest upward trend in the expression of all Hoxd complex genes at 24 and 48 hpf in the mutant line. However, these changes were not statistically significant.
• Bioinformatic Insights:
  • Conservation: The core CR2 and CR3 sequences are nearly perfectly conserved across vertebrates, whereas the deleted 626-bp subregion shows patchy conservation (strong in mammals, absent in Xenopus and zebrafish).
  • Epigenetics: The deleted subregion is uniformly repressed by Polycomb group proteins (H3K27me3) across cell types. In contrast, the retained core CR2 exhibits a "bivalent" or poised promoter signature in specific cell lines and active enhancer marks in embryonic limbs.
  • Transcription: Strand-specific RNA-ISH detected transcription of CR2 and CR3 in the same direction as Hoxd genes (opposite to evx2). Signals were observed as early as the 64-cell stage (suggesting maternal deposition) and became restricted to the anterior region (eye, brain, anterior somites) by 24 hpf.

Significance and Claims
The paper claims that the core ultraconserved sequences within CR2 and CR3 are developmentally indispensable, as their complete removal is lethal or prevents germline transmission. This finding argues against the notion that ultraconservation is phenotypically inert.

However, the study also demonstrates that the flanking, less-conserved subregions of these elements (specifically the 626-bp region of CR2) are dispensable for gross morphological development in zebrafish. The authors propose that the lack of a severe phenotype in the CR2∆626 line, despite the deletion of a conserved element, may be due to:

1. Lineage Specificity: The deleted subregion may have acquired essential regulatory roles only in higher vertebrates (e.g., for digit formation in mammals), which are absent in teleost fish.
2. Functional Redundancy: The intact core CR2 and CR3, along with broader genetic buffering networks, may compensate for the loss of the subregion.
3. Quantitative vs. Qualitative Effects: The deleted subregion may normally contribute a weak, constitutive repressive input to Hoxd transcription. Its removal results in a modest, non-significant upregulation that does not cross the threshold for overt morphological defects.

The authors conclude that within a single ultraconserved region, functional constraints are concentrated in the core, while flanking sequences may be dispensable in certain contexts. They also highlight the potential role of CR-derived non-coding RNAs in Hoxd regulation, suggesting a possible self-tethering mechanism where transcripts bind back to their genomic loci to tune transcriptional output, though this remains speculative and untested. The study underscores that the functional impact of UCE deletions depends critically on which part of the element is removed and the specific developmental context.