Gene conversion causing human inherited disease: evidence for involvement of non-B-DNA-forming sequences and recombination-promoting motifs in DNA breakage and repair

Hereditary Transthyretin-Related Amyloidosis: Genetic Heterogeneity and Early Personalized Gene Therapy

Point mutations of the transthyretin (TTR) gene are related with hereditary amyloidosis (hATTR). The number of people affected by this rare disease is only partially estimated. The real impact of somatic mosaicism and other genetic factors on expressivity, complexity, progression, and transmission of the disease should be better investigated. The relevance of this rare disease is increasing and many efforts have been made to improve the time to diagnosis and to estimate the real number of cases in endemic and non-endemic areas. In this context, somatic mosaicism should be better investigated to explain the complexity of the heterogeneity of the hATTR clinical features, to better estimate the number of new cases, and to focus on early and personalized gene therapy. Gene therapy can potentially improve the living conditions of affected individuals and is one of the central goals in research on amyloidosis related to the TTR gene, with the advantage of overcoming liver transplantation as the sole treatment for hATTR disease.

Distinct sequence features underlie microdeletions and gross deletions in the human genome

Microdeletions and gross deletions are important causes (~20%) of human inherited disease and their genomic locations are strongly influenced by the local DNA sequence environment. This notwithstanding, no study has systematically examined their underlying generative mechanisms. Here, we obtained 42,098 pathogenic microdeletions and gross deletions from the Human Gene Mutation Database (HGMD) that together form a continuum of germline deletions ranging in size from 1 to 28,394,429 bp. We analyzed the DNA sequence within 1 kb of the breakpoint junctions and found that the frequencies of non‐B DNA‐forming repeats, GC‐content, and the presence of seven of 78 specific sequence motifs in the vicinity of pathogenic deletions correlated with deletion length for deletions of length ≤30 bp. Further, we found that the presence of DR, GQ, and STR repeats is important for the formation of longer deletions (>30 bp) but not for the formation of shorter deletions (≤30 bp) while significantly (χ², p < 2E−16) more microhomologies were identified flanking short deletions than long deletions (length >30 bp). We provide evidence to support a functional distinction between microdeletions and gross deletions. Finally, we propose that a deletion length cut‐off of 25–30 bp may serve as an objective means to functionally distinguish microdeletions from gross deletions.

Non-B DNA: a major contributor to small- and large-scale variation in nucleotide substitution frequencies across the genome

Approximately 13% of the human genome can fold into non-canonical (non-B) DNA structures (e.g. G-quadruplexes, Z-DNA, etc.), which have been implicated in vital cellular processes. Non-B DNA also hinders replication, increasing errors and facilitating mutagenesis, yet its contribution to genome-wide variation in mutation rates remains unexplored. Here, we conducted a comprehensive analysis of nucleotide substitution frequencies at non-B DNA loci within noncoding, non-repetitive genome regions, their ±2 kb flanking regions, and 1-Megabase windows, using human-orangutan divergence and human single-nucleotide polymorphisms. Functional data analysis at single-base resolution demonstrated that substitution frequencies are usually elevated at non-B DNA, with patterns specific to each non-B DNA type. Mirror, direct and inverted repeats have higher substitution frequencies in spacers than in repeat arms, whereas G-quadruplexes, particularly stable ones, have higher substitution frequencies in loops than in stems. Several non-B DNA types also affect substitution frequencies in their flanking regions. Finally, non-B DNA explains more variation than any other predictor in multiple regression models for diversity or divergence at 1-Megabase scale. Thus, non-B DNA substantially contributes to variation in substitution frequencies at small and large scales. Our results highlight the role of non-B DNA in germline mutagenesis with implications to evolution and genetic diseases.

Centromeres under Pressure: Evolutionary Innovation in Conflict with Conserved Function

Centromeres are essential genetic elements that enable spindle microtubule attachment for chromosome segregation during mitosis and meiosis. While this function is preserved across species, centromeres display an array of dynamic features, including: (1) rapidly evolving DNA; (2) wide evolutionary diversity in size, shape and organization; (3) evidence of mutational processes to generate homogenized repetitive arrays that characterize centromeres in several species; (4) tolerance to changes in position, as in the case of neocentromeres; and (5) intrinsic fragility derived by sequence composition and secondary DNA structures. Centromere drive underlies rapid centromere DNA evolution due to the "selfish" pursuit to bias meiotic transmission and promote the propagation of stronger centromeres. Yet, the origins of other dynamic features of centromeres remain unclear. Here, we review our current understanding of centromere evolution and plasticity. We also detail the mutagenic processes proposed to shape the divergent genetic nature of centromeres. Changes to centromeres are not simply evolutionary relics, but ongoing shifts that on one side promote centromere flexibility, but on the other can undermine centromere integrity and function with potential pathological implications such as genome instability.

Recombination rates along the entire Epstein Barr virus genome display a highly heterogeneous landscape

Epstein Barr virus (EBV) has a large DNA genome assumed to be stable, but also subject to mutational processes such as nucleotide substitution and recombination, the latter explored to a lesser extent. Moreover, differences in the extent of recombination events across herpes sub-families were recently reported. Given the relevance of recombination in viral evolution and its possible impact in pathogenesis, we aimed to fully characterize and quantify its extension in all available EBV complete genome by assessing global and local recombination rate values (⍴/bp). Our results provide the first EBV recombination map based on recombination rates assessment, both at a global and gene by gene level, where the mean value for the entire genome was 0.035 (HPDI 0.020-0.062) ⍴/bp. We quantified how this evolutionary process changes along the EBV genome, and proved it to be non-homogeneous, since regulatory regions depicted the lowest recombination rate values while repetitive regions the highest signal. Moreover, GC content rich regions seem not to be linked to high recombination rates as previously reported. At an intragenic level, four genes (EBNA3C, EBNA3B, BRRF2 and BBLF2-BBLF3) presented a recombination rate above genome average. We specifically quantified the signal strength among different recombination-initiators previously described features and concluded that those which elicited the greatest amount of changes in ⍴/bp, TGGAG and CCCAG, were two well characterized recombination inducing motifs in eukaryotic cells. Strikingly, although TGGAG was not the most frequently detected DNA motif across the EBV genome (697 hits), it still induced a significantly greater proportion of initiation events (0.025 events/hits) than other more represented motifs, p-value = 0.04; one tailed proportion test. Present results support the idea that diversity and evolution of herpesviruses are impacted by mechanisms, such as recombination, which extends beyond the usual consideration of point mutations.

Copy number variation of functional RBMY1 is associated with sperm motility: an azoospermia factor-linked candidate for asthenozoospermia

STUDY QUESTION

What is the influence of copy number variation (CNV) in functional RNA binding motif protein Y-linked family 1 (RBMY1) on spermatogenic phenotypes?

SUMMARY ANSWER

The RBMY1 functional copy dosage is positively correlated with sperm motility, and dosage insufficiency is an independent risk factor for asthenozoospermia.

WHAT IS KNOWN ALREADY

RBMY1, a multi-copy gene expressed exclusively in the adult testis, is one of the most important candidates for male infertility in the azoospermia factor (AZF) region of the Y-chromosome. RBMY1 encodes an RNA-binding protein that serves as a pre-mRNA splicing regulator during spermatogenesis, and male mice deficient in Rbmy are sterile.

STUDY DESIGN, SIZE, DURATION

A total of 3127 adult males were recruited from 2009 to 2016; of this group, the dosage of RBMY1 functional copy were investigated in 486 fertile males. In the remaining 2641 males with known spermatogenesis status, 1070 Y-chromosome haplogroup (Y-hg) O3* or O3e carriers without chromosomal aberration or known AZF structure mutations responsible for spermatogenic impairment, including 506 men with normozoospermia and 564 men with oligozoospermia or/and asthenozoospermia, were screened, and the RBMY1 functional copy dosage and copy conversion were determined to explore their associations with sperm phenotypes. The correlation between RBMY1 dosage and its mRNA level or RBMY1 protein level and the correlation between sperm RBMY1 level and motility were analysed in 15 testis tissue samples and eight semen samples. Ten additional semen samples were used to confirm the subcellular localization of RBMY1 in individual sperm.

PARTICIPANTS/MATERIALS, SETTING, METHODS

All the Han volunteers donating whole blood, semen and testis tissue were from southwest China. RBMY1 copy number, copy conversion, mRNA/protein amount and protein location in sperm were detected using the AccuCopy® assay method, paralog ratio test, quantitative PCR, western blotting and immunofluorescence staining methods, respectively.

MAIN RESULTS AND THE ROLE OF CHANCE

This study identified Y-hg-independent CNV of functional RBMY1 in the enrolled population. A difference in the distribution of RBMY1 copy number was observed between the group with normal sperm motility and the group with asthenozoospermia. A positive correlation between the RBMY1 copy dosage and sperm motility was identified, and the males with fewer than six copies of RBMY1 showed an elevated risk for asthenozoospermia relative to those with six RBMY1 copies, the most common dosage in the population. The RBMY1 copy dosage was positively correlated with its mRNA and protein level in the testis. Sperm with high motility were found to carry more RBMY1 protein than those with relatively low motility. The RBMY1 protein was confirmed to predominantly localize in the neck and mid-piece region of sperm as well as the principal piece of the sperm tail. Our population study completes a chain of evidence suggesting that RBMY1 influences the susceptibility of males to asthenozoospermia by modulating sperm motility.

LIMITATIONS REASONS FOR CAUTION

High sequence similarity between the RBMY1 functional copies and a large number of pseudogenes potentially reduces the accuracy of the copy number detection. The mechanism underlying the CNV in RBMY1 is still unclear, and the effect of the structural variations in the RBMY1 copy cluster on the copy dosage of other protein-coding genes located in the region cannot be excluded, which may potentially bias our observations.

WIDER IMPLICATIONS OF THE FINDINGS

Asthenozoospermia is a multi-factor complex disease with a limited number of proven susceptibility genes. This study identified a novel genomic candidate independently contributing to the condition, enriching our understanding of the role of AZF-linked genes in male reproduction. Our finding provides insight into the physiological and pathological characteristics of RBMY1 in terms of sperm motility, supplies persuasive evidence of the significance of RBMY1 copy number analysis in the clinical counselling of male infertility resulting from asthenozoospermia.

STUDY FUNDING/COMPETING INTEREST(S)

This work was funded by the National Natural Science Foundation of China (Nos. 81370748 and 30971598). The authors have no conflicts of interest.

Heterogeneity of the Epstein-Barr Virus (EBV) Major Internal Repeat Reveals Evolutionary Mechanisms of EBV and a Functional Defect in the Prototype EBV Strain B95-8

Epstein-Barr virus (EBV) is a ubiquitous pathogen of humans that can cause several types of lymphoma and carcinoma. Like other herpesviruses, EBV has diversified through both coevolution with its host and genetic exchange between virus strains. Sequence analysis of the EBV genome is unusually challenging because of the large number and lengths of repeat regions within the virus. Here we describe the sequence assembly and analysis of the large internal repeat 1 of EBV (IR1; also known as the BamW repeats) for more than 70 strains. The diversity of the latency protein EBV nuclear antigen leader protein (EBNA-LP) resides predominantly within the exons downstream of IR1. The integrity of the putative BWRF1 open reading frame (ORF) is retained in over 80% of strains, and deletions truncating IR1 always spare BWRF1. Conserved regions include the IR1 latency promoter (Wp) and one zone upstream of and two within BWRF1. IR1 is heterogeneous in 70% of strains, and this heterogeneity arises from sequence exchange between strains as well as from spontaneous mutation, with interstrain recombination being more common in tumor-derived viruses. This genetic exchange often incorporates regions of <1 kb, and allelic gene conversion changes the frequency of small regions within the repeat but not close to the flanks. These observations suggest that IR1-and, by extension, EBV-diversifies through both recombination and breakpoint repair, while concerted evolution of IR1 is driven by gene conversion of small regions. Finally, the prototype EBV strain B95-8 contains four nonconsensus variants within a single IR1 repeat unit, including a stop codon in the EBNA-LP gene. Repairing IR1 improves EBNA-LP levels and the quality of transformation by the B95-8 bacterial artificial chromosome (BAC).IMPORTANCE Epstein-Barr virus (EBV) infects the majority of the world population but causes illness in only a small minority of people. Nevertheless, over 1% of cancers worldwide are attributable to EBV. Recent sequencing projects investigating virus diversity to see if different strains have different disease impacts have excluded regions of repeating sequence, as they are more technically challenging. Here we analyze the sequence of the largest repeat in EBV (IR1). We first characterized the variations in protein sequences encoded across IR1. In studying variations within the repeat of each strain, we identified a mutation in the main laboratory strain of EBV that impairs virus function, and we suggest that tumor-associated viruses may be more likely to contain DNA mixed from two strains. The patterns of this mixing suggest that sequences can spread between strains (and also within the repeat) by copying sequence from another strain (or repeat unit) to repair DNA damage.

Oyster RNA-seq Data Support the Development of Malacoherpesviridae Genomics

The family of double-stranded DNA (dsDNA) Malacoherpesviridae includes viruses able to infect marine mollusks and detrimental for worldwide aquaculture production. Due to fast-occurring mortality and a lack of permissive cell lines, the available data on the few known Malacoherpesviridae provide only partial support for the study of molecular virus features, life cycle, and evolutionary history. Following thorough data mining of bivalve and gastropod RNA-seq experiments, we used more than five million Malacoherpesviridae reads to improve the annotation of viral genomes and to characterize viral InDels, nucleotide stretches, and SNPs. Both genome and protein domain analyses confirmed the evolutionary diversification and gene uniqueness of known Malacoherpesviridae. However, the presence of Malacoherpesviridae-like sequences integrated within genomes of phylogenetically distant invertebrates indicates broad diffusion of these viruses and indicates the need for confirmatory investigations. The manifest co-occurrence of OsHV-1 genotype variants in single RNA-seq samples of Crassostrea gigas provide further support for the Malacoherpesviridae diversification. In addition to simple sequence motifs inter-punctuating viral ORFs, recombination-inducing sequences were found to be enriched in the OsHV-1 and AbHV1-AUS genomes. Finally, the highly correlated expression of most viral ORFs in multiple oyster samples is consistent with the burst of viral proteins during the lytic phase.

A novel homozygous complex deletion in CFTR caused cystic fibrosis in a Chinese patient

Cystic fibrosis (CF) is the most frequent lethal genetic disorder among Caucasians, but is considered to be a very rare disease in Chinese population. Here, we present an 11-year-old Chinese CF patient with disseminated bronchiectasis and salty sweat, for whom exon sequencing followed by multiplex ligation-dependent probe amplification analysis of the CFTR gene was applied for mutation screening. A homozygous deletion involving exon 20 of CFTR was observed in the patient's genome. Molecular characterization of the breakpoints indicated that both alleles of this locus had an identical novel complex rearrangement (c.3140-454_c.3367+249del931ins13, p.R1048_G1123del), leading to an in-frame removal of 76 amino acid residues in the second transmembrane domains of the CFTR protein. Although a same haplotype containing this complex rearrangement was observed on both of the maternal and paternal alleles, the parents denied any blood relationship as far as they know. Genome-wide homozygosity mapping was performed through SNP microarray and only a single homozygous interval of ~14.1 Mb at chromosome 7 containing the CFTR gene was observed, indicating the possible origin of the deletion from a common ancestor many generations ago. This study expands the mutation spectrum of CFTR in patients of Chinese origin and further emphasizes the necessity of MLPA analysis in mutation screening for CF patients.

Analysis of Multiple Families With Single Individuals Affected by Pseudohypoparathyroidism Type Ib (PHP1B) Reveals Only One Novel Maternally Inherited GNAS Deletion

Proximal tubular resistance to parathyroid hormone (PTH) resulting in hypocalcemia and hyperphosphatemia are preeminent abnormalities in pseudohypoparathyroidism type Ib (PHP1B), but resistance toward other hormones as well as variable features of Albright's Hereditary Osteodystrophy (AHO) can occur also. Genomic DNA from PHP1B patients shows epigenetic changes at one or multiple differentially methylated regions (DMRs) within GNAS, the gene encoding Gαs and splice variants thereof. In the autosomal dominant disease variant, these methylation abnormalities are caused by deletions in STX16 or GNAS on the maternal allele. The molecular defect(s) leading to sporadic PHP1B (sporPHP1B) remains in most cases unknown and we therefore analyzed 60 sporPHP1B patients and available family members by microsatellite markers, single nucleotide polymorphisms (SNPs), multiplex ligation-dependent probe amplification (MLPA), and methylation-specific MLPA (MS-MLPA). All investigated cases revealed broad GNAS methylation changes, but no evidence for inheritance of two paternal chromosome 20q alleles. Some patients with partial epigenetic modifications in DNA from peripheral blood cells showed more complete GNAS methylation changes when testing their immortalized lymphoblastoid cells. Analysis of siblings and children of sporPHP1B patients provided no evidence for an abnormal mineral ion regulation and no changes in GNAS methylation. Only one patient revealed, based on MLPA and microsatellite analyses, evidence for an allelic loss, which resulted in the discovery of two adjacent, maternally inherited deletions (37,597 and 1427 bp, respectively) that remove the area between GNAS antisense exons 3 and 5, including exon NESP. Our findings thus emphasize that the region comprising antisense exons 3 and 4 is required for establishing all maternal GNAS methylation imprints. The genetic defect(s) leading in sporPHP1B to epigenetic GNAS changes and thus PTH-resistance remains unknown, but it seems unlikely that this disease variant is caused by heterozygous inherited or de novo mutations involving GNAS.

Concurrent nucleotide substitution mutations in the human genome are characterized by a significantly decreased transition/transversion ratio

There is accumulating evidence that the number of multiple-nucleotide substitutions (MNS) occurring in closely spaced sites in eukaryotic genomes is significantly higher than would be predicted from the random accumulation of independently generated single-nucleotide substitutions (SNS). Although this excess can in principle be accounted for by the concept of transient hypermutability, a general mutational signature of concurrent MNS mutations has not so far been evident. Employing a dataset (N = 449) of "concurrent" double MNS mutations causing human inherited disease, we have identified just such a mutational signature: concurrently generated double MNS mutations exhibit a >twofold lower transition/transversion ratio (termed RTs/Tv ) than independently generated de novo SNS mutations (<0.80 vs. 2.10; P = 2.69 × 10(-14) ). We replicated this novel finding through a similar analysis employing two double MNS variant datasets with differing abundances of concurrent events (150,521 variants with both substitutions on the same haplotypic lineage vs. 94,875 variants whose component substitutions were on different haplotypic lineages) plus 5,430,874 SNS variants, all being derived from the whole-genome sequencing of seven Chinese individuals. Evaluation of the newly observed mutational signature in diverse contexts provides solid support for the postulated role of translesion synthesis DNA polymerases in transient hypermutability.

A Role for Non-B DNA Forming Sequences in Mediating Microlesions Causing Human Inherited Disease

Missense/nonsense mutations and microdeletions/microinsertions (<21 bp) represent ∼ 76% of all mutations causing human inherited disease, and their occurrence has been associated with sequence motifs (direct, inverted, and mirror repeats; G-quartets) capable of adopting non-B DNA structures. We found that a significant proportion (∼ 21%) of both microdeletions and microinsertions occur within direct repeats, and are explicable by slipped misalignment. A novel mutational mechanism, DNA triplex formation followed by DNA repair, may explain ∼ 5% of microdeletions and microinsertions at mirror repeats. Further, G-quartets, direct, and inverted repeats also appear to play a prominent role in mediating missense mutations, whereas only direct and inverted repeats mediate nonsense mutations. We suggest a mutational mechanism involving slipped strand mispairing, slipped structure formation, and DNA repair, to explain ∼ 15% of missense and ∼ 12% of nonsense mutations yielding perfect direct repeats from imperfect repeats, or the extension of existing direct repeats. Similar proportions of missense and nonsense mutations were explicable by hairpin/loop formation and DNA repair, yielding perfect inverted repeats from imperfect repeats. We also propose a model for single base-pair substitution based on one-electron oxidation reactions at G-quadruplex DNA. Overall, the proposed mechanisms provide support for a role for non-B DNA structures in human gene mutagenesis.

High-resolution breakpoint analysis provides evidence for the sequence-directed nature of genome rearrangements in hereditary disorders

Although most of the pertinent data on the sequence-directed processes leading to genome rearrangements (GRs) have come from studies on somatic tissues, little is known about GRs in the germ line of patients with hereditary disorders. This study aims at identifying DNA motifs and higher order structures of genome architecture, which can result in losses and gains of genetic material in the germ line. We first identified candidate motifs by studying 112 pathogenic germ-line GRs in hereditary colorectal cancer patients, and subsequently created an algorithm, termed recombination type ratio, which correctly predicts the propensity of rearrangements with respect to homologous versus nonhomologous recombination events.

Potential non-B DNA regions in the human genome are associated with higher rates of nucleotide mutation and expression variation

While individual non-B DNA structures have been shown to impact gene expression, their broad regulatory role remains elusive. We utilized genomic variants and expression quantitative trait loci (eQTL) data to analyze genome-wide variation propensities of potential non-B DNA regions and their relation to gene expression. Independent of genomic location, these regions were enriched in nucleotide variants. Our results are consistent with previously observed mutagenic properties of these regions and counter a previous study concluding that G-quadruplex regions have a reduced frequency of variants. While such mutagenicity might undermine functionality of these elements, we identified in potential non-B DNA regions a signature of negative selection. Yet, we found a depletion of eQTL-associated variants in potential non-B DNA regions, opposite to what might be expected from their proposed regulatory role. However, we also observed that genes downstream of potential non-B DNA regions showed higher expression variation between individuals. This coupling between mutagenicity and tolerance for expression variability of downstream genes may be a result of evolutionary adaptation, which allows reconciling mutagenicity of non-B DNA structures with their location in functionally important regions and their potential regulatory role.

Characterization of a new V gene replacement in the absence of activation-induced cytidine deaminase and its contribution to human B-cell receptor diversity

In B cells, B-cell receptor (BCR) immunoglobulin revision is a common route for modifying unwanted antibody specificities via a mechanism called VH replacement. This in vivo process, mostly affecting heavy-chain rearrangement, involves the replacement of all or part of a previously rearranged IGHV gene with another germline IGHV gene located upstream. Two different mechanisms of IGHV replacement have been reported: type 1, involving the recombination activating genes complex and requiring a framework region 3 internal recombination signal; and type 2, involving an unidentified mechanism different from that of type 1. In the case of light-chain loci, BCR immunoglobulin editing ensures that a second V-J rearrangement occurs. This helps to maintain tolerance, by generating a novel BCR with a new antigenic specificity. We report that human B cells can, surprisingly, undergo type 2 replacement associated with κ light-chain rearrangements. The de novo IGKV-IGKJ products result from the partial replacement of a previously rearranged IGKV gene by a new germline IGKV gene, in-frame and without deletion or addition of nucleotides. There are wrcy/rgyw motifs at the 'IGKV donor-IGKV recipient chimera junction' as described for type 2 IGHV replacement, but activation-induced cytidine deaminase (AID) expression was not detected. This unusual mechanism of homologous recombination seems to be a variant of gene conversion-like recombination, which does not require AID. The recombination phenomenon described here provides new insight into immunoglobulin locus recombination and BCR immunoglobulin repertoire diversity.

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.

DHX9 helicase is involved in preventing genomic instability induced by alternatively structured DNA in human cells

Sequences that have the capacity to adopt alternative (i.e. non-B) DNA structures in the human genome have been implicated in stimulating genomic instability. Previously, we found that a naturally occurring intra-molecular triplex (H-DNA) caused genetic instability in mammals largely in the form of DNA double-strand breaks. Thus, it is of interest to determine the mechanism(s) involved in processing H-DNA. Recently, we demonstrated that human DHX9 helicase preferentially unwinds inter-molecular triplex DNA in vitro. Herein, we used a mutation-reporter system containing H-DNA to examine the relevance of DHX9 activity on naturally occurring H-DNA structures in human cells. We found that H-DNA significantly increased mutagenesis in small-interfering siRNA-treated, DHX9-depleted cells, affecting mostly deletions. Moreover, DHX9 associated with H-DNA in the context of supercoiled plasmids. To further investigate the role of DHX9 in the recognition/processing of H-DNA, we performed binding assays in vitro and chromatin immunoprecipitation assays in U2OS cells. DHX9 recognized H-DNA, as evidenced by its binding to the H-DNA structure and enrichment at the H-DNA region compared with a control region in human cells. These composite data implicate DHX9 in processing H-DNA structures in vivo and support its role in the overall maintenance of genomic stability at sites of alternatively structured DNA.

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.

Patterns and mutational signatures of tandem base substitutions causing human inherited disease

Tandem base substitutions (TBSs) are multiple mutations that comprise two or more contiguous nucleotide substitutions without any net gain or loss of bases. They have recently become recognized as a distinct category of human genomic variant. However, their role in causing human inherited disease so far has not been studied methodically. Here, using data from the Human Gene Mutation Database (http://www.hgmd.org), we identified 477 events to be TBSs (doublets, 448; triplets, 16; and quadruplets to octuplets, 13). A comprehensive sequence pattern and context analysis implied the likely fundamental importance of translesion synthesis (TLS) DNA polymerases in generating these diverse TBSs but revealed that TLS polymerases may operate differently in generating TBSs of ≤ 3 bases (bypass of endogenous DNA lesions) than those of ≥ 4 bases (serial replication slippage). Moreover, GC was found to be the most frequently affected dinucleotide with GC/GC>AA/TT being the most frequent double TBS. Comparison with cancer genome mutational spectra allowed us to conclude that human germline TBSs arise predominantly through the action of endogenous mechanisms of mutagenesis rather than through exposure to exogenous mutagens. Finally, the rates of double and triple TBSs were estimated to be 0.2-1.2 × 10(-10) and 0.8-4.8 × 10(-12) per base per generation, respectively.

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.

An evolutionary view of the mechanism for immune and genome diversity

An ortholog of activation-induced cytidine deaminase (AID) was, evolutionarily, the first enzyme to generate acquired immune diversity by catalyzing gene conversion and probably somatic hypermutation (SHM). AID began to mediate class switch recombination (CSR) only after the evolution of frogs. Recent studies revealed that the mechanisms for generating immune and genetic diversity share several critical features. Meiotic recombination, V(D)J recombination, CSR, and SHM all require H3K4 trimethyl histone modification to specify the target DNA. Genetic instability related to dinucleotide or triplet repeats depends on DNA cleavage by topoisomerase 1, which also initiates DNA cleavage in both SHM and CSR. These similarities suggest that AID hijacked the basic mechanism for genome instability when AID evolved in jawless fish. Thus, the risk of introducing genome instability into nonimmunoglobulin loci is unavoidable but tolerable compared with the advantage conferred on the host of being protected against pathogens by the enormous Ig diversification.

Evolutionary comparison of the mechanism of DNA cleavage with respect to immune diversity and genomic instability

It is generally assumed that the genetic mechanism for immune diversity is unique and distinct from that for general genome diversity, in part because of the high efficiency and strict regulation of immune diversity. This expectation was partially met by the discovery of RAG1 and -2, which catalyze V(D)J recombination to generate the immune repertoire of B and T lymphocyte receptors. RAG1 and -2 were later shown to be derived from a transposon. On the other hand, activation-induced cytidine deaminase (AID), which mediates both somatic hypermutation (SHM) and the class-switch recombination (CSR) of the immunoglobulin genes, evolved earlier than RAG1 and -2 in jawless vertebrates. This review compares immune diversity and general genome diversity from an evolutionary perspective, shedding light on the roles of DNA-cleaving enzymes and target recognition markers. This comparison revealed that AID-mediated SHM and CSR share the cleaving enzyme topoisomerase 1 with transcription-associated mutation (TAM) and triplet contraction, which is involved in many genetic diseases. These genome-altering events appear to target DNA with non-B structure, which is induced by the inefficient correction of the excessive supercoiling that is caused by active transcription. Furthermore, an epigenetic modification on chromatin (histone H3K4 trimethylation) is used as a mark for DNA cleavage sites in meiotic recombination, V(D)J recombination, CSR, and SHM. We conclude that acquired immune diversity evolved via the appearance of an AID orthologue that utilized a preexisting mechanism for genomic instability, such as TAM.

Gene Conversion-Like Events in the Diversification of Human Rearranged IGHV3-23*01 Gene Sequences

Gene conversion (GCV), a mechanism mediated by activation-induced cytidine deaminase (AID) is well established as a mechanism of immunoglobulin diversification in a few species. However, definitive evidence of GCV-like events in human immunoglobulin genes is scarce. The lack of evidence of GCV in human rearranged immunoglobulin gene sequences is puzzling given the presence of highly similar germline donors and the presence of all the enzymatic machinery required for GCV. In this study, we undertook a computational analysis of rearranged IGHV3-23(*)01 gene sequences from common variable immunodeficiency (CVID) patients, AID-deficient patients, and healthy individuals to survey "GCV-like" activities. We analyzed rearranged IGHV3-23(*)01 gene sequences obtained from total PBMC RNA and single-cell polymerase chain reaction of individual B cell lysates. Our search identified strong evidence of GCV-like activity. We observed that GCV-like tracts are flanked by AID hotspot motifs. Structural modeling of IGHV3-23(*)01 gene sequence revealed that hypermutable bases flanking GCV-like tracts are in the single stranded DNA (ssDNA) of stable stem-loop structures (SLSs). ssDNA is inherently fragile and also an optimal target for AID. We speculate that GCV could have been initiated by the targeting of hypermutable bases in ssDNA state in stable SLSs, plausibly by AID. We have observed that the frequency of GCV-like events is significantly higher in rearranged IGHV3-23-(*)01 sequences from healthy individuals compared to that of CVID patients. We did not observe GCV-like events in rearranged IGHV3-23-(*)01 sequences from AID-deficient patients. GCV, unlike somatic hypermutation (SHM), can result in multiple base substitutions that can alter many amino acids. The extensive changes in antibody affinity by GCV-like events would be instrumental in protecting humans against pathogens that diversify their genome by antigenic shift.

Multiple cellular mechanisms prevent chromosomal rearrangements involving repetitive DNA

Repetitive DNA is present in the eukaryotic genome in the form of segmental duplications, tandem and interspersed repeats, and satellites. Repetitive sequences can be beneficial by serving specific cellular functions (e.g. centromeric and telomeric DNA) and by providing a rapid means for adaptive evolution. However, such elements are also substrates for deleterious chromosomal rearrangements that affect fitness and promote human disease. Recent studies analyzing the role of nuclear organization in DNA repair and factors that suppress non-allelic homologous recombination (NAHR) have provided insights into how genome stability is maintained in eukaryotes. In this review, we outline the types of repetitive sequences seen in eukaryotic genomes and how recombination mechanisms are regulated at the DNA sequence, cell organization, chromatin structure, and cell cycle control levels to prevent chromosomal rearrangements involving these sequences.

Interlocus gene conversion events introduce deleterious mutations into at least 1% of human genes associated with inherited disease

Establishing the molecular basis of DNA mutations that cause inherited disease is of fundamental importance to understanding the origin, nature, and clinical sequelae of genetic disorders in humans. The majority of disease-associated mutations constitute single-base substitutions and short deletions and/or insertions resulting from DNA replication errors and the repair of damaged bases. However, pathological mutations can also be introduced by nonreciprocal recombination events between paralogous sequences, a phenomenon known as interlocus gene conversion (IGC). IGC events have thus far been linked to pathology in more than 20 human genes. However, the large number of duplicated gene sequences in the human genome implies that many more disease-associated mutations could originate via IGC. Here, we have used a genome-wide computational approach to identify disease-associated mutations derived from IGC events. Our approach revealed hundreds of known pathological mutations that could have been caused by IGC. Further, we identified several dozen high-confidence cases of inherited disease mutations resulting from IGC in ∼1% of all genes analyzed. About half of the donor sequences associated with such mutations are functional paralogous genes, suggesting that epistatic interactions or differential expression patterns will determine the impact upon fitness of specific substitutions between duplicated genes. In addition, we identified thousands of hitherto undescribed and potentially deleterious mutations that could arise via IGC. Our findings reveal the extent of the impact of interlocus gene conversion upon the spectrum of human inherited disease.

Evolution of siglec-11 and siglec-16 genes in hominins

We previously reported a human-specific gene conversion of SIGLEC11 by an adjacent paralogous pseudogene (SIGLEC16P), generating a uniquely human form of the Siglec-11 protein, which is expressed in the human brain. Here, we show that Siglec-11 is expressed exclusively in microglia in all human brains studied-a finding of potential relevance to brain evolution, as microglia modulate neuronal survival, and Siglec-11 recruits SHP-1, a tyrosine phosphatase that modulates microglial biology. Following the recent finding of a functional SIGLEC16 allele in human populations, further analysis of the human SIGLEC11 and SIGLEC16/P sequences revealed an unusual series of gene conversion events between two loci. Two tandem and likely simultaneous gene conversions occurred from SIGLEC16P to SIGLEC11 with a potentially deleterious intervening short segment happening to be excluded. One of the conversion events also changed the 5' untranslated sequence, altering predicted transcription factor binding sites. Both of the gene conversions have been dated to ~1-1.2 Ma, after the emergence of the genus Homo, but prior to the emergence of the common ancestor of Denisovans and modern humans about 800,000 years ago, thus suggesting involvement in later stages of hominin brain evolution. In keeping with this, recombinant soluble Siglec-11 binds ligands in the human brain. We also address a second-round more recent gene conversion from SIGLEC11 to SIGLEC16, with the latter showing an allele frequency of ~0.1-0.3 in a worldwide population study. Initial pseudogenization of SIGLEC16 was estimated to occur at least 3 Ma, which thus preceded the gene conversion of SIGLEC11 by SIGLEC16P. As gene conversion usually disrupts the converted gene, the fact that ORFs of hSIGLEC11 and hSIGLEC16 have been maintained after an unusual series of very complex gene conversion events suggests that these events may have been subject to hominin-specific selection forces.

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.

The consequences of structural genomic alterations in humans: genomic disorders, genomic instability and cancer

Over the last decade or so, sophisticated technological advances in array-based genomics have firmly established the contribution of structural alterations in the human genome to a variety of complex developmental disorders, and also to diseases such as cancer. In fact, multiple 'novel' disorders have been identified as a direct consequence of these advances. Our understanding of the molecular events leading to the generation of these structural alterations is also expanding. Many of the models proposed to explain these complex rearrangements involve DNA breakage and the coordinated action of DNA replication, repair and recombination machinery. Here, and within the context of Genomic Disorders, we will briefly overview the principal models currently invoked to explain these chromosomal rearrangements, including Non-Allelic Homologous Recombination (NAHR), Fork Stalling Template Switching (FoSTeS), Microhomology Mediated Break-Induced Repair (MMBIR) and Breakage-fusion-bridge cycle (BFB). We will also discuss an unanticipated consequence of certain copy number variations (CNVs) whereby the CNVs potentially compromise fundamental processes controlling genomic stability including DNA replication and the DNA damage response. We will illustrate these using specific examples including Genomic Disorders (DiGeorge/Veleocardiofacial syndrome, HSA21 segmental aneuploidy and rec (3) syndrome) and cell-based model systems. Finally, we will review some of the recent exciting developments surrounding specific CNVs and their contribution to cancer development as well as the latest model for cancer genome rearrangement; 'chromothripsis'.

TaqMan-based real-time polymerase chain reaction for detection of FUT2 copy number variations: identification of novel Alu-mediated deletion

BACKGROUND

The human FUT2 locus, which encodes a secretor-type α(1,2)fucosyltransferase, is known to be highly polymorphic. In addition to many single-nucleotide polymorphisms, three recombination alleles with a deletion of complete or partial FUT2 coding region have been reported.

STUDY DESIGN AND METHODS

To detect copy number variations (CNVs) of the FUT2 gene including three recombinant alleles by a high-throughput system, we developed a triplex TaqMan real-time polymerase chain reaction (PCR) method. The relative number of copies of two regions of the FUT2 gene, the 5' flanking (FUT2-5') and FUT2-promoter (Prom) regions, were determined by comparing the number of threshold cycles (Ct) to those of the albumin gene (ALB) as the internal control (ΔCt).

RESULTS

The mean 2(-ΔΔCt) values (FUT2-5'/ALB or Prom/ALB) obtained from 237 samples with known FUT2 copy numbers clearly differentiated two nonoverlapping intervals that corresponded to the one-copy-number samples ranging from 0.42 to 0.59 and two-copy-number samples ranging from 0.81 to 1.19; no FUT2-5' signal for recombination alleles was detected in homozygotes. Using this assay, we detected an individual in a Chinese population with a loss of one copy of the FUT2-5' region resulting from a novel Alu-mediated FUT2 deletion (approx. 4 kb).

CONCLUSIONS

The TaqMan real-time PCR method was able to detect the number of copies of FUT2 and distinguish different kinds of known CNVs. This system is robust, fast, and suitable for high-throughput analysis.

The Rate and Tract Length of Gene Conversion between Duplicated Genes

Interlocus gene conversion occurs such that a certain length of DNA fragment is non-reciprocally transferred (copied and pasted) between paralogous regions. To understand the rate and tract length of gene conversion, there are two major approaches. One is based on mutation-accumulation experiments, and the other uses natural DNA sequence variation. In this review, we overview the two major approaches and discuss their advantages and disadvantages. In addition, to demonstrate the importance of statistical analysis of empirical and evolutionary data for estimating tract length, we apply a maximum likelihood method to several data sets.

Non-B DNA-forming sequences and WRN deficiency independently increase the frequency of base substitution in human cells

Although alternative DNA secondary structures (non-B DNA) can induce genomic rearrangements, their associated mutational spectra remain largely unknown. The helicase activity of WRN, which is absent in the human progeroid Werner syndrome, is thought to counteract this genomic instability. We determined non-B DNA-induced mutation frequencies and spectra in human U2OS osteosarcoma cells and assessed the role of WRN in isogenic knockdown (WRN-KD) cells using a supF gene mutation reporter system flanked by triplex- or Z-DNA-forming sequences. Although both non-B DNA and WRN-KD served to increase the mutation frequency, the increase afforded by WRN-KD was independent of DNA structure despite the fact that purified WRN helicase was found to resolve these structures in vitro. In U2OS cells, ∼70% of mutations comprised single-base substitutions, mostly at G·C base-pairs, with the remaining ∼30% being microdeletions. The number of mutations at G·C base-pairs in the context of NGNN/NNCN sequences correlated well with predicted free energies of base stacking and ionization potentials, suggesting a possible origin via oxidation reactions involving electron loss and subsequent electron transfer (hole migration) between neighboring bases. A set of ∼40,000 somatic mutations at G·C base pairs identified in a lung cancer genome exhibited similar correlations, implying that hole migration may also be involved. We conclude that alternative DNA conformations, WRN deficiency and lung tumorigenesis may all serve to increase the mutation rate by promoting, through diverse pathways, oxidation reactions that perturb the electron orbitals of neighboring bases. It follows that such "hole migration" is likely to play a much more widespread role in mutagenesis than previously anticipated.

Comparative analysis of germline and somatic microlesion mutational spectra in 17 human tumor suppressor genes

Mutations associated with tumorigenesis may either arise somatically or can be inherited through the germline. We performed a comparison of somatic, germline, shared (found in both soma and germline) and somatic recurrent mutational spectra for 17 human tumor suppressor genes, which focused upon missense single base-pair substitutions and microdeletions/microinsertions. Somatic and germline mutational spectra were similar in relation to C.G>T.A transitions but differed with respect to the frequency of A.T>G.C, A.T>T.A, and C.G>A.T substitutions. Shared missense mutations were characterized by higher mutability rates, greater physicochemical differences between wild-type and mutant residues, and a tendency to occur in evolutionarily conserved residues and within CpG/CpHpG oligonucleotides. Mononucleotide runs (≥4 bp) were identified as hotspots for shared microdeletions/microinsertions. Both germline and somatic microdeletions/microinsertions were found to be significantly overrepresented within the "indel-hotspot" motif, GTAAGT. Using a naïve Bayes' classifier trained to discriminate between five missense mutation groups, 63% of mutations in our dataset were on average correctly recognized. Applying this classifier to an independent dataset of probable driver mutations, we concluded that ∼50% of these somatic missense mutations possess features consistent with their being either shared or recurrent, suggesting that a disproportionate number of such lesions are likely to be drivers of tumorigenesis.

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.

Gene conversion in human genetic disease

Gene conversion is a specific type of homologous recombination that involves the unidirectional transfer of genetic material from a ‘donor’ sequence to a highly homologous ‘acceptor’. We have recently reviewed the molecular mechanisms underlying gene conversion, explored the key part that this process has played in fashioning extant human genes, and performed a meta-analysis of gene-conversion events known to have caused human genetic disease. Here we shall briefly summarize some of the latest developments in the study of pathogenic gene conversion events, including (i) the emerging idea of minimal efficient sequence homology (MESH) for homologous recombination, (ii) the local DNA sequence features that appear to predispose to gene conversion, (iii) a mechanistic comparison of gene conversion and transient hypermutability, and (iv) recently reported examples of pathogenic gene conversion events.

Telomere healing following DNA polymerase arrest-induced breakages is likely the main mechanism generating chromosome 4p terminal deletions

Constitutional developmental disorders are frequently caused by terminal chromosomal deletions. The mechanisms and/or architectural features that might underlie those chromosome breakages remain largely unexplored. Because telomeres are the vital DNA protein complexes stabilizing linear chromosomes against chromosome degradation, fusion, and incomplete replication, those terminal-deleted chromosomes acquired new telomeres either by telomere healing or by telomere capture. To unravel the mechanisms leading to chromosomal breakage and healing, we sequenced nine chromosome 4p terminal deletion boundaries. A computational analysis of the breakpoint flanking region, including 12 previously published pure terminal breakage sites, was performed in order to identify architectural features that might be involved in this process. All terminal 4p truncations were likely stabilized by telomerase-mediated telomere healing. In the majority of breakpoints multiple genetic elements have a potential to induce secondary structures and an enrichment in replication stalling site motifs were identified. These findings suggest DNA replication stalling-induced chromosome breakage during early development is the first mechanistic step leading toward terminal deletion syndromes.

dbCRID: a database of chromosomal rearrangements in human diseases

Chromosomal rearrangement (CR) events result from abnormal breaking and rejoining of the DNA molecules, or from crossing-over between repetitive DNA sequences, and they are involved in many tumor and non-tumor diseases. Investigations of disease-associated CR events can not only lead to important discoveries about DNA breakage and repair mechanisms, but also offer important clues about the pathologic causes and the diagnostic/therapeutic targets of these diseases. We have developed a database of Chromosomal Rearrangements In Diseases (dbCRID, http://dbCRID.biolead.org), a comprehensive database of human CR events and their associated diseases. For each reported CR event, dbCRID documents the type of the event, the disease or symptoms associated, and--when possible--detailed information about the CR event including precise breakpoint positions, junction sequences, genes and gene regions disrupted and experimental techniques applied to discover/analyze the CR event. With 2643 records of disease-associated CR events curated from 1172 original studies, dbCRID is a comprehensive and dynamic resource useful for studying DNA breakage and repair mechanisms, and for analyzing the genetic basis of human tumor and non-tumor diseases.

Chimeric CYP21A1P/CYP21A2 genes identified in Czech patients with congenital adrenal hyperplasia

Congenital adrenal hyperplasia (CAH) comprises a group of autosomal recessive disorders caused by an enzymatic deficiency which impairs the biosynthesis of cortisol and, in the majority of severe cases, also the biosynthesis of aldosterone. Approximately 95% of all CAH cases are caused by mutations in the steroid 21-hydroxylase gene (CYP21A2). The CYP21A2 gene and its inactive pseudogene (CYP21A1P) are located within the HLA class III region of the major histocompatibility complex (MHC) locus on chromosome 6p21.3. In this study, we describe chimeric CYP21A1P/CYP21A2 genes detected in our patients with 21-hydroxylase deficiency (21OHD). Chimeric CYP21A1P/CYP21A2 genes were present in 171 out of 508 mutated CYP21A2 alleles (33.8%). We detected four types of chimeric CYP21A1P/CYP21A2 genes: three of them have been described previously as CH-1, CH-3, CH-4, and one type is novel. The novel chimeric gene, termed CH-7, was detected in 21.4% of the mutant alleles. Possible causes of CYP21A1P/CYP21A2 formation are associated with 1) high recombination rate in the MHC locus, 2) high recombination rate between highly homologous genes and pseudogenes in the CYP21 gene area, and 3) the existence of chi-like sequences and repetitive minisatellite consensus sequences in CYP21A2 and CYP21A1P which play a role in promoting genetic recombination.

A novel third type of recurrent NF1 microdeletion mediated by nonallelic homologous recombination between LRRC37B-containing low-copy repeats in 17q11.2

Large microdeletions encompassing the neurofibromatosis type-1 (NF1) gene and its flanking regions at 17q11.2 belong to the group of genomic disorders caused by aberrant recombination between segmental duplications. The most common NF1 microdeletions (type-1) span 1.4-Mb and have breakpoints located within NF1-REPs A and C, low-copy repeats (LCRs) containing LRRC37-core duplicons. We have identified a novel type of recurrent NF1 deletion mediated by nonallelic homologous recombination (NAHR) between the highly homologous NF1-REPs B and C. The breakpoints of these approximately 1.0-Mb ("type-3") NF1 deletions were characterized at the DNA sequence level in three unrelated patients. Recombination regions, spanning 275, 180, and 109-bp, respectively, were identified within the LRRC37B-P paralogues of NF1-REPs B and C, and were found to contain sequences capable of non-B DNA formation. Both LCRs contain LRRC37-core duplicons, abundant and highly dynamic sequences in the human genome. NAHR between LRRC37-containing LCRs at 17q21.31 is known to have mediated the 970-kb polymorphic inversions of the MAPT-locus that occurred independently in different primate species, but also underlies the syndromes associated with recurrent 17q21.31 microdeletions and reciprocal microduplications. The novel NF1 microdeletions reported here provide further evidence for the unusually high recombinogenic potential of LRRC37-containing LCRs in the human genome.

Complete ascertainment of intragenic copy number mutations (CNMs) in the CFTR gene and its implications for CNM formation at other autosomal loci

Over the last 20 years since the discovery of the cystic fibrosis transmembrane conductance regulator (CFTR) gene, more than 1,600 different putatively pathological CFTR mutations have been identified. Until now, however, copy number mutations (CNMs) involving the CFTR gene have not been methodically analyzed, resulting almost certainly in the underascertainment of CFTR gene duplications compared with deletions. Here, high-resolution array comparative genomic hybridization (averaging one interrogating probe every 95 bp) was used to analyze the entire length of the CFTR gene (189 kb) in 233 cystic fibrosis chromosomes lacking conventional mutations. We succeeded in identifying five duplication CNMs that would otherwise have been refractory to analysis. Based upon findings from this and other studies, we propose that deletion and duplication CNMs in the human autosomal genome are likely to be generated in the proportion of approximately 2-3:1. We further postulate that intragenic gene duplication CNMs in other disease loci may have been routinely underascertained. Finally, our analysis of +/-20 bp flanking each of the 40 CFTR breakpoints characterized at the DNA sequence level provide support for the emerging concept that non-B DNA conformations in combination with specific sequence motifs predispose to both recurring and nonrecurring genomic rearrangements.

Biological roles of translin and translin-associated factor-X: RNA metabolism comes to the fore

Translin, and its binding partner protein TRAX (translin-associated factor-X) are a paralogous pair of conserved proteins, which have been implicated in a broad spectrum of biological activities, including cell growth regulation, mRNA processing, spermatogenesis, neuronal development/function, genome stability regulation and carcinogenesis, although their precise role in some of these processes remains unclear. Furthermore, translin (with or without TRAX) has nucleic-acid-binding activity and it is apparent that controlling nucleic acid metabolism and distribution are central to the biological role(s) of this protein and its partner TRAX. More recently, translin and TRAX have together been identified as enhancer components of an RNAi (RNA interference) pathway in at least one organism and this might provide critical insight into the biological roles of this enigmatic partnership. In the present review we discuss the biological and the biochemical properties of these proteins that indicate that they play a central and important role in eukaryotic cell biology.

Genomic rearrangements in inherited disease and cancer

Genomic rearrangements in inherited disease and cancer involve gross alterations of chromosomes or large chromosomal regions and can take the form of deletions, duplications, insertions, inversions or translocations. The characterization of a considerable number of rearrangement breakpoints has now been accomplished at the nucleotide sequence level, thereby providing an invaluable resource for the detailed study of the mutational mechanisms which underlie genomic recombination events. A better understanding of these mutational mechanisms is vital for improving the design of mutation detection strategies. At least five categories of mutational mechanism are known to give rise to genomic rearrangements: (i) homologous recombination including non-allelic homologous recombination (NAHR), gene conversion, single strand annealing (SSA) and break-induced replication (BIR), (ii) non-homologous end joining (NHEJ), (iii) microhomology-mediated replication-dependent recombination (MMRDR), (iv) long interspersed element-1 (LINE-1 or L1)-mediated retrotransposition and (v) telomere healing. Focussing on the first three of these general mechanisms, we compare and contrast their hallmark characteristics, and discuss the role of various local DNA sequence features (e.g. recombination-promoting motifs, repetitive sequences and sequences capable of non-B DNA formation) in mediating the recombination events that underlie gross genomic rearrangements. Finally, we explore how studies both at the level of the gene (using the neurofibromatosis type-1 gene as an example) and the whole genome (using data derived from cancer genome sequencing studies) are shaping our understanding of the impact of genomic rearrangements as a cause of human genetic disease.

The AZFc region of the Y chromosome: at the crossroads between genetic diversity and male infertility

BACKGROUND

The three azoospermia factor (AZF) regions of the Y chromosome represent genomic niches for spermatogenesis genes. Yet, the most distal region, AZFc, is a major generator of large-scale variation in the human genome. Determining to what extent this variability affects spermatogenesis is a highly contentious topic in human reproduction.

METHODS

In this review, an extensive characterization of the molecular mechanisms responsible for AZFc genotypical variation is undertaken. Such data are complemented with the assessment of the clinical consequences for male fertility imputable to the different AZFc variants. For this, a critical re-evaluation of 23 association studies was performed in order to extract unifying conclusions by curtailing methodological heterogeneities.

RESULTS

Intrachromosomal homologous recombination mechanisms, either crossover or non-crossover based, are the main drivers for AZFc genetic diversity. In particular, rearrangements affecting gene dosage are the most likely to introduce phenotypical disruptions in the spermatogenic profile. In the specific cases of partial AZFc deletions, both the actual existence and the severity of the spermatogenic defect are dependent on the evolutionary background of the Y chromosome.

CONCLUSIONS

AZFc is one of the most genetically dynamic regions in the human genome. This property may serve as counter against the genetic degeneracy associated with the lack of a meiotic partner. However, such strategy comes at a price: some rearrangements represent a risk factor or a de-facto causative agent of spermatogenic disruption. Interestingly, this precarious balance is modulated, among other yet unknown factors, by the evolutionary history of the Y chromosome.

Closely spaced multiple mutations as potential signatures of transient hypermutability in human genes

Data from diverse organisms suggests that transient hypermutability is a general mutational mechanism with the potential to generate multiple synchronous mutations, a phenomenon probably best exemplified by closely spaced multiple mutations (CSMMs). Here we have attempted to extend the concept of transient hypermutability from somatic cells to the germline, using human inherited disease-causing multiple mutations as a model system. Employing stringent criteria for data inclusion, we have retrospectively identified numerous potential examples of pathogenic CSMMs that exhibit marked similarities to the CSMMs reported in other systems. These examples include (1) eight multiple mutations, each comprising three or more components within a sequence tract of <100 bp; (2) three possible instances of "mutation showers"; and (3) numerous highly informative "homocoordinate" mutations. Using the proportion of CpG substitution as a crude indicator of the relative likelihood of transient hypermutability, we present evidence to suggest that CSMMs comprising at least one pair of mutations separated by < or =100 bp may constitute signatures of transient hypermutability in human genes. Although this analysis extends the generality of the concept of transient hypermutability and provides new insights into what may be considered a novel mechanism of mutagenesis underlying human inherited disease, it has raised serious concerns regarding current practices in mutation screening.