Human genomic deletions mediated by recombination between Alu elements

Breakpoint characterization of large deletions in EXT1 or EXT2 in 10 multiple osteochondromas families

Background

Osteochondromas (cartilage-capped bone tumors) are by far the most commonly treated of all primary benign bone tumors (50%). In 15% of cases, these tumors occur in the context of a hereditary syndrome called multiple osteochondromas (MO), an autosomal dominant skeletal disorder characterized by the formation of multiple cartilage-capped bone tumors at children's metaphyses. MO is caused by various mutations in EXT1 or EXT2, whereby large genomic deletions (single-or multi-exonic) are responsible for up to 8% of MO-cases.

Methods

Here we report on the first molecular characterization of ten large EXT1- and EXT2-deletions in MO-patients. Deletions were initially indentified using MLPA or FISH analysis and were subsequently characterized using an MO-specific tiling path array, allele-specific PCR-amplification and sequencing analysis.

Results

Within the set of ten large deletions, the deleted regions ranged from 2.7 to 260 kb. One EXT2 exon 8 deletion was found to be recurrent. All breakpoints were located outside the coding exons of EXT1 and EXT2. Non-allelic homologous recombination (NAHR) mediated by Alu-sequences, microhomology mediated replication dependent recombination (MMRDR) and non-homologous end-joining (NHEJ) were hypothesized as the causal mechanisms in different deletions.

Conclusions

Molecular characterization of EXT1- and EXT2-deletion breakpoints in MO-patients indicates that NAHR between Alu-sequences as well as NHEJ are causal and that the majority of these deletions are nonrecurring. These observations emphasize once more the huge genetic variability which is characteristic for MO. To our knowledge, this is the first study characterizing large genomic deletions in EXT1 and EXT2.

Mobile DNA and the TE-Thrust hypothesis: supporting evidence from the primates

Transposable elements (TEs) are increasingly being recognized as powerful facilitators of evolution. We propose the TE-Thrust hypothesis to encompass TE-facilitated processes by which genomes self-engineer coding, regulatory, karyotypic or other genetic changes. Although TEs are occasionally harmful to some individuals, genomic dynamism caused by TEs can be very beneficial to lineages. This can result in differential survival and differential fecundity of lineages. Lineages with an abundant and suitable repertoire of TEs have enhanced evolutionary potential and, if all else is equal, tend to be fecund, resulting in species-rich adaptive radiations, and/or they tend to undergo major evolutionary transitions. Many other mechanisms of genomic change are also important in evolution, and whether the evolutionary potential of TE-Thrust is realized is heavily dependent on environmental and ecological factors. The large contribution of TEs to evolutionary innovation is particularly well documented in the primate lineage. In this paper, we review numerous cases of beneficial TE-caused modifications to the genomes of higher primates, which strongly support our TE-Thrust hypothesis.

Genomic structure and expression analyses of the PYGM gene in the thoroughbred horse

Muscle glycogen Phosphorylase (PYGM) has been shown to catalyze the degradation of glycogen to glucose-1-phosphate. The PYGM gene can contribute to providing energy to the body by disassembling the glycogen in muscle. Here, we analyzed the genomic structure and expression of the PYGM gene in the thoroughbred horse. The PYGM gene, containing several transposable elements (MIRs, LINEs, and MERs), was highly conserved in mammalian genomes. In order to understand the expression of the horse PYGM gene, we performed quantitative RT-PCR using 11 thoroughbred horse tissue samples. The horse PYGM gene was broadly expressed in all tissues tested. In particular, the highest expression of the horse PYGM gene was observed in skeletal muscle tissue relative to the other tissues. Interestingly, the horse PYGM gene contains fewer mobile elements than its human ortholog, resulting in an increase in the structural stability of the PYGM gene sequence. This study provides insights into the genomic structure of the horse PYGM gene that may be useful in future studies of its association with exercise capability.

Alu distribution and mutation types of cancer genes

Background

Alu elements are the most abundant retrotransposable elements comprising ~11% of the human genome. Many studies have highlighted the role that Alu elements have in genetic instability and how their contribution to the assortment of mutagenic events can lead to cancer. As of yet, little has been done to quantitatively assess the association between Alu distribution and genes that are causally implicated in oncogenesis.

Results

We have investigated the effect of various Alu densities on the mutation type based classifications of cancer genes. In order to establish the direct relationship between Alus and the cancer genes of interest, genome wide Alu-related densities were measured using genes rather than the sliding windows of fixed length as the units. Several novel genomic features, such as the density of the adjacent Alu pairs and the number of Alu-Exon-Alu triplets, were developed in order to extend the investigation via the multivariate statistical analysis toward more advanced biological insight. In addition, we characterized the genome-wide intron Alu distribution with a mixture model that distinguished genes containing Alu elements from those with no Alus, and evaluated the gene-level effect of the 5'-TTAAAA motif associated with Alu insertion sites using a two-step regression analysis method.

Conclusions

The study resulted in several novel findings worthy of further investigation. They include: (1) Recessive cancer genes (tumor suppressor genes) are enriched with Alu elements (p < 0.01) compared to dominant cancer genes (oncogenes) and the entire set of genes in the human genome; (2) Alu-related genomic features can be used to cluster cancer genes into biological meaningful groups; (3) The retention of exon Alus has been restricted in the human genome development, and an upper limit to the chromosome-level exon Alu densities is suggested by the distribution profile; (4) For the genes with at least one intron Alu repeat in individual chromosomes, the intron Alu densities can be well fitted by a Gamma distribution; (5) The effect of the 5'-TTAAAA motif on Alu densities varies across different chromosomes.

Isolation of MECP2-null Rett Syndrome patient hiPS cells and isogenic controls through X-chromosome inactivation

Rett syndrome (RTT) is a neurodevelopmental autism spectrum disorder that affects girls due primarily to mutations in the gene encoding methyl-CpG binding protein 2 (MECP2). The majority of RTT patients carry missense and nonsense mutations leading to a hypomorphic MECP2, while null mutations leading to the complete absence of a functional protein are rare. MECP2 is an X-linked gene subject to random X-chromosome inactivation resulting in mosaic expression of mutant MECP2. The lack of human brain tissue motivates the need for alternative human cellular models to study RTT. Here we report the characterization of a MECP2 mutation in a classic female RTT patient involving rearrangements that remove exons 3 and 4 creating a functionally null mutation. To generate human neuron models of RTT, we isolated human induced pluripotent stem (hiPS) cells from RTT patient fibroblasts. RTT-hiPS cells retained the MECP2 mutation, are pluripotent and fully reprogrammed, and retained an inactive X-chromosome in a nonrandom pattern. Taking advantage of the latter characteristic, we obtained a pair of isogenic wild-type and mutant MECP2 expressing RTT-hiPS cell lines that retained this MECP2 expression pattern upon differentiation into neurons. Phenotypic analysis of mutant RTT-hiPS cell-derived neurons demonstrated a reduction in soma size compared with the isogenic control RTT-hiPS cell-derived neurons from the same RTT patient. Analysis of isogenic control and mutant hiPS cell-derived neurons represents a promising source for understanding the pathogenesis of RTT and the role of MECP2 in human neurons.

Hereditary hemorrhagic telangiectasia: two distinct ENG deletions in one family

Hereditary hemorrhagic telangiectasia (HHT) is an autosomal dominant disorder characterized by aberrant vascular development. Mutations in endoglin (ENG) or activin A receptor type II-like 1 (ACVRL1) account for around 90% of HHT patients, 10% of those are large deletions or duplications. We report here the first observation of two distinct, large ENG deletions segregating in one pedigree. An ENG exon 4-7 deletion was observed in a patient with HHT. This deletion was identified in several affected family members. However, some affected family members had an ENG exon 3 deletion instead. These deletions were detected by multiplex ligation-dependent probe amplification and confirmed by mRNA sequencing and an oligo-CGH array. Linkage analysis revealed that one individual with the exon 3 deletion inherited the same chromosome from his mother who has the exon 4-7 deletion. This finding has important clinical implications because it shows that targeted family-specific mutation analysis for exon deletions could have led to the misdiagnosis of some affected family members.

A human genome structural variation sequencing resource reveals insights into mutational mechanisms

Understanding the prevailing mutational mechanisms responsible for human genome structural variation requires uniformity in the discovery of allelic variants and precision in terms of breakpoint delineation. We develop a resource based on capillary end sequencing of 13.8 million fosmid clones from 17 human genomes and characterize the complete sequence of 1054 large structural variants corresponding to 589 deletions, 384 insertions, and 81 inversions. We analyze the 2081 breakpoint junctions and infer potential mechanism of origin. Three mechanisms account for the bulk of germline structural variation: microhomology-mediated processes involving short (2-20 bp) stretches of sequence (28%), nonallelic homologous recombination (22%), and L1 retrotransposition (19%). The high quality and long-range continuity of the sequence reveals more complex mutational mechanisms, including repeat-mediated inversions and gene conversion, that are most often missed by other methods, such as comparative genomic hybridization, single nucleotide polymorphism microarrays, and next-generation sequencing.

RISCI--Repeat Induced Sequence Changes Identifier: a comprehensive, comparative genomics-based, in silico subtractive hybridization pipeline to identify repeat induced sequence changes in closely related genomes

Background -

The availability of multiple whole genome sequences has facilitated in silico identification of fixed and polymorphic transposable elements (TE). Whereas polymorphic loci serve as makers for phylogenetic and forensic analysis, fixed species-specific transposon insertions, when compared to orthologous loci in other closely related species, may give insights into their evolutionary significance. Besides, TE insertions are not isolated events and are frequently associated with subtle sequence changes concurrent with insertion or post insertion. These include duplication of target site, 3' and 5' flank transduction, deletion of the target locus, 5' truncation or partial deletion and inversion of the transposon, and post insertion changes like inter or intra element recombination, disruption etc. Although such changes have been studied independently, no automated platform to identify differential transposon insertions and the associated array of sequence changes in genomes of the same or closely related species is available till date. To this end, we have designed RISCI - 'Repeat Induced Sequence Changes Identifier' - a comprehensive, comparative genomics-based, in silico subtractive hybridization pipeline to identify differential transposon insertions and associated sequence changes using specific alignment signatures, which may then be examined for their downstream effects.

Results -

We showcase the utility of RISCI by comparing full length and truncated L1HS and AluYa5 retrotransposons in the reference human genome with the chimpanzee genome and the alternate human assemblies (Celera and HuRef). Comparison of the reference human genome with alternate human assemblies using RISCI predicts 14 novel polymorphisms in full length L1HS, 24 in truncated L1HS and 140 novel polymorphisms in AluYa5 insertions, besides several insertion and post insertion changes. We present comparison with two previous studies to show that RISCI predictions are broadly in agreement with earlier reports. We also demonstrate its versatility by comparing various strains of Mycobacterium tuberculosis for IS 6100 insertion polymorphism.

Conclusions -

RISCI combines comparative genomics with subtractive hybridization, inferring changes only when exclusive to one of the two genomes being compared. The pipeline is generic and may be applied to most transposons and to any two or more genomes sharing high sequence similarity. Such comparisons, when performed on a larger scale, may pull out a few critical events, which may have seeded the divergence between the two species under comparison.

A common molecular mechanism underlies two phenotypically distinct 17p13.1 microdeletion syndromes

DNA copy-number variations (CNVs) underlie many neuropsychiatric conditions, but they have been less studied in cancer. We report the association of a 17p13.1 CNV, childhood-onset developmental delay (DD), and cancer. Through a screen of over 4000 patients with diverse diagnoses, we identified eight probands harboring microdeletions at TP53 (17p13.1). We used a purpose-built high-resolution array with 93.75% breakpoint accuracy to fine map these microdeletions. Four patients were found to have a common phenotype including DD, hypotonia, and hand and foot abnormalities, constituting a unique syndrome. Notably, these patients were not affected with cancer. Moreover, none of the TP53-deletion patients affected with cancer (n = 4) had neurocognitive impairments. DD patients have larger deletions, which encompass but do not disrupt TP53, whereas cancer-affected patients harbor CNVs with at least one breakpoint within TP53. Most 17p13.1 deletions arise by Alu-mediated nonallelic homologous recombination. Furthermore, we identify a critical genomic region associated with DD and containing six underexpressed genes. We conclude that, although they overlap, 17p13.1 CNVs are associated with distinct phenotypes depending on the position of the breakpoint with respect to TP53. Further, detailed characterization of breakpoints revealed a common formation signature. Future studies should consider whether other loci in the genome also give rise to phenotypically distinct disorders by means of a common mechanism, resulting in a similar formation signature.

GJB1/Connexin 32 whole gene deletions in patients with X-linked Charcot-Marie-Tooth disease

The X-linked form of Charcot-Marie-Tooth disease (CMTX) is the second most common form of this genetically heterogeneous inherited peripheral neuropathy. CMT1X is caused by mutations in the GJB1 gene. Most of the mutations causative for CMT1X are missense mutations. In addition, a few disease causative nonsense mutations and frameshift deletions that lead to truncated forms of the protein have also been reported to be associated with CMT1X. Previously, there have been reports of patients with deletions of the coding sequence of GJB1; however, the size and breakpoints of these deletions were not assessed. Here, we report five patients with deletions that range in size from 12.2 to 48.3 kb and that completely eliminate the entire coding sequence of the GJB1 gene, resulting in a null allele for this locus. Analyses of the breakpoints of these deletions showed that they are nonrecurrent and that they can be generated by different mechanisms. In addition to PMP22, GJB1 is the second CMT gene for which both point mutations and genomic rearrangements can cause a neuropathy phenotype, stressing the importance of CMT as a genomic disorder.

Are we Genomic Mosaics? Variations of the Genome of Somatic Cells can Contribute to Diversify our Phenotypes

Theoretical and experimental evidences support the hypothesis that the genomes and the epigenomes may be different in the somatic cells of complex organisms. In the genome, the differences range from single base substitutions to chromosome number; in the epigenome, they entail multiple postsynthetic modifications of the chromatin. Somatic genome variations (SGV) may accumulate during development in response both to genetic programs, which may differ from tissue to tissue, and to environmental stimuli, which are often undetected and generally irreproducible. SGV may jeopardize physiological cellular functions, but also create novel coding and regulatory sequences, to be exposed to intraorganismal Darwinian selection. Genomes acknowledged as comparatively poor in genes, such as humans', could thus increase their pristine informational endowment. A better understanding of SGV will contribute to basic issues such as the "nature vs nurture" dualism and the inheritance of acquired characters. On the applied side, they may explain the low yield of cloning via somatic cell nuclear transfer, provide clues to some of the problems associated with transdifferentiation, and interfere with individual DNA analysis. SGV may be unique in the different cells types and in the different developmental stages, and thus explain the several hundred gaps persisting in the human genomes "completed" so far. They may compound the variations associated to our epigenomes and make of each of us an "(epi)genomic" mosaic. An ensuing paradigm is the possibility that a single genome (the ephemeral one assembled at fertilization) has the capacity to generate several different brains in response to different environments.

What drives recombination hotspots to repeat DNA in humans?

Recombination between homologous, but non-allelic, stretches of DNA such as gene families, segmental duplications and repeat elements is an important source of mutation. In humans, recent studies have identified short DNA motifs that both determine the location of 40 per cent of meiotic cross-over hotspots and are significantly enriched at the breakpoints of recurrent non-allelic homologous recombination (NAHR) syndromes. Unexpectedly, the most highly penetrant form of the motif occurs on the background of an inactive repeat element family (THE1 elements) and the motif also has strong recombinogenic activity on currently active element families including Alu and LINE2 elements. Analysis of genetic variation among members of these repeat families indicates an important role for NAHR in their evolution. Given the potential for double-strand breaks within repeat DNA to cause pathological rearrangement, the association between repeats and hotspots is surprising. Here we consider possible explanations for why selection acting against NAHR has not eliminated hotspots from repeat DNA including mechanistic constraints, possible benefits to repeat DNA from recruiting hotspots and rapid evolution of the recombination machinery. I suggest that rapid evolution of hotspot motifs may, surprisingly, tend to favour sequences present in repeat DNA and outline the data required to differentiate between hypotheses.

Breakpoint determination of 15 large deletions in Peutz-Jeghers subjects

The Peutz-Jeghers Syndrome (PJS) is an autosomal dominant polyposis disorder with increased risk of multiple cancers. STK11/LKB1 (hereafter named STK11) germline mutations account for the large majority of PJS cases whereas large deletions account for about 30% of the cases. We report here the first thorough molecular characterization of 15 large deletions identified in a cohort of 51 clinically well-characterized PJS patients. The deletions were identified by MLPA analysis and characterized by custom CGH-array and quantitative PCR to define their boundaries. The deletions, ranging from 2.9 to 180 kb, removed one or more loci contiguous to the STK11 gene in six patients, while partial STK11 gene deletions were present in the remaining nine cases. By means of DNA sequencing, we were able to precisely characterize the breakpoints in each case. Of the 30 breakpoints, 16 were located in Alu elements, revealing non-allelic homologous recombination (NAHR) as the putative mechanism for the deletions of the STK11 gene, which lays in a region with high Alu density. In the remaining cases, other mechanisms could be hypothesized, such as microhomology-mediated end-joining (MMEJ) or non-homologous end-joining (NHEJ). In conclusion we here demonstrated the non-random occurrence of large deletions associated with PJS. All our patients had a classical PJS phenotype, which shows that haploinsufficiency for SBNO2, C19orf26, ATP5D, MIDN, C19orf23, CIRBP, C19orf24,and EFNA2, does not apparently affect their clinical phenotype.

Mechanisms for nonrecurrent genomic rearrangements associated with CMT1A or HNPP: rare CNVs as a cause for missing heritability

Genomic rearrangements involving the peripheral myelin protein gene (PMP22) in human chromosome 17p12 are associated with neuropathy: duplications cause Charcot-Marie-Tooth disease type 1A (CMT1A), whereas deletions lead to hereditary neuropathy with liability to pressure palsies (HNPP). Our previous studies showed that >99% of these rearrangements are recurrent and mediated by nonallelic homologous recombination (NAHR). Rare copy number variations (CNVs) generated by nonrecurrent rearrangements also exist in 17p12, but their underlying mechanisms are not well understood. We investigated 21 subjects with rare CNVs associated with CMT1A or HNPP by oligonucleotide-based comparative genomic hybridization microarrays and breakpoint sequence analyses, and we identified 17 unique CNVs, including two genomic deletions, ten genomic duplications, two complex rearrangements, and three small exonic deletions. Each of these CNVs includes either the entire PMP22 gene, or exon(s) only, or ultraconserved potential regulatory sequences upstream of PMP22, further supporting the contention that PMP22 is the critical gene mediating the neuropathy phenotypes associated with 17p12 rearrangements. Breakpoint sequence analysis reveals that, different from the predominant NAHR mechanism in recurrent rearrangement, various molecular mechanisms, including nonhomologous end joining, Alu-Alu-mediated recombination, and replication-based mechanisms (e.g., FoSTeS and/or MMBIR), can generate nonrecurrent 17p12 rearrangements associated with neuropathy. We document a multitude of ways in which gene function can be altered by CNVs. Given the characteristics, including small size, structural complexity, and location outside of coding regions, of selected rare CNVs, their identification remains a challenge for genome analysis. Rare CNVs may potentially represent an important portion of "missing heritability" for human diseases.

Four different ways of alternative transcripts generation mechanism in ADRA1A gene

The ADRA1A (Alpha-1-adrenergic receptor) gene is one of the members of G protein-coupled receptor superfamily. Alternative splicing of this gene was known to generate four transcript variants which code four isoforms with various C-terminal regions. In this study, we conducted expression analysis and evolutionary characterization of alternative transcripts of the ADRA1A gene. In total, 10 alternative transcripts were identified using experimental approaches and in silico analysis. Among them, 6 alternative transcripts (T1, T2, T3, T4, T4-1, and T8) were validated by RT-PCR approaches and sequencing procedures. From the alternative splicing analysis, it could be assumed that there were three different alternative transcripts generation mechanisms and unknown mechanism. First one is the integration event of three different TEs (AluSc, L1MC5, and MIR3) as seen on the last exons of T3, T4, T4-1, T5, T6, and T7 transcripts. The second mechanism is a differential promoter usage on T8. The third one is a substitution event of the 3' splicing site during the primate evolution on T3. The last one is an unknown mechanism which was identified in the T4-1 transcript. Therefore, alternative transcripts of human ADRA1A gene occurred by four different ways, such as integration of TEs, differential promoter usage, substitution of splicing sites, and unknown mechanism.

[Functional analysis of retroviral endogenous inserts in the human genome evolution]

Retroelements, mobile elements produced in DNA by reverse transcription, comprise about 40% of the human genome. A small part of these elements appeared in the genome quite recently after the divergence of humans and chimpanzees had occurred. Evolutionarily young retroelements are represented by the members of four groups, SVA, Alu, L1, and the endogenous HERV-K (HML-2) virus. These retroelements could play a functional role in the course of the molecular evolution of human DNA. We comprehensively studied the contribution of human-specific endogenous viruses (hsERV) to the structural modifications and regulation of the human genome. We found that hsERV presented in 134 copies occupied about 330 000 bp of human DNA. They added to genomic sequences the copies of 50 functional retroviral genes as well as 134 potential promoters and enhancers, 50% of which are located in the regions adjacent to known genes, and 22% in gene introns. At least 67% of these elements are human-specific promoters in vivo. hsERV viruses regulate the activity of known protein-encoding genes by means of RNA interference, function as enhancers, and provide new polyadenylation signals for mRNA.

Recombination-associated sequence homogenization of neighboring Alu elements: signature of nonallelic gene conversion

Recently, researchers have begun to recognize that, in order to establish neutral models for disease association and evolutionary genomics studies, it is crucial to have a clear understanding of the genomic impact of nonallelic gene conversion. Drawing on previous successes in characterizing this phenomenon over protein-coding gene families, we undertook a computational analysis of neighboring Alu sequences in the genome scale. For this purpose, we developed adjusted comutation rate (aCMR), a novel statistical method measuring the excess number of identical point mutations shared by adjacent Alu sequences, vis-à-vis random pairs. Using aCMR, we uncovered a remarkable genome-wide sequence homogenization of neighboring Alus, with the strongest signal observed in the pseudoautosomal regions of the X and Y chromosomes. The magnitude of sequence homogenization between Alu pairs is greater with shorter interlocus distance, higher sequence identity, and parallel orientation. Moreover, shared substitutions show a strong directionality toward GC nucleotides, with multiple substitutions tending to cluster within the Alu sequence. Taken together, these observed recombination-associated sequence homogenization patterns are best explained by frequent ubiquitous gene conversion events between neighboring Alus. We believe that these observations help to illuminate the nature and impact of the enigmatic phenomenon of gene conversion.

Identification of a Tc1-like transposon integration site in the genome of the flounder (Platichthys flesus): a novel use of an inverse PCR method

The inverse PCR method has been developed and applied employed for the identification of the integration sites of the Tc1-like transposons in the genome of the flounder, Platichthys flesus. One Tc1-like insertion instance was recognized and characterized, demonstrating an efficiency of the method for determining of transposon integration sites. The similarity of the sequence flanking transposon (SFT) to reverse transcriptase sequences (RVT) was demonstrated. It is likely that the insertion took place within currently degenerated LINE (long interspersed nuclear elements) retrotransposon.

Homozygous deletion of the STK11/LKB1 locus and the generation of novel fusion transcripts in cervical cancer cells

The STK11/LKB1 gene encodes a ubiquitously expressed serine/threonine kinase that is mutated in multiple sporadic cancers including non-small cell lung carcinomas, pancreatic cancers, and melanomas. LKB1 plays a role in multiple cellular functions including cell growth, cell cycle progression, metabolism, cell polarity, and migration. To date, only a limited number of studies have assessed the status of LKB1 in cervical cancers. Herein, we investigate DNA methylation, DNA mutation, and transcription at the LKB1 locus in cervical cancer cell lines. We identified homozygous deletions of 25-85kb in the HeLa and SiHa cell lines. Deletion breakpoint analysis in HeLa cells revealed that the deletion resulted from an Alu-recombination-mediated deletion (ARMD) and generated a novel LKB1 fusion transcript driven by an uncharacterized CpG island promoter located approximately 11kb upstream of LKB1. Although the homozygous deletion in SiHa cells removes the entire LKB1 gene and portions of the neighboring genes SBNO2 and c19orf26, this deletion also generates a fusion transcript driven by the c19orf26 promoter and composed of both c19orf26 and SBNO2 sequences. Further analyses of public gene expression and mutation databases suggest that LKB1 and its neighboring genes are frequently dysregulated in primary cervical cancers. Thus, homozygous deletions affecting LKB1 in cervical cancers may generate multiple fusion transcripts involving LKB1, SBNO2, and c19orf26.

A mobile threat to genome stability: The impact of non-LTR retrotransposons upon the human genome

It is now commonly agreed that the human genome is not the stable entity originally presumed. Deletions, duplications, inversions, and insertions are common, and contribute significantly to genomic structural variations (SVs). Their collective impact generates much of the inter-individual genomic diversity observed among humans. Not only do these variations change the structure of the genome; they may also have functional implications, e.g. altered gene expression. Some SVs have been identified as the cause of genetic disorders, including cancer predisposition. Cancer cells are notorious for their genomic instability, and often show genomic rearrangements at the microscopic and submicroscopic level to which transposable elements (TEs) contribute. Here, we review the role of TEs in genome instability, with particular focus on non-LTR retrotransposons. Currently, three non-LTR retrotransposon families - long interspersed element 1 (L1), SVA (short interspersed element (SINE-R), variable number of tandem repeats (VNTR), and Alu), and Alu (a SINE) elements - mobilize in the human genome, and cause genomic instability through both insertion- and post-insertion-based mutagenesis. Due to the abundance and high sequence identity of TEs, they frequently mislead the homologous recombination repair pathway into non-allelic homologous recombination, causing deletions, duplications, and inversions. While less comprehensively studied, non-LTR retrotransposon insertions and TE-mediated rearrangements are probably more common in cancer cells than in healthy tissue. This may be at least partially attributed to the commonly seen global hypomethylation as well as general epigenetic dysfunction of cancer cells. Where possible, we provide examples that impact cancer predisposition and/or development.

Identification of uncommon recurrent Potocki-Lupski syndrome-associated duplications and the distribution of rearrangement types and mechanisms in PTLS

Nonallelic homologous recombination (NAHR) can mediate recurrent rearrangements in the human genome and cause genomic disorders. Smith-Magenis syndrome (SMS) and Potocki-Lupski syndrome (PTLS) are genomic disorders associated with a 3.7 Mb deletion and its reciprocal duplication in 17p11.2, respectively. In addition to these common recurrent rearrangements, an uncommon recurrent 5 Mb SMS-associated deletion has been identified. However, its reciprocal duplication predicted by the NAHR mechanism had not been identified. Here we report the molecular assays on 74 subjects with PTLS-associated duplications, 35 of whom are newly investigated. By both oligonucleotide-based comparative genomic hybridization and recombination hot spot analyses, we identified two cases of the predicted 5 Mb uncommon recurrent PTLS-associated duplication. Interestingly, the crossovers occur in proximity to a recently delineated allelic homologous recombination (AHR) hot spot-associated sequence motif, further documenting the common hot spot features shared between NAHR and AHR. An additional eight subjects with nonrecurrent PTLS duplications were identified. The smallest region of overlap (SRO) for all of the 74 PTLS duplications examined is narrowed to a 125 kb interval containing only RAI1, a gene recently further implicated in autism. Sequence complexities consistent with DNA replication-based mechanisms were identified in four of eight (50%) newly identified nonrecurrent PTLS duplications. Our findings of the uncommon recurrent PTLS-associated duplication at a relative prevalence reflecting the de novo mutation rate and the distribution of 17p11.2 duplication types in PTLS reveal insights into both the contributions of new mutations and the different underlying mechanisms that generate genomic rearrangements causing genomic disorders.

The (r)evolution of SINE versus LINE distributions in primate genomes: sex chromosomes are important

The densities of transposable elements (TEs) in the human genome display substantial variation both within individual chromosomes and among chromosome types (autosomes and the two sex chromosomes). Finding an explanation for this variability has been challenging, especially in light of genome landscapes unique to the sex chromosomes. Here, using a multiple regression framework, we investigate primate Alu and L1 densities shaped by regional genome features and location on a particular chromosome type. As a result of our analysis, first, we build statistical models explaining up to 79% and 44% of variation in Alu and L1 element density, respectively. Second, we analyze sex chromosome versus autosome TE densities corrected for regional genomic effects. We discover that sex-chromosome bias in Alu and L1 distributions not only persists after accounting for these effects, but even presents differences in patterns, confirming preferential Alu integration in the male germline, yet likely integration of L1s in both male and female germlines or in early embryogenesis. Additionally, our models reveal that local base composition (measured by GC content and density of L1 target sites) and natural selection (inferred via density of most conserved elements) are significant to predicting densities of L1s. Interestingly, measurements of local double-stranded breaks (a 13-mer associated with genome instability) strongly correlate with densities of Alu elements; little evidence was found for the role of recombination-driven deletion in driving TE distributions over evolutionary time. Thus, Alu and L1 densities have been influenced by the combination of distinct local genome landscapes and the unique evolutionary dynamics of sex chromosomes.

Unique functions of repetitive transcriptomes

Repetitive sequences occupy a huge fraction of essentially every eukaryotic genome. Repetitive sequences cover more than 50% of mammalian genomic DNAs, whereas gene exons and protein-coding sequences occupy only ~3% and 1%, respectively. Numerous genomic repeats include genes themselves. They generally encode "selfish" proteins necessary for the proliferation of transposable elements (TEs) in the host genome. The major part of evolutionary "older" TEs accumulated mutations over time and fails to encode functional proteins. However, repeats have important functions also on the RNA level. Repetitive transcripts may serve as multifunctional RNAs by participating in the antisense regulation of gene activity and by competing with the host-encoded transcripts for cellular factors. In addition, genomic repeats include regulatory sequences like promoters, enhancers, splice sites, polyadenylation signals, and insulators, which actively reshape cellular transcriptomes. TE expression is tightly controlled by the host cells, and some mechanisms of this regulation were recently decoded. Finally, capacity of TEs to proliferate in the host genome led to the development of multiple biotechnological applications.

Laboratory methods for the analysis of primate mobile elements

Mobile elements represent a unique and powerful set of tools for understanding the variation in a genome. Methods exist not only to utilize the polymorphisms among and within taxa to various ends but also to investigate the mechanism through which mobilization occurs. The number of methods to accomplish these ends is ever growing. Here, we present several protocols designed to assay mobile element-based variation within and among individual genomes.

Computational methods for the analysis of primate mobile elements

Transposable elements (TE), defined as discrete pieces of DNA that can move from one site to another site in genomes, represent significant components of eukaryotic genomes, including primates. Comparative genome-wide analyses have revealed the considerable structural and functional impact of TE families on primate genomes. Insights into these questions have come in part from the development of computational methods that allow detailed and reliable identification, annotation, and evolutionary analyses of the many TE families that populate primate genomes. Here, we present an overview of these computational methods and describe efficient data mining strategies for providing a comprehensive picture of TE biology in newly available genome sequences.

Mechanisms of copy number variation and hybrid gene formation in the KIR immune gene complex

The fine-scale structure of the majority of copy number variation (CNV) regions remains unknown. The killer immunoglobulin receptor (KIR) gene complex exhibits significant CNV. The evolutionary plasticity of the KIRs and their broad biomedical relevance makes it important to understand how these immune receptors evolve. In this paper, we describe haplotype re-arrangement creating novel loci at the KIR complex. We completely sequenced, after fosmid cloning, two rare contracted haplotypes. Evidence of frequent hybrid KIR genes in samples from many populations suggested that re-arrangements may be frequent and selectively advantageous. We propose mechanisms for formation of novel hybrid KIR genes, facilitated by protrusive non-B DNA structures at transposon recombination sites. The heightened propensity to generate novel hybrid KIR receptors may provide a proactive evolutionary measure, to militate against pathogen evasion or subversion. We propose that CNV in KIR is an evolutionary strategy, which KIR typing for disease association must take into account.

Retroelements and their impact on genome evolution and functioning

Retroelements comprise a considerable fraction of eukaryotic genomes. Since their initial discovery by Barbara McClintock in maize DNA, retroelements have been found in genomes of almost all organisms. First considered as a "junk DNA" or genomic parasites, they were shown to influence genome functioning and to promote genetic innovations. For this reason, they were suggested as an important creative force in the genome evolution and adaptation of an organism to altered environmental conditions. In this review, we summarize the up-to-date knowledge of different ways of retroelement involvement in structural and functional evolution of genes and genomes, as well as the mechanisms generated by cells to control their retrotransposition.

Alu repeats increase local recombination rates

Background

Recombination rates vary widely across the human genome, but little of that variation is correlated with known DNA sequence features. The genome contains more than one million Alu mobile element insertions, and these insertions have been implicated in non-homologous recombination, modulation of DNA methylation, and transcriptional regulation. If individual Alu insertions have even modest effects on local recombination rates, they could collectively have a significant impact on the pattern of linkage disequilibrium in the human genome and on the evolution of the Alu family itself.

Results

We carried out sequencing, SNP identification, and SNP genotyping around 19 AluY insertion loci in 347 individuals sampled from diverse populations, then used the SNP genotypes to estimate local recombination rates around the AluY loci. The loci and SNPs were chosen so as to minimize other factors (such as SNP ascertainment bias and SNP density) that could influence recombination rate estimates. We detected a significant increase in recombination rate within ~2 kb of the AluY insertions in our African population sample. To test this observation against a larger set of AluY insertions, we applied our locus- and SNP-selection design and analyses to the HapMap Phase II data. In that data set, we observed a significantly increased recombination rate near AluY insertions in both the CEU and YRI populations.

Conclusion

We show that the presence of a fixed AluY insertion is significantly predictive of an elevated local recombination rate within 2 kb of the insertion, independent of other known predictors. The magnitude of this effect, approximately a 6% increase, is comparable to the effects of some recombinogenic DNA sequence motifs identified via their association with recombination hot spots.

Reading between the LINEs to see into the past

Transposable elements (TEs) are an important source of genome diversity and play a crucial role in genome evolution. A recent study by Zhao et al. describes novel patterns of TE diversification in the genome of the extinct mammoth Mammuthus primigenius. Analysis of Mammuthus has provided a unique genome landscape, a pivotal species for understanding TEs and genome evolution and hints at the diversity we verge on discovering by expanding our taxonomic sampling among genomes. Strategies based on this work might also revolutionize investigations of the interface between TE dynamics and genome diversity.

The impact of retrotransposons on human genome evolution

Their ability to move within genomes gives transposable elements an intrinsic propensity to affect genome evolution. Non-long terminal repeat (LTR) retrotransposons--including LINE-1, Alu and SVA elements--have proliferated over the past 80 million years of primate evolution and now account for approximately one-third of the human genome. In this Review, we focus on this major class of elements and discuss the many ways that they affect the human genome: from generating insertion mutations and genomic instability to altering gene expression and contributing to genetic innovation. Increasingly detailed analyses of human and other primate genomes are revealing the scale and complexity of the past and current contributions of non-LTR retrotransposons to genomic change in the human lineage.

Origins and functional impact of copy number variation in the human genome

Structural variations of DNA greater than 1 kilobase in size account for most bases that vary among human genomes, but are still relatively under-ascertained. Here we use tiling oligonucleotide microarrays, comprising 42 million probes, to generate a comprehensive map of 11,700 copy number variations (CNVs) greater than 443 base pairs, of which most (8,599) have been validated independently. For 4,978 of these CNVs, we generated reference genotypes from 450 individuals of European, African or East Asian ancestry. The predominant mutational mechanisms differ among CNV size classes. Retrotransposition has duplicated and inserted some coding and non-coding DNA segments randomly around the genome. Furthermore, by correlation with known trait-associated single nucleotide polymorphisms (SNPs), we identified 30 loci with CNVs that are candidates for influencing disease susceptibility. Despite this, having assessed the completeness of our map and the patterns of linkage disequilibrium between CNVs and SNPs, we conclude that, for complex traits, the heritability void left by genome-wide association studies will not be accounted for by common CNVs.

Molecular mechanisms underlying large genomic deletions in ornithine transcarbamylase (OTC) gene

Ornithine transcarbamylase deficiency (OTCD) is an X-linked urea cycle error causing hyperammonemia and orotic aciduria. Clinical diagnosis is generally confirmed by mutation detection. However, in approximately 20% of the patients, no mutation is found by conventional mutation-searching strategies, which fail to detect deletions spanning at least a whole exon, large rearrangements, or mutations at non-coding regions. To detect large deletions or duplications, we have applied the multiplex ligation-dependent probe amplification (MLPA) methodology to three OTCD patients (two females and one male). MLPA revealed copy number alterations of OTC exons in all of them. The two females were found to be heterozygous for deletions of either exon 2 or exons 6-9, and the male was confirmed to lack all OTC exons. Females' characterization of the deletion breakpoints by long polymerase chain reaction and sequencing revealed the mutations c.78-3544_217-129del5921 and c.541-600_1005 + 1880del10862 corresponding to exon 2 and exon 6-9 deletions, respectively. Examination of the deletion-flanking regions suggests that exon 2 deletion probably resulted from replication slippage facilitated by a secondary structure formed by two inverted Alu repeats, whereas an Alu-Alu homologous recombination was probably responsible for the exon 6-9 deletion. This work contributes to the identification of novel disease-causing mutations in OTCD and increases the knowledge on possible mutational mechanisms generating deletions in OTC.

Evolution of serum albumin intron-1 is shaped by a 5' truncated non-long terminal repeat retrotransposon in western Palearctic water frogs (Neobatrachia)

A 5' truncated non-LTR CR1-like retrotransposon, named RanaCR1, was identified in the serum albumin intron-1 (SAI-1) of at least seven species of western Palearctic water frogs (WPWF). Based on sequence similarity of the carboxy-terminal region (CTR) of ORF2 and/or the highly conserved 3' untranslated region (3' UTR), RanaCR1-like elements occur also in the genome of Xenopus tropicalis and Rana temporaria. Unlike other CR1 elements, RanaCR1 contains a CA microsatellite in its 3' UTR. The low nucleotide diversity of the 3' UTR compared to the CTR and to SAI-1 suggests that this region still plays a role in WPWF, either as a structure-stabilizing element, or within a species-specific transcriptional network. Length variation of water frog SAI-1 sequences is caused by deletions that extend in some cases beyond the 5' or 3' ends of RanaCR1, probably a result of selection for structural and functional stability of the primary transcript. The impact of RanaCR1 on SAI-1 evolution is also indicated by the significant negative correlation between the length of both SAI-1 and RanaCR1 and the percentage GC content of RanaCR1. Both SAI-1 and RanaCR1 sequences support the sister group relationship of R. perezi and R. saharica, which are placed in the phylogenetic tree at a basal position, the sister clade to other water frog taxa. It also supports the monophyly of the R. lessonae group; of Anatolian water frogs (R. cf. bedriagae), which are not conspecific with R. bedriagae, and of the European ridibunda group. Within the ridibunda clade, Greek frogs are clearly separated, supporting the hypothesis that Balkan water frogs represent a distinct species. Frogs from Atyrau (Kazakhstan), the type locality of R. ridibunda, were heterozygous for a ridibunda and a cf. bedriagae specific allele.

Understanding what determines the frequency and pattern of human germline mutations

Surprising findings about human germline mutation have come from applying new technologies to detect rare mutations in germline DNA, from analysing DNA sequence divergence between humans and closely related species, and from investigating human polymorphic variation. In this Review we discuss how these approaches affect our current understanding of the roles of sex, age, mutation hot spots, germline selection and genomic factors in determining human nucleotide substitution mutation patterns and frequencies. To enhance our understanding of mutation and disease, more extensive molecular data on the human germ line with regard to mutation origin, DNA repair, epigenetic status and the effect of newly arisen mutations on gamete development are needed.

Complex human chromosomal and genomic rearrangements

Copy number variation (CNV) is a major source of genetic variation among humans. In addition to existing as benign polymorphisms, CNVs can also convey clinical phenotypes, including genomic disorders, sporadic diseases and complex human traits. CNV results from genomic rearrangements that can represent simple deletion or duplication of a genomic segment, or be more complex. Complex chromosomal rearrangements (CCRs) have been known for some time but their mechanisms have remained elusive. Recent technology advances and high-resolution human genome analyses have revealed that complex genomic rearrangements can account for a large fraction of non-recurrent rearrangements at a given locus. Various mechanisms, most of which are DNA-replication-based, for example fork stalling and template switching (FoSTeS) and microhomology-mediated break-induced replication (MMBIR), have been proposed for generating such complex genomic rearrangements and are probably responsible for CCR.

Mapping insertions, deletions and SNPs on Venter's chromosomes

BACKGROUND

The very recent availability of fully sequenced individual human genomes is a major revolution in biology which is certainly going to provide new insights into genetic diseases and genomic rearrangements.

RESULTS

We mapped the insertions, deletions and SNPs (single nucleotide polymorphisms) that are present in Craig Venter's genome, more precisely on chromosomes 17 to 22, and compared them with the human reference genome hg17. Our results show that insertions and deletions are almost absent in L1 and generally scarce in L2 isochore families (GC-poor L1+L2 isochores represent slightly over half of the human genome), whereas they increase in GC-rich isochores, largely paralleling the densities of genes, retroviral integrations and Alu sequences. The distributions of insertions/deletions are in striking contrast with those of SNPs which exhibit almost the same density across all isochore families with, however, a trend for lower concentrations in gene-rich regions.

CONCLUSIONS

Our study strongly suggests that the distribution of insertions/deletions is due to the structure of chromatin which is mostly open in gene-rich, GC-rich isochores, and largely closed in gene-poor, GC-poor isochores. The different distributions of insertions/deletions and SNPs are clearly related to the two different responsible mechanisms, namely recombination and point mutations.

The DNA replication FoSTeS/MMBIR mechanism can generate genomic, genic and exonic complex rearrangements in humans

We recently proposed a DNA replication-based mechanism of fork stalling and template switching (FoSTeS) to explain the complex genomic rearrangements associated with a dysmyelinating central nervous system disorder in humans. The FoSTeS mechanism has been further generalized and molecular mechanistic details have been provided in the microhomology-mediated break-induced replication (MMBIR) model that may underlie many structural variations in genomes from all domains of life. Here we provide evidence that human genomic rearrangements ranging in size from several megabases to a few hundred base pairs can be generated by FoSTeS/MMBIR. Furthermore, we show that FoSTeS/MMBIR-mediated rearrangements can occur mitotically and can result in duplication or triplication of individual genes or even rearrangements of single exons. The FoSTeS/MMBIR mechanism can explain both the gene duplication-divergence hypothesis and exon shuffling, suggesting an important role in both genome and single-gene evolution.

Diagnostic strategies for autosomal dominant acute porphyrias: retrospective analysis of 467 unrelated patients referred for mutational analysis of the HMBS, CPOX, or PPOX gene

BACKGROUND

Clinically indistinguishable attacks of acute porphyria occur in acute intermittent porphyria (AIP), hereditary coproporphyria (HCP), and variegate porphyria (VP). There are few evidence-based diagnostic strategies for these disorders.

METHODS

The diagnostic sensitivity of mutation detection was determined by sequencing and gene-dosage analysis to search for mutations in 467 sequentially referred, unrelated patients. The diagnostic accuracy of plasma fluorescence scanning, fecal porphyrin analysis, and porphobilinogen deaminase (PBGD) assay was assessed in mutation-positive patients (AIP, 260 patients; VP, 152 patients; HCP, 31 patients).

RESULTS

Sensitivities (95% CI) for mutation detection were as follows: AIP, 98.1% (95.6%-99.2%); HCP, 96.9% (84.3%-99.5%); VP, 100% (95.7%-100%). We identified 5 large deletions in the HMBS gene (hydroxymethylbilane synthase) and one in the CPOX gene (coproporphyrinogen oxidase). The plasma fluorescence scan was positive more often in VP (99% of patients) than in AIP (68%) or HCP (29%). The wavelength of the fluorescence emission peak and the fecal coproporphyrin isomer ratio had high diagnostic specificity and sensitivity for differentiating between AIP, HCP, and VP. DNA analysis followed by PBGD assay in mutation-negative patients had greater diagnostic accuracy for AIP than either test alone.

CONCLUSIONS

When PBG excretion is increased, 2 investigations (plasma fluorescence scanning, the coproporphyrin isomer ratio) are sufficient, with rare exceptions, to identify the type of acute porphyria. When the results of PBG, 5-aminolevulinate, and porphyrin analyses are within reference intervals and clinical suspicion that a past illness was caused by an acute porphyria remains high, mutation analysis of the HMBS gene followed by PBGD assay is an effective strategy for diagnosis or exclusion of AIP.

Mobile elements create structural variation: analysis of a complete human genome

Structural variants (SVs) are common in the human genome. Because approximately half of the human genome consists of repetitive, transposable DNA sequences, it is plausible that these elements play an important role in generating SVs in humans. Sequencing of the diploid genome of one individual human (HuRef) affords us the opportunity to assess, for the first time, the impact of mobile elements on SVs in an individual in a thorough and unbiased fashion. In this study, we systematically evaluated more than 8000 SVs to identify mobile element-associated SVs as small as 100 bp and specific to the HuRef genome. Combining computational and experimental analyses, we identified and validated 706 mobile element insertion events (including Alu, L1, SVA elements, and nonclassical insertions), which added more than 305 kb of new DNA sequence to the HuRef genome compared with the Human Genome Project (HGP) reference sequence (hg18). We also identified 140 mobile element-associated deletions, which removed approximately 126 kb of sequence from the HuRef genome. Overall, approximately 10% of the HuRef-specific indels larger than 100 bp are caused by mobile element-associated events. More than one-third of the insertion/deletion events occurred in genic regions, and new Alu insertions occurred in exons of three human genes. Based on the number of insertions and the estimated time to the most recent common ancestor of HuRef and the HGP reference genome, we estimated the Alu, L1, and SVA retrotransposition rates to be one in 21 births, 212 births, and 916 births, respectively. This study presents the first comprehensive analysis of mobile element-related structural variants in the complete DNA sequence of an individual and demonstrates that mobile elements play an important role in generating inter-individual structural variation.

Human-specific modulation of transcriptional activity provided by endogenous retroviral insertions

Many phenotypic differences exist between Homo sapiens and its closest relatives, chimpanzees, and these differences can arise as a result of variations in the regulation of certain genes common to these closely related species. Human-specific endogenous retroviruses (HERVs) and their solitary long terminal repeats (LTRs) are probable candidates for such a role due to the presence of regulatory elements, such as enhancers, promoters, splice sites, and polyadenylation signals. In this study we show for the first time that HERVs can participate in the specific antisense regulation of human gene expression owing to their LTR promoter activity. We found that two HERV LTRs situated in the introns of genes SLC4A8 (for sodium bicarbonate cotransporter) and IFT172 (for intraflagellar transport protein 172) in the antisense orientation serve in vivo as promoters for generating RNAs complementary to the exons of enclosing genes. The antisense transcripts formed from LTR promoter were shown to decrease the mRNA level of the corresponding genes. The human-specific regulation of these genes suggests their involvement in the evolutionary process.

Copy number variation in human health, disease, and evolution

Copy number variation (CNV) is a source of genetic diversity in humans. Numerous CNVs are being identified with various genome analysis platforms, including array comparative genomic hybridization (aCGH), single nucleotide polymorphism (SNP) genotyping platforms, and next-generation sequencing. CNV formation occurs by both recombination-based and replication-based mechanisms and de novo locus-specific mutation rates appear much higher for CNVs than for SNPs. By various molecular mechanisms, including gene dosage, gene disruption, gene fusion, position effects, etc., CNVs can cause Mendelian or sporadic traits, or be associated with complex diseases. However, CNV can also represent benign polymorphic variants. CNVs, especially gene duplication and exon shuffling, can be a predominant mechanism driving gene and genome evolution.

Chromosomal inversions between human and chimpanzee lineages caused by retrotransposons

The long interspersed element-1 (LINE-1 or L1) and Alu elements are the most abundant mobile elements comprising 21% and 11% of the human genome, respectively. Since the divergence of human and chimpanzee lineages, these elements have vigorously created chromosomal rearrangements causing genomic difference between humans and chimpanzees by either increasing or decreasing the size of genome. Here, we report an exotic mechanism, retrotransposon recombination-mediated inversion (RRMI), that usually does not alter the amount of genomic material present. Through the comparison of the human and chimpanzee draft genome sequences, we identified 252 inversions whose respective inversion junctions can clearly be characterized. Our results suggest that L1 and Alu elements cause chromosomal inversions by either forming a secondary structure or providing a fragile site for double-strand breaks. The detailed analysis of the inversion breakpoints showed that L1 and Alu elements are responsible for at least 44% of the 252 inversion loci between human and chimpanzee lineages, including 49 RRMI loci. Among them, three RRMI loci inverted exonic regions in known genes, which implicates this mechanism in generating the genomic and phenotypic differences between human and chimpanzee lineages. This study is the first comprehensive analysis of mobile element bases inversion breakpoints between human and chimpanzee lineages, and highlights their role in primate genome evolution.