Structure and evolution of cereal genomes

Alternative splicing of flowering time gene FT is associated with halving of time to flowering in coconut

Coconut palm has two distinct types-"tall" and "dwarf"-which differ morphologically. Tall coconut varieties need 8-10 years to start flowering, while dwarf coconut varieties only require 3-5 years. We compared seedling and reproductive stage transcriptomes for both coconut types to determine potential molecular mechanisms underlying control of flowering time in coconut. Several key genes in the photoperiod pathway were differentially expressed between seedling and reproductive leaf samples in both tall and dwarf coconut. These genes included suppressor of overexpression of constans (SOC1), flowering locus T (FT), and Apetala 1 (AP1). Alternative splicing analysis of genes in the photoperiod pathway further revealed that the FT gene produces different transcripts in tall compared to dwarf coconut. The shorter alternative splice variant of FT [which included a 6 bp deletion, alternative 3' splicing sites (A3SS)] was found to be exclusively present in dwarf coconut varieties but absent in most tall coconut varieties. Our results provide a valuable information resource as well as suggesting a probable mechanism for differentiation of flowering time onset in coconut, providing a target for future breeding work in accelerating time to flowering in this crop species.

DArT Markers Effectively Target Gene Space in the Rye Genome

Large genome size and complexity hamper considerably the genomics research in relevant species. Rye (Secale cereale L.) has one of the largest genomes among cereal crops and repetitive sequences account for over 90% of its length. Diversity Arrays Technology is a high-throughput genotyping method, in which a preferential sampling of gene-rich regions is achieved through the use of methylation sensitive restriction enzymes. We obtained sequences of 6,177 rye DArT markers and following a redundancy analysis assembled them into 3,737 non-redundant sequences, which were then used in homology searches against five Pooideae sequence sets. In total 515 DArT sequences could be incorporated into publicly available rye genome zippers providing a starting point for the integration of DArT- and transcript-based genomics resources in rye. Using Blast2Go pipeline we attributed putative gene functions to 1101 (29.4%) of the non-redundant DArT marker sequences, including 132 sequences with putative disease resistance-related functions, which were found to be preferentially located in the 4RL and 6RL chromosomes. Comparative analysis based on the DArT sequences revealed obvious inconsistencies between two recently published high density consensus maps of rye. Furthermore we demonstrated that DArT marker sequences can be a source of SSR polymorphisms. Obtained data demonstrate that DArT markers effectively target gene space in the large, complex, and repetitive rye genome. Through the annotation of putative gene functions and the alignment of DArT sequences relative to reference genomes we obtained information, that will complement the results of the studies, where DArT genotyping was deployed, by simplifying the gene ontology and microcolinearity based identification of candidate genes.

The DYT1-interacting proteins bHLH010, bHLH089 and bHLH091 are redundantly required for Arabidopsis anther development and transcriptome

The anther is the male reproductive organ of flowering plants, and the Arabidopsis bHLH transcription factors encoded by DYSFUNCTIONAL TAPETUM1 (DYT1) and ABORTED MICROSPORE (AMS) are required for control of the complex transcriptional networks regulating anther development. Knowledge of the mechanisms by which the bHLH proteins affect this diverse gene expression is quite limited. We examine here three recently duplicated Arabidopsis bHLH genes, bHLH010, bHLH089 and bHLH091, using evolutionary, genetic, morphological and transcriptomic approaches, and uncover their redundant functions in anther development. These three genes are relatively highly expressed in the tapetum of the Arabidopsis anther; single mutants at each of the bHLH010, bHLH089 and bHLH091 loci are developmentally normal, but the various double and triple combinations progressively exhibit increasingly defective anther phenotypes (abnormal tapetum morphology, delayed callose degeneration, and aborted pollen development), indicating their redundant functions in male fertility. Further transcriptomic and molecular analyses suggest that these three proteins act slightly later than DYT1, and also form protein complexes with DYT1, subsequently affecting the correct expression of many DYT1 target genes in the anther development transcriptional network. This study demonstrated that bHLH010, bHLH089 and bHLH091 together are important for the normal transcriptome of the developing Arabidopsis anther, possibly by forming a feed-forward loop with DYT1.

Identification of HSP20 gene family in wheat and barley and their differential expression profiling under heat stress

Small heat shock proteins (sHSPs) are chaperones that play an important role in various developmental, biotic and abiotic stresses. The sHSP family possess a conserved domain of approximately 80 to 100 amino acids called alpha-crystalline domain (ACD), flanked by N- and C-terminal regions. Search for complete proteomes and expressed sequenced tag (EST) database of wheat and barley using Hidden Markov Model and BLAST algorithm was conducted. Here, we report genome-wide identification and characterization of 27 newly TaHSP20 candidate genes in wheat and 13 HvHSP20 in barley, describing structures, phylogenetic relationships, conserved protein motifs, and expression patterns. The structural analysis highlights that this gene family possesses a conserved ACD region at the C-terminal. Detailed pattern analysis of HSP20 revealed presence of P-G doublet and I/V/L-X-I/V/L motif that helps in oligomerization. Identification of conserved motif sequences of wheat and barley HSP20 strongly supported their identity as sHSP families. This study illustrates for the first time 3D model prediction of full-length wheat HSP20 (TaHSP20) protein and ACD region. Digital expression analysis was also carried out in order to reveal a widespread distribution of the sHSP family genes at various developmental stages of wheat and barley. In addition, five selected transcripts of both wheat and barley were validated for their expression profile under 35 °C and 42 °C heat stress conditions. Results indicate up-regulation of all the transcripts under heat stress condition except TaCBM38894 candidate, which showed down-regulation in wheat.

Comparative analysis of syntenic genes in grass genomes reveals accelerated rates of gene structure and coding sequence evolution in polyploid wheat

Cycles of whole-genome duplication (WGD) and diploidization are hallmarks of eukaryotic genome evolution and speciation. Polyploid wheat (Triticum aestivum) has had a massive increase in genome size largely due to recent WGDs. How these processes may impact the dynamics of gene evolution was studied by comparing the patterns of gene structure changes, alternative splicing (AS), and codon substitution rates among wheat and model grass genomes. In orthologous gene sets, significantly more acquired and lost exonic sequences were detected in wheat than in model grasses. In wheat, 35% of these gene structure rearrangements resulted in frame-shift mutations and premature termination codons. An increased codon mutation rate in the wheat lineage compared with Brachypodium distachyon was found for 17% of orthologs. The discovery of premature termination codons in 38% of expressed genes was consistent with ongoing pseudogenization of the wheat genome. The rates of AS within the individual wheat subgenomes (21%-25%) were similar to diploid plants. However, we uncovered a high level of AS pattern divergence between the duplicated homeologous copies of genes. Our results are consistent with the accelerated accumulation of AS isoforms, nonsynonymous mutations, and gene structure rearrangements in the wheat lineage, likely due to genetic redundancy created by WGDs. Whereas these processes mostly contribute to the degeneration of a duplicated genome and its diploidization, they have the potential to facilitate the origin of new functional variations, which, upon selection in the evolutionary lineage, may play an important role in the origin of novel traits.

Bioenergy grass feedstock: current options and prospects for trait improvement using emerging genetic, genomic, and systems biology toolkits

For lignocellulosic bioenergy to become a viable alternative to traditional energy production methods, rapid increases in conversion efficiency and biomass yield must be achieved. Increased productivity in bioenergy production can be achieved through concomitant gains in processing efficiency as well as genetic improvement of feedstock that have the potential for bioenergy production at an industrial scale. The purpose of this review is to explore the genetic and genomic resource landscape for the improvement of a specific bioenergy feedstock group, the C4 bioenergy grasses. First, bioenergy grass feedstock traits relevant to biochemical conversion are examined. Then we outline genetic resources available bioenergy grasses for mapping bioenergy traits to DNA markers and genes. This is followed by a discussion of genomic tools and how they can be applied to understanding bioenergy grass feedstock trait genetic mechanisms leading to further improvement opportunities.

A major locus qS12, located in a duplicated segment of chromosome 12, causes spikelet sterility in an indica-japonica rice hybrid

Chromosome segment duplications are integral in genome evolution by providing a source for the origin of new genes. In the rice genome, besides an ancient polyploidy event known in the rice common ancestor, it had been identified that there was a special segmental duplication involving chromosomes 11 and 12, but the biological role of this duplication remains unknown. In this study, by using a set of chromosome segment substitution lines (CSSLs) and near isogenic lines (NILs) derived from the indica cultivar 9311 and japonica cultivar Nipponbare, a major QTL (qS12) resulting in hybrid male sterility was mapped within ~400 kb region adjacent to the special duplicated segment on the short arm of chromosome 12. Compared to the japonica cultivar Nipponbare, the two sides of the qS12 candidate region were inverted in the indica cultivar 9311. Among 47 of the 111 rice genotypes evaluated by molecular markers, the inverted sides were detected, and found completely homologous to indica cultivar 9311. These results suggested that the two inverted sides protect the sequence in the qS12 regions from recombination. On the short-arm of chromosome 12, two QTLs S-e and S25, in addition to qS12, were previously detected as a distinct segregation distortion and pollen semi-sterility loci. We propose these three hybrid sterility loci are the same locus, and the duplicated segment on chromosome 12 may play a prominent role in diversification, i.e., sub-speciation of cultivated rice.

Genetic dissection of the biotic stress response using a genome-scale gene network for rice

Rice is a staple food for one-half the world's population and a model for other monocotyledonous species. Thus, efficient approaches for identifying key genes controlling simple or complex traits in rice have important biological, agricultural, and economic consequences. Here, we report on the construction of RiceNet, an experimentally tested genome-scale gene network for a monocotyledonous species. Many different datasets, derived from five different organisms including plants, animals, yeast, and humans, were evaluated, and 24 of the most useful were integrated into a statistical framework that allowed for the prediction of functional linkages between pairs of genes. Genes could be linked to traits by using guilt-by-association, predicting gene attributes on the basis of network neighbors. We applied RiceNet to an important agronomic trait, the biotic stress response. Using network guilt-by-association followed by focused protein-protein interaction assays, we identified and validated, in planta, two positive regulators, LOC_Os01g70580 (now Regulator of XA21; ROX1) and LOC_Os02g21510 (ROX2), and one negative regulator, LOC_Os06g12530 (ROX3). These proteins control resistance mediated by rice XA21, a pattern recognition receptor. We also showed that RiceNet can accurately predict gene function in another major monocotyledonous crop species, maize. RiceNet thus enables the identification of genes regulating important crop traits, facilitating engineering of pathways critical to crop productivity.

Genome-wide identification and expression profiling of auxin response factor (ARF) gene family in maize

Background

Auxin signaling is vital for plant growth and development, and plays important role in apical dominance, tropic response, lateral root formation, vascular differentiation, embryo patterning and shoot elongation. Auxin Response Factors (ARFs) are the transcription factors that regulate the expression of auxin responsive genes. The ARF genes are represented by a large multigene family in plants. The first draft of full maize genome assembly has recently been released, however, to our knowledge, the ARF gene family from maize (ZmARF genes) has not been characterized in detail.

Results

In this study, 31 maize (Zea mays L.) genes that encode ARF proteins were identified in maize genome. It was shown that maize ARF genes fall into related sister pairs and chromosomal mapping revealed that duplication of ZmARFs was associated with the chromosomal block duplications. As expected, duplication of some ZmARFs showed a conserved intron/exon structure, whereas some others were more divergent, suggesting the possibility of functional diversification for these genes. Out of these 31 ZmARF genes, 14 possess auxin-responsive element in their promoter region, among which 7 appear to show small or negligible response to exogenous auxin. The 18 ZmARF genes were predicted to be the potential targets of small RNAs. Transgenic analysis revealed that increased miR167 level could cause degradation of transcripts of six potential targets (ZmARF3, 9, 16, 18, 22 and 30). The expressions of maize ARF genes are responsive to exogenous auxin treatment. Dynamic expression patterns of ZmARF genes were observed in different stages of embryo development.

Conclusions

Maize ARF gene family is expanded (31 genes) as compared to Arabidopsis (23 genes) and rice (25 genes). The expression of these genes in maize is regulated by auxin and small RNAs. Dynamic expression patterns of ZmARF genes in embryo at different stages were detected which suggest that maize ARF genes may be involved in seed development and germination.

Seventy million years of concerted evolution of a homoeologous chromosome pair, in parallel, in major Poaceae lineages

Whole genome duplication ~70 million years ago provided raw material for Poaceae (grass) diversification. Comparison of rice (Oryza sativa), sorghum (Sorghum bicolor), maize (Zea mays), and Brachypodium distachyon genomes revealed that one paleo-duplicated chromosome pair has experienced very different evolution than all the others. For tens of millions of years, the two chromosomes have experienced illegitimate recombination that has been temporally restricted in a stepwise manner, producing structural stratification in the chromosomes. These strata formed independently in different grass lineages, with their similarities (low sequence divergence between paleo-duplicated genes) preserved in parallel for millions of years since the divergence of these lineages. The pericentromeric region of this homeologous chromosome pair accounts for two-thirds of the gene content differences between the modern chromosomes. Both intriguing and perplexing is a distal chromosomal region with the greatest DNA similarity between surviving duplicated genes but also with the highest concentration of lineage-specific gene pairs found anywhere in these genomes and with a significantly elevated gene evolutionary rate. Intragenomic similarity near this chromosomal terminus may be important in hom(e)ologous chromosome pairing. Chromosome structural stratification, together with enrichment of autoimmune response-related (nucleotide binding site-leucine-rich repeat) genes and accelerated DNA rearrangement and gene loss, confer a striking resemblance of this grass chromosome pair to the sex chromosomes of other taxa.

The evolutionary history of PDR in Brachypodium distachyon polyploids

The ATP-binding cassette transporter genes include the pleiotropic drug resistance (PDR) family found only in fungi and plants. These transporters transport toxic compounds across biological membranes. Here, we investigated the evolution of the PDR1 gene in Brachypodium distachyon, a widely distributed temperate grass species that belongs to the Poaceae (Gramineae) family, which also contains the domesticated cereal crops. Because this species has multiple ploidy levels, investigating PDR1 evolution in B. distachyon will offer insights into the formation and evolution of polyploidy. From 23 B. distachyon ecotypes, 39 PDR1 homologs were identified. All ecotypes had either one or two PDR1 copies. Based on restriction site analysis, the PDR1 homologs were classified as E or H type. All but one diploid and tetraploid ecotypes had only a single H type PDR1. All but one hexaploid ecotypes had both an E and a H type PDR1. Phylogenetic analysis revealed that each type formed a well-supported cluster. The two PDR1 types appeared to evolve differently. These different evolutionary patterns could indicate a difference in age between the two types or might indicate different mutation rates or selection pressures on the two types. The phylogenetic analysis also revealed that the hexaploid ecotypes shared a genomic origin for their E type PDR1, but there were multiple origins for hexaploid H type PDR1 homologs. Overall, the results suggest that tetraploid and hexaploid might be misnomers in B. distachyon and suggest a complex polyploidization history during B. distachyon evolution.

Use of chromosome walking in discovery of single-nucleotide polymorphism in noncoding regions of a candidate actin gene in Pinus radiata

Untranslated regions (UTRs) of eukaryotic mRNAs play crucial roles in post-transcriptional regulation of gene expression via the modulation of nucleocytoplasmic mRNA transport, translation efficiency, subcellular localization, and message stability. Single-nucleotide polymorphisms (SNPs) in UTRs of a candidate gene may also change the post-transcriptional regulation of a gene or function by nucleotide mutation. For species that have not been entirely sequenced genomically, new methods need to be devised to discover SNPs in noncoding regions of candidate genes. In this study, based on the expressed sequence tag (EST) of Pinus radiata (Monterey pine), we obtained all the sequences of UTRs of the actin gene by using a chromosome walking method. We also detected all the SNPs in and around the coding region of the actin gene. In this way, the full genomic sequence (2154 bp) of the actin gene was identified, including the 5’UTR, introns, the coding sequence, and the 3’UTR. PCR amplification and DNA fragment sequencing from 200 unrelated P. radiata trees revealed a total of 21 SNPs in the actin gene, of which 3 were located in the 5’UTR, 3 in the introns, 10 in the coding sequence, and 5 in the 3’UTR. We show that chromosome walking can be used for obtaining the sequence of UTRs, and then, based on this sequence, to discover SNPs in the noncoding regions of candidate genes from this species without an entire genomic sequence.

The evolution of Brassica napus FLOWERING LOCUS T paralogues in the context of inverted chromosomal duplication blocks

Background

The gene FLOWERING LOCUS T (FT) and its orthologues play a central role in the integration of flowering signals within Arabidopsis and other diverse species. Multiple copies of FT, with different cis-intronic sequence, exist and appear to operate harmoniously within polyploid crop species such as Brassica napus (AACC), a member of the same plant family as Arabidopsis.

Results

We have identified six BnFT paralogues from the genome of B. napus and mapped them to six distinct regions, each of which is homologous to a common ancestral block (E) of Arabidopsis chromosome 1. Four of the six regions were present within inverted duplicated regions of chromosomes A7 and C6. The coding sequences of BnFT paralogues showed 92-99% identities to each other and 85-87% identity with that of Arabidopsis. However, two of the paralogues on chromosomes A2 and C2, BnA2.FT and BnC2.FT, were found to lack the distinctive CArG box that is located within intron 1 that has been shown in Arabidopsis to be the binding site for theFLC protein. Three BnFT paralogues (BnA2.FT, BnC6.FT.a and BnC6.FT.b) were associated with two major QTL clusters for flowering time. One of the QTLs encompassing two BnFT paralogues (BnC6.FT.a and BnC6.FT.b) on chromosome C6 was resolved further using near isogenic lines, specific alleles of which were both shown to promote flowering. Association analysis of the three BnFT paralogues across 55 cultivars of B. napus showed that the alleles detected in the original parents of the mapping population used to detect QTL (NY7 and Tapidor) were ubiquitous amongst spring and winter type cultivars of rapeseed. It was inferred that the ancestral FT homologues in Brassica evolved from two distinct copies, one of which was duplicated along with inversion of the associated chromosomal segment prior to the divergence of B. rapa (AA) and B. oleracea (CC). At least ten such inverted duplicated blocks (IDBs) were identified covering a quarter of the whole B. napus genome.

Conclusion

Six orthologues of Arabidopsis FT were identified and mapped in the genome of B. napus which sheds new light on the evolution of paralogues in polyploidy species. The allelic variation of BnFT paralogues results in functional differences affecting flowering time between winter and spring type cultivars of oilseed Brassica. The prevalent inverted duplicated blocks, two of which were located by four of the six BnFT paralogues, contributed to gene duplications and might represent predominant pathway of evolution in Brassica.

Genome-wide comparative analysis of the Brassica rapa gene space reveals genome shrinkage and differential loss of duplicated genes after whole genome triplication

Euchromatic regions of the Brassica rapa genome were sequenced and mapped onto the corresponding regions in the Arabidopsis thaliana genome.

Background

Brassica rapa is one of the most economically important vegetable crops worldwide. Owing to its agronomic importance and phylogenetic position, B. rapa provides a crucial reference to understand polyploidy-related crop genome evolution. The high degree of sequence identity and remarkably conserved genome structure between Arabidopsis and Brassica genomes enables comparative tiling sequencing using Arabidopsis sequences as references to select the counterpart regions in B. rapa, which is a strong challenge of structural and comparative crop genomics.

Results

We assembled 65.8 megabase-pairs of non-redundant euchromatic sequence of B. rapa and compared this sequence to the Arabidopsis genome to investigate chromosomal relationships, macrosynteny blocks, and microsynteny within blocks. The triplicated B. rapa genome contains only approximately twice the number of genes as in Arabidopsis because of genome shrinkage. Genome comparisons suggest that B. rapa has a distinct organization of ancestral genome blocks as a result of recent whole genome triplication followed by a unique diploidization process. A lack of the most recent whole genome duplication (3R) event in the B. rapa genome, atypical of other Brassica genomes, may account for the emergence of B. rapa from the Brassica progenitor around 8 million years ago.

Conclusions

This work demonstrates the potential of using comparative tiling sequencing for genome analysis of crop species. Based on a comparative analysis of the B. rapa sequences and the Arabidopsis genome, it appears that polyploidy and chromosomal diploidization are ongoing processes that collectively stabilize the B. rapa genome and facilitate its evolution.

Genome comparisons reveal a dominant mechanism of chromosome number reduction in grasses and accelerated genome evolution in Triticeae

Single-nucleotide polymorphism was used in the construction of an expressed sequence tag map of Aegilops tauschii, the diploid source of the wheat D genome. Comparisons of the map with the rice and sorghum genome sequences revealed 50 inversions and translocations; 2, 8, and 40 were assigned respectively to the rice, sorghum, and Ae. tauschii lineages, showing greatly accelerated genome evolution in the large Triticeae genomes. The reduction of the basic chromosome number from 12 to 7 in the Triticeae has taken place by a process during which an entire chromosome is inserted by its telomeres into a break in the centromeric region of another chromosome. The original centromere-telomere polarity of the chromosome arms is maintained in the new chromosome. An intrachromosomal telomere-telomere fusion resulting in a pericentric translocation of a chromosome segment or an entire arm accompanied or preceded the chromosome insertion in some instances. Insertional dysploidy has been recorded in three grass subfamilies and appears to be the dominant mechanism of basic chromosome number reduction in grasses. A total of 64% and 66% of Ae. tauschii genes were syntenic with sorghum and rice genes, respectively. Synteny was reduced in the vicinity of the termini of modern Ae. tauschii chromosomes but not in the vicinity of the ancient termini embedded in the Ae. tauschii chromosomes, suggesting that the dependence of synteny erosion on gene location along the centromere-telomere axis either evolved recently in the Triticeae phylogenetic lineage or its evolution was recently accelerated.

Development and functional assessment of EST-derived 2RL-specific markers for 2BS.2RL translocations

ESTs-derived markers are useful for comparative genomic analysis and can also serve as phenotype-linked functional markers. Here, we report the development of EST-derived 2RL-specific markers and the evaluation of the possibility of functional assessment of markers tagging 2RL, which carries Hessian fly resistance genes (loci). To identify transcripts specific to 2RL, unigene sequences in combination with wheat progenitor genomes were used. Total 275 contigs mapped to the long arms of homoeologous group 2 chromosomes were downloaded. To obtain a cluster corresponding to each of the wheat 275 contigs, unigene sequences of wheat, rice, barley, and rye were pooled for cross-species clusters. Out of 275 clusters examined, it was possible to design 112 cross-species primer pairs for genome-specific amplifications. Out of 112 cross-species primer pairs, 45 primer pairs (40%) produced amplicons from at least one species (three wheat progenitors or rye). Among the 45 contigs, 73% were associated with one of known functions and 82% of the contigs associated with known functions were also associated with one of the GO categories. On the basis of the oligonucleotide sequence alignment of each of 45 genome-specific amplifications, 21 amplifications (47%) were suitable for designing RR genome-specific primers, which are specific to translocated rye chromatin 2RL. Six primer pairs (13%) successfully produced amplicons in the 2BS.2RL translocation lines and not in the non-2RLs. Functional assessment of one of the 2RL-specific markers, NSFT03P2_Contig4445, was performed on Hessian fly infested NILs. Under Hessian fly infestation, significantly high expression of a gene tagged by a 2RL-specific marker (NSFT03P2_Contig4445) was observed 1 day after infestation. EST-derived 2RL-specific marker development from this study provides a basis for the development of ESTs-derived markers for detecting wheat-rye translocations. In addition, these markers could be employed in elucidating functional analysis of genes on 2RL.

Origin and diversification of land plant CC-type glutaredoxins

Glutaredoxins (GRXs) are ubiquitous glutathione-dependent oxidoreductase enzymes implicated in redox homeostasis, particularly oxidative stress response. Three major classes of GRX genes exist, the CPYC, CGFS classes are present in all pro- and eukaryote species, whereas the CC-type class GRXs are specific to land plants. In the basal land plant Physcomitrella patens, only two CC-type GRXs are present, compared with 21 in Arabidopsis. In contrast, sizes of the CPYC and CGFS classes remained rather similar throughout plant evolution, raising the interesting question as to when the CC-type GRXs first originated and how and why they expanded during land plant evolution. Recent evidence suggests that CC-type GRXs may have been recruited during evolution into diverse plant-specific functions of flower development (ROXY1, ROXY2) and pathogenesis response (ROXY19/GRX480). In the present study, GRX genes from the genomes of a range of green algae and evolutionarily diverse land plant species were identified; Ostreococcus, Micromonas, Volvox, Selaginella, Vitis, Sorghum, and Brachypodium. Previously identified sequences from Chlamydomonas, Physcomitrella, Oryza, Arabidopsis, and Populus were integrated to generate a more comprehensive understanding of the forces behind the evolution of various GRX classes. The analysis indicates that the CC-type GRXs probably arose by diversification from the CPYC class, at a time coinciding with colonization of land by plants. This strong differential expansion of the CC-type class occurred exclusively in the angiosperms, mainly through paleopolyploidy duplication events shortly after the monocot–eudicot split, and more recently through multiple tandem duplications that occurred independently in five investigated angiosperm lineages. The presented data suggest that following duplications, subfunctionalization, and subsequent neofunctionalization likely facilitated the sequestration of land plant-specific CC-type GRXs into novel functions such as development and pathogenesis response.

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.

Structure and expression analysis of rice paleo duplications

Having a well-known history of genome duplication, rice is a good model for studying structural and functional evolution of paleo duplications. Improved sequence alignment criteria were used to characterize 10 major chromosome-to-chromosome duplication relationships associated with 1440 paralogous pairs, covering 47.8% of the rice genome, with 12.6% of genes that are conserved within sister blocks. Using a micro-array experiment, a genome-wide expression map has been produced, in which 2382 genes show significant differences of expression in root, leaf and grain. By integrating both structural (1440 paralogous pairs) and functional information (2382 differentially expressed genes), we identified 115 paralogous gene pairs for which at least one copy is differentially expressed in one of the three tissues. A vast majority of the 115 paralogous gene pairs have been neofunctionalized or subfunctionalized as 88%, 89% and 96% of duplicates, respectively, expressed in grain, leaf and root show distinct expression patterns. On the basis of a Gene Ontology analysis, we have identified and characterized the gene families that have been structurally and functionally preferentially retained in the duplication showing that the vast majority (>85%) of duplicated have been either lost or have been subfunctionalized or neofunctionalized during 50–70 million years of evolution.

Quantitative trait loci (QTL) analysis for rice grain width and fine mapping of an identified QTL allele gw-5 in a recombination hotspot region on chromosome 5

Rice grain width and shape play a crucial role in determining grain quality and yield. The genetic basis of rice grain width was dissected into six additive quantitative trait loci (QTL) and 11 pairs of epistatic QTL using an F(7) recombinant inbred line (RIL) population derived from a single cross between Asominori (japonica) and IR24 (indica). QTL by environment interactions were evaluated in four environments. Chromosome segment substitution lines (CSSLs) harboring the six additive effect QTL were used to evaluate gene action across eight environments. A major, stable QTL, qGW-5, consistently decreased rice grain width in both the Asominori/IR24 RIL and CSSL populations with the genetic background Asominori. By investigating the distorted segregation of phenotypic values of rice grain width and genotypes of molecular markers in BC(4)F(2) and BC(4)F(3) populations, qGW-5 was dissected into a single recessive gene, gw-5, which controlled both grain width and length-width ratio. gw-5 was narrowed down to a 49.7-kb genomic region with high recombination frequencies on chromosome 5 using 6781 BC(4)F(2) individuals and 10 newly developed simple sequence repeat markers. Our results provide a basis for map-based cloning of the gw-5 gene and for marker-aided gene/QTL pyramiding in rice quality breeding.

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.

The ploidy effects in plant gene expression: progress, problems and prospects

Polyploidy and haploid are widely employed in the studies of genetics and evolution, and great progress has been made in these fields, inspiring the enthusiasm of scientists to explore the ploidy effects in gene expression. In this paper, we review the gene expression and its regulation in polyploids, especially in autopolyploids. We summarize some limitations in previous reports on polyploidy gene expression and its regulation, especially the limitations in the research materials. We propose an idea to create homologous ploidy series with twin-seedlings and to employ high-throughput techniques to investigate the polyploidy transcriptome and its regulation.

Insights into the Musa genome: syntenic relationships to rice and between Musa species

Background

Musa species (Zingiberaceae, Zingiberales) including bananas and plantains are collectively the fourth most important crop in developing countries. Knowledge concerning Musa genome structure and the origin of distinct cultivars has greatly increased over the last few years. Until now, however, no large-scale analyses of Musa genomic sequence have been conducted. This study compares genomic sequence in two Musa species with orthologous regions in the rice genome.

Results

We produced 1.4 Mb of Musa sequence from 13 BAC clones, annotated and analyzed them along with 4 previously sequenced BACs. The 443 predicted genes revealed that Zingiberales genes share GC content and distribution characteristics with eudicot and Poaceae genomes. Comparison with rice revealed microsynteny regions that have persisted since the divergence of the Commelinid orders Poales and Zingiberales at least 117 Mya. The previously hypothesized large-scale duplication event in the common ancestor of major cereal lineages within the Poaceae was verified. The divergence time distributions for Musa-Zingiber (Zingiberaceae, Zingiberales) orthologs and paralogs provide strong evidence for a large-scale duplication event in the Musa lineage after its divergence from the Zingiberaceae approximately 61 Mya. Comparisons of genomic regions from M. acuminata and M. balbisiana revealed highly conserved genome structure, and indicated that these genomes diverged circa 4.6 Mya.

Conclusion

These results point to the utility of comparative analyses between distantly-related monocot species such as rice and Musa for improving our understanding of monocot genome evolution. Sequencing the genome of M. acuminata would provide a strong foundation for comparative genomics in the monocots. In addition a genome sequence would aid genomic and genetic analyses of cultivated Musa polyploid genotypes in research aimed at localizing and cloning genes controlling important agronomic traits for breeding purposes.

Identification and characterization of shared duplications between rice and wheat provide new insight into grass genome evolution

The grass family comprises the most important cereal crops and is a good system for studying, with comparative genomics, mechanisms of evolution, speciation, and domestication. Here, we identified and characterized the evolution of shared duplications in the rice (Oryza sativa) and wheat (Triticum aestivum) genomes by comparing 42,654 rice gene sequences with 6426 mapped wheat ESTs using improved sequence alignment criteria and statistical analysis. Intraspecific comparisons identified 29 interchromosomal duplications covering 72% of the rice genome and 10 duplication blocks covering 67.5% of the wheat genome. Using the same methodology, we assessed orthologous relationships between the two genomes and detected 13 blocks of colinearity that represent 83.1 and 90.4% of the rice and wheat genomes, respectively. Integration of the intraspecific duplications data with colinearity relationships revealed seven duplicated segments conserved at orthologous positions. A detailed analysis of the length, composition, and divergence time of these duplications and comparisons with sorghum (Sorghum bicolor) and maize (Zea mays) indicated common and lineage-specific patterns of conservation between the different genomes. This allowed us to propose a model in which the grass genomes have evolved from a common ancestor with a basic number of five chromosomes through a series of whole genome and segmental duplications, chromosome fusions, and translocations.

LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons

Long terminal repeat retrotransposons (LTR elements) are ubiquitous eukaryotic transposable elements. They play important roles in the evolution of genes and genomes. Ever-growing amount of genomic sequences of many organisms present a great challenge to fast identifying them. That is the first and indispensable step to study their structure, distribution, functions and other biological impacts. However, until today, tools for efficient LTR retrotransposon discovery are very limited. Thus, we developed LTR_FINDER web server. Given DNA sequences, it predicts locations and structure of full-length LTR retrotransposons accurately by considering common structural features. LTR_FINDER is a system capable of scanning large-scale sequences rapidly and the first web server for ab initio LTR retrotransposon finding. We illustrate its usage and performance on the genome of Saccharomyces cerevisiae. The web server is freely accessible at http://tlife.fudan.edu.cn/ltr_finder/.

The FLOWERING LOCUS T-like gene family in barley (Hordeum vulgare)

The FLOWERING LOCUS T (FT) gene plays a central role in integrating flowering signals in Arabidopsis because its expression is regulated antagonistically by the photoperiod and vernalization pathways. FT belongs to a family of six genes characterized by a phosphatidylethanolamine-binding protein (PEBP) domain. In rice (Oryza sativa), 19 PEBP genes were previously described, 13 of which are FT-like genes. Five FT-like genes were found in barley (Hordeum vulgare). HvFT1, HvFT2, HvFT3, and HvFT4 were highly homologous to OsFTL2 (the Hd3a QTL), OsFTL1, OsFTL10, and OsFTL12, respectively, and this relationship was supported by comparative mapping. No rice equivalent was found for HvFT5. HvFT1 was highly expressed under long-day (inductive) conditions at the time of the morphological switch of the shoot apex from vegetative to reproductive growth. HvFT2 and HvFT4 were expressed later in development. HvFT1 was therefore identified as the main barley FT-like gene involved in the switch to flowering. Mapping of HvFT genes suggests that they provide important sources of flowering-time variation in barley. HvFTI was a candidate for VRN-H3, a dominant mutation giving precocious flowering, while HvFT3 was a candidate for Ppd-H2, a major QTL affecting flowering time in short days.

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.

Functional classification, genomic organization, putatively cis-acting regulatory elements, and relationship to quantitative trait loci, of sorghum genes with rhizome-enriched expression

Rhizomes are organs of fundamental importance to plant competitiveness and invasiveness. We have identified genes expressed at substantially higher levels in rhizomes than other plant parts, and explored their functional categorization, genomic organization, regulatory motifs, and association with quantitative trait loci (QTLs) conferring rhizomatousness. The finding that genes with rhizome-enriched expression are distributed across a wide range of functional categories suggests some degree of specialization of individual members of many gene families in rhizomatous plants. A disproportionate share of genes with rhizome-enriched expression was implicated in secondary and hormone metabolism, and abiotic stimuli and development. A high frequency of unknown-function genes reflects our still limited knowledge of this plant organ. A putative oligosaccharyl transferase showed the highest degree of rhizome-specific expression, with several transcriptional or regulatory protein complex factors also showing high (but lesser) degrees of specificity. Inferred by the upstream sequences of their putative rice (Oryza sativa) homologs, sorghum (Sorghum bicolor) genes that were relatively highly expressed in rhizome tip tissues were enriched for cis-element motifs, including the pyrimidine box, TATCCA box, and CAREs box, implicating the gibberellins in regulation of many rhizome-specific genes. From cDNA clones showing rhizome-enriched expression, expressed sequence tags forming 455 contigs were plotted on the rice genome and aligned to QTL likelihood intervals for ratooning and rhizomatous traits in rice and sorghum. Highly expressed rhizome genes were somewhat enriched in QTL likelihood intervals for rhizomatousness or ratooning, with specific candidates including some of the most rhizome-specific genes. Some rhizomatousness and ratooning QTLs were shown to be potentially related to one another as a result of ancient duplication, suggesting long-term functional conservation of the underlying genes. Insight into genes and pathways that influence rhizome growth set the stage for genetic and/or exogenous manipulation of rhizomatousness, and for further dissection of the molecular evolution of rhizomatousness.

Integration of hybridization-based markers (overgos) into physical maps for comparative and evolutionary explorations in the genus Oryza and in Sorghum

BACKGROUND

With the completion of the genome sequence for rice (Oryza sativa L.), the focus of rice genomics research has shifted to the comparison of the rice genome with genomes of other species for gene cloning, breeding, and evolutionary studies. The genus Oryza includes 23 species that shared a common ancestor 8-10 million years ago making this an ideal model for investigations into the processes underlying domestication, as many of the Oryza species are still undergoing domestication. This study integrates high-throughput, hybridization-based markers with BAC end sequence and fingerprint data to construct physical maps of rice chromosome 1 orthologues in two wild Oryza species. Similar studies were undertaken in Sorghum bicolor, a species which diverged from cultivated rice 40-50 million years ago.

RESULTS

Overgo markers, in conjunction with fingerprint and BAC end sequence data, were used to build sequence-ready BAC contigs for two wild Oryza species. The markers drove contig merges to construct physical maps syntenic to rice chromosome 1 in the wild species and provided evidence for at least one rearrangement on chromosome 1 of the O. sativa versus Oryza officinalis comparative map. When rice overgos were aligned to available S. bicolor sequence, 29% of the overgos aligned with three or fewer mismatches; of these, 41% gave positive hybridization signals. Overgo hybridization patterns supported colinearity of loci in regions of sorghum chromosome 3 and rice chromosome 1 and suggested that a possible genomic inversion occurred in this syntenic region in one of the two genomes after the divergence of S. bicolor and O. sativa.

CONCLUSION

The results of this study emphasize the importance of identifying conserved sequences in the reference sequence when designing overgo probes in order for those probes to hybridize successfully in distantly related species. As interspecific markers, overgos can be used successfully to construct physical maps in species which diverged less than 8 million years ago, and can be used in a more limited fashion to examine colinearity among species which diverged as much as 40 million years ago. Additionally, overgos are able to provide evidence of genomic rearrangements in comparative physical mapping studies.

Biochemical characterization of the major sorghum grain peroxidase

The major cationic peroxidase in sorghum grain (SPC4) , which is ubiquitously present in all sorghum varieties was purified to apparent homogeneity, and found to be a highly basic protein (pI approximately 11). MS analysis showed that SPC4 consists of two glycoforms with molecular masses of 34,227 and 35,629 Da and it contains a type-b heme. Chemical deglycosylation allowed to estimate sugar contents of 3.0% and 6.7% (w/w) in glycoform I and II, respectively, and a mass of the apoprotein of 33,246 Da. High performance anion exchange chromatography allowed to determine the carbohydrate constituents of the polysaccharide chains. The N-terminal sequence of SPC4 is not blocked by pyroglutamate. MS analysis showed that six peptides, including the N-terminal sequence of SPC4 matched with the predicted tryptic peptides of gene indice TC102191 of sorghum chromosome 1, indicating that TC102191 codes for the N-terminal part of the sequence of SPC4, including a signal peptide of 31 amino acids. The N-terminal fragment of SPC4 (213 amino acids) has a high sequence identity with barley BP1 (85%), rice Prx23 (90%), wheat WSP1 (82%) and maize peroxidase (58%), indicative for a common ancestor. SPC4 is activated by calcium ions. Ca2+ binding increased the protein conformational stability by raising the melting temperature (Tm) from 67 to 82 degrees C. SPC4 catalyzed the oxidation of a wide range of aromatic substrates, being catalytically more efficient with hydroxycinnamates than with tyrosine derivatives. In spite of the conserved active sites, SPC4 differs from BP1 in being active with aromatic compounds above pH 5.

A comparative genomics strategy for targeted discovery of single-nucleotide polymorphisms and conserved-noncoding sequences in orphan crops

Completed genome sequences provide templates for the design of genome analysis tools in orphan species lacking sequence information. To demonstrate this principle, we designed 384 PCR primer pairs to conserved exonic regions flanking introns, using Sorghum/Pennisetum expressed sequence tag alignments to the Oryza genome. Conserved-intron scanning primers (CISPs) amplified single-copy loci at 37% to 80% success rates in taxa that sample much of the approximately 50-million years of Poaceae divergence. While the conserved nature of exons fostered cross-taxon amplification, the lesser evolutionary constraints on introns enhanced single-nucleotide polymorphism detection. For example, in eight rice (Oryza sativa) genotypes, polymorphism averaged 12.1 per kb in introns but only 3.6 per kb in exons. Curiously, among 124 CISPs evaluated across Oryza, Sorghum, Pennisetum, Cynodon, Eragrostis, Zea, Triticum, and Hordeum, 23 (18.5%) seemed to be subject to rigid intron size constraints that were independent of per-nucleotide DNA sequence variation. Furthermore, we identified 487 conserved-noncoding sequence motifs in 129 CISP loci. A large CISP set (6,062 primer pairs, amplifying introns from 1,676 genes) designed using an automated pipeline showed generally higher abundance in recombinogenic than in nonrecombinogenic regions of the rice genome, thus providing relatively even distribution along genetic maps. CISPs are an effective means to explore poorly characterized genomes for both DNA polymorphism and noncoding sequence conservation on a genome-wide or candidate gene basis, and also provide anchor points for comparative genomics across a diverse range of species.

Leafing through the genomes of our major crop plants: strategies for capturing unique information

Crop plants not only have economic significance, but also comprise important botanical models for evolution and development. This is reflected by the recent increase in the percentage of publicly available sequence data that are derived from angiosperms. Further genome sequencing of the major crop plants will offer new learning opportunities, but their large, repetitive, and often polyploid genomes present challenges. Reduced-representation approaches - such as EST sequencing, methyl filtration and Cot-based cloning and sequencing - provide increased efficiency in extracting key information from crop genomes without full-genome sequencing. Combining these methods with phylogenetically stratified sampling to allow comparative genomic approaches has the potential to further accelerate progress in angiosperm genomics.

Computational identification and systematic analysis of the ACR gene family in Oryza sativa

Based on sequence similarity search and domain detection, nine ACT domain repeat protein-coding genes (the "ACR" genes) in rice were identified, which were mainly distributed on the chromosomes 2, 3, 4, and 8. An InterPro database search indicated that four copies of the ACT domain linearly occupied the entire polypeptide. The first three ACT domains were linked by two different sequences. However, the fourth ACT domain was extremely close to ACT3. Gene structure comparisons showed large differences in exon numbers, from three to eight, among members of the rice ACR gene family. In addition, it appeared that gene duplication might be operative when the compositions of exons and introns were analyzed. Phylogenetic analysis divided the ACR gene family into five distinct groups, and this division was generally according to the expression patterns of the ACR genes. The Arabidopsis and rice ACR proteins were clustered across together, suggesting that these ACR genes might originate from an ancient common ancestor. Notably, the identification of orthologues and paralogues would be useful for rice gene functional annotation.

Grains of knowledge: genomics of model cereals

The economic and scientific importance of the cereals has motivated a rich history of research into their genetics, development, and evolution. The nearly completed sequence of the rice genome is emblematic of a transition to high-throughput genomics and computational biology that has also pervaded study of many other cereals. The relatively close (ca. <50 million years old) relationships among morphologically diverse cereals native to environments that sample much of global geographic diversity make the cereals particularly attractive for comparative studies of plant genome evolution. Extensive germplasm resources, largely a byproduct of their economic importance, together with growing collections of defined mutants, provide foundations for a host of post-genomic studies to shed more light on the relationship between sequence and function in this important group. Using the rapidly growing capabilities of several informatics resources, genomic data from model cereals are likely to be leveraged tremendously in the study and improvement of a wide range of crop plants that sustain much of the world's population, including many which still lack primary genomic resources.

Buffering of crucial functions by paleologous duplicated genes may contribute cyclicality to angiosperm genome duplication

Genome duplication followed by massive gene loss has permanently shaped the genomes of many higher eukaryotes, particularly angiosperms. It has long been believed that a primary advantage of genome duplication is the opportunity for the evolution of genes with new functions by modification of duplicated genes. If so, then patterns of genetic diversity among strains within taxa might reveal footprints of selection that are consistent with this advantage. Contrary to classical predictions that duplicated genes may be relatively free to acquire unique functionality, we find among both Arabidopsis ecotypes and Oryza subspecies that SNPs encode less radical amino acid changes in genes for which there exists a duplicated copy at a "paleologous" locus than in "singleton" genes. Preferential retention of duplicated genes encoding long complex proteins and their unexpectedly slow divergence (perhaps because of homogenization) suggest that a primary advantage of retaining duplicated paleologs may be the buffering of crucial functions. Functional buffering and functional divergence may represent extremes in the spectrum of duplicated gene fates. Functional buffering may be especially important during "genomic turmoil" immediately after genome duplication but continues to act approximately 60 million years later, and its gradual deterioration may contribute cyclicality to genome duplication in some lineages.

Genome-wide comparative analysis of the IQD gene families in Arabidopsis thaliana and Oryza sativa

BACKGROUND

Calcium signaling plays a prominent role in plants for coordinating a wide range of developmental processes and responses to environmental cues. Stimulus-specific generation of intracellular calcium transients, decoding of calcium signatures, and transformation of the signal into cellular responses are integral modules of the transduction process. Several hundred proteins with functions in calcium signaling circuits have been identified, and the number of downstream targets of calcium sensors is expected to increase. We previously identified a novel, calmodulin-binding nuclear protein, IQD1, which stimulates glucosinolate accumulation and plant defense in Arabidopsis thaliana. Here, we present a comparative genome-wide analysis of a new class of putative calmodulin target proteins in Arabidopsis and rice.

RESULTS

We identified and analyzed 33 and 29 IQD1-like genes in Arabidopsis thaliana and Oryza sativa, respectively. The encoded IQD proteins contain a plant-specific domain of 67 conserved amino acid residues, referred to as the IQ67 domain, which is characterized by a unique and repetitive arrangement of three different calmodulin recruitment motifs, known as the IQ, 1-5-10, and 1-8-14 motifs. We demonstrated calmodulin binding for IQD20, the smallest IQD protein in Arabidopsis, which consists of a C-terminal IQ67 domain and a short N-terminal extension. A striking feature of IQD proteins is the high isoelectric point (approximately 10.3) and frequency of serine residues (approximately 11%). We compared the Arabidopsis and rice IQD gene families in terms of gene structure, chromosome location, predicted protein properties and motifs, phylogenetic relationships, and evolutionary history. The existence of an IQD-like gene in bryophytes suggests that IQD proteins are an ancient family of calmodulin-binding proteins and arose during the early evolution of land plants.

CONCLUSION

Comparative phylogenetic analyses indicate that the major IQD gene lineages originated before the monocot-eudicot divergence. The extant IQD loci in Arabidopsis primarily resulted from segmental duplication and reflect preferential retention of paralogous genes, which is characteristic for proteins with regulatory functions. Interaction of IQD1 and IQD20 with calmodulin and the presence of predicted calmodulin binding sites in all IQD family members suggest that IQD proteins are a new class of calmodulin targets. The basic isoelectric point of IQD proteins and their frequently predicted nuclear localization suggest that IQD proteins link calcium signaling pathways to the regulation of gene expression. Our comparative genomics analysis of IQD genes and encoded proteins in two model plant species provides the first step towards the functional dissection of this emerging family of putative calmodulin targets.

Comprehensive molecular cytogenetic analysis of sorghum genome architecture: distribution of euchromatin, heterochromatin, genes and recombination in comparison to rice

Cytogenetic maps of sorghum chromosomes 3-7, 9, and 10 were constructed on the basis of the fluorescence in situ hybridization (FISH) of approximately 18-30 BAC probes mapped across each of these chromosomes. Distal regions of euchromatin and pericentromeric regions of heterochromatin were delimited for all 10 sorghum chromosomes and their DNA content quantified. Euchromatic DNA spans approximately 50% of the sorghum genome, ranging from approximately 60% of chromosome 1 (SBI-01) to approximately 33% of chromosome 7 (SBI-07). This portion of the sorghum genome is predicted to encode approximately 70% of the sorghum genes ( approximately 1 gene model/12.3 kbp), assuming that rice and sorghum encode a similar number of genes. Heterochromatin spans approximately 411 Mbp of the sorghum genome, a region characterized by a approximately 34-fold lower rate of recombination and approximately 3-fold lower gene density compared to euchromatic DNA. The sorghum and rice genomes exhibit a high degree of macrocolinearity; however, the sorghum genome is approximately 2-fold larger than the rice genome. The distal euchromatic regions of sorghum chromosomes 3-7 and 10 are approximately 1.8-fold larger overall and exhibit an approximately 1.5-fold lower average rate of recombination than the colinear regions of the homeologous rice chromosomes. By contrast, the pericentromeric heterochromatic regions of these chromosomes are on average approximately 3.6-fold larger in sorghum and recombination is suppressed approximately 15-fold compared to the colinear regions of rice chromosomes.

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.

Phylogenomic analysis of the PEBP gene family in cereals

The TFL1 and FT genes, which are key genes in the control of flowering time in Arabidopsis thaliana, belong to a small multigene family characterized by a specific phosphatidylethanolamine-binding protein domain, termed the PEBP gene family. Several PEBP genes are found in dicots and monocots, and act on the control of flowering time. We investigated the evolution of the PEBP gene family in cereals. First, taking advantage of the complete rice genome sequence and EST databases, we found 19 PEBP genes in this species, 6 of which were not previously described. Ten genes correspond to five pairs of paralogs mapped on known duplicated regions of the rice genome. Phylogenetic analysis of Arabidopsis and rice genes indicates that the PEBP gene family consists of three main homology classes (the so-called TFL1-LIKE, MFT-LIKE, and FT-LIKE subfamilies), in which gene duplication and/or loss occurred independently in Arabidopsis and rice. Second, phylogenetic analyses of genomic and EST sequences from five cereal species indicate that the three subfamilies of PEBP genes have been conserved in cereals. The tree structure suggests that the ancestral grass genome had at least two MFT-like genes, two TFL1-like genes, and eight FT-like genes. A phylogenomic approach leads to some hypotheses about conservation of gene function within the subfamilies.

How can we use genomics to improve cereals with rice as a reference genome?

Rice serves as a model crop for cereal genomics. The availability of complete genome sequences, together with various genomic resources available for both rice and Arabidopsis, have revolutionized our understanding of the genetic make-up of crop plants. Both macrocolinearity revealed by comparative mapping and microcolinearity revealed by sequence comparisons among the grasses indicate that sequencing and functional analysis of the rice genome will have a significant impact on other cereals in terms of both genomic studies and crop improvement. The availability of mutants, introgression libraries, and advanced transformation techniques make functional genomics in rice and other cereals more manageable than ever before. A wide array of genetic markers, including anchor markers for comparative mapping, SSRs and SNPs are widely used in genetic mapping, germplasm evaluation and marker assisted selection. An integrated database that combines genome information for rice and other cereals is key to the effective utilization of all genomics resources for cereal improvement. To maximize the potential of genomics for plant breeding, experiments must be further miniaturized and costs must be reduced. Many techniques, including targeted gene disruption or allele substitution, insertional mutagenesis, RNA interference and homologous recombination, need to be refined before they can be widely used in functional genomic analysis and plant breeding.

A pilot study of transcription unit analysis in rice using oligonucleotide tiling-path microarray

As the international efforts to sequence the rice genome are completed, an immediate challenge and opportunity is to comprehensively and accurately define all transcription units in the rice genome. Here we describe a strategy of using high-density oligonucleotide tiling-path microarrays to map transcription of the japonica rice genome. In a pilot experiment to test this approach, one array representing the reverse strand of the last 11.2 Mb sequence of chromosome 10 was analyzed in detail based on a mathematical model developed in this study. Analysis of the array data detected 77% of the reference gene models in a mixture of four RNA populations. Moreover, significant transcriptional activities were found in many of the previously annotated intergenic regions. These preliminary results demonstrate the utility of genome tiling microarrays in evaluating annotated rice gene models and in identifying novel transcription units that will facilitate rice genome annotation.

Comparative genomics of Gossypium and Arabidopsis: unraveling the consequences of both ancient and recent polyploidy

Both ancient and recent polyploidy, together with post-polyploidization loss of many duplicated gene copies, complicates angiosperm comparative genomics. To explore an approach by which these challenges might be mitigated, genetic maps of extant diploid and tetraploid cottons (Gossypium spp.) were used to infer the approximate order of 3016 loci along the chromosomes of their hypothetical common ancestor. The inferred Gossypium gene order corresponded more closely than the original maps did to a similarly inferred ancestral gene order predating an independent paleopolyploidization (alpha) in Arabidopsis. At least 59% of the cotton map and 53% of the Arabidopsis transcriptome showed correspondence in multilocus gene arrangements based on one or both of two software packages (CrimeStatII, FISH). Genomic regions in which chromosome structural rearrangement has been rapid (obscuring gene order correspondence) have also been subject to greater divergence of individual gene sequences. About 26%-44% of corresponding regions involved multiple Arabidopsis or cotton chromosomes, in some cases consistent with known, more ancient, duplications. The genomic distributions of multiple-locus probes provided early insight into the consequences for chromosome structure of an ancient large-scale duplication in cotton. Inferences that mitigate the consequences of ancient duplications improve leveraging of genomic information for model organisms in the study of more complex genomes.

Extensive de Novo genomic variation in rice induced by introgression from wild rice (Zizania latifolia Griseb.)

To study the possible impact of alien introgression on a recipient plant genome, we examined >6000 unbiased genomic loci of three stable rice recombinant inbred lines (RILs) derived from intergeneric hybridization between rice (cv. Matsumae) and a wild relative (Zizania latifolia Griseb.) followed by successive selfing. Results from amplified fragment length polymorphism (AFLP) analysis showed that, whereas the introgressed Zizania DNA comprised <0.1% of the genome content in the RILs, extensive and genome-wide de novo variations occurred in up to 30% of the analyzed loci for all three lines studied. The AFLP-detected changes were validated by DNA gel-blot hybridization and/or sequence analysis of genomic loci corresponding to a subset of the differentiating AFLP fragments. A BLAST analysis revealed that the genomic variations occurred in diverse sequences, including protein-coding genes, transposable elements, and sequences of unknown functions. Pairwise sequence comparison of selected loci between a RIL and its rice parent showed that the variations represented either base substitutions or small insertion/deletions. Genome variations were detected in all 12 rice chromosomes, although their distribution was uneven both among and within chromosomes. Taken together, our results imply that even cryptic alien introgression can be highly mutagenic to a recipient plant genome.

Phylogenomic analysis of the receptor-like proteins of rice and Arabidopsis

The tomato (Lycopersicon esculentum) Cf-9 resistance gene encodes the first characterized member of the plant receptor-like protein (RLP) family. Other RLPs such as CLAVATA2 and TOO MANY MOUTHS are known to regulate development. The domain structure of RLPs consists of extracellular leucine-rich repeats, a transmembrane helix, and a short cytoplasmic region. Here, we identify 90 RLPs in rice (Oryza sativa) and compare them with functionally characterized RLPs from different plant species and with 56 Arabidopsis (Arabidopsis thaliana) RLPs, including the downy mildew resistance protein RPP27. Many RLPs cluster into four distinct superclades, three of which include RLPs known to be involved in plant defense. Sequence comparisons reveal diagnostic amino acid residues that may specify different molecular functions in different RLP subtypes. This analysis of rice RLPs thus identified at least 73 candidate resistance genes and four genes potentially involved in development. Due to the synteny between rice and other Gramineae, this analysis should provide valuable tools for experimental studies in rice and other cereals.

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

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

Large intraspecific haplotype variability at the Rph7 locus results from rapid and recent divergence in the barley genome

To study genome evolution and diversity in barley (Hordeum vulgare), we have sequenced and compared more than 300 kb of sequence spanning the Rph7 leaf rust disease resistance gene in two barley cultivars. Colinearity was restricted to five genic and two intergenic regions representing <35% of the two sequences. In each interval separating the seven conserved regions, the number and type of repetitive elements were completely different between the two homologous sequences, and a single gene was absent in one cultivar. In both cultivars, the nonconserved regions consisted of approximately 53% repetitive sequences mainly represented by long-terminal repeat retrotransposons that have inserted <1 million years ago. PCR-based analysis of intergenic regions at the Rph7 locus and at three other independent loci in 41 H. vulgare lines indicated large haplotype variability in the cultivated barley gene pool. Together, our data indicate rapid and recent divergence at homologous loci in the genome of H. vulgare, possibly providing the molecular mechanism for the generation of high diversity in the barley gene pool. Finally, comparative analysis of the gene composition in barley, wheat (Triticum aestivum), rice (Oryza sativa), and sorghum (Sorghum bicolor) suggested massive gene movements at the Rph7 locus in the Triticeae lineage.

Two ancient rounds of polyploidy in rice genome

An ancient genome duplication (PPP1) that predates divergence of the cereals has recently been recognized. We report here another potentially older large-scale duplication (PPP2) event that predates monocot-dicot divergence in the genome of rice (Oryza sativa L.), as inferred from the age distribution of pairs of duplicate genes based on recent genome data for rice. Our results suggest that paleopolyploidy was widespread and played an important role in the evolution of rice.