Segmental duplication density decrease with distance to human-mouse breaks of synteny

Mappability and read length

Power-law distributions are the main functional form for the distribution of repeat size and repeat copy number in the human genome. When the genome is broken into fragments for sequencing, the limited size of fragments and reads may prevent an unique alignment of repeat sequences to the reference sequence. Repeats in the human genome can be as long as 10⁴ bases, or 10⁵ − 10⁶ bases when allowing for mismatches between repeat units. Sequence reads from these regions are therefore unmappable when the read length is in the range of 10³ bases. With a read length of 1000 bases, slightly more than 1% of the assembled genome, and slightly less than 1% of the 1 kb reads, are unmappable, excluding the unassembled portion of the human genome (8% in GRCh37/hg19). The slow decay (long tail) of the power-law function implies a diminishing return in converting unmappable regions/reads to become mappable with the increase of the read length, with the understanding that increasing read length will always move toward the direction of 100% mappability.

Mouse segmental duplication and copy number variation

Detailed analyses of the clone-based genome assembly reveal that the recent duplication content of mouse (4.94%) is now comparable to that of human (5.5%), in contrast to previous estimates from the whole-genome shotgun sequence assembly. However, the architecture of mouse and human genomes differs markedly: most mouse duplications are organized into discrete clusters of tandem duplications that show depletion of genes and transcripts and enrichment of long interspersed nuclear element (LINE) and long terminal repeat (LTR) retroposons. We assessed copy number variation of the C57BL/6J duplicated regions within 15 mouse strains previously used for genetic association studies, sequencing and the Mouse Phenome Project. We determined that over 60% of these base pairs are polymorphic among the strains (on average, there was 20 Mb of copy-number-variable DNA between different mouse strains). Our data suggest that different mouse strains show comparable, if not greater, copy number polymorphism when compared to human; however, such variation is more locally restricted. We show large and complex patterns of interstrain copy number variation restricted to large gene families associated with spermatogenesis, pregnancy, viviparity, pheromone signaling and immune response.

DupMasker: a tool for annotating primate segmental duplications

Segmental duplications (SDs) play an important role in genome rearrangement, evolution, and the copy-number variation (CNV) of primate genomes. Such sequences are difficult to detect, a priori, because they share no defining sequence features that distinguish them from unique portions of the genome. Current sequence annotation of segmental duplications requires computationally intensive, genome-wide self-comparisons that cannot be easily implemented on new data sets. Based on the successful implementation of RepeatMasker, we developed a new genome annotation tool, DupMasker. The program uses a library of nonredundant consensus sequences of human segmental duplications, wherein a majority of the ancestral origins have been determined based on comparisons to mammalian outgroup genomes. Using DupMasker, new human and nonhuman primate (NHP) sequences may be readily queried to provide details on the origin and degree of sequence identity of each duplicon. This program can be applied to delineate the order and orientation of duplicons within complex duplication blocks and used to characterize structural variation differences between sequenced human haplotypes. We predict this tool will be valuable in the annotation of large-insert sequence clones, allowing putative unique and duplicated regions of the genomes to be annotated prior to whole genome assembly comparisons.

Human and mouse protein-noncoding snoRNA host genes with dissimilar nucleotide sequences show chromosomal synteny

snoRNAs are small protein-noncoding RNAs essential for pre-rRNA processing and ribosome biogenesis, and are encoded intronically in host genes (HGs) that are either protein coding or noncoding. mRNAs of protein-noncoding HGs differ in their nucleotide sequences among species. Although the reason for such sequential divergence has not been well explained, we present evidence here that such structurally different HGs have evolved from a common ancestral gene. We first identified two novel protein-noncoding HGs (mU50HG-a and mU50HG-b) that intronically encode a mouse ortholog of a human snoRNA, hU50. The sequences of mU50HG mRNA differed from that of hU50HG. However, a chromosome mapping study revealed that mU50HG is located at 9E3-1, the murine segment syntenic to human 6q15, where hU50HG is located. Synteny is a phenomenon whereby gene orthologs are arranged in the same order at equivalent chromosomal loci in different species; synteny between two species means it is highly likely that the genes have evolved from a common ancestral gene. We then extended this mapping study to other protein-noncoding snoRNA-HGs, and found again that they are syntenic, implying that they have evolved from genes of common ancestral species. Furthermore, on these syntenic segments, exons of adjacent protein-coding genes were found to be far better conserved than those of noncoding HGs, suggesting that the exons of protein-noncoding snoRNA-HGs have been much more fragile during evolution.