A chromosomal rearrangement hotspot can be identified from population genetic variation and is coincident with a hotspot for allelic recombination

AGAP duplicons associate with structural diversity at Chromosome 10q11.22

The 10q11.22 chromosomal region is a duplication-rich interval of the human genome and one of the last to be fully assembled. It carries copy number-variable genes associated with intellectual disability, bipolar disorder, and obesity. In this study, we characterized the structural diversity at this locus by analyzing 64 haploid assemblies produced by the Human Pangenome Reference Consortium. We identified 11 alternative haplotypes that differ in the copy number and/or orientation of large genomic segments, ranging from hundreds of kilobase pairs (kbp) to over one megabase pair (Mbp). We uncovered a 2.4 Mbp size difference between the shortest and longest haplotypes. Breakpoint analysis revealed that genomic instability results from nonallelic homologous recombination between segmental duplication (SD) pairs with varying similarity (94.4%-99.6%). Nonetheless, these pairs generally recombine at positions where their identity is higher (>99.6%). Recurrent inversions occur with different breakpoints within the same inverted SD pair. Inversion polymorphisms shuffle the entire SD arrangement, creating new predispositions to copy-number variations. The SD architecture is associated with a catarrhine-specific subgroup of the AGAP gene family, which likely triggered the accumulation of SDs at this locus over the past 25 million years of human evolution. Our results reveal extensive structural diversity and genomic instability at the 10q11.22 locus, and expand the general understanding of the mutational mechanisms behind SD-mediated rearrangements.

Inverted triplications formed by iterative template switches generate structural variant diversity at genomic disorder loci

The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a complex genomic rearrangement (CGR). Although it has been identified as an important pathogenic DNA mutation signature in genomic disorders and cancer genomes, its architecture remains unresolved. Here, we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the DNA of 24 patients identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted structural variant (SV) haplotypes. Using a combination of short-read genome sequencing (GS), long-read GS, optical genome mapping, and single-cell DNA template strand sequencing (strand-seq), the haplotype structure was resolved in 18 samples. The point of template switching in 4 samples was shown to be a segment of ∼2.2-5.5 kb of 100% nucleotide similarity within inverted repeat pairs. These data provide experimental evidence that inverted low-copy repeats act as recombinant substrates. This type of CGR can result in multiple conformers generating diverse SV haplotypes in susceptible dosage-sensitive loci.

Break-induced replication underlies formation of inverted triplications and generates unexpected diversity in haplotype structures

Background

The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a type of complex genomic rearrangement (CGR) hypothesized to result from replicative repair of DNA due to replication fork collapse. It is often mediated by a pair of inverted low-copy repeats (LCR) followed by iterative template switches resulting in at least two breakpoint junctions in cis . Although it has been identified as an important mutation signature of pathogenicity for genomic disorders and cancer genomes, its architecture remains unresolved and is predicted to display at least four structural variation (SV) haplotypes.

Results

Here we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the genomic DNA of 24 patients with neurodevelopmental disorders identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted SV haplotypes. Using a combination of short-read genome sequencing (GS), long- read GS, optical genome mapping and StrandSeq the haplotype structure was resolved in 18 samples. This approach refined the point of template switching between inverted LCRs in 4 samples revealing a DNA segment of ∼2.2-5.5 kb of 100% nucleotide similarity. A prediction model was developed to infer the LCR used to mediate the non-allelic homology repair.

Conclusions

These data provide experimental evidence supporting the hypothesis that inverted LCRs act as a recombinant substrate in replication-based repair mechanisms. Such inverted repeats are particularly relevant for formation of copy-number associated inversions, including the DUP-TRP/INV-DUP structures. Moreover, this type of CGR can result in multiple conformers which contributes to generate diverse SV haplotypes in susceptible loci .

CRISPR/Cas9-induced gene conversion between ATAD3 paralogs

Paralogs and pseudogenes are abundant within the human genome, and can mediate non-allelic homologous recombination (NAHR) or gene conversion events. The ATAD3 locus contains three paralogs situated in tandem, and is therefore prone to NAHR-mediated deletions and duplications associated with severe neurological phenotypes. To study this locus further, we aimed to generate biallelic loss-of-function variants in ATAD3A by CRISPR/Cas9 genome editing. Unexpectedly, two of the generated clones underwent gene conversion, as evidenced by replacement of the targeted sequence of ATAD3A by a donor sequence from its paralog ATAD3B. We highlight the complexity of CRISPR/Cas9 design, end-product formation, and recombination repair mechanisms for CRISPR/Cas9 delivery as a nucleic acid molecular therapy when targeting genes that have paralogs or pseudogenes, and advocate meticulous evaluation of resultant clones in model organisms. In addition, we suggest that endogenous gene conversion may be used to repair missense variants in genes with paralogs or pseudogenes.

Classification of NF1 microdeletions and its importance for establishing genotype/phenotype correlations in patients with NF1 microdeletions

An estimated 5–11% of patients with neurofibromatosis type-1 (NF1) harbour large deletions encompassing the NF1 gene and flanking regions. These NF1 microdeletions are subclassified into type 1, 2, 3 and atypical deletions which are distinguishable from each other by their extent and by the number of genes included within the deletion regions as well as the frequency of mosaicism with normal cells. Most common are type-1 NF1 deletions which encompass 1.4-Mb and 14 protein-coding genes. Type-1 deletions are frequently associated with overgrowth, global developmental delay, cognitive disability and dysmorphic facial features which are uncommon in patients with intragenic pathogenic NF1 gene variants. Further, patients with type-1 NF1 deletions frequently exhibit high numbers of neurofibromas and have an increased risk of malignant peripheral nerve sheath tumours. Genes located within the type-1 NF1 microdeletion interval and co-deleted with NF1 are likely to act as modifiers responsible for the severe disease phenotype in patients with NF1 microdeletions, thereby causing the NF1 microdeletion syndrome. Genotype/phenotype correlations in patients with NF1 microdeletions of different lengths are important to identify such modifier genes. However, these correlations are critically dependent upon the accurate characterization of the deletions in terms of their extent. In this review, we outline the utility as well as the shortcomings of multiplex ligation-dependent probe amplification (MLPA) to classify the different types of NF1 microdeletion and indicate the importance of high-resolution microarray analysis for correct classification, a necessary precondition to identify those genes responsible for the NF1 microdeletion syndrome.

Assisted reproduction mediated resurrection of a feline model for Chediak-Higashi syndrome caused by a large duplication in LYST

Chediak-Higashi Syndrome (CHS) is a well-characterized, autosomal recessively inherited lysosomal disease caused by mutations in lysosomal trafficking regulator (LYST). The feline model for CHS was originally maintained for ~20 years. However, the colonies were disbanded and the CHS cat model was lost to the research community before the causative mutation was identified. To resurrect the cat model, semen was collected and cryopreserved from a lone, fertile,  CHS carrier male. Using cryopreserved semen, laparoscopic oviductal artificial insemination was performed on three queens, two queens produced 11 viable kittens. To identify the causative mutation, a fibroblast cell line, derived from an affected cat from the original colony, was whole genome sequenced. Visual inspection of the sequence data identified a candidate causal variant as a ~20 kb tandem duplication within LYST, spanning exons 30 through to 38 (NM_001290242.1:c.8347-2422_9548 + 1749dup). PCR genotyping of the produced offspring demonstrated three individuals inherited the mutant allele from the CHS carrier male. This study demonstrated the successful use of cryopreservation and assisted reproduction to maintain and resurrect biomedical models and has defined the variant causing Chediak-Higashi syndrome in the domestic cat.

Into the bloom: Molecular response of pelagic tunicates to fluctuating food availability

The planktonic tunicates appendicularians and thaliaceans are highly efficient filter feeders on a wide range of prey size including bacteria and have shorter generation times than any other marine grazers. These traits allow some tunicate species to reach high population densities and ensure their success in a favourable environment. However, there are still few studies focusing on which genes and gene pathways are associated with responses of pelagic tunicates to environmental variability. Herein, we present the effect of food availability increase on tunicate community and gene expression at the Marquesas Islands (South-East Pacific Ocean). By using data from the Tara Oceans expedition, we show that changes in phytoplankton density and composition trigger the success of a dominant larvacean species (an undescribed appendicularian). Transcriptional signature to the autotroph bloom suggests key functions in specific physiological processes, i.e., energy metabolism, muscle contraction, membrane trafficking, and proteostasis. The relative abundance of reverse transcription-related Pfams was lower at bloom conditions, suggesting a link with adaptive genetic diversity in tunicates in natural ecosystems. Downstream of the bloom, pelagic tunicates were outcompeted by copepods. Our work represents the first metaomics study of the biological effects of phytoplankton bloom on a key zooplankton taxon.

Predicting human genes susceptible to genomic instability associated with Alu/Alu-mediated rearrangements

Alu elements, the short interspersed element numbering more than 1 million copies per human genome, can mediate the formation of copy number variants (CNVs) between substrate pairs. These Alu/Alu-mediated rearrangements (AAMRs) can result in pathogenic variants that cause diseases. To investigate the impact of AAMR on gene variation and human health, we first characterized Alus that are involved in mediating CNVs (CNV-Alus) and observed that these Alus tend to be evolutionarily younger. We then computationally generated, with the assistance of a supercomputer, a test data set consisting of 78 million Alu pairs and predicted ∼18% of them are potentially susceptible to AAMR. We further determined the relative risk of AAMR in 12,074 OMIM genes using the count of predicted CNV-Alu pairs and experimentally validated the predictions with 89 samples selected by correlating predicted hotspots with a database of CNVs identified by clinical chromosomal microarrays (CMAs) on the genomes of approximately 54,000 subjects. We fine-mapped 47 duplications, 40 deletions, and two complex rearrangements and examined a total of 52 breakpoint junctions of simple CNVs. Overall, 94% of the candidate breakpoints were at least partially Alu mediated. We successfully predicted all (100%) of Alu pairs that mediated deletions (n = 21) and achieved an 87% positive predictive value overall when including AAMR-generated deletions and duplications. We provided a tool, AluAluCNVpredictor, for assessing AAMR hotspots and their role in human disease. These results demonstrate the utility of our predictive model and provide insights into the genomic features and molecular mechanisms underlying AAMR.

The coexistence of copy number variations (CNVs) and single nucleotide polymorphisms (SNPs) at a locus can result in distorted calculations of the significance in associating SNPs to disease

With the recent advance in genome-wide association studies (GWAS), disease-associated single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) have been extensively reported. Accordingly, the issue of incorrect identification of recombination events that can induce the distortion of multi-allelic or hemizygous variants has received more attention. However, the potential distorted calculation bias or significance of a detected association in a GWAS due to the coexistence of CNVs and SNPs in the same genomic region may remain under-recognized. Here we performed the association study within a congenital scoliosis (CS) cohort whose genetic etiology was recently elucidated as a compound inheritance model, including mostly one rare variant deletion CNV null allele and one common variant non-coding hypomorphic haplotype of the TBX6 gene. We demonstrated that the existence of a deletion in TBX6 led to an overestimation of the contribution of the SNPs on the hypomorphic allele. Furthermore, we generalized a model to explain the calculation bias, or distorted significance calculation for an association study, that can be 'induced' by CNVs at a locus. Meanwhile, overlapping between the disease-associated SNPs from published GWAS and common CNVs (overlap 10%) and pathogenic/likely pathogenic CNVs (overlap 99.69%) was significantly higher than the random distribution (p < 1 × 10^-6 and p = 0.034, respectively), indicating that such co-existence of CNV and SNV alleles might generally influence data interpretation and potential outcomes of a GWAS. We also verified and assessed the influence of colocalizing CNVs to the detection sensitivity of disease-associated SNP variant alleles in another adolescent idiopathic scoliosis (AIS) genome-wide association study. We proposed that detecting co-existent CNVs when evaluating the association signals between SNPs and disease traits could improve genetic model analyses and better integrate GWAS with robust Mendelian principles.

Pronounced maternal parent-of-origin bias for type-1 NF1 microdeletions

Neurofibromatosis type 1 (NF1) is caused, in 4.7-11% of cases, by large deletions encompassing the NF1 gene and its flanking regions within 17q11.2. Different types of large NF1 deletion occur which are distinguishable by their breakpoint location and underlying mutational mechanism. Most common are the type-1 NF1 deletions of 1.4 Mb which exhibit recurrent breakpoints caused by nonallelic homologous recombination (NAHR), also termed unequal crossover. Here, we analyzed 37 unrelated families of patients with de novo type-1 NF1 deletions by means of short tandem repeat (STR) profiling to determine the parental origin of the deletions. We observed that 33 of the 37 type-1 deletions were of maternal origin (89.2% of cases; p < 0.0001). Analysis of the patients' siblings indicated that, in 14 informative cases, ten (71.4%) deletions resulted from interchromosomal unequal crossover during meiosis I. Our findings indicate a strong maternal parent-of-origin bias for type-1 NF1 deletions. A similarly pronounced maternal transmission bias has been reported for recurrent copy number variants (CNVs) within 16p11.2 associated with autism, but not so far for any other NAHR-mediated pathogenic CNVs. Region-specific genomic features are likely to be responsible for the maternal bias in the origin of both the 16p11.2 CNVs and type-1 NF1 deletions.

The impact of recombination on human mutation load and disease

Recombination promotes genomic integrity among cells and tissues through double-strand break repair, and is critical for gamete formation and fertility through a strict regulation of the molecular mechanisms associated with proper chromosomal disjunction. In humans, congenital defects and recurrent structural abnormalities can be attributed to aberrant meiotic recombination. Moreover, mutations affecting genes involved in recombination pathways are directly linked to pathologies including infertility and cancer. Recombination is among the most prominent mechanism shaping genome variation, and is associated with not only the structuring of genomic variability, but is also tightly linked with the purging of deleterious mutations from populations. Together, these observations highlight the multiple roles of recombination in human genetics: its ability to act as a major force of evolution, its molecular potential to maintain genome repair and integrity in cell division and its mutagenic cost impacting disease evolution.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Consideration of the haplotype diversity at nonallelic homologous recombination hotspots improves the precision of rearrangement breakpoint identification

Precise characterization of nonallelic homologous recombination (NAHR) breakpoints is key to identifying those features that influence NAHR frequency. Until now, analysis of NAHR-mediated rearrangements has generally been performed by comparison of the breakpoint-spanning sequences with the human genome reference sequence. We show here that the haplotype diversity of NAHR hotspots may interfere with breakpoint-mapping. We studied the transmitting parents of individuals with germline type-1 NF1 deletions mediated by NAHR within the paralogous recombination site 1 (PRS1) or paralogous recombination site 2 (PRS2) hotspots. Several parental wild-type PRS1 and PRS2 haplotypes were identified that exhibited considerable sequence differences with respect to the reference sequence, which also affected the number of predicted PRDM9-binding sites. Sequence comparisons between the parental wild-type PRS1 or PRS2 haplotypes and the deletion breakpoint-spanning sequences from the patients (method #2) turned out to be an accurate means to assign NF1 deletion breakpoints and proved superior to crude reference sequence comparisons that neglect to consider haplotype diversity (method #1). The mean length of the deletion breakpoint regions assigned by method #2 was 269-bp in contrast to 502-bp by method #1. Our findings imply that paralog-specific haplotype diversity of NAHR hotspots (such as PRS2) and population-specific haplotype diversity must be taken into account in order to accurately ascertain NAHR-mediated rearrangement breakpoints.

The next generation of population-based spinal muscular atrophy carrier screening: comprehensive pan-ethnic SMN1 copy-number and sequence variant analysis by massively parallel sequencing

PURPOSE

To investigate pan-ethnic SMN1 copy-number and sequence variation by hybridization-based target enrichment coupled with massively parallel sequencing or next-generation sequencing (NGS).

METHODS

NGS reads aligned to SMN1 and SMN2 exon 7 were quantified to determine the total combined copy number of SMN1 and SMN2. The ratio of SMN1 to SMN2 was calculated based on a single-nucleotide difference that distinguishes the two genes. SMN1 copy-number results were compared between the NGS and quantitative polymerase chain reaction and/or multiplex ligation-dependent probe amplification. The NGS data set was also queried for the g.27134T>G single-nucleotide polymorphism (SNP) and other SMN1 sequence pathogenic variants.

RESULTS

The sensitivity of the test to detect spinal muscular atrophy (SMA) carriers with one copy of SMN1 was 100% (95% confidence interval (CI): 95.9-100%; n = 90) and specificity was 99.6% (95% CI: 99.4-99.7%; n = 6,648). Detection of the g.27134T>G SNP by NGS was 100% concordant with an restriction fragment-length polymorphism method (n = 493). Ten single-nucleotide variants in SMN1 were detectable by NGS and confirmed by gene-specific amplicon-based sequencing. This comprehensive approach yielded SMA carrier detection rates of 90.3-95.0% in five ethnic groups studied.

CONCLUSION

We have developed a novel, comprehensive SMN1 copy-number and sequence variant analysis method by NGS that demonstrated improved SMA carrier detection rates across the entire population examined.Genet Med advance online publication 19 January 2017.

SRBreak: A Read-Depth and Split-Read Framework to Identify Breakpoints of Different Events Inside Simple Copy-Number Variable Regions

Copy-number variation (CNV) has been associated with increased risk of complex diseases. High-throughput sequencing (HTS) technologies facilitate the detection of copy-number variable regions (CNVRs) and their breakpoints. This helps in understanding genome structure as well as their evolution process. Various approaches have been proposed for detecting CNV breakpoints, but currently it is still challenging for tools based on a single analysis method to identify breakpoints of CNVs. It has been shown, however, that pipelines which integrate multiple approaches are able to report more reliable breakpoints. Here, based on HTS data, we have developed a pipeline to identify approximate breakpoints (±10 bp) relating to different ancestral events within a specific CNVR. The pipeline combines read-depth and split-read information to infer breakpoints, using information from multiple samples to allow an imputation approach to be taken. The main steps involve using a normal mixture model to cluster samples into different groups, followed by simple kernel-based approaches to maximize information obtained from read-depth and split-read approaches, after which common breakpoints of groups are inferred. The pipeline uses split-read information directly from CIGAR strings of BAM files, without using a re-alignment step. On simulated data sets, it was able to report breakpoints for very low-coverage samples including those for which only single-end reads were available. When applied to three loci from existing human resequencing data sets (NEGR1, LCE3, IRGM) the pipeline obtained good concordance with results from the 1000 Genomes Project (92, 100, and 82%, respectively). The package is available at https://github.com/hoangtn/SRBreak, and also as a docker-based application at https://registry.hub.docker.com/u/hoangtn/srbreak/.

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.

Mechanisms of germ line genome instability

During meiosis, numerous DNA double-strand breaks (DSBs) are formed as part of the normal developmental program. This seemingly destructive behavior is necessary for successful meiosis, since repair of the DSBs through homologous recombination (HR) helps to produce physical links between the homologous chromosomes essential for correct chromosome segregation later in meiosis. However, DSB formation at such a massive scale also introduces opportunities to generate gross chromosomal rearrangements. In this review, we explore ways in which meiotic DSBs can result in such genomic alterations.

Multiallelic Positions in the Human Genome: Challenges for Genetic Analyses

As the amount of human genomic sequence available from personal genomes and exomes has increased, so too has the observation of genomic positions having two or more alternative alleles, so-called multiallelic sites. For portions of the haploid genome that are present in more than one copy, including segmental duplications, variation at such multisite variant positions becomes even more complex. Despite the frequency of multiallelic variants, a number of commonly used resources and tools in genomic research and diagnostics do not support these multiallelic variants all together or require special modifications. Here, we explore the frequency of multiallelic sites in large samples with whole exome sequencing and discuss potential outcomes of failing to account for multiple variant alleles. We also briefly discuss some commonly utilized resources that fully support multiallelic sites.

Fine mapping of meiotic NAHR-associated crossovers causing large NF1 deletions

Large deletions encompassing the NF1 gene and its flanking regions belong to the group of genomic disorders caused by copy number changes that are mediated by the local genomic architecture. Although nonallelic homologous recombination (NAHR) is known to be a major mutational mechanism underlying such genomic copy number changes, the sequence determinants of NAHR location and frequency are still poorly understood since few high-resolution mapping studies of NAHR hotspots have been performed to date. Here, we have characterized two NAHR hotspots, PRS1 and PRS2, separated by 20 kb and located within the low-copy repeats NF1-REPa and NF1-REPc, which flank the human NF1 gene region. High-resolution mapping of the crossover sites identified in 78 type 1 NF1 deletions mediated by NAHR indicated that PRS2 is a much stronger NAHR hotspot than PRS1 since 80% of these deletions exhibited crossovers within PRS2, whereas 20% had crossovers within PRS1. The identification of the most common strand exchange regions of these 78 deletions served to demarcate the cores of the PRS1 and PRS2 hotspots encompassing 1026 and 1976 bp, respectively. Several sequence features were identified that may influence hotspot intensity and direct the positional preference of NAHR to the hotspot cores. These features include regions of perfect sequence identity encompassing 700 bp at the hotspot core, the presence of PRDM9 binding sites perfectly matching the consensus motif for the most common PRDM9 variant, specific pre-existing patterns of histone modification and open chromatin conformations that are likely to facilitate PRDM9 binding.

Structural variation mutagenesis of the human genome: Impact on disease and evolution

Watson-Crick base-pair changes, or single-nucleotide variants (SNV), have long been known as a source of mutations. However, the extent to which DNA structural variation, including duplication and deletion copy number variants (CNV) and copy number neutral inversions and translocations, contribute to human genome variation and disease has been appreciated only recently. Moreover, the potential complexity of structural variants (SV) was not envisioned; thus, the frequency of complex genomic rearrangements and how such events form remained a mystery. The concept of genomic disorders, diseases due to genomic rearrangements and not sequence-based changes for which genomic architecture incite genomic instability, delineated a new category of conditions distinct from chromosomal syndromes and single-gene Mendelian diseases. Nevertheless, it is the mechanistic understanding of CNV/SV formation that has promoted further understanding of human biology and disease and provided insights into human genome and gene evolution. Environ. Mol. Mutagen. 56:419-436, 2015. © 2015 Wiley Periodicals, Inc.

Complex genomic rearrangements at the PLP1 locus include triplication and quadruplication

Inverted repeats (IRs) can facilitate structural variation as crucibles of genomic rearrangement. Complex duplication-inverted triplication-duplication (DUP-TRP/INV-DUP) rearrangements that contain breakpoint junctions within IRs have been recently associated with both MECP2 duplication syndrome (MIM#300260) and Pelizaeus-Merzbacher disease (PMD, MIM#312080). We investigated 17 unrelated PMD subjects with copy number gains at the PLP1 locus including triplication and quadruplication of specific genomic intervals-16/17 were found to have a DUP-TRP/INV-DUP rearrangement product. An IR distal to PLP1 facilitates DUP-TRP/INV-DUP formation as well as an inversion structural variation found frequently amongst normal individuals. We show that a homology-or homeology-driven replicative mechanism of DNA repair can apparently mediate template switches within stretches of microhomology. Moreover, we provide evidence that quadruplication and potentially higher order amplification of a genomic interval can occur in a manner consistent with rolling circle amplification as predicted by the microhomology-mediated break induced replication (MMBIR) model.

Bayesian inference of shared recombination hotspots between humans and chimpanzees

Recombination generates variation and facilitates evolution. Recombination (or lack thereof) also contributes to human genetic disease. Methods for mapping genes influencing complex genetic diseases via association rely on linkage disequilibrium (LD) in human populations, which is influenced by rates of recombination across the genome. Comparative population genomic analyses of recombination using related primate species can identify factors influencing rates of recombination in humans. Such studies can indicate how variable hotspots for recombination may be both among individuals (or populations) and over evolutionary timescales. Previous studies have suggested that locations of recombination hotspots are not conserved between humans and chimpanzees. We made use of the data sets from recent resequencing projects and applied a Bayesian method for identifying hotspots and estimating recombination rates. We also reanalyzed SNP data sets for regions with known hotspots in humans using samples from the human and chimpanzee. The Bayes factors (BF) of shared recombination hotspots between human and chimpanzee across regions were obtained. Based on the analysis of the aligned regions of human chromosome 21, locations where the two species show evidence of shared recombination hotspots (with high BFs) were identified. Interestingly, previous comparative studies of human and chimpanzee that focused on the known human recombination hotspots within the β-globin and HLA regions did not find overlapping of hotspots. Our results show high BFs of shared hotspots at locations within both regions, and the estimated locations of shared hotspots overlap with the locations of human recombination hotspots obtained from sperm-typing studies.

The genetics of microdeletion and microduplication syndromes: an update

Chromosomal abnormalities, including microdeletions and microduplications, have long been associated with abnormal developmental outcomes. Early discoveries relied on a common clinical presentation and the ability to detect chromosomal abnormalities by standard karyotype analysis or specific assays such as fluorescence in situ hybridization. Over the past decade, the development of novel genomic technologies has allowed more comprehensive, unbiased discovery of microdeletions and microduplications throughout the human genome. The ability to quickly interrogate large cohorts using chromosome microarrays and, more recently, next-generation sequencing has led to the rapid discovery of novel microdeletions and microduplications associated with disease, including very rare but clinically significant rearrangements. In addition, the observation that some microdeletions are associated with risk for several neurodevelopmental disorders contributes to our understanding of shared genetic susceptibility for such disorders. Here, we review current knowledge of microdeletion/duplication syndromes, with a particular focus on recurrent rearrangement syndromes.

Population-specific differences in gene conversion patterns between human SUZ12 and SUZ12P are indicative of the dynamic nature of interparalog gene conversion

Nonallelic homologous gene conversion (NAHGC) resulting from interparalog recombination without crossover represents an important influence on the evolution of duplicated sequences in the human genome. In 17q11.2, different paralogous sequences mediate large NF1 deletions by nonallelic homologous recombination with crossover (NAHR). Among these paralogs are SUZ12 and its pseudogene SUZ12P which harbour the breakpoints of type-2 (1.2-Mb) NF1 deletions. Such deletions are caused predominantly by mitotic NAHR since somatic mosaicism with normal cells is evident in most patients. Investigating whether SUZ12 and SUZ12P have also been involved in NAHGC, we observed gene conversion tracts between these paralogs in both Africans (AFR) and Europeans (EUR). Since germline type-2 NF1 deletions resulting from meiotic NAHR are very rare, the vast majority of the gene conversion tracts in SUZ12 and SUZ12P are likely to have resulted from mitotic recombination during premeiotic cell divisions of germ cells. A higher number of gene conversion tracts were noted within SUZ12 and SUZ12P in AFR as compared to EUR. Further, the distinctive signature of NAHGC (a high number of SNPs per paralog and a high number of shared SNPs between paralogs), a characteristic of many actively recombining paralogs, was observed in both SUZ12 and SUZ12P but only in AFR and not in EUR. A novel polymorphic 2.3-kb deletion in SUZ12P was identified which exhibited a high allele frequency in EUR. We postulate that this interparalog structural difference, together with low allelic recombination rates, could have caused a reduction in NAHGC between SUZ12 and SUZ12P during human evolution.

Analysis of crossover breakpoints yields new insights into the nature of the gene conversion events associated with large NF1 deletions mediated by nonallelic homologous recombination

Large NF1 deletions are mediated by nonallelic homologous recombination (NAHR). An in-depth analysis of gene conversion operating in the breakpoint-flanking regions of large NF1 deletions was performed to investigate whether the rate of discontinuous gene conversion during NAHR with crossover is increased, as has been previously noted in NAHR-mediated rearrangements. All 20 germline type-1 NF1 deletions analyzed were mediated by NAHR associated with continuous gene conversion within the breakpoint-flanking regions. Continuous gene conversion was also observed in 31/32 type-2 NF1 deletions investigated. In contrast to the meiotic type-1 NF1 deletions, type-2 NF1 deletions are predominantly of post-zygotic origin. Our findings therefore imply that the mitotic as well as the meiotic NAHR intermediates of large NF1 deletions are processed by long-patch mismatch repair (MMR), thereby ensuring gene conversion tract continuity instead of the discontinuous gene conversion that is characteristic of short-patch repair. However, the single type-2 NF1 deletion not exhibiting continuous gene conversion was processed without MMR, yielding two different deletion-bearing chromosomes, which were distinguishable in terms of their breakpoint positions. Our findings indicate that MMR failure during NAHR, followed by post-meiotic/mitotic segregation, has the potential to give rise to somatic mosaicism in human genomic rearrangements by generating breakpoint heterogeneity.

Interaction-based evolution: how natural selection and nonrandom mutation work together

BACKGROUND

The modern evolutionary synthesis leaves unresolved some of the most fundamental, long-standing questions in evolutionary biology: What is the role of sex in evolution? How does complex adaptation evolve? How can selection operate effectively on genetic interactions? More recently, the molecular biology and genomics revolutions have raised a host of critical new questions, through empirical findings that the modern synthesis fails to explain: for example, the discovery of de novo genes; the immense constructive role of transposable elements in evolution; genetic variance and biochemical activity that go far beyond what traditional natural selection can maintain; perplexing cases of molecular parallelism; and more.

PRESENTATION OF THE HYPOTHESIS

Here I address these questions from a unified perspective, by means of a new mechanistic view of evolution that offers a novel connection between selection on the phenotype and genetic evolutionary change (while relying, like the traditional theory, on natural selection as the only source of feedback on the fit between an organism and its environment). I hypothesize that the mutation that is of relevance for the evolution of complex adaptation-while not Lamarckian, or "directed" to increase fitness-is not random, but is instead the outcome of a complex and continually evolving biological process that combines information from multiple loci into one. This allows selection on a fleeting combination of interacting alleles at different loci to have a hereditary effect according to the combination's fitness.

TESTING AND IMPLICATIONS OF THE HYPOTHESIS

This proposed mechanism addresses the problem of how beneficial genetic interactions can evolve under selection, and also offers an intuitive explanation for the role of sex in evolution, which focuses on sex as the generator of genetic combinations. Importantly, it also implies that genetic variation that has appeared neutral through the lens of traditional theory can actually experience selection on interactions and thus has a much greater adaptive potential than previously considered. Empirical evidence for the proposed mechanism from both molecular evolution and evolution at the organismal level is discussed, and multiple predictions are offered by which it may be tested.

REVIEWERS

This article was reviewed by Nigel Goldenfeld (nominated by Eugene V. Koonin), Jürgen Brosius and W. Ford Doolittle.

NAHR-mediated copy-number variants in a clinical population: mechanistic insights into both genomic disorders and Mendelizing traits

We delineated and analyzed directly oriented paralogous low-copy repeats (DP-LCRs) in the most recent version of the human haploid reference genome. The computationally defined DP-LCRs were cross-referenced with our chromosomal microarray analysis (CMA) database of 25,144 patients subjected to genome-wide assays. This computationally guided approach to the empirically derived large data set allowed us to investigate genomic rearrangement relative frequencies and identify new loci for recurrent nonallelic homologous recombination (NAHR)-mediated copy-number variants (CNVs). The most commonly observed recurrent CNVs were NPHP1 duplications (233), CHRNA7 duplications (175), and 22q11.21 deletions (DiGeorge/velocardiofacial syndrome, 166). In the ∼25% of CMA cases for which parental studies were available, we identified 190 de novo recurrent CNVs. In this group, the most frequently observed events were deletions of 22q11.21 (48), 16p11.2 (autism, 34), and 7q11.23 (Williams-Beuren syndrome, 11). Several features of DP-LCRs, including length, distance between NAHR substrate elements, DNA sequence identity (fraction matching), GC content, and concentration of the homologous recombination (HR) hot spot motif 5'-CCNCCNTNNCCNC-3', correlate with the frequencies of the recurrent CNVs events. Four novel adjacent DP-LCR-flanked and NAHR-prone regions, involving 2q12.2q13, were elucidated in association with novel genomic disorders. Our study quantitates genome architectural features responsible for NAHR-mediated genomic instability and further elucidates the role of NAHR in human disease.

The role of gene conversion in preserving rearrangement hotspots in the human genome

Hotspots of non-allelic homologous recombination (NAHR) have a crucial role in creating genetic diversity and are also associated with dozens of genomic disorders. Recent studies suggest that many human NAHR hotspots have been preserved throughout the evolution of primates. NAHR hotspots are likely to remain active as long as the segmental duplications (SDs) promoting NAHR retain sufficient similarity. Here, we propose an evolutionary model of SDs that incorporates the effect of gene conversion and compare it with a null model that assumes SDs evolve independently without gene conversion. The gene conversion model predicts a much longer lifespan of NAHR hotspots compared with the null model. We show that the literature on copy number variants (CNVs) and genomic disorders, and also the results of additional analysis of CNVs, are all more consistent with the gene conversion model.

Duplication and retention biases of essential and non-essential genes revealed by systematic knockdown analyses

When a duplicate gene has no apparent loss-of-function phenotype, it is commonly considered that the phenotype has been masked as a result of functional redundancy with the remaining paralog. This is supported by indirect evidence showing that multi-copy genes show loss-of-function phenotypes less often than single-copy genes and by direct tests of phenotype masking using select gene sets. Here we take a systematic genome-wide RNA interference approach to assess phenotype masking in paralog pairs in the Caenorhabditis elegans genome. Remarkably, in contrast to expectations, we find that phenotype masking makes only a minor contribution to the low knockdown phenotype rate for duplicate genes. Instead, we find that non-essential genes are highly over-represented among duplicates, leading to a low observed loss-of-function phenotype rate. We further find that duplicate pairs derived from essential and non-essential genes have contrasting evolutionary dynamics: whereas non-essential genes are both more often successfully duplicated (fixed) and lost, essential genes are less often duplicated but upon successful duplication are maintained over longer periods. We expect the fundamental evolutionary duplication dynamics presented here to be broadly applicable.

Duplicate genes occur in all organisms. It has been found that mutations in duplicate genes cause defects much less often than when single copy genes are mutated. It is widely believed that this is due to functional redundancy—that is, the two genes can carry out similar functions so that the non-mutated duplicate gene can cover for or “mask” the phenotype of the mutation in the first duplicate. To determine whether this hypothesis is true, it is necessary to test systematically whether defects indeed occur in the organism when both duplicate genes are inhibited. We have for the first time carried out such an analysis in a multicellular organism, the nematode Caenorhabditis elegans. In contrast to expectations, we observed that when both copies of duplicate genes are inhibited deleterious effects are very rare. We show that this is because duplicate genes are much more often non-essential compared to genes where there is only a single copy. Non-essential genes are also lost from the genome much more often than essential genes. However, when essential genes are duplicated, they remain present in the genome over longer periods. Our results give a framework to explain the evolutionary dynamics of duplications in the genome.

Differential relationship of DNA replication timing to different forms of human mutation and variation

Human genetic variation is distributed nonrandomly across the genome, though the principles governing its distribution are only partially known. DNA replication creates opportunities for mutation, and the timing of DNA replication correlates with the density of SNPs across the human genome. To enable deeper investigation of how DNA replication timing relates to human mutation and variation, we generated a high-resolution map of the human genome's replication timing program and analyzed its relationship to point mutations, copy number variations, and the meiotic recombination hotspots utilized by males and females. DNA replication timing associated with point mutations far more strongly than predicted from earlier analyses and showed a stronger relationship to transversion than transition mutations. Structural mutations arising from recombination-based mechanisms and recombination hotspots used more extensively by females were enriched in early-replicating parts of the genome, though these relationships appeared to relate more strongly to the genomic distribution of causative sequence features. These results indicate differential and sex-specific relationship of DNA replication timing to different forms of mutation and recombination.

A common copy-number breakpoint of ERBB2 amplification in breast cancer colocalizes with a complex block of segmental duplications

Introduction

Segmental duplications (low-copy repeats) are the recently duplicated genomic segments in the human genome that display nearly identical (> 90%) sequences and account for about 5% of euchromatic regions. In germline, duplicated segments mediate nonallelic homologous recombination and thus cause both non-disease-causing copy-number variants and genomic disorders. To what extent duplicated segments play a role in somatic DNA rearrangements in cancer remains elusive. Duplicated segments often cluster and form genomic blocks enriched with both direct and inverted repeats (complex genomic regions). Such complex regions could be fragile and play a mechanistic role in the amplification of the ERBB2 gene in breast tumors, because repeated sequences are known to initiate gene amplification in model systems.

Methods

We conducted polymerase chain reaction (PCR)-based assays for primary breast tumors and analyzed publically available array-comparative genomic hybridization data to map a common copy-number breakpoint in ERBB2-amplified primary breast tumors. We further used molecular, bioinformatics, and population-genetics approaches to define duplication contents, structural variants, and haplotypes within the common breakpoint.

Results

We found a large (> 300-kb) block of duplicated segments that was colocalized with a common-copy number breakpoint for ERBB2 amplification. The breakpoint that potentially initiated ERBB2 amplification localized in a region 1.5 megabases (Mb) on the telomeric side of ERBB2. The region is very complex, with extensive duplications of KRTAP genes, structural variants, and, as a result, a paucity of single-nucleotide polymorphism (SNP) markers. Duplicated segments are varied in size and degree of sequence homology, indicating that duplications have occurred recurrently during genome evolution.

Conclusions

Amplification of the ERBB2 gene in breast tumors is potentially initiated by a complex region that has unusual genomic features and thus requires rigorous, labor-intensive investigation. The haplotypes we provide could be useful to identify the potential association between the complex region and ERBB2 amplification.

Identification of recurrent type-2 NF1 microdeletions reveals a mitotic nonallelic homologous recombination hotspot underlying a human genomic disorder

Nonallelic homologous recombination (NAHR) is one of the major mechanisms underlying copy number variation in the human genome. Although several disease-associated meiotic NAHR breakpoints have been analyzed in great detail, hotspots for mitotic NAHR are not well characterized. Type-2 NF1 microdeletions, which are predominantly of postzygotic origin, constitute a highly informative model with which to investigate the features of mitotic NAHR. Here, a custom-designed MLPA- and PCR-based approach was used to identify 23 novel NAHR-mediated type-2 NF1 deletions. Breakpoint analysis of these 23 type-2 deletions, together with 17 NAHR-mediated type-2 deletions identified previously, revealed that the breakpoints are nonuniformly distributed within the paralogous SUZ12 and SUZ12P sequences. Further, the analysis of this large group of type-2 deletions revealed breakpoint recurrence within short segments (ranging in size from 57 to 253-bp) as well as the existence of a novel NAHR hotspot of 1.9-kb (termed PRS4). This hotspot harbored 20% (8/40) of the type-2 deletion breakpoints and contains the 253-bp recurrent breakpoint region BR6 in which four independent type-2 deletion breakpoints were identified. Our findings indicate that a combination of an open chromatin conformation and short non-B DNA-forming repeats may predispose to recurrent mitotic NAHR events between SUZ12 and its pseudogene.

Evolutionary history of copy-number-variable locus for the low-affinity Fcγ receptor: mutation rate, autoimmune disease, and the legacy of helminth infection

Both sequence variation and copy-number variation (CNV) of the genes encoding receptors for immunoglobulin G (Fcγ receptors) have been genetically and functionally associated with a number of autoimmune diseases. However, the molecular nature and evolutionary context of this variation is unknown. Here, we describe the structure of the CNV, estimate its mutation rate and diversity, and place it in the context of the known functional alloantigen variation of these genes. Deletion of Fcγ receptor IIIB, associated with systemic lupus erythematosus, is a result of independent nonallelic homologous recombination events with a frequency of approximately 0.1%. We also show that pathogen diversity, in particular helminth diversity, has played a critical role in shaping the functional variation at these genes both between mammalian species and between human populations. Positively selected amino acids are involved in the interaction with IgG and include some amino acids that are known polymorphic alloantigens in humans. This supports a genetic contribution to the hygiene hypothesis, which states that past evolution in the context of helminth diversity has left humans with an array of susceptibility alleles for autoimmune disease in the context of a helminth-free environment. This approach shows the link between pathogens and autoimmune disease at the genetic level and provides a strategy for interrogating the genetic variation underlying autoimmune-disease risk and infectious-disease susceptibility.

Small rare recurrent deletions and reciprocal duplications in 2q21.1, including brain-specific ARHGEF4 and GPR148

We have identified a rare small (~450 kb unique sequence) recurrent deletion in a previously linked attention-deficit hyperactivity disorder (ADHD) locus at 2q21.1 in five unrelated families with developmental delay (DD)/intellectual disability (ID), ADHD, epilepsy and other neurobehavioral abnormalities from 17 035 samples referred for clinical chromosomal microarray analysis. Additionally, a DECIPHER (http://decipher.sanger.ac.uk) patient 2311 was found to have the same deletion and presented with aggressive behavior. The deletion was not found in either six control groups consisting of 13 999 healthy individuals or in the DGV database. We have also identified reciprocal duplications in five unrelated families with autism, developmental delay (DD), seizures and ADHD. This genomic region is flanked by large, complex low-copy repeats (LCRs) with directly oriented subunits of ~109 kb in size that have 97.7% DNA sequence identity. We sequenced the deletion breakpoints within the directly oriented paralogous subunits of the flanking LCR clusters, demonstrating non-allelic homologous recombination as a mechanism of formation. The rearranged segment harbors five genes: GPR148, FAM123C, ARHGEF4, FAM168B and PLEKHB2. Expression of ARHGEF4 (Rho guanine nucleotide exchange factor 4) is restricted to the brain and may regulate the actin cytoskeletal network, cell morphology and migration, and neuronal function. GPR148 encodes a G-protein-coupled receptor protein expressed in the brain and testes. We suggest that small rare recurrent deletion of 2q21.1 is pathogenic for DD/ID, ADHD, epilepsy and other neurobehavioral abnormalities and, because of its small size, low frequency and more severe phenotype might have been missed in other previous genome-wide screening studies using single-nucleotide polymorphism analyses.

Mechanisms for recurrent and complex human genomic rearrangements

During the last two decades, the importance of human genome copy number variation (CNV) in disease has become widely recognized. However, much is not understood about underlying mechanisms. We show how, although model organism research guides molecular understanding, important insights are gained from study of the wealth of information available in the clinic. We describe progress in explaining nonallelic homologous recombination (NAHR), a major cause of copy number change occurring when control of allelic recombination fails, highlight the growing importance of replicative mechanisms to explain complex events, and describe progress in understanding extreme chromosome reorganization (chromothripsis). Both nonhomologous end-joining and aberrant replication have significant roles in chromothripsis. As we study CNV, the processes underlying human genome evolution are revealed.

Segmental duplication, microinversion, and gene loss associated with a complex inversion breakpoint region in Drosophila

Chromosomal inversions are usually portrayed as simple two-breakpoint rearrangements changing gene order but not gene number or structure. However, increasing evidence suggests that inversion breakpoints may often have a complex structure and entail gene duplications with potential functional consequences. Here, we used a combination of different techniques to investigate the breakpoint structure and the functional consequences of a complex rearrangement fixed in Drosophila buzzatii and comprising two tandemly arranged inversions sharing the middle breakpoint: 2m and 2n. By comparing the sequence in the breakpoint regions between D. buzzatii (inverted chromosome) and D. mojavensis (noninverted chromosome), we corroborate the breakpoint reuse at the molecular level and infer that inversion 2m was associated with a duplication of a ~13 kb segment and likely generated by staggered breaks plus repair by nonhomologous end joining. The duplicated segment contained the gene CG4673, involved in nuclear transport, and its two nested genes CG5071 and CG5079. Interestingly, we found that other than the inversion and the associated duplication, both breakpoints suffered additional rearrangements, that is, the proximal breakpoint experienced a microinversion event associated at both ends with a 121-bp long duplication that contains a promoter. As a consequence of all these different rearrangements, CG5079 has been lost from the genome, CG5071 is now a single copy nonnested gene, and CG4673 has a transcript ~9 kb shorter and seems to have acquired a more complex gene regulation. Our results illustrate the complex effects of chromosomal rearrangements and highlight the need of complementing genomic approaches with detailed sequence-level and functional analyses of breakpoint regions if we are to fully understand genome structure, function, and evolutionary dynamics.

Characterization of the nonallelic homologous recombination hotspot PRS3 associated with type-3 NF1 deletions

Nonallelic homologous recombination (NAHR) is the major mechanism underlying recurrent genomic rearrangements, including the large deletions at 17q11.2 that cause neurofibromatosis type 1 (NF1). Here, we identify a novel NAHR hotspot, responsible for type-3 NF1 deletions that span 1.0 Mb. Breakpoint clustering within this 1-kb hotspot, termed PRS3, was noted in 10 of 11 known type-3 NF1 deletions. PRS3 is located within the LRRC37B pseudogene of the NF1-REPb and NF1-REPc low-copy repeats. In contrast to other previously characterized NAHR hotspots, PRS3 has not developed on a preexisting allelic homologous recombination hotspot. Furthermore, the variation pattern of PRS3 and its flanking regions is unusual since only NF1-REPc (and not NF1-REPb) is characterized by a high single nucleotide polymorphism (SNP) frequency, suggestive of unidirectional sequence transfer via nonallelic homologous gene conversion (NAHGC). By contrast, the previously described intense NAHR hotspots within the CMT1A-REPs, and the PRS1 and PRS2 hotspots underlying type-1 NF1 deletions, experience frequent bidirectional sequence transfer. PRS3 within NF1-REPc was also found to be involved in NAHGC with the LRRC37B gene, the progenitor locus of the LRRC37B-P duplicons, as indicated by the presence of shared SNPs between these loci. PRS3 therefore represents a weak (and probably evolutionarily rather young) NAHR hotspot with unique properties.

On the sequence-directed nature of human gene mutation: the role of genomic architecture and the local DNA sequence environment in mediating gene mutations underlying human inherited disease

Different types of human gene mutation may vary in size, from structural variants (SVs) to single base-pair substitutions, but what they all have in common is that their nature, size and location are often determined either by specific characteristics of the local DNA sequence environment or by higher order features of the genomic architecture. The human genome is now recognized to contain "pervasive architectural flaws" in that certain DNA sequences are inherently mutation prone by virtue of their base composition, sequence repetitivity and/or epigenetic modification. Here, we explore how the nature, location and frequency of different types of mutation causing inherited disease are shaped in large part, and often in remarkably predictable ways, by the local DNA sequence environment. The mutability of a given gene or genomic region may also be influenced indirectly by a variety of noncanonical (non-B) secondary structures whose formation is facilitated by the underlying DNA sequence. Since these non-B DNA structures can interfere with subsequent DNA replication and repair and may serve to increase mutation frequencies in generalized fashion (i.e., both in the context of subtle mutations and SVs), they have the potential to serve as a unifying concept in studies of mutational mechanisms underlying human inherited disease.

Characterization of copy number-stable regions in the human genome

In the past few years the number of copy number variants (CNVs) identified in the human genome has increased significantly, but our understanding of the functional impact of CNVs is still limited. Clinically significant variations cannot easily be distinguished from benign, complicating interpretation of patient data. Multiple studies have focused on analysis of regions that vary in copy number in specific disorders. Here we use the opposite strategy and focus our analysis on regions that never seem to vary in the general population, hypothesizing that these are copy number stable because variations within them are deleterious. Our results show that copy number stable regions are characterized by correlation with a number of genomic features, allowing us to define a list of genomic regions that are dosage sensitive in humans. We find that these dosage-sensitive regions show significant overlap with de novo CNVs identified in patients with intellectual disability or autism. There is also a significant association between copy number stable regions and rare inherited variants in autism patients, but not in controls. Based on this predictive power, we propose that copy number stable regions can be used to complement maps of known CNVs to facilitate interpretation of patient data.

Copy number gain at Xp22.31 includes complex duplication rearrangements and recurrent triplications

Genomic instability is a feature of the human Xp22.31 region wherein deletions are associated with X-linked ichthyosis, mental retardation and attention deficit hyperactivity disorder. A putative homologous recombination hotspot motif is enriched in low copy repeats that mediate recurrent deletion at this locus. To date, few efforts have focused on copy number gain at Xp22.31. However, clinical testing revealed a high incidence of duplication of Xp22.31 in subjects ascertained and referred with neurobehavioral phenotypes. We systematically studied 61 unrelated subjects with rearrangements revealing gain in copy number, using multiple molecular assays. We detected not only the anticipated recurrent and simple nonrecurrent duplications, but also unexpectedly identified recurrent triplications and other complex rearrangements. Breakpoint analyses enabled us to surmise the mechanisms for many of these rearrangements. The clinical significance of the recurrent duplications and triplications were assessed using different approaches. We cannot find any evidence to support pathogenicity of the Xp22.31 duplication. However, our data suggest that the Xp22.31 duplication may serve as a risk factor for abnormal phenotypes. Our findings highlight the need for more robust Xp22.31 triplication detection in that such further gain may be more penetrant than the duplications. Our findings reveal the distribution of different mechanisms for genomic duplication rearrangements at a given locus, and provide insights into aspects of strand exchange events between paralogous sequences in the human genome.

Intrachromosomal mitotic nonallelic homologous recombination is the major molecular mechanism underlying type-2 NF1 deletions

Nonallelic homologous recombination (NAHR) is responsible for the recurrent rearrangements that give rise to genomic disorders. Although meiotic NAHR has been investigated in multiple contexts, much less is known about mitotic NAHR despite its importance for tumorigenesis. Because type-2 NF1 microdeletions frequently result from mitotic NAHR, they represent a good model in which to investigate the features of mitotic NAHR. We have used microsatellite analysis and SNP arrays to distinguish between the various alternative recombinational possibilities, thereby ascertaining that 17 of 18 type-2 NF1 deletions, with breakpoints in the SUZ12 gene and its highly homologous pseudogene, originated via intrachromosomal recombination. This high proportion of intrachromosomal NAHR causing somatic type-2 NF1 deletions contrasts with the interchromosomal origin of germline type-1 NF1 microdeletions, whose breakpoints are located within the NF1-REPs (low-copy repeats located adjacent to the SUZ12 sequences). Further, meiotic NAHR causing type-1 NF1 deletions occurs within recombination hotspots characterized by high GC-content and DNA duplex stability, whereas the type-2 breakpoints associated with the mitotic NAHR events investigated here do not cluster within hotspots and are located within regions of significantly lower GC-content and DNA stability. Our findings therefore point to fundamental mechanistic differences between the determinants of mitotic and meiotic NAHR.

Structural variation of the human genome: mechanisms, assays, and role in male infertility

Genomic disorders are defined as diseases caused by rearrangements of the genome incited by a genomic architecture that conveys instability. Y-chromosome related dysfunctions such as male infertility are frequently associated with gross DNA rearrangements resulting from its peculiar genomic architecture. The Y-chromosome has evolved into a highly specialized chromosome to perform male functions, mainly spermatogenesis. Direct and inverted repeats, some of them palindromes with highly identical nucleotide sequences that can form DNA cruciform structures, characterize the genomic structure of the Y-chromosome long arm. Some particular Y chromosome genomic deletions can cause spermatogenic failure likely because of removal of one or more transcriptional units with a potential role in spermatogenesis. We describe mechanisms underlying the formation of human genomic rearrangements on autosomes and review Y-chromosome deletions associated with male infertility.

Observation and prediction of recurrent human translocations mediated by NAHR between nonhomologous chromosomes

Four unrelated families with the same unbalanced translocation der(4)t(4;11)(p16.2;p15.4) were analyzed. Both of the breakpoint regions in 4p16.2 and 11p15.4 were narrowed to large ∼359-kb and ∼215-kb low-copy repeat (LCR) clusters, respectively, by aCGH and SNP array analyses. DNA sequencing enabled mapping the breakpoints of one translocation to 24 bp within interchromosomal paralogous LCRs of ∼130 kb in length and 94.7% DNA sequence identity located in olfactory receptor gene clusters, indicating nonallelic homologous recombination (NAHR) as the mechanism for translocation formation. To investigate the potential involvement of interchromosomal LCRs in recurrent chromosomal translocation formation, we performed computational genome-wide analyses and identified 1143 interchromosomal LCR substrate pairs, >5 kb in size and sharing >94% sequence identity that can potentially mediate chromosomal translocations. Additional evidence for interchromosomal NAHR mediated translocation formation was provided by sequencing the breakpoints of another recurrent translocation, der(8)t(8;12)(p23.1;p13.31). The NAHR sites were mapped within 55 bp in ∼7.8-kb paralogous subunits of 95.3% sequence identity located in the ∼579-kb (chr 8) and ∼287-kb (chr 12) LCR clusters. We demonstrate that NAHR mediates recurrent constitutional translocations t(4;11) and t(8;12) and potentially many other interchromosomal translocations throughout the human genome. Furthermore, we provide a computationally determined genome-wide "recurrent translocation map."