Colinearity and its exceptions in orthologous adh regions of maize and sorghum

Development of Genomic Simple Sequence Repeat Markers for Evaluating Resources of Armillaira ostoyae and Their Transferability to Armillaira gallica

In this study, we aimed to develop simple sequence repeat (SSR) markers for evaluating resources in Armillaria ostoyae and examine their transferability to Armillaria gallica, related species. SSR markers were developed using the released A. ostoyae whole-genome sequence (GenBank assembly accession: GCA_900157425.1). The SSR regions were analyzed using the MISA (MIcroSAtellite identification tool) program. A total of 2319 SSR loci consisting of 922 (39.76%) mononucleotide, 763 (32.90%) trinucleotide, and 517 (22.29%) dinucleotide motifs were identified. Marker design involved an arbitrary choice of 150 SSR loci, considering motif abundance. A total of 22 strains of A. ostoyae were analyzed using the developed markers, and 105 markers were successfully amplified. The mean values of major allele frequency, number of alleles, expected heterozygosity, observed heterozygosity, and polymorphism information content (PIC) values were approximately 5.89, 5.4, 0.541, 0.255, and 0.504, respectively. A. gallica was analyzed, and 52 markers (49.5%) were successfully amplified to evaluate the transferability of the developed SSR markers. When these markers were used, the mean values of major allele frequency, number of alleles, expected heterozygosity, observed heterozygosity, and PIC were calculated to be approximately 0.615, 4.3, 0.517, 0.133, and 0.502, respectively. In conclusion, SSR markers were developed using the genome of A. ostoyae, and some of these markers exhibited transferability to A. gallica. These results can be used for resource evaluation of A. ostoyae and A. gallica.

cis-Regulatory Elements in Plant Development, Adaptation, and Evolution

cis-Regulatory elements encode the genomic blueprints that ensure the proper spatiotemporal patterning of gene expression necessary for appropriate development and responses to the environment. Accumulating evidence implicates changes to gene expression as a major source of phenotypic novelty in eukaryotes, including acute phenotypes such as disease and cancer in mammals. Moreover, genetic and epigenetic variation affecting cis-regulatory sequences over longer evolutionary timescales has become a recurring theme in studies of morphological divergence and local adaptation. Here, we discuss the functions of and methods used to identify various classes of cis-regulatory elements, as well as their role in plant development and response to the environment. We highlight opportunities to exploit cis-regulatory variants underlying plant development and environmental responses for crop improvement efforts. Although a comprehensive understanding of cis-regulatory mechanisms in plants has lagged behind that in animals, we showcase several breakthrough findings that have profoundly influenced plant biology and shaped the overall understanding of transcriptional regulation in eukaryotes.

Dissection of Structural Reorganization of Wheat 5B Chromosome Associated With Interspecies Recombination Suppression

Chromosomal rearrangements that lead to recombination suppression can have a significant impact on speciation, and they are also important for breeding. The regions of recombination suppression in wheat chromosome 5B were identified based on comparisons of the 5B map of a cross between the Chinese Spring (CS) variety of hexaploid wheat and CS-5Bdic (genotype CS with 5B substituted with its homologue from tetraploid Triticum dicoccoides) with several 5B maps of tetraploid and hexaploid wheat. In total, two regions were selected in which recombination suppression occurred in cross CS × CS-5Bdic when compared with other maps: one on the short arm, 5BS_RS, limited by markers BS00009810/BS00022336, and the second on the long arm, 5BL_RS, between markers Ra_c10633_2155 and BS00087043. The regions marked as 5BS_RS and 5BL_RS, with lengths of 5 Mb and 3.6 Mb, respectively, were mined from the 5B pseudomolecule of CS and compared to the homoeologous regions (7.6 and 3.8 Mb, respectively) of the 5B pseudomolecule of Zavitan (T. dicoccoides). It was shown that, in the case of 5BS_RS, the local heterochromatin islands determined by the satellite DNA (119.2) and transposable element arrays, as well as the dissimilarity caused by large insertions/deletions (chromosome rearrangements) between 5BSs aestivum/dicoccoides, are likely the key determinants of recombination suppression in the region. Two major and two minor segments with significant loss of similarity were recognized within the 5BL_RS region. It was shown that the loss of similarity, which can lead to suppression of recombination in the 5BL_RS region, is caused by chromosomal rearrangements, driven by the activity of mobile genetic elements (both DNA transposons and long terminal repeat retrotransposons) and their divergence during evolution. It was noted that the regions marked as 5BS_RS and 5BL_RS are associated with chromosomal rearrangements identified earlier by С-banding analysis of intraspecific polymorphism of tetraploid emmer wheat. The revealed divergence in 5BS_RS and 5BL_RS may be a consequence of interspecific hybridization, plant genetic adaptation, or both.

The Dynamism of Transposon Methylation for Plant Development and Stress Adaptation

Plant development processes are regulated by epigenetic alterations that shape nuclear structure, gene expression, and phenotypic plasticity; these alterations can provide the plant with protection from environmental stresses. During plant growth and development, these processes play a significant role in regulating gene expression to remodel chromatin structure. These epigenetic alterations are mainly regulated by transposable elements (TEs) whose abundance in plant genomes results in their interaction with genomes. Thus, TEs are the main source of epigenetic changes and form a substantial part of the plant genome. Furthermore, TEs can be activated under stress conditions, and activated elements cause mutagenic effects and substantial genetic variability. This introduces novel gene functions and structural variation in the insertion sites and primarily contributes to epigenetic modifications. Altogether, these modifications indirectly or directly provide the ability to withstand environmental stresses. In recent years, many studies have shown that TE methylation plays a major role in the evolution of the plant genome through epigenetic process that regulate gene imprinting, thereby upholding genome stability. The induced genetic rearrangements and insertions of mobile genetic elements in regions of active euchromatin contribute to genome alteration, leading to genomic stress. These TE-mediated epigenetic modifications lead to phenotypic diversity, genetic variation, and environmental stress tolerance. Thus, TE methylation is essential for plant evolution and stress adaptation, and TEs hold a relevant military position in the plant genome. High-throughput techniques have greatly advanced the understanding of TE-mediated gene expression and its associations with genome methylation and suggest that controlled mobilization of TEs could be used for crop breeding. However, development application in this area has been limited, and an integrated view of TE function and subsequent processes is lacking. In this review, we explore the enormous diversity and likely functions of the TE repertoire in adaptive evolution and discuss some recent examples of how TEs impact gene expression in plant development and stress adaptation.

DCT4-A New Member of the Dicarboxylate Transporter Family in C4 Grasses

Malate transport shuttles atmospheric carbon into the Calvin–Benson cycle during NADP-ME C₄ photosynthesis. Previous characterizations of several plant dicarboxylate transporters (DCT) showed that they efficiently exchange malate across membranes. Here, we identify and characterize a previously unknown member of the DCT family, DCT4, in Sorghum bicolor. We show that SbDCT4 exchanges malate across membranes and its expression pattern is consistent with a role in malate transport during C₄ photosynthesis. SbDCT4 is not syntenic to the characterized photosynthetic gene ZmDCT2, and an ortholog is not detectable in the maize reference genome. We found that the expression patterns of DCT family genes in the leaves of Zea mays, and S. bicolor varied by cell type. Our results suggest that subfunctionalization, of members of the DCT family, for the transport of malate into the bundle sheath plastids, occurred during the process of independent recurrent evolution of C₄ photosynthesis in grasses of the PACMAD clade. We also show that this subfunctionalization is lineage independent. Our results challenge the dogma that key C₄ genes must be orthologues of one another among C₄ species, and shed new light on the evolution of C₄ photosynthesis.

Small RNA-Mediated De Novo Silencing of Ac/Ds Transposons Is Initiated by Alternative Transposition in Maize

Although transposable elements (TEs) comprise a major fraction of many higher eukaryotic genomes, most TEs are silenced by host defense mechanisms. The means by which otherwise active TEs are recognized and silenced remains poorly understood. Here we analyzed two independent cases of spontaneous silencing of the active maize Ac/Ds transposon system. This silencing is initiated by alternative transposition, a type of aberrant transposition event that engages the termini of two nearby separate TEs. Alternative transposition during DNA replication can generate Composite Insertions that contain inverted duplications of the transposon sequences. We show that the inverted duplications of two Composite Insertions are transcribed to produce double-stranded RNAs that trigger the production of two distinct classes of small interfering RNAs: a 24-nt class complementary to the TE terminal inverted repeats and noncoding subterminal regions, and a 21- to 22-nt class corresponding to the TE transcribed regions. Plants containing these small interfering RNA-generating Composite Insertions exhibit decreased levels of Ac transcript and heritable repression of Ac/Ds transposition. Further, we demonstrate that Composite Insertions can heritably silence otherwise active elements in trans This study documents the first case of transposon silencing induced by alternative transposition and may represent a general initiating mechanism for silencing of DNA transposons.

Comparative Genomics, Whole-Genome Re-sequencing and Expression Profile Analysis of Nucleobase:Cation Symporter 2 (NCS2) Genes in Maize

Nucleobase:cation symporter 2 (NCS2) proteins are important for the transport of free nucleobases, participating in diverse plant growth and developmental processes, as well as response to abiotic stress. To date, a comprehensive analysis of the NCS2 gene family has not been performed in maize. In this study, we conducted a comparative genomics analysis of NCS2 genes in 28 plant species, ranging from aquatic algae to land plants, concentrating mainly on maize. Gene duplication events contributed to the expansion of NCS2 genes from lower aquatic plants to higher angiosperms, and whole-genome/segmental and single-gene duplication events were responsible for the expansion of the maize NCS2 gene family. Phylogenetic construction showed three NCS2 subfamilies, I, II, and III. According to homology-based relationships, members of subfamily I are NCS2/AzgA-like genes, whereas those in subfamilies II and III are NCS2/NATs. Moreover, subfamily I exhibited ancient origins. A motif compositional analysis showed that one symbolic motif (motif 4) of the NCS2/NAT genes was absent in subfamily I. In maize, three NCS2/AzgA-like and 21 NCS2/NAT genes were identified, and purifying selection influenced the duplication of maize NCS2 genes. Additionally, a population genetic analysis of NCS2 genes revealed that ZmNCS2-21 showed the greatest diversity between the 78 inbred and 22 wild surveyed maize populations. An expression profile analysis using transcriptome data and quantitative real-time PCR revealed that NCS2 genes in maize are involved in diverse developmental processes and responses to abiotic stresses, including abscisic acid, salt (NaCl), polyethylene glycol, and low (4°C) and high (42°C) temperatures. ZmNCS2 genes with relatively close relationships had similar expression patterns, strongly indicating functional redundancy. Finally, ZmNCS2-16 and ZmNCS2-23 localize in the plasma membrane, which confirmed their predicted membrane structures. These results provide a foundation for future studies regarding the functions of ZmNCS2 proteins, particularly those with potentially important roles in plant responses to abiotic stresses.

A comparative gene analysis with rice identified orthologous group II HKT genes and their association with Na(+) concentration in bread wheat

BACKGROUND

Although the HKT transporter genes ascertain some of the key determinants of crop salt tolerance mechanisms, the diversity and functional role of group II HKT genes are not clearly understood in bread wheat. The advanced knowledge on rice HKT and whole genome sequence was, therefore, used in comparative gene analysis to identify orthologous wheat group II HKT genes and their role in trait variation under different saline environments.

RESULTS

The four group II HKTs in rice identified two orthologous gene families from bread wheat, including the known TaHKT2;1 gene family and a new distinctly different gene family designated as TaHKT2;2. A single copy of TaHKT2;2 was found on each homeologous chromosome arm 7AL, 7BL and 7DL and each gene was expressed in leaf blade, sheath and root tissues under non-stressed and at 200 mM salt stressed conditions. The proteins encoded by genes of the TaHKT2;2 family revealed more than 93% amino acid sequence identity but ≤52% amino acid identity compared to the proteins encoded by TaHKT2;1 family. Specifically, variations in known critical domains predicted functional differences between the two protein families. Similar to orthologous rice genes on chromosome 6L, TaHKT2;1 and TaHKT2;2 genes were located approximately 3 kb apart on wheat chromosomes 7AL, 7BL and 7DL, forming a static syntenic block in the two species. The chromosomal region on 7AL containing TaHKT2;1 7AL-1 co-located with QTL for shoot Na(+) concentration and yield in some saline environments.

CONCLUSION

The differences in copy number, genes sequences and encoded proteins between TaHKT2;2 homeologous genes and other group II HKT gene families within and across species likely reflect functional diversity for ion selectivity and transport in plants. Evidence indicated that neither TaHKT2;2 nor TaHKT2;1 were associated with primary root Na(+) uptake but TaHKT2;1 may be associated with trait variation for Na(+) exclusion and yield in some but not all saline environments.

Alternative Transposition Generates New Chimeric Genes and Segmental Duplications at the Maize p1 Locus

The maize Ac/Ds transposon family was the first transposable element system identified and characterized by Barbara McClintock. Ac/Ds transposons belong to the hAT family of class II DNA transposons. We and others have shown that Ac/Ds elements can undergo a process of alternative transposition in which the Ac/Ds transposase acts on the termini of two separate, nearby transposons. Because these termini are present in different elements, alternative transposition can generate a variety of genome alterations such as inversions, duplications, deletions, and translocations. Moreover, Ac/Ds elements transpose preferentially into genic regions, suggesting that structural changes arising from alternative transposition may potentially generate chimeric genes at the rearrangement breakpoints. Here we identified and characterized 11 independent cases of gene fusion induced by Ac alternative transposition. In each case, a functional chimeric gene was created by fusion of two linked, paralogous genes; moreover, each event was associated with duplication of the ∼70-kb segment located between the two paralogs. An extant gene in the maize B73 genome that contains an internal duplication apparently generated by an alternative transposition event was also identified. Our study demonstrates that alternative transposition-induced duplications may be a source for spontaneous creation of diverse genome structures and novel genes in maize.

Genetic diversity of maize germplasm assessed by retrotransposon-based markers

Maize is one of the most important crops and also a model for grass genome research. Transposable elements comprise over 78% of the maize genome and their ability to generate new copies makes them good potential markers. Interretrotransposon-amplified polymorphism (IRAP) and retrotransposon microsatellite amplified polymorphism (REMAP) protocols were used for the first time in maize to study the genetic variability between maize cultivars. Ten PCR primers were selected based on a systematic analysis of the sequence conservation in the extremities of different high copy number transposable elements, whereas one primer was chosen based on a microsatellite sequence. Of the 16 primer combinations tested, 14 produced polymorphic bands. These markers were used to identify genetic similarity among 20 maize cultivars selected by their different kernel oil content. Genetic similarity analysis was performed based on the polymorphic band profiles and dendrograms were developed by the unweighted pair-group method with arithmetic averages. Clustering technique revealed that samples were grouped into three clusters that differed in their kernel oil content and size, and in their relative embryo size. In the current investigation, there is evidence that IRAP/REMAP may be useful as markers in maize.

Genomic resources for gene discovery, functional genome annotation, and evolutionary studies of maize and its close relatives

Maize is one of the most important food crops and a key model for genetics and developmental biology. A genetically anchored and high-quality draft genome sequence of maize inbred B73 has been obtained to serve as a reference sequence. To facilitate evolutionary studies in maize and its close relatives, much like the Oryza Map Alignment Project (OMAP) (www.OMAP.org) bacterial artificial chromosome (BAC) resource did for the rice community, we constructed BAC libraries for maize inbred lines Zheng58, Chang7-2, and Mo17 and maize wild relatives Zea mays ssp. parviglumis and Tripsacum dactyloides. Furthermore, to extend functional genomic studies to maize and sorghum, we also constructed binary BAC (BIBAC) libraries for the maize inbred B73 and the sorghum landrace Nengsi-1. The BAC/BIBAC vectors facilitate transfer of large intact DNA inserts from BAC clones to the BIBAC vector and functional complementation of large DNA fragments. These seven Zea Map Alignment Project (ZMAP) BAC/BIBAC libraries have average insert sizes ranging from 92 to 148 kb, organellar DNA from 0.17 to 2.3%, empty vector rates between 0.35 and 5.56%, and genome equivalents of 4.7- to 8.4-fold. The usefulness of the Parviglumis and Tripsacum BAC libraries was demonstrated by mapping clones to the reference genome. Novel genes and alleles present in these ZMAP libraries can now be used for functional complementation studies and positional or homology-based cloning of genes for translational genomics.

A high density physical map of chromosome 1BL supports evolutionary studies, map-based cloning and sequencing in wheat

BACKGROUND

As for other major crops, achieving a complete wheat genome sequence is essential for the application of genomics to breeding new and improved varieties. To overcome the complexities of the large, highly repetitive and hexaploid wheat genome, the International Wheat Genome Sequencing Consortium established a chromosome-based strategy that was validated by the construction of the physical map of chromosome 3B. Here, we present improved strategies for the construction of highly integrated and ordered wheat physical maps, using chromosome 1BL as a template, and illustrate their potential for evolutionary studies and map-based cloning.

RESULTS

Using a combination of novel high throughput marker assays and an assembly program, we developed a high quality physical map representing 93% of wheat chromosome 1BL, anchored and ordered with 5,489 markers including 1,161 genes. Analysis of the gene space organization and evolution revealed that gene distribution and conservation along the chromosome results from the superimposition of the ancestral grass and recent wheat evolutionary patterns, leading to a peak of synteny in the central part of the chromosome arm and an increased density of non-collinear genes towards the telomere. With a density of about 11 markers per Mb, the 1BL physical map provides 916 markers, including 193 genes, for fine mapping the 40 QTLs mapped on this chromosome.

CONCLUSIONS

Here, we demonstrate that high marker density physical maps can be developed in complex genomes such as wheat to accelerate map-based cloning, gain new insights into genome evolution, and provide a foundation for reference sequencing.

New insights into nested long terminal repeat retrotransposons in Brassica species

Long terminal repeat (LTR) retrotransposons, one of the foremost types of transposons, continually change or modify gene function and reorganize the genome through bursts of dramatic proliferation. Many LTR-TEs preferentially insert within other LTR-TEs, but the cause and evolutionary significance of these nested LTR-TEs are not well understood. In this study, a total of 1.52Gb of Brassica sequence containing 2020 bacterial artificial chromosomes (BACs) was scanned, and six bacterial artificial chromosome (BAC) clones with extremely nested LTR-TEs (LTR-TEs density: 7.24/kb) were selected for further analysis. The majority of the LTR-TEs in four of the six BACs were found to be derived from the rapid proliferation of retrotransposons originating within the BAC regions, with only a few LTR-TEs originating from the proliferation and insertion of retrotransposons from outside the BAC regions approximately 5-23Mya. LTR-TEs also preferably inserted into TA-rich repeat regions. Gene prediction by Genescan identified 207 genes in the 0.84Mb of total BAC sequences. Only a few genes (3/207) could be matched to the Brassica expressed sequence tag (EST) database, indicating that most genes were inactive after retrotransposon insertion. Five of the six BACs were putatively centromeric. Hence, nested LTR-TEs in centromere regions are rapidly duplicated, repeatedly inserted, and act to suppress activity of genes and to reshuffle the structure of the centromeric sequences. Our results suggest that LTR-TEs burst and proliferate on a local scale to create nested LTR-TE regions, and that these nested LTR-TEs play a role in the formation of centromeres.

Characterization of transposable elements in the ectomycorrhizal fungus Laccaria bicolor

Background

The publicly available Laccaria bicolor genome sequence has provided a considerable genomic resource allowing systematic identification of transposable elements (TEs) in this symbiotic ectomycorrhizal fungus. Using a TE-specific annotation pipeline we have characterized and analyzed TEs in the L. bicolor S238N-H82 genome.

Methodology/Principal Findings

TEs occupy 24% of the 60 Mb L. bicolor genome and represent 25,787 full-length and partial copy elements distributed within 171 families. The most abundant elements were the Copia-like. TEs are not randomly distributed across the genome, but are tightly nested or clustered. The majority of TEs exhibits signs of ancient transposition except some intact copies of terminal inverted repeats (TIRS), long terminal repeats (LTRs) and a large retrotransposon derivative (LARD) element. There were three main periods of TE expansion in L. bicolor: the first from 57 to 10 Mya, the second from 5 to 1 Mya and the most recent from 0.5 Mya ago until now. LTR retrotransposons are closely related to retrotransposons found in another basidiomycete, Coprinopsis cinerea.

Conclusions

This analysis 1) represents an initial characterization of TEs in the L. bicolor genome, 2) contributes to improve genome annotation and a greater understanding of the role TEs played in genome organization and evolution and 3) provides a valuable resource for future research on the genome evolution within the Laccaria genus.

Revolver and superior: novel transposon-like gene families of the plant kingdom

High-throughput sequencing of eukaryotic genomes has revived interest in the structure and function of repetitive genomic sequences, previously referred to as junk DNA. Repetitive sequences, including transposable elements, are now believed to play a significant role in genomic differentiation and evolution. Some are also expressed as regulatory noncoding RNAs. Vast DNA databases exist for higher eukaryotes; however, with the exception of homologues of known repetitive-sequence-families and transposable elements, most repetitive elements still need to be annotated. Revolver and Superior, both discovered in the Triticeae, are novel classes of transposon-like genes and major components of large cereal genomes. Revolver was isolated from rye via genome subtraction of sequences common to rye and wheat. Superior was isolated from rye by cleavage with EcoO109I, the recognition sites of which consist of a 5'- PuGGNCCPy-3' multi-sequence. Revolver is 2929-3041 bp long with an inverted repeat sequence on each end. The Superior family elements are 1292-1432 bp in length, with divergent 5' regions, indicating the presence of considerable structural diversity. Revolver and Superior are transcriptionally active elements; Revolver harbors a single gene consisting of three exons and two introns, encoding a protein of 139 amino acid residues. Revolver variants range in size from 2665 bp to 4269 bp, with some variants lacking the 5' region, indicating structural diversity around the first exon. Revolver and Superior are dispersed across all seven chromosomes of rye. Revolver has existed since the diploid progenitor of wheat, and has been amplified or lost in several species during the evolution of the Triticeae. This article reviews the recently discovered Revolver and Superior families of plant transposons, which do not share identity with any known autonomous transposable elements or repetitive elements from any living species.

Collinearity-based marker mining for the fine mapping of Pm6, a powdery mildew resistance gene in wheat

The genome sequences of rice (Oryza sativa L.) and Brachypodium distachyon and the comprehensive Triticeae EST (Expressed Sequence Tag) resources provide invaluable information for comparative genomics analysis. The powdery mildew resistance gene, Pm6, which was introgressed into common wheat from Triticum timopheevii, was previously mapped to the wheat chromosome bin of 2BL [fraction length (FL) 0.50-1.00] with limited DNA markers. In this study, we saturated the Pm6 locus in wheat using the collinearity-based markers by extensively exploiting these genomic resources. All wheat ESTs located in the bin 2BL FL 0.50-1.00 and their corresponding orthologous genes on rice chromosome 4 were firstly used to develop STS (Sequence Tagged Site) markers. Those identified markers that flanked the Pm6 locus were then used to identify the collinear regions in the genomes of rice and Brachypodium. Triticeae ESTs with orthologous genes in these collinear regions were further used to develop new conserved markers for the fine mapping of Pm6. Using two F(2) populations derived from crosses of IGVI-465 × Prins and IGVI-466 × Prins, we mapped a total of 29 markers to the Pm6 locus. Among them, 14 markers were co-segregated with Pm6 in the IGVI-466/Prins population. Comparative genome analysis showed that the collinear region of the 29 linked markers covers a ~5.6-Mb region in chromosome 5L of Brachypodium and a ~6.0-Mb region in chromosome 4L of rice. The marker order is conserved between rice and Brachypodium, but re-arrangements are present in wheat. Comparative mapping in the two populations showed that two conserved markers (CINAU123 and CINAU127) flanked the Pm6 locus, and an LRR-receptor-like protein kinase cluster was identified in the collinear regions of Brachypodium and rice. This putative resistance gene cluster provides a potential target site for further fine mapping and cloning of Pm6. Moreover, the newly developed conserved markers closely linked to Pm6 can be used for the marker-assisted selection (MAS) of Pm6 in wheat breeding programs.

Genome size and transposable element content as determined by high-throughput sequencing in maize and Zea luxurians

The genome of maize (Zea mays ssp. mays) consists mostly of transposable elements (TEs) and varies in size among lines. This variation extends to other species in the genus Zea: although maize and Zea luxurians diverged only ∼140,000 years ago, their genomes differ in size by ∼50%. We used paired-end Illumina sequencing to evaluate the potential contribution of TEs to the genome size difference between these two species. We aligned the reads both to a filtered gene set and to an exemplar database of unique repeats representing 1,514 TE families; ∼85% of reads mapped against TE repeats in both species. The relative contribution of TE families to the B73 genome was highly correlated with previous estimates, suggesting that reliable estimates of TE content can be obtained from short high-throughput sequencing reads, even at low coverage. Because we used paired-end reads, we could assess whether a TE was near a gene by determining if one paired read mapped to a TE and the second read mapped to a gene. Using this method, Class 2 DNA elements were found significantly more often in genic regions than Class 1 RNA elements, but Class 1 elements were found more often near other TEs. Overall, we found that both Class 1 and 2 TE families account for ∼70% of the genome size difference between B73 and luxurians. Interestingly, the relative abundance of TE families was conserved between species (r = 0.97), suggesting genome-wide control of TE content rather than family-specific effects.

The history and disposition of transposable elements in polyploid Gossypium

Transposable elements (TEs) are a major component of plant genomes. It is of particular interest to explore the potential activation of TE proliferation, especially in hybrids and polyploids, which often are associated with rapid genomic and epigenetic restructuring. Here we explore the consequences of genomic merger and doubling on copia and gypsy-like Gorge3 long terminal repeat (LTR) retrotransposons as well as on non-LTR long interspersed nuclear elements (LINEs) in allotetraploid cotton, Gossypium hirsutum. Using phylogenetic and quantitative methods, we describe the composition and genomic origin of TEs in polyploid Gossypium. In addition, we present information on ancient and recent transposition activities of the three TE types and demonstrate the absence of an impressive proliferation of TEs following polyploidization in Gossypium. Further, we provide evidence for present-day transcription of LINEs, a relatively minor component of Gossypium genomes, whereas the more abundant LTR retrotransposons display limited expression and only under stressed conditions.

Rice structural variation: a comparative analysis of structural variation between rice and three of its closest relatives in the genus Oryza

Rapid progress in comparative genomics among the grasses has revealed similar gene content and order despite exceptional differences in chromosome size and number. Large- and small-scale genomic variations are of particular interest, especially among cultivated and wild species, as they encode rapidly evolving features that may be important in adaptation to particular environments. We present a genome-wide study of intermediate-sized structural variation (SV) among rice (Oryza sativa) and three of its closest relatives in the genus Oryza (Oryza nivara, Oryza rufipogon and Oryza glaberrima). We computationally identified regional expansions, contractions and inversions in the Oryza species genomes relative to O. sativa by combining data from paired-end clone alignments to the O. sativa reference genome and physical maps. A subset of the computational predictions was validated using a new approach for BAC size determination. The result was a confirmed catalog of 674 expansions (25-38 Mb) and 611 (4-19 Mb) contractions, and 140 putative inversions (14-19 Mb) between the three Oryza species and O. sativa. In the expanded regions unique to O. sativa we found enrichment in transposable elements (TEs): long terminal repeats (LTRs) were randomly located across the chromosomes, and their insertion times corresponded to the date of the A genome radiation. Also, rice-expanded regions contained an over-representation of single-copy genes related to defense factors in the environment. This catalog of confirmed SV in reference to O. sativa provides an entry point for future research in genome evolution, speciation, domestication and novel gene discovery.

Structural and functional divergence of a 1-Mb duplicated region in the soybean (Glycine max) genome and comparison to an orthologous region from Phaseolus vulgaris

Soybean (Glycine max) has undergone at least two rounds of polyploidization, resulting in a paleopolyploid genome that is a mosaic of homoeologous regions. To determine the structural and functional impact of these duplications, we sequenced two ~1-Mb homoeologous regions of soybean, Gm8 and Gm15, derived from the most recent ~13 million year duplication event and the orthologous region from common bean (Phaseolus vulgaris), Pv5. We observed inversions leading to major structural variation and a bias between the two chromosome segments as Gm15 experienced more gene movement (gene retention rate of 81% in Gm15 versus 91% in Gm8) and a nearly twofold increase in the deletion of long terminal repeat (LTR) retrotransposons via solo LTR formation. Functional analyses of Gm15 and Gm8 revealed decreases in gene expression and synonymous substitution rates for Gm15, for instance, a 38% increase in transcript levels from Gm8 relative to Gm15. Transcriptional divergence of homoeologs was found based on expression patterns among seven tissues and developmental stages. Our results indicate asymmetric evolution between homoeologous regions of soybean as evidenced by structural changes and expression variances of homoeologous genes.

Effective isolation of retrotransposons and repetitive DNA families from the wheat genome

New classes of repetitive DNA elements were effectively identified by isolating small fragments of the elements from the wheat genome. A wheat A genome library was constructed from Triticum monococcum by degenerate cleavage with EcoO109I, the recognition sites of which consisted of 5'-PuGGNCCPy-3' multi-sequences. Three novel repetitive sequences pTm6, pTm69 and pTm58 derived from the A genome were screened and tested for high copy number using a blotting approach. pTm6 showed identity with integrase domains of the barley Ty1-Copia-retrotransposon BARE-1 and pTm58 showed similarity to the barley Ty3-gypsy-like retrotransposon Romani. pTm69, however, constituted a tandem array with useful genomic specificities, but did not share any identity with known repetitive elements. This study also sought to isolate wheat D-genome-specific repetitive elements regardless of the level of methylation, by genomic subtraction. Total genomic DNA of Aegilops tauschii was cleaved into short fragments with a methylation-insensitive 4 bp cutter, MboI, and then common DNA sequences between Ae. tauschii and Triticum turgidum were subtracted by annealing with excess T. turgidum genomic DNA. The D genome repetitive sequence pAt1 was isolated and used to identify an additional novel repetitive sequence family from wheat bacterial artificial chromosomes with a size range of 1 395-1 850 bp. The methods successfully led pathfinding of two unique repetitive families.

Orthologous comparisons of the Hd1 region across genera reveal Hd1 gene lability within diploid Oryza species and disruptions to microsynteny in Sorghum

Heading date is one of the most important quantitative traits responsible for the domestication of rice. We compared a 155-kb reference segment of the Oryza sativa ssp. japonica cv. Nipponbare genome surrounding Hd1, a major heading date gene in rice, with orthologous regions from nine diploid Oryza species that diverged over a relatively short time frame (∼16 My) to study sequence evolution around a domestication locus. The orthologous Hd1 region from Sorghum bicolor was included to compare and contrast the evolution in a more distant relative of rice. Consistent with other observations at the adh1/adh2, monoculm1, and sh2/a1 loci in grass species, we found high gene colinearity in the Hd1 region amidst size differences that were lineage specific and long terminal repeat retrotransposon driven. Unexpectedly, the Hd1 gene was deleted in O. glaberrima, whereas the O. rufipogon and O. punctata copies had degenerative mutations, suggesting that other heading date loci might compensate for the loss or nonfunctionality of Hd1 in these species. Compared with the japonica Hd1 region, the orthologous region in sorghum exhibited micro-rearrangements including gene translocations, seven additional genes, and a gene triplication and truncation event predating the divergence from Oryza.

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.

A highly conserved gene island of three genes on chromosome 3B of hexaploid wheat: diverse gene function and genomic structure maintained in a tightly linked block

BACKGROUND

The complexity of the wheat genome has resulted from waves of retrotransposable element insertions. Gene deletions and disruptions generated by the fast replacement of repetitive elements in wheat have resulted in disruption of colinearity at a micro (sub-megabase) level among the cereals. In view of genomic changes that are possible within a given time span, conservation of genes between species tends to imply an important functional or regional constraint that does not permit a change in genomic structure. The ctg1034 contig completed in this paper was initially studied because it was assigned to the Sr2 resistance locus region, but detailed mapping studies subsequently assigned it to the long arm of 3B and revealed its unusual features.

RESULTS

BAC shotgun sequencing of the hexaploid wheat (Triticum aestivum cv. Chinese Spring) genome has been used to assemble a group of 15 wheat BACs from the chromosome 3B physical map FPC contig ctg1034 into a 783,553 bp genomic sequence. This ctg1034 sequence was annotated for biological features such as genes and transposable elements. A three-gene island was identified among >80% repetitive DNA sequence. Using bioinformatics analysis there were no observable similarity in their gene functions. The ctg1034 gene island also displayed complete conservation of gene order and orientation with syntenic gene islands found in publicly available genome sequences of Brachypodium distachyon, Oryza sativa, Sorghum bicolor and Zea mays, even though the intergenic space and introns were divergent.

CONCLUSION

We propose that ctg1034 is located within the heterochromatic C-band region of deletion bin 3BL7 based on the identification of heterochromatic tandem repeats and presence of significant matches to chromodomain-containing gypsy LTR retrotransposable elements. We also speculate that this location, among other highly repetitive sequences, may account for the relative stability in gene order and orientation within the gene island.Sequence data from this article have been deposited with the GenBank Data Libraries under accession no. GQ422824.

Exceptional diversity, non-random distribution, and rapid evolution of retroelements in the B73 maize genome

Recent comprehensive sequence analysis of the maize genome now permits detailed discovery and description of all transposable elements (TEs) in this complex nuclear environment. Reiteratively optimized structural and homology criteria were used in the computer-assisted search for retroelements, TEs that transpose by reverse transcription of an RNA intermediate, with the final results verified by manual inspection. Retroelements were found to occupy the majority (>75%) of the nuclear genome in maize inbred B73. Unprecedented genetic diversity was discovered in the long terminal repeat (LTR) retrotransposon class of retroelements, with >400 families (>350 newly discovered) contributing >31,000 intact elements. The two other classes of retroelements, SINEs (four families) and LINEs (at least 30 families), were observed to contribute 1,991 and ∼35,000 copies, respectively, or a combined ∼1% of the B73 nuclear genome. With regard to fully intact elements, median copy numbers for all retroelement families in maize was 2 because >250 LTR retrotransposon families contained only one or two intact members that could be detected in the B73 draft sequence. The majority, perhaps all, of the investigated retroelement families exhibited non-random dispersal across the maize genome, with LINEs, SINEs, and many low-copy-number LTR retrotransposons exhibiting a bias for accumulation in gene-rich regions. In contrast, most (but not all) medium- and high-copy-number LTR retrotransposons were found to preferentially accumulate in gene-poor regions like pericentromeric heterochromatin, while a few high-copy-number families exhibited the opposite bias. Regions of the genome with the highest LTR retrotransposon density contained the lowest LTR retrotransposon diversity. These results indicate that the maize genome provides a great number of different niches for the survival and procreation of a great variety of retroelements that have evolved to differentially occupy and exploit this genomic diversity.

Although TEs are a major component of all studied plant genomes, and are the most significant contributors to genome structure and evolution in almost all eukaryotes that have been investigated, their properties and reasons for existence are not well understood in any eukaryotic genome. In order to begin a comprehensive study of TE contributions to the structure, function, and evolution of both genes and genomes, we first identified all of the TEs in maize and then investigated whether there were non-random patterns in their dispersal. We used homology and TE structure criteria in an effort to discover all of the retroelements in the recently sequenced genome from maize inbred B73. We found that the retroelements are incredibly diverse in maize, with many hundreds of families that show different insertion and/or retention specificities across the maize chromosomes. Most of these element families are present in low copy numbers and had been missed by previous searches that relied on a high-copy-number criterion. Different element families exhibited very different biases for accumulation across the chromosomes, indicating that they can detect and utilize many different chromatin environments.

Detailed analysis of a contiguous 22-Mb region of the maize genome

Most of our understanding of plant genome structure and evolution has come from the careful annotation of small (e.g., 100 kb) sequenced genomic regions or from automated annotation of complete genome sequences. Here, we sequenced and carefully annotated a contiguous 22 Mb region of maize chromosome 4 using an improved pseudomolecule for annotation. The sequence segment was comprehensively ordered, oriented, and confirmed using the maize optical map. Nearly 84% of the sequence is composed of transposable elements (TEs) that are mostly nested within each other, of which most families are low-copy. We identified 544 gene models using multiple levels of evidence, as well as five miRNA genes. Gene fragments, many captured by TEs, are prevalent within this region. Elimination of gene redundancy from a tetraploid maize ancestor that originated a few million years ago is responsible in this region for most disruptions of synteny with sorghum and rice. Consistent with other sub-genomic analyses in maize, small RNA mapping showed that many small RNAs match TEs and that most TEs match small RNAs. These results, performed on approximately 1% of the maize genome, demonstrate the feasibility of refining the B73 RefGen_v1 genome assembly by incorporating optical map, high-resolution genetic map, and comparative genomic data sets. Such improvements, along with those of gene and repeat annotation, will serve to promote future functional genomic and phylogenomic research in maize and other grasses.

From dwarves to giants? Plant height manipulation for biomass yield

The increasing demand for lignocellulosic biomass for the production of biofuels provides value to vegetative plant tissue and leads to a paradigm shift for optimizing plant architecture in bioenergy crops. Plant height (PHT) is among the most important biomass yield components and is the focus of this review, with emphasis on the energy grasses maize (Zea mays) and sorghum (Sorghum bicolor). We discuss the scientific advances in the identification of PHT quantitative trait loci (QTLs) and the understanding of pathways and genes controlling PHT, especially gibberellins and brassinosteroids. We consider pleiotropic effects of QTLs or genes affecting PHT on other agronomically important traits and, finally, we discuss strategies for applying this knowledge to the improvement of dual-purpose or dedicated bioenergy crops.

Structural analysis of 83-kb genomic DNA from Thellungiella halophila: sequence features and microcolinearity between salt cress and Arabidopsis thaliana

Salt cress (Thellungiella halophila) has become a desirable plant model for molecular analysis of the mechanisms of salt tolerance. Analysis of its physiological action and expressed EST has resulted in better understanding. However, less is known about its genomic features. Here we determined a continuous sequence approximately 83 kb from a salt cress BAC clone, providing the first insight into the genomic feature for this species. The gene density is approximately one gene per 3.6 kb in this sequence. Many types of repetitive sequences are present in this salt cress sequence, including LTR retroelements, DNA transposons and a number of simple sequence repeats. Comparison of sequence similarity indicated that salt cress shares a close relationship with Arabidopsis. Extensive conservation and high-level microcolinearity were uncovered for both genomes. Our study also indicated that genomic DNA alternations (involving chromosome inversion, sequence loss and gene translocation) contributed to the genomic discrepancies between salt cress and Arabidopsis.

Comparative sequence analysis of MONOCULM1-orthologous regions in 14 Oryza genomes

Comparative genomics is a powerful tool to decipher gene and genome evolution. Placing multiple genome comparisons in a phylogenetic context improves the sensitivity of evolutionary inferences. In the genus Oryza, this comparative approach can be used to investigate gene function, genome evolution, domestication, polyploidy, and ecological adaptation. A large genomic region surrounding the MONOCULM1 (MOC1) locus was chosen for study in 14 Oryza species, including 10 diploids and 4 allotetraploids. Sequencing and annotation of 18 bacterial artificial chromosome clones for these species revealed highly conserved gene colinearity and structure in the MOC1 region. Since the Oryza radiation about 14 Mya, differences in transposon amplification appear to be responsible for the different current sizes of the Oryza genomes. In the MOC1 region, transposons were only conserved between genomes of the same type (e.g., AA or BB). In addition to the conserved gene content, several apparent genes have been generated de novo or uniquely retained in the AA lineage. Two different 3-gene segments have been inserted into the MOC1 region of O. coarctata (KK) or O. sativa by unknown mechanism(s). Large and apparently noncoding sequences flanking the MOC1 gene were observed to be under strong purifying selection. The allotetraploids Oryza alta and Oryza minuta were found to be products of recent polyploidization, less than 1.6 and 0.4 Mya, respectively. In allotetraploids, pseudogenization of duplicated genes was common, caused by large deletions, small frame-shifting insertions/deletions, or nonsense mutations.

Dynamic evolution of oryza genomes is revealed by comparative genomic analysis of a genus-wide vertical data set

Oryza (23 species; 10 genome types) contains the world's most important food crop - rice. Although the rice genome serves as an essential tool for biological research, little is known about the evolution of the other Oryza genome types. They contain a historical record of genomic changes that led to diversification of this genus around the world as well as an untapped reservoir of agriculturally important traits. To investigate the evolution of the collective Oryza genome, we sequenced and compared nine orthologous genomic regions encompassing the Adh1-Adh2 genes (from six diploid genome types) with the rice reference sequence. Our analysis revealed the architectural complexities and dynamic evolution of this region that have occurred over the past approximately 15 million years. Of the 46 intact genes and four pseudogenes in the japonica genome, 38 (76%) fell into eight multigene families. Analysis of the evolutionary history of each family revealed independent and lineage-specific gain and loss of gene family members as frequent causes of synteny disruption. Transposable elements were shown to mediate massive replacement of intergenic space (>95%), gene disruption, and gene/gene fragment movement. Three cases of long-range structural variation (inversions/deletions) spanning several hundred kilobases were identified that contributed significantly to genome diversification.

A single nucleotide polymorphism (SNP839) in the adh1 reference gene affects the quantitation of genetically modified maize (Zea mays L.)

The real-time PCR methods recommended in the European Union for the quantitation of genetically modified (GM) maize events NK603, GA21, and MON 863 measure the number of copies of the GM event in relation to those of the maize-specific adh1 reference gene. The study reported here revealed that the targeted 70 base pair adh1 region exhibits a single nucleotide polymorphism (SNP839) that hampers the binding of the reverse primer used in the adh1 detection method. Partial fragments of the adh1-A and adh1-F allele were cloned. By allele-specific real-time PCR, it was shown that SNP839 corresponds to a common allelic polymorphism in maize. As a result, the quantitation of the GM maize events mentioned is positively or negatively biased, depending on the adh1 genotype of sample and calibrant. Therefore, it is proposed to revise the quantitative detection methods for NK603, GA21, and MON 863 maize.

Molecular characterization of the waxy locus in sorghum

A comparison of approximately 4.5 kb of nucleotide sequence from the waxy locus (the granule-bound starch synthase I [GBSS I] locus) from a waxy line, BTxARG1, and a non-waxy line, QL39, revealed an extremely high level of sequence conservation. Among a total of 24 nucleotide differences and 9 indels, only 2 nucleotide changes resulted in altered amino acid residues. Protein folding prediction software suggested that one of the amino acid changes (Glu to His) may result in an altered protein structure, which may explain the apparently inactive GBSS I present in BTxARG1. This SNP was not found in the second waxy line, RTx2907, which does not produce GBSS I, and no other SNPs or indels were found in the approximately 4 kb of sequence obtained from RTx2907. Using one indel, the waxy locus was mapped to sorghum chromosome SBI-10, which is syntenous to maize chromosome 9; the waxy locus has been mapped to this maize chromosome. The distribution of indels in a diverse set of sorghum germplasm suggested that there are two broad types of non-waxy GBSS I alleles, each type comprising several alleles, and that the two waxy alleles in BTxARG1 and RTx2907 have evolved from one of the non-waxy allele types. The Glu/His polymorphism was found only in BTxARG1 and derived lines and has potential as a perfect marker for the BTxARG1 source of the waxy allele at the GBSS I locus. The indels correctly predicted the non-waxy phenotype in approximately 65% of diverse sorghum germplasm. The indels co-segregated perfectly with phenotype in two sorghum populations derived from crosses between a waxy and a non-waxy sorghum line, correctly identifying heterozygous lines. Thus, these indel markers or sequence-based SNP markers can be used to follow waxy alleles in sorghum breeding programs in selected pedigrees.

Contrasting patterns of transposable-element insertion polymorphism and nucleotide diversity in autotetraploid and allotetraploid Arabidopsis species

It has been hypothesized that polyploidy permits the proliferation of transposable elements, due to both the masking of deleterious recessive mutations and the breakdown of host silencing mechanisms. We investigated the patterns of insertion polymorphism of an Ac-like transposable element and nucleotide diversity at 18 gene fragments in the allotetraploid Arabidopsis suecica and the autotetraploid A. arenosa. All identified insertions were fixed in A. suecica, and many were clearly inherited from the parental species A. thaliana or A. arenosa. These results are inconsistent with a rapid increase in transposition associated with hybrid breakdown but support the evidence from nucleotide polymorphism patterns of a recent single origin of this species leading to genomewide fixations of transposable elements. In contrast, most insertions were segregating at very low frequencies in A. arenosa samples, showing a significant departure from neutrality in favor of purifying selection, even when we account for population subdivision inferred from sequence variation. Patterns of nucleotide variation at reference genes are consistent with the TE results, showing evidence for higher effective population sizes in A. arenosa than in related diploid taxa but a near complete population bottleneck associated with the origins of A. suecica.

Use of transferable Nicotiana tabacum L. microsatellite markers for investigating genetic diversity in the genus Nicotiana

The recent development of microsatellite markers for tobacco, Nicotiana tabacum L., may be valuable for genetic studies within the genus Nicotiana. The first objective was to evaluate transferability of 100 N. tabacum microsatellite primer combinations to 5 diploid species closely related to tobacco. The number of primer combinations that amplified scorable bands in these species ranged from 42 to 56. Additional objectives were to assess levels of genetic diversity amongst available accessions of diploid relatives closely related to tobacco (species of sections Sylvestres and Tomentosae), and to evaluate the efficacy of microsatellite markers for establishing species relationships in comparison with existing phylogenetic reconstructions. A subset of 46 primer combinations was therefore used to genotype 3 synthetic tobaccos and an expanded collection of 51 Nicotiana accessions representing 15 species. The average genetic similarity for 7 diverse accessions of tobacco was greater than the average similarity for N. otophora accessions, but lower than the average genetic similarities for N. sylvestris, N. tomentosa, N. kawakamii, and N. tomentosiformis accessions. A microsatellite-based phylogenetic tree was largely congruent with taxonomic representations based on morphological, cytological, and molecular observations. Results will be useful for selection of parents for creation of diploid mapping populations and for germplasm introgression activities.

Revolver is a new class of transposon-like gene composing the triticeae genome

Revolver discovered in the Triticeae plant is a novel class of transposon-like gene and a major component of the large cereal genome. An 89 bp segment of Revolver that is enriched in the genome of rye was isolated by deleting the DNA sequences common to rye and wheat. The entire structure of Revolver was determined by using rye genomic clones, which were screened by the 89 bp probe. Revolver consists of 2929-3041 bp with an inverted repeated sequence on each end and is dispersed through all seven chromosomes of the rye genome. Revolver is transcriptionally active, and the isolated full-length cDNA (726 bp) reveals that Revolver harbors a single gene consisting of three exons (342, 88, and 296 bp) and two introns (750 and 1237 bp), and encodes 139 amino acid residues of protein, which shows similarity to some transcriptional regulators. Revolver variants ranging from 2665 to 4269 bp, in which 5' regions were destructed, indicate structural diversities around the first exon. Revolver does not share identity with any known class I or class II autonomous transposable elements of any living species. DNA blot analysis of Triticeae plants shows that Revolver has existed since the diploid progenitor of wheat, and has been amplified or lost in several species during the evolution of the Triticeae.

The 172-kb genomic DNA region of the O. rufipogon yld1.1 locus: comparative sequence analysis with O. sativa ssp. japonica and O. sativa ssp. indica

Common wild rice (Oryza rufipogon) plays an important role by contributing to modern rice breeding. In this paper, we report the sequence and analysis of a 172-kb genomic DNA region of wild rice around the RM5 locus, which is associated with the yield QTL yld1.1. Comparative sequence analysis between orthologous RM5 regions from Oryza sativa ssp. japonica, O. sativa ssp. indica and O. rufipogon revealed a high level of conserved synteny in the content, homology, structure, orientation, and physical distance of all 14 predicted genes. Twelve of the putative genes were supported by matches to proteins with known function, whereas two were predicted by homology to rice and other plant expressed sequence tags or complementary DNAs. The remarkably high level of conservation found in coding, intronic and intergenic regions may indicate high evolutionary selection on the RM5 region. Although our analysis has not defined which gene(s) determine the yld1.1 phenotype, allelic variation and the insertion of transposable elements, among other nucleotide changes, represent potential variation responsible for the yield QTL. However, as suggested previously, two putative receptor-like protein kinase genes remain the key suspects for yld1.1.

Transposable elements in Coffea (Gentianales: Rubiacea) transcripts and their role in the origin of protein diversity in flowering plants

Transposable elements are major components of plant genomes and they influence their evolution, acting as recombination hot spots, acquiring specific cell functions or becoming part of protein-coding regions. The latter is the subject of the present analysis. This study is a report on the annotation of transposable elements (TEs) in expressed sequences of Coffea arabica, Coffea canephora and Coffea racemosa, showing the occurrence of 383 ESTs and 142 unigenes with TE fragments in these three Coffea species. Based on selected unigenes, it was possible to suggest 26 putative proteins with TE-cassette insertions, demonstrating a likely contribution to protein variability. The genes for two of those proteins, the fertility restorer (FR) and the pyrophosphate-dependent phosphofructokinase (PPi-PFKs) genes, were selected for evaluating the impact of TE-cassettes on host gene evolution of other plant genomes (Arabidopsis thaliana, Oryza sativa and Populus trichocarpa). This survey allowed identifying a FR gene in O. sativa harboring multiple insertions of LTR retrotransposons that originated new exons, which however does not necessarily mean a case of molecular domestication. A possible transduction event of a fragment of the PPi-PFK beta-subunit gene mediated by Helitron ATREPX1 in Arabidopsis thaliana was also highlighted.

TEnest: automated chronological annotation and visualization of nested plant transposable elements

Organisms with a high density of transposable elements (TEs) exhibit nesting, with subsequent repeats found inside previously inserted elements. Nesting splits the sequence structure of TEs and makes annotation of repetitive areas challenging. We present TEnest, a repeat identification and display tool made specifically for highly repetitive genomes. TEnest identifies repetitive sequences and reconstructs separated sections to provide full-length repeats and, for long-terminal repeat (LTR) retrotransposons, calculates age since insertion based on LTR divergence. TEnest provides a chronological insertion display to give an accurate visual representation of TE integration history showing timeline, location, and families of each TE identified, thus creating a framework from which evolutionary comparisons can be made among various regions of the genome. A database of repeats has been developed for maize (Zea mays), rice (Oryza sativa), wheat (Triticum aestivum), and barley (Hordeum vulgare) to illustrate the potential of TEnest software. All currently finished maize bacterial artificial chromosomes totaling 29.3 Mb were analyzed with TEnest to provide a characterization of the repeat insertions. Sixty-seven percent of the maize genome was found to be made up of TEs; of these, 95% are LTR retrotransposons. The rate of solo LTR formation is shown to be dissimilar across retrotransposon families. Phylogenetic analysis of TE families reveals specific events of extreme TE proliferation, which may explain the high quantities of certain TE families found throughout the maize genome. The TEnest software package is available for use on PlantGDB under the tools section (http://www.plantgdb.org/prj/TE_nest/TE_nest.html); the source code is available from (http://wiselab.org).

A high-throughput Agrobacterium-mediated transformation system for the grass model species Brachypodium distachyon L

In the ongoing process of developing Brachypodium distachyon as a model plant for temperate cereals and forage grasses, we have developed a high-throughput Agrobacterium-mediated transformation system for a diploid accession. Embryogenic callus, derived from immature embryos of the accession BDR018, were transformed with Agrobacterium tumefaciens strain AGL1 carrying two T-DNA plasmids, pDM805 and pWBV-Ds-Ubi-bar-Ds. Transient and stable transformation efficiencies were optimised by varying the pre-cultivation period, which had a strong effect on stable transformation efficiency. On average 55% of 17-day-old calli co-inoculated with Agrobacterium regenerated stable transgenic plants. Stable transformation frequencies of up to 80%, which to our knowledge is the highest transformation efficiency reported in graminaceous species, were observed. In a study of 177 transgenic lines transformed with pDM805, all of the regenerated transgenic lines were resistant to BASTA, while the gusA gene was expressed in 88% of the transgenic lines. Southern blot analysis revealed that 35% of the tested plants had a single T-DNA integration. Segregation analysis performed on progenies of ten selected T(0) plants indicated simple Mendelian inheritance of the two transgenes. Furthermore, the presence of two selection marker genes, bar and hpt, on the T-DNA of pWBV-Ds-Ubi-bar-Ds allowed us to characterize the developed transformation protocol with respect to full-length integration rate. Even when not selected for, full-length integration occurred in 97% of the transformants when using bialaphos as selection agent.

Evolutionary dynamics of an ancient retrotransposon family provides insights into evolution of genome size in the genus Oryza

Long terminal repeat (LTR) retrotransposons constitute a significant portion of most eukaryote genomes and can dramatically change genome size and organization. Although LTR retrotransposon content variation is well documented, the dynamics of genomic flux caused by their activity are poorly understood on an evolutionary time scale. This is primarily because of the lack of an experimental system composed of closely related species whose divergence times are within the limits of the ability to detect ancestrally related retrotransposons. The genus Oryza, with 24 species, ten genome types, different ploidy levels and over threefold genome size variation, constitutes an ideal experimental system to explore genus-level transposon dynamics. Here we present data on the discovery and characterization of an LTR retrotransposon family named RWG in the genus Oryza. Comparative analysis of transposon content (approximately 20 to 27,000 copies) and transpositional history of this family across the genus revealed a broad spectrum of independent and lineage-specific changes that have implications for the evolution of genome size and organization. In particular, we provide evidence that the basal GG genome of Oryza (O. granulata) has expanded by nearly 25% by a burst of the RWG lineage Gran3 subsequent to speciation. Finally we describe the recent evolutionary origin of Dasheng, a large retrotransposon derivative of the RWG family, specifically found in the A, B and C genome lineages of Oryza.

Differential impact of retrotransposon populations on the genome of allotetraploid tobacco (Nicotiana tabacum)

LTR-retrotransposons contribute substantially to the structural diversity of plant genomes. Recent models of genome evolution suggest that retrotransposon amplification is offset by removal of retrotransposon sequences, leading to a turnover of retrotransposon populations. While bursts of amplification have been documented, it is not known whether removal of retrotransposon sequences occurs continuously, or is triggered by specific stimuli over short evolutionary periods. In this work, we have characterized the evolutionary dynamics of four populations of copia-type retrotransposons in allotetraploid tobacco (Nicotiana tabacum) and its two diploid progenitors Nicotiana sylvestris and Nicotiana tomentosiformis. We have used SSAP (Sequence-Specific Amplification Polymorphism) to evaluate the contribution retrotransposons have made to the diversity of tobacco and its diploid progenitor species, to quantify the contribution each diploid progenitor has made to tobacco's retrotransposon populations, and to estimate losses or amplifications of retrotransposon sequences subsequent to tobacco's formation. Our results show that the tobacco genome derives from a turnover of retrotransposon sequences with removals concomitant with new insertions. We have detected unique behaviour specific to each retrotransposon population, with differences likely reflecting distinct evolutionary histories and activities of particular elements. Our results indicate that the retrotransposon content of a given plant species is strongly influenced by the host evolutionary history, with periods of rapid turnover of retrotransposon sequences stimulated by allopolyploidy.

New insights into Oryza genome evolution: high gene colinearity and differential retrotransposon amplification

An approximately 247-kb genomic region from FF genome of wild rice Oryza brachyantha, possessing the smallest Oryza genome, was compared to the orthologous approximately 450-kb region from AA genome, O. sativa L. ssp. japonica. 37 of 38 genes in the orthologous regions are shared between japonica and O. brachyantha. Analyses of nucleotide substitution in coding regions suggest the two genomes diverged approximately 10 million years ago. Comparisons of transposable elements (TEs) reveal that the density of DNA TEs in O. brachyantha is comparable to O. sativa; however, the density of RNA TEs is dramatically lower. The genomic fraction of RNA TEs in japonica is two times greater than in O. brachyantha. Differences, particularly in RNA TEs, in this region and in BAC end sequences from five wild and two cultivated Oryza species explain major genome size differences between sativa and brachyantha. Gene expression analyses of three ObDREB1 genes in the sequenced region indicate orthologous genes retain similar expression patterns following cold stress. Our results demonstrate that size and number of RNA TEs play a major role in genomic differentiation and evolution in Oryza. Additionally, distantly related O. brachyantha shares colinearity with O. sativa, offering opportunities to use comparative genomics to explore the genetic diversity of wild species to improve cultivated rice.

Analysis of EST sequences suggests recent origin of allotetraploid colonial and creeping bentgrasses

Advances in plant genomics have permitted the analysis of several members of the grass family, including the major domesticated species, and provided new insights into the evolution of the major crops on earth. Two members, colonial bentgrass (Agrostis capillaris L.) and creeping bentgrass (A. stolonifera L.) have only recently been domesticated and provide an interesting case of polyploidy and comparison to crops that have undergone human selection for thousands of years. As an initial step of characterizing these genomes, we have sampled roughly 10% of their gene content, thereby also serving as a starting point for the construction of their physical and genetic maps. Sampling mRNA from plants subjected to environmental stress showed a remarkable increase in transcription of transposable elements. Both colonial and creeping bentgrass are allotetraploids and are considered to have one genome in common, designated the A2 genome. Analysis of conserved genes present among the ESTs suggests the colonial and creeping bentgrass A2 genomes diverged from a common ancestor approximately 2.2 million years ago (MYA), thereby providing an enhanced evolutionary zoom in respect to the origin of maize, which formed 4.8 MYA, and tetraploid wheat, which formed only 0.5 MYA and is the progenitor of domesticated hexaploid wheat.

Orthologous comparison in a gene-rich region among grasses reveals stability in the sugarcane polyploid genome

Modern sugarcane (Saccharum spp.) is an important grass that contributes 60% of the raw sugar produced worldwide and has a high biofuel production potential. It was created about a century ago through hybridization of two highly polyploid species, namely S. officinarum and S. spontaneum. We investigated genome dynamics in this highly polyploid context by analyzing two homoeologous sequences (97 and 126 kb) in a region that has already been studied in several cereals. Our findings indicate that the two Saccharum species diverged 1.5-2 million years ago from one another and 8-9 million years ago from sorghum. The two sugarcane homoeologous haplotypes show perfect colinearity as well as high gene structure conservation. Apart from the insertion of a few retrotransposable elements, high homology was also observed for the non-transcribed regions. Relative to sorghum, the sugarcane sequences displayed colinearity, with the exception of two genes present only in sorghum, and striking homology in most non-coding parts of the genome. The gene distribution highlighted high synteny and colinearity with rice, and partial colinearity with each homoeologous maize region, which became perfect when the sequences were combined. The haplotypes observed in sugarcane may thus closely represent the ancestral Andropogoneae haplotype. This analysis of sugarcane haplotype organization at the sequence level suggests that the high ploidy in sugarcane did not induce generalized reshaping of its genome, thus challenging the idea that polyploidy quickly induces generalized rearrangement of genomes. These results also confirm the view that sorghum is the model of choice for sugarcane.