Duplication and DNA segmental loss in the rice genome: implications for diploidization

* Large-scale duplication events have been recently uncovered in the rice genome, but different interpretations were proposed regarding the extent of the duplications. * Through analysing the 370 Mb genome sequences assembled into 12 chromosomes of Oryza sativa subspecies indica, we detected 10 duplicated blocks on all 12 chromosomes that contained 47% of the total predicted genes. Based on the phylogenetic analysis, we inferred that this was a result of a genome duplication that occurred c. 70 million years ago, supporting the polyploidy origin of the rice genome. In addition, a segmental duplication was also identified involving chromosomes 11 and 12, which occurred c. 5 million years ago. * Following the duplications, there have been large-scale chromosomal rearrangements and deletions. About 30-65% of duplicated genes were lost shortly after the duplications, leading to a rapid diploidization. * Together with other lines of evidence, we propose that polyploidization is still an ongoing process in grasses of polyploidy origins.     

Molecular Evidence for Asymmetric Evolution of Sister Duplicated Blocks after Cereal Polyploidy

Polyploidy (genome duplication) is thought to have contributed to the evolution of the eukaryotic genome, but complex genome structures and massive gene loss during evolution has complicated detection of these ancestral duplication events. The major factors determining the fate of duplicated genes are currently unclear, as are the processes by which duplicated genes evolve after polyploidy. Fine-scale analysis between homologous regions may allow us to better understand post-polyploidy evolution. Here, using gene-by-gene and gene-by-genome strategies, we identified the S5 region and four homologous regions within the japonica genome. Additional phylogenomic analyses of the comparable duplicated blocks indicate that four successive duplication events gave rise to these five regions, allowing us to propose a model for this local chromosomal evolution. According to this model, gene loss may play a major role in post-duplication genetic evolution at the segmental level. Moreover, we found molecular evidence that one of the sister duplicated blocks experienced more gene loss and a more rapid evolution subsequent to two recent duplication events. Given that these two recent duplication events were likely involved in polyploidy, this asymmetric evolution (gene loss and gene divergence) may be one possible mechanism accounting for the diploidization at the segmental level. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s11103-005-4414-1     

Segmental duplications and copy-number variation in the human genome

The human genome contains numerous blocks of highly homologous duplicated sequence. This higher-order architecture provides a substrate for recombination and recurrent chromosomal rearrangement associated with genomic disease. However, an assessment of the role of segmental duplications in normal variation has not yet been made. On the basis of the duplication architecture of the human genome, we defined a set of 130 potential rearrangement hotspots and constructed a targeted bacterial artificial chromosome (BAC) microarray (with 2,194 BACs) to assess copy-number variation in these regions by array comparative genomic hybridization. Using our segmental duplication BAC microarray, we screened a panel of 47 normal individuals, who represented populations from four continents, and we identified 119 regions of copy-number polymorphism (CNP), 73 of which were previously unreported. We observed an equal frequency of duplications and deletions, as well as a 4-fold enrichment of CNPs within hotspot regions, compared with control BACs (P < .000001), which suggests that segmental duplications are a major catalyst of large-scale variation in the human genome. Importantly, segmental duplications themselves were also significantly enriched >4-fold within regions of CNP. Almost without exception, CNPs were not confined to a single population, suggesting that these either are recurrent events, having occurred independently in multiple founders, or were present in early human populations. Our study demonstrates that segmental duplications define hotspots of chromosomal rearrangement, likely acting as mediators of normal variation as well as genomic disease, and it suggests that the consideration of genomic architecture can significantly improve the ascertainment of large-scale rearrangements. Our specialized segmental duplication BAC microarray and associated database of structural polymorphisms will provide an important resource for the future characterization of human genomic disorders.     

Divergence of genes encoding non-specific lipid transfer proteins in the poaceae family

The genes encoding non-specific lipid transfer proteins (nsLTPs), members of a small multigene family, show a complex pattern of expressional regulation, suggesting that some diversification may have resulted from changes in their expression after duplication. In this study, the evolution of nsLTP genes within the Poaceae family was characterized via a survey of the pseudogenes and unigenes encoding the nsLTP in rice pseudomolecules and the NCBI unigene database. nsLTP-rich regions were detected in the distal portions of rice chromosomes 11 and 12; these may have resulted from the most recent large segmental duplication in the rice genome. Two independent tandem duplications were shown to occur within the nsLTP-rich regions of rice. The genomic distribution of the nsLTP genes in the rice genome differs from that in wheat. This may be attributed to gene migration, chromosomal rearrangement, and/or differential gene loss. The genomic distribution pattern of nsLTP genes in the Poaceae family points to the existence of some differences among cereal nsLTP genes, all of which diverged from an ancient gene. The unigenes encoding nsLTPs in each cereal species are clustered into five groups. The somewhat different distribution of nsLTP-encoding EST clones between the groups across cereal species imply that independent duplication(s) followed by subfunctionalization (and/or neofunctionalization) of the nsLTP gene family in each species occurred during speciation.     

Genome change in wheat observed through the structure and expression of α/β-gliadin genes

To better understand genome structure and the expression of α/β-gliadin multigenes in hexaploid wheat, bacterial artificial chromosome (BAC) clones containing α/β-gliadin genes from the three loci, Gli-A2, Gli-B2, and Gli-D2, were screened. Based on their restriction fragment patterns, we selected five BAC clones, namely, two clones for Gli-A2, two clones for Gli-B2, and one clone for Gli-D2, to fully sequence. Approximately 200 kb was sequenced for each locus. In total, twelve α/β-gliadin intact genes and four pseudogenes were found, and retrotransposons or other transposons existed in each BAC clone. Dot-plot analysis revealed the pattern of genome segmental duplication within each BAC. We calculated time since duplication of each set of α/β-gliadin genes and insertion of retrotransposons. Duplication of all adjacent genes within the same BAC clone took place before or after allotetrapolyploidization, but duplication of certain genes occurred before diploid differentiation of wheat species. Retrotransposons were also inserted before and after the segmental duplication events. Furthermore, translocation of α/β-gliadin genes from chromosomes 1 to 6 apparently occurred before the diversification of various wheat genomes. Duplication of genome segments containing α/β-gliadin genes and retrotransposons were brought about through unequal crossing-over or saltatory replication and α/β-gliadin genes per se were duplicated without any recombination events. Out of twelve intact α/β-gliadin genes detected from their sequences, nine were expressed, although their patterns of expression were distinct. Since they have similar cis-elements and promoter structures, the mechanisms underlying their distinct gene expression and possible applications are discussed.     

Genome-level evolution of resistance genes in Arabidopsis thaliana

Pathogen resistance genes represent some of the most abundant and diverse gene families found within plant genomes. However, evolutionary mechanisms generating resistance gene diversity at the genome level are not well understood. We used the complete Arabidopsis thaliana genome sequence to show that most duplication of individual NBS-LRR sequences occurs at close physical proximity to the parent sequence and generates clusters of closely related NBS-LRR sequences. Deploying the statistical strength of phylogeographic approaches and using chromosomal location as a proxy for spatial location, we show that apparent duplication of NBS-LRR genes to ectopic chromosomal locations is largely the consequence of segmental chromosome duplication and rearrangement, rather than the independent duplication of individual sequences. Although accounting for a smaller fraction of NBS-LRR gene duplications, segmental chromosome duplication and rearrangement events have a large impact on the evolution of this multigene family. Intergenic exchange is dramatically lower between NBS-LRR sequences located in different chromosome regions as compared to exchange between sequences within the same chromosome region. Consequently, once translocated to new chromosome locations, NBS-LRR gene copies have a greater likelihood of escaping intergenic exchange and adopting new functions than do gene copies located within the same chromosomal region. We propose an evolutionary model that relates processes of genome evolution to mechanisms of evolution for the large, diverse, NBS-LRR gene family.     

Theoretical and practical advances in genome halving

MOTIVATION: Duplication of an organism's entire genome is a rare but spectacular event, enabling the rapid emergence of multiple new gene functions. Over time, the parallel linkage of duplicated genes across chromosomes may be disrupted by reciprocal translocations, while the intra-chromosomal order of genes may be shuffled by inversions and transpositions. Some duplicate genes may evolve unrecognizably or be deleted. As a consequence, the only detectable signature of an ancient duplication event in a modern genome may be the presence of various chromosomal segments containing parallel paralogous genes, with each segment appearing exactly twice in the genome. The problem of reconstructing the linkage structure of an ancestral genome before duplication is known as genome halving with unordered chromosomes. RESULTS: In this paper, we derive a new upper bound on the genome halving distance that is tighter than the best known, and a new lower bound that is almost always tighter than the best known. We also define the notion of genome halving diameter, and obtain both upper and lower bounds for it. Our tighter bounds on genome halving distance yield a new algorithm for reconstructing an ancestral duplicated genome. We create a software package GenomeHalving based on this new algorithm and test it on the yeast genome, identifying a sequence of translocations for halving the yeast genome that is shorter than previously conjectured possible.     

Patterns of segmental duplication in the human genome

We analyzed the completed human genome for recent segmental duplications (size > or = 1 kb and sequence similarity > or = 90%). We found that approximately 4% of the genome is covered by duplications and that the extent of segmental duplication varies from 1% to 14% among the 24 chromosomes. Intrachromosomal duplication is more frequent than interchromosomal duplication in 15 chromosomes. The duplication frequencies in pericentromeric and subtelomeric regions are greater than the genome average by approximately threefold and fourfold. We examined factors that may affect the frequency of duplication in a region. Within individual chromosomes, the duplication frequency shows little correlation with local gene density, repeat density, recombination rate, and GC content, except chromosomes 7 and Y. For the entire genome, the duplication frequency is correlated with each of the above factors. Based on known genes and Ensembl genes, the proportion of duplications containing complete genes is 3.4% and 10.7%, respectively. The proportion of duplications containing genes is higher in intrachromosomal than in interchromosomal duplications, and duplications containing genes have a higher sequence similarity and tend to be longer than duplications containing no genes. Our simulation suggests that many duplications containing genes have been selectively maintained in the genome.     

杨树基因组中染色体基因丰度分化及EST序列对基因覆盖度的检测(英文)

高等植物基因组中,大部分序列为非表达序列,基因序列所占的比例很小,了解基因在基因组中的分布是研究基因组结构的一个重要方面。在美国能源部资助下,一个毛果杨无性系的基因组测序已经完成并对公众发布。杨树全基因组序列的完成,为我们了解林木基因组中基因的分布提供了一个特例。在本文中,我们利用泊松分析对杨树基因组中基因在各个染色体上的密度进行了检测,结果表明杨树基因组中各条染色体的基因含量存在显著差异。杨树全基因组测序项目揭示现代杨树基因组起源于一次古全基因组复制事件(称为杨柳科基因组复制),所以杨树基因组不同染色体间存在很大的同源复制片段。但是我们的研究显示,杨树基因组中大多数高度同源的染色体上基因的密度与染色体间的同源性没有明显关系,这说明杨柳科全基因组复制事件后,各个高度同源染色体上的基因发生了流失,且基因流失的速率是不一样的。同时本文还对近九万条毛果杨EST序列进行了比对分析,结果显示这些EST序列覆盖的基因仅占杨树基因组中基因总数的16.8%左右。EST测序虽然是发现基因的一个重要手段,但小规模EST测序对基因的覆盖度很低,所以小规模EST测序的应用价值是有限的。     

Incongruent evolution of chromosomal size in rice

To investigate genome size evolution, it is usually informative to compare closely related species that vary dramatically in genome size. A whole genome duplication (polyploidy) that occurred in rice (Oryza sativa) about 70 million years ago has been well documented based on current genome sequencing. The presence of three distinct duplicate blocks from the polyploidy, of which one duplicated segment in a block is intact (no sequencing gap) and less than half the length of its syntenic duplicate segment, provided an excellent opportunity for elucidating the causes of their size variation during the post-polyploid time. The results indicated that incongruent patterns (shrunken, balanced and inflated) of chromosomal size evolution occurred in the three duplicate blocks, spanning over 30 Mb among chromosomes 2, 3, 6, 7, and 10, with an average of 20.3% for each. DNA sequences of chromosomes 2 and 3 appeared to had become as short as about half of their initial sequence lengths, chromosomes 6 and 7 had remained basically balanced, and chromosome 10 had become dramatically enlarged (approximately 70%). The size difference between duplicate segments of rice was mainly caused by variations in non-repetitive DNA loss. Amplification of long terminal repeat retrotransposons also played an important role. Moreover, a relationship seems to exist between the chromosomal size differences and the nonhomologous combination in corresponding regions in the rice genome. These findings help shed light on the evolutionary mechanism of genomic sequence variation after polyploidy and genome size evolution.     

Inferences on the genome structure of progenitor maize through comparative analysis of rice, maize and the domesticated panicoids.

Corn and rice genetic linkage map alignments were extended and refined by the addition of 262 new, reciprocally mapped maize cDNA loci. Twenty chromosomal rearrangements were identified in maize relative to rice and these included telomeric fusions between rice linkage groups, nested insertion of rice linkage groups, intrachromosomal inversions, and a nonreciprocal translocation. Maize genome evolution was inferred relative to other species within the Panicoideae and a progenitor maize genome with eight linkage groups was proposed. Conservation of composite linkage groups indicates that the tetrasomic state arose during maize evolution either from duplication of one progenitor corn genome (autoploidy) or from a cross between species that shared the composite linkages observed in modern maize (alloploidy). New evidence of a quadruplicated homeologous segment on maize chromosomes 2 and 10, and 3 and 4, corresponded to the internally duplicated region on rice chromosomes 11 and 12 and suggested that this duplication in the rice genome predated the divergence of the Panicoideae and Oryzoideae subfamilies. Charting of the macroevolutionary steps leading to the modern maize genome clarifies the interpretation of intercladal comparative maps and facilitates alignments and genomic cross-referencing of genes and phenotypes among grass family members.     

Chromosomal breakpoint mapping by arrayCGH using flow-sorted chromosomes

Despite the recent completion of the human genome project, the mapping of disease-related chromosomal translocation breakpoints and genes has remained laborious. Here, we describe a novel and rapid procedure to map such translocation breakpoints using flow-sorted chromosomes in combination with array-based comparative genomic hybridization (arrayCGH). To test the feasibility of this approach, we used a t(12;15)(q13;q25)-positive cell line with known breakpoint positions as a model. The derivative 12 chromosomes were flow-sorted, labeled, and hybridized to a genome-wide array containing 3648 well-characterized human genomic clones. The exact locations of the breakpoints on both chromosome 12 and 15 could be determined in a single hybridization experiment. In addition, we have tested the minimal amount of material necessary to perform these experiments and show that it is possible to obtain highly reliable profiles using as little as 10,000 flow-sorted chromosomes.     

Genome duplication events and functional reduction of ploidy levels in sturgeon (Acipenser, Huso and Scaphirhynchus).

Sturgeon (order Acipenserformes) provide an ideal taxonomic context for examination of genome duplication events. Multiple levels of ploidy exist among these fish. In a novel microsatellite approach, data from 962 fish from 20 sturgeon species were used for analysis of ploidy in sturgeon. Allele numbers in a sample of individuals were assessed at six microsatellite loci. Species with approximately 120 chromosomes are classified as functional diploid species, species with approximately 250 chromosomes as functional tetraploid species, and with approximately 500 chromosomes as functional octaploids. A molecular phylogeny of the sturgeon was determined on the basis of sequences of the entire mitochondrial cytochrome b gene. By mapping the estimated levels of ploidy on this proposed phylogeny we demonstrate that (I) polyploidization events independently occurred in the acipenseriform radiation; (II) the process of functional genome reduction is nearly finished in species with approximately 120 chromosomes and more active in species with approximately 250 chromosomes and approximately 500 chromosomes; and (III) species with approximately 250 and approximately 500 chromosomes arose more recently than those with approximately 120 chromosomes. These results suggest that gene silencing, chromosomal rearrangements, and transposition events played an important role in the acipenseriform genome formation. Furthermore, this phylogeny is broadly consistent with previous hypotheses but reveals a highly supported oceanic (Atlantic-Pacific) subdivision within the Acipenser/Huso complex.     

A rhesus macaque radiation hybrid map and comparative analysis with the human genome

The genomes of nonhuman primates are powerful references for better understanding the recent evolution of the human genome. Here we compare the order of 802 genomic markers mapped in a rhesus macaque (Macaca mulatta) radiation hybrid panel with the human genome, allowing for nearly complete cross-reference to the human genome at an average resolution of 3.5 Mb. At least 23 large-scale chromosomal rearrangements, mostly inversions, are needed to explain the changes in marker order between human and macaque. Analysis of the breakpoints flanking inverted chromosomal segments and estimation of their duplication divergence dates provide additional evidence implicating segmental duplications as a major mechanism of chromosomal rearrangement in recent primate evolution.     

Detection of gene duplications and block duplications in eukaryotic genomes

Several eukaryotic genomes have been completely sequenced and this provides an opportunity to investigate the extent and characteristics (e.g., single gene duplication, block duplication, etc.) of gene duplication in a genome. Detecting duplicate genes in a genome, however, is not a simple problem because of several complications such as domain shuffling, the existence of isoforms derived from alternative splicing, and annotational errors in the databases. We describe a method for overcoming these difficulties and the extents of gene duplication in the genomes of Drosophila melanogaster, Caenorhabditis elegans, and yeast inferred from this method. We also describe a method for detecting block duplications in a genome. Application of this method showed that block duplication is a common phenomenon in both yeast and nematode. The patterns of block duplication in the two species are, however, markedly different. Yeast shows much more extensive block duplication than nematode, with some chromosomes having more than 40% of the duplications derived from block duplications. Moreover, in yeast the majority of block duplications occurred between chromosomes, while in nematode most block duplications occurred within chromosomes.     

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.     

Directed Synthesis of a Segmental Chromosomal Transposition: An Approach to the Study of Chromosomes Lethal to the Gametophyte Generation of Maize

Because of the haploid nature of the gametophyte generation of plants, most mutations that are lethal or detrimental to the gametophytes cannot be recovered. Our laboratory is currently developing several techniques to overcome this situation. In this paper, a procedure is described to generate directed segmental chromosomal transpositions. The method involves recovery of recombinants between reciprocal translocation overlaps such that one region of the genome is inserted into a nonhomologous chromosome in a predetermined and directed manner. This duplicated segment then could serve to cover deficiencies or mutations, lethal to the gametophytes, in the region from whence it originated. The manipulation of segmental chromosomal transpositions for analyzing mutants lethal or detrimental to the gametophyte generation is discussed. The procedure to generate transpositions, the translocations between normal A and supernumerary B chromosomes that generate deficiencies in the male gametes, the r-X1 chromosome that generates deficiencies in the female gametes and other techniques available in maize form a system to analyze gametophyte lethal mutations.     

Construction of an SfiI macrorestriction map of the Candida albicans genome.

The opportunistic fungal pathogen, Candida albicans, is diploid as usually isolated and has no apparent sexual cycle. Genetic analysis has therefore been very difficult. Molecular genetics has yielded important information in the past few years, but it too is hampered by the lack of a good genetic map. Using the well-characterized strain 1006 and strain WO-1, which undergoes the white-opaque phenotypic transition, we have developed a genomic restriction map of C. albicans with the enzyme SfiI. There are approximately 34 SfiI restriction sites in the C. albicans genome. Restriction fragments were separated by pulsed-field electrophoresis and were assigned to chromosomes by hybridization of complete and partial digests with known chromosome-specific probes as well as by digestion of isolated chromosomes. Telomeric fragments were identified by hybridization with a telomere-specific probe (C. Sadhu, M.J. McEachern, E.P. Rustchenko-Bulgac, J. Schmid, D.R. Soll, and J.B. Hicks, J. Bacteriol. 173:842-850, 1991). WO-1 differs from 1006 in that it has undergone three reciprocal chromosomal translocations. Analysis of the translocation products indicates that each translocation has occurred at or near an SfiI site; thus, the SfiI fragments from the two strains are similar or identical. The tendency for translocation to occur at or near SfiI sites may be related to the repeated sequence RPS 1, which contains four such sites and could provide homology for ectopic pairing and crossing over. The genome size of both strains is about 16 to 17 megabases, in good agreement with previous determinations.     

Genome Resilience and Prevalence of Segmental Duplications Following Fast Neutron Irradiation of Soybean

Fast neutron radiation has been used as a mutagen to develop extensive mutant collections. However, the genome-wide structural consequences of fast neutron radiation are not well understood. Here, we examine the genome-wide structural variants observed among 264 soybean [Glycine max (L.) Merrill] plants sampled from a large fast neutron-mutagenized population. While deletion rates were similar to previous reports, surprisingly high rates of segmental duplication were also found throughout the genome. Duplication coverage extended across entire chromosomes and often prevailed at chromosome ends. High-throughput resequencing analysis of selected mutants resolved specific chromosomal events, including the rearrangement junctions for a large deletion, a tandem duplication, and a translocation. Genetic mapping associated a large deletion on chromosome 10 with a quantitative change in seed composition for one mutant. A tandem duplication event, located on chromosome 17 in a second mutant, was found to cosegregate with a short petiole mutant phenotype, and thus may serve as an example of a morphological change attributable to a DNA copy number gain. Overall, this study provides insight into the resilience of the soybean genome, the patterns of structural variation resulting from fast neutron mutagenesis, and the utility of fast neutron-irradiated mutants as a source of novel genetic losses and gains.