Pericentromeric regions of soybean (Glycine max L. Merr.) chromosomes consist of retroelements and tandemly repeated DNA and are structurally and evolutionarily labile

Chromosome and Genome Diversity in the Genus Trifolium (Fabaceae)

Trifolium L. is an economically important genus that is characterized by variable karyotypes relating to its ploidy level and basic chromosome numbers. The advent of genomic resources combined with molecular cytogenetics provides an opportunity to develop our understanding of plant genomes in general. Here, we summarize the current state of knowledge on Trifolium genomes and chromosomes and review methodologies using molecular markers that have contributed to Trifolium research. We discuss possible future applications of cytogenetic methods in research on the Trifolium genome and chromosomes.

Plasticity in Triticeae centromere DNA sequences: a wheat × tall wheatgrass (decaploid) model

Centromeres mediate chromosome attachment to microtubules and maintain the integrity of chromosomes for proper segregation of the sister chromatids during cell division. Advances in the assembly of Triticeae genome sequences combined with the capacity to recover hybrid species derived from very distantly related species provides potential experimental systems for linking retrotransposon amplification and repositioning of centromeres via non-mendelian inheritance in partial amphiploid breeds. The decaploid tall wheatgrass (Thinopyrum ponticum) is one of the most successfully used perennial species in wheat breeding for generating translocation lines with valuable agronomic traits. We found that wheat centromere retrotransposons CRW and Quinta widely occur within the tall wheatgrass genome. In addition, one of the genome donors to Th. ponticum, Pseudoroegneria stipifolia (StSt), has been shown to have Abigail and a satellite repeat, CentSt. We also found two other centromeric retrotransposons, Abia and CL135 in Th. ponticum by ChIP-seq. Examination of partial amphiploid lines that were generated in the 1970s demonstrated extensive modification in centromere sequences using CentSt, Abigail and Abia as probes. We also detected that St-genome chromosomes were more enriched with Abigail and CentSt, whereas E-genome chromosomes were enriched with CRW and Quinta in tall wheatgrass and its closer relatives. It can be concluded that bursts of transposition of retrotransposons and repositioning of centromeres via non-mendelian segregation are common in partial amphiploids derived from interspecific hybrids. Practically speaking, our study reveals that the existence of homologous centromere functional sequences in both a donor and its receptor can substantially contribute to the successful transfer of alien genes into crop species. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/sra/SRR9089557; https://www.ncbi.nlm.nih.gov/sra/SRR9089558; https://www.ncbi.nlm.nih.gov/sra/SRR9089559; https://www.ncbi.nlm.nih.gov/sra/SRR9089560; https://www.ncbi.nlm.nih.gov/sra/SRR9089561; https://www.ncbi.nlm.nih.gov/sra/SRR9089562; https://www.ncbi.nlm.nih.gov/sra/SRR9089563; https://www.ncbi.nlm.nih.gov/sra/SRR9089564; https://www.ncbi.nlm.nih.gov/nuccore/MK999394; https://www.ncbi.nlm.nih.gov/nuccore/MK999395; https://www.ncbi.nlm.nih.gov/nuccore/MK999396.

Genome-wide nucleotide patterns and potential mechanisms of genome divergence following domestication in maize and soybean

BACKGROUND

Plant domestication provides a unique model to study genome evolution. Many studies have been conducted to examine genes, genetic diversity, genome structure, and epigenome changes associated with domestication. Interestingly, domesticated accessions have significantly higher [A] and [T] values across genome-wide polymorphic sites than accessions sampled from the corresponding progenitor species. However, the relative contributions of different genomic regions to this genome divergence pattern and underlying mechanisms have not been well characterized.

RESULTS

Here, we investigate the genome-wide base-composition patterns by analyzing millions of SNPs segregating among 100 accessions from a teosinte-maize comparison set and among 302 accessions from a wild-domesticated soybean comparison set. We show that non-genic part of the genome has a greater contribution than genic SNPs to the [AT]-increase observed between wild and domesticated accessions in maize and soybean. The separation between wild and domesticated accessions in [AT] values is significantly enlarged in non-genic and pericentromeric regions. Motif frequency and sequence context analyses show the motifs (PyCG) related to solar-UV signature are enriched in these regions, particularly when they are methylated. Additional analysis using population-private SNPs also implicates the role of these motifs in relatively recent mutations. With base-composition across polymorphic sites as a genome phenotype, genome scans identify a set of putative candidate genes involved in UV damage repair pathways.

CONCLUSIONS

The [AT]-increase is more pronounced in genomic regions that are non-genic, pericentromeric, transposable elements; methylated; and with low recombination. Our findings establish important links among UV radiation, mutation, DNA repair, methylation, and genome evolution.

Organization and evolution of two repetitive sequences, 18-24J and 12-13P, in the genome of Chenopodium (Amaranthaceae)

The abundance and chromosomal organization of two repetitive sequences named 12-13P and 18-24J were analyzed in 24 diploid and nine polyploid species of Chenopodium s.l., with special attention to Chenopodium s.s. Both sequences were predominantly present in species of Chenopodium s.s.; however, differences in the amplification levels were observed among the species. The 12-13P repeat was highly amplified in all of the analyzed Eurasian species, whereas the American diploids showed a marked variation in the amplification levels. The 12-13P repeat contains a tandemly arranged 40 bp minisatellite element forming a large proportion of the genome of Chenopodium (up to 3.5%). FISH revealed its localization to the pericentromeric regions of the chromosomes. The chromosomal distribution of 12-13P delivered additional chromosomal marker for B-genome diploids. The 18-24J repeat showed a dispersed organization in all of the chromosomes of the analyzed diploid species and the Eurasian tetraploids. In the American allotetraploids (C. quinoa, C. berlandieri) and Eurasian allohexaploids (e.g., C. album) very intense hybridization signals of 18-24J were observed only on 18 chromosomes that belong to the B subgenome of these polyploids. Combined cytogenetic and molecular analyses suggests that reorganization of these two repeats accompanied the diversification and speciation of diploid (especially A genome) and polyploid species of Chenopodium s.s.

Assembly and annotation of a draft genome sequence for Glycine latifolia, a perennial wild relative of soybean

Glycine latifolia (Benth.) Newell & Hymowitz (2n = 40), one of the 27 wild perennial relatives of soybean, possesses genetic diversity and agronomically favorable traits that are lacking in soybean. Here, we report the 939-Mb draft genome assembly of G. latifolia (PI 559298) using exclusively linked-reads sequenced from a single Chromium library. We organized scaffolds into 20 chromosome-scale pseudomolecules utilizing two genetic maps and the Glycine max (L.) Merr. genome sequence. High copy numbers of putative 91-bp centromere-specific tandem repeats were observed in consecutive blocks within predicted pericentromeric regions on several pseudomolecules. No 92-bp putative centromeric repeats, which are abundant in G. max, were detected in G. latifolia or Glycine tomentella. Annotation of the assembled genome and subsequent filtering yielded a high confidence gene set of 54 475 protein-coding loci. In comparative analysis with five legume species, genes related to defense responses were significantly overrepresented in Glycine-specific orthologous gene families. A total of 304 putative nucleotide-binding site (NBS)-leucine-rich-repeat (LRR) genes were identified in this genome assembly. Different from other legume species, we observed a scarcity of TIR-NBS-LRR genes in G. latifolia. The G. latifolia genome was also predicted to contain genes encoding 367 LRR-receptor-like kinases, a family of proteins involved in basal defense responses and responses to abiotic stress. The genome sequence and annotation of G. latifolia provides a valuable source of alternative alleles and novel genes to facilitate soybean improvement. This study also highlights the efficacy and cost-effectiveness of the application of Chromium linked-reads in diploid plant genome de novo assembly.

In Situ Hybridization in Rice (Oryza sativa)

Fluorescence in situ hybridization (FISH) is widely used in cytogenetics to determine the localization of DNA sequences on target chromosomes, to provide visible information regarding the physical position of DNA sequences, to determine the abundance and distribution of repetitive sequences that comprise a large proportion of genomes, and to determine the relative chromosome positions of multiple sequences in physical mapping. By mapping on extended chromatin fibers, fiber-FISH can be used to determine the structure and organization of genes or DNA sequences with a high resolution (to a few kilobases). The protocols described here will provide procedures of FISH on metaphase chromosomes and extended chromatin fibers of rice (Oryza sativa). © 2016 by John Wiley & Sons, Inc.

Highly distinct chromosomal structures in cowpea (Vigna unguiculata), as revealed by molecular cytogenetic analysis

Cowpea (Vigna unguiculata (L.) Walp) is an important legume, particularly in developing countries. However, little is known about its genome or chromosome structure. We used molecular cytogenetics to characterize the structure of pachytene chromosomes to advance our knowledge of chromosome and genome organization of cowpea. Our data showed that cowpea has highly distinct chromosomal structures that are cytologically visible as brightly DAPI-stained heterochromatic regions. Analysis of the repetitive fraction of the cowpea genome present at centromeric and pericentromeric regions confirmed that two retrotransposons are major components of pericentromeric regions and that a 455-bp tandem repeat is found at seven out of 11 centromere pairs in cowpea. These repeats likely evolved after the divergence of cowpea from common bean and form chromosomal structure unique to cowpea. The integration of cowpea genetic and physical chromosome maps reveals potential regions of suppressed recombination due to condensed heterochromatin and a lack of pairing in a few chromosomal termini. This study provides fundamental knowledge on cowpea chromosome structure and molecular cytogenetics tools for further chromosome studies.

Distribution of new satellites and simple sequence repeats in annual and perennial Glycine species

The repeat sequences occupied more than 50 % of soybean genome. In order to understand where these repeat sequences distributed in soybean genome and its related Glycine species, we examined three new repeat sequences-soybean repeat sequence (SBRS1, SBRS2 and SBRS3), some nonspecific repeat sequences and 45S rDNA on several Glycine species, including annual and perennial accessions in this study. In the annual species, G. soja, signals for SBRS1 and ATT repeat can be found on each chromosome in GG genome, but those for SBRS2 and SBRS3 were located at three specific loci. In perennial Glycine species, these three SBR repeat frequently co-localized with 45S rDNA, two major 45S rDNA loci were found in all tetraploid species. However, an extra minor locus was found in one accession of the G. pescadrensis (Tab074), but not in another accession (Tab004). We demonstrate that some repetitive sequences are present in all Glycine species used in the study, but the abundancy is different in annual or perennial species. We suggest this study may provide additional information in investigations of the phylogeny in the Glycine species.

Punctuated distribution of recombination hotspots and demarcation of pericentromeric regions in Phaseolus vulgaris L

High density genetic maps are a reliable tool for genetic dissection of complex plant traits. Mapping resolution is often hampered by the variable crossover and non-crossover events occurring across the genome, with pericentromeric regions (pCENR) showing highly suppressed recombination rates. The efficiency of linkage mapping can further be improved by characterizing and understanding the distribution of recombinational activity along individual chromosomes. In order to evaluate the genome wide recombination rate in common beans (Phaseolus vulgaris L.) we developed a SNP-based linkage map using the genotype-by-sequencing approach with a 188 recombinant inbred line family generated from an inter gene pool cross (Andean x Mesoamerican). We identified 1,112 SNPs that were subsequently used to construct a robust linkage map with 11 groups, comprising 513 recombinationally unique marker loci spanning 943 cM (LOD 3.0). Comparative analysis showed that the linkage map spanned >95% of the physical map, indicating that the map is almost saturated. Evaluation of genome-wide recombination rate indicated that at least 45% of the genome is highly recombinationally suppressed, and allowed us to estimate locations of pCENRs. We observed an average recombination rate of 0.25 cM/Mb in pCENRs as compared to the rest of genome that showed 3.72 cM/Mb. However, several hot spots of recombination were also detected with recombination rates reaching as high as 34 cM/Mb. Hotspots were mostly found towards the end of chromosomes, which also happened to be gene-rich regions. Analyzing relationships between linkage and physical map indicated a punctuated distribution of recombinational hot spots across the genome.

The wild side of a major crop: soybean's perennial cousins from Down Under

The accumulation of over 30 years of basic research on the biology, genetic variation, and evolution of the wild perennial relatives of soybean (Glycine max) provides a foundation to improve cultivated soybean. The cultivated soybean and its wild progenitor, G. soja, have a center of origin in eastern Asia and are the only two species in the annual subgenus Soja. Systematic and evolutionary studies of the ca. 30 perennial species of subgenus Glycine, native to Australia, have benefited from the availability of the G. max genomic sequence. The perennial species harbor many traits of interest to soybean breeders, among them resistance to major soybean pathogens such as cyst nematode and leaf rust. New species in the Australian subgenus continue to be described, due to the collection of new material and to insights gleaned through systematic studies of accessions in germplasm collections. Ongoing studies in perennial species focus on genomic regions that contain genes for key traits relevant to soybean breeding. These comparisons also include the homoeologous regions that are the result of polyploidy in the common ancestor of all Glycine species. Subgenus Glycine includes a complex of recently formed allopolyploids that are the focus of studies aimed at elucidating genomic, transcriptomic, physiological, taxonomic, morphological, developmental, and ecological processes related to polyploid evolution. Here we review what has been learned over the past 30 years and outline ongoing work on photosynthesis, nitrogen fixation, and floral biology, much of it drawing on new technologies and resources.

Identification of high-quality single-nucleotide polymorphisms in Glycine latifolia using a heterologous reference genome sequence

Like many widely cultivated crops, soybean [Glycine max (L.) Merr.] has a relatively narrow genetic base, while its perennial distant relatives in the subgenus Glycine Willd. are more genetically diverse and display desirable traits not present in cultivated soybean. To identify single-nucleotide polymorphisms (SNPs) between a pair of G. latifolia accessions that were resistant or susceptible to Sclerotinia sclerotiorum (Lib.) de Bary, reduced-representations of DNAs from each accession were sequenced. Approximately 30 % of the 36 million 100-nt reads produced from each of the two G. latifolia accessions aligned primarily to gene-rich euchromatic regions on the distal arms of G. max chromosomes. Because a genome sequence was not available for G. latifolia, the G. max genome sequence was used as a reference to identify 9,303 G. latifolia SNPs that aligned to unique positions in the G. max genome with at least 98 % identity and no insertions and deletions. To validate a subset of the SNPs, nine TaqMan and 384 GoldenGate allele-specific G. latifolia SNP assays were designed and analyzed in F2 G. latifolia populations derived from G. latifolia plant introductions (PI) 559298 and 559300. All nine TaqMan markers and 91 % of the 291 polymorphic GoldenGate markers segregated in a 1:2:1 ratio. Genetic linkage maps were assembled for G. latifolia, nine of which were uninterrupted and nearly collinear with the homoeologous G. max chromosomes. These results made use of a heterologous reference genome sequence to identify more than 9,000 informative high-quality SNPs for G. latifolia, a subset of which was used to generate the first genetic maps for any perennial Glycine species.

Hidden magicians of genome evolution

Transposable elements (TEs) represent genome's dynamic component, causing mutations and genetic variations. Transposable elements can invade eukaryotic genomes in a short span; these are silenced by homology-dependent gene silencing and some functional parts of silenced elements are utilized to perform novel cellular functions. However, during the past two decades, major interest has been focused on the positive contribution of these elements in the evolution of genomes. The interaction between mobile DNAs and their host genomes are quite diverse, ranging from modifications of gene structure to alterations in general genome architecture and can be regarded as hidden magicians in shaping evolution of genomes. Some of the prominent examples that impressively demonstrate the beneficial impact of TEs on host biology over evolutionary time include their role in structure and functions of eukaryotic genomes.

Levels of DNA methylation and histone methylation and acetylation change in root tip cells of soybean seedlings grown at different temperatures

In order to check whether changes in DNA and histone modifications occur in the nuclei of root tip cells of soybean seedlings grown 1) under control conditions (25 °C), 2) subjected to chilling stress (10 °C) and 3) recovered (25 °C) after chilling, measurements of fluorescence intensity with the use of antibodies to heterochromatin as well as to euchromatin markers were carried out. Moreover, the number and sizes of chromocentres were analyzed. The studies showed that during chilling stress the fluorescence intensity for the markers characteristic of heterochromatin increased while for the markers of euchromatin decreased in comparison to the control. After the recovery the converse situation was observed, i.e. increase in fluorescence intensity for euchromatin markers and decrease in heterochromatin markers. The number of chromocentres remained unchanged in the nuclei of all three studied variants. However, differences in the sizes of chromocentres were observed - the highest number of big chromocentres and simultaneously the lowest number of small chromocentres were in the nuclei of stressed plants. Conversely - in the nuclei of recovered plants there were the lowest number of big chromocentres and the highest number of small ones. The treatment of seedlings with the inhibitors of DNA methylation (5-aza-dC) and histone deacetylation (NaBu) also caused changes in fluorescence intensity and chromocentre sizes in soybean nuclei. These results suggest that DNA and histone modification patterns can be altered in soybean nuclei by different growth temperatures and by appropriate inhibitors influencing epigenetic chromatic modifications.

Genomic organization and dynamics of repetitive DNA sequences in representatives of three Fagaceae genera

Oaks, chestnuts, and beeches are economically important species of the Fagaceae. To understand the relationship between these members of this family, a deep knowledge of their genome composition and organization is needed. In this work, we have isolated and characterized several AFLP fragments obtained from Quercus rotundifolia Lam. through homology searches in available databases. Genomic polymorphisms involving some of these sequences were evaluated in two species of Quercus , one of Castanea , and one of Fagus with specific primers. Comparative FISH analysis with generated sequences was performed in interphase nuclei of the four species, and the co-immunolocalization of 5-methylcytosine was also studied. Some of the sequences isolated proved to be genus-specific, while others were present in all the genera. Retroelements, either gypsy-like of the Tat/Athila clade or copia-like, are well represented, and most are dispersed in euchromatic regions of these species with no DNA methylation associated, pointing to an interspersed arrangement of these retroelements with potential gene-rich regions. A particular gypsy-sequence is dispersed in oaks and chestnut nuclei, but its confinement to chromocenters in beech evidences genome restructuring events during evolution of Fagaceae. Several sequences generated in this study proved to be good tools to comparatively study Fagaceae genome organization.

Fluorescence in situ hybridization-based karyotyping of soybean translocation lines

Soybean (Glycine max [L.] Merr.) is a major crop species and, therefore, a major target of genomic and genetic research. However, in contrast to other plant species, relatively few chromosomal aberrations have been identified and characterized in soybean. This is due in part to the difficulty of cytogenetic analysis of its small, morphologically homogeneous chromosomes. The recent development of a fluorescence in situ hybridization -based karyotyping system for soybean has enabled our characterization of most of the chromosomal translocation lines identified to date. Utilizing genetic data from existing translocation studies in soybean, we identified the chromosomes and approximate breakpoints involved in five translocation lines.

A fluorescence in situ hybridization system for karyotyping soybean

The development of a universal soybean (Glycine max [L.] Merr.) cytogenetic map that associates classical genetic linkage groups, molecular linkage groups, and a sequence-based physical map with the karyotype has been impeded due to the soybean chromosomes themselves, which are small and morphologically homogeneous. To overcome this obstacle, we screened soybean repetitive DNA to develop a cocktail of fluorescent in situ hybridization (FISH) probes that could differentially label mitotic chromosomes in root tip preparations. We used genetically anchored BAC clones both to identify individual chromosomes in metaphase spreads and to complete a FISH-based karyotyping cocktail that permitted simultaneous identification of all 20 chromosome pairs. We applied these karyotyping tools to wild soybean, G. soja Sieb. and Zucc., which represents a large gene pool of potentially agronomically valuable traits. These studies led to the identification and characterization of a reciprocal chromosome translocation between chromosomes 11 and 13 in two accessions of wild soybean. The data confirm that this translocation is widespread in G. soja accessions and likely accounts for the semi-sterility found in some G. soja by G. max crosses.

Recent Insights into Mechanisms of Genome Size Change in Plants

Genome sizes vary considerably across all eukaryotes and even among closely related species. The genesis and evolutionary dynamics of that variation have generated considerable interest, as have the patterns of variation themselves. Here we review recent developments in our understanding of genome size evolution in plants, drawing attention to the higher order processes that can influence the mechanisms generating changing genome size.

Functional centromeres in soybean include two distinct tandem repeats and a retrotransposon

The centromere as a kinetochore assembly site is fundamental to the partitioning of genetic material during cell division. In order to determine the functional centromeres of soybean, we characterized the soybean centromere-specific histone H3 (GmCENH3) protein and developed an antibody against the N-terminal end. Using this antibody, we cloned centromere-associated DNA sequences by chromatin immunoprecipitation. Our analyses indicate that soybean centromeres are composed of two distinct satellite repeats (GmCent-1 and GmCent-4) and retrotransposon-related sequences (GmCR). The possible allopolyploid origin of the soybean genome is discussed in view of the centromeric satellite sequences present.

FIDEL-a retrovirus-like retrotransposon and its distinct evolutionary histories in the A- and B-genome components of cultivated peanut

In this paper, we describe a Ty3-gypsy retrotransposon from allotetraploid peanut (Arachis hypogaea) and its putative diploid ancestors Arachis duranensis (A-genome) and Arachis ipaënsis (B-genome). The consensus sequence is 11,223 bp. The element, named FIDEL (Fairly long Inter-Dispersed Euchromatic LTR retrotransposon), is more frequent in the A- than in the B-genome, with copy numbers of about 3,000 (±950, A. duranensis), 820 (±480, A. ipaënsis), and 3,900 (±1,500, A. hypogaea) per haploid genome. Phylogenetic analysis of reverse transcriptase sequences showed distinct evolution of FIDEL in the ancestor species. Fluorescent in situ hybridization revealed disperse distribution in euchromatin and absence from centromeres, telomeric regions, and the nucleolar organizer region. Using paired sequences from bacterial artificial chromosomes, we showed that elements appear less likely to insert near conserved ancestral genes than near the fast evolving disease resistance gene homologs. Within the Ty3-gypsy elements, FIDEL is most closely related with the Athila/Calypso group of retrovirus-like retrotransposons. Putative transmembrane domains were identified, supporting the presence of a vestigial envelope gene. The results emphasize the importance of FIDEL in the evolution and divergence of different Arachis genomes and also may serve as an example of the role of retrotransposons in the evolution of legume genomes in general.

Genome sequence of the palaeopolyploid soybean

Soybean (Glycine max) is one of the most important crop plants for seed protein and oil content, and for its capacity to fix atmospheric nitrogen through symbioses with soil-borne microorganisms. We sequenced the 1.1-gigabase genome by a whole-genome shotgun approach and integrated it with physical and high-density genetic maps to create a chromosome-scale draft sequence assembly. We predict 46,430 protein-coding genes, 70% more than Arabidopsis and similar to the poplar genome which, like soybean, is an ancient polyploid (palaeopolyploid). About 78% of the predicted genes occur in chromosome ends, which comprise less than one-half of the genome but account for nearly all of the genetic recombination. Genome duplications occurred at approximately 59 and 13 million years ago, resulting in a highly duplicated genome with nearly 75% of the genes present in multiple copies. The two duplication events were followed by gene diversification and loss, and numerous chromosome rearrangements. An accurate soybean genome sequence will facilitate the identification of the genetic basis of many soybean traits, and accelerate the creation of improved soybean varieties.

Genome sequence of the palaeopolyploid soybean

Soybean (Glycine max) is one of the most important crop plants for seed protein and oil content, and for its capacity to fix atmospheric nitrogen through symbioses with soil-borne microorganisms. We sequenced the 1.1-gigabase genome by a whole-genome shotgun approach and integrated it with physical and high-density genetic maps to create a chromosome-scale draft sequence assembly. We predict 46,430 protein-coding genes, 70% more than Arabidopsis and similar to the poplar genome which, like soybean, is an ancient polyploid (palaeopolyploid). About 78% of the predicted genes occur in chromosome ends, which comprise less than one-half of the genome but account for nearly all of the genetic recombination. Genome duplications occurred at approximately 59 and 13 million years ago, resulting in a highly duplicated genome with nearly 75% of the genes present in multiple copies. The two duplication events were followed by gene diversification and loss, and numerous chromosome rearrangements. An accurate soybean genome sequence will facilitate the identification of the genetic basis of many soybean traits, and accelerate the creation of improved soybean varieties.

Hypervariable 3' UTR region of plant LTR-retrotransposons as a source of novel satellite repeats

The repetitive sequence PisTR-A has an unusual organization in the pea (Pisum sativum) genome, being present both as short dispersed repeats as well as long arrays of tandemly arranged satellite DNA. Cloning, sequencing and FISH analysis of both PisTR-A variants revealed that the former occurs in the genome embedded within the sequence of Ty3/gypsy-like Ogre elements, whereas the latter forms homogenized arrays of satellite repeats at several genomic loci. The Ogre elements carry the PisTR-A sequences in their 3' untranslated region (UTR) separating the gag-pol region from the 3' LTR. This region was found to be highly variable among pea Ogre elements, and includes a number of other tandem repeats along with or instead of PisTR-A. Bioinformatic analysis of LTR-retrotransposons mined from available plant genomic sequence data revealed that the frequent occurrence of variable tandem repeats within 3' UTRs is a typical feature of the Tat lineage of plant retrotransposons. Comparison of these repeats to known plant satellite sequences uncovered two other instances of satellites with sequence similarity to a Tat-like retrotransposon 3' UTR regions. These observations suggest that some retrotransposons may significantly contribute to satellite DNA evolution by generating a library of short repeat arrays that can subsequently be dispersed through the genome and eventually further amplified and homogenized into novel satellite repeats.

Dynamic rearrangements determine genome organization and useful traits in soybean

Soybean (Glycine max) is a paleopolyploid whose genome has gone through at least two rounds of polyploidy and subsequent diploidization events. Several studies have investigated the changes in genome structure produced by the relatively recent polyploidy event, but little is known about the ancient polyploidy due to the high frequency of gene loss after duplication. Our previous study, regarding a region responsible for bacterial leaf pustule, reported two homeologous Rxp regions produced by the recent whole-genome duplication event. In this study, we identified the full set of four homeologous Rxp regions (ranging from 1.96 to 4.60 Mb) derived from both the recent and ancient polyploidy events, and this supports the quadruplicated structure of the soybean genome. Among the predicted genes on chromosome 17 (linkage group D2), 71% of them were conserved in a recently duplicated region, while 21% and 24% of duplicated genes were retained in two homeologous regions formed by the ancient polyploidy. Furthermore, comparative analysis showed a 2:1 relationship between soybean and Medicago truncatula, since M. truncatula did not undergo the recent polyploidy event that soybean did. Unlike soybean, M. truncatula homeologous regions were highly fractionated and their synteny did not exist, revealing different rates of diploidization process between the two species. Our data show that extensive synteny remained in the four homeologous regions in soybean, even though the soybean genome experienced dynamic genome rearrangements following paleopolyploidy events. Moreover, multiple Rxp quantitative trait loci on different soybean chromosomes actually comprise homeologous regions produced by two rounds of polyploidy events.

Molecular and chromosomal evidence for allopolyploidy in soybean

Recent studies have documented that the soybean (Glycine max) genome has undergone two rounds of large-scale genome and/or segmental duplication. To shed light on the timing and nature of these duplication events, we characterized and analyzed two subfamilies of high-copy centromeric satellite repeats, CentGm-1 and CentGm-2, using a combination of computational and molecular cytogenetic approaches. These two subfamilies of satellite repeats mark distinct subsets of soybean centromeres and, in at least one case, a pair of homologs, suggesting their origins from an allopolyploid event. The satellite monomers of each subfamily are arranged in large tandem arrays, and intermingled monomers of the two subfamilies were not detected by fluorescence in situ hybridization on extended DNA fibers nor at the sequence level. This indicates that there has been little recombination and homogenization of satellite DNA between these two sets of centromeres. These satellite repeats are also present in Glycine soja, the proposed wild progenitor of soybean, but could not be detected in any other relatives of soybean examined in this study, suggesting the rapid divergence of the centromeric satellite DNA within the Glycine genus. Together, these observations provide direct evidence, at molecular and chromosomal levels, in support of the hypothesis that the soybean genome has experienced a recent allopolyploidization event.

Divergence in centromere structure distinguishes related genomes in Coix lacryma-jobi and its wild relative

Knowledge about the composition and structure of centromeres is critical for understanding how centromeres perform their functional roles. Here, we report the sequences of one centromere-associated bacterial artificial chromosome clone from a Coix lacryma-jobi library. Two Ty3/gypsy-class retrotransposons, centromeric retrotransposon of C. lacryma-jobi (CRC) and peri-centromeric retrotransposon of C. lacryma-jobi, and a (peri)centromere-specific tandem repeat with a unit length of 153 bp were identified. The CRC is highly homologous to centromere-specific retrotransposons reported in grass species. An 80-bp DNA region in the 153-bp satellite repeat was found to be conserved to centromeric satellite repeats from maize, rice, and pearl millet. Fluorescence in situ hybridization showed that the three repetitive sequences were located in (peri-)centromeric regions of both C. lacryma-jobi and Coix aquatica. However, the 153-bp satellite repeat was only detected on 20 out of the 30 chromosomes in C. aquatica. Immunostaining with an antibody against rice CENH3 indicates that the 153-bp satellite repeat and CRC might be both the major components for functional centromeres, but not all the 153-bp satellite repeats or CRC sequences are associated with CENH3. The evolution of centromeric repeats of C. lacryma-jobi during the polyploidization was discussed.

Differential accumulation of retroelements and diversification of NB-LRR disease resistance genes in duplicated regions following polyploidy in the ancestor of soybean

The genomes of most, if not all, flowering plants have undergone whole genome duplication events during their evolution. The impact of such polyploidy events is poorly understood, as is the fate of most duplicated genes. We sequenced an approximately 1 million-bp region in soybean (Glycine max) centered on the Rpg1-b disease resistance gene and compared this region with a region duplicated 10 to 14 million years ago. These two regions were also compared with homologous regions in several related legume species (a second soybean genotype, Glycine tomentella, Phaseolus vulgaris, and Medicago truncatula), which enabled us to determine how each of the duplicated regions (homoeologues) in soybean has changed following polyploidy. The biggest change was in retroelement content, with homoeologue 2 having expanded to 3-fold the size of homoeologue 1. Despite this accumulation of retroelements, over 77% of the duplicated low-copy genes have been retained in the same order and appear to be functional. This finding contrasts with recent analyses of the maize (Zea mays) genome, in which only about one-third of duplicated genes appear to have been retained over a similar time period. Fluorescent in situ hybridization revealed that the homoeologue 2 region is located very near a centromere. Thus, pericentromeric localization, per se, does not result in a high rate of gene inactivation, despite greatly accelerated retrotransposon accumulation. In contrast to low-copy genes, nucleotide-binding-leucine-rich repeat disease resistance gene clusters have undergone dramatic species/homoeologue-specific duplications and losses, with some evidence for partitioning of subfamilies between homoeologues.

The soybean-Phytophthora resistance locus Rps1-k encompasses coiled coil-nucleotide binding-leucine rich repeat-like genes and repetitive sequences

BACKGROUND

A series of Rps (resistance to Pytophthora sojae) genes have been protecting soybean from the root and stem rot disease caused by the Oomycete pathogen, Phytophthora sojae. Five Rps genes were mapped to the Rps1 locus located near the 28 cM map position on molecular linkage group N of the composite genetic soybean map. Among these five genes, Rps1-k was introgressed from the cultivar, Kingwa. Rps1-k has been providing stable and broad-spectrum Phytophthora resistance in the major soybean-producing regions of the United States. Rps1-k has been mapped and isolated. More than one functional Rps1-k gene was identified from the Rps1-k locus. The clustering feature at the Rps1-k locus might have facilitated the expansion of Rps1-k gene numbers and the generation of new recognition specificities. The Rps1-k region was sequenced to understand the possible evolutionary steps that shaped the generation of Phytophthora resistance genes in soybean.

RESULTS

Here the analyses of sequences of three overlapping BAC clones containing the 184,111 bp Rps1-k region are reported. A shotgun sequencing strategy was applied in sequencing the BAC contig. Sequence analysis predicted a few full-length genes including two Rps1-k genes, Rps1-k-1 and Rps1-k-2. Previously reported Rps1-k-3 from this genomic region 1 was evolved through intramolecular recombination between Rps1-k-1 and Rps1-k-2 in Escherichia coli. The majority of the predicted genes are truncated and therefore most likely they are nonfunctional. A member of a highly abundant retroelement, SIRE1, was identified from the Rps1-k region. The Rps1-k region is primarily composed of repetitive sequences. Sixteen simple repeat and 63 tandem repeat sequences were identified from the locus.

CONCLUSION

These data indicate that the Rps1 locus is located in a gene-poor region. The abundance of repetitive sequences in the Rps1-k region suggested that the location of this locus is in or near a heterochromatic region. Poor recombination frequencies combined with presence of two functional Rps genes at this locus has been providing stable Phytophthora resistance in soybean.

Gene duplication and paleopolyploidy in soybean and the implications for whole genome sequencing

BACKGROUND

Soybean, Glycine max (L.) Merr., is a well documented paleopolyploid. What remains relatively under characterized is the level of sequence identity in retained homeologous regions of the genome. Recently, the Department of Energy Joint Genome Institute and United States Department of Agriculture jointly announced the sequencing of the soybean genome. One of the initial concerns is to what extent sequence identity in homeologous regions would have on whole genome shotgun sequence assembly.

RESULTS

Seventeen BACs representing approximately 2.03 Mb were sequenced as representative potential homeologous regions from the soybean genome. Genetic mapping of each BAC shows that 11 of the 20 chromosomes are represented. Sequence comparisons between homeologous BACs shows that the soybean genome is a mosaic of retained paleopolyploid regions. Some regions appear to be highly conserved while other regions have diverged significantly. Large-scale "batch" reassembly of all 17 BACs combined showed that even the most homeologous BACs with upwards of 95% sequence identity resolve into their respective homeologous sequences. Potential assembly errors were generated by tandemly duplicated pentatricopeptide repeat containing genes and long simple sequence repeats. Analysis of a whole-genome shotgun assembly of 80,000 randomly chosen JGI-DOE sequence traces reveals some new soybean-specific repeat sequences.

CONCLUSION

This analysis investigated both the structure of the paleopolyploid soybean genome and the potential effects retained homeology will have on assembling the whole genome shotgun sequence. Based upon these results, homeologous regions similar to those characterized here will not cause major assembly issues.

Gene duplication and paleopolyploidy in soybean and the implications for whole genome sequencing

BACKGROUND

Soybean, Glycine max (L.) Merr., is a well documented paleopolyploid. What remains relatively under characterized is the level of sequence identity in retained homeologous regions of the genome. Recently, the Department of Energy Joint Genome Institute and United States Department of Agriculture jointly announced the sequencing of the soybean genome. One of the initial concerns is to what extent sequence identity in homeologous regions would have on whole genome shotgun sequence assembly.

RESULTS

Seventeen BACs representing approximately 2.03 Mb were sequenced as representative potential homeologous regions from the soybean genome. Genetic mapping of each BAC shows that 11 of the 20 chromosomes are represented. Sequence comparisons between homeologous BACs shows that the soybean genome is a mosaic of retained paleopolyploid regions. Some regions appear to be highly conserved while other regions have diverged significantly. Large-scale "batch" reassembly of all 17 BACs combined showed that even the most homeologous BACs with upwards of 95% sequence identity resolve into their respective homeologous sequences. Potential assembly errors were generated by tandemly duplicated pentatricopeptide repeat containing genes and long simple sequence repeats. Analysis of a whole-genome shotgun assembly of 80,000 randomly chosen JGI-DOE sequence traces reveals some new soybean-specific repeat sequences.

CONCLUSION

This analysis investigated both the structure of the paleopolyploid soybean genome and the potential effects retained homeology will have on assembling the whole genome shotgun sequence. Based upon these results, homeologous regions similar to those characterized here will not cause major assembly issues.

Reinterpreting pericentromeric heterochromatin

In fission yeast, pericentromeric heterochromatin is directly responsible for the sister chromatid cohesion that assures accurate chromosome segregation. In plants, however, heterochromatin and chromosome segregation appear to be largely unrelated: chromosome transmission is impaired by mutations in cohesion but not by mutations that affect heterochromatin formation. We argue that the formation of pericentromeric heterochromatin is primarily a response to constraints on chromosome mechanics that disfavor the transmission of recombination events in pericentromeric regions. This effect allows pericentromeres to expand to enormous sizes by the accumulation of transposons and through large-scale insertions and inversions. Although sister chromatid cohesion is spatially limited to pericentromeric regions at mitosis and meiosis II, the cohesive domains appear to be defined independently of heterochromatin. The available data from plants suggest that sister chromatid cohesion is marked by histone phosphorylation and mediated by Aurora kinases.

Sequence conservation of homeologous bacterial artificial chromosomes and transcription of homeologous genes in soybean (Glycine max L. Merr.)

The paleopolyploid soybean genome was investigated by sequencing homeologous BAC clones anchored by duplicate N-hydroxycinnamoyl/benzoyltransferase (HCBT) genes. The homeologous BACs were genetically mapped to linkage groups C1 and C2. Annotation of the 173,747- and 98,760-bp BACs showed that gene conservation in both order and orientation is high between homeologous regions with only a single gene insertion/deletion and local tandem duplications differing between the regions. The nucleotide sequence conservation extends into intergenic regions as well, probably due to conserved regulatory sequences. Most of the homeologs appear to have a role in either transcription/DNA binding or cellular signaling, suggesting a potential preference for retention of duplicate genes with these functions. Reverse transcriptase-PCR analysis of homeologs showed that in the tissues sampled, most homeologs have not diverged greatly in their transcription profiles. However, four cases of changes in transcription were identified, primarily in the HCBT gene cluster. Because a mapped locus corresponds to a soybean cyst nematode (SCN) QTL, the potential role of HCBT genes in response to SCN is discussed. These results are the first sequenced-based analysis of homeologous BACs in soybean, a diploidized paleopolyploid.

Tuareg, a novel miniature-inverted repeat family of pearl millet (Pennisetum glaucum) related to the PIF superfamily of maize

Miniature-inverted repeat transposable elements (MITEs) are abundantly repeated in plant genomes and are especially found in genic regions where they could contribute regulatory elements for gene expression. We describe with molecular and cytological tools the first MITE family reported in pearl millet: Tuareg. It was initially detected in the pearl millet ortholog of Teosinte-branched1, an important developmental gene involved in the domestication of maize. The Tuareg family was amplified recently in the pearl millet genome and elements were found more abundant in wild than in domesticated plants. We found that they shared similarity in their terminal repeats with the previously described mPIF MITEs and that they are also present in other Pennisetum species, in maize and more distantly related grasses. The Tuareg family may be part of MITEs activated by PIF-like transposases and it could have been mobile since pearl millet domestication.

Current status and the future of fluorescence in situ hybridization (FISH) in plant genome research

Fluorescence in situ hybridization (FISH), which allows direct mapping of DNA sequences on chromosomes, has become the most important technique in plant molecular cytogenetics research. Repetitive DNA sequence can generate unique FISH patterns on individual chromosomes for karyotyping and phylogenetic analysis. FISH on meiotic pachytene chromosomes coupled with digital imaging systems has become an efficient method to develop physical maps in plant species. FISH on extended DNA fibers provides a high-resolution mapping approach to analyze large DNA molecules and to characterize large genomic loci. FISH-based physical mapping provides a valuable complementary approach in genome sequencing and map-based cloning research. We expect that FISH will continue to play an important role in relating DNA sequence information to chromosome biology. FISH coupled with immunoassays will be increasingly used to study features of chromatin at the cytological level that control expression and regulation of genes.

Survey sequencing of soybean elucidates the genome structure, composition and identifies novel repeats

In order to expand our knowledge of the soybean genome and to create a useful DNA repeat sequence database, over 24 000 DNA fragments from a soybean [Glycine max (L.) Merr.] cv. Williams 82 genomic shotgun library were sequenced. Additional sequences came from over 29 000 bacterial artificial chromosome (BAC) end sequences derived from a BstI library of the cv. Williams 82 genome. Analysis of these sequences identified 348 different DNA repeats, many of which appear to be novel. To extend the utility of the work, a pilot study was also conducted using methylation filtration to estimate the hypomethylated, soybean gene space. A comparison between 8366 sequences obtained from a filtered library and 23 788 from an unfiltered library indicate a gene-enrichment of ~3.2-fold in the hypomethylated sequences. Given the 1.1-Gb soybean genome, our analysis predicts a ~343-Mb hypomethylated, gene-rich space.

Chromosome-level homeology in paleopolyploid soybean (Glycine max) revealed through integration of genetic and chromosome maps

Soybean has 20 chromosome pairs that are derived from at least two rounds of genomewide duplication or polyploidy events although, cytogenetically, soybean behaves like a diploid and has disomic inheritance for most loci. Genetically anchored genomic clones were used as probes for fluorescence in situ hybridization (FISH) to determine the level of postpolyploid chromosomal rearrangements and to integrate the genetic and physical maps to (1) assign linkage groups to specific chromosomes, (2) assess chromosomal structure, and (3) determine the distribution of recombination along the length of a chromosome. FISH mapping of seven putatively gene-rich BACs from linkage group L (chromosome 19) revealed that most of the genetic map correlates to the highly euchromatic long arm and that there is extensive homeology with another chromosome pair, although colinearity of some loci does appear to be disrupted. Moreover, mapping of BACs containing high-copy sequences revealed sequestration of high-copy repeats to the pericentromeric regions of this chromosome. Taken together, these data present a model of chromosome structure in a highly duplicated but diploidized eukaryote, soybean.