Mosaic organization of orthologous sequences in grass genomes

Major impacts of widespread structural variation on sorghum

Genetic diversity is critical to crop breeding and improvement, and dissection of the genomic variation underlying agronomic traits can both assist breeding and give insight into basic biological mechanisms. Although recent genome analyses in plants reveal many structural variants (SVs), most current studies of crop genetic variation are dominated by single-nucleotide polymorphisms (SNPs). The extent of the impact of SVs on global trait variation, as well as their utility in genome-wide selection, is not yet understood. In this study, we built an SV data set based on whole-genome resequencing of diverse sorghum lines (n = 363), validated the correlation of photoperiod sensitivity and variety type, and identified SV hotspots underlying the divergent evolution of cellulosic and sweet sorghum. In addition, we showed the complementary contribution of SVs for heritability of traits related to sorghum adaptation. Importantly, inclusion of SV polymorphisms in association studies revealed genotype-phenotype associations not observed with SNPs alone. Three-way genome-wide association studies (GWAS) based on whole-genome SNP, SV, and integrated SNP + SV data sets showed substantial associations between SVs and sorghum traits. The addition of SVs to GWAS substantially increased heritability estimates for some traits, indicating their important contribution to functional allelic variation at the genome level. Our discovery of the widespread impacts of SVs on heritable gene expression variation could render a plausible mechanism for their disproportionate impact on phenotypic variation. This study expands our knowledge of SVs and emphasizes the extensive impacts of SVs on sorghum.

The regulation of zein biosynthesis in maize endosperm

We review the current knowledge regarding the regulation of zein storage proteins biosynthesis and protein body formation, which are crucial processes for the successful accumulation of nutrients in maize kernels. Storage proteins in the seeds of crops in the grass family (Poaceae) are a major source of dietary protein for humans. In maize (Zea mays), proteins are the second largest nutrient component in the kernels, accounting for ~ 10% of the kernel weight. Over half of the storage proteins in maize kernels are zeins, which lack two essential amino acids, lysine and tryptophan. This deficiency limits the use of maize proteins in the food and feed industries. Zeins are encoded by a large super-gene family. During endosperm development, zeins accumulate in protein bodies, which are derived from the rough endoplasmic reticulum. In recent years, our knowledge of the pathways of zein biosynthesis and their deposition within the endosperm has been greatly expanded. In this review, we summarize the current understanding of zeins, including the genes encoding these proteins, their expression patterns and transcriptional regulation, the process of protein body formation, and other biological processes affecting zein accumulation.

Structure, Organization, and Expression of the Alpha Prolamin Multigenic Family Bring New Insights into the Evolutionary Relationships among Grasses

Prolamins are the major seed storage proteins of grasses. In maize and related species, prolamins are classified into α-, β-, γ-, and δ-subclasses by their solubility properties. α-prolamins are encoded by multigene families and have a secondary structure that consists of tandem α-helix repeats. Maize has two α-prolamin subclasses, namely the 19 and 22 kDa subclasses that contain nine and 10 α-helix repeats, respectively. Here, we present an evolutionary study based on the structure, organization, and expression of α-prolamins in maize, sugarcane, sorghum, and coix. True 22 kDa subclasses containing 10 repeats are conserved in all four species, but true 19 kDa subclasses containing nine repeats are found only in maize and sugarcane. We discovered a 19 kDa-like α-coixin that, as in sorghum, is encoded by few genes. These data suggest that a 19 kDa progenitor present in the ancestor common to maize, coix, sorghum, and sugarcane was preserved at low copy number in coix and sorghum, while amplified into multigene family architecture in maize and sugarcane. The expression profiling of α-prolamins, verified by two-dimensional gels, showed highly conserved multispot composition for the 19 kDa α-prolamins in maize and sugarcane. Coix and sorghum did not present true 19 kDa α-prolamin spots. Our data show remarkable similarity between maize and sugarcane 19 kDa α-prolamins regarding both gene structure and expression. Since the multigene architecture of 19 kDa α-canein appeared after sugarcane diverged from sorghum, our data suggest that maize and sugarcane might have acquired the multigene family encoding these storage proteins from a common ancestor.

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.

QTL for fibre-related traits in grain × sweet sorghum as a tool for the enhancement of sorghum as a biomass crop

Compared to maize and temperate grasses, sorghum has received less attention in terms of improving cell wall components. The objectives of this study were to identify quantitative trait loci (QTL) with main effects, epistatic and pleiotropic effects along with QTL × environment (QE) interactions controlling fibre-related traits in sorghum. Neutral detergent fibre (NDF), acid detergent fibre (ADF), acid detergent lignin (ADL), cellulose, hemicellulose, fresh leaf mass, stripped stalk mass, dry stalk mass, fresh biomass and dry biomass were analysed from a population of 188 grain × sweet sorghum recombinant inbred lines. A genetic map consisting of 157 DNA markers was constructed, and QTL were detected using composite interval mapping (CIM). CIM detected more than 5 additive QTL per trait explaining 7.1-24.7% of the phenotypic variation. Abundant co-localization of these QTL was observed across all chromosomes, and the highest cluster was identified on chromosome 6. Searching for candidate genes using the confidence interval of our QTL clusters reveals that these clusters might comprise a set of genes that are tightly linked. Some QTL showed multiple effects; however, the allele for each trait was favouring the parent with the increasing effect. QE interactions were observed for QTL showing multiple effects. Additive × additive interaction was observed for 7 out of 10 traits, indicating the importance of epistatic analysis. However, the phenotypic variation explained by digenic interactions was lower compared to the individual QTL. Our results indicate that various genetic components contribute to fibre-related traits and should be considered during the enhancement of sorghum for lignocellulosic biomass.

Recent insertion of a 52-kb mitochondrial DNA segment in the wheat lineage

The assembly of a 1.3-Mb size region of the wheat genome has provided the opportunity to study a recent nuclear mitochondrial DNA insertion (NUMT). In the present study, we have studied two bacterial artificial chromosomes (BACs) and characterized a 52-kb NUMT segment from the tetraploid and hexaploid wheat BAC libraries. The conserved orthologous NUMT regions from tetraploid and hexaploid wheat Langdon and Chinese Spring shared identical gene haplotypes even though mutations (insertions, deletions, and substitutions) had occurred. The 52-kb NUMT was present in hexaploid variety Chinese Spring, but absent in variety Hope, by sequence comparison of their corresponding region. Amplifying the NUMT junctions using a set of the wheat materials including diploid, tetraploid, and hexaploid lines showed that none of the diploid wheat carried the region and only some tetraploid and hexaploid wheat were positive for the NUMT. Age estimation of the NUMT displayed the mean ages of Langdon NUMT and Chinese Spring NUMT to be 378,000 and 416,000 years ago, respectively. Reverse transcription PCR and sequencing of the nad7 gene showed 28 C → U RNA editing sites and four partial editing sites, as expected for mitochondrial DNA expression. Specific SNPs discriminated between cDNA from the nucleus and the mitochondria and suggested that the nuclear copy was not expressed. The mitochondrial DNA studied was inserted into the genome quite recently within the wheat lineage and gave rise to the non-coding nuclear nad7 gene. The NUMT segment could be lost and acquired frequently during the wheat evolution.

Construction of a Coix BAC library and isolation of the 22 kDa α-coixin gene cluster

Coix lacryma-jobi L. (Coix) is a close relative of maize and is considered a valuable genetic resource for crop improvement. Here we report the construction of the first Coix bacterial artificial chromosome (BAC) library using accession PI 324059. This BAC library contains about 230 400 clones with an average insert size of 113 kb, has low organellar DNA contamination, and provides 16.3-fold coverage of the genome. The library was stored in 12 × 96 pools that could be screened with a PCR protocol. Library screening was performed for the 22 kDa α-coixin gene family. A total of 57 positive pools were identified, and single clones were isolated from 19 of these pools. Based on DNA fingerprinting and Southern blot analysis, these 19 BAC clones form a single contig of about 340 kb in length, indicating that the 22 kDa α-coixin genes occur in a cluster. These results demonstrated the suitability of this BAC library for gene isolation and comparative genomics studies of the Coix genome.

The amplification and evolution of orthologous 22-kDa α-prolamin tandemly arrayed genes in coix, sorghum and maize genomes

Tandemly arrayed genes (TAGs) account for about one-third of the duplicated genes in eukaryotic genomes. They provide raw genetic material for biological evolution, and play important roles in genome evolution. The 22-kDa prolamin genes in cereal genomes represent typical TAG organization, and provide the good material to investigate gene amplification of TAGs in closely related grass genomes. Here, we isolated and sequenced the Coix 22-kDa prolamin (coixin) gene cluster (283 kb), and carried out a comparative analysis with orthologous 22-kDa prolamin gene clusters from maize and sorghum. The 22-kDa prolamin gene clusters descended from orthologous ancestor genes, but underwent independent gene amplification paths after the separation of these species, therefore varied dramatically in sequence and organization. Our analysis indicated that the gene amplification model of 22-kDa prolamin gene clusters can be divided into three major stages. In the first stage, rare gene duplications occurred from the ancestor gene copy accidentally. In the second stage, rounds of gene amplification occurred by unequal crossing over to form tandem gene array(s). In the third stage, gene array was further diverged by other genomic activities, such as transposon insertions, segmental rearrangements, etc. Unlike their highly conserved sequences, the amplified 22-kDa prolamin genes diverged rapidly at their expression capacities and expression levels. Such processes had no apparent correlation to age or order of amplified genes within TAG cluster, suggesting a fast evolving nature of TAGs after gene amplification. These results provided insights into the amplification and evolution of TAG families in grasses.

Ancestral grass karyotype reconstruction unravels new mechanisms of genome shuffling as a source of plant evolution

The comparison of the chromosome numbers of today's species with common reconstructed paleo-ancestors has led to intense speculation of how chromosomes have been rearranged over time in mammals. However, similar studies in plants with respect to genome evolution as well as molecular mechanisms leading to mosaic synteny blocks have been lacking due to relevant examples of evolutionary zooms from genomic sequences. Such studies require genomes of species that belong to the same family but are diverged to fall into different subfamilies. Our most important crops belong to the family of the grasses, where a number of genomes have now been sequenced. Based on detailed paleogenomics, using inference from n = 5-12 grass ancestral karyotypes (AGKs) in terms of gene content and order, we delineated sequence intervals comprising a complete set of junction break points of orthologous regions from rice, maize, sorghum, and Brachypodium genomes, representing three different subfamilies and different polyploidization events. By focusing on these sequence intervals, we could show that the chromosome number variation/reduction from the n = 12 common paleo-ancestor was driven by nonrandom centric double-strand break repair events. It appeared that the centromeric/telomeric illegitimate recombination between nonhomologous chromosomes led to nested chromosome fusions (NCFs) and synteny break points (SBPs). When intervals comprising NCFs were compared in their structure, we concluded that SBPs (1) were meiotic recombination hotspots, (2) corresponded to high sequence turnover loci through repeat invasion, and (3) might be considered as hotspots of evolutionary novelty that could act as a reservoir for producing adaptive phenotypes.

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.

Determination of the in vivo distribution of nuclear matrix attachment regions using a polymerase chain reaction-based assay in Arabidopsis thaliana

Matrix attachment regions (MARs) are the regions on genomic DNA that are attached to the nuclear matrix in eukaryotes. Previous in vitro and in silico MAR analyses have shown that MARs distribute at average intervals of about 5 kb on the Arabidopsis thaliana genome. However, the in vivo evidence for the distribution of MARs in A. thaliana is lacking. Therefore, we have used a polymerase chain reaction (PCR)-based method to investigate the in vivo locations of MARs across an 80 kb region of A. thaliana genome. This assay indicated that the average interval of MARs within this region is 4.7 kb (range 1 to 11 kb), well consistent with the previous in vitro and in silico MAR studies. This result suggests that average size of the chromatin loop in A. thaliana is smaller when compared with the other eukaryotes, in which the sizes are known to vary in the range from 9 to 100 kb. However, we found that the number of genes per chromatin loop (1-3 genes) in A. thaliana is similar to those found in other eukaryotes. Furthermore, as in animals' MARs, DNase I hypersensitive sites were also found in the MARs end-region in A. thaliana. Our results suggest that basic organization of chromatin loop in A. thaliana is similar to other eukaryotes in the view that it contains a few genes, and that the average size of chromatin loop in eukaryotes is possibly determined by genome structure, such as gene density and average gene size.

Amplification of prolamin storage protein genes in different subfamilies of the Poaceae

Prolamins are seed storage proteins in cereals and represent an important source of essential amino acids for feed and food. Genes encoding these proteins resulted from dispersed and tandem amplification. While previous studies have concentrated on protein sequences from different grass species, we now can add a new perspective to their relationships by asking how their genes are shared by ancestry and copied in different lineages of the same family of species. These differences are derived from alignment of chromosomal regions, where collinearity is used to identify prolamin genes in syntenic positions, also called orthologous gene copies. New or paralogous gene copies are inserted in tandem or new locations of the same genome. More importantly, one can detect the loss of older genes. We analyzed chromosomal intervals containing prolamin genes from rice, sorghum, wheat, barley, and Brachypodium, representing different subfamilies of the Poaceae. The Poaceae commonly known as the grasses includes three major subfamilies, the Ehrhartoideae (rice), Pooideae (wheat, barley, and Brachypodium), and Panicoideae (millets, maize, sorghum, and switchgrass). Based on chromosomal position and sequence divergence, it becomes possible to infer the order of gene amplification events. Furthermore, the loss of older genes in different subfamilies seems to permit a faster pace of divergence of paralogous genes. Change in protein structure affects their physical properties, subcellular location, and amino acid composition. On the other hand, regulatory sequence elements and corresponding transcriptional activators of new gene copies are more conserved than coding sequences, consistent with the tissue-specific expression of these genes.

Reconstruction of monocotelydoneous proto-chromosomes reveals faster evolution in plants than in animals

Paleogenomics seeks to reconstruct ancestral genomes from the genes of today's species. The characterization of paleo-duplications represented by 11,737 orthologs and 4,382 paralogs identified in five species belonging to three of the agronomically most important subfamilies of grasses, that is, Ehrhartoideae (rice) Panicoideae (sorghum, maize), and Pooideae (wheat, barley), permitted us to propose a model for an ancestral genome with a minimal size of 33.6 Mb structured in five proto-chromosomes containing at least 9,138 predicted proto-genes. It appears that only four major evolutionary shuffling events (alpha, beta, gamma, and delta) explain the divergence of these five cereal genomes during their evolution from a common paleo-ancestor. Comparative analysis of ancestral gene function with rice as a reference indicated that five categories of genes were preferentially modified during evolution. Furthermore, alignments between the five grass proto-chromosomes and the recently identified seven eudicot proto-chromosomes indicated that additional very active episodes of genome rearrangements and gene mobility occurred during angiosperm evolution. If one compares the pace of primate evolution of 90 million years (233 species) to 60 million years of the Poaceae (10,000 species), change in chromosome structure through speciation has accelerated significantly in plants.

Comparative genomic analysis of C4 photosynthetic pathway evolution in grasses

BACKGROUND

Sorghum is the first C4 plant and the second grass with a full genome sequence available. This makes it possible to perform a whole-genome-level exploration of C4 pathway evolution by comparing key photosynthetic enzyme genes in sorghum, maize (C4) and rice (C3), and to investigate a long-standing hypothesis that a reservoir of duplicated genes is a prerequisite for the evolution of C4 photosynthesis from a C3 progenitor.

RESULTS

We show that both whole-genome and individual gene duplication have contributed to the evolution of C4 photosynthesis. The C4 gene isoforms show differential duplicability, with some C4 genes being recruited from whole genome duplication duplicates by multiple modes of functional innovation. The sorghum and maize carbonic anhydrase genes display a novel mode of new gene formation, with recursive tandem duplication and gene fusion accompanied by adaptive evolution to produce C4 genes with one to three functional units. Other C4 enzymes in sorghum and maize also show evidence of adaptive evolution, though differing in level and mode. Intriguingly, a phosphoenolpyruvate carboxylase gene in the C3 plant rice has also been evolving rapidly and shows evidence of adaptive evolution, although lacking key mutations that are characteristic of C4 metabolism. We also found evidence that both gene redundancy and alternative splicing may have sheltered the evolution of new function.

CONCLUSIONS

Gene duplication followed by functional innovation is common to evolution of most but not all C4 genes. The apparently long time-lag between the availability of duplicates for recruitment into C4 and the appearance of C4 grasses, together with the heterogeneity of origins of C4 genes, suggests that there may have been a long transition process before the establishment of C4 photosynthesis.

Change of gene structure and function by non-homologous end-joining, homologous recombination, and transposition of DNA

An important objective in genome research is to relate genome structure to gene function. Sequence comparisons among orthologous and paralogous genes and their allelic variants can reveal sequences of functional significance. Here, we describe a 379-kb region on chromosome 1 of maize that enables us to reconstruct chromosome breakage, transposition, non-homologous end-joining, and homologous recombination events. Such a high-density composition of various mechanisms in a small chromosomal interval exemplifies the evolution of gene regulation and allelic diversity in general. It also illustrates the evolutionary pace of changes in plants, where many of the above mechanisms are of somatic origin. In contrast to animals, somatic alterations can easily be transmitted through meiosis because the germline in plants is contiguous to somatic tissue, permitting the recovery of such chromosomal rearrangements. The analyzed region contains the P1-wr allele, a variant of the genetically well-defined p1 gene, which encodes a Myb-like transcriptional activator in maize. The P1-wr allele consists of eleven nearly perfect P1-wr 12-kb repeats that are arranged in a tandem head-to-tail array. Although a technical challenge to sequence such a structure by shotgun sequencing, we overcame this problem by subcloning each repeat and ordering them based on nucleotide variations. These polymorphisms were also critical for recombination and expression analysis in presence and absence of the trans-acting epigenetic factor Ufo1. Interestingly, chimeras of the p1 and p2 genes, p2/p1 and p1/p2, are framing the P1-wr cluster. Reconstruction of sequence amplification steps at the p locus showed the evolution from a single Myb-homolog to the multi-gene P1-wr cluster. It also demonstrates how non-homologous end-joining can create novel gene fusions. Comparisons to orthologous regions in sorghum and rice also indicate a greater instability of the maize genome, probably due to diploidization following allotetraploidization.

The fractionated orthology of Bs2 and Rx/Gpa2 supports shared synteny of disease resistance in the Solanaceae

Comparative genomics provides a powerful tool for the identification of genes that encode traits shared between crop plants and model organisms. Pathogen resistance conferred by plant R genes of the nucleotide-binding-leucine-rich-repeat (NB-LRR) class is one such trait with great agricultural importance that occupies a critical position in understanding fundamental processes of pathogen detection and coevolution. The proposed rapid rearrangement of R genes in genome evolution would make comparative approaches tenuous. Here, we test the hypothesis that orthology is predictive of R-gene genomic location in the Solanaceae using the pepper R gene Bs2. Homologs of Bs2 were compared in terms of sequence and gene and protein architecture. Comparative mapping demonstrated that Bs2 shared macrosynteny with R genes that best fit criteria determined to be its orthologs. Analysis of the genomic sequence encompassing solanaceous R genes revealed the magnitude of transposon insertions and local duplications that resulted in the expansion of the Bs2 intron to 27 kb and the frequently detected duplications of the 5'-end of R genes. However, these duplications did not impact protein expression or function in transient assays. Taken together, our results support a conservation of synteny for NB-LRR genes and further show that their distribution in the genome has been consistent with global rearrangements.

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.

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.

Sequencing and comparative analysis of a conserved syntenic segment in the Solanaceae

Comparative genomics is a powerful tool for gaining insight into genomic function and evolution. However, in plants, sequence data that would enable detailed comparisons of both coding and noncoding regions have been limited in availability. Here we report the generation and analysis of sequences for an unduplicated conserved syntenic segment (CSS) in the genomes of five members of the agriculturally important plant family Solanaceae. This CSS includes a 105-kb region of tomato chromosome 2 and orthologous regions of the potato, eggplant, pepper, and petunia genomes. With a total neutral divergence of 0.73-0.78 substitutions/site, these sequences are similar enough that most noncoding regions can be aligned, yet divergent enough to be informative about evolutionary dynamics and selective pressures. The CSS contains 17 distinct genes with generally conserved order and orientation, but with numerous small-scale differences between species. Our analysis indicates that the last common ancestor of these species lived approximately 27-36 million years ago, that more than one-third of short genomic segments (5-15 bp) are under selection, and that more than two-thirds of selected bases fall in noncoding regions. In addition, we identify genes under positive selection and analyze hundreds of conserved noncoding elements. This analysis provides a window into 30 million years of plant evolution in the absence of polyploidization.

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.

Development of genic-microsatellite markers for sorghum staygreen QTL using a comparative genomic approach with rice

The already available comprehensive genome sequence information of model crops along with the transcriptomic resource from other crops provides an excellent opportunity for comparative genome analysis. We studied the synteny between each of the four major sorghum staygreen quantitative trait loci (QTL) regions with that in the rice genome and attempted to increase marker density around the QTL with genic-microsatellites from the sorghum transcriptomic resource using the rice genome as template. For each of the sorghum QTL regions, the reported RFLP markers were compiled, used for sequence similarity searches against the rice genome which identified syntenous regions on rice chromosome 1 for Stg1 and Stg2 QTL, on chromosome 9 for Stg3 QTL, and on chromosome 11 for Stg4 QTL. Using the Gramene genome browsing tool, 869 non-redundant sorghum expressed sequence tags (ESTs) were selected and 50 genic-microsatellites (18, 12, 15, and 5, for Stg1, Stg2, Stg3, and Stg4 QTL, respectively) could be developed. We could experimentally establish synteny of the Stg1, Stg2, Stg3, and Stg4 QTL regions with that of the rice genome by mapping ten polymorphic genic-microsatellite markers (20%) to the positions of the staygreen QTL. The simple strategy demonstrated in the present study could readily be extrapolated to other cereals of the Poaceae family. The markers developed in this study provide a basis for the isolation of genes underling these QTL using an association study or map-based gene isolation approach, and create an additional option for MAS of the staygreen trait in sorghum.

Contrasted microcolinearity and gene evolution within a homoeologous region of wheat and barley species

We study here the evolution of genes located in the same physical locus using the recently sequenced Ha locus in seven wheat genomes in diploid, tetraploid, and hexaploid species and compared them with barley and rice orthologous regions. We investigated both the conservation of microcolinearity and the molecular evolution of genes, including coding and noncoding sequences. Microcolinearity is restricted to two groups of genes (Unknown gene-2, VAMP, BGGP, Gsp-1, and Unknown gene-8 surrounded by several copies of ATPase), almost conserved in rice and barley, but in a different relative position. Highly conserved genes between wheat and rice run along with genes harboring different copy numbers and highly variable sequences between close wheat genomes. The coding sequence evolution appeared to be submitted to heterogeneous selective pressure and intronic sequences analysis revealed that the molecular clock hypothesis is violated in most cases.

Genome size and sequence composition of moso bamboo: a comparative study

Moso bamboo (Phyllostachys pubescens) is one of the world's most important bamboo species. It has the largest area of all planted bamboo--over two-thirds of the total bamboo forest area--and the highest economic value in China. Moso bamboo is a tetraploid (4x=48) and a special member of the grasses family. Although several genomes have been sequenced or are being sequenced in the grasses family, we know little about the genome of the bambusoids (bamboos). In this study, the moso bamboo genome size was estimated to be about 2034 Mb by flow cytometry (FCM), using maize (cv. B73) and rice (cv. Nipponbare) as internal references. The rice genome has been sequenced and the maize genome is being sequenced. We found that the size of the moso bamboo genome was similar to that of maize but significantly larger than that of rice. To determine whether the bamboo genome had a high proportion of repeat elements, similar to that of the maize genome, approximately 1000 genome survey sequences (GSS) were generated. Sequence analysis showed that the proportion of repeat elements was 23.3% for the bamboo genome, which is significantly lower than that of the maize genome (65.7%). The bamboo repeat elements were mainly Gypsy/DIRS1 and Ty1/Copia LTR retrotransposons (14.7%), with a few DNA transposons. However, more genomic sequences are needed to confirm the above results due to several factors, such as the limitation of our GSS data. This study is the first to investigate sequence composition of the bamboo genome. Our results are valuable for future genome research of moso and other bamboos.

Comparative analysis of the 100 kb region containing the Pi-k(h) locus between indica and japonica rice lines

We have recently cloned a pathogen inducible blast resistance gene Pi-k^h from the indica rice line Tetep using a positional cloning approach. In this study, we carried out structural organization analysis of the Pi-k^h locus in both indica and japonica rice lines. A 100 kb region containing 50 kb upstream and 50 kb downstream sequences flanking to the Pi-k^h locus was selected for the investigation. A total of 16 genes in indica and 15 genes in japonica were predicted and annotated in this region. The average GC content of indica and japonica genes in this region was 53.15% and 49.3%, respectively. Both indica and japonica sequences were polymorphic for simple sequence repeats having mono-, di-, tri-, tetra-, and pentanucleotides. Sequence analysis of the specific blast resistant Pi-k^h allele of Tetep and the susceptible Pi-k^h allele of the japonica rice line Nipponbare showed differences in the number and distribution of motifs involved in phosphorylation, resulting in the resistance phenotype in Tetep.

Large-scale, lineage-specific expansion of a bric-a-brac/tramtrack/broad complex ubiquitin-ligase gene family in rice

Selective ubiquitination of proteins is directed by diverse families of ubiquitin-protein ligases (or E3s) in plants. One important type uses Cullin-3 as a scaffold to assemble multisubunit E3 complexes containing one of a multitude of bric-a-brac/tramtrack/broad complex (BTB) proteins that function as substrate recognition factors. We previously described the 80-member BTB gene superfamily in Arabidopsis thaliana. Here, we describe the complete BTB superfamily in rice (Oryza sativa spp japonica cv Nipponbare) that contains 149 BTB domain-encoding genes and 43 putative pseudogenes. Amino acid sequence comparisons of the rice and Arabidopsis superfamilies revealed a near equal repertoire of putative substrate recognition module types. However, phylogenetic comparisons detected numerous gene duplication and/or loss events since the rice and Arabidopsis BTB lineages split, suggesting possible functional specialization within individual BTB families. In particular, a major expansion and diversification of a subset of BTB proteins containing Meprin and TRAF homology (MATH) substrate recognition sites was evident in rice and other monocots that likely occurred following the monocot/dicot split. The MATH domain of a subset appears to have evolved significantly faster than those in a smaller core subset that predates flowering plants, suggesting that the substrate recognition module in many monocot MATH-BTB E3s are diversifying to ubiquitinate a set of substrates that are themselves rapidly changing. Intriguing possibilities include pathogen proteins attempting to avoid inactivation by the monocot host.

Large-scale, lineage-specific expansion of a bric-a-brac/tramtrack/broad complex ubiquitin-ligase gene family in rice

Selective ubiquitination of proteins is directed by diverse families of ubiquitin-protein ligases (or E3s) in plants. One important type uses Cullin-3 as a scaffold to assemble multisubunit E3 complexes containing one of a multitude of bric-a-brac/tramtrack/broad complex (BTB) proteins that function as substrate recognition factors. We previously described the 80-member BTB gene superfamily in Arabidopsis thaliana. Here, we describe the complete BTB superfamily in rice (Oryza sativa spp japonica cv Nipponbare) that contains 149 BTB domain-encoding genes and 43 putative pseudogenes. Amino acid sequence comparisons of the rice and Arabidopsis superfamilies revealed a near equal repertoire of putative substrate recognition module types. However, phylogenetic comparisons detected numerous gene duplication and/or loss events since the rice and Arabidopsis BTB lineages split, suggesting possible functional specialization within individual BTB families. In particular, a major expansion and diversification of a subset of BTB proteins containing Meprin and TRAF homology (MATH) substrate recognition sites was evident in rice and other monocots that likely occurred following the monocot/dicot split. The MATH domain of a subset appears to have evolved significantly faster than those in a smaller core subset that predates flowering plants, suggesting that the substrate recognition module in many monocot MATH-BTB E3s are diversifying to ubiquitinate a set of substrates that are themselves rapidly changing. Intriguing possibilities include pathogen proteins attempting to avoid inactivation by the monocot host.

Comparative analyses between lolium/festuca introgression lines and rice reveal the major fraction of functionally annotated gene models is located in recombination-poor/very recombination-poor regions of the genome

Publication of the rice genome sequence has allowed an in-depth analysis of genome organization in a model monocot plant species. This has provided a powerful tool for genome analysis in large-genome unsequenced agriculturally important monocot species such as wheat, barley, rye, Lolium, etc. Previous data have indicated that the majority of genes in large-genome monocots are located toward the ends of chromosomes in gene-rich regions that undergo high frequencies of recombination. Here we demonstrate that a substantial component of the coding sequences in monocots is localized proximally in regions of very low and even negligible recombination frequencies. The implications of our findings are that during domestication of monocot plant species selection has concentrated on genes located in the terminal regions of chromosomes within areas of high recombination frequency. Thus a large proportion of the genetic variation available for selection of superior plant genotypes has not been exploited. In addition our findings raise the possibility of the evolutionary development of large supergene complexes that confer a selective advantage to the individual.

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.

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.

Transcriptome analysis reveals season-specific rbcS gene expression profiles in diploid perennial ryegrass (Lolium perenne L.)

Perennial ryegrass (Lolium perenne L.) is a major grass species used for forage and turf throughout the world, and gains by conventional breeding have reached a plateau. Perennial ryegrass is an outcrossing, self-incompatible diploid (2n = 2x = 14) with a relatively large genome (4067 Mbp/diploid genome; Evans, G.M., Rees, H., Snell, C.L. and Sun, S. (1972) The relation between nuclear DNA amount and the duration of the mitotic cycle. Chrom. Today, 3, 24-31). Using tissues sourced from active pastures during the peak of the autumn, winter, spring and summer seasons, we analysed the ryegrass transcriptome employing a Serial Analysis of Gene Expression (SAGE) protocol, with the dual goals of understanding the seasonal changes in perennial ryegrass gene expression and enhancing our ability to select genes for genetic manipulation. A total of 159,002 14-mer SAGE tags was sequenced and mapped to the perennial ryegrass DNA database, comprising methyl-filtered (GeneThresher) and expressed sequence tag (EST) sequences. The analysis of 14,559 unique SAGE tags, which were present more than once in our SAGE library, revealed 964, 1331, 346 and 131 exclusive transcripts to autumn, winter, spring and summer, respectively. Intriguingly, our analysis of the SAGE tags revealed season-specific expression profiles for the small subunit of ribulose-1,5-bisphosphate carboxylase (Rubisco), LprbcS. The transcript level for LprbcS was highest in spring, and then decreased gradually between summer and winter. Five different copies of LprbcS were revealed in ryegrass, with one possibly producing splice variant transcripts. Two highly expressed LprbcS genes were reported, one of which was not active in autumn. Another LprbcS gene showed an inverse expression profile to the autumn inactive LprbcS in a manner to compensate the expression level.

Comparative sequence and genetic analyses of asparagus BACs reveal no microsynteny with onion or rice

The Poales (includes the grasses) and Asparagales [includes onion (Allium cepa L.) and asparagus (Asparagus officinalis L.)] are the two most economically important monocot orders. The Poales are a member of the commelinoid monocots, a group of orders sister to the Asparagales. Comparative genomic analyses have revealed a high degree of synteny among the grasses; however, it is not known if this synteny extends to other major monocot groups such as the Asparagales. Although we previously reported no evidence for synteny at the recombinational level between onion and rice, microsynteny may exist across shorter genomic regions in the grasses and Asparagales. We sequenced nine asparagus BACs to reveal physically linked genic-like sequences and determined their most similar positions in the onion and rice genomes. Four of the asparagus BACs were selected using molecular markers tightly linked to the sex-determining M locus on chromosome 5 of asparagus. These BACs possessed only two putative coding regions and had long tracts of degenerated retroviral elements and transposons. Five asparagus BACs were selected after hybridization of three onion cDNAs that mapped to three different onion chromosomes. Genic-like sequences that were physically linked on the cDNA-selected BACs or genetically linked on the M-linked BACs showed significant similarities (e < -20) to expressed sequences on different rice chromosomes, revealing no evidence for microsynteny between asparagus and rice across these regions. Genic-like sequences that were linked in asparagus were used to identify highly similar (e < -20) expressed sequence tags (ESTs) of onion. These onion ESTs mapped to different onion chromosomes and no relationship was observed between physical or genetic linkages in asparagus and genetic linkages in onion. These results further indicate that synteny among grass genomes does not extend to a sister order in the monocots and that asparagus may not be an appropriate smaller genome model for plants in the Asparagales with enormous nuclear genomes.

Molecular genetic characterization of the Lr34/Yr18 slow rusting resistance gene region in wheat

Wheat expressed sequence tags (wESTs) were identified in a genomic interval predicted to span the Lr34/Yr18 slow rusting region on chromosome 7DS and that corresponded to genes located in the syntenic region of rice chromosome 6 (between 2.02 and 2.38 Mb). A subset of the wESTs was also used to identify corresponding bacterial artificial chromosome (BAC) clones from the diploid D genome of wheat (Aegilops tauschii). Conservation and deviation of micro-colinearity within blocks of genes were found in the D genome BACs relative to the orthologous sequences in rice. Extensive RFLP analysis using the wEST derived clones as probes on a panel of wheat genetic stocks with or without Lr34/Yr18 revealed monomorphic patterns as the norm in this region of the wheat genome. A similar pattern was observed with single nucleotide polymorphism analysis on a subset of the wEST derived clones and subclones from corresponding D genome BACs. One exception was a wEST derived clone that produced a consistent RFLP pattern that distinguished the Lr34/Yr18 genetic stocks and well-established cultivars known either to possess or lack Lr34/Yr18. Conversion of the RFLP to a codominant sequence tagged site (csLV34) revealed a bi-allelic locus, where a variant size of 79 bp insertion in an intron sequence was associated with lines or cultivars that lacked Lr34/Yr18. This association with Lr34/Yr18 was validated in wheat cultivars from diverse backgrounds. Genetic linkage between csLV34 and Lr34/Yr18 was estimated at 0.4 cM.

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.

Uneven chromosome contraction and expansion in the maize genome

Maize (Zea mays or corn), both a major food source and an important cytogenetic model, evolved from a tetraploid that arose about 4.8 million years ago (Mya). As a result, maize has extensive duplicated regions within its genome. We have sequenced the two copies of one such region, generating 7.8 Mb of sequence spanning 17.4 cM of the short arm of chromosome 1 and 6.6 Mb (25.6 cM) from the long arm of chromosome 9. Rice, which did not undergo a similar whole genome duplication event, has only one orthologous region (4.9 Mb) on the short arm of chromosome 3, and can be used as reference for the maize homoeologous regions. Alignment of the three regions allowed identification of syntenic blocks, and indicated that the maize regions have undergone differential contraction in genic and intergenic regions and expansion by the insertion of retrotransposable elements. Approximately 9% of the predicted genes in each duplicated region are completely missing in the rice genome, and almost 20% have moved to other genomic locations. Predicted genes within these regions tend to be larger in maize than in rice, primarily because of the presence of predicted genes in maize with larger introns. Interestingly, the general gene methylation patterns in the maize homoeologous regions do not appear to have changed with contraction or expansion of their chromosomes. In addition, no differences in methylation of single genes and tandemly repeated gene copies have been detected. These results, therefore, provide new insights into the diploidization of polyploid species.

Single-copy genes define a conserved order between rice and wheat for understanding differences caused by duplication, deletion, and transposition of genes

The high-quality rice genome sequence is serving as a reference for comparative genome analysis in crop plants, especially cereals. However, early comparisons with bread wheat showed complex patterns of conserved synteny (gene content) and colinearity (gene order). Here, we show the presence of ancient duplicated segments in the progenitor of wheat, which were first identified in the rice genome. We also show that single-copy (SC) rice genes, those representing unique matches with wheat expressed sequence tag (EST) unigene contigs in the whole rice genome, show more than twice the proportion of genes mapping to syntenic wheat chromosome as compared to the multicopy (MC) or duplicated rice genes. While 58.7% of the 1,244 mapped SC rice genes were located in single syntenic wheat chromosome groups, the remaining 41.3% were distributed randomly to the other six non-syntenic wheat groups. This could only be explained by a background dispersal of genes in the genome through transposition or other unknown mechanism. The breakdown of rice-wheat synteny due to such transpositions was much greater near the wheat centromeres. Furthermore, the SC rice genes revealed a conserved primordial gene order that gives clues to the origin of rice and wheat chromosomes from a common ancestor through polyploidy, aneuploidy, centromeric fusions, and translocations. Apart from the bin-mapped wheat EST contigs, we also compared 56,298 predicted rice genes with 39,813 wheat EST contigs assembled from 409,765 EST sequences and identified 7,241 SC rice gene homologs of wheat. Based on the conserved colinearity of 1,063 mapped SC rice genes across the bins of individual wheat chromosomes, we predicted the wheat bin location of 6,178 unmapped SC rice gene homologs and validated the location of 213 of these in the telomeric bins of 21 wheat chromosomes with 35.4% initial success. This opens up the possibility of directed mapping of a large number of conserved SC rice gene homologs in wheat. Overall, only 46.4% of these SC genes code for proteins with known functional domains; the remaining 53.6% have unknown function, and hence, represent an important, but yet, under explored category of genes.

Organization and variability of the maize genome

With a size approximating that of the human genome, the maize genome is about to become the largest plant genome yet sequenced. Contributing to that size are a whole-genome duplication event and a retrotransposition explosion that produced a large amount of repetitive DNA. This DNA is greatly under-represented in cDNA collections, so analysis of the maize transcriptome has been an expedient way of assessing the gene content of maize. Over 2 million maize cDNA sequences are now available, making maize the third most widely studied organism, behind mouse and man. To date, the sequencing of large-sized DNA clones has been largely driven by the genetic interests of different investigators. The recent construction of a physical map that is anchored to the genetic map will aid immensely in the maize genome-sequencing effort. However, studies showing that the repetitive DNA component is highly polymorphic among maize inbred lines point to the need to sample vertically a few specific regions of the genome to evaluate the extent and importance of this variability.

Molecular cytogenetics and DNA sequence analysis of an apomixis-linked BAC in Paspalum simplex reveal a non pericentromere location and partial microcolinearity with rice

Apomixis in plants is a form of clonal reproduction through seeds. A BAC clone linked to apomictic reproduction in Paspalum simplex was used to locate the apomixis locus on meiotic chromosome preparations. Fluorescent in situ hybridisation revealed the existence of a single locus embedded in a heterochromatin-poor region not adjacent to the centromere. We report here for the first time information regarding the sequencing of a large DNA clone from the apomixis locus. The presence of two genes whose rice homologs were mapped on the telomeric part of the long arm of rice chromosome 12 confirmed the strong synteny between the apomixis locus of P. simplex with the related area of the rice genome at the map level. Comparative analysis of this region with rice as representative of a sexual species revealed large-scale rearrangements due to transposable elements and small-scale rearrangements due to deletions and single point mutations. Both types of rearrangements induced the loss of coding capacity of large portions of the "apomictic" genes compared to their rice homologs. Our results are discussed in relation to the use of rice genome data for positional cloning of apomixis genes and to the possible role of rearranged supernumerary genes in the apomictic process of P. simplex.

The maize genome as a model for efficient sequence analysis of large plant genomes

The genomes of flowering plants vary in size from about 0.1 to over 100 gigabase pairs (Gbp), mostly because of polyploidy and variation in the abundance of repetitive elements in intergenic regions. High-quality sequences of the relatively small genomes of Arabidopsis (0.14 Gbp) and rice (0.4 Gbp) have now been largely completed. The sequencing of plant genomes that have a more representative size (the mean for flowering plant genomes is 5.6 Gbp) has been seen as a daunting task, partly because of their size and partly because of the numerous highly conserved repeats. Nevertheless, creative strategies and powerful new tools have been generated recently in the plant genetics community, so that sequencing large plant genomes is now a realistic possibility. Maize (2.4-2.7 Gbp) will be the first gigabase-size plant genome to be sequenced using these novel approaches. Pilot studies on maize indicate that the new gene-enrichment, gene-finishing and gene-orientation technologies are efficient, robust and comprehensive. These strategies will succeed in sequencing the gene-space of large genome plants, and in locating all of these genes and adjacent sequences on the genetic and physical maps.

ramosa2 encodes a LATERAL ORGAN BOUNDARY domain protein that determines the fate of stem cells in branch meristems of maize

Genetic control of grass inflorescence architecture is critical given that cereal seeds provide most of the world's food. Seeds are borne on axillary branches, which arise from groups of stem cells in axils of leaves and whose branching patterns dictate most of the variation in plant form. Normal maize (Zea mays) ears are unbranched, and tassels have long branches only at their base. The ramosa2 (ra2) mutant of maize has increased branching with short branches replaced by long, indeterminate ones. ra2 was cloned by chromosome walking and shown to encode a LATERAL ORGAN BOUNDARY domain transcription factor. ra2 is transiently expressed in a group of cells that predicts the position of axillary meristem formation in inflorescences. Expression in different mutant backgrounds places ra2 upstream of other genes that regulate branch formation. The early expression of ra2 suggests that it functions in the patterning of stem cells in axillary meristems. Alignment of ra2-like sequences reveals a grass-specific domain in the C terminus that is not found in Arabidopsis thaliana. The ra2-dm allele suggests this domain is required for transcriptional activation of ra1. The ra2 expression pattern is conserved in rice (Oryza sativa), barley (Hordeum vulgare), sorghum (Sorghum bicolor), and maize, suggesting that ra2 is critical for shaping the initial steps of grass inflorescence architecture.

Structure and architecture of the maize genome

Maize (Zea mays or corn) plays many varied and important roles in society. It is not only an important experimental model plant, but also a major livestock feed crop and a significant source of industrial products such as sweeteners and ethanol. In this study we report the systematic analysis of contiguous sequences of the maize genome. We selected 100 random regions averaging 144 kb in size, representing about 0.6% of the genome, and generated a high-quality dataset for sequence analysis. This sampling contains 330 annotated genes, 91% of which are supported by expressed sequence tag data from maize and other cereal species. Genes averaged 4 kb in size with five exons, although the largest was over 59 kb with 31 exons. Gene density varied over a wide range from 0.5 to 10.7 genes per 100 kb and genes did not appear to cluster significantly. The total repetitive element content we observed (66%) was slightly higher than previous whole-genome estimates (58%-63%) and consisted almost exclusively of retroelements. The vast majority of genes can be aligned to at least one sequence read derived from gene-enrichment procedures, but only about 30% are fully covered. Our results indicate that much of the increase in genome size of maize relative to rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana) is attributable to an increase in number of both repetitive elements and genes.

Complex microcolinearity among wheat, rice, and barley revealed by fine mapping of the genomic region harboring a major QTL for resistance to Fusarium head blight in wheat

A major quantitative trait locus (QTL), Qfhs.ndsu-3BS, for resistance to Fusarium head blight (FHB) in wheat has been identified and verified by several research groups. The objectives of this study were to construct a fine genetic map of this QTL region and to examine microcolinearity in the QTL region among wheat, rice, and barley. Two simple sequence repeat (SSR) markers (Xgwm533 and Xgwm493) flanking this QTL were used to screen for recombinants in a population of 3,156 plants derived from a single F(7) plant heterozygous for the Qfhs.ndsu-3BS region. A total of 382 recombinants were identified, and they were genotyped with two more SSR markers and eight sequence-tagged site (STS) markers. A fine genetic map of the Qfhs.ndsu-3BS region was constructed and spanned 6.3 cM. Based on replicated evaluations of homozygous recombinant lines for Type II FHB resistance, Qfhs.ndsu-3BS, redesignated as Fhb1, was placed into a 1.2-cM marker interval flanked by STS3B-189 and STS3B-206. Primers of STS markers were designed from wheat expressed sequence tags homologous to each of six barley genes expected to be located near this QTL region. A comparison of the wheat fine genetic map and physical maps of rice and barley revealed inversions and insertions/deletions. This suggests a complex microcolinearity among wheat, rice, and barley in this QTL region.

High-resolution mapping of the barley leaf rust resistance gene Rph5 using barley expressed sequence tags (ESTs) and synteny with rice

The rapidly growing expressed sequence tag (EST) resources of species representing the Poacea family and availability of comprehensive sequence information for the rice (Oryza sativa) genome create an excellent opportunity for comparative genome analysis. Extensive synteny between rice chromosome 1 and barley (Hordeum vulgare L.) chromosome 3 has proven extremely useful for saturation mapping of chromosomal regions containing target genes of large-genome barley with conserved orthologous genes from the syntenic regions of the rice genome. Rph5 is a gene conferring resistance to the barley leaf rust pathogen Puccinia hordei. It was mapped to chromosome 3HS, which is syntenic with rice chromosome 1S. The objective of this study was to increase marker density within the sub-centimorgan region around Rph5, using sequence-tagged site (STS) markers that were developed based on barley ESTs syntenic to the phage (P1)-derived artificial chromosome (PAC) clones comprising the distal region of rice chromosome 1S. Five rice PAC clones were used as queries in a blastn search to screen 375,187 barley ESTs. Ninety-four non-redundant EST sequences were identified from the EST database and used as templates to design 174 pairs of primer combinations. As a result, 9 barley EST-based STS markers were incorporated into the 'Bowman' x 'Magnif 102' high-resolution map of the Rph5 region. More importantly, six markers, including five EST-derived STS sequences, were found to co-segregate with Rph5. The results of this study demonstrate the usefulness of rice genomic resources for efficient deployment of barley ESTs for marker saturation of targeted barley genomic regions.

Genetic mapping of expressed sequences in onion and in silico comparisons with rice show scant colinearity

The Poales (which include the grasses) and Asparagales [which include onion (Allium cepa L.) and other Allium species] are the two most economically important monocot orders. Enormous genomic resources have been developed for the grasses; however, their applicability to other major monocot groups, such as the Asparagales, is unclear. Expressed sequence tags (ESTs) from onion that showed significant similarities (80% similarity over at least 70% of the sequence) to single positions in the rice genome were selected. One hundred new genetic markers developed from these ESTs were added to the intraspecific map derived from the BYG15-23xAC43 segregating family, producing 14 linkage groups encompassing 1,907 cM at LOD 4. Onion linkage groups were assigned to chromosomes using alien addition lines of Allium fistulosum L. carrying single onion chromosomes. Visual comparisons of genetic linkage in onion with physical linkage in rice revealed scant colinearity; however, short regions of colinearity could be identified. Our results demonstrate that the grasses may not be appropriate genomic models for other major monocot groups such as the Asparagales; this will make it necessary to develop genomic resources for these important plants.

Chromosome evolution in eukaryotes: a multi-kingdom perspective

In eukaryotes, chromosomal rearrangements, such as inversions, translocations and duplications, are common and range from part of a gene to hundreds of genes. Lineage-specific patterns are also seen: translocations are rare in dipteran flies, and angiosperm genomes seem prone to polyploidization. In most eukaryotes, there is a strong association between rearrangement breakpoints and repeat sequences. Current data suggest that some repeats promoted rearrangements via non-allelic homologous recombination, for others the association might not be causal but reflects the instability of particular genomic regions. Rearrangement polymorphisms in eukaryotes are correlated with phenotypic differences, so are thought to confer varying fitness in different habitats. Some seem to be under positive selection because they either trap favorable allele combinations together or alter the expression of nearby genes. There is little evidence that chromosomal rearrangements cause speciation, but they probably intensify reproductive isolation between species that have formed by another route.

Fertility restorer locus Rf1 [corrected] of sorghum (Sorghum bicolor L.) encodes a pentatricopeptide repeat protein not present in the colinear region of rice chromosome 12

With an aim to clone the sorghum fertility restorer gene Rf1, a high-resolution genetic and physical map of the locus was constructed. The Rf1 locus was resolved to a 32-kb region spanning four open reading frames: a plasma membrane Ca(2+)-ATPase, a cyclin D-1, an unknown protein, and a pentatricopeptide repeat (PPR13) gene family member. An approximately 19-kb region spanning the cyclin D-1 and unknown protein genes was completely conserved between sterile and fertile plants as was the sequence spanning the coding region of the Ca(2+)-ATPase. In contrast, 19 sequence polymorphisms were located in an approximately 7-kb region spanning PPR13, and all markers cosegregated with the fertility restoration phenotype. PPR13 was predicted to encode a mitochondrial-targeted protein containing a single exon with 14 PPR repeats, and the protein is classified as an E-type PPR subfamily member. To permit sequence-based comparison of the sorghum and rice genomes in the Rf1 region, 0.53 Mb of sorghum chromosome 8 was sequenced and compared to the colinear region of rice chromosome 12. Genome comparison revealed a mosaic pattern of colinearity with an approximately 275-kb gene-poor region with little gene conservation and an adjacent, approximately 245-kb gene-rice region that is more highly conserved between rice and sorghum. Despite being located in a region of high gene conservation, sorghum PPR13 was not located in a colinear position on rice chromosome 12. The present results suggest that sorghum PPR13 represents a potential candidate for the sorghum Rf1 gene, and its presence in the sorghum genome indicates a single gene transposition event subsequent to the divergence of rice and sorghum ancestors.

Fertility restorer locus Rf1 [corrected] of sorghum (Sorghum bicolor L.) encodes a pentatricopeptide repeat protein not present in the colinear region of rice chromosome 12

With an aim to clone the sorghum fertility restorer gene Rf1, a high-resolution genetic and physical map of the locus was constructed. The Rf1 locus was resolved to a 32-kb region spanning four open reading frames: a plasma membrane Ca(2+)-ATPase, a cyclin D-1, an unknown protein, and a pentatricopeptide repeat (PPR13) gene family member. An approximately 19-kb region spanning the cyclin D-1 and unknown protein genes was completely conserved between sterile and fertile plants as was the sequence spanning the coding region of the Ca(2+)-ATPase. In contrast, 19 sequence polymorphisms were located in an approximately 7-kb region spanning PPR13, and all markers cosegregated with the fertility restoration phenotype. PPR13 was predicted to encode a mitochondrial-targeted protein containing a single exon with 14 PPR repeats, and the protein is classified as an E-type PPR subfamily member. To permit sequence-based comparison of the sorghum and rice genomes in the Rf1 region, 0.53 Mb of sorghum chromosome 8 was sequenced and compared to the colinear region of rice chromosome 12. Genome comparison revealed a mosaic pattern of colinearity with an approximately 275-kb gene-poor region with little gene conservation and an adjacent, approximately 245-kb gene-rice region that is more highly conserved between rice and sorghum. Despite being located in a region of high gene conservation, sorghum PPR13 was not located in a colinear position on rice chromosome 12. The present results suggest that sorghum PPR13 represents a potential candidate for the sorghum Rf1 gene, and its presence in the sorghum genome indicates a single gene transposition event subsequent to the divergence of rice and sorghum ancestors.