Genome-wide intraspecific DNA-sequence variations in rice

Maroon Rice Genomic Diversity Reflects 350 Years of Colonial History

Maroons in Suriname and French Guiana descend from enslaved Africans who escaped the plantations during colonial times. Maroon farmers still cultivate a large diversity of rice, their oldest staple crop. The oral history and written records of Maroons by colonial authorities provide contrasting perspectives on the origins of Maroon rice. Here, we analyzed the genomic ancestry of 136 newly sequenced Maroon rice varieties and found seven genomic groups that differ in their geographical associations. We interpreted these findings in light of ethnobotanical and archival investigations to reconstruct the historical contexts associated with the introduction of rice varieties to the Guianas. We found that two rice groups trace to West Africa, which we propose are linked to the transatlantic slave trade (c. 1526 to 1825). We posit that the Maroon rice stock additionally contains varieties that derive from rice introduced by indentured laborers from Java (1890 onwards), USA rice breeders (1932 onwards), and Hmong refugees who fled the Vietnam War (1991). Furthermore, on the Maroon fields, we found rice types never documented before that were derived from crosses. Overall, our results demonstrate that the Maroon farmers prioritize maintenance of a high stock diversity, which we posit reflects the expertise they inherited from their (African) ancestors. Ignored by agricultural modernization initiatives, Maroon farmers today are custodians of a unique cultural heritage. Notably, the genomic findings underline many Maroon stories about their past. We anticipate that a similar study approach can be applied to other heirloom crops of (Indigenous) communities that may have preserved their history on their farms to reconstruct, acknowledge, and honor the past.

DNAMarkMaker: streamlining ARMS and CAPS marker development from resequencing data with NGS short reads

DNA markers serve as essential tools in breeding selection and genetic analysis. However, developing DNA markers can be time-consuming and labor-intensive due to the need to identify polymorphisms between cultivars/lines and to design suitable primers. To address these challenges, we have developed DNAMarkMaker, a tool designed to automate the process of primer design for Amplification Refractory Mutation System (ARMS) and Cleaved Amplified Polymorphic Sequences (CAPS) markers, utilizing resequencing data. One key feature of DNAMarkMaker is its user-friendly graphical user interface (GUI), ensuring its accessibility and ease of use, even for researchers not well-versed in bioinformatics. We confirmed DNAMarkMaker's applicability by developing DNA markers for rice, potato, and turnip-each representing distinct genome structures: homozygous diploid, heterozygous autotetraploid, and heterozygous diploid, respectively. DNAMarkMaker will contribute to the rapid and efficient development of DNA markers, accelerating breeding and genetic analysis in various crops.

A Comprehensive Genome-Wide Investigation of the Cytochrome 71 (OsCYP71) Gene Family: Revealing the Impact of Promoter and Gene Variants (Ser33Leu) of OsCYP71P6 on Yield-Related Traits in Indica Rice (Oryza sativa L.)

The cytochrome P450 (CYP450) gene family plays a critical role in plant growth and developmental processes, nutrition, and detoxification of xenobiotics in plants. In the present research, a comprehensive set of 105 OsCYP71 family genes was pinpointed within the genome of indica rice. These genes were categorized into twelve distinct subfamilies, where members within the same subgroup exhibited comparable gene structures and conserved motifs. In addition, 105 OsCYP71 genes were distributed across 11 chromosomes, and 36 pairs of OsCYP71 involved in gene duplication events. Within the promoter region of OsCYP71, there exists an extensive array of cis-elements that are associated with light responsiveness, hormonal regulation, and stress-related signaling. Further, transcriptome profiling revealed that a majority of the genes exhibited responsiveness to hormones and were activated across diverse tissues and developmental stages in rice. The OsCYP71P6 gene is involved in insect resistance, senescence, and yield-related traits in rice. Hence, understanding the association between OsCYP71P6 genetic variants and yield-related traits in rice varieties could provide novel insights for rice improvement. Through the utilization of linear regression models, a total of eight promoters were identified, and a specific gene variant (Ser33Leu) within OsCYP71P6 was found to be linked to spikelet fertility. Additionally, different alleles of the OsCYP71P6 gene identified through in/dels polymorphism in 131 rice varieties were validated for their allelic effects on yield-related traits. Furthermore, the single-plant yield, spikelet number, panicle length, panicle weight, and unfilled grain per panicle for the OsCYP71P6-1 promoter insertion variant were found to contribute 20.19%, 13.65%, 5.637%, 8.79%, and 36.86% more than the deletion variant, respectively. These findings establish a robust groundwork for delving deeper into the functions of OsCYP71-family genes across a range of biological processes. Moreover, these findings provide evidence that allelic variation in the promoter and amino acid substitution of Ser33Leu in the OsCYP71P6 gene could potentially impact traits related to rice yield. Therefore, the identified promoter variants in the OsCYP71P6 gene could be harnessed to amplify rice yields.

A New Way of Rice Breeding: Polyploid Rice Breeding

Polyploid rice, first discovered by Japanese scientist Eiiti Nakamori in 1933, has a history of nearly 90 years. In the following years, polyploid rice studies have mainly focused on innovations in breeding theory, induction technology and the creation of new germplasm, the analysis of agronomic traits and nutritional components, the study of gametophyte development and reproduction characteristics, DNA methylation modification and gene expression regulation, distant hybridization and utilization among subspecies, species and genomes. In recent years, PMeS lines and neo-tetraploid rice lines with stable high seed setting rate characteristics have been successively selected, breaking through the bottleneck of low seed setting rate of polyploid rice. Following, a series of theoretical and applied studies on high seed setting rate tetraploid rice were carried out. This has pushed research on polyploid rice to a new stage, opening new prospects for polyploid rice breeding.

A Two-in-One Strategy: Target and Nontarget Site Mechanisms Both Play Important Role in IMI-Resistant Weedy Rice

The introduction of Clearfield technology allows the use of imidazolinone (IMI) herbicides to control weedy rice. Imidazolinone herbicides stop the acetolactate synthase (ALS) enzyme from synthesizing branched-chain amino acids, resulting in the death of the plant. Since the launch of Clearfield technology in Malaysia in 2010, many farmers have replaced traditional cultivars with Clearfield (CL) rice lines (MR220-CL1 and MR220-CL2). This technology was initially effective; however, in recent years, local farmers have reported the reduced efficacy of IMI herbicides in controlling the spread of weedy rice. Under IMI herbicide treatment, in previous weedy rice studies, the target-site resistance (TSR) mechanism of the ALS gene has been suggested as a key factor conferring herbicide resistance. In our study, a combination of ALS gene sequencing, enzyme colorimetric assay, and a genome-wide association study (GWAS) highlighted that a non-target-site resistance (NTSR) can be an alternative molecular mechanism in IMI-resistant weedy rice. This is supported by a series of evidence, including a weak correlation between single nucleotide polymorphisms (SNPs) within the ALS exonic region and ALS enzyme activity. Our findings suggest that the adaptability of weedy rice in Clearfield rice fields can be more complicated than previously found in other rice strains.

Variations in the Structure and Evolution of Rice WRKY Genes in Indica and Japonica Genotypes and their Co-expression Network in Mediating Disease Resistance

WRKY transcription factor (TF) family regulates many functions in plant growth and development and also during biotic and abiotic stress. In this study, 101 WRKY TF gene models in indica and japonica rice were used to conduct evolutionary analysis, gene structure analysis, and motif composition. Co-expression analysis was carried out first by selecting the differentially expressing genes that showed a significant change in response to the pathogens from Rice Oligonucleotide Array Database (ROAD). About 82 genes showed responses to infection by Magnaporthe oryzae or Xanthomonas oryzae pv. oryzae. Co-expression gene network was constructed using direct neighborhood and context associated inbuilt mode in RiceNetv2 tool. Only 41 genes showed interaction with 2299 non-WRKY genes. Variations exist in the structure and evolution of WRKY genes among indica and japonica genotypes which have important implications in their differential roles including disease resistance. WRKY genes mediate a complex networking and co-express along with other WRKY and non-WRKY genes to mediate resistance against fungal and bacterial pathogens in rice.

Genomic variation in 3,010 diverse accessions of Asian cultivated rice

Here we analyse genetic variation, population structure and diversity among 3,010 diverse Asian cultivated rice (Oryza sativa L.) genomes from the 3,000 Rice Genomes Project. Our results are consistent with the five major groups previously recognized, but also suggest several unreported subpopulations that correlate with geographic location. We identified 29 million single nucleotide polymorphisms, 2.4 million small indels and over 90,000 structural variations that contribute to within- and between-population variation. Using pan-genome analyses, we identified more than 10,000 novel full-length protein-coding genes and a high number of presence-absence variations. The complex patterns of introgression observed in domestication genes are consistent with multiple independent rice domestication events. The public availability of data from the 3,000 Rice Genomes Project provides a resource for rice genomics research and breeding.

Genomic variation in 3,010 diverse accessions of Asian cultivated rice

Here we analyse genetic variation, population structure and diversity among 3,010 diverse Asian cultivated rice (Oryza sativa L.) genomes from the 3,000 Rice Genomes Project. Our results are consistent with the five major groups previously recognized, but also suggest several unreported subpopulations that correlate with geographic location. We identified 29 million single nucleotide polymorphisms, 2.4 million small indels and over 90,000 structural variations that contribute to within- and between-population variation. Using pan-genome analyses, we identified more than 10,000 novel full-length protein-coding genes and a high number of presence–absence variations. The complex patterns of introgression observed in domestication genes are consistent with multiple independent rice domestication events. The public availability of data from the 3,000 Rice Genomes Project provides a resource for rice genomics research and breeding.

Analyses of genetic variation and population structure based on over 3,000 cultivated rice (Oryza sativa) genomes reveal subpopulations that correlate with geographic location and patterns of introgression consistent with multiple rice domestication events.

Towards plant pangenomics

As an increasing number of genome sequences become available for a wide range of species, there is a growing understanding that the genome of a single individual is insufficient to represent the gene diversity within a whole species. Many studies examine the sequence diversity within genes, and this allelic variation is an important source of phenotypic variation which can be selected for by man or nature. However, the significant gene presence/absence variation that has been observed within species and the impact of this variation on traits is only now being studied in detail. The sum of the genes for a species is termed the pangenome, and the determination and characterization of the pangenome is a requirement to understand variation within a species. In this review, we explore the current progress in pangenomics as well as methods and approaches for the characterization of pangenomes for a wide range of plant species.

Whole genome de novo assemblies of three divergent strains of rice, Oryza sativa, document novel gene space of aus and indica

BACKGROUND

The use of high throughput genome-sequencing technologies has uncovered a large extent of structural variation in eukaryotic genomes that makes important contributions to genomic diversity and phenotypic variation. When the genomes of different strains of a given organism are compared, whole genome resequencing data are typically aligned to an established reference sequence. However, when the reference differs in significant structural ways from the individuals under study, the analysis is often incomplete or inaccurate.

RESULTS

Here, we use rice as a model to demonstrate how improvements in sequencing and assembly technology allow rapid and inexpensive de novo assembly of next generation sequence data into high-quality assemblies that can be directly compared using whole genome alignment to provide an unbiased assessment. Using this approach, we are able to accurately assess the "pan-genome" of three divergent rice varieties and document several megabases of each genome absent in the other two.

CONCLUSIONS

Many of the genome-specific loci are annotated to contain genes, reflecting the potential for new biological properties that would be missed by standard reference-mapping approaches. We further provide a detailed analysis of several loci associated with agriculturally important traits, including the S5 hybrid sterility locus, the Sub1 submergence tolerance locus, the LRK gene cluster associated with improved yield, and the Pup1 cluster associated with phosphorus deficiency, illustrating the utility of our approach for biological discovery. All of the data and software are openly available to support further breeding and functional studies of rice and other species.

Whole genome de novo assemblies of three divergent strains of rice, Oryza sativa, document novel gene space of aus and indica

Background

The use of high throughput genome-sequencing technologies has uncovered a large extent of structural variation in eukaryotic genomes that makes important contributions to genomic diversity and phenotypic variation. When the genomes of different strains of a given organism are compared, whole genome resequencing data are typically aligned to an established reference sequence. However, when the reference differs in significant structural ways from the individuals under study, the analysis is often incomplete or inaccurate.

Results

Here, we use rice as a model to demonstrate how improvements in sequencing and assembly technology allow rapid and inexpensive de novo assembly of next generation sequence data into high-quality assemblies that can be directly compared using whole genome alignment to provide an unbiased assessment. Using this approach, we are able to accurately assess the ‘pan-genome’ of three divergent rice varieties and document several megabases of each genome absent in the other two.

Conclusions

Many of the genome-specific loci are annotated to contain genes, reflecting the potential for new biological properties that would be missed by standard reference-mapping approaches. We further provide a detailed analysis of several loci associated with agriculturally important traits, including the S5 hybrid sterility locus, the Sub1 submergence tolerance locus, the LRK gene cluster associated with improved yield, and the Pup1 cluster associated with phosphorus deficiency, illustrating the utility of our approach for biological discovery. All of the data and software are openly available to support further breeding and functional studies of rice and other species.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0506-z) contains supplementary material, which is available to authorized users.

A strategy for genome-wide identification of gene based polymorphisms in rice reveals non-synonymous variation and functional genotypic markers

The genetic diversity of plants has traditionally been employed to improve crop plants to suit human needs, and in the future feed the increasing population and protect crops from environmental stresses and climate change. Genome-wide sequencing is a reality and can be used to make association to crop traits to be utilized by high-throughput marker based selection methods. This study describes a strategy of using next generation sequencing (NGS) data from the rice genome to make comparisons to the high-quality reference genome, identify functional polymorphisms within genes that might result in function changes and be used to study correlations to traits and employed in genetic mapping. We analyzed the NGS data of Oryza sativa ssp indica cv. G4 covering 241 Mb with ∼20X coverage and compared to the reference genome of Oryza sativa ssp. japonica to describe the genome-wide distribution of gene-based single nucleotide polymorphisms (SNPs). The analysis shows that the 63% covered genome consists of 1.6 million SNPs with 6.9 SNPs/Kb, and including 80,146 insertions and 92,655 deletions (INDELs) genome-wide. There are a total of 1,139,801 intergenic SNPs, 295,136 SNPs in intronic/non-coding regions, 195,098 in coding regions, 23,242 SNPs at the five-prime (5′) UTR regions and 22,686 SNPs at the three-prime (3′) UTR region. SNP variation was found in 40,761 gene loci, which include 75,262 synonymous and 119,836 non-synonymous changes, and functional reading frame changes through 3,886 inducing STOP-codon (isSNP) and 729 preventing STOP-codon (psSNP) variation. There are quickly evolving 194 high SNP hotspot genes (>100 SNPs/gene), and 1,513 out of 2,458 transcription factors displaying 2,294 non-synonymous SNPs that can be a major source of phenotypic diversity within the species. All data is searchable at https://plantstress-pereira.uark.edu/oryza2. We envision that this strategy will be useful for the identification of genes for crop traits and molecular breeding of rice cultivars.

Genome-wide association analyses provide genetic and biochemical insights into natural variation in rice metabolism

Plant metabolites are important to world food security in terms of maintaining sustainable yield and providing food with enriched phytonutrients. Here we report comprehensive profiling of 840 metabolites and a further metabolic genome-wide association study based on ∼6.4 million SNPs obtained from 529 diverse accessions of Oryza sativa. We identified hundreds of common variants influencing numerous secondary metabolites with large effects at high resolution. We observed substantial heterogeneity in the natural variation of metabolites and their underlying genetic architectures among different subspecies of rice. Data mining identified 36 candidate genes modulating levels of metabolites that are of potential physiological and nutritional importance. As a proof of concept, we functionally identified or annotated five candidate genes influencing metabolic traits. Our study provides insights into the genetic and biochemical bases of rice metabolome variation and can be used as a powerful complementary tool to classical phenotypic trait mapping for rice improvement.

Global genomic diversity of Oryza sativa varieties revealed by comparative physical mapping

Bacterial artificial chromosome (BAC) physical maps embedding a large number of BAC end sequences (BESs) were generated for Oryza sativa ssp. indica varieties Minghui 63 (MH63) and Zhenshan 97 (ZS97) and were compared with the genome sequences of O. sativa spp. japonica cv. Nipponbare and O. sativa ssp. indica cv. 93-11. The comparisons exhibited substantial diversities in terms of large structural variations and small substitutions and indels. Genome-wide BAC-sized and contig-sized structural variations were detected, and the shared variations were analyzed. In the expansion regions of the Nipponbare reference sequence, in comparison to the MH63 and ZS97 physical maps, as well as to the previously constructed 93-11 physical map, the amounts and types of the repeat contents, and the outputs of gene ontology analysis, were significantly different from those of the whole genome. Using the physical maps of four wild Oryza species from OMAP (http://www.omap.org) as a control, we detected many conserved and divergent regions related to the evolution process of O. sativa. Between the BESs of MH63 and ZS97 and the two reference sequences, a total of 1532 polymorphic simple sequence repeats (SSRs), 71,383 SNPs, 1767 multiple nucleotide polymorphisms, 6340 insertions, and 9137 deletions were identified. This study provides independent whole-genome resources for intra- and intersubspecies comparisons and functional genomics studies in O. sativa. Both the comparative physical maps and the GBrowse, which integrated the QTL and molecular markers from GRAMENE (http://www.gramene.org) with our physical maps and analysis results, are open to the public through our Web site (http://gresource.hzau.edu.cn/resource/resource.html).

Genomic resources for functional analyses of the rice genome

With the availability of the rice genome sequence, rice research communities are entering a new era of plant functional genomics. The last decade has seen rapid worldwide progress on establishing platforms for rice functional genomic research. These platforms offer practical toolkits and genomic resources for high-throughput identification of genes and pathways. In this review, we summarize available genomic resources for functional analyses of the rice genome. These genomic resources include high-quality bacterial artificial chromosome libraries, large-scale expression sequence tags, full-length cDNA collections, large amounts of data on global expression profiles, various mutant libraries and integrated bioinformatics databases. We not only present the current status of genomic resources but also discuss their usage in elucidating gene functions of the rice genome.

Understanding reproductive isolation based on the rice model

Reproductive isolation is both an indicator of speciation and a mechanism for maintaining species identity. Here we review the progress in studies of hybrid sterility in rice to illustrate the present understanding of the molecular and evolutionary mechanisms underlying reproductive isolation. Findings from molecular characterization of genes controlling hybrid sterility can be summarized with three evolutionary genetic models. The parallel divergence model features duplicated loci generated by genome evolution; in this model, the gametes abort when the two copies of loss-of-function mutants meet in hybrids. In the sequential divergence model, mutations of two linked loci occur sequentially in one lineage, and negative interaction between the ancestral and nascent alleles of different genes causes incompatibility. The parallel-sequential divergence model involves three tightly linked loci, exemplified by a killer-protector system formed of mutations in two steps. We discuss the significance of such findings and their implications for crop improvement.

SNPs in stress-responsive rice genes: validation, genotyping, functional relevance and population structure

Background

Single nucleotide polymorphism (SNP) validation and large-scale genotyping are required to maximize the use of DNA sequence variation and determine the functional relevance of candidate genes for complex stress tolerance traits through genetic association in rice. We used the bead array platform-based Illumina GoldenGate assay to validate and genotype SNPs in a select set of stress-responsive genes to understand their functional relevance and study the population structure in rice.

Results

Of the 384 putative SNPs assayed, we successfully validated and genotyped 362 (94.3%). Of these 325 (84.6%) showed polymorphism among the 91 rice genotypes examined. Physical distribution, degree of allele sharing, admixtures and introgression, and amino acid replacement of SNPs in 263 abiotic and 62 biotic stress-responsive genes provided clues for identification and targeted mapping of trait-associated genomic regions. We assessed the functional and adaptive significance of validated SNPs in a set of contrasting drought tolerant upland and sensitive lowland rice genotypes by correlating their allelic variation with amino acid sequence alterations in catalytic domains and three-dimensional secondary protein structure encoded by stress-responsive genes. We found a strong genetic association among SNPs in the nine stress-responsive genes with upland and lowland ecological adaptation. Higher nucleotide diversity was observed in indica accessions compared with other rice sub-populations based on different population genetic parameters. The inferred ancestry of 16% among rice genotypes was derived from admixed populations with the maximum between upland aus and wild Oryza species.

Conclusions

SNPs validated in biotic and abiotic stress-responsive rice genes can be used in association analyses to identify candidate genes and develop functional markers for stress tolerance in rice.

Dynamic intra-japonica subspecies variation and resource application

We constructed a physical map of O. sativa ssp. japonica cv. ZH11 and compared it and its random sample sequences with the Nipponbare RefSeq derived from the same subspecies. This comparison showed that the two japonica genomes were highly syntenic but revealed substantial differences in terms of structural variations, rates of substitutions and indels, and transposable element content. For example, contractions/expansions as large as 450 kb and repeat sequences that were present in high copy numbers only in ZH11 were detected. In tri-alignment regions using the indica variety 93-11 sequence as an outgroup, we found that: (1) the substitution rates of the two japonica-indica inter-subspecies comparison combinations were close but almost a magnitude higher than the substitution rate between the japonica rice varieties ZH11 and Nipponbare; (2) of the substitutions found between ZH11 and Nipponbare, 47.2% occurred in ZH11 and 52.6% in Nipponbare; (3) of the indels found between ZH11 and Nipponbare, the indels that occurred in ZH11 were 15.8 times of those in Nipponbare. Of the indels that occurred in ZH11, 75.67% were insertions and 24.33% deletions. Of the indels that occurred in Nipponbare, 48.23% were insertions and 51.77% were deletions. The ZH11 comparative map covered four Nipponbare physical gaps, detected assembly errors in the Nipponbare sequence, and was integrated with the FSTs of a large ZH11 T-DNA insertion mutant library. ZH11 BAC clones can be browsed, searched, and obtained at our website, http://GResource.hzau.edu.cn.

Genetic architecture of aluminum tolerance in rice (Oryza sativa) determined through genome-wide association analysis and QTL mapping

Aluminum (Al) toxicity is a primary limitation to crop productivity on acid soils, and rice has been demonstrated to be significantly more Al tolerant than other cereal crops. However, the mechanisms of rice Al tolerance are largely unknown, and no genes underlying natural variation have been reported. We screened 383 diverse rice accessions, conducted a genome-wide association (GWA) study, and conducted QTL mapping in two bi-parental populations using three estimates of Al tolerance based on root growth. Subpopulation structure explained 57% of the phenotypic variation, and the mean Al tolerance in Japonica was twice that of Indica. Forty-eight regions associated with Al tolerance were identified by GWA analysis, most of which were subpopulation-specific. Four of these regions co-localized with a priori candidate genes, and two highly significant regions co-localized with previously identified QTLs. Three regions corresponding to induced Al-sensitive rice mutants (ART1, STAR2, Nrat1) were identified through bi-parental QTL mapping or GWA to be involved in natural variation for Al tolerance. Haplotype analysis around the Nrat1 gene identified susceptible and tolerant haplotypes explaining 40% of the Al tolerance variation within the aus subpopulation, and sequence analysis of Nrat1 identified a trio of non-synonymous mutations predictive of Al sensitivity in our diversity panel. GWA analysis discovered more phenotype-genotype associations and provided higher resolution, but QTL mapping identified critical rare and/or subpopulation-specific alleles not detected by GWA analysis. Mapping using Indica/Japonica populations identified QTLs associated with transgressive variation where alleles from a susceptible aus or indica parent enhanced Al tolerance in a tolerant Japonica background. This work supports the hypothesis that selectively introgressing alleles across subpopulations is an efficient approach for trait enhancement in plant breeding programs and demonstrates the fundamental importance of subpopulation in interpreting and manipulating the genetics of complex traits in rice.

Complex evolution of S5, a major reproductive barrier regulator, in the cultivated rice Oryza sativa and its wild relatives

• The hybrid sterility gene S5 comprises three types of alleles in cultivated rice. Such tri-allelic system provided a unique opportunity to study the molecular bases of evolutionary changes underlying reproductive isolation in plants. • We analysed the sequence diversity and evolutionary history of S5 in 138 Oryza accessions. We also examined the effect of the two functional variations (C819A and C1412T) in determining hybrid sterility by transformation. • Nineteen haplotypes were identified, which were classified into the indica-like, the japonica-like and the wide-compatibility gene (WCG)-like group, according to the sequence features of the tri-allelic system. The origin and evolutionary course of the three allelic groups were investigated, thus confirming the independent origins of indica and japonica subspecies. There were perfect associations between C819A and C1412T in the rice germplasm assayed, and the combination of C819 and C1412 was required for hybrid sterility. Evidence of positive selection in the WCG-like alleles suggested that they might have been favored by selection for higher compatibility in hybrids. • The complex evolution of S5 revealed the counteractive function of the three allelic groups at the species level. S5 might perform an important primary function in an evolutionary scale, and hybrid sterility acts as a 'byproduct' of this speciation gene.

Exceptional lability of a genomic complex in rice and its close relatives revealed by interspecific and intraspecific comparison and population analysis

Background

Extensive DNA rearrangement of genic colinearity, as revealed by comparison of orthologous genomic regions, has been shown to be a general concept describing evolutionary dynamics of plant genomes. However, the nature, timing, lineages and adaptation of local genomic rearrangement in closely related species (e.g., within a genus) and haplotype variation of genomic rearrangement within populations have not been well documented.

Results

We previously identified a hotspot for genic rearrangement and transposon accumulation in the Orp region of Asian rice (Oryza sativa, AA) by comparison with its orthologous region in sorghum. Here, we report the comparative analysis of this region with its orthologous regions in the wild progenitor species (O. nivara, AA) of Asian rice and African rice (O. glaberrima) using the BB genome Oryza species (O. punctata) as an outgroup, and investigation of transposon insertion sites and a segmental inversion event in the AA genomes at the population level. We found that Orp region was primarily and recently expanded in the Asian rice species O. sativa and O. nivara. LTR-retrotransposons shared by the three AA-genomic regions have been fixed in all the 94 varieties that represent different populations of the AA-genome species/subspecies, indicating their adaptive role in genome differentiation. However, LTR-retrotransposons unique to either O. nivara or O. sativa regions exhibited dramatic haplotype variation regarding their presence or absence between or within populations/subpopulations.

Conclusions

The LTR-retrotransposon insertion hotspot in the Orp region was formed recently, independently and concurrently in different AA-genome species, and that the genic rearrangements detected in different species appear to be differentially triggered by transposable elements. This region is located near the end of the short arm of chromosome 8 and contains a high proportion of LTR-retrotransposons similar to observed in the centromeric region of this same chromosome, and thus may represent a genomic region that has recently switched from euchromatic to heterochromatic states. The haplotype variation of LTR-retrotransposon insertions within this region reveals substantial admixture among various subpopulations as established by molecular markers at the whole genome level, and can be used to develop retrotransposon junction markers for simple and rapid classification of O. sativa germplasm.

A global analysis of QTLs for expression variations in rice shoots at the early seedling stage

Analyses of quantitative trait loci (QTLs) for expression levels (eQTLs) of genes reveal a genetic relationship between expression variation and the regulator, thus unlocking information for identifying the regulatory network. Oligo-nucleotide expression microarrays hybridized with RNA can simultaneously provide data for molecular markers and transcript abundance. In this study, we used an Affymetrix GeneChip Rice Genome Array to analyze eQTLs in rice shoots at 72 h after germination from 110 recombinant inbred lines (RILs) derived from a cross between Zhenshan 97 and Minghui 63. In total, 1632 single-feature polymorphisms (SFPs) plus 23 PCR markers were identified and placed into 601 recombinant bins, spanning 1459 cM in length, which were used as markers to genotype the RILs. We obtained 16,372 expression traits (e-traits) each with at least one eQTL, resulting in 26,051 eQTLs in total, including both cis- and trans-eQTLs. We also identified 171 eQTL hot spots in the rice genome, each of which controls transcript variations of many e-traits. Gene ontology analysis revealed an enrichment of certain functional categories of genes in some of the eQTL hot spots. In particular, eQTLs for e-traits involving the DNA metabolic process was significantly enriched in several eQTL hot spots on chromosomes 3, 5 and 10. Several e-traits co-localizing with cis-eQTLs showed significant correlations with hundreds of e-traits, indicating possible co-regulation. We also detected correlations between QTLs for shoot dry weight and eQTLs, revealing possible candidate genes for the trait. These results provided clues for the identification and characterization of the regulatory network in the whole genome at the transcriptional level.

Fine definition of the pedigree haplotypes of closely related rice cultivars by means of genome-wide discovery of single-nucleotide polymorphisms

Background

To create useful gene combinations in crop breeding, it is necessary to clarify the dynamics of the genome composition created by breeding practices. A large quantity of single-nucleotide polymorphism (SNP) data is required to permit discrimination of chromosome segments among modern cultivars, which are genetically related. Here, we used a high-throughput sequencer to conduct whole-genome sequencing of an elite Japanese rice cultivar, Koshihikari, which is closely related to Nipponbare, whose genome sequencing has been completed. Then we designed a high-throughput typing array based on the SNP information by comparison of the two sequences. Finally, we applied this array to analyze historical representative rice cultivars to understand the dynamics of their genome composition.

Results

The total 5.89-Gb sequence for Koshihikari, equivalent to 15.7× the entire rice genome, was mapped using the Pseudomolecules 4.0 database for Nipponbare. The resultant Koshihikari genome sequence corresponded to 80.1% of the Nipponbare sequence and led to the identification of 67 051 SNPs. A high-throughput typing array consisting of 1917 SNP sites distributed throughout the genome was designed to genotype 151 representative Japanese cultivars that have been grown during the past 150 years. We could identify the ancestral origin of the pedigree haplotypes in 60.9% of the Koshihikari genome and 18 consensus haplotype blocks which are inherited from traditional landraces to current improved varieties. Moreover, it was predicted that modern breeding practices have generally decreased genetic diversity

Conclusions

Detection of genome-wide SNPs by both high-throughput sequencer and typing array made it possible to evaluate genomic composition of genetically related rice varieties. With the aid of their pedigree information, we clarified the dynamics of chromosome recombination during the historical rice breeding process. We also found several genomic regions decreasing genetic diversity which might be caused by a recent human selection in rice breeding. The definition of pedigree haplotypes by means of genome-wide SNPs will facilitate next-generation breeding of rice and other crops.

Genomics and bioinformatics resources for crop improvement

Recent remarkable innovations in platforms for omics-based research and application development provide crucial resources to promote research in model and applied plant species. A combinatorial approach using multiple omics platforms and integration of their outcomes is now an effective strategy for clarifying molecular systems integral to improving plant productivity. Furthermore, promotion of comparative genomics among model and applied plants allows us to grasp the biological properties of each species and to accelerate gene discovery and functional analyses of genes. Bioinformatics platforms and their associated databases are also essential for the effective design of approaches making the best use of genomic resources, including resource integration. We review recent advances in research platforms and resources in plant omics together with related databases and advances in technology.

Dissection of a QTL reveals an adaptive, interacting gene complex associated with transgressive variation for flowering time in rice

A days to heading QTL (dth1.1) located on the short arm of rice chromosome 1 was sub-divided into eight sub-introgression lines (SILs) to analyze the genetic basis of transgressive variation for flowering time. Each SIL contained one or more introgression(s) from O. rufipogon in the genetic background of the elite Oryza sativa cultivar, Jefferson. Each introgression was defined at high resolution using molecular markers and those in the dth1.1 region were associated with the presence of one or more flowering time genes (GI, SOC1, FT-L8, EMF1, and PNZIP). SILs and controls were evaluated for flowering time under both short- and long-day growing conditions. Under short-day lengths, lines with introgressions carrying combinations of linked flowering time genes (GI/SOC1, SOC1/FT-L8, GI/SOC1/FT-L8 and EMF1/PNZIP) from the late parent, O. rufipogon, flowered earlier than the recurrent parent, Jefferson, while recombinant lines carrying smaller introgressions marked by the presence of GI, SOC1, EMF1 or PNZIP alone no longer flowered early. Under long-day length, lines carrying SOC1/FT-L8, SOC1 or PNZIP flowered early, while those carrying GI or EMF1 delayed flowering. Across all experiments and in the field, only SIL_SOC1/FT-L8 was consistently early. A preliminary yield evaluation indicated that the transgressive early flowering observed in several of the SILs was also associated with a measurable and positive effect on yield. These SILs represent a new source of variation that can be used in breeding programs to manipulate flowering time in rice cultivars without the reduction in yield that is often associated with early maturing phenotypes.

Gene expression profiles deciphering rice phenotypic variation between Nipponbare (Japonica) and 93-11 (Indica) during oxidative stress

Rice is a very important food staple that feeds more than half the world's population. Two major Asian cultivated rice (Oryza sativa L.) subspecies, japonica and indica, show significant phenotypic variation in their stress responses. However, the molecular mechanisms underlying this phenotypic variation are still largely unknown. A common link among different stresses is that they produce an oxidative burst and result in an increase of reactive oxygen species (ROS). In this study, methyl viologen (MV) as a ROS agent was applied to investigate the rice oxidative stress response. We observed that 93-11 (indica) seedlings exhibited leaf senescence with severe lesions under MV treatment compared to Nipponbare (japonica). Whole-genome microarray experiments were conducted, and 1,062 probe sets were identified with gene expression level polymorphisms between the two rice cultivars in addition to differential expression under MV treatment, which were assigned as Core Intersectional Probesets (CIPs). These CIPs were analyzed by gene ontology (GO) and highlighted with enrichment GO terms related to toxin and oxidative stress responses as well as other responses. These GO term-enriched genes of the CIPs include glutathine S-transferases (GSTs), P450, plant defense genes, and secondary metabolism related genes such as chalcone synthase (CHS). Further insertion/deletion (InDel) and regulatory element analyses for these identified CIPs suggested that there may be some eQTL hotspots related to oxidative stress in the rice genome, such as GST genes encoded on chromosome 10. In addition, we identified a group of marker genes individuating the japonica and indica subspecies. In summary, we developed a new strategy combining biological experiments and data mining to study the possible molecular mechanism of phenotypic variation during oxidative stress between Nipponbare and 93-11. This study will aid in the analysis of the molecular basis of quantitative traits.

Do genetic recombination and gene density shape the pattern of DNA elimination in rice long terminal repeat retrotransposons?

In flowering plants, the accumulation of small deletions through unequal homologous recombination (UR) and illegitimate recombination (IR) is proposed to be the major process counteracting genome expansion, which is caused primarily by the periodic amplification of long terminal repeat retrotransposons (LTR-RTs). However, the full suite of evolutionary forces that govern the gain or loss of transposable elements (TEs) and their distribution within a genome remains unclear. Here, we investigated the distribution and structural variation of LTR-RTs in relation to the rates of local genetic recombination (GR) and gene densities in the rice (Oryza sativa) genome. Our data revealed a positive correlation between GR rates and gene densities and negative correlations between LTR-RT densities and both GR and gene densities. The data also indicate a tendency for LTR-RT elements and fragments to be shorter in regions with higher GR rates; the size reduction of LTR-RTs appears to be achieved primarily through solo LTR formation by UR. Comparison of indica and japonica rice revealed patterns and frequencies of LTR-RT gain and loss within different evolutionary timeframes. Different LTR-RT families exhibited variable distribution patterns and structural changes, but overall LTR-RT compositions and genes were organized according to the GR gradients of the genome. Further investigation of non-LTR-RTs and DNA transposons revealed a negative correlation between gene densities and the abundance of DNA transposons and a weak correlation between GR rates and the abundance of long interspersed nuclear elements (LINEs)/short interspersed nuclear elements (SINEs). Together, these observations suggest that GR and gene density play important roles in shaping the dynamic structure of the rice genome.

SNP genotyping by DNA photoligation: application to SNP detection of genes from food crops

We describe a simple and inexpensive single-nucleotide polymorphism (SNP) typing method, using DNA photoligation with 5-carboxyvinyl-2'-deoxyuridine and two fluorophores. This SNP-typing method facilitates qualitative determination of genes from indica and japonica rice, and showed a high degree of single nucleotide specificity up to 10 000. This method can be used in the SNP typing of actual genomic DNA samples from food crops.

Analysis of intraspecies diversity in wheat and barley genomes identifies breakpoints of ancient haplotypes and provides insight into the structure of diploid and hexaploid triticeae gene pools

A large number of wheat (Triticum aestivum) and barley (Hordeum vulgare) varieties have evolved in agricultural ecosystems since domestication. Because of the large, repetitive genomes of these Triticeae crops, sequence information is limited and molecular differences between modern varieties are poorly understood. To study intraspecies genomic diversity, we compared large genomic sequences at the Lr34 locus of the wheat varieties Chinese Spring, Renan, and Glenlea, and diploid wheat Aegilops tauschii. Additionally, we compared the barley loci Vrs1 and Rym4 of the varieties Morex, Cebada Capa, and Haruna Nijo. Molecular dating showed that the wheat D genome haplotypes diverged only a few thousand years ago, while some barley and Ae. tauschii haplotypes diverged more than 500,000 years ago. This suggests gene flow from wild barley relatives after domestication, whereas this was rare or absent in the D genome of hexaploid wheat. In some segments, the compared haplotypes were very similar to each other, but for two varieties each at the Rym4 and Lr34 loci, sequence conservation showed a breakpoint that separates a highly conserved from a less conserved segment. We interpret this as recombination breakpoints of two ancient haplotypes, indicating that the Triticeae genomes are a heterogeneous and variable mosaic of haplotype fragments. Analysis of insertions and deletions showed that large events caused by transposable element insertions, illegitimate recombination, or unequal crossing over were relatively rare. Most insertions and deletions were small and caused by template slippage in short homopolymers of only a few base pairs in size. Such frequent polymorphisms could be exploited for future molecular marker development.

RICD: a rice indica cDNA database resource for rice functional genomics

BACKGROUND

The Oryza sativa L. indica subspecies is the most widely cultivated rice. During the last few years, we have collected over 20,000 putative full-length cDNAs and over 40,000 ESTs isolated from various cDNA libraries of two indica varieties Guangluai 4 and Minghui 63. A database of the rice indica cDNAs was therefore built to provide a comprehensive web data source for searching and retrieving the indica cDNA clones.

RESULTS

Rice Indica cDNA Database (RICD) is an online MySQL-PHP driven database with a user-friendly web interface. It allows investigators to query the cDNA clones by keyword, genome position, nucleotide or protein sequence, and putative function. It also provides a series of information, including sequences, protein domain annotations, similarity search results, SNPs and InDels information, and hyperlinks to gene annotation in both The Rice Annotation Project Database (RAP-DB) and The TIGR Rice Genome Annotation Resource, expression atlas in RiceGE and variation report in Gramene of each cDNA.

CONCLUSION

The online rice indica cDNA database provides cDNA resource with comprehensive information to researchers for functional analysis of indica subspecies and for comparative genomics. The RICD database is available through our website http://www.ncgr.ac.cn/ricd.

Genome-wide analysis of transposon insertion polymorphisms reveals intraspecific variation in cultivated rice

Insertions and precise eliminations of transposable elements generated numerous transposon insertion polymorphisms (TIPs) in rice (Oryza sativa). We observed that TIPs represent more than 50% of large insertions and deletions (>100 bp) in the rice genome. Using a comparative genomic approach, we identified 2,041 TIPs between the genomes of two cultivars, japonica Nipponbare and indica 93-11. We also identified 691 TIPs between Nipponbare and indica Guangluai 4 in the 23-Mb collinear regions of chromosome 4. Among them, retrotransposon-based insertion polymorphisms were used to reveal the evolutionary relationships of these three cultivars. Our conservative estimates suggest that the TIPs generated approximately 14% of the genomic DNA sequence differences between subspecies indica and japonica. It was also found that more than 10% of TIPs were located in expressed gene regions, representing an important source of genetic variation. Transcript evidence implies that these TIPs induced a series of genetic differences between two subspecies, including interrupting host genes, creating different expression forms, drastically changing intron length, and affecting expression levels of adjacent genes. These analyses provide genome-wide insights into evolutionary history and genetic variation of rice.

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.

Rice yellow mottle virus stress responsive genes from susceptible and tolerant rice genotypes

BACKGROUND

The effects of viral infection involve concomitant plant gene variations and cellular changes. A simple system is required to assess the complexity of host responses to viral infection. The genome of the Rice yellow mottle virus (RYMV) is a single-stranded RNA with a simple organisation. It is the most well-known monocotyledon virus model. Several studies on its biology, structure and phylogeography have provided a suitable background for further genetic studies. 12 rice chromosome sequences are now available and provide strong support for genomic studies, particularly physical mapping and gene identification.

RESULTS

The present data, obtained through the cDNA-AFLP technique, demonstrate differential responses to RYMV of two different rice cultivars, i.e. susceptible IR64 (Oryza sativa indica), and partially resistant Azucena (O. s. japonica). This RNA profiling provides a new original dataset that will enable us to gain greater insight into the RYMV/rice interaction and the specificity of the host response. Using the SIM4 subroutine, we took the intron/exon structure of the gene into account and mapped 281 RYMV stress responsive (RSR) transcripts on 12 rice chromosomes corresponding to 234 RSR genes. We also mapped previously identified deregulated proteins and genes involved in partial resistance and thus constructed the first global physical map of the RYMV/rice interaction. RSR transcripts on rice chromosomes 4 and 10 were found to be not randomly distributed. Seven genes were identified in the susceptible and partially resistant cultivars, and transcripts were colocalized for these seven genes in both cultivars. During virus infection, many concomitant plant gene expression changes may be associated with host changes caused by the infection process, general stress or defence responses. We noted that some genes (e.g. ABC transporters) were regulated throughout the kinetics of infection and differentiated susceptible and partially resistant hosts.

CONCLUSION

We enhanced the first RYMV/rice interaction map by combining information from the present study and previous studies on proteins and ESTs regulated during RYMV infection, thus providing a more comprehensive view on genes related to plant responses. This combined map provides a new tool for exploring molecular mechanisms underlying the RYMV/rice interaction.

Case study for identification of potentially indel-caused alternative expression isoforms in the rice subspecies japonica and indica by integrative genome analysis

Alternative splicing (AS) is one of the most significant components of the functional complexity of the eukaryote genome, increasing protein diversity, creating isoforms, and affecting mRNA stability. Recently, whole genome sequences and large microarray data sets have become available, making data integration feasible and allowing the study of the possible regulatory mechanism of AS in rice (Oryza sativa) by erecting and testing hypotheses before doing bench studies. We have developed a new strategy and have identified 215 rice genes with alternative expression isoforms related to insertion and deletion (indel) between subspecies indica and subspecies japonica. We did a case study for alternative expression isoforms of the rice peroxidase gene LOC_Os06g48030 to investigate possible mechanisms by which indels caused alternative splicing between the indica and the japonica varieties by mining of array data together with validation by RT-PCR and genome sequencing analysis. Multiple poly(A) signals were detected in the specific indel region for LOC_Os06g48030. We present a new methodology to promote more discoveries of potentially indel-caused AS genes in rice, which may serve as the foundation for research into the regulatory mechanism of alternative expression isoforms between subspecies.

The number of genes having different alleles between rice cultivars estimated by SNP analysis

Identification of single nucleotide polymorphisms (SNPs) in a large number of genes will enable estimation of the number of genes having different alleles in a population. In the present study, SNPs between 21 rice cultivars including 17 Japanese cultivars, one upland rice, and three indica cultivars were analyzed by PCR-RF-SSCP. PCR-RF-SSCP analysis was found to be a more efficient method for detecting SNPs than mismatch-cleavage analysis, though both PCR-RF-SSCP and mismatch-cleavage are useful for screening SNPs. The number of DNA fragments showing polymorphism between Japanese cultivars was 134 in the 1,036 genes analyzed. In 137 genes, 638 DNA polymorphisms were identified. Out of 52 genes having polymorphisms in the exons, one had a frame-shift mutation, three had polymorphism causing amino acid insertions or deletions, and 16 genes had missense polymorphisms. The number of genes having frame-shift mutations and missense polymorphisms between the 17 Japanese cultivars was estimated to be 41 and 677 on average, respectively, and those between japonica and indica to be 425 and 6,977, respectively. Chromosomal regions of cultivars selected in rice breeding processes were identified by SNP analysis of genes.

A collection of 10,096 indica rice full-length cDNAs reveals highly expressed sequence divergence between Oryza sativa indica and japonica subspecies

Relatively few indica rice full-length cDNAs were available to aid in the annotation of rice genes. The data presented here described the sequencing and analysis of 10,096 full-length cDNAs from Oryza sativa subspecies indica Guangluai 4. Of them, 9,029 matched rice genomic sequences in publicly-available databases, and 1,200 were identified as new rice genes. Comparison with the knowledge-based Oryza Molecular Biological Encyclopedia japonica cDNA collection indicated that 3,316 (41.6%) of the 7,965 indica-japonica cDNA pairs showed no distinct variations at protein level (2,117 indica-japonica cDNA pairs showed fully identical and 1,199 indica-japonica cDNA pairs showed no frame shift). Moreover, 3,645 (45.8%) of the indica-japonica pairs showed substantial differences at the protein level due to single nucleotide polymorphisms (SNPs), insertions or deletions, and sequence-segment variations between indica and japonica subspecies. Further experimental verifications using PCR screening and quantitative reverse transcriptional PCR revealed unique transcripts for indica subspecies. Comparative analysis also showed that most of rice genes were evolved under purifying selection. These variations might distinguish the phenotypic changes of the two cultivated rice subspecies indica and japonica. Analysis of these cDNAs extends known rice genes and identifies new ones in rice.

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.

Rice functional genomics research in China

Rice functional genomics is a scientific approach that seeks to identify and define the function of rice genes, and uncover when and how genes work together to produce phenotypic traits. Rapid progress in rice genome sequencing has facilitated research in rice functional genomics in China. The Ministry of Science and Technology of China has funded two major rice functional genomics research programmes for building up the infrastructures of the functional genomics study such as developing rice functional genomics tools and resources. The programmes were also aimed at cloning and functional analyses of a number of genes controlling important agronomic traits from rice. National and international collaborations on rice functional genomics study are accelerating rice gene discovery and application.

Highly asymmetric rice genomes

Background

Individuals in the same species are assumed to share the same genomic set. However, it is not unusual to find an orthologous gene only in small subset of the species, and recent genomic studies suggest that structural rearrangements are very frequent between genomes in the same species. Two recently sequenced rice genomes Oryza sativa L. var. Nipponbare and O. sativa L. var. 93-11 provide an opportunity to systematically investigate the extent of the gene repertoire polymorphism, even though the genomic data of 93-11 derived from whole-short-gun sequencing is not yet as complete as that of Nipponbare.

Results

We compared gene contents and the genomic locations between two rice genomes. Our conservative estimates suggest that at least 10% of the genes in the genomes were either under presence/absence polymorphism (5.2%) or asymmetrically located between genomes (4.7%). The proportion of these "asymmetric genes" varied largely among gene groups, in which disease resistance (R) genes and the RLK kinase gene group had 11.6 and 7.8 times higher proportion of asymmetric genes than housekeeping genes (Myb and MADS). The significant difference in the proportion of asymmetric genes among gene groups suggests that natural selection is responsible for maintaining genomic asymmetry. On the other hand, the nucleotide diversity in 17 R genes under presence/absence polymorphism was generally low (average nucleotide diversity = 0.0051).

Conclusion

The genomic symmetry was disrupted by 10% of asymmetric genes, which could cause genetic variation through more unequal crossing over, because these genes had no allelic counterparts to pair and then they were free to pair with homologues at non-allelic loci, during meiosis in heterozygotes. It might be a consequence of diversifying selection that increased the structural divergence among genomes, and of purifying selection that decreased nucleotide divergence in each R gene locus.

Comparative physical mapping between Oryza sativa (AA genome type) and O. punctata (BB genome type)

A comparative physical map of the AA genome (Oryza sativa) and the BB genome (O. punctata) was constructed by aligning a physical map of O. punctata, deduced from 63,942 BAC end sequences (BESs) and 34,224 fingerprints, onto the O. sativa genome sequence. The level of conservation of each chromosome between the two species was determined by calculating a ratio of BES alignments. The alignment result suggests more divergence of intergenic and repeat regions in comparison to gene-rich regions. Further, this characteristic enabled localization of heterochromatic and euchromatic regions for each chromosome of both species. The alignment identified 16 locations containing expansions, contractions, inversions, and transpositions. By aligning 40% of the punctata BES on the map, 87% of the punctata FPC map covered 98% of the O. sativa genome sequence. The genome size of O. punctata was estimated to be 8% larger than that of O. sativa with individual chromosome differences of 1.5-16.5%. The sum of expansions and contractions observed in regions >500 kb were similar, suggesting that most of the contractions/expansions contributing to the genome size difference between the two species are small, thus preserving the macro-collinearity between these species, which diverged approximately 2 million years ago.

Biased distribution of microsatellite motifs in the rice genome

Microsatellites are useful tools to study the extent of divergence between two taxonomic groups that show high sequence similarity. We have compared microsatellite distribution to illustrate genetic variation between the two rice genomes, Oryza sativa L. ssp. indica and Oryza sativa L. ssp. japonica. Microsatellite distribution proved to be non random as certain regions of very high microsatellite density have been identified. Microsatellite density in the subspecies japonica was computed marginally higher than in the subspecies indica in the genomic regions compared between the two subspecies. Unexpectedly high microsatellite densities were observed in 5'-untranslated regions of genes. These regions also displayed a clear motif bias. Some of the longest microsatellite repeats were found in intron sequences. Frequency, as well as motif bias was also noted with respect to the association of microsatellites with transposable elements. Microsatellite mutability values were exemplarily estimated for 90 loci by aligning the microsatellite containing regions between the two genomes. Poor rates of finding an orthologue corresponded with high microsatellite mutability in rice. These insights are likely to play a significant role in selecting microsatellite loci to be used in molecular breeding and studying evolutionary dynamics of the two subspecies.

Differentiation of the two rice subspecies indica and japonica: a Gene Ontology perspective

Cultivated rice (Oryza sativa) is comprised of two subspecies: japonica and indica. Polymorphism levels between putative homologues were determined for genes whose japonica homologue had been classified into functional categories using the Gene Ontology (GO) system. Genes were partitioned into below-average and above-average polymorphism groups, and then the set of genes having each GO term was checked for the randomness of its distribution into these polymorphism groups using a series of False Discovery Rate (FDR) tests. The robustness of the conclusions was enhanced by employing different cutoff values and sequence samplings in the FDR tests. Significant nonrandom polymorphism distributions were found for protein-coding sequences in many GO categories. In contrast, a random distribution for nearly all GO terms was seen with intron sequences. These results were extended by measuring the nonsynonymous to synonymous codon usage ratio (dN/dS) using a permutation test, which showed that some above-average polymorphism GO categories also had a high proportion of genes with a dN/dS ratio greater than one, suggesting positive selection on these GO categories during indica-japonica differentiation. An analysis of predominant gene names in the significant GO categories divided them into four functional classes: production of defense-related compounds, cell wall, cell signaling, and transcription factors.

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.