Uneven chromosome contraction and expansion in the maize genome

Production and processing of siRNA precursor transcripts from the highly repetitive maize genome

Mutations affecting the maintenance of heritable epigenetic states in maize identify multiple RNA–directed DNA methylation (RdDM) factors including RMR1, a novel member of a plant-specific clade of Snf2-related proteins. Here we show that RMR1 is necessary for the accumulation of a majority of 24 nt small RNAs, including those derived from Long-Terminal Repeat (LTR) retrotransposons, the most common repetitive feature in the maize genome. A genetic analysis of DNA transposon repression indicates that RMR1 acts upstream of the RNA–dependent RNA polymerase, RDR2 (MOP1). Surprisingly, we show that non-polyadenylated transcripts from a sampling of LTR retrotransposons are lost in both rmr1 and rdr2 mutants. In contrast, plants deficient for RNA Polymerase IV (Pol IV) function show an increase in polyadenylated LTR RNA transcripts. These findings support a model in which Pol IV functions independently of the small RNA accumulation facilitated by RMR1 and RDR2 and support that a loss of Pol IV leads to RNA Polymerase II–based transcription. Additionally, the lack of changes in general genome homeostasis in rmr1 mutants, despite the global loss of 24 nt small RNAs, challenges the perceived roles of siRNAs in maintaining functional heterochromatin in the genomes of outcrossing grass species.

Most eukaryotic genomes are divided into two functional classes of regulation: the euchromatic and the heterochromatic. Heterochromatic regions, often composed of potentially deleterious transposons and retrotransposons, are typically viewed as “silent” or not transcribed. Paradoxically, evidence from multiple organisms indicates that heterochromatic regions must be transcribed to maintain a heterochromatic character. In plants, specialized RNA polymerase complexes are thought to specifically process repetitive regions of the genome into small RNA molecules that facilitate maintenance of a heterochromatic environment. We investigated the role of this specialized polymerase pathway in maintaining maize genome homeostasis with particular focus on RMR1, a novel protein related to a family of DNA repair proteins, whose function in modifying repetitive regions of the genome is unknown. We find most small RNA generation is dependent on RMR1, which appears to function downstream of the specialized polymerase, RNA polymerase IV. However, we provide evidence that the function of RNA polymerase IV is not disrupted by the absence of small RNA generation. Our results suggest the division of the plant genome into euchromatin and heterochromatin is maintained by template competition between the specialized plant polymerases and canonical RNA polymerase II, and not by the subsequent generation of small RNA molecules.

Advances in maize genomics and their value for enhancing genetic gains from breeding

Maize is an important crop for food, feed, forage, and fuel across tropical and temperate areas of the world. Diversity studies at genetic, molecular, and functional levels have revealed that, tropical maize germplasm, landraces, and wild relatives harbor a significantly wider range of genetic variation. Among all types of markers, SNP markers are increasingly the marker-of-choice for all genomics applications in maize breeding. Genetic mapping has been developed through conventional linkage mapping and more recently through linkage disequilibrium-based association analyses. Maize genome sequencing, initially focused on gene-rich regions, now aims for the availability of complete genome sequence. Conventional insertion mutation-based cloning has been complemented recently by EST- and map-based cloning. Transgenics and nutritional genomics are rapidly advancing fields targeting important agronomic traits including pest resistance and grain quality. Substantial advances have been made in methodologies for genomics-assisted breeding, enhancing progress in yield as well as abiotic and biotic stress resistances. Various genomic databases and informatics tools have been developed, among which MaizeGDB is the most developed and widely used by the maize research community. In the future, more emphasis should be given to the development of tools and strategic germplasm resources for more effective molecular breeding of tropical maize products.

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

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

The 'inner circle' of the cereal genomes

Early marker-based macrocolinearity studies between the grass genomes led to arranging their chromosomes into concentric 'crop circles' of synteny blocks that initially consisted of 30 rice-independent linkage groups representing the ancestral cereal genome structure. Recently, increased marker density and genome sequencing of several cereal genomes allowed the characterization of intragenomic duplications and their integration with intergenomic colinearity data to identify paleo-duplications and propose a model for the evolution of the grass genomes from a common ancestor. On the basis of these data an 'inner circle' comprising five ancestral chromosomes was defined providing a new reference for the grass chromosomes and new insights into their ancestral relationships and origin, as well as an efficient tool to design cross-genome markers for genetic studies.

Methylation-sensitive linking libraries enhance gene-enriched sequencing of complex genomes and map DNA methylation domains

Background

Many plant genomes are resistant to whole-genome assembly due to an abundance of repetitive sequence, leading to the development of gene-rich sequencing techniques. Two such techniques are hypomethylated partial restriction (HMPR) and methylation spanning linker libraries (MSLL). These libraries differ from other gene-rich datasets in having larger insert sizes, and the MSLL clones are designed to provide reads localized to "epigenetic boundaries" where methylation begins or ends.

Results

A large-scale study in maize generated 40,299 HMPR sequences and 80,723 MSLL sequences, including MSLL clones exceeding 100 kb. The paired end reads of MSLL and HMPR clones were shown to be effective in linking existing gene-rich sequences into scaffolds. In addition, it was shown that the MSLL clones can be used for anchoring these scaffolds to a BAC-based physical map. The MSLL end reads effectively identified epigenetic boundaries, as indicated by their preferential alignment to regions upstream and downstream from annotated genes. The ability to precisely map long stretches of fully methylated DNA sequence is a unique outcome of MSLL analysis, and was also shown to provide evidence for errors in gene identification. MSLL clones were observed to be significantly more repeat-rich in their interiors than in their end reads, confirming the correlation between methylation and retroelement content. Both MSLL and HMPR reads were found to be substantially gene-enriched, with the SalI MSLL libraries being the most highly enriched (31% align to an EST contig), while the HMPR clones exhibited exceptional depletion of repetitive DNA (to ~11%). These two techniques were compared with other gene-enrichment methods, and shown to be complementary.

Conclusion

MSLL technology provides an unparalleled approach for mapping the epigenetic status of repetitive blocks and for identifying sequences mis-identified as genes. Although the types and natures of epigenetic boundaries are barely understood at this time, MSLL technology flags both approximate boundaries and methylated genes that deserve additional investigation. MSLL and HMPR sequences provide a valuable resource for maize genome annotation, and are a uniquely valuable complement to any plant genome sequencing project. In order to make these results fully accessible to the community, a web display was developed that shows the alignment of MSLL, HMPR, and other gene-rich sequences to the BACs; this display is continually updated with the latest ESTs and BAC sequences.

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.

Distinct size distribution of endogeneous siRNAs in maize: Evidence from deep sequencing in the mop1-1 mutant

Small RNAs from plants are known to be highly complex and abundant, with this complexity proportional to genome size. Most endogenous siRNAs in Arabidopsis are dependent on RNA-DEPENDENT RNA POLYMERASE 2 (RDR2) for their biogenesis. Recent work has demonstrated that the maize MEDIATOR OF PARAMUTATION1 (mop1) gene is a predicted ortholog of RDR2. The mop1 gene is required for establishment of paramutation and maintenance of transcriptional silencing of transposons and transgenes, suggesting the potential involvement of small RNAs. We analyzed small RNAs in wild-type maize and in the isogenic mop1-1 loss-of-function mutant by using Illumina's sequencing-by-synthesis (SBS) technology, which allowed us to characterize the complement of maize small RNAs to considerable depth. Similar to rdr2 in Arabidopsis, in mop1-1, the 24-nucleotide (nt) endogenous heterochromatic short-interfering siRNAs were dramatically reduced, resulting in an enrichment of miRNAs and transacting siRNAs. In contrast to the Arabidopsis rdr2 mutant, the mop1-1 plants retained a highly abundant heterochromatic approximately 22-nt class of small RNAs, suggesting a second mechanism for heterochromatic siRNA production. The enrichment of miRNAs and loss of 24-nt heterochromatic siRNAs in mop1-1 should be advantageous for miRNA discovery as the maize genome becomes more fully sequenced.

Organization of the prolamin gene family provides insight into the evolution of the maize genome and gene duplications in grass species

Zea mays, commonly known as corn, is perhaps the most greatly produced crop in terms of tonnage and a major food, feed, and biofuel resource. Here we analyzed its prolamin gene family, encoding the major seed storage proteins, as a model for gene evolution by syntenic alignments with sorghum and rice, two genomes that have been sequenced recently. Because a high-density gene map has been constructed for maize inbred B73, all prolamin gene copies can be identified in their chromosomal context. Alignment of respective chromosomal regions of these species via conserved genes allow us to identify the pedigree of prolamin gene copies in space and time. Its youngest and largest gene family, the alpha prolamins, arose about 22-26 million years ago (Mya) after the split of the Panicoideae (including maize, sorghum, and millet) from the Pooideae (including wheat, barley, and oats) and Oryzoideae (rice). The first dispersal of alpha prolamin gene copies occurred before the split of the progenitors of maize and sorghum about 11.9 Mya. One of the two progenitors of maize gained a new alpha zein locus, absent in the other lineage, to form a nonduplicated locus in maize after allotetraplodization about 4.8 Mya. But dispersed copies gave rise to tandem duplications through uneven expansion and gene silencing of this gene family in maize and sorghum, possibly because of maize's greater recombination and mutation rates resulting from its diploidization process. Interestingly, new gene loci in maize represent junctions of ancestral chromosome fragments and sites of new centromeres in sorghum and rice.

Diverged copies of the seed regulatory Opaque-2 gene by a segmental duplication in the progenitor genome of rice, sorghum, and maize

Comparative analyses of the sequence of entire genomes have shown that gene duplications, chromosomal segmental duplications, or even whole genome duplications (WGD) have played prominent roles in the evolution of many eukaryotic species. Here, we used the ancient duplication of a well known transcription factor in maize, encoded by the Opaque-2 (O2) locus, to examine the general features of divergences of chromosomal segmental duplications in a lineage-specific manner. We took advantage of contiguous chromosomal sequence information in rice (Oryza sativa, Nipponbare), sorghum (Sorghum bicolor, Btx623), and maize (Zea mays, B73) that were aligned by conserved gene order (synteny). This analysis showed that the maize O2 locus is contained within a 1.25 million base-pair (Mb) segment on chromosome 7, which was duplicated approximately 56 million years ago (mya) before the split of rice and maize 50 mya. The duplicated region on chromosome 1 is only half the size and contains the maize OHP gene, which does not restore the o2 mutation although it encodes a protein with the same DNA and protein binding properties in endosperm. The segmental duplication is not only found in rice, but also in sorghum, which split from maize 11.9 mya. A detailed analysis of the duplicated regions provided examples for complex rearrangements including deletions, duplications, conversions, inversions, and translocations. Furthermore, the rice and sorghum genomes appeared to be more stable than the maize genome, probably because maize underwent allotetraploidization and then diploidization.

Involvement of the MADS-box gene ZMM4 in floral induction and inflorescence development in maize

The switch from vegetative to reproductive growth is marked by the termination of vegetative development and the adoption of floral identity by the shoot apical meristem (SAM). This process is called the floral transition. To elucidate the molecular determinants involved in this process, we performed genome-wide RNA expression profiling on maize (Zea mays) shoot apices at vegetative and early reproductive stages using massively parallel signature sequencing technology. Profiling revealed significant up-regulation of two maize MADS-box (ZMM) genes, ZMM4 and ZMM15, after the floral transition. ZMM4 and ZMM15 map to duplicated regions on chromosomes 1 and 5 and are linked to neighboring MADS-box genes ZMM24 and ZMM31, respectively. This gene order is syntenic with the vernalization1 locus responsible for floral induction in winter wheat (Triticum monococcum) and similar loci in other cereals. Analyses of temporal and spatial expression patterns indicated that the duplicated pairs ZMM4-ZMM24 and ZMM15-ZMM31 are coordinately activated after the floral transition in early developing inflorescences. More detailed analyses revealed ZMM4 expression initiates in leaf primordia of vegetative shoot apices and later increases within elongating meristems acquiring inflorescence identity. Expression analysis in late flowering mutants positioned all four genes downstream of the floral activators indeterminate1 (id1) and delayed flowering1 (dlf1). Overexpression of ZMM4 leads to early flowering in transgenic maize and suppresses the late flowering phenotype of both the id1 and dlf1 mutations. Our results suggest ZMM4 may play roles in both floral induction and inflorescence development.

Physical and genetic structure of the maize genome reflects its complex evolutionary history

Maize (Zea mays L.) is one of the most important cereal crops and a model for the study of genetics, evolution, and domestication. To better understand maize genome organization and to build a framework for genome sequencing, we constructed a sequence-ready fingerprinted contig-based physical map that covers 93.5% of the genome, of which 86.1% is aligned to the genetic map. The fingerprinted contig map contains 25,908 genic markers that enabled us to align nearly 73% of the anchored maize genome to the rice genome. The distribution pattern of expressed sequence tags correlates to that of recombination. In collinear regions, 1 kb in rice corresponds to an average of 3.2 kb in maize, yet maize has a 6-fold genome size expansion. This can be explained by the fact that most rice regions correspond to two regions in maize as a result of its recent polyploid origin. Inversions account for the majority of chromosome structural variations during subsequent maize diploidization. We also find clear evidence of ancient genome duplication predating the divergence of the progenitors of maize and rice. Reconstructing the paleoethnobotany of the maize genome indicates that the progenitors of modern maize contained ten chromosomes.

As a cash crop and a model biological system, maize is of great public interest. To facilitate maize molecular breeding and its basic biology research, we built a high-resolution physical map with two different fingerprinting methods on the same set of bacterial artificial chromosome clones. The physical map was integrated to a high-density genetic map and further serves as a framework for the maize genome-sequencing project. Comparative genomics showed that the euchromatic regions between rice and maize are very conserved. Physically we delimited these conserved regions and thus detected many genome rearrangements. We defined extensively the duplication blocks within the maize genome. These blocks allowed us to reconstruct the chromosomes of the maize progenitor. We detected that maize genome has experienced two rounds of genome duplications, an ancient one before maize–rice divergence and a recent one after tetraploidization.

Development and evaluation of a high-throughput, low-cost genotyping platform based on oligonucleotide microarrays in rice

BACKGROUND

We report the development of a microarray platform for rapid and cost-effective genetic mapping, and its evaluation using rice as a model. In contrast to methods employing whole-genome tiling microarrays for genotyping, our method is based on low-cost spotted microarray production, focusing only on known polymorphic features.

RESULTS

We have produced a genotyping microarray for rice, comprising 880 single feature polymorphism (SFP) elements derived from insertions/deletions identified by aligning genomic sequences of the japonica cultivar Nipponbare and the indica cultivar 93-11. The SFPs were experimentally verified by hybridization with labeled genomic DNA prepared from the two cultivars. Using the genotyping microarrays, we found high levels of polymorphism across diverse rice accessions, and were able to classify all five subpopulations of rice with high bootstrap support. The microarrays were used for mapping of a gene conferring resistance to Magnaporthe grisea, the causative organism of rice blast disease, by quantitative genotyping of samples from a recombinant inbred line population pooled by phenotype.

CONCLUSION

We anticipate this microarray-based genotyping platform, based on its low cost-per-sample, to be particularly useful in applications requiring whole-genome molecular marker coverage across large numbers of individuals.

Genetics of intrinsic postzygotic isolation in a circumpolar plant species, Draba nivalis (Brassicaceae)

Sterility barriers, ranging from incomplete to fully developed, were recently demonstrated within taxonomic species of the genus Draba, suggesting the existence of numerous, cryptic biological species. Because these taxa are predominately selfers and of Pleistocene origin, it was concluded that hybrid sterility evolved quickly and possibly by genetic drift. Here we used genetic mapping and QTL analyses to determine the genetic basis of hybrid sterility between geographically distant populations of one of these taxonomic species, Draba nivalis. Fifty microsatellite loci were mapped, and QTL analyses identified five loci underlying seed fertility and two underlying pollen fertility. Four of five seed fertility QTLs reduced fertility in heterozygotes, an observation most consistent with drift-based fixation of underdominant sterility loci. However, several nuclear-nuclear interactions were also found, including two that acted like reciprocal translocations with lowest fitness in double heterozygotes, and two that had a pattern of fitness consistent with Bateson-Dobzhansky-Muller incompatibilities. In contrast, pollen fertility QTLs exhibited additive inheritance, with lowest fertility associated with the paternal allele, a pattern of inheritance suggestive of cytonuclear incompatibilities. The results imply that multiple genetic mechanisms underlie the rapid evolution of reproductive barriers in Draba.

Maize genome structure variation: interplay between retrotransposon polymorphisms and genic recombination

Although maize (Zea mays) retrotransposons are recombinationally inert, the highly polymorphic structure of maize haplotypes raises questions regarding the local effect of intergenic retrotransposons on recombination. To examine this effect, we compared recombination in the same genetic interval with and without a large retrotransposon cluster. We used three different bz1 locus haplotypes, McC, B73, and W22, in the same genetic background. We analyzed recombination between the bz1 and stc1 markers in heterozygotes that differ by the presence and absence of a 26-kb intergenic retrotransposon cluster. To facilitate the genetic screen, we used Ds and Ac markers that allowed us to identify recombinants by their seed pigmentation. We sequenced 239 recombination junctions and assigned them to a single nucleotide polymorphism-delimited interval in the region. The genetic distance between the markers was twofold smaller in the presence of the retrotransposon cluster. The reduction was seen in bz1 and stc1, but no recombination occurred in the highly polymorphic intergenic region of either heterozygote. Recombination within genes shuffled flanking retrotransposon clusters, creating new chimeric haplotypes and either contracting or expanding the physical distance between markers. Our findings imply that haplotype structure will profoundly affect the correlation between genetic and physical distance for the same interval in maize.

A genomic and expression compendium of the expanded PEBP gene family from maize

The phosphatidylethanolamine-binding proteins (PEBPs) represent an ancient protein family found across the biosphere. In animals they are known to act as kinase and serine protease inhibitors controlling cell growth and differentiation. In plants the most extensively studied PEBP genes, the Arabidopsis (Arabidopsis thaliana) FLOWERING LOCUS T (FT) and TERMINAL FLOWER1 (TFL1) genes, function, respectively, as a promoter and a repressor of the floral transition. Twenty-five maize (Zea mays) genes that encode PEBP-like proteins, likely the entire gene family, were identified and named Zea mays CENTRORADIALIS (ZCN), after the first described plant PEBP gene from Antirrhinum. The maize family is expanded relative to eudicots (typically six to eight genes) and rice (Oryza sativa; 19 genes). Genomic structures, map locations, and syntenous relationships with rice were determined for 24 of the maize ZCN genes. Phylogenetic analysis assigned the maize ZCN proteins to three major subfamilies: TFL1-like (six members), MOTHER OF FT AND TFL1-like (three), and FT-like (15). Expression analysis demonstrated transcription for at least 21 ZCN genes, many with developmentally specific patterns and some having alternatively spliced transcripts. Expression patterns and protein structural analysis identified maize candidates likely having conserved gene function of TFL1. Expression patterns and interaction of the ZCN8 protein with the floral activator DLF1 in the yeast (Saccharomyces cerevisiae) two-hybrid assay strongly supports that ZCN8 plays an orthologous FT function in maize. The expression of other ZCN genes in roots, kernels, and flowers implies their involvement in diverse developmental processes.

A transgenomic cytogenetic sorghum (Sorghum propinquum) bacterial artificial chromosome fluorescence in situ hybridization map of maize (Zea mays L.) pachytene chromosome 9, evidence for regions of genome hyperexpansion

A cytogenetic FISH map of maize pachytene-stage chromosome 9 was produced with 32 maize marker-selected sorghum BACs as probes. The genetically mapped markers used are distributed along the linkage maps at an average spacing of 5 cM. Each locus was mapped by means of multicolor direct FISH with a fluorescently labeled probe mix containing a whole-chromosome paint, a single sorghum BAC clone, and the centromeric sequence, CentC. A maize-chromosome-addition line of oat was used for bright unambiguous identification of the maize 9 fiber within pachytene chromosome spreads. The locations of the sorghum BAC-FISH signals were determined, and each new cytogenetic locus was assigned a centiMcClintock position on the short (9S) or long (9L) arm. Nearly all of the markers appeared in the same order on linkage and cytogenetic maps but at different relative positions on the two. The CentC FISH signal was localized between cdo17 (at 9L.03) and tda66 (at 9S.03). Several regions of genome hyperexpansion on maize chromosome 9 were found by comparative analysis of relative marker spacing in maize and sorghum. This transgenomic cytogenetic FISH map creates anchors between various maps of maize and sorghum and creates additional tools and information for understanding the structure and evolution of the maize genome.

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.

Rapid evolution and complex structural organization in genomic regions harboring multiple prolamin genes in the polyploid wheat genome

Genes encoding wheat prolamins belong to complicated multi-gene families in the wheat genome. To understand the structural complexity of storage protein loci, we sequenced and analyzed orthologous regions containing both gliadin and LMW-glutenin genes from the A and B genomes of a tetraploid wheat species, Triticum turgidum ssp. durum. Despite their physical proximity to one another, the gliadin genes and LMW-glutenin genes are organized quite differently. The gliadin genes are found to be more clustered than the LMW-glutenin genes which are separated from each other by much larger distances. The separation of the LMW-glutenin genes is the result of both the insertion of large blocks of repetitive DNA owing to the rapid amplification of retrotransposons and the presence of genetic loci interspersed between them. Sequence comparisons of the orthologous regions reveal that gene movement could be one of the major factors contributing to the violation of microcolinearity between the homoeologous A and B genomes in wheat. The rapid sequence rearrangements and differential insertion of repetitive DNA has caused the gene islands to be not conserved in compared regions. In addition, we demonstrated that the i-type LMW-glutenin originated from a deletion of 33-bps in the 5' coding region of the m-type gene. Our results show that multiple rounds of segmental duplication of prolamin genes have driven the amplification of the omega-gliadin genes in the region; such segmental duplication could greatly increase the repetitive DNA content in the genome depending on the amount of repetitive DNA present in the original duplicate region.

Evolutionary history and positional shift of a rice centromere

Rice centromere 8 was previously proposed to be an "immature" centromere that recently arose from a genic region. Our comparative genomics analysis indicates that Cen8 was formed at its current location at least 7-9 million years ago and was physically shifted by a more recent inversion of a segment spanning centromeric and pericentromeric regions.

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.

Sequence-indexed mutations in maize using the UniformMu transposon-tagging population

BACKGROUND

Gene knockouts are a critical resource for functional genomics. In Arabidopsis, comprehensive knockout collections were generated by amplifying and sequencing genomic DNA flanking insertion mutants. These Flanking Sequence Tags (FSTs) map each mutant to a specific locus within the genome. In maize, FSTs have been generated using DNA transposons. Transposable elements can generate unstable insertions that are difficult to analyze for simple knockout phenotypes. Transposons can also generate somatic insertions that fail to segregate in subsequent generations.

RESULTS

Transposon insertion sites from 106 UniformMu FSTs were tested for inheritance by locus-specific PCR. We confirmed 89% of the FSTs to be germinal transposon insertions. We found no evidence for somatic insertions within the 11% of insertion sites that were not confirmed. Instead, this subset of insertion sites had errors in locus-specific primer design due to incomplete or low-quality genomic sequences. The locus-specific PCR assays identified a knockout of a 6-phosphogluconate dehydrogenase gene that co-segregates with a seed mutant phenotype. The mutant phenotype linked to this knockout generates novel hypotheses about the role for the plastid-localized oxidative pentose phosphate pathway during grain-fill.

CONCLUSION

We show that FSTs from the UniformMu population identify stable, germinal insertion sites in maize. Moreover, we show that these sequence-indexed mutations can be readily used for reverse genetic analysis. We conclude from these data that the current collection of 1,882 non-redundant insertion sites from UniformMu provide a genome-wide resource for reverse genetics.

The maize (Zea mays L.) RTCS gene encodes a LOB domain protein that is a key regulator of embryonic seminal and post-embryonic shoot-borne root initiation

Maize has a complex root system composed of different root types formed during different stages of development. The rtcs (rootless concerning crown and seminal roots) mutant is impaired in the initiation of the embryonic seminal roots and the post-embryonic shoot-borne root system. The primary root of the mutant shows a reduced gravitropic response, while its elongation, lateral root density and reaction to exogenously applied auxin is not affected. We report here the map-based cloning of the RTCS gene which encodes a 25.5 kDa LOB domain protein located on chromosome 1S. The RTCS gene has been duplicated during evolution. The RTCS-LIKE (RTCL) gene displays 72% sequence identity on the protein level. Both genes are preferentially expressed in roots. Expression of RTCS in coleoptilar nodes is confined to emerging shoot-borne root primordia. Sequence analyses of the RTCS and RTCL upstream genomic regions identified auxin response elements. Reverse transcriptase-PCR revealed that both genes are auxin induced. Microsynteny analyses between maize and rice genomes revealed co-linearity of 14 genes in the RTCS region. We conclude from our data that RTCS and RTCL are auxin-responsive genes involved in the early events that lead to the initiation and maintenance of seminal and shoot-borne root primordia formation.

MIPSPlantsDB--plant database resource for integrative and comparative plant genome research

Genome-oriented plant research delivers rapidly increasing amount of plant genome data. Comprehensive and structured information resources are required to structure and communicate genome and associated analytical data for model organisms as well as for crops. The increase in available plant genomic data enables powerful comparative analysis and integrative approaches. PlantsDB aims to provide data and information resources for individual plant species and in addition to build a platform for integrative and comparative plant genome research. PlantsDB is constituted from genome databases for Arabidopsis, Medicago, Lotus, rice, maize and tomato. Complementary data resources for cis elements, repetive elements and extensive cross-species comparisons are implemented. The PlantsDB portal can be reached at .

Recent proliferation and translocation of pollen group 1 allergen genes in the maize genome

The dominant allergenic components of grass pollen are known by immunologists as group 1 allergens. These constitute a set of closely related proteins from the beta-expansin family and have been shown to have cell wall-loosening activity. Group 1 allergens may facilitate the penetration of pollen tubes through the grass stigma and style. In maize (Zea mays), group 1 allergens are divided into two classes, A and B. We have identified 15 genes encoding group 1 allergens in maize, 11 genes in class A and four genes in class B, as well as seven pseudogenes. The genes in class A can be divided by sequence relatedness into two complexes, whereas the genes in class B constitute a single complex. Most of the genes identified are represented in pollen-specific expressed sequence tag libraries and are under purifying selection, despite the presence of multiple copies that are nearly identical. Group 1 allergen genes are clustered in at least six different genomic locations. The single class B location and one of the class A locations show synteny with the rice (Oryza sativa) regions where orthologous genes are found. Both classes are expressed at high levels in mature pollen but at low levels in immature flowers. The set of genes encoding maize group 1 allergens is more complex than originally anticipated. If this situation is common in grasses, it may account for the large number of protein variants, or group 1 isoallergens, identified previously in turf grass pollen by immunologists.

Comparison of orthologous loci from small grass genomes Brachypodium and rice: implications for wheat genomics and grass genome annotation

Brachypodium sylvaticum and Brachypodium distachyon were recently proposed as new model plants because of their small genomes and their phylogenetic position between rice and Triticeae crops. We sequenced a 371-kb region in B. sylvaticum, the largest genomic sequence available so far from this species, providing quantitative data on gene conservation, collinearity and phylogeny. We compared it with orthologous regions from rice and wheat. Brachypodium and wheat show perfect macro-collinearity of genetic markers, whereas rice contains an approximately 220-kb inversion. Rice contains almost twice as many genes as Brachypodium in the region studied, whereas wheat has about 40% more. Through comparative annotation, we identified alternative transcripts and improved the annotation for several rice genes, indicating that approximately 15% of rice genes might require re-annotation. Surprisingly, our data suggest that 10-15% of functional sequences in small grass genomes may not encode any proteins. From available genomic and expressed sequence tag sequences, we estimated Brachypodium to have diverged from wheat about 35-40 Mya, significantly more recently than the divergence of rice and wheat. However, our data also indicate that orthologous regions from Brachypodium and wheat differ considerably in gene content, thus the Brachypodium genome sequence probably cannot replace genomic studies in the large Triticeae genomes.

Maize haplotype with a helitron-amplified cytidine deaminase gene copy

Background

Genetic maps are based on recombination of orthologous gene sequences between different strains of the same species. Therefore, it was unexpected to find extensive non-collinearity of genes between different inbred strains of maize. Interestingly, disruption of gene collinearity can be caused among others by a rolling circle-type copy and paste mechanism facilitated by Helitrons. However, understanding the role of this type of gene amplification has been hampered by the lack of finding intact gene sequences within Helitrons.

Results

By aligning two haplotypes of the z1C1 locus of maize we found a Helitron that contains two genes, one encoding a putative cytidine deaminase and one a hypothetical protein with part of a 40S ribosomal protein. The cytidine deaminase gene, called ZmCDA3, has been copied from the ZmCDA1 gene on maize chromosome 7 about 4.5 million years ago (mya) after maize was formed by whole-genome duplication from two progenitors. Inbred lines contain gene copies of both progenitors, the ZmCDA1 and ZmCDA2 genes. Both genes diverged when the progenitors of maize split and are derived from the same progenitor as the rice OsCDA1 gene. The ZmCDA1 and ZmCDA2 genes are both transcribed in leaf and seed tissue, but transcripts of the paralogous ZmCDA3 gene have not been found yet. Based on their protein structure the maize CDA genes encode a nucleoside deaminase that is found in bacterial systems and is distinct from the mammalian RNA and/or DNA modifying enzymes.

Conclusion

The conservation of a paralogous gene sequence encoding a cytidine deaminase gene over 4.5 million years suggests that Helitrons could add functional gene sequences to new chromosomal positions and thereby create new haplotypes. However, the function of such paralogous gene copies cannot be essential because they are not present in all maize strains. However, it is interesting to note that maize hybrids can outperform their inbred parents. Therefore, certain haplotypes may function only in combination with other haplotypes or under specialized environmental conditions.