Characterization of the past and current duplication activities in the human 22q11.2 region

Resolving the 22q11.2 deletion using CTLR-Seq reveals chromosomal rearrangement mechanisms and individual variance in breakpoints

We developed a generally applicable method, CRISPR/Cas9-targeted long-read sequencing (CTLR-Seq), to resolve, haplotype-specifically, the large and complex regions in the human genome that had been previously impenetrable to sequencing analysis, such as large segmental duplications (SegDups) and their associated genome rearrangements. CTLR-Seq combines in vitro Cas9-mediated cutting of the genome and pulse-field gel electrophoresis to isolate intact large (i.e., up to 2,000 kb) genomic regions that encompass previously unresolvable genomic sequences. These targets are then sequenced (amplification-free) at high on-target coverage using long-read sequencing, allowing for their complete sequence assembly. We applied CTLR-Seq to the SegDup-mediated rearrangements that constitute the boundaries of, and give rise to, the 22q11.2 Deletion Syndrome (22q11DS), the most common human microdeletion disorder. We then performed de novo assembly to resolve, at base-pair resolution, the full sequence rearrangements and exact chromosomal breakpoints of 22q11.2DS (including all common subtypes). Across multiple patients, we found a high degree of variability for both the rearranged SegDup sequences and the exact chromosomal breakpoint locations, which coincide with various transposons within the 22q11.2 SegDups, suggesting that 22q11DS can be driven by transposon-mediated genome recombination. Guided by CTLR-Seq results from two 22q11DS patients, we performed three-dimensional chromosomal folding analysis for the 22q11.2 SegDups from patient-derived neurons and astrocytes and found chromosome interactions anchored within the SegDups to be both cell type-specific and patient-specific. Lastly, we demonstrated that CTLR-Seq enables cell-type specific analysis of DNA methylation patterns within the deletion haplotype of 22q11DS.

22q11.2 Low Copy Repeats Expanded in the Human Lineage

Segmental duplications or low copy repeats (LCRs) constitute duplicated regions interspersed in the human genome, currently neglected in standard analyses due to their extreme complexity. Recent functional studies have indicated the potential of genes within LCRs in synaptogenesis, neuronal migration, and neocortical expansion in the human lineage. One of the regions with the highest proportion of duplicated sequence is the 22q11.2 locus, carrying eight LCRs (LCR22-A until LCR22-H), and rearrangements between them cause the 22q11.2 deletion syndrome. The LCR22-A block was recently reported to be hypervariable in the human population. It remains unknown whether this variability also exists in non-human primates, since research is strongly hampered by the presence of sequence gaps in the human and non-human primate reference genomes. To chart the LCR22 haplotypes and the associated inter- and intra-species variability, we de novo assembled the region in non-human primates by a combination of optical mapping techniques. A minimal and likely ancient haplotype is present in the chimpanzee, bonobo, and rhesus monkey without intra-species variation. In addition, the optical maps identified assembly errors and closed gaps in the orthologous chromosome 22 reference sequences. These findings indicate the LCR22 expansion to be unique to the human population, which might indicate involvement of the region in human evolution and adaptation. Those maps will enable LCR22-specific functional studies and investigate potential associations with the phenotypic variability in the 22q11.2 deletion syndrome.

The Transposable Element Environment of Human Genes Differs According to Their Duplication Status and Essentiality

Transposable elements (TEs) are major components of eukaryotic genomes and represent approximately 45% of the human genome. TEs can be important sources of novelty in genomes and there is increasing evidence that TEs contribute to the evolution of gene regulation in mammals. Gene duplication is an evolutionary mechanism that also provides new genetic material and opportunities to acquire new functions. To investigate how duplicated genes are maintained in genomes, here, we explored the TE environment of duplicated and singleton genes. We found that singleton genes have more short-interspersed nuclear elements and DNA transposons in their vicinity than duplicated genes, whereas long-interspersed nuclear elements and long-terminal repeat retrotransposons have accumulated more near duplicated genes. We also discovered that this result is highly associated with the degree of essentiality of the genes with an unexpected accumulation of short-interspersed nuclear elements and DNA transposons around the more-essential genes. Our results underline the importance of taking into account the TE environment of genes to better understand how duplicated genes are maintained in genomes.

Formation of human long intergenic non-coding RNA genes, pseudogenes, and protein genes: Ancestral sequences are key players

Pathways leading to formation of non-coding RNA and protein genes are varied and complex. We report finding a conserved repeat sequence present in human and chimpanzee genomes that appears to have originated from a common primate ancestor. This sequence is repeatedly copied in human chromosome 22 (chr22) low copy repeats (LCR22) or segmental duplications and forms twenty-one different genes, which include the human long intergenic non-coding RNA (lincRNA) family FAM230, a newly discovered lincRNA gene family termed conserved long intergenic non-coding RNAs (clincRNA), pseudogene families, as well as the gamma-glutamyltransferase (GGT) protein gene family and the RNA pseudogenes that originate from GGT sequences. Of particular interest are the GGT5 and USP18 protein genes that appear to have formed from an homologous repeat sequence that also forms the clincRNA gene family. The data point to ancestral DNA sequences, conserved through evolution and duplicated in humans by chromosomal repeat sequences that may serve as functional genomic elements in the development of diverse genes.

Atypical chromosome 22q11.2 deletions are complex rearrangements and have different mechanistic origins

The majority (99%) of individuals with 22q11.2 deletion syndrome (22q11.2DS) have a deletion that is caused by non-allelic homologous recombination between two of four low copy repeat clusters on chromosome 22q11.2 (LCR22s). However, in a small subset of patients, atypical deletions are observed with at least one deletion breakpoint within unique sequence between the LCR22s. The position of the chromosome breakpoints and the mechanisms driving those atypical deletions remain poorly studied. Our large-scale, whole genome sequencing study of >1500 subjects with 22q11.2DS identified six unrelated individuals with atypical deletions of different types. Using a combination of whole genome sequencing data and fiber-fluorescence in situ hybridization, we mapped the rearranged alleles in these subjects. In four of them, the distal breakpoints mapped within one of the LCR22s and we found that the deletions likely occurred by replication-based mechanisms. Interestingly, in two of them, an inversion probably preceded inter-chromosomal 'allelic' homologous recombination between differently oriented LCR22-D alleles. Inversion associated allelic homologous recombination (AHR) may well be a common mechanism driving (atypical) deletions on 22q11.2.

Molecular genetics of 22q11.2 deletion syndrome

The 22q11.2 deletion syndrome (22q11.2DS) is a congenital malformation and neuropsychiatric disorder caused by meiotic chromosome rearrangements. One of the goals of this review is to summarize the current state of basic research studies of 22q11.2DS. It highlights efforts to understand the mechanisms responsible for the 22q11.2 deletion that occurs in meiosis. This mechanism involves the four sets of low copy repeats (LCR22) that are dispersed in the 22q11.2 region and the deletion is mediated by nonallelic homologous recombination events. This review also highlights selected genes mapping to the 22q11.2 region that may contribute to the typical clinical findings associated with the disorder and explain that mutations in genes on the remaining allele can uncover rare recessive conditions. Another important aspect of 22q11.2DS is the existence of phenotypic heterogeneity. While some patients are mildly affected, others have severe medical, cognitive, and/or psychiatric challenges. Variability may be due in part to the presence of genetic modifiers. This review discusses current genome-wide efforts to identify such modifiers that could shed light on molecular pathways required for normal human development, cognition or behavior.

The 22q11 low copy repeats are characterized by unprecedented size and structural variability

Low copy repeats (LCRs) are recognized as a significant source of genomic instability, driving genome variability and evolution. The Chromosome 22 LCRs (LCR22s) mediate nonallelic homologous recombination (NAHR) leading to the 22q11 deletion syndrome (22q11DS). However, LCR22s are among the most complex regions in the genome, and their structure remains unresolved. The difficulty in generating accurate maps of LCR22s has also hindered localization of the deletion end points in 22q11DS patients. Using fiber FISH and Bionano optical mapping, we assembled LCR22 alleles in 187 cell lines. Our analysis uncovered an unprecedented level of variation in LCR22s, including LCR22A alleles ranging in size from 250 to 2000 kb. Further, the incidence of various LCR22 alleles varied within different populations. Additionally, the analysis of LCR22s in 22q11DS patients and their parents enabled further refinement of the rearrangement site within LCR22A and -D, which flank the 22q11 deletion. The NAHR site was localized to a 160-kb paralog shared between the LCR22A and -D in seven 22q11DS patients. Thus, we present the most comprehensive map of LCR22 variation to date. This will greatly facilitate the investigation of the role of LCR variation as a driver of 22q11 rearrangements and the phenotypic variability among 22q11DS patients.

Formation of a Family of Long Intergenic Noncoding RNA Genes with an Embedded Translocation Breakpoint Motif in Human Chromosomal Low Copy Repeats of 22q11.2-Some Surprises and Questions

A family of long intergenic noncoding RNA (lincRNA) genes, FAM230 is formed via gene sequence duplication, specifically in human chromosomal low copy repeats (LCR) or segmental duplications. This is the first group of lincRNA genes known to be formed by segmental duplications and is consistent with current views of evolution and the creation of new genes via DNA low copy repeats. It appears to be an efficient way to form multiple lincRNA genes. But as these genes are in a critical chromosomal region with respect to the incidence of abnormal translocations and resulting genetic abnormalities, the 22q11.2 region, and also carry a translocation breakpoint motif, several intriguing questions arise concerning the presence and function of the translocation breakpoint sequence in RNA genes situated in LCR22s.

A family of long intergenic non-coding RNA genes in human chromosomal region 22q11.2 carry a DNA translocation breakpoint/AT-rich sequence

FAM230C, a long intergenic non-coding RNA (lincRNA) gene in human chromosome 13 (chr13) is a member of lincRNA genes termed family with sequence similarity 230. An analysis using bioinformatics search tools and alignment programs was undertaken to determine properties of FAM230C and its related genes. Results reveal that the DNA translocation element, the Translocation Breakpoint Type A (TBTA) sequence, which consists of satellite DNA, Alu elements, and AT-rich sequences is embedded in the FAM230C gene. Eight lincRNA genes related to FAM230C also carry the TBTA sequences. These genes were formed from a large segment of the 3’ half of the FAM230C sequence duplicated in chr22, and are specifically in regions of low copy repeats (LCR22)s, in or close to the 22q.11.2 region. 22q11.2 is a chromosomal segment that undergoes a high rate of DNA translocation and is prone to genetic deletions. FAM230C-related genes present in other chromosomes do not carry the TBTA motif and were formed from the 5’ half region of the FAM230C sequence. These findings identify a high specificity in lincRNA gene formation by gene sequence duplication in different chromosomes.

Nested Inversion Polymorphisms Predispose Chromosome 22q11.2 to Meiotic Rearrangements

Inversion polymorphisms between low-copy repeats (LCRs) might predispose chromosomes to meiotic non-allelic homologous recombination (NAHR) events and thus lead to genomic disorders. However, for the 22q11.2 deletion syndrome (22q11.2DS), the most common genomic disorder, no such inversions have been uncovered as of yet. Using fiber-FISH, we demonstrate that parents transmitting the de novo 3 Mb LCR22A-D 22q11.2 deletion, the reciprocal duplication, and the smaller 1.5 Mb LCR22A-B 22q11.2 deletion carry inversions of LCR22B-D or LCR22C-D. Hence, the inversions predispose chromosome 22q11.2 to meiotic rearrangements and increase the individual risk for transmitting rearrangements. Interestingly, the inversions are nested or flanking rather than coinciding with the deletion or duplication sizes. This finding raises the possibility that inversions are a prerequisite not only for 22q11.2 rearrangements but also for all NAHR-mediated genomic disorders.

An exploratory study of predisposing genetic factors for DiGeorge/velocardiofacial syndrome

DiGeorge/velocardiofacial syndrome (DGS/VCFS) is a disorder caused by a 22q11.2 deletion mediated by non-allelic homologous recombination (NAHR) between low-copy repeats (LCRs). We have evaluated the role of LCR22 genomic architecture and PRDM9 variants as DGS/VCFS predisposing factors. We applied FISH using fosmid probes on chromatin fibers to analyze the number of tandem repeat blocks in LCR22 in two DGS/VCFS fathers-of-origin with proven 22q11.2 NAHR susceptibility. Results revealed copy number variations (CNVs) of L9 and K3 fosmids in these individuals compared to controls. The total number of L9 and K3 copies was also characterized using droplet digital PCR (ddPCR). Although we were unable to confirm variations, we detected an additional L9 amplicon corresponding to a pseudogene. Moreover, none of the eight DGS/VCFS parents-of-origin was heterozygote for the inv(22)(q11.2) haplotype. PRDM9 sequencing showed equivalent allelic distributions between DGS/VCFS parents-of-origin and controls, although a new PRDM9 allele (L50) was identified in one case. Our results support the hypothesis that LCR22s variations influences 22q11.2 NAHR events, however further studies are needed to confirm this association and clarify the contribution of pseudogenes and rare PDRM9 alleles to NAHR susceptibility.

Variant discovery and breakpoint region prediction for studying the human 22q11.2 deletion using BAC clone and whole genome sequencing analysis

Velo-cardio-facial syndrome/DiGeorge syndrome/22q11.2 deletion syndrome (22q11.2DS) is caused by meiotic non-allelic homologous recombination events between flanking low copy repeats termed LCR22A and LCR22D, resulting in a 3 million base pair (Mb) deletion. Due to their complex structure, large size and high sequence identity, genetic variation within LCR22s among different individuals has not been well characterized. In this study, we sequenced 13 BAC clones derived from LCR22A/D and aligned them with 15 previously available BAC sequences to create a new genetic variation map. The thousands of variants identified by this analysis were not uniformly distributed in the two LCR22s. Moreover, shared single nucleotide variants between LCR22A and LCR22D were enriched in the Breakpoint Cluster Region pseudogene (BCRP) block, suggesting the existence of a possible recombination hotspot there. Interestingly, breakpoints for atypical 22q11.2 rearrangements have previously been located to BCRPs To further explore this finding, we carried out in-depth analyses of whole genome sequence (WGS) data from two unrelated probands harbouring a de novo 3Mb 22q11.2 deletion and their normal parents. By focusing primarily on WGS reads uniquely mapped to LCR22A, using the variation map from our BAC analysis to help resolve allele ambiguity, and by performing PCR analysis, we infer that the deletion breakpoints were most likely located near or within the BCRP module. In summary, we found a high degree of sequence variation in LCR22A and LCR22D and a potential recombination breakpoint near or within the BCRP block, providing a starting point for future breakpoint mapping using additional trios.

Modeling a model: Mouse genetics, 22q11.2 Deletion Syndrome, and disorders of cortical circuit development

Understanding the developmental etiology of autistic spectrum disorders, attention deficit/hyperactivity disorder and schizophrenia remains a major challenge for establishing new diagnostic and therapeutic approaches to these common, difficult-to-treat diseases that compromise neural circuits in the cerebral cortex. One aspect of this challenge is the breadth and overlap of ASD, ADHD, and SCZ deficits; another is the complexity of mutations associated with each, and a third is the difficulty of analyzing disrupted development in at-risk or affected human fetuses. The identification of distinct genetic syndromes that include behavioral deficits similar to those in ASD, ADHC and SCZ provides a critical starting point for meeting this challenge. We summarize clinical and behavioral impairments in children and adults with one such genetic syndrome, the 22q11.2 Deletion Syndrome, routinely called 22q11DS, caused by micro-deletions of between 1.5 and 3.0 MB on human chromosome 22. Among many syndromic features, including cardiovascular and craniofacial anomalies, 22q11DS patients have a high incidence of brain structural, functional, and behavioral deficits that reflect cerebral cortical dysfunction and fall within the spectrum that defines ASD, ADHD, and SCZ. We show that developmental pathogenesis underlying this apparent genetic "model" syndrome in patients can be defined and analyzed mechanistically using genomically accurate mouse models of the deletion that causes 22q11DS. We conclude that "modeling a model", in this case 22q11DS as a model for idiopathic ASD, ADHD and SCZ, as well as other behavioral disorders like anxiety frequently seen in 22q11DS patients, in genetically engineered mice provides a foundation for understanding the causes and improving diagnosis and therapy for these disorders of cortical circuit development.

22q11.2 Microduplication in a patient with 19p13.12-13.13 deletion

BACKGROUND

The chromosome 22q11.2 region microduplication has been described in patients with variable phenotypes. Here we present a 3-month-old girl with both 22q11.2 microduplication and 19p13.12-13.13 deletion. The presence of both genomic imbalances in one patient has not been previously reported in literature.

METHODS

A routine G-banding karyotype analysis was performed using peripheral lymphocytes. Chromosome microarray analysis (CMA) was done using Affymetrix CytoScan™ HD array.

RESULTS

The result of karyotyping showed that the patient is 46,XX,t(12;19)(q24.3;p13.1), but CMA detected a 2.8Mb microduplication within the region 22q11.2 (chr22: 18,648,866-21,465,659) and a 1.2Mb deletion on the chromosome 19at band p13.12-p13.13 (chr19: 13,107,938-14,337,347) in her genome, while no abnormalities were identified on 12q24.3. The 3-month-old girl presented with microcephaly, cleft palate, low set and retroverted ears, and facial dysmorphism which consisted of the following: a long narrow face, widely spaced eyes, downslanting palpebral fissures, broad nasal base, short philtrum, thin upper lip, and micro/retrognathia. She also had a congenital right pulmonary artery sling and tracheal stenosis and suffered from significant hypotonia and partial bilateral mixed hearing loss.

CONCLUSIONS

We report a case of 22q11.2 duplication syndrome with 19p13.12-13.13 deletion. Synergistic effect from the two genomic imbalances is likely responsible for the complicated clinical features observed in this patient.

The role of retrotransposons in gene family expansions: insights from the mouse Abp gene family

Background

Retrotransposons have been suggested to provide a substrate for non-allelic homologous recombination (NAHR) and thereby promote gene family expansion. Their precise role, however, is controversial. Here we ask whether retrotransposons contributed to the recent expansions of the Androgen-binding protein (Abp) gene families that occurred independently in the mouse and rat genomes.

Results

Using dot plot analysis, we found that the most recent duplication in the Abp region of the mouse genome is flanked by L1Md_T elements. Analysis of the sequence of these elements revealed breakpoints that are the relicts of the recombination that caused the duplication, confirming that the duplication arose as a result of NAHR using L1 elements as substrates. L1 and ERVII retrotransposons are considerably denser in the Abp regions than in one Mb flanking regions, while other repeat types are depleted in the Abp regions compared to flanking regions. L1 retrotransposons preferentially accumulated in the Abp gene regions after lineage separation and roughly followed the pattern of Abp gene expansion. By contrast, the proportion of shared vs. lineage-specific ERVII repeats in the Abp region resembles the rest of the genome.

Conclusions

We confirmed the role of L1 repeats in Abp gene duplication with the identification of recombinant L1Md_T elements at the edges of the most recent mouse Abp gene duplication. High densities of L1 and ERVII repeats were found in the Abp gene region with abrupt transitions at the region boundaries, suggesting that their higher densities are tightly associated with Abp gene duplication. We observed that the major accumulation of L1 elements occurred after the split of the mouse and rat lineages and that there is a striking overlap between the timing of L1 accumulation and expansion of the Abp gene family in the mouse genome. Establishing a link between the accumulation of L1 elements and the expansion of the Abp gene family and identification of an NAHR-related breakpoint in the most recent duplication are the main contributions of our study.

22q11.2 Microduplication syndrome and epilepsy with continuous spikes and waves during sleep (CSWS). A case report and review of the literature

Chromosome 22q11.2 microduplication syndrome is characterized by a variable and usually mild phenotype and by incomplete penetrance. Neurological features of the syndrome may entail intellectual or learning disability, motor delay, and other neurodevelopmental disorders. However, seizures or abnormal EEG are reported in a few cases. We describe a 6-year-old girl with microduplication of chromosome 22q11.2 and epilepsy with continuous spikes and waves during sleep (CSWS). Her behavioral disorder, characterized by hyperactivity, impulsiveness, attention deficit, and aggressiveness, became progressively evident a few months after epilepsy onset, suggesting a link with the interictal epileptic activity characterizing CSWS. We hypothesize that, at least in some cases, the neurodevelopmental deficit seen in the 22q11.2 microduplication syndrome could be the consequence of a disorder of cerebral electrogenesis, suggesting the need for an EEG recording in affected individuals. Moreover, an array-CGH analysis should be performed in all individuals with cryptogenic epilepsy and CSWS.

High-resolution fish on DNA fibers for low-copy repeats genome architecture studies

Low-copy repeats (LCRs) constitute 5% of the human genome. LCRs act as substrates for non-allelic homologous recombination (NAHR) leading to genomic structural variation. The aim of this study was to assess the potential of Fiber-FISH for LCRs direct visualization to support investigations of genome architecture within these challenging genomic regions. We describe a set of Fiber-FISH experiments designed for the study of the LCR22-2. This LCR is involved in recurrent reorganizations causing different genomic disorders. Four fosmid clones covering the entire length of the LCR22-2 and two single-copy BAC-clones, delimiting the LCR22-2 proximally and distally, were selected. The probes were hybridized in different multiple color combinations on DNA fibers from two karyotypically normal cell lines. We were able to identify three distinct structural haplotypes characterized by differences in copy-number and arrangement of the LCR22-2 genes and pseudogenes. Our results show that Multicolor Fiber-FISH is a viable methodological approach for the analysis of genome organization within complex LCR regions.

Next generation quantitative genetics in plants

Most characteristics in living organisms show continuous variation, which suggests that they are controlled by multiple genes. Quantitative trait loci (QTL) analysis can identify the genes underlying continuous traits by establishing associations between genetic markers and observed phenotypic variation in a segregating population. The new high-throughput sequencing (HTS) technologies greatly facilitate QTL analysis by providing genetic markers at genome-wide resolution in any species without previous knowledge of its genome. In addition HTS serves to quantify molecular phenotypes, which aids to identify the loci responsible for QTLs and to understand the mechanisms underlying diversity. The constant improvements in price, experimental protocols, computational pipelines, and statistical frameworks are making feasible the use of HTS for any research group interested in quantitative genetics. In this review I discuss the application of HTS for molecular marker discovery, population genotyping, and expression profiling in QTL analysis.

Unmasking of a Recessive SCARF2 Mutation by a 22q11.12 de novo Deletion in a Patient with Van den Ende-Gupta Syndrome

Van den Ende-Gupta syndrome (VDEGS) is a congenital condition characterized by craniofacial and skeletal manifestations, specifically blepharophimosis, malar and maxillary hypoplasia, distinctive nose, arachnocamptodactyly, and long slender bones of the hands and feet. To date, only 24 patients have been described. It is generally thought that the syndrome is transmitted by an autosomal recessive mode of inheritance, although evidence for genetic heterogeneity has recently been presented. We report on a girl followed from birth up to 3 years of life with a set of peculiar minor anomalies, arachnocamptodactyly of hands and feet, characteristic of VDEGS in association with a 22q11.12 deletion. Recently, the VDEGS gene was mapped to the DiGeorge syndrome region on 22q11.2, and homozygous mutations in the SCARF2 gene were identified. We now report the first patient with VDEGS due to compound heterozygosity for the common 22q11.2 microdeletion and a hemizygous SCARF2 splice site mutation.