Non-recurrent 17p11.2 deletions are generated by homologous and non-homologous mechanisms

Several recurrent common chromosomal deletion and duplication breakpoints have been localized to large, highly homologous, low-copy repeats (LCRs). The mechanism responsible for these rearrangements, viz., non-allelic homologous recombination between LCR copies, has been well established. However, fewer studies have examined the mechanisms responsible for non-recurrent rearrangements with non-homologous breakpoint regions. Here, we have analyzed four uncommon deletions of 17p11.2, involving the Smith–Magenis syndrome region. Using somatic cell hybrid lines created from patient lymphoblasts, we have utilized a strategy based on the polymerase chain reaction to refine the deletion breakpoints and to obtain sequence data at the deletion junction. Our analyses have revealed that two of the four deletions are a product of Alu/Alu recombination, whereas the remaining two deletions result from a non-homologous end-joining mechanism. Of the breakpoints studied, three of eight are located in LCRs, and five of eight are within repetitive elements, including Alu and MER5B sequences. These findings suggest that higher-order genomic architecture, such as LCRs, and smaller repetitive sequences, such as Alu elements, can mediate chromosomal deletions via homologous and non-homologous mechanisms. These data further implicate homologous recombination as the predominant mechanism of deletion formation in this genomic interval.     

Genomic Evolution of Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) Isolates Revealed by Deep Sequencing

Most studies on PRRSV evolution have been limited to a particular region of the viral genome. A thorough genome-wide understanding of the impact of different mechanisms on shaping PRRSV genetic diversity is still lacking. To this end, deep sequencing was used to obtain genomic sequences of a diverse set of 16 isolates from a region of Hong Kong with a complex PRRSV epidemiological record. Genome assemblies and phylogenetic typing indicated the co-circulation of strains of both genotypes (type 1and type 2) with varying Nsp2 deletion patterns and distinct evolutionary lineages (“High Fever”-like and local endemic type). Recombination analyses revealed genomic breakpoints in structural and non-structural regions of genomes of both genotypes with evidence of many recombination events originating from common ancestors. Additionally, the high fold of coverage per nucleotide allowed the characterization of minor variants arising from the quasispecies of each strain. Overall, 0.56–2.83% of sites were found to be polymorphic with respect to cognate consensus genomes. The distribution of minor variants across each genome was not uniform indicating the influence of selective forces. Proportion of variants capable of causing an amino acid change in their respective codons ranged between 25–67% with many predicted to be non-deleterious. Low frequency deletion variants were also detected providing one possible mechanism for their sudden emergence as cited in previous reports.     

Gross deletions involving IGHM, BTK, or Artemis: a model for genomic lesions mediated by transposable elements

Most genetic disruptions underlying human disease are microlesions, whereas gross lesions are rare with gross deletions being most frequently found (6%). Similar observations have been made in primary immunodeficiency genes, such as BTK, but for unknown reasons the IGHM and DCLRE1C (Artemis) gene defects frequently represent gross deletions ( approximately 60%). We characterized the gross deletion breakpoints in IGHM-, BTK-, and Artemis-deficient patients. The IGHM deletion breakpoints did not show involvement of recombination signal sequences or immunoglobulin switch regions. Instead, five IGHM, eight BTK, and five unique Artemis breakpoints were located in or near sequences derived from transposable elements (TE). The breakpoints of four out of five disrupted Artemis alleles were located in highly homologous regions, similar to Ig subclass deficiencies and Vh deletion polymorphisms. Nevertheless, these observations suggest a role for TEs in mediating gross deletions. The identified gross deletion breakpoints were mostly located in TE subclasses that were specifically overrepresented in the involved gene as compared to the average in the human genome. This concerned both long (LINE1) and short (Alu, MIR) interspersed elements, as well as LTR retrotransposons (ERV). Furthermore, a high total TE content (>40%) was associated with an increased frequency of gross deletions. Both findings were further investigated and confirmed in a total set of 20 genes disrupted in human disease. Thus, to our knowledge for the first time, we provide evidence that a high TE content, irrespective of the type of element, results in the increased incidence of gross deletions as gene disruption underlying human disease.     

Uncommon Alu-mediated NF1 microdeletion with a breakpoint inside the NF1 gene

Neurofibromatosis type 1 (NF1) microdeletion syndrome is caused by haploinsufficiency of the NF1 gene and of gene(s) located in adjacent flanking regions. Most of the NF1 deletions originate by nonallelic homologous recombination between repeated sequences (REP-P and -M) mapped to 17q11.2, while a few uncommon deletions show unusual breakpoints. We characterized an uncommon 1.5-Mb deletion of an NF1 patient displaying a mild phenotype. We applied high-resolution FISH analysis allowing us to obtain the sequence of the first junction fragment of an uncommon deletion showing the telomeric breakpoint inside the IVS23a of the NF1 gene. Sequence analysis of the centromeric and telomeric boundaries revealed that the breakpoints were present in the AluJb and AluSx regions, respectively, showing 85% homology. The centromeric breakpoint is localized inside a chi-like element; a few copies of this sequence are also located very close to both breakpoints. The in silico analysis of the breakpoint intervals, aimed at identifying consensus sequences of several motifs usually involved in deletions and translocations, suggests that Alu sequences, probably associated with the chi-like element, might be the only recombinogenic motif directly mediating this large deletion.     

Analysis of 22 deletion breakpoints in dystrophin intron 49

Over 60% of Duchenne and Becker muscular dystrophies are caused by deletions spanning tens or hundreds of kilobases in the dystrophin gene. The molecular mechanisms underlying the loss of DNA at this genomic locus are not yet understood. By studying the distribution of deletion breakpoints at the genomic level, we have previously shown that intron 49 exhibits a higher relative density of breakpoints than most dystrophin introns. To determine whether the mechanisms leading to deletions in this intron preferentially involve specific sequence elements, we sublocalized 22 deletion endpoints along its length by a polymerase-chain-reaction-based approach and, in particular, analyzed the nucleotide sequences of five deletion junctions. Deletion breakpoints were homogeneously distributed throughout the intron length, and no extensive homology was observed between the sequences adjacent to each breakpoint. However, a short sequence able to curve the DNA molecule was found at or near three breakpoint junctions.     

<Emphasis Type="Italic">Mariner</Emphasis>-like elements in <Emphasis Type="Italic">Rhynchosciara americana</Emphasis> (Sciaridae) genome: molecular and cytological aspects

Two mariner-like elements, Ramar1 and Ramar2, are described in the genome of Rhynchosciara americana, whose nucleotide consensus sequences were derived from multiple defective copies containing deletions, frame shifts and stop codons. Ramar1 contains several conserved amino acid blocks which were identified, including a specific D,D(34)D signature motif. Ramar2 is a defective mariner-like element, which contains a deletion overlapping in most of the internal region of the transposase ORF while its extremities remain intact. Predicted transposase sequences demonstrated that Ramar1 and Ramar2 phylogenetically present high identity to mariner-like elements of mauritiana subfamily. Southern blot analysis indicated that Ramar1 is widely represented in the genome of Rhynchosciara americana. In situ hybridizations showed Ramar1 localized in several chromosome regions, mainly in pericentromeric heterochromatin and their boundaries, while Ramar2 appeared as a single band in chromosome A.     

Long-Range Genomic Enrichment,Sequencing, and Assembly to Determine Unknown Sequences Flanking a Known microRNA

Conserved plant microRNAs (miRNAs) modulate important biological processes but little is known about conserved cis-regulatory elements (CREs) surrounding MIRNA genes. We developed a solution-based targeted genomic enrichment methodology to capture, enrich, and sequence flanking genomic regions surrounding conserved MIRNA genes with a locked-nucleic acid (LNA)-modified, biotinylated probe complementary to the mature miRNA sequence. Genomic DNA bound by the probe is captured by streptavidin-coated magnetic beads, amplified, sequenced and assembled de novo to obtain genomic DNA sequences flanking MIRNA locus of interest. We demonstrate the sensitivity and specificity of this enrichment methodology in Arabidopsis thaliana to enrich targeted regions spanning 10–20 kb surrounding known MIR166 and MIR165 loci. Assembly of the sequencing reads successfully recovered all targeted loci. While further optimization for larger, more complex genomes is needed, this method may enable determination of flanking genomic DNA sequence surrounding a known core (like a conserved mature miRNA) from multiple species that currently don''t have a full genome assembly available.     

Structural characterization of Brachypodium genome and its syntenic relationship with rice and wheat

Brachypodium distachyon (Brachypodium) has been recently recognized as an emerging model system for both comparative and functional genomics in grass species. In this study, 55,221 repeat masked Brachypodium BAC end sequences (BES) were used for comparative analysis against the 12 rice pseudomolecules. The analysis revealed that ~26.4% of BES have significant matches with the rice genome and 82.4% of the matches were homologous to known genes. Further analysis of paired-end BES and ~1.0 Mb sequences from nine selected BACs proved to be useful in revealing conserved regions and regions that have undergone considerable genomic changes. Differential gene amplification, insertions/deletions and inversions appeared to be the common evolutionary events that caused variations of microcolinearity at different orthologous genomic regions. It was found that ~17% of genes in the two genomes are not colinear in the orthologous regions. Analysis of BAC sequences also revealed higher gene density (~9 kb/gene) and lower repeat DNA content (~13.1%) in Brachypodium when compared to the orthologous rice regions, consistent with the smaller size of the Brachypodium genome. The 119 annotated Brachypodium genes were BLASTN compared against the wheat EST database and deletion bin mapped wheat ESTs. About 77% of the genes retrieved significant matches in the EST database, while 9.2% matched to the bin mapped ESTs. In some cases, genes in single Brachypodium BACs matched to multiple ESTs that were mapped to the same deletion bins, suggesting that the Brachypodium genome will be useful for ordering wheat ESTs within the deletion bins and developing specific markers at targeted regions in the wheat genome.     

Resolving the breakpoints of the 17q21.31 microdeletion syndrome with next-generation sequencing

Recurrent deletions have been associated with numerous diseases and genomic disorders. Few, however, have been resolved at the molecular level because their breakpoints often occur in highly copy-number-polymorphic duplicated sequences. We present an approach that uses a combination of somatic cell hybrids, array comparative genomic hybridization, and the specificity of next-generation sequencing to determine breakpoints that occur within segmental duplications. Applying our technique to the 17q21.31 microdeletion syndrome, we used genome sequencing to determine copy-number-variant breakpoints in three deletion-bearing individuals with molecular resolution. For two cases, we observed breakpoints consistent with nonallelic homologous recombination involving only H2 chromosomal haplotypes, as expected. Molecular resolution revealed that the breakpoints occurred at different locations within a 145 kbp segment of >99% identity and disrupt KANSL1 (previously known as KANSL1). In the remaining case, we found that unequal crossover occurred interchromosomally between the H1 and H2 haplotypes and that this event was mediated by a homologous sequence that was once again missing from the human reference. Interestingly, the breakpoints mapped preferentially to gaps in the current reference genome assembly, which we resolved in this study. Our method provides a strategy for the identification of breakpoints within complex regions of the genome harboring high-identity and copy-number-polymorphic segmental duplication. The approach should become particularly useful as high-quality alternate reference sequences become available and genome sequencing of individuals'' DNA becomes more routine.     

两个糙皮侧耳菌株的比较基因组学研究

随着糙皮侧耳基因组测序的完成,基因组数据的深入挖掘和利用成为研究重点。本研究基于较高质量的糙皮侧耳CCMSSC00389-1和CCMSSC03989-1基因组图谱,开展了种内比较基因组学分析。CCMSSC00389-1和CCMSSC03989-1的11条染色体序列表现出良好的共线性,且分别含有89.92%和91.68%的保守基因。对菌株特有基因的GO富集分析表明,两菌株各自分化出一些独特的调控方式或途径。CCMSSC00389-1和CCMSSC03989-1之间共存在931 542个单核苷酸多态性位点,231 654个插入/缺失和9 221个结构变异。对与遗传变异重叠基因的GO富集分析表明,碳水化合物降解、物质运输/催化和调控/蛋白活性相关的基因分别容易发生SNP、In/Del和SV变异。CCMSSC00389-1中50.326kb缺失序列的断点两侧含有双链断裂修复（DSB）信号,推测DSB介导50.326kb序列的缺失。     

Identification of transposable element-mediated deletions in 27 Korean individuals based on whole genome sequencing data

The human genome has various genomic structural variations such as insertion/deletions between human individuals. These structural variations have led to genomic fluidity and rearrangements in individuals and populations. To investigate Korean-specific structural genomic variations, we performed next generation sequencing with 30× mean coverage from 27 Korean individuals using illumina-HiSeq 2000 platform. We collected a total of 119 deletion loci as transposable element-mediated Korean-specific deletion (KSD) candidates. Of the 119 loci, 35 were filtered out due to computational overlapping regions. A total of 78 loci were validated by PCR amplification with 27 Korean individuals and 80 human individuals from four different populations. We confirmed deletion breakpoints of the 78 loci using Sanger sequencing. We also investigated different deletion mechanisms based on sequencing alignment analysis. We found at least one KSD locus in 80 human individual panel. It has not been previously reported in human genomes. Here, for the first time, we report transposable element-mediated KSD study based on whole genome sequencing data of 27 Korean.     

Breakpoint junctions of chromosome 9p deletions in two human glioma cell lines.

Interstitial deletions of the short arm of chromosome 9 are associated with glioma, acute lymphoblastic leukemia, melanoma, mesothelioma, lung cancer, and bladder cancer. The distal breakpoints of the deletions (in relation to the centromere) in 14 glioma and leukemia cell lines have been mapped within the 400 kb IFN gene cluster located at band 9p21. To obtain information about the mechanism of these deletions, we have isolated and analyzed the nucleotide sequences at the breakpoint junctions in two glioma-derived cell lines. The A1235 cell line has a complex rearrangement of chromosome 9, including a deletion and an inversion that results in two breakpoint junctions. Both breakpoints of the distal inversion junction occurred within AT-rich regions. In the A172 cell line, a tandem heptamer repeat was found on either side of the deletion breakpoint junction. The distal breakpoint occurred 5' of IFNA2; the 256 bp sequenced from the proximal side of the breakpoint revealed 95% homology to long interspersed nuclear elements. One- and two-base-pair overlaps were observed at these junctions. The possible role of sequence overlaps, and repetitive sequences, in the rearrangement is discussed.     

Genome Rearrangements Detected by SNP Microarrays in Individuals with Intellectual Disability Referred with Possible Williams Syndrome

Background

Intellectual disability (ID) affects 2–3% of the population and may occur with or without multiple congenital anomalies (MCA) or other medical conditions. Established genetic syndromes and visible chromosome abnormalities account for a substantial percentage of ID diagnoses, although for ∼50% the molecular etiology is unknown. Individuals with features suggestive of various syndromes but lacking their associated genetic anomalies pose a formidable clinical challenge. With the advent of microarray techniques, submicroscopic genome alterations not associated with known syndromes are emerging as a significant cause of ID and MCA.

Methodology/Principal Findings

High-density SNP microarrays were used to determine genome wide copy number in 42 individuals: 7 with confirmed alterations in the WS region but atypical clinical phenotypes, 31 with ID and/or MCA, and 4 controls. One individual from the first group had the most telomeric gene in the WS critical region deleted along with 2 Mb of flanking sequence. A second person had the classic WS deletion and a rearrangement on chromosome 5p within the Cri du Chat syndrome (OMIM:123450) region. Six individuals from the ID/MCA group had large rearrangements (3 deletions, 3 duplications), one of whom had a large inversion associated with a deletion that was not detected by the SNP arrays.

Conclusions/Significance

Combining SNP microarray analyses and qPCR allowed us to clone and sequence 21 deletion breakpoints in individuals with atypical deletions in the WS region and/or ID or MCA. Comparison of these breakpoints to databases of genomic variation revealed that 52% occurred in regions harboring structural variants in the general population. For two probands the genomic alterations were flanked by segmental duplications, which frequently mediate recurrent genome rearrangements; these may represent new genomic disorders. While SNP arrays and related technologies can identify potentially pathogenic deletions and duplications, obtaining sequence information from the breakpoints frequently provides additional information.     

Structural analysis of aberrant chromosomes that occur spontaneously in diploid Saccharomyces cerevisiae: retrotransposon Ty1 plays a crucial role in chromosomal rearrangements.

The structural analysis of aberrant chromosomes is important for our understanding of the molecular mechanisms underlying chromosomal rearrangements. We have identified a number of diploid Saccharomyces cerevisiae clones that have undergone loss of heterozygosity (LOH) leading to functional inactivation of the hemizygous URA3 marker placed on the right arm of chromosome III. Aberrant-sized chromosomes derived from chromosome III were detected in approximately 8% of LOH clones. Here, we have analyzed the structure of the aberrant chromosomes in 45 LOH clones with a PCR-based method that determines the ploidy of a series of loci on chromosome III. The alterations included various deletions and amplifications. Sequencing of the junctions revealed that all the breakpoints had been made within repeat sequences in the yeast genome, namely, MAT-HMR, which resulted in intrachromosomal deletion, and retrotransposon Ty1 elements, which were involved in various translocations. Although the translocations involved different breakpoints on different chromosomes, all breakpoints were exclusively within Ty1 elements. Some of the resulting Ty1 elements left at the breakpoints had a complex construction that indicated the involvement of other Ty1 elements not present at the parental breakpoints. These indicate that Ty1 elements are crucially involved in the generation of chromosomal rearrangements in diploid yeast cells.