Complex human chromosomal and genomic rearrangements

Extrachromosomal circular DNAs in prostate adenocarcinoma: global characterizations and a novel prediction model

Background

The role of focal amplifications and extrachromosomal circular DNA (eccDNA) is still uncertain in prostate adenocarcinoma (PRAD). Here, we first mapped the global characterizations of eccDNA and then investigate the characterization of eccDNA-amplified key differentially expressed encoded genes (eKDEGs) in the progression, immune response and immunotherapy of PRAD.

Methods

Circular_seq was used in conjunction with the TCGA-PRAD transcriptome dataset to sequence, annotate, and filter for eccDNA-amplified differentially expressed coding genes (eDEGs) in PRAD and para-cancerous normal prostate tissues. Afterwards, risk models were created and eKDEGs linked to the PRAD prognosis were identified using Cox and Lasso regression analysis. The immune microenvironment of the risk model was quantified using a variety of immunological algorithms, which also identified its characteristics with regard to immunotherapy, immune response, and immune infiltration.

Results

In this research, there was no significant difference in the size, type, and chromosomal distribution of eccDNA in PRAD and para-cancerous normal prostate tissues. However, 4,290 differentially expressed eccDNAs were identified and 1,981 coding genes were amplified. Following that, 499 eDEGs were tested in conjunction with the transcriptome dataset from TCGA-PRAD. By using Cox and Lasso regression techniques, ZNF330 and PITPNM3 were identified as eKDEGs of PRAD, and a new PRAD risk model was conducted based on this. Survival analysis showed that the high-risk group of this model was associated with poor prognosis and validated in external data. Immune infiltration analysis showed that the model risks affected immune cell infiltration in PRAD, not only mediating changes in immune cell function, but also correlating with immunophenotyping. Furthermore, the high-risk group was negatively associated with anti-CTLA-4/anti-PD-1 response and mutational burden. In addition, Tumor Immune Dysfunction and Exclusion analyses showed that high-risk group was more prone to immune escape. Drug sensitivity analyses identified 10 drugs, which were instructive for PRAD treatment.

Conclusion

ZNF330 and PITPNM are the eKDEGs for PRAD, which can be used as potential new prognostic markers. The two-factor combined risk model can effectively assess the survival and prognosis of PRAD patients, but also can predict the different responses of immunotherapy to PRAD patients, which may provide new ideas for PRAD immunotherapy.

Genomic Balancing Act: deciphering DNA rearrangements in the complex chromosomal aberration involving 5p15.2, 2q31.1, and 18q21.32

Despite extensive research into the genetic underpinnings of neurodevelopmental disorders (NDD), many clinical cases remain unresolved. We studied a female proband with a NDD, mildly dysmorphic facial features, and brain stem hypoplasia on neuroimaging. Comprehensive genomic analyses revealed a terminal 5p loss and a terminal 18q gain in the proband while a diploid copy number for chromosomes 5 and 18 in both parents. Genomic investigations in the proband identified an unbalanced translocation t(5;18) with additional genetic material from chromosome 2 (2q31.3) inserted at the breakpoint, pointing to a complex chromosomal rearrangement (CCR) involving 5p15.2, 2q31.3, and 18q21.32. Breakpoint junction analyses enabled by long-read genome sequencing unveiled the presence of four distinct junctions in the father, who is a carrier of a balanced CCR. The proband inherited from the father both the abnormal chromosome 5 resulting in segmental aneusomies of chr5 (loss) and chr18 (gain) and a der(2) homologue. Evidences suggest a chromoplexy mechanism for this CCR derivation, involving double-strand breaks (DSBs) repaired by non-homologous end joining (NHEJ) or alternative end joining (alt-EJ). The complexity of the CCR and the segregation of homologues elucidate the genetic model for this family. This study demonstrates the importance of combining multiple genomic technologies to uncover genetic causes of complex neurodevelopmental syndromes and to better understand genetic disease mechanisms.

Long-read sequencing and optical mapping generates near T2T assemblies that resolves a centromeric translocation

Long-read genome sequencing (lrGS) is a promising method in genetic diagnostics. Here we investigate the potential of lrGS to detect a disease-associated chromosomal translocation between 17p13 and the 19 centromere. We constructed two sets of phased and non-phased de novo assemblies; (i) based on lrGS only and (ii) hybrid assemblies combining lrGS with optical mapping using lrGS reads with a median coverage of 34X. Variant calling detected both structural variants (SVs) and small variants and the accuracy of the small variant calling was compared with those called with short-read genome sequencing (srGS). The de novo and hybrid assemblies had high quality and contiguity with N50 of 62.85 Mb, enabling a near telomere to telomere assembly with less than a 100 contigs per haplotype. Notably, we successfully identified the centromeric breakpoint of the translocation. A concordance of 92% was observed when comparing small variant calling between srGS and lrGS. In summary, our findings underscore the remarkable potential of lrGS as a comprehensive and accurate solution for the analysis of SVs and small variants. Thus, lrGS could replace a large battery of genetic tests that were used for the diagnosis of a single symptomatic translocation carrier, highlighting the potential of lrGS in the realm of digital karyotyping.

Structural rearrangements as a recurrent pathogenic mechanism for SETBP1 haploinsufficiency

Chromosomal structural rearrangements consist of anomalies in genomic architecture that may or may not be associated with genetic material gain and loss. Evaluating the precise breakpoint is crucial from a diagnostic point of view, highlighting possible gene disruption and addressing to appropriate genotype-phenotype association. Structural rearrangements can either occur randomly within the genome or present with a recurrence, mainly due to peculiar genomic features of the surrounding regions. We report about three non-related individuals, harboring chromosomal structural rearrangements interrupting SETBP1, leading to gene haploinsufficiency. Two out of them resulted negative to Chromosomal Microarray Analysis (CMA), being the rearrangement balanced at a microarray resolution. The third one, presenting with a complex three-chromosome rearrangement, had been previously diagnosed with SETBP1 haploinsufficiency due to a partial gene deletion at one of the chromosomal breakpoints. We thoroughly characterized the rearrangements by means of Optical Genome Mapping (OGM) and Whole Genome Sequencing (WGS), providing details about the involved sequences and the underlying mechanisms. We propose structural variants as a recurrent event in SETBP1 haploinsufficiency, which may be overlooked by laboratory routine genomic analyses (CMA and Whole Exome Sequencing) or only partially determined when associated with genomic losses at breakpoints. We finally introduce a possible role of SETBP1 in a Noonan-like phenotype.

A unifying model that explains the origins of human inverted copy number variants

With the release of the telomere-to-telomere human genome sequence and the availability of both long-read sequencing and optical genome mapping techniques, the identification of copy number variants (CNVs) and other structural variants is providing new insights into human genetic disease. Different mechanisms have been proposed to account for the novel junctions in these complex architectures, including aberrant forms of DNA replication, non-allelic homologous recombination, and various pathways that repair DNA breaks. Here, we have focused on a set of structural variants that include an inverted segment and propose that they share a common initiating event: an inverted triplication with long, unstable palindromic junctions. The secondary rearrangement of these palindromes gives rise to the various forms of inverted structural variants. We postulate that this same mechanism (ODIRA: origin-dependent inverted-repeat amplification) that creates the inverted CNVs in inherited syndromes also generates the palindromes found in cancers.

Recent gene selection and drug resistance underscore clinical adaptation across Candida species

Understanding how microbial pathogens adapt to treatments, humans and clinical environments is key to infer mechanisms of virulence, transmission and drug resistance. This may help improve therapies and diagnostics for infections with a poor prognosis, such as those caused by fungal pathogens, including Candida. Here we analysed genomic variants across approximately 2,000 isolates from six Candida species (C. glabrata, C. auris, C. albicans, C. tropicalis, C. parapsilosis and C. orthopsilosis) and identified genes under recent selection, suggesting a highly complex clinical adaptation. These involve species-specific and convergently affected adaptive mechanisms, such as adhesion. Using convergence-based genome-wide association studies we identified known drivers of drug resistance alongside potentially novel players. Finally, our analyses reveal an important role of structural variants and suggest an unexpected involvement of (para)sexual recombination in the spread of resistance. Our results provide insights on how opportunistic pathogens adapt to human-related environments and unearth candidate genes that deserve future attention.

Complex Chromosomal Rearrangement Involving Chromosomes 10 and 11, Accompanied by Two Adjacent 11p14.1p13 and 11p13p12 Deletions, Identified in a Patient with WAGR Syndrome

Three years ago, our patient, at that time a 16-month-old boy, was discovered to have bilateral kidney lesions with a giant tumor in the right kidney. Chemotherapy and bilateral nephron-sparing surgery (NSS) for Wilms tumor with nephroblastomatosis was carried out. The patient also had eye affection, including glaucoma, eye enlargement, megalocornea, severe corneal swelling and opacity, complete aniridia, and nystagmus. The diagnosis of WAGR syndrome was suspected. De novo complex chromosomal rearrangement with balanced translocation t(10,11)(p15;p13) and a pericentric inversion inv(11)(p13q12), accompanied by two adjacent 11p14.1p13 and 11p13p12 deletions, were identified. Deletions are raised through the complex molecular mechanism of two subsequent rearrangements affecting chromosomes 11 and 10. WAGR syndrome diagnosis was clinically and molecularly confirmed, highlighting the necessity of comprehensive genetic testing in patients with congenital aniridia and/or WAGR syndrome.

Genomic structural variation: A complex but important driver of human evolution

Structural variants (SVs)-including duplications, deletions, and inversions of DNA-can have significant genomic and functional impacts but are technically difficult to identify and assay compared with single-nucleotide variants. With the aid of new genomic technologies, it has become clear that SVs account for significant differences across and within species. This phenomenon is particularly well-documented for humans and other primates due to the wealth of sequence data available. In great apes, SVs affect a larger number of nucleotides than single-nucleotide variants, with many identified SVs exhibiting population and species specificity. In this review, we highlight the importance of SVs in human evolution by (1) how they have shaped great ape genomes resulting in sensitized regions associated with traits and diseases, (2) their impact on gene functions and regulation, which subsequently has played a role in natural selection, and (3) the role of gene duplications in human brain evolution. We further discuss how to incorporate SVs in research, including the strengths and limitations of various genomic approaches. Finally, we propose future considerations in integrating existing data and biospecimens with the ever-expanding SV compendium propelled by biotechnology advancements.

Applications of advanced technologies for detecting genomic structural variation

Chromosomal structural variation (SV) encompasses a heterogenous class of genetic variants that exerts strong influences on human health and disease. Despite their importance, many structural variants (SVs) have remained poorly characterized at even a basic level, a discrepancy predicated upon the technical limitations of prior genomic assays. However, recent advances in genomic technology can identify and localize SVs accurately, opening new questions regarding SV risk factors and their impacts in humans. Here, we first define and classify human SVs and their generative mechanisms, highlighting characteristics leveraged by various SV assays. We next examine the first-ever gapless assembly of the human genome and the technical process of assembling it, which required third-generation sequencing technologies to resolve structurally complex loci. The new portions of that "telomere-to-telomere" and subsequent pangenome assemblies highlight aspects of SV biology likely to develop in the near-term. We consider the strengths and limitations of the most promising new SV technologies and when they or longstanding approaches are best suited to meeting salient goals in the study of human SV in population-scale genomics research, clinical, and public health contexts. It is a watershed time in our understanding of human SV when new approaches are expected to fundamentally change genomic applications.

Identification of complex and cryptic chromosomal rearrangements by optical genome mapping

BACKGROUND

Optical genome mapping (OGM) has developed into a highly promising method for detecting structural variants (SVs) in human genomes. Complex chromosomal rearrangements (CCRs) and cryptic translocations are rare events that are considered difficult to detect by routine cytogenetic methods. In this study, OGM was applied to delineate the precise chromosomal rearrangements in three cases with uncertain or unconfirmed CCRs detected by conventional karyotyping and one case with a cryptic translocation suggested by fetal chromosomal microarray analysis (CMA).

RESULTS

In the three cases with CCRs, OGM not only confirmed or revised the original karyotyping results but also refined the precise chromosomal structures. In the case with a suspected translocation not detected by karyotyping, OGM efficiently identified the cryptic translocation and defined the genomic breakpoints with relatively high accuracy.

CONCLUSIONS

Our study confirmed OGM as a robust alternative approach to karyotyping for the detection of chromosomal structural rearrangements, including CCRs and cryptic translocations.

Quantitative trait locus for calving traits on Bos taurus autosome 18 in Holstein cattle is embedded in a complex genomic region

Although the quantitative trait locus (QTL) on chromosome 18 (BTA18) associated with paternal calving ease and stillbirth in Holstein Friesian cattle and its cross has been known for over 20 years, to our knowledge, the exact causal genetic sequence has yet escaped identification. The aim of this study was to re-examine the region of the published QTL on BTA18 and to investigate the possible reasons behind this elusiveness. For this purpose, we carried out a combined linkage disequilibrium and linkage analysis using genotyping data of 2,697 German Holstein Friesian (HF) animals and subsequent whole-genome sequencing (WGS) data analyses and genome assembly of HF samples. We confirmed the known QTL in the 95% confidence interval of 1.089 Mbp between 58.34 and 59.43 Mbp on BTA18. Additionally, these 4 SNPs in the near-perfect linkage disequilibrium with the QTL haplotype were identified: rs381577268 (on 57,816,137 bp, C/T), rs381878735 (on 59,574,329 bp, A/T), rs464221818 (on 59,329,176 bp, C/T), and rs472502785 (on 59,345,689 bp, T/C). Search for the causal mutation using short and long-read sequences, and methylation data of the BTA18 QTL region did not reveal any candidates though. The assembly showed problems in the region, as well as an abundance of segmental duplications within and around the region. Taking the QTL of BTA18 in Holstein cattle as an example, the data presented in this study comprehensively characterize the genomic features that could also be relevant for other such elusive QTL in various other cattle breeds and livestock species as well.

Investigation of the mechanism of copy number variations involving the α-globin gene cluster on chromosome 16: two case reports and literature review

Thalassemia is one of the most common single-gene disorder worldwide. An important genetic cause of thalassemia is copy number variations (CNVs) in the α-globin gene cluster. However, there is no unified summary and discussion on the detailed information and mechanisms of these CNVs. In this study, two novel CNVs, a tandem duplication (αααα¹⁵⁹) and deletion (--²⁵⁹), were identified in two Chinese families with thalassemia patients, according to the results of hematologic analysis, routine genetic testing for thalassemia, multiplex ligation-dependent probe amplification (MLPA), next-generation sequencing (NGS) and other molecular methods. Co-inherited with β^CD41-42 mutation and --^SEA deletion separately, αααα¹⁵⁹ and --²⁵⁹ resulted in a patient with β-thalassemia intermedia and a lethal fetus with Hb Bart's hydrops fetalis syndrome, respectively. Next, a literature review was performed to summarize all known CNVs involving the α-globin gene cluster. The molecular structure characteristics of these CNVs were analyzed and the possible mechanism was explored. It is the first time to analyze the generation mechanism of genome arrangements in the α-globin gene cluster systematically.

Complex genomic rearrangements: an underestimated cause of rare diseases

Complex genomic rearrangements (CGRs) are known contributors to disease but are often missed during routine genetic screening. Identifying CGRs requires (i) identifying copy number variants (CNVs) concurrently with inversions, (ii) phasing multiple breakpoint junctions incis, as well as (iii) detecting and resolving structural variants (SVs) within repeats. We demonstrate how combining cytogenetics and new sequencing methodologies is being successfully applied to gain insights into the genomic architecture of CGRs. In addition, we review CGR patterns and molecular features revealed by studying constitutional genomic disorders. These data offer invaluable lessons to individuals interested in investigating CGRs, evaluating their clinical relevance and frequency, as well as assessing their impact(s) on rare genetic diseases.

A Complex Genomic Rearrangement Resulting in Loss of Function of SCN1A and SCN2A in a Patient with Severe Developmental and Epileptic Encephalopathy

Complex genomic rearrangements (CGRs) are structural variants arising from two or more chromosomal breaks, which are challenging to characterize by conventional or molecular cytogenetic analysis (karyotype and FISH). The integrated approach of standard and genomic techniques, including optical genome mapping (OGM) and genome sequencing, is crucial for disclosing and characterizing cryptic chromosomal rearrangements at high resolutions. We report on a patient with a complex developmental and epileptic encephalopathy in which karyotype analysis showed a de novo balanced translocation involving the long arms of chromosomes 2 and 18. Microarray analysis detected a 194 Kb microdeletion at 2q24.3 involving the SCN2A gene, which was considered the likely translocation breakpoint on chromosome 2. However, OGM redefined the translocation breakpoints by disclosing a paracentric inversion at 2q24.3 disrupting SCN1A. This combined genomic high-resolution approach allowed a fine characterization of the CGR, which involves two different chromosomes with four breakpoints. The patient’s phenotype resulted from the concomitant loss of function of SCN1A and SCN2A.

Extrachromosomal circular DNA: Current status and future prospects

Extrachromosomal circular DNA (eccDNA) is a double-stranded DNA molecule found in various organisms, including humans. In the past few decades, the research on eccDNA has mainly focused on cancers and their associated diseases. Advancements in modern omics technologies have reinvigorated research on eccDNA and shed light on the role of these molecules in a range of diseases and normal cell phenotypes. In this review, we first summarize the formation of eccDNA and its modes of action in eukaryotic cells. We then outline eccDNA as a disease biomarker and reveal its regulatory mechanism. We finally discuss the future prospects of eccDNA, including basic research and clinical application. Thus, with the deepening of understanding and exploration of eccDNAs, they hold great promise in future biomedical research and clinical translational application.

Biology in balance: human diploid genome integrity, gene dosage, and genomic medicine

The path to completion of the functional annotation of the haploid human genome reference build, exploration of the clan genomics hypothesis, understanding human gene and genome functional biology, and gene genome and organismal evolution, is in reach.

Computational Approaches for Understanding Sequence Variation Effects on the 3D Genome Architecture

Decoding how genomic sequence and its variations affect 3D genome architecture is indispensable for understanding the genetic architecture of various traits and diseases. The 3D genome organization can be significantly altered by genome variations and in turn impact the function of the genomic sequence. Techniques for measuring the 3D genome architecture across spatial scales have opened up new possibilities for understanding how the 3D genome depends upon the genomic sequence and how it can be altered by sequence variations. Computational methods have become instrumental in analyzing and modeling the sequence effects on 3D genome architecture, and recent development in deep learning sequence models have opened up new opportunities for studying the interplay between sequence variations and the 3D genome. In this review, we focus on computational approaches for both the detection and modeling of sequence variation effects on the 3D genome, and we discuss the opportunities presented by these approaches. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Novel gross deletion at the LHX4 gene locus in a child with growth hormone deficiency

OBJECTIVE

To identify and characterize a novel deletion at the LHX4 gene locus in a proband with growth hormone deficiency (GHD).

METHODS

Long range polymerase chain reaction (PCR) amplification was used to confirm the suspected deletion and to identify the rough locations of the end points. Sanger sequencing was carried out to identify the exact end points of the deletion.

RESULTS

Suspected deletion was confirmed via long range PCR amplification. Sanger sequencing identified the end points of the deletion within three nucleotide repeat sequences ("CTT"). The total length of the deleted segment was 12 127 base pairs and it includes complete exon 5 and exon 6 of the LHX4 gene. Therefore the homeodomain motif coded by exons 4 and 5, might be affected.

CONCLUSION

We have identified a novel deletion that spans exon 5 and exon 6 of the LHX4 gene that could have occurred via microhomology mediated non-recurrent rearrangement. The deletion characterized does not appear to have been reported before. To our knowledge this novel deletion is the first identified LHX4 variant from Sri Lanka and it explains the phenotype of the proband characterized by growth hormone deficiency, hypoplastic anterior pituitary and subsequent deficiency of thyroid stimulating hormone and adrenocorticotropic hormone (ACTH).

Xq27.1 palindrome mediated interchromosomal insertion likely causes familial congenital bilateral laryngeal abductor paralysis (Plott syndrome)

Bilateral laryngeal abductor paralysis is a rare entity and the second most common cause of stridor in newborns. So far, no conclusive genetic or chromosomal aberration has been reported for X-linked isolated bilateral vocal cord paralysis, also referred to as Plott syndrome. Via whole genome sequencing (WGS), we identified a complex interchromosomal insertion in a large family with seven affected males. The 404 kb inserted fragment originates from chromosome 10q21.3, contains no genes and is inserted inversionally into the intergenic chromosomal region Xq27.1, 82 kb centromeric to the nearest gene SOX3. The patterns found at the breakpoint junctions resemble typical characteristics that arise in replication-based mechanisms with long-distance template switching. Non protein-coding insertions into the same genomic region have been described to result in different phenotypes, indicating that the phenotypic outcome likely depends on the introduction of regulatory elements. In conclusion, our data adds Plott syndrome as another entity, likely caused by the insertion of non-coding DNA into the intergenic chromosomal region Xq27.1. In this regard, we demonstrate the importance of WGS as a powerful diagnostic test in unsolved genetic diseases, as this genomic rearrangement has not been detected by current first-line diagnostic tests, i.e., exome sequencing and chromosomal microarray analysis.

Stimulation of adaptive gene amplification by origin firing under replication fork constraint

Adaptive mutations can cause drug resistance in cancers and pathogens, and increase the tolerance of agricultural pests and diseases to chemical treatment. When and how adaptive mutations form is often hard to discern, but we have shown that adaptive copy number amplification of the copper resistance gene CUP1 occurs in response to environmental copper due to CUP1 transcriptional activation. Here we dissect the mechanism by which CUP1 transcription in budding yeast stimulates copy number variation (CNV). We show that transcriptionally stimulated CNV requires TREX-2 and Mediator, such that cells lacking TREX-2 or Mediator respond normally to copper but cannot acquire increased resistance. Mediator and TREX-2 can cause replication stress by tethering transcribed loci to nuclear pores, a process known as gene gating, and transcription at the CUP1 locus causes a TREX-2-dependent accumulation of replication forks indicative of replication fork stalling. TREX-2-dependent CUP1 gene amplification occurs by a Rad52 and Rad51-mediated homologous recombination mechanism that is enhanced by histone H3K56 acetylation and repressed by Pol32 and Pif1. CUP1 amplification is also critically dependent on late-firing replication origins present in the CUP1 repeats, and mutations that remove or inactivate these origins strongly suppress the acquisition of copper resistance. We propose that replicative stress imposed by nuclear pore association causes replication bubbles from these origins to collapse soon after activation, leaving a tract of H3K56-acetylated chromatin that promotes secondary recombination events during elongation after replication fork re-start events. The capacity for inefficient replication origins to promote copy number variation renders certain genomic regions more fragile than others, and therefore more likely to undergo adaptive evolution through de novo gene amplification.

De novo mutation rates at the single-mutation resolution in a human HBB gene-region associated with adaptation and genetic disease

Although it is known that the mutation rate varies across the genome, previous estimates were based on averaging across various numbers of positions. Here, we describe a method to measure the origination rates of target mutations at target base positions and apply it to a 6-bp region in the human hemoglobin subunit beta (HBB) gene and to the identical, paralogous hemoglobin subunit delta (HBD) region in sperm cells from both African and European donors. The HBB region of interest (ROI) includes the site of the hemoglobin S (HbS) mutation, which protects against malaria, is common in Africa, and has served as a classic example of adaptation by random mutation and natural selection. We found a significant correspondence between de novo mutation rates and past observations of alleles in carriers, showing that mutation rates vary substantially in a mutation-specific manner that contributes to the site frequency spectrum. We also found that the overall point mutation rate is significantly higher in Africans than in Europeans in the HBB region studied. Finally, the rate of the 20A→T mutation, called the “HbS mutation” when it appears in HBB, is significantly higher than expected from the genome-wide average for this mutation type. Nine instances were observed in the African HBB ROI, where it is of adaptive significance, representing at least three independent originations; no instances were observed elsewhere. Further studies will be needed to examine mutation rates at the single-mutation resolution across these and other loci and organisms and to uncover the molecular mechanisms responsible.

RNA-driven JAZF1-SUZ12 gene fusion in human endometrial stromal cells

Oncogenic fusion genes as the result of chromosomal rearrangements are important for understanding genome instability in cancer cells and developing useful cancer therapies. To date, the mechanisms that create such oncogenic fusion genes are poorly understood. Previously we reported an unappreciated RNA-driven mechanism in human prostate cells in which the expression of chimeric RNA induces specified gene fusions in a sequence-dependent manner. One fundamental question yet to be addressed is whether such RNA-driven gene fusion mechanism is generalizable, or rather, a special case restricted to prostate cells. In this report, we demonstrated that the expression of designed chimeric RNAs in human endometrial stromal cells leads to the formation of JAZF1-SUZ12, a cancer fusion gene commonly found in low-grade endometrial stromal sarcomas. The process is specified by the sequence of chimeric RNA involved and inhibited by estrogen or progesterone. Furthermore, it is the antisense rather than sense chimeric RNAs that effectively drive JAZF1-SUZ12 gene fusion. The induced fusion gene is validated both at the RNA and the genomic DNA level. The ability of designed chimeric RNAs to drive and recapitulate the formation of JAZF1-SUZ12 gene fusion in endometrial cells represents another independent case of RNA-driven gene fusion, suggesting that RNA-driven genomic recombination is a permissible mechanism in mammalian cells. The results could have fundamental implications in the role of RNA in genome stability, and provide important insight in early disease mechanisms related to the formation of cancer fusion genes.

Fusion genes resulting from chromosomal translocations are important for understanding cancer mechanisms and developing anti-cancer therapies. Fusion gene are presumed to occur prior to fusion RNA expression. However, studies have reported the presence of fusion RNAs in individuals who were negative for chromosomal translocations. The observation, that fusion RNA could be present prior to fusion gene formation, raises the provocative hypothesis that fusion RNA, or any cellular RNA with sequence compositions resembling that of fusion RNA, could act as a template to mediate genomic rearrangement which leads to the final gene fusion. In this report, we demonstrated that the expression of designed chimeric RNAs in human endometrial stromal cells leads to the formation of JAZF1-SUZ12, a cancer fusion gene found in endometrial stromal sarcomas. The process is specified by the sequence of chimeric RNA involved and inhibited by estrogen or progesterone. Furthermore, it is the antisense rather than sense chimeric RNAs that effectively drive JAZF1-SUZ12 gene fusion. The results could have fundamental implications in the role of RNA in mammalian genome stability, provide important insight in early disease mechanism, as well as developing gene editing technology via mechanisms native to mammalian cells.

An Apparently Balanced Complex Chromosome Rearrangement Involving Seven Breaks and Four Chromosomes in a Healthy Female and Segregation/Recombination in Her Affected Son

Duplication of the distal 1q and 4p segments are both characterized by the presence of intellectual disability/neurodevelopmental delay and dysmorphisms. Here, we describe a male with a complex chromosome rearrangement (CCR) presenting with overlapping clinical findings between these 2 syndromes. In order to better characterize this CCR, classical karyotyping, FISH, and chromosomal microarray analysis were performed on material from the patient and his parents, which revealed an unbalanced karyotype with duplications at 1q41q43 and 4p15.2p14 in the proband. The rearrangements, which were derived from a maternal balanced karyotype, included an insertion of a segment from the long to the short arm of chromosome 1, a balanced translocation involving chromosomes 14 and 18, and an insertion of a segment from the short arm of chromosome 4 into the derived chromosome 14. This study aimed to better define the clinical history and prognosis of a patient with this rare category of chromosomal aberration. Our results suggest that the frequency of CCR in the general population may be underestimated; when balanced, they may not have a phenotypic effect. Moreover, they emphasize the need for cytogenetic techniques complementary to chromosomal microarray for proper genetic counseling.

Reducing variability of breast cancer subtype predictors by grounding deep learning models in prior knowledge

Deep learning neural networks have improved performance in many cancer informatics problems, including breast cancer subtype classification. However, many networks experience underspecificationwheremultiplecombinationsofparametersachievesimilarperformance, bothin training and validation. Additionally, certain parameter combinations may perform poorly when the test distribution differs from the training distribution. Embedding prior knowledge from the literature may address this issue by boosting predictive models that provide crucial, in-depth information about a given disease. Breast cancer research provides a wealth of such knowledge, particularly in the form of subtype biomarkers and genetic signatures. In this study, we draw on past research on breast cancer subtype biomarkers, label propagation, and neural graph machines to present a novel methodology for embedding knowledge into machine learning systems. We embed prior knowledge into the loss function in the form of inter-subject distances derived from a well-known published breast cancer signature. Our results show that this methodology reduces predictor variability on state-of-the-art deep learning architectures and increases predictor consistency leading to improved interpretation. We find that pathway enrichment analysis is more consistent after embedding knowledge. This novel method applies to a broad range of existing studies and predictive models. Our method moves the traditional synthesis of predictive models from an arbitrary assignment of weights to genes toward a more biologically meaningful approach of incorporating knowledge.

Levodopa-responsive dystonia caused by biallelic PRKN exon inversion invisible to exome sequencing

Biallelic pathogenic variants in PRKN (PARK2), encoding the E3 ubiquitin ligase parkin, lead to early-onset Parkinson's disease. Structural variants, including duplications or deletions, are common in PRKN due to their location within the fragile site FRA6E. These variants are readily detectable by copy number variation analysis. We studied four siblings with levodopa-responsive dystonia by exome sequencing followed by genome sequencing. Affected individuals developed juvenile levodopa-responsive dystonia with subsequent appearance of parkinsonism and motor fluctuations that improved by subthalamic stimulation. Exome sequencing and copy number variation analysis were not diagnostic, yet revealed a shared homozygous block including PRKN. Genome sequencing revealed an inversion within PRKN, with intronic breakpoints flanking exon 5. Breakpoint junction analysis implicated non-homologous end joining and possibly replicative mechanisms as the repair pathways involved. Analysis of cDNA indicated skipping of exon 5 (84 bp) that was replaced by 93 bp of retained intronic sequence, preserving the reading frame yet altering a significant number of residues. Balanced copy number inversions in PRKN are associated with a severe phenotype. Such structural variants, undetected by exome analysis and by copy number variation analysis, should be considered in the relevant clinical setting. These findings raise the possibility that PRKN structural variants are more common than currently estimated.

Double-strand breaks induce short-scale DNA replication and damage amplification in the fully grown mouse oocytes

Break-induced replication (BIR) is essential for the repair of DNA double-strand breaks (DSBs) with single ends. DSBs-induced microhomology-mediated BIR (mmBIR) and template-switching can increase the risk of complex genome rearrangement. In addition, DSBs can also induce the multi-invasion-mediated DSB amplification. The mmBIR-induced genomic rearrangement has been identified in cancer cells and patients with rare diseases. However, when and how mmBIR is initiated have not been fully and deeply studied. Furthermore, it is not well understood about the conditions for initiation of multi-invasion-mediated DSB amplification. In the G2 phase oocyte of mouse, we identified a type of short-scale BIR (ssBIR) using the DNA replication indicator 5-ethynyl-2'-deoxyuridine (EdU). These ssBIRs could only be induced in the fully grown oocytes but not the growing oocytes. If the DSB oocytes were treated with Rad51 or Chek1/2 inhibitors, both EdU signals and DSB marker γH2A.X foci would decrease. In addition, the DNA polymerase inhibitor Aphidicolin could inhibit the ssBIR and another inhibitor ddATP could reduce the number of γH2A.X foci in the DSB oocytes. In conclusion, our results showed that DNA DSBs in the fully grown oocytes can initiate ssBIR and be amplified by Rad51 or DNA replication.

Study of complex structural variations of X-linked deafness-2 based on single-molecule sequencing

X-linked deafness-2 (DFNX2) is cochlear incomplete partition type III (IP-III), one of inner ear malformations characterized by an abnormally wide opening in the bone separating the basal turn of the cochlea from the internal auditory canal, fixation of the stapes and cerebrospinal fluid (CSF) gusher upon stapedectomy or cochleostomy. The causative gene of DFNX2 was POU3F4. To investigate the genetic causes of DFNX2 and compare the efficiency of different sequencing methods, 12 unrelated patients were enrolled in the present study. Targeted next-generation sequencing (NGS) and long-read sequencing were used to analyze the genetic etiology of DFNX2. Six variants of POU3F4 were identified in this cohort by NGS. Three patients with a negative diagnosis based on NGS were enrolled in further long-read sequencing. Two of them were all found to carry structural variations (SVs) on chromosome X, consisting of an 870-kb deletion (DEL) at upstream of POU3F4 and an 8-Mb inversion (INV). The 870-kb DEL may have been due to non-homologous end joining (NHEJ), while non-allelic homologous recombination (NAHR) within a single chromatid may have accounted for the 8-Mb INV. Common POU3F4 mutations in DFNX2 included point mutations, small insertions and deletions (INDELs), and exon mutations, which can be detected by Sanger sequencing and NGS. Single-molecule long-read sequencing constitutes an additional and valuable method for accurate detection of pathogenic SVs in IP-III patients with negative NGS results.

Congenital chromoanagenesis in the routine postnatal chromosomal microarray analyses

Chromosomal microarray analyses (CMA) have greatly increased both the yield and diagnostic accuracy of postnatal analysis; it has been used as a first-tier cytogenetic test in patients with intellectual disability, autism spectrum disorder, and multiple congenital abnormalities. During the last 15 years, we performed CMA in approximately 8,000 patients with neurodevelopmental and/or congenital disorders, of which 13 (0.16%) genetically catastrophic complex chromosomal rearrangements were identified. These ultrarare rearrangements showed clustering of breakpoints, characteristic of chromoanagenesis events. Al1 13 complex events display underlying formation mechanisms, originating either by a synchronization of the shattering of clustered chromosome regions in which regional asynchrony of DNA replication may be one of the main causes of disruption. We provide an overview of the copy number profiling in these patients. Although several previous studies have suggested that chromoanagenesis is often a genetic disease source in postnatal diagnostic screening, due to either the challenge of clinical interpretation of these complex rearrangements or the limitation of microarray resolution relative to the small size and complexity of chromogenic induced chromosome abnormalities, bringing further attention and to study its occurrence in the clinical setting is extremely important.

Nanopore Sequencing Indicates That Tandem Amplification of Chromosome 20q11.21 in Human Pluripotent Stem Cells Is Driven by Break-Induced Replication

Copy number variants (CNVs) are genomic rearrangements implicated in numerous congenital and acquired diseases, including cancer. The appearance of culture-acquired CNVs in human pluripotent stem cells (PSCs) has prompted concerns for their use in regenerative medicine. A particular problem in PSC is the frequent occurrence of CNVs in the q11.21 region of chromosome 20. However, the exact mechanism of origin of this amplicon remains elusive due to the difficulty in delineating its sequence and breakpoints. Here, we have addressed this problem using long-read Nanopore sequencing of two examples of this CNV, present as duplication and as triplication. In both cases, the CNVs were arranged in a head-to-tail orientation, with microhomology sequences flanking or overlapping the proximal and distal breakpoints. These breakpoint signatures point to a mechanism of microhomology-mediated break-induced replication in CNV formation, with surrounding Alu sequences likely contributing to the instability of this genomic region.

Long-read genome sequencing for the molecular diagnosis of neurodevelopmental disorders

Exome and genome sequencing have proven to be effective tools for the diagnosis of neurodevelopmental disorders (NDDs), but large fractions of NDDs cannot be attributed to currently detectable genetic variation. This is likely, at least in part, a result of the fact that many genetic variants are difficult or impossible to detect through typical short-read sequencing approaches. Here, we describe a genomic analysis using Pacific Biosciences circular consensus sequencing (CCS) reads, which are both long (>10 kb) and accurate (>99% bp accuracy). We used CCS on six proband-parent trios with NDDs that were unexplained despite extensive testing, including genome sequencing with short reads. We identified variants and created de novo assemblies in each trio, with global metrics indicating these datasets are more accurate and comprehensive than those provided by short-read data. In one proband, we identified a likely pathogenic (LP), de novo L1-mediated insertion in CDKL5 that results in duplication of exon 3, leading to a frameshift. In a second proband, we identified multiple large de novo structural variants, including insertion-translocations affecting DGKB and MLLT3, which we show disrupt MLLT3 transcript levels. We consider this extensive structural variation likely pathogenic. The breadth and quality of variant detection, coupled to finding variants of clinical and research interest in two of six probands with unexplained NDDs, support the hypothesis that long-read genome sequencing can substantially improve rare disease genetic discovery rates.

Complex genomic alterations and intellectual disability: an interpretative challenge

BACKGROUND

Complex chromosomal rearrangements (CCRs) are structural rearrangements involving more than three chromosomes or having more than two breaks; approximately 70% are not associated with any clinical phenotype. Here, we describe a CCR segregating in a two-generation family.

METHOD

A 4-year-old male was evaluated for developmental delay, mild intellectual disability and epicanthus. Karyotype, fluorescence in situ hybridisation (FISH) analysis and array comparative genomic hybridisation (aCGH) analysis were performed on the patient and of all family members.

RESULT

Array CGH analysis of the proband detected two non-contiguous genomic gains of chromosome 2 at bands q32.3q33.2 and bands q36.1q36.3. Both karyotype and FISH analysis revealed a recombinant chromosome 2 with a direct insertion of regions q32.3q33.2 and q36.1q36.3 into region q12. Both of these regions were also present in their original location. Karyotype and FISH analysis of the father revealed a de novo direct insertion of regions q32.3q33.2 and q36.1q36.3 into region q12. Moreover, a de novo balanced translocation involving the q arm of the same chromosome 2 and the p arm of chromosome 10 was observed in the father of the proband. The single nucleotide polymorphism (SNP) array analysis and haplotype reconstruction confirmed the paternal origin of the duplications. Karyotype, FISH analysis and array CGH analysis of other family members were all normal.

CONCLUSION

This report underlines the importance of using different methods to correctly evaluate the origin and the structure of CCRs in order to provide an appropriate management of the patients and a good estimation of the reproductive risk of the family.

Chromosome structural variation in tumorigenesis: mechanisms of formation and carcinogenesis

With the rapid development of next-generation sequencing technology, chromosome structural variation has gradually gained increased clinical significance in tumorigenesis. However, the molecular mechanism(s) underlying this structural variation remain poorly understood. A search of the literature shows that a three-dimensional chromatin state plays a vital role in inducing structural variation and in the gene expression profiles in tumorigenesis. Structural variants may result in changes in copy number or deletions of coding sequences, as well as the perturbation of structural chromatin features, especially topological domains, and disruption of interactions between genes and their regulatory elements. This review focuses recent work aiming at elucidating how structural variations develop and misregulate oncogenes and tumor suppressors, to provide general insights into tumor formation mechanisms and to provide potential targets for future anticancer therapies.

Breakpoint junction analysis for complex genomic rearrangements with the caldera volcano-like pattern

Chromosomal triplications can be classified into recurrent and nonrecurrent triplications. Most of the nonrecurrent triplications are embedded in duplicated segments, and duplication-inverted triplication-duplication (DUP-TRP/INV-DUP) has been established as one of the mechanisms of triplication. This study aimed to reveal the underlying mechanism of the TRP-DUP-TRP pattern of chromosomal aberrations, in which the appearance of moving averages obtained through array-based comparative genomic hybridization analysis is similar to the shadows of the caldera volcano-like pattern, which were first identified in two patients with neurodevelopmental disabilities. For this purpose, whole-genome sequencing using long-read Nanopore sequencing was carried out to confirm breakpoint junctions. Custom array analysis and Sanger sequencing were also used to detect all breakpoint junctions. As a result, the TRP-DUP-TRP pattern consisted of only two patterns of breakpoint junctions in both patients. In patient 1, microhomologies were identified in breakpoint junctions. In patient 2, more complex architectures with insertional segments were identified. Thus, replication-based mechanisms were considered as a mechanism of the TRP-DUP-TRP pattern.

Replication stress at microsatellites causes DNA double-strand breaks and break-induced replication

Short tandemly repeated DNA sequences, termed microsatellites, are abundant in the human genome. These microsatellites exhibit length instability and susceptibility to DNA double-strand breaks (DSBs) due to their tendency to form stable non-B DNA structures. Replication-dependent microsatellite DSBs are linked to genome instability signatures in human developmental diseases and cancers. To probe the causes and consequences of microsatellite DSBs, we designed a dual-fluorescence reporter system to detect DSBs at expanded (CTG/CAG) _n and polypurine/polypyrimidine (Pu/Py) mirror repeat structures alongside the c-myc replication origin integrated at a single ectopic chromosomal site. Restriction cleavage near the (CTG/CAG)₁₀₀ microsatellite leads to homology-directed single-strand annealing between flanking AluY elements and reporter gene deletion that can be detected by flow cytometry. However, in the absence of restriction cleavage, endogenous and exogenous replication stressors induce DSBs at the (CTG/CAG)₁₀₀ and Pu/Py microsatellites. DSBs map to a narrow region at the downstream edge of the (CTG)₁₀₀ lagging-strand template. (CTG/CAG) _n chromosome fragility is repeat length-dependent, whereas instability at the (Pu/Py) microsatellites depends on replication polarity. Strikingly, restriction-generated DSBs and replication-dependent DSBs are not repaired by the same mechanism. Knockdown of DNA damage response proteins increases (Rad18, polymerase (Pol) η, Pol κ) or decreases (Mus81) the sensitivity of the (CTG/CAG)₁₀₀ microsatellites to replication stress. Replication stress and DSBs at the ectopic (CTG/CAG)₁₀₀ microsatellite lead to break-induced replication and high-frequency mutagenesis at a flanking thymidine kinase gene. Our results show that non-B structure-prone microsatellites are susceptible to replication-dependent DSBs that cause genome instability.

Whole Genome Low-Coverage Sequencing Concurrently Detecting Copy Number Variations and Their Underlying Complex Chromosomal Rearrangements by Systematic Breakpoint Mapping in Intellectual Deficiency/Developmental Delay Patients

Simple copy number variations (CNVs) detected by chromosomal microarray (CMA) can result from complex structural changes. Therefore, it is necessary to characterize potential structural changes that cause pathogenic CNVs. We applied whole-genome low-coverage sequencing (WGLCS) to concurrently detect pathogenic CNVs and their associated chromosomal rearrangements in 15 patients. All the patients had an average of 2–3 pathogenic CNVs involving 1–2 chromosomes. WGLCS identified all the 34 pathogenic CNVs found by microarray. By identifying chimeric read pairs, WGLCS mapped 70 breakpoints in these patients, of which 47 were finely mapped at the nucleotide level and confirmed by subsequent PCR amplification and Sanger sequencing of the junction fragments. In 15 patients, structural rearrangements were defined at molecular level in 13 patients. In 13 patients, WGLCS reveal no additional results in two patients. In another 11 patients, WGLCS revealed new breakpoints or finely mapped the genes disrupted by breakpoints or 1–6 bp microhomology and/or short insertion (4–70 bp) in the breakpoints junctions. However, structural changes in the other two patients still remained unclear after WGLCS was performed. The structural alteration identified in the 13 patients could be divided into the following categories: (1) interstitial inverted duplication with concomitant terminal deletion (inv dup del) (P1,P4,P9,P11); (2) the product of pericentric inversion (P5); (3) ring chromosome (P8); (4) interstitial duplication and/or triplication (P6, P7); and (5) +der(22)t(11;22) (P2,P15); (6) complex structural rearrangements (P3,P12,P14). WGLCS displayed the ability to discover CNVs and define breakpoints and its disrupted genes and its surrounding sequences in one experiment at base-pair-resolution, which help us to learn more about the mechanisms of formation of observed genomic rearrangements, and in which DNA replicative/repair mechanism might contribute to the formation of complex rearrangements in 11 patients. Clear karyotype at molecular level could help provide an accurate evaluation of recurrent risk and guide prenatal diagnosis or reproductive planning.

The involvement of U-type dicentric chromosomes in the formation of terminal deletions with or without adjacent inverted duplications

An inverted duplication with a terminal deletion (inv-dup-del) is one of the complex constitutional structural rearrangements that can occur in a chromosome. Although breakages of dicentric chromosome have been suggested, the precise mechanism of this is yet to be fully understood. In our present study, we investigated the genomic structure of 10 inv-dup-del cases to elucidate this mechanism. Two recurrent 8p inv-dup-del cases harbored a large copy-number-neutral region between the duplication and deletion in common. Although the other non-recurrent cases did not appear to have this copy-number-neutral region, refined sequencing analysis identified that they contained a small intervening region at the junction between the inverted and non-inverted segment. The size of this small intervening region ranged from 1741 to 3728 bp. Combined with a presence of microhomology at the junction, a resolution of the replication fork stalling through template switching within the same replication fork is suggested. We further observed two cases with mosaicism of the dicentric chromosome and various structural rearrangements related to the dicentric chromosome. Refined analysis allowed us to identify different breakpoints on the same chromosome in the same case, implicating multiple rounds of U-type formation and its breakage. From these results, we propose that a replication-based mechanism generates unstable dicentric chromosomes and that their breakage leads to the formation of inv-dup-dels and other related derivative chromosomes.

Whole-genome sequencing reveals complex chromosome rearrangement disrupting NIPBL in infant with Cornelia de Lange syndrome

Clinical laboratory diagnostic evaluation of the genomes of children with suspected genetic disorders, including chromosomal microarray and exome sequencing, cannot detect copy number neutral genomic rearrangements such as inversions, balanced translocations, and complex chromosomal rearrangements (CCRs). We describe an infant with a clinical diagnosis of Cornelia de Lange syndrome (CdLS) in whom chromosome analysis revealed a de novo complex balanced translocation, 46,XY,t(5;7;6)(q11.2;q32;q13)dn. Subsequent molecular characterization by whole-genome sequencing (WGS) identified 23 breakpoints, delineating segments derived from four chromosomes (5;6;7;21) in ancestral or inverted orientation. One of the breakpoints disrupted a known CdLS gene, NIPBL. Further investigation revealed paternal origin of the CCR allele, clustering of the breakpoint junctions, and molecular repair signatures suggestive of a single catastrophic event. Notably, very short DNA segments (25 and 41 bp) were included in the reassembled chromosomes, lending additional support that the DNA repair machinery can detect and repair such segments. Interestingly, there was an independent paternally derived miniscule complex rearrangement, possibly predisposing to subsequent genomic instability. In conclusion, we report a CCR causing a monogenic Mendelian disorder, urging WGS analysis of similar unsolved cases with suspected Mendelian disorders. Breakpoint analysis allowed for identification of the underlying molecular diagnosis and implicated chromoanagenesis in CCR formation.

Structural variant identification and characterization

Structural variant (SV) differences between human genomes can cause germline and mosaic disease as well as inter-individual variation. De-regulation of accurate DNA repair and genomic surveillance mechanisms results in a large number of SVs in cancer. Analysis of the DNA sequences at SV breakpoints can help identify pathways of mutagenesis and regions of the genome that are more susceptible to rearrangement. Large-scale SV analyses have been enabled by high-throughput genome-level sequencing on humans in the past decade. These studies have shed light on the mechanisms and prevalence of complex genomic rearrangements. Recent advancements in both sequencing and other mapping technologies as well as calling algorithms for detection of genomic rearrangements have helped propel SV detection into population-scale studies, and have begun to elucidate previously inaccessible regions of the genome. Here, we discuss the genomic organization of simple and complex SVs, the molecular mechanisms of their formation, and various ways to detect them. We also introduce methods for characterizing SVs and their consequences on human genomes.

The rare Alus element-mediated chimerism of multiple de novo complex rearrangement sequences in GAN result in giant axonal neuropathy

Giant axonal neuropathy (GAN) is a rare and grievous autosomal recessive neurodegenerative disease due to loss-of-function mutation in GAN. However, the chimerism of complex rearrangement sequences of GAN has not been reported so far, and the mechanism for its complex rearrangements remains to be determined. We identified a family with clinical symptoms of GAN and aimed to reveal a genetic cause underlying this disease. By whole-exome sequencing in the patient we identified a novel homozygous frameshift mutation with 1 bp deletion (c.27delC) in GAN. However, when analyzed the patient's genomic DNA (gDNA) by quantitative real-time PCR and breakpoint DNA sequencing, we found the chimerism of multiple complex rearrangement sequences encompassing exon 1 of GAN in the patient's genome. The microhomology and localization of the breakpoint indicated that they may be caused by Alu repeat elements. We also found that the mRNA expression level of GAN in patient's lymphocyte was decreased, confirming the pathogenicity of these mutations. Our study is the first reported on many complex rearrangement sequences mosaic in GAN mediated by Alu element. The patient here is not a simple homozygous frameshift mutation, but a compound heterozygous paternal c.27delC mutation and the chimerism of multiple de novo complex rearrangement sequences in GAN. Our results may also provide new insights into the formation and pathogenicity of complex rearrangement in GAN, and may be helpful to genetic counseling and genetic testing. It also enriches the Alu-mediated disease-associated database which are important for correct clinical interpretation.

Distinct patterns of complex rearrangements and a mutational signature of microhomeology are frequently observed in PLP1 copy number gain structural variants

BACKGROUND

We investigated the features of the genomic rearrangements in a cohort of 50 male individuals with proteolipid protein 1 (PLP1) copy number gain events who were ascertained with Pelizaeus-Merzbacher disease (PMD; MIM: 312080). We then compared our new data to previous structural variant mutagenesis studies involving the Xq22 region of the human genome. The aggregate data from 159 sequenced join-points (discontinuous sequences in the reference genome that are joined during the rearrangement process) were studied. Analysis of these data from 150 individuals enabled the spectrum and relative distribution of the underlying genomic mutational signatures to be delineated.

METHODS

Genomic rearrangements in PMD individuals with PLP1 copy number gain events were investigated by high-density customized array or clinical chromosomal microarray analysis and breakpoint junction sequence analysis.

RESULTS

High-density customized array showed that the majority of cases (33/50; ~ 66%) present with single duplications, although complex genomic rearrangements (CGRs) are also frequent (17/50; ~ 34%). Breakpoint mapping to nucleotide resolution revealed further previously unknown structural and sequence complexities, even in single duplications. Meta-analysis of all studied rearrangements that occur at the PLP1 locus showed that single duplications were found in ~ 54% of individuals and that, among all CGR cases, triplication flanked by duplications is the most frequent CGR array CGH pattern observed. Importantly, in ~ 32% of join-points, there is evidence for a mutational signature of microhomeology (highly similar yet imperfect sequence matches).

CONCLUSIONS

These data reveal a high frequency of CGRs at the PLP1 locus and support the assertion that replication-based mechanisms are prominent contributors to the formation of CGRs at Xq22. We propose that microhomeology can facilitate template switching, by stabilizing strand annealing of the primer using W-C base complementarity, and is a mutational signature for replicative repair.

A 14q distal chromoanagenesis elucidated by whole genome sequencing

Chromoanagenesis represents an extreme form of genomic rearrangements involving multiple breaks occurring on a single or multiple chromosomes. It has been recently described in both acquired and rare constitutional genetic disorders. Constitutional chromoanagenesis events could lead to abnormal phenotypes including developmental delay and congenital anomalies, and have also been implicated in some specific syndromic disorders. We report the case of a girl presenting with growth retardation, hypotonia, microcephaly, dysmorphic features, coloboma, and hypoplastic corpus callosum. Karyotype showed a de novo structurally abnormal chromosome 14q31qter region. Molecular characterization using SNP-array revealed a complex unbalanced rearrangement in 14q31.1-q32.2, on the paternal chromosome 14, including thirteen interstitial deletions ranging from 33 kb to 1.56 Mb in size, with a total of 4.1 Mb in size, thus suggesting that a single event like chromoanagenesis occurred. To our knowledge, this is one of the first case of 14q distal deletion due to a germline chromoanagenesis. Genome sequencing allowed the characterization of 50 breakpoints, leading to interruption of 10 genes including YY1 which fit with the patient's phenotype. This precise genotyping of breaking junction allowed better definition of genotype-phenotype correlations.

Megabase Length Hypermutation Accompanies Human Structural Variation at 17p11.2

DNA rearrangements resulting in human genome structural variants (SVs) are caused by diverse mutational mechanisms. We used long- and short-read sequencing technologies to investigate end products of de novo chromosome 17p11.2 rearrangements and query the molecular mechanisms underlying both recurrent and non-recurrent events. Evidence for an increased rate of clustered single-nucleotide variant (SNV) mutation in cis with non-recurrent rearrangements was found. Indel and SNV formation are associated with both copy-number gains and losses of 17p11.2, occur up to ∼1 Mb away from the breakpoint junctions, and favor C > G transversion substitutions; results suggest that single-stranded DNA is formed during the genesis of the SV and provide compelling support for a microhomology-mediated break-induced replication (MMBIR) mechanism for SV formation. Our data show an additional mutational burden of MMBIR consisting of hypermutation confined to the locus and manifesting as SNVs and indels predominantly within genes.

Characterization and evolutionary dynamics of complex regions in eukaryotic genomes

Complex regions in eukaryotic genomes are typically characterized by duplications of chromosomal stretches that often include one or more genes repeated in a tandem array or in relatively close proximity. Nevertheless, the repetitive nature of these regions, together with the often high sequence identity among repeats, have made complex regions particularly recalcitrant to proper molecular characterization, often being misassembled or completely absent in genome assemblies. This limitation has prevented accurate functional and evolutionary analyses of these regions. This is becoming increasingly relevant as evidence continues to support a central role for complex genomic regions in explaining human disease, developmental innovations, and ecological adaptations across phyla. With the advent of long-read sequencing technologies and suitable assemblers, the development of algorithms that can accommodate sample heterozygosity, and the adoption of a pangenomic-like view of these regions, accurate reconstructions of complex regions are now within reach. These reconstructions will finally allow for accurate functional and evolutionary studies of complex genomic regions, underlying the generation of genotype-phenotype maps of unprecedented resolution.

Chromoanagenesis: cataclysms behind complex chromosomal rearrangements

Background

During the last decade, genome sequencing projects in cancer genomes as well as in patients with congenital diseases and healthy individuals have led to the identification of new types of massive chromosomal rearrangements arising during single chaotic cellular events. These unanticipated catastrophic phenomenon are termed chromothripsis, chromoanasynthesis and chromoplexis., and are grouped under the name of “chromoanagenesis”.

Results

For each process, several specific features have been described, allowing each phenomenon to be distinguished from each other and to understand its mechanism of formation and to better understand its aetiology. Thus, chromothripsis derives from chromosome shattering followed by the random restitching of chromosomal fragments with low copy-number change whereas chromoanasynthesis results from erroneous DNA replication of a chromosome through serial fork stalling and template switching with variable copy-number gains, and chromoplexy refers to the occurrence of multiple inter-and intra-chromosomal translocations and deletions with little or no copy-number alterations in prostate cancer. Cumulating data and experimental models have shown that chromothripsis and chromoanasynthesis may essentially result from lagging chromosome encapsulated in micronuclei or telomere attrition and end-to-end telomere fusion.

Conclusion

The concept of chromanagenesis has provided new insight into the aetiology of complex structural rearrangements, the connection between defective cell cycle progression and genomic instability, and the complexity of cancer evolution. Increasing reported chromoanagenesis events suggest that these chaotic mechanisms are probably much more frequent than anticipated.

Exonic duplication of the OTC gene by a complex rearrangement that likely occurred via a replication-based mechanism: a case report

BACKGROUND

Ornithine transcarbamylase deficiency (OTCD) is an X-linked recessive disorder involving a defect in the urea cycle caused by OTC gene mutations. Although a total of 417 disease-causing mutations in OTC have been reported, structural abnormalities in this gene are rare. We here describe a female OTCD case caused by an exonic duplication of the OTC gene (exons 1-6).

CASE PRESENTATION

A 23-year-old woman with late-onset OTCD diagnosed by biochemical testing was subjected to subsequent genetic testing. Sanger sequencing revealed no pathogenic mutation throughout the coding exons of the OTC gene, but multiplex ligation-dependent probe amplification (MLPA) revealed duplication of exons 1-6. Further genetic analyses revealed an inversion of duplicated exon 1 and a tandem duplication of exons 2-6. Each of the junctions of the inversion harbored a microhomology and non-templated microinsertion, respectively, suggesting a replication-based mechanism. The duplication was also of de novo origin but segregation analysis indicated that it took place in the paternal chromosome.

CONCLUSION

We report the first OTCD case harboring an exonic duplication in the OTC gene. The functional defects caused by this anomaly were determined via structural analysis of its complex rearrangements.

RNA-mediated gene fusion in mammalian cells

One of the hallmarks of cancer is the formation of oncogenic fusion genes as a result of chromosomal translocations. Fusion genes are presumed to form before fusion RNA expression. However, studies have reported the presence of fusion RNAs in individuals who were negative for chromosomal translocations. These observations give rise to "the cart before the horse" hypothesis, in which the genesis of a fusion RNA precedes the fusion gene. The fusion RNA then guides the genomic rearrangements that ultimately result in a gene fusion. However, RNA-mediated genomic rearrangements in mammalian cells have never been demonstrated. Here we provide evidence that expression of a chimeric RNA drives formation of a specified gene fusion via genomic rearrangement in mammalian cells. The process is: (i) specified by the sequence of chimeric RNA involved, (ii) facilitated by physiological hormone levels, (iii) permissible regardless of intrachromosomal (TMPRSS2-ERG) or interchromosomal (TMPRSS2-ETV1) fusion, and (iv) can occur in normal cells before malignant transformation. We demonstrate that, contrary to "the cart before the horse" model, it is the antisense rather than sense chimeric RNAs that effectively drive gene fusion, and that this disparity can be explained by transcriptional conflict. Furthermore, we identified an endogenous RNA AZI1 that functions as the "initiator" RNA to induce TMPRSS2-ERG fusion. RNA-driven gene fusion demonstrated in this report provides important insight in early disease mechanisms, and could have fundamental implications in the biology of mammalian genome stability, as well as gene-editing technology via mechanisms native to mammalian cells.

Flanking complex copy number variants in the same family formed through unequal crossing-over during meiosis

Two phenomena that have been described in germline complex genomic rearrangements (CGRs) formation are chromothripsis and chromoanasynthesis, characterized by distinct features such as the orientation and copy number of the involved fragments. Herein we present different CGRs on chromosome 5p in a mother and her daughter that through unequal crossing-over during meiosis has evolved from a chromothriptic rearrangement in the mother into another complex rearrangement in her daughter involving both deletions and duplications. Initially, both rearrangements were classified as simple copy number variants, but follow-up studies using whole-genome sequencing revealed a much more complex nature of both rearrangements and enabled us to decipher the biological process involved in the formation of the rearrangement found in the daughter. In conclusion, these two cases highlight the need of analyzing the inheritance patterns of CGRs, and provide an example of a disease-causing CGR formed through multiple genetic events.

The 3D genome: From fundamental principles to disease and cancer

In higher eukaryotes, the three-dimensional (3D) organization of the genome is intimately related to numerous key biological functions including gene expression, DNA repair and DNA replication regulations. Alteration of this 3D organization is detrimental to the organism and can give rise to a broad range of diseases such as cancers. Here, we review recent advances in the field. We first describe how the genome is packed in 3D to form chromosome territories, compartments and domains. We also give an overview of the recent techniques that allow to map the genome in 3D up to the kilobase resolution. We then discuss potential mechanisms by which genome misfolding can affect proper gene expression by distal enhancers, and how the 3D genome influences the formation of genomic rearrangements.

Novel applications of array comparative genomic hybridization in molecular diagnostics

INTRODUCTION

In 2004, the implementation of array comparative genomic hybridization (array comparative genome hybridization [CGH]) into clinical practice marked a new milestone for genetic diagnosis. Array CGH and single-nucleotide polymorphism (SNP) arrays enable genome-wide detection of copy number changes in a high resolution, and therefore microarray has been recognized as the first-tier test for patients with intellectual disability or multiple congenital anomalies, and has also been applied prenatally for detection of clinically relevant copy number variations in the fetus. Area covered: In this review, the authors summarize the evolution of array CGH technology from their diagnostic laboratory, highlighting exonic SNP arrays developed in the past decade which detect small intragenic copy number changes as well as large DNA segments for the region of heterozygosity. The applications of array CGH to human diseases with different modes of inheritance with the emphasis on autosomal recessive disorders are discussed. Expert commentary: An exonic array is a powerful and most efficient clinical tool in detecting genome wide small copy number variants in both dominant and recessive disorders. However, whole-genome sequencing may become the single integrated platform for detection of copy number changes, single-nucleotide changes as well as balanced chromosomal rearrangements in the near future.