A fine physical map of the rice chromosome 4

Twenty years of rice genomics research: From sequencing and functional genomics to quantitative genomics

Since the completion of the rice genome sequencing project in 2005, we have entered the era of rice genomics, which is still in its ascendancy. Rice genomics studies can be classified into three stages: structural genomics, functional genomics, and quantitative genomics. Structural genomics refers primarily to genome sequencing for the construction of a complete map of rice genome sequence. This is fundamental for rice genetics and molecular biology research. Functional genomics aims to decode the functions of rice genes. Quantitative genomics is large-scale sequence- and statistics-based research to define the quantitative traits and genetic features of rice populations. Rice genomics has been a transformative influence on rice biological research and contributes significantly to rice breeding, making rice a good model plant for studying crop sciences.

Chilling tolerance in rice: Past and present

Rice is generally sensitive to chilling stress, which seriously affects growth and yield. Since early in the last century, considerable efforts have been made to understand the physiological and molecular mechanisms underlying the response to chilling stress and improve rice chilling tolerance. Here, we review the research trends and advances in this field. The phenotypic and biochemical changes caused by cold stress and the physiological explanations are briefly summarized. Using published data from the past 20 years, we reviewed the past progress and important techniques in the identification of quantitative trait loci (QTL), novel genes, and cellular pathways involved in rice chilling tolerance. The advent of novel technologies has significantly advanced studies of cold tolerance, and the characterization of QTLs, key genes, and molecular modules have sped up molecular design breeding for cold tolerance in rice varieties. In addition to gene function studies based on overexpression or artificially generated mutants, elucidating natural allelic variation in specific backgrounds is emerging as a novel approach for the study of cold tolerance in rice, and the superior alleles identified using this approach can directly facilitate breeding.

Engineering meiotic recombination pathways in rice

In the last 15 years, outstanding progress has been made in understanding the function of meiotic genes in the model dicot and monocot plants Arabidopsis and rice (Oryza sativa L.), respectively. This knowledge allowed to modulate meiotic recombination in Arabidopsis and, more recently, in rice. For instance, the overall frequency of crossovers (COs) has been stimulated 2.3- and 3.2-fold through the inactivation of the rice FANCM and RECQ4 DNA helicases, respectively, two genes involved in the repair of DNA double-strand breaks (DSBs) as noncrossovers (NCOs) of the Class II crossover pathway. Differently, the programmed induction of DSBs and COs at desired sites is currently explored by guiding the SPO11-1 topoisomerase-like transesterase, initiating meiotic recombination in all eukaryotes, to specific target regions of the rice genome. Furthermore, the inactivation of 3 meiosis-specific genes, namely PAIR1, OsREC8 and OsOSD1, in the Mitosis instead of Meiosis (MiMe) mutant turned rice meiosis into mitosis, thereby abolishing recombination and achieving the first component of apomixis, apomeiosis. The successful translation of Arabidopsis results into a crop further allowed the implementation of two breakthrough strategies that triggered parthenogenesis from the MiMe unreduced clonal egg cell and completed the second component of diplosporous apomixis. Here, we review the most recent advances in and future prospects of the manipulation of meiotic recombination in rice and potentially other major crops, all essential for global food security.

Development of a quantitative pachytene chromosome map and its unification with somatic chromosome and linkage maps of rice (Oryza sativa L.)

A quantitative pachytene chromosome map of rice (Oryza sativa L.) was developed using imaging methods. The map depicts not only distribution patterns of chromomeres specific to pachytene chromosomes, but also the higher order information of chromosomal structures, such as heterochromatin (condensed regions), euchromatin (decondensed regions), the primary constrictions (centromeres), and the secondary constriction (nucleolar organizing regions, NOR). These features were image analyzed and quantitatively mapped onto the map by Chromosome Image Analyzing System ver. 4.0 (CHIAS IV). Correlation between H3K9me2, an epigenetic marker and formation and/or maintenance of heterochromatin, thus was, clearly visualized. Then the pachytene chromosome map was unified with the existing somatic chromosome and linkage maps by physically mapping common DNA markers among them, such as a rice A genome specific tandem repeat sequence (TrsA), 5S and 45S ribosomal RNA genes, five bacterial artificial chromosome (BAC) clones, four P1 bacteriophage artificial chromosome (PAC) clones using multicolor fluorescence in situ hybridization (FISH). Detailed comparison between the locations of the DNA probes on the pachytene chromosomes using multicolor FISH, and the linkage map enabled determination of the chromosome number and short/long arms of individual pachytene chromosomes using the chromosome number and arm assignment designated for the linkage map. As a result, the quantitative pachytene chromosome map was unified with two other major rice chromosome maps representing somatic prometaphase chromosomes and genetic linkages. In conclusion, the unification of the three rice maps serves as an indispensable basic information, not only for an in-depth comparison between genetic and chromosomal data, but also for practical breeding programs.

RelocaTE2: a high resolution transposable element insertion site mapping tool for population resequencing

Background

Transposable element (TE) polymorphisms are important components of population genetic variation. The functional impacts of TEs in gene regulation and generating genetic diversity have been observed in multiple species, but the frequency and magnitude of TE variation is under appreciated. Inexpensive and deep sequencing technology has made it affordable to apply population genetic methods to whole genomes with methods that identify single nucleotide and insertion/deletion polymorphisms. However, identifying TE polymorphisms, particularly transposition events or non-reference insertion sites can be challenging due to the repetitive nature of these sequences, which hamper both the sensitivity and specificity of analysis tools.

Methods

We have developed the tool RelocaTE2 for identification of TE insertion sites at high sensitivity and specificity. RelocaTE2 searches for known TE sequences in whole genome sequencing reads from second generation sequencing platforms such as Illumina. These sequence reads are used as seeds to pinpoint chromosome locations where TEs have transposed. RelocaTE2 detects target site duplication (TSD) of TE insertions allowing it to report TE polymorphism loci with single base pair precision.

Results and Discussion

The performance of RelocaTE2 is evaluated using both simulated and real sequence data. RelocaTE2 demonstrate high level of sensitivity and specificity, particularly when the sequence coverage is not shallow. In comparison to other tools tested, RelocaTE2 achieves the best balance between sensitivity and specificity. In particular, RelocaTE2 performs best in prediction of TSDs for TE insertions. Even in highly repetitive regions, such as those tested on rice chromosome 4, RelocaTE2 is able to report up to 95% of simulated TE insertions with less than 0.1% false positive rate using 10-fold genome coverage resequencing data. RelocaTE2 provides a robust solution to identify TE insertion sites and can be incorporated into analysis workflows in support of describing the complete genotype from light coverage genome sequencing.

Development of Oryza sativa L. by Oryza punctata Kotschy ex Steud. monosomic addition lines with high value traits by interspecific hybridization

This paper describes the development of monosomic alien addition and disomic introgression lines through a cross between autotetraploid indica rice and Oryza punctata toward tapping valuable traits for rice improvement. Oryza punctata is a distantly related wild Oryza species having BB genome with untapped genetic resources for rice improvement. Low crossability between the cultivated O. sativa and O. punctata restricts the success of transferring many desirable traits into cultivated rice. Artificially induced autotetraploids of an elite breeding line, IR31917-45-3-2, were produced and crossed with O. punctata. Allotriploid F1 plants were backcrossed to IR31917-45-3-2 and generated progenies with extra chromosomes from O. punctata. Twenty BC1F1 and 59 BC2F1 plants were produced with chromosome numbers ranging from 24 (2n) to 29 (2n + 5) and 2n (24) to 26 (2n + 2), respectively. Eleven monosomic alien addition lines (MAALs) were characterized morphologically and cytologically and designated as MAAL 1-12. MAALs were genotyped using O. punctata genome-specific molecular markers and detected chromosome segments inherited from O. punctata. O. punctata introgressions across all the chromosomes of O. sativa were identified except for chromosome 8. The most frequent introgressions were observed in chromosomes 4, 6, 10, and 11, which could be the recombination hotspots between A and B genomes. Some of the qualitative traits such as black hull, purple coleoptile base, purple stigma, long awn, and short grain size from O. punctata were inherited in some disomic introgression lines (DILs). Several DILs inherited genes from O. punctata conferring resistance to brown planthopper, green leafhopper, and diseases such as bacterial blight and blast. This is the first report on successful gene transfer from O. punctata into O. sativa.

Highly distinct chromosomal structures in cowpea (Vigna unguiculata), as revealed by molecular cytogenetic analysis

Cowpea (Vigna unguiculata (L.) Walp) is an important legume, particularly in developing countries. However, little is known about its genome or chromosome structure. We used molecular cytogenetics to characterize the structure of pachytene chromosomes to advance our knowledge of chromosome and genome organization of cowpea. Our data showed that cowpea has highly distinct chromosomal structures that are cytologically visible as brightly DAPI-stained heterochromatic regions. Analysis of the repetitive fraction of the cowpea genome present at centromeric and pericentromeric regions confirmed that two retrotransposons are major components of pericentromeric regions and that a 455-bp tandem repeat is found at seven out of 11 centromere pairs in cowpea. These repeats likely evolved after the divergence of cowpea from common bean and form chromosomal structure unique to cowpea. The integration of cowpea genetic and physical chromosome maps reveals potential regions of suppressed recombination due to condensed heterochromatin and a lack of pairing in a few chromosomal termini. This study provides fundamental knowledge on cowpea chromosome structure and molecular cytogenetics tools for further chromosome studies.

Construction of high resolution genetic linkage maps to improve the soybean genome sequence assembly Glyma1.01

BACKGROUND

A landmark in soybean research, Glyma1.01, the first whole genome sequence of variety Williams 82 (Glycine max L. Merr.) was completed in 2010 and is widely used. However, because the assembly was primarily built based on the linkage maps constructed with a limited number of markers and recombinant inbred lines (RILs), the assembled sequence, especially in some genomic regions with sparse numbers of anchoring markers, needs to be improved. Molecular markers are being used by researchers in the soybean community, however, with the updating of the Glyma1.01 build based on the high-resolution linkage maps resulting from this research, the genome positions of these markers need to be mapped.

RESULTS

Two high density genetic linkage maps were constructed based on 21,478 single nucleotide polymorphism loci mapped in the Williams 82 x G. soja (Sieb. & Zucc.) PI479752 population with 1083 RILs and 11,922 loci mapped in the Essex x Williams 82 population with 922 RILs. There were 37 regions or single markers where marker order in the two populations was in agreement but was not consistent with the physical position in the Glyma1.01 build. In addition, 28 previously unanchored scaffolds were positioned. Map data were used to identify false joins in the Glyma1.01 assembly and the corresponding scaffolds were broken and reassembled to the new assembly, Wm82.a2.v1. Based upon the plots of the genetic on physical distance of the loci, the euchromatic and heterochromatic regions along each chromosome in the new assembly were delimited. Genomic positions of the commonly used markers contained in BARCSOYSSR_1.0 database and the SoySNP50K BeadChip were updated based upon the Wm82.a2.v1 assembly.

CONCLUSIONS

The information will facilitate the study of recombination hot spots in the soybean genome, identification of genes or quantitative trait loci controlling yield, seed quality and resistance to biotic or abiotic stresses as well as other genetic or genomic research.

BAC and RNA sequencing reveal the brown planthopper resistance gene BPH15 in a recombination cold spot that mediates a unique defense mechanism

BACKGROUND

Brown planthopper (BPH, Nilaparvata lugens Stål), is the most destructive phloem-feeding insect pest of rice (Oryza sativa). The BPH-resistance gene BPH15 has been proved to be effective in controlling the pest and widely applied in rice breeding programs. Nevertheless, molecular mechanism of the resistance remain unclear. In this study, we narrowed down the position of BPH15 on chromosome 4 and investigated the transcriptome of BPH15 rice after BPH attacked.

RESULTS

We analyzed 13,000 BC2F2 plants of cross between susceptible rice TN1 and the recombinant inbred line RI93 that carrying the BPH15 gene from original resistant donor B5. BPH15 was mapped to a 0.0269 cM region on chromosome 4, which is 210-kb in the reference genome of Nipponbare. Sequencing bacterial artificial chromosome (BAC) clones that span the BPH15 region revealed that the physical size of BPH15 region in resistant rice B5 is 580-kb, much bigger than the corresponding region in the reference genome of Nipponbare. There were 87 predicted genes in the BPH15 region in resistant rice. The expression profiles of predicted genes were analyzed. Four jacalin-related lectin proteins genes and one LRR protein gene were found constitutively expressed in resistant parent and considered the candidate genes of BPH15. The transcriptomes of resistant BPH15 introgression line and the susceptible recipient line were analyzed using high-throughput RNA sequencing. In total, 2,914 differentially expressed genes (DEGs) were identified. BPH-responsive transcript profiles were distinct between resistant and susceptible plants and between the early stage (6 h after infestation, HAI) and late stage (48 HAI). The key defense mechanism was related to jasmonate signaling, ethylene signaling, receptor kinase, MAPK cascades, Ca(2+) signaling, PR genes, transcription factors, and protein posttranslational modifications.

CONCLUSIONS

Our work combined BAC and RNA sequencing to identify candidate genes of BPH15 and revealed the resistance mechanism that it mediated. These results increase our understanding of plant-insect interactions and can be used to protect against this destructive agricultural pest.

BAC and RNA sequencing reveal the brown planthopper resistance gene BPH15 in a recombination cold spot that mediates a unique defense mechanism

BACKGROUND

Brown planthopper (BPH, Nilaparvata lugens Stål), is the most destructive phloem-feeding insect pest of rice (Oryza sativa). The BPH-resistance gene BPH15 has been proved to be effective in controlling the pest and widely applied in rice breeding programs. Nevertheless, molecular mechanism of the resistance remain unclear. In this study, we narrowed down the position of BPH15 on chromosome 4 and investigated the transcriptome of BPH15 rice after BPH attacked.

RESULTS

We analyzed 13,000 BC2F2 plants of cross between susceptible rice TN1 and the recombinant inbred line RI93 that carrying the BPH15 gene from original resistant donor B5. BPH15 was mapped to a 0.0269 cM region on chromosome 4, which is 210-kb in the reference genome of Nipponbare. Sequencing bacterial artificial chromosome (BAC) clones that span the BPH15 region revealed that the physical size of BPH15 region in resistant rice B5 is 580-kb, much bigger than the corresponding region in the reference genome of Nipponbare. There were 87 predicted genes in the BPH15 region in resistant rice. The expression profiles of predicted genes were analyzed. Four jacalin-related lectin proteins genes and one LRR protein gene were found constitutively expressed in resistant parent and considered the candidate genes of BPH15. The transcriptomes of resistant BPH15 introgression line and the susceptible recipient line were analyzed using high-throughput RNA sequencing. In total, 2,914 differentially expressed genes (DEGs) were identified. BPH-responsive transcript profiles were distinct between resistant and susceptible plants and between the early stage (6 h after infestation, HAI) and late stage (48 HAI). The key defense mechanism was related to jasmonate signaling, ethylene signaling, receptor kinase, MAPK cascades, Ca(2+) signaling, PR genes, transcription factors, and protein posttranslational modifications.

CONCLUSIONS

Our work combined BAC and RNA sequencing to identify candidate genes of BPH15 and revealed the resistance mechanism that it mediated. These results increase our understanding of plant-insect interactions and can be used to protect against this destructive agricultural pest.

A genome-wide association study of seed protein and oil content in soybean

BACKGROUND

Association analysis is an alternative to conventional family-based methods to detect the location of gene(s) or quantitative trait loci (QTL) and provides relatively high resolution in terms of defining the genome position of a gene or QTL. Seed protein and oil concentration are quantitative traits which are determined by the interaction among many genes with small to moderate genetic effects and their interaction with the environment. In this study, a genome-wide association study (GWAS) was performed to identify quantitative trait loci (QTL) controlling seed protein and oil concentration in 298 soybean germplasm accessions exhibiting a wide range of seed protein and oil content.

RESULTS

A total of 55,159 single nucleotide polymorphisms (SNPs) were genotyped using various methods including Illumina Infinium and GoldenGate assays and 31,954 markers with minor allele frequency >0.10 were used to estimate linkage disequilibrium (LD) in heterochromatic and euchromatic regions. In euchromatic regions, the mean LD (r2) rapidly declined to 0.2 within 360 Kbp, whereas the mean LD declined to 0.2 at 9,600 Kbp in heterochromatic regions. The GWAS results identified 40 SNPs in 17 different genomic regions significantly associated with seed protein. Of these, the five SNPs with the highest associations and seven adjacent SNPs were located in the 27.6-30.0 Mbp region of Gm20. A major seed protein QTL has been previously mapped to the same location and potential candidate genes have recently been identified in this region. The GWAS results also detected 25 SNPs in 13 different genomic regions associated with seed oil. Of these markers, seven SNPs had a significant association with both protein and oil.

CONCLUSIONS

This research indicated that GWAS not only identified most of the previously reported QTL controlling seed protein and oil, but also resulted in narrower genomic regions than the regions reported as containing these QTL. The narrower GWAS-defined genome regions will allow more precise marker-assisted allele selection and will expedite positional cloning of the causal gene(s).

Bridging the genotyping gap: using genotyping by sequencing (GBS) to add high-density SNP markers and new value to traditional bi-parental mapping and breeding populations

Genotyping by sequencing (GBS) is the latest application of next-generation sequencing protocols for the purposes of discovering and genotyping SNPs in a variety of crop species and populations. Unlike other high-density genotyping technologies which have mainly been applied to general interest "reference" genomes, the low cost of GBS makes it an attractive means of saturating mapping and breeding populations with a high density of SNP markers. One barrier to the widespread use of GBS has been the difficulty of the bioinformatics analysis as the approach is accompanied by a high number of erroneous SNP calls which are not easily diagnosed or corrected. In this study, we use a 384-plex GBS protocol to add 30,984 markers to an indica (IR64) × japonica (Azucena) mapping population consisting of 176 recombinant inbred lines of rice (Oryza sativa) and we release our imputation and error correction pipeline to address initial GBS data sparsity and error, and streamline the process of adding SNPs to RIL populations. Using the final imputed and corrected dataset of 30,984 markers, we were able to map recombination hot and cold spots and regions of segregation distortion across the genome with a high degree of accuracy, thus identifying regions of the genome containing putative sterility loci. We mapped QTL for leaf width and aluminum tolerance, and were able to identify additional QTL for both phenotypes when using the full set of 30,984 SNPs that were not identified using a subset of only 1,464 SNPs, including a previously unreported QTL for aluminum tolerance located directly within a recombination hotspot on chromosome 1. These results suggest that adding a high density of SNP markers to a mapping or breeding population through GBS has a great value for numerous applications in rice breeding and genetics research.

Development and evaluation of SoySNP50K, a high-density genotyping array for soybean

The objective of this research was to identify single nucleotide polymorphisms (SNPs) and to develop an Illumina Infinium BeadChip that contained over 50,000 SNPs from soybean (Glycine max L. Merr.). A total of 498,921,777 reads 35-45 bp in length were obtained from DNA sequence analysis of reduced representation libraries from several soybean accessions which included six cultivated and two wild soybean (G. soja Sieb. et Zucc.) genotypes. These reads were mapped to the soybean whole genome sequence and 209,903 SNPs were identified. After applying several filters, a total of 146,161 of the 209,903 SNPs were determined to be ideal candidates for Illumina Infinium II BeadChip design. To equalize the distance between selected SNPs, increase assay success rate, and minimize the number of SNPs with low minor allele frequency, an iteration algorithm based on a selection index was developed and used to select 60,800 SNPs for Infinium BeadChip design. Of the 60,800 SNPs, 50,701 were targeted to euchromatic regions and 10,000 to heterochromatic regions of the 20 soybean chromosomes. In addition, 99 SNPs were targeted to unanchored sequence scaffolds. Of the 60,800 SNPs, a total of 52,041 passed Illumina's manufacturing phase to produce the SoySNP50K iSelect BeadChip. Validation of the SoySNP50K chip with 96 landrace genotypes, 96 elite cultivars and 96 wild soybean accessions showed that 47,337 SNPs were polymorphic and generated successful SNP allele calls. In addition, 40,841 of the 47,337 SNPs (86%) had minor allele frequencies ≥ 10% among the landraces, elite cultivars and the wild soybean accessions. A total of 620 and 42 candidate regions which may be associated with domestication and recent selection were identified, respectively. The SoySNP50K iSelect SNP beadchip will be a powerful tool for characterizing soybean genetic diversity and linkage disequilibrium, and for constructing high resolution linkage maps to improve the soybean whole genome sequence assembly.

Ac-Ds solutions for rice insertion mutagenesis

Rice is the model plant for monocotyledons. Since the completion of the high-quality sequence of its genome, the international community is deploying efforts to identify the function of the 30-40,000 nontransposable element genes of rice. These efforts comprise the creation of large collections of rice mutants accessible to the international scientific community. In addition to loss of function mutants, insertion mutagenesis using Agrobacterium-mediated transformation and engineered mobile elements allows the identification of genes through enhancer or gene trapping or activation tagging. The maize transposable element Ac-Ds is known to be active in rice since the early 1990s and it does not interfere with endogenous rice transposons. This is the guaranty that induced mutation obtained with the mobility of the Ds element will be stable when the source of Ac transposase is removed from the mutated genome. In this paper, we describe single- or double-component T-DNA constructs that have been used to introduce a functional Ac-Ds system in rice and allowed the generation and selection of different type of Ds insertion mutations in the rice genome.

Genomic analysis of wild tomato introgressions determining metabolism- and yield-associated traits

With the aim of determining the genetic basis of metabolic regulation in tomato fruit, we constructed a detailed physical map of genomic regions spanning previously described metabolic quantitative trait loci of a Solanum pennellii introgression line population. Two genomic libraries from S. pennellii were screened with 104 colocated markers from five selected genomic regions, and a total of 614 bacterial artificial chromosome (BAC)/cosmids were identified as seed clones. Integration of sequence data with the genetic and physical maps of Solanum lycopersicum facilitated the anchoring of 374 of these BAC/cosmid clones. The analysis of this information resulted in a genome-wide map of a nondomesticated plant species and covers 10% of the physical distance of the selected regions corresponding to approximately 1% of the wild tomato genome. Comparative analyses revealed that S. pennellii and domesticated tomato genomes can be considered as largely colinear. A total of 1,238,705 bp from both BAC/cosmid ends and nine large insert clones were sequenced, annotated, and functionally categorized. The sequence data allowed the evaluation of the level of polymorphism between the wild and cultivated tomato species. An exhaustive microsynteny analysis allowed us to estimate the divergence date of S. pennellii and S. lycopersicum at 2.7 million years ago. The combined results serve as a reference for comparative studies both at the macrosyntenic and microsyntenic levels. They also provide a valuable tool for fine-mapping of quantitative trait loci in tomato. Furthermore, they will contribute to a deeper understanding of the regulatory factors underpinning metabolism and hence defining crop chemical composition.

Selection on grain shattering genes and rates of rice domestication

Molecular cloning of major quantitative trait loci (QTLs) responsible for the reduction of rice grain shattering, a hallmark of cereal domestication, provided opportunities for in-depth investigation of domestication processes. Here, we studied nucleotide variation at the shattering loci, sh4 and qSH1, for cultivated rice, Oryza sativa ssp. indica and Oryza sativa ssp. japonica, and the wild progenitors, Oryza nivara andOryza rufipogon. The nonshattering sh4 allele was fixed in all rice cultivars, with levels of sequence polymorphism significantly reduced in both indica and japonica cultivars relative to the wild progenitors. The sh4 phylogeny together with the neutrality tests and coalescent simulations suggested that sh4 had a single origin and was fixed by artificial selection during the domestication of rice. Selection on qSH1 was not detected in indica and remained unclear in japonica. Selection on sh4 could be strong enough to have driven its fixation in a population of cultivated rice within a period of c. 100 yr. The slow fixation of the nonshattering phenotype observed at the archeological sites might be a result of relatively weak selection on mutations other than sh4 in early rice cultivation. The fixation of sh4 could have been achieved later through strong selection for the optimal phenotype.

Assignment of genetic linkage maps to diploid Solanum tuberosum pachytene chromosomes by BAC-FISH technology

A cytogenetic map has been developed for diploid potato (Solanum tuberosum), in which the arms of the 12 potato bivalents can be identified in pachytene complements using multicolor fluorescence in situ hybridization (FISH) with a set of 60 genetically anchored bacterial artificial chromosome (BAC) clones from the RHPOTKEY BAC library. This diagnostic set of selected BACs (five per chromosome) hybridizes to euchromatic regions and corresponds to well-defined loci in the ultradense genetic map, and with these probes a new detailed and reliable pachytene karyotype could be established. Chromosome size has been estimated both from microscopic length measurements and from 4′,6-diamidino-2-phenylindole fluorescence-based DNA content measurements. In both approaches, chromosome 1 is the largest (100–115 Mb) and chromosome 11 the smallest (49–53 Mb). Detailed measurements of mega-base-pair to micrometer ratios have been obtained for chromosome 5, with average values of 1.07 Mb/μm for euchromatin and 3.67 Mb/μm for heterochromatin. In addition, our FISH results helped to solve two discrepancies in the potato genetic map related to chromosomes 8 and 12. Finally, we discuss the significance of the potato cytogenetic map for extended FISH studies in potato and related Solanaceae, which will be especially beneficial for the potato genome-sequencing project.

High potential of a transposon mPing as a marker system in japonica x japonica cross in rice

Although quantitative traits loci (QTL) analysis has been widely performed to isolate agronomically important genes, it has been difficult to obtain molecular markers between individuals with similar phenotypes (assortative mating). Recently, the miniature inverted-repeat transposable element mPing was shown to be active in the japonica strain Gimbozu EG4 where it had accumulated more than 1000 copies. In contrast, most other japonicas, including Nipponbare, have 50 or fewer mPing insertions in their genome. In this study we have exploited the polymorphism of mPing insertion sites to generate 150 PCR markers in a cross between the closely related japonicas, Nipponbare × Gimbozu (EG4). These new markers were distributed in genic regions of the whole genome and showed significantly higher polymorphism (150 of 183) than all other molecular markers tested including short sequence repeat markers (46 of 661). In addition, we performed QTL analysis with these markers using recombinant inbred lines derived from Nipponbare × Gimbozu EG4, and successfully mapped a locus involved in heading date on the short arm of chromosome 6. Moreover, we could easily map two novel loci involved in the culm length on the short arms of chromosomes 3 and 10.

High-resolution chromosome mapping of BACs using multi-colour FISH and pooled-BAC FISH as a backbone for sequencing tomato chromosome 6

Within the framework of the International Solanaceae Genome Project, the genome of tomato (Solanum lycopersicum) is currently being sequenced. We follow a 'BAC-by-BAC' approach that aims to deliver high-quality sequences of the euchromatin part of the tomato genome. BACs are selected from various libraries of the tomato genome on the basis of markers from the F2.2000 linkage map. Prior to sequencing, we validated the precise physical location of the selected BACs on the chromosomes by five-colour high-resolution fluorescent in situ hybridization (FISH) mapping. This paper describes the strategies and results of cytogenetic mapping for chromosome 6 using 75 seed BACs for FISH on pachytene complements. The cytogenetic map obtained showed discrepancies between the actual chromosomal positions of these BACs and their markers on the linkage group. These discrepancies were most notable in the pericentromere heterochromatin, thus confirming previously described suppression of cross-over recombination in that region. In a so called pooled-BAC FISH, we hybridized all seed BACs simultaneously and found a few large gaps in the euchromatin parts of the long arm that are still devoid of seed BACs and are too large for coverage by expanding BAC contigs. Combining FISH with pooled BACs and newly recruited seed BACs will thus aid in efficient targeting of novel seed BACs into these areas. Finally, we established the occurrence of repetitive DNA in heterochromatin/euchromatin borders by combining BAC FISH with hybridization of a labelled repetitive DNA fraction (Cot-100). This strategy provides an excellent means to establish the borders between euchromatin and heterochromatin in this chromosome.

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.

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.

Validation of rice genome sequence by optical mapping

Background

Rice feeds much of the world, and possesses the simplest genome analyzed to date within the grass family, making it an economically relevant model system for other cereal crops. Although the rice genome is sequenced, validation and gap closing efforts require purely independent means for accurate finishing of sequence build data.

Results

To facilitate ongoing sequencing finishing and validation efforts, we have constructed a whole-genome SwaI optical restriction map of the rice genome. The physical map consists of 14 contigs, covering 12 chromosomes, with a total genome size of 382.17 Mb; this value is about 11% smaller than original estimates. 9 of the 14 optical map contigs are without gaps, covering chromosomes 1, 2, 3, 4, 5, 7, 8 10, and 12 in their entirety – including centromeres and telomeres. Alignments between optical and in silico restriction maps constructed from IRGSP (International Rice Genome Sequencing Project) and TIGR (The Institute for Genomic Research) genome sequence sources are comprehensive and informative, evidenced by map coverage across virtually all published gaps, discovery of new ones, and characterization of sequence misassemblies; all totalling ~14 Mb. Furthermore, since optical maps are ordered restriction maps, identified discordances are pinpointed on a reliable physical scaffold providing an independent resource for closure of gaps and rectification of misassemblies.

Conclusion

Analysis of sequence and optical mapping data effectively validates genome sequence assemblies constructed from large, repeat-rich genomes. Given this conclusion we envision new applications of such single molecule analysis that will merge advantages offered by high-resolution optical maps with inexpensive, but short sequence reads generated by emerging sequencing platforms. Lastly, map construction techniques presented here points the way to new types of comparative genome analysis that would focus on discernment of structural differences revealed by optical maps constructed from a broad range of rice subspecies and varieties.

A plant-transformation-competent BIBAC/BAC-based map of rice for functional analysis and genetic engineering of its genomic sequence

Sequencing of the rice genome has provided a platform for functional genomics research of rice and other cereal species. However, multiple approaches are needed to determine the functions of its genes and sequences and to use the genome sequencing results for genetic improvement of cereal crops. Here, we report a plant-transformation-competent, binary bacterial artificial chromosome (BIBAC) and bacterial artificial chromosome (BAC) based map of rice to facilitate these studies. The map was constructed from 20 835 BIBAC and BAC clones, and consisted of 579 overlapping BIBAC/BAC contigs. To facilitate functional analysis of chromosome 8 genomic sequence and cloning of the genes and QTLs mapped to the chromosome, we anchored the chromosomal contigs to the existing rice genetic maps. The chromosomal map consists of 11 contigs, 59 genetic markers, and 36 sequence tagged sites, spanning a total of ca. 38 Mb in physical length. Comparative analysis between the genetic and physical maps of chromosome 8 showed that there are 3 "hot" and 2 "cold" spots of genetic recombination along the chromosomal arms in addition to the "cold spot" in the centromeric region, suggesting that the sequence component contents of a chromosome may affect its local genetic recombination frequencies. Because of its plant transformability, the BIBAC/BAC map could provide a platform for functional analysis of the rice genome sequence and effective use of the sequencing results for gene and QTL cloning and molecular breeding.

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.

Rice as a model for centromere and heterochromatin research

Rice (Oryza sativa) has become an important model plant species in numerous research projects involving genome, molecular and evolutionary biology. In this review we describe the reasons why rice provides an excellent model system for centromere and heterochromatin research. In most multicellular eukaryotes, centromeres and heterochromatic domains contain long arrays of repetitive DNA elements that are recalcitrant to DNA sequencing. In contrast, three rice centromeres and the majority of the cytologically defined heterochromatin in the rice genome have been sequenced to high quality, providing an unparalleled resource compared to other model multicellular eukaryotes. Most importantly, active genes have been discovered in the functional domains of several rice centromeres. The centromeric genes and sequence resources provide an unprecedented opportunity to study function and evolution of centromeres and centromere-associated genes.

A 1,000-loci transcript map of the barley genome: new anchoring points for integrative grass genomics

An integrated barley transcript map (consensus map) comprising 1,032 expressed sequence tag (EST)-based markers (total 1,055 loci: 607 RFLP, 190 SSR, and 258 SNP), and 200 anchor markers from previously published data, has been generated by mapping in three doubled haploid (DH) populations. Between 107 and 179 EST-based markers were allocated to the seven individual barley linkage groups. The map covers 1118.3 cM with individual linkage groups ranging from 130 cM (chromosome 4H) to 199 cM (chromosome 3H), yielding an average marker interval distance of 0.9 cM. 475 EST-based markers showed a syntenic organisation to known colinear linkage groups of the rice genome, providing an extended insight into the status of barley/rice genome colinearity as well as ancient genome duplications predating the divergence of rice and barley. The presented barley transcript map is a valuable resource for targeted marker saturation and identification of candidate genes at agronomically important loci. It provides new anchor points for detailed studies in comparative grass genomics and will support future attempts towards the integration of genetic and physical mapping information.

High-resolution single-copy gene fluorescence in situ hybridization and its use in the construction of a cytogenetic map of maize chromosome 9

High-resolution cytogenetic maps provide important biological information on genome organization and function, as they correlate genetic distance with cytological structures, and are an invaluable complement to physical sequence data. The most direct way to generate a cytogenetic map is to localize genetically mapped genes onto chromosomes by fluorescence in situ hybridization (FISH). Detection of single-copy genes on plant chromosomes has been difficult. In this study, we developed a squash FISH procedure allowing successful detection of single-copy genes on maize (Zea mays) pachytene chromosomes. Using this method, the shortest probe that can be detected is 3.1 kb, and two sequences separated by approximately 100 kb can be resolved. To show the robust nature of this protocol, we localized nine genetically mapped single-copy genes on chromosome 9 in one FISH experiment. Integration of existing information from genetic maps and the BAC contig-based physical map with the cytological structure of chromosome 9 provides a comprehensive cross-referenced cytogenetic map and shows the dramatic reduction of recombination in the pericentromeric heterochromatic region. To establish a feasible mapping system for maize, we also developed a probe cocktail for unambiguous identification of the 10 maize pachytene chromosomes. These results provide a starting point toward constructing a high-resolution integrated cytogenetic map of maize.

Caught red-handed: Rc encodes a basic helix-loop-helix protein conditioning red pericarp in rice

Rc is a domestication-related gene required for red pericarp in rice (Oryza sativa). The red grain color is ubiquitous among the wild ancestors of O. sativa, in which it is closely associated with seed shattering and dormancy. Rc encodes a basic helix-loop-helix (bHLH) protein that was fine-mapped to an 18.5-kb region on rice chromosome 7 using a cross between Oryza rufipogon (red pericarp) and O. sativa cv Jefferson (white pericarp). Sequencing of the alleles from both mapping parents as well as from two independent genetic stocks of Rc revealed that the dominant red allele differed from the recessive white allele by a 14-bp deletion within exon 6 that knocked out the bHLH domain of the protein. A premature stop codon was identified in the second mutant stock that had a light red pericarp. RT-PCR experiments confirmed that the Rc gene was expressed in both red- and white-grained rice but that a shortened transcript was present in white varieties. Phylogenetic analysis, supported by comparative mapping in rice and maize (Zea mays), showed that Rc, a positive regulator of proanthocyanidin, is orthologous with INTENSIFIER1, a negative regulator of anthocyanin production in maize, and is not in the same clade as rice bHLH anthocyanin regulators.

An active DNA transposon nDart causing leaf variegation and mutable dwarfism and its related elements in rice

While characterized mutable alleles caused by DNA transposons have been abundant in maize since the discovery of Dissociation conferring variegation by Barbara McClintock, only a few mutable alleles have been described in rice even though the rice genome contains various transposons. Here, we show that a spontaneous mutable virescent allele, pyl-v, is caused by the disruption of the nuclear-coded essential chloroplast protease gene, OsClpP5, due to insertion of a 607-bp non-autonomous DNA transposon, non-autonomous DNA-based active rice transposon one (nDart1), belonging to the hAT superfamily. The transposition of nDart1 can be induced by crossing with a line containing an autonomous element, aDart, and stabilized by segregating out of aDart. We also identified a novel mutable dwarf allele thl-m caused by an insertion of nDart1. The japonica cultivar Nipponbare carries no aDart, although it contains epigenetically silenced Dart element(s), which can be activated by 5-azacytidine. Nipponbare bears four subgroups of about 3.6-kb Dart-like sequences, three of which contain potential transposase genes, and around 3.6-kb elements without an apparent transposase gene, as well as three subgroups of about 0.6-kb nDart1-related elements that are all internal deletions of the Dart-like sequences. Both nDart1 and 3.6-kb Dart-like elements were also present in indica varieties 93-11 and Kasalath. nDart1 appears to be the most active mutagen among nDart1-related elements contributing to generating natural variations. A candidate for an autonomous element, aDart, and a possible application of nDart1 for transposon tagging are discussed.

A fine physical map of the rice chromosome 5

A fine physical map of the rice (Oryza sativa spp. Japonica var. Nipponbare) chromosome 5 with bacterial artificial chromosome (BAC) and PI-derived artificial chromosome (PAC) clones was constructed through integration of 280 sequenced BAC/PAC clones and 232 sequence tagged site/expressed sequence tag markers with the use of fingerprinted contig data of the Nipponbare genome. This map consists of five contigs covering 99% of the estimated chromosome size (30.08 Mb). The four physical gaps were estimated at 30 and 20 kb for gaps 1-3 and gap 4, respectively. We have submitted 42.2-Mb sequences with 29.8 Mb of nonoverlapping sequences to public databases. BAC clones corresponding to telomere and centromere regions were confirmed by BAC-fluorescence in situ hybridization (FISH) on a pachytene chromosome. The genetically centromeric region at 54.6 cM was covered by a minimum tiling path spanning 2.1 Mb with no physical gaps. The precise position of the centromere was revealed by using three overlapping BAC/PACs for approximately 150 kb. In addition, FISH results revealed uneven chromatin condensation around the centromeric region at the pachytene stage. This map is of use for positional cloning and further characterization of the rice functional genomics.

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

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

A tiling microarray expression analysis of rice chromosome 4 suggests a chromosome-level regulation of transcription

The complete genome sequence of cultivated rice (Oryza sativa) provides an unprecedented opportunity to understand the biology of this model cereal. An essential and necessary step in this effort is the determination of the coding information and expression patterns of each sequenced chromosome. Here, we report an analysis of the transcriptional activity of rice chromosome 4 using a tiling path microarray based on PCR-generated genomic DNA fragments. Six representative rice organ types were examined using this microarray to catalog the transcribed regions of rice chromosome 4 and to reveal organ- and developmental stage-specific transcription patterns. This analysis provided expression support for 82% of the gene models in the chromosome. Transcriptional activities in 1643 nonannotated regions were also detected. Comparison with cytologically defined chromatin features indicated that in juvenile-stage rice the euchromatic region is more actively transcribed than is the transposon-rich heterochromatic portion of the chromosome. Interestingly, increased transcription of transposon-related gene models in certain heterochromatic regions was observed in mature-stage rice organs and in suspension-cultured cells. These results suggest a close correlation between transcriptional activity and chromosome organization and the developmental regulation of transcription activity at the chromosome level.

High-resolution mapping of theSandZloci ofPhalaris coerulescens

Self incompatibility (SI) in Phalaris coerulescens is gametophytically determined by two unlinked multi allelic loci (S and Z). Neither the S nor Z genes have yet been cloned. As part of a map-based cloning strategy, high-resolution maps of the S and Z regions were generated from distorted segregating populations using RFLP probes from wheat, barley, oat, and Phalaris. The S locus was delimited to 0.26 cM with two boundary markers (Xwg811 and Xpsr168) and cosegregated with Xbm2 and Xbcd762. Xbcd266 was the closest marker linked to Z (0.9 cM). A high level of colinearity in the S and Z regions was found in both self-incompatible and -compatible species. The S locus was localized to the subcentromere region of chromosome 1 and the Z locus to the long arm end of chromosome 2. Several rice BAC clones orthologous to the S and Z locus regions were identified. This opens the possibility of using the rice genome sequence data to generate more closely linked markers and identify SI candidate genes. These results add further support to the conservation of gene order in the S and Z regions of the grass genomes.Key words: Phalaris coerulescens, self-incompatibility, distorted segregation, mapping, map-based cloning, synteny mapping.

Structure and evolution of cereal genomes

The cereal species, of central importance to our diet, began to diverge 50-70 million years ago. For the past few thousand years, these species have undergone largely parallel selection regimes associated with domestication and improvement. The rice genome sequence provides a platform for organizing information about diverse cereals, and together with genetic maps and sequence samples from other cereals is yielding new insights into both the shared and the independent dimensions of cereal evolution. New data and population-based approaches are identifying genes that have been involved in cereal improvement. Reduced-representation sequencing promises to accelerate gene discovery in many large-genome cereals, and to better link the under-explored genomes of 'orphan' cereals with state-of-the-art knowledge.

Construction and characterization of a bacterial artificial chromosome (BAC) library of hexaploid wheat (Triticum aestivum L.) and validation of genome coverage using locus-specific primers

A BAC library of hexaploid wheat was constructed using the spring wheat cultivar Triticum aestivum L. 'Glenlea'. Fresh shoot tissue from 7- to 10-day-old seedlings was used to obtain HMW DNA. The library was constructed using the HindIII site of pIndigoBAC-5 and the BamHI site of pIndigoBAC-5 and pECBAC1. A total of 12 ligations were used to construct the entire library, which contains over 650 000 clones. Ninety-six percent of the clones had inserts. The insert size ranged from 5 to 189 kb with an average of 79 kb. The entire library was gridded onto 24 high-density filters using a 5 x 5 array. A subset of these membranes was hybridized with two intergenic chloroplast probes and the percentage of clones containing chloroplast DNA (cpDNA) was calculated to be 2.2%. The genome coverage was estimated to be 3.1 x haploid genome equivalents, giving a 95.3% probability of identifying a clone corresponding to any wheat DNA sequence. BAC pools were constructed and screened using markers targeting the Glu-B1 locus (1BL), the hardness loci (5AS, 5BS, 5DS), the leaf rust resistance locus Lr1 (5DL), and the major fusarium head blight QTL locus located on 3BS. These markers were either locus-specific amplicons or microsatellites. A total of 49 BAC clones were identified for 14 markers giving an average of 3.5 clones/marker, thereby corroborating the estimated 3.1x genome coverage. An example using the gene encoding the HMW glutenin Bx7 is illustrated.

Development of genome-wide DNA polymorphism database for map-based cloning of rice genes

DNA polymorphism is the basis to develop molecular markers that are widely used in genetic mapping today. A genome-wide rice (Oryza sativa) DNA polymorphism database has been constructed in this work using the genomes of Nipponbare, a cultivar of japonica, and 93-11, a cultivar of indica. This database contains 1,703,176 single nucleotide polymorphisms (SNPs) and 479,406 Insertion/Deletions (InDels), approximately one SNP every 268 bp and one InDel every 953 bp in rice genome. Both SNPs and InDels in the database were experimentally validated. Of 109 randomly selected SNPs, 107 SNPs (98.2%) are accurate. PCR analysis indicated that 90% (97 of 108) of InDels in the database could be used as molecular markers, and 68% to 89% of the 97 InDel markers have polymorphisms between other indica cultivars (Guang-lu-ai 4 and Long-te-pu B) and japonica cultivars (Zhong-hua 11 and 9522). This suggests that this database can be used not only for Nipponbare and 93-11, but also for other japonica and indica cultivars. While validating InDel polymorphisms in the database, a set of InDel markers with each chromosome 3 to 5 marker was developed. These markers are inexpensive and easy to use, and can be used for any combination of japonica and indica cultivars used in this work. This rice DNA polymorphism database will be a valuable resource and important tool for map-based cloning of rice gene, as well as in other various research on rice (http://shenghuan.shnu.edu.cn/ricemarker).

Structural features of the rice chromosome 4 centromere

A complete sequence of a chromosome centromere is necessary for fully understanding centromere function. We reported the sequence structures of the first complete rice chromosome centromere through sequencing a large insert bacterial artificial chromosome clone-based contig, which covered the rice chromosome 4 centromere. Complete sequencing of the 124-kb rice chromosome 4 centromere revealed that it consisted of 18 tracts of 379 tandemly arrayed repeats known as CentO and a total of 19 centromeric retroelements (CRs) but no unique sequences were detected. Four tracts, composed of 65 CentO repeats, were located in the opposite orientation, and 18 CentO tracts were flanked by 19 retroelements. The CRs were classified into four types, and the type I retroelements appeared to be more specific to rice centromeres. The preferential insert of the CRs among CentO repeats indicated that the centromere-specific retroelements may contribute to centromere expansion during evolution. The presence of three intact retrotransposons in the centromere suggests that they may be responsible for functional centromere initiation through a transcription-mediated mechanism.

High-resolution physical mapping in Pennisetum squamulatum reveals extensive chromosomal heteromorphism of the genomic region associated with apomixis

Gametophytic apomixis is asexual reproduction as a consequence of parthenogenetic development of a chromosomally unreduced egg. The trait leads to the production of embryos with a maternal genotype, i.e. progeny are clones of the maternal plant. The application of the trait in agriculture could be a tremendous tool for crop improvement through conventional and nonconventional breeding methods. Unfortunately, there are no major crops that reproduce by apomixis, and interspecific hybridization with wild relatives has not yet resulted in commercially viable germplasm. Pennisetum squamulatum is an aposporous apomict from which the gene(s) for apomixis has been transferred to sexual pearl millet by backcrossing. Twelve molecular markers that are linked with apomixis coexist in a tight linkage block called the apospory-specific genomic region (ASGR), and several of these markers have been shown to be hemizygous in the polyploid genome of P. squamulatum. High resolution genetic mapping of these markers has not been possible because of low recombination in this region of the genome. We now show the physical arrangement of bacterial artificial chromosomes containing apomixis-linked molecular markers by high resolution fluorescence in situ hybridization on pachytene chromosomes. The size of the ASGR, currently defined as the entire hemizygous region that hybridizes with apomixis-linked bacterial artificial chromosomes, was estimated on pachytene and mitotic chromosomes to be approximately 50 Mbp (a quarter of the chromosome). The ASGR includes highly repetitive sequences from an Opie-2-like retrotransposon family that are particularly abundant in this region of the genome.

Structural and functional analysis of rice genome

Rice is an excellent system for plant genomics as it represents a modest size genome of 430 Mb. It feeds more than half the population of the world. Draft sequences of the rice genome, derived by whole-genome shotgun approach at relatively low coverage (4-6 X), were published and the International Rice Genome Sequencing Project (IRGSP) declared high quality (>10 X), genetically anchored, phase 2 level sequence in 2002. In addition, phase 3 level finished sequence of chromosomes 1, 4 and 10 (out of 12 chromosomes of rice) has already been reported by scientists from IRGSP consortium. Various estimates of genes in rice place the number at >50,000. Already, over 28,000 full-length cDNAs have been sequenced, most of which map to genetically anchored genome sequence. Such information is very useful in revealing novel features of macro- and micro-level synteny of rice genome with other cereals. Microarray analysis is unraveling the identity of rice genes expressing in temporal and spatial manner and should help target candidate genes useful for improving traits of agronomic importance. Simultaneously, functional analysis of rice genome has been initiated by marker-based characterization of useful genes and employing functional knock-outs created by mutation or gene tagging. Integration of this enormous information is expected to catalyze tremendous activity on basic and applied aspects of rice genomics.

Transcript abundance of rml1, encoding a putative GT1-like factor in rice, is up-regulated by Magnaporthe grisea and down-regulated by light

We isolated and sequenced both genomic DNA and cDNA clones, which encoded a putative GT1-like protein with 385 amino acids, from cultivated rice (Oryza sativa ssp. indica). This protein shows significant amino acid sequence similarities with trihelix DNA-binding GT-1a/B2F and GT-1 factors that were identified in dicot plants. Northern blotting analysis indicated that the transcript of the rice GT-1 factor in seedling was up-regulated by the rice blast fungus Magnaporthe grisea, down-regulated by various continuous light conditions and expressed rhythmically in light/dark cycles. This GT1-like factor gene was therefore designated as rml1 (rice gene regulated by M. grisea and light). The putative RML1 protein, encoded by this single copy gene, is thus identified as a new member of the plant-specific GT family of transcription factors in rice.

Sequence variations of simple sequence repeats on chromosome-4 in two subspecies of the Asian cultivated rice

Computational screening of the chromosome-4 sequence of the rice cultivar Nipponbare ( Oryza sativa L. japonica) revealed 1,844 tandem simple sequence repeats (SSRs) or microsatellites with SSR motifs >/=20 bp and repeated unit length of 1-6 base pairs. Thus SSRs occur once in every 18.8 kb, on the average, on the chromosome with one SSR per 23.8 kb and 16 kb on the short and long arms, respectively. No SSR was detected in the core region of the centromere. Poly(AT)(n) repeats represented the most abundant and length polymorphic class of SSRs on the chromosome, but it did not occur in the exons. GC-rich trinucleotide repeats were most abundant in the coding regions, representing 71.69% of the SSRs identified in the exons. Two hundred and twenty four SSRs were associated with the repetitive DNA sequences, most of them were poly(AT)(n) tracts. Sequence variations of SSRs between two cultivars, representing the two subspecies of the Asian cultivated rice indica and japonica, were identified, revealing that divergence and convergence of the two subspecies could be traced by the analysis of SSRs. These results provide a great opportunity for SSR-based marker development and comparative genome analysis of the two subspecies of the Asian cultivated rice.

Establishing an efficient Ac/Ds tagging system in rice: large-scale analysis of Ds flanking sequences

A two-element Activator/Dissociation (Ac/Ds) gene trap system was successfully established in rice (Oryza sativa ssp. japonica cv. Nipponbare) to generate a collection of stable, unlinked and single-copy Ds transposants. The germinal transposition frequency of Ds was estimated as an average of 51% by analyzing 4413 families. Study of Ds transposition pattern in siblings revealed that 79% had at least two different insertions, suggesting late transposition during rice development. Analysis of 2057 Ds flanking sequences showed that 88% of them were unique, whereas the rest within T-DNA. The insertions were distributed randomly throughout the genome; however, there was a bias toward chromosomes 4 and 7, which had two times as many insertions as that expected. A hot spot for Ds insertions was identified on chromosome 7 within a 40-kbp region. One-third of Ds flanking sequences was homologous to either proteins or rice expressed sequence tags (ESTs), confirming a preference for Ds transposition into coding regions. Analysis of 200 Ds lines on chromosome 1 revealed that 72% insertions were found in genic region. Anchoring of more than 800 insertions to yeast artificial chromosome (YAC)-based EST map showed that Ds transposes preferentially into regions rich in expressed sequences. High germinal transposition frequency and independent transpositions among siblings show that the efficiency of this system is suitable for large-scale transposon mutagenesis in rice.

Physical maps and recombination frequency of six rice chromosomes

We constructed physical maps of rice chromosomes 1, 2, and 6-9 with P1-derived artificial chromosome (PAC) and bacterial artificial chromosome (BAC) clones. These maps, with only 20 gaps, cover more than 97% of the predicted length of the six chromosomes. We submitted a total of 193 Mbp of non-overlapping sequences to public databases. We analyzed the DNA sequences of 1316 genetic markers and six centromere-specific repeats to facilitate characterization of chromosomal recombination frequency and of the genomic composition and structure of the centromeric regions. We found marked changes in the relative recombination rate along the length of each chromosome. Chromosomal recombination at the centromere core and surrounding regions on the six chromosomes was completely suppressed. These regions have a total physical length of about 23 Mbp, corresponding to 11.4% of the entire size of the six chromosomes. Chromosome 6 has the longest quiescent region, with about 5.6 Mbp, followed by chromosome 8, with quiescent region about half this size. Repetitive sequences accounted for at least 40% of the total genomic sequence on the partly sequenced centromeric region of chromosome 1. Rice CentO satellite DNA is arrayed in clusters and is closely associated with the presence of Centromeric Retrotransposon of Rice (CRR)- and RIce RetroElement 7 (RIRE7)-like retroelement sequences. We also detected relatively small coldspot regions outside the centromeric region; their repetitive content and gene density were similar to those of regions with normal recombination rates. Sequence analysis of these regions suggests that either the amount or the organization patterns of repetitive sequences may play a role in the inactivation of recombination.

Fine genetic mapping of a gene required for Rice yellow mottle virus cell-to-cell movement

The very high resistance to Rice yellow mottle virus observed in the two rice varieties Gigante ( Oryza sativa) and Tog 5681 ( O. glaberrima) is monogenic and recessive. Bulked segregant analysis was carried out to identify AFLP markers linked to the resistance gene. Mapping of PCR-specific markers, CAPS and microsatellite markers on 429 individuals of an IR64 x Gigante F(2) population pinpointed this resistance gene on the long arm of chromosome 4 in a 3.7-cM interval spanned by PCR markers. These markers also flanked the resistance gene of the O. glaberrima accession Tog 5681 and confirmed previous allelism tests. The rarity of this recessive natural resistance was in line with a resistance mechanism model based on point mutations of a host component required for cell-to-cell movement of the virus. Preliminary data on the genetic divergence between the two cultivated rice species in the vicinity of the resistance locus suggested that two different resistance alleles are present in Gigante and Tog 5681. A large set of recombinants is now available to envisage physical mapping and cloning of the gene.

Genome-wide intraspecific DNA-sequence variations in rice

Genome-wide comparative analysis of the DNA sequences of two major cultivated rice subspecies, Oryza sativa L. ssp indica and Oryza sativa L. ssp japonica, have revealed their extensive microcolinearity in gene order and content. However, deviations from colinearity are frequent owing to insertions or deletions. Intraspecific sequence polymorphisms commonly occur in both coding and non-coding regions. These variations often affect gene structures and may contribute to intraspecific phenotypic adaptations.

Towards an accurate sequence of the rice genome

Several more- or less-elaborated rice genome sequences have been produced recently using different strategies. It has become possible to compare them and to unravel the major features of the rice genome in terms of nucleotide composition, repeats, gene content and variability. It has also become possible to compare the rice and Arabidopsis genomes and to evaluate rice as a model genome.

Sequence and analysis of rice chromosome 4

Rice is the principal food for over half of the population of the world. With its genome size of 430 megabase pairs (Mb), the cultivated rice species Oryza sativa is a model plant for genome research. Here we report the sequence analysis of chromosome 4 of O. sativa, one of the first two rice chromosomes to be sequenced completely. The finished sequence spans 34.6 Mb and represents 97.3% of the chromosome. In addition, we report the longest known sequence for a plant centromere, a completely sequenced contig of 1.16 Mb corresponding to the centromeric region of chromosome 4. We predict 4,658 protein coding genes and 70 transfer RNA genes. A total of 1,681 predicted genes match available unique rice expressed sequence tags. Transposable elements have a pronounced bias towards the euchromatic regions, indicating a close correlation of their distributions to genes along the chromosome. Comparative genome analysis between cultivated rice subspecies shows that there is an overall syntenic relationship between the chromosomes and divergence at the level of single-nucleotide polymorphisms and insertions and deletions. By contrast, there is little conservation in gene order between rice and Arabidopsis.