Long-read genome sequencing identifies cryptic structural variants in congenital aniridia cases.

BACKGROUND: Haploinsufficiency of the transcription factor PAX6 is the main cause of congenital aniridia, a genetic disorder characterized by iris and foveal hypoplasia. 11p13 microdeletions altering PAX6 or its downstream regulatory region (DRR) are present in about 25% of patients; however, only a few complex rearrangements have been described to date. Here, we performed nanopore-based whole-genome sequencing to assess the presence of cryptic structural variants (SVs) on the only two unsolved "PAX6-negative" cases from a cohort of 110 patients with congenital aniridia after unsuccessfully short-read sequencing approaches. RESULTS: Long-read sequencing (LRS) unveiled balanced chromosomal rearrangements affecting the PAX6 locus at 11p13 in these two patients and allowed nucleotide-level breakpoint analysis. First, we identified a cryptic 4.9 Mb de novo inversion disrupting intron 7 of PAX6, further verified by targeted polymerase chain reaction amplification and sequencing and FISH-based cytogenetic analysis. Furthermore, LRS was decisive in correctly mapping a t(6;11) balanced translocation cytogenetically detected in a second proband with congenital aniridia and considered non-causal 15 years ago. LRS resolved that the breakpoint on chromosome 11 was indeed located at 11p13, disrupting the DNase I hypersensitive site 2 enhancer within the DRR of PAX6, 161 Kb from the causal gene. Patient-derived RNA expression analysis demonstrated PAX6 haploinsufficiency, thus supporting that the 11p13 breakpoint led to a positional effect by cleaving crucial enhancers for PAX6 transactivation. LRS analysis was also critical for mapping the exact breakpoint on chromosome 6 to the highly repetitive centromeric region at 6p11.1. CONCLUSIONS: In both cases, the LRS-based identified SVs have been deemed the hidden pathogenic cause of congenital aniridia. Our study underscores the limitations of traditional short-read sequencing in uncovering pathogenic SVs affecting low-complexity regions of the genome and the value of LRS in providing insight into hidden sources of variation in rare genetic diseases.

Opposite T3 Response of ACTG1-FOS Subnetwork Differentiate Tailfin Fate in Xenopus Tadpole and Post-hatching Axolotl.

Amphibian post-embryonic development and Thyroid Hormones (TH) signaling are deeply and intimately connected. In anuran amphibians, TH induce the spectacular and complex process known as metamorphosis. In paedomorphic salamanders, at similar development time, raising levels of TH fail to induce proper metamorphosis, as many "larval" tissues (e.g., gills, tailfin) are maintained. Why does the same evolutionary conserved signaling pathway leads to alternative phenotypes? We used a combination of developmental endocrinology, functional genomics and network biology to compare the transcriptional response of tailfin to TH, in the post-hatching paedormorphic Axolotl salamander and Xenopus tadpoles. We also provide a technological framework that efficiently reduces large lists of regulated genes down to a few genes of interest, which is well-suited to dissect endocrine regulations. We first show that Axolotl tailfin undergoes a strong and robust TH-dependent transcriptional response at post embryonic transition, despite the lack of visible anatomical changes. We next show that Fos and Actg1, which structure a single and dense subnetwork of cellular sensors and regulators, display opposite regulation between the two species. We finally show that TH treatments and natural variations of TH levels follow similar transcriptional dynamics. We suggest that, at the molecular level, tailfin fate correlates with the alternative transcriptional states of an fos-actg1 sub-network, which also includes transcription factors and regulators of cell fate. We propose that this subnetwork is one of the molecular switches governing the initiation of distinct TH responses, with transcriptional programs conducting alternative tailfin fate (maintenance vs. resorption) 2 weeks post-hatching.

Developmental and cancer-associated plasticity of DNA replication preferentially targets GC-poor, lowly expressed and late-replicating regions.

The spatiotemporal program of metazoan DNA replication is regulated during development and altered in cancers. We have generated novel OK-seq, Repli-seq and RNA-seq data to compare the DNA replication and gene expression programs of twelve cancer and non-cancer human cell types. Changes in replication fork directionality (RFD) determined by OK-seq are widespread but more frequent within GC-poor isochores and largely disconnected from transcription changes. Cancer cell RFD profiles cluster with non-cancer cells of similar developmental origin but not with different cancer types. Importantly, recurrent RFD changes are detected in specific tumour progression pathways. Using a model for establishment and early progression of chronic myeloid leukemia (CML), we identify 1027 replication initiation zones (IZs) that progressively change efficiency during long-term expression of the BCR-ABL1 oncogene, being twice more often downregulated than upregulated. Prolonged expression of BCR-ABL1 results in targeting of new IZs and accentuation of previous efficiency changes. Targeted IZs are predominantly located in GC-poor, late replicating gene deserts and frequently silenced in late CML. Prolonged expression of BCR-ABL1 results in massive deletion of GC-poor, late replicating DNA sequences enriched in origin silencing events. We conclude that BCR-ABL1 expression progressively affects replication and stability of GC-poor, late-replicating regions during CML progression.