A comparison of rice chloroplast genomes

Using high quality sequence reads extracted from our whole genome shotgun repository, we assembled two chloroplast genome sequences from two rice (Oryza sativa) varieties, one from 93-11 (a typical indica variety) and the other from PA64S (an indica-like variety with maternal origin of japonica), which are both parental varieties of the super-hybrid rice, LYP9. Based on the patterns of high sequence coverage, we partitioned chloroplast sequence variations into two classes, intravarietal and intersubspecific polymorphisms. Intravarietal polymorphisms refer to variations within 93-11 or PA64S. Intersubspecific polymorphisms were identified by comparing the major genotypes of the two subspecies represented by 93-11 and PA64S, respectively. Some of the minor genotypes occurring as intravarietal polymorphisms in one variety existed as major genotypes in the other subspecific variety, thus giving rise to intersubspecific polymorphisms. In our study, we found that the intersubspecific variations of 93-11 (indica) and PA64S (japonica) chloroplast genomes consisted of 72 single nucleotide polymorphisms and 27 insertions or deletions. The intersubspecific polymorphism rates between 93-11 and PA64S were 0.05% for single nucleotide polymorphisms and 0.02% for insertions or deletions, nearly 8 and 10 times lower than their respective nuclear genomes. Based on the total number of nucleotide substitutions between the two chloroplast genomes, we dated the divergence of indica and japonica chloroplast genomes as occurring approximately 86,000 to 200,000 years ago.     

The Organelle Genomes of Hassawi Rice (Oryza sativa L.) and Its Hybrid in Saudi Arabia: Genome Variation, Rearrangement, and Origins

Hassawi rice (Oryza sativa L.) is a landrace adapted to the climate of Saudi Arabia, characterized by its strong resistance to soil salinity and drought. Using high quality sequencing reads extracted from raw data of a whole genome sequencing project, we assembled both chloroplast (cp) and mitochondrial (mt) genomes of the wild-type Hassawi rice (Hassawi-1) and its dwarf hybrid (Hassawi-2). We discovered 16 InDels (insertions and deletions) but no SNP (single nucleotide polymorphism) is present between the two Hassawi cp genomes. We identified 48 InDels and 26 SNPs in the two Hassawi mt genomes and a new type of sequence variation, termed reverse complementary variation (RCV) in the rice cp genomes. There are two and four RCVs identified in Hassawi-1 when compared to 93-11 (indica) and Nipponbare (japonica), respectively. Microsatellite sequence analysis showed there are more SSRs in the genic regions of both cp and mt genomes in the Hassawi rice than in the other rice varieties. There are also large repeats in the Hassawi mt genomes, with the longest length of 96,168 bp and 96,165 bp in Hassawi-1 and Hassawi-2, respectively. We believe that frequent DNA rearrangement in the Hassawi mt and cp genomes indicate ongoing dynamic processes to reach genetic stability under strong environmental pressures. Based on sequence variation analysis and the breeding history, we suggest that both Hassawi-1 and Hassawi-2 originated from the Indonesian variety Peta since genetic diversity between the two Hassawi cultivars is very low albeit an unknown historic origin of the wild-type Hassawi rice.     

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).     

Chloroplast genome variation in upland and lowland switchgrass

Switchgrass (Panicum virgatum L.) exists at multiple ploidies and two phenotypically distinct ecotypes. To facilitate interploidal comparisons and to understand the extent of sequence variation within existing breeding pools, two complete switchgrass chloroplast genomes were sequenced from individuals representative of the upland and lowland ecotypes. The results demonstrated a very high degree of conservation in gene content and order with other sequenced plastid genomes. The lowland ecotype reference sequence (Kanlow Lin1) was 139,677 base pairs while the upland sequence (Summer Lin2) was 139,619 base pairs. Alignments between the lowland reference sequence and short-read sequence data from existing sequence datasets identified as either upland or lowland confirmed known polymorphisms and indicated the presence of other differences. Insertions and deletions principally occurred near stretches of homopolymer simple sequence repeats in intergenic regions while most Single Nucleotide Polymorphisms (SNPs) occurred in intergenic regions and introns within the single copy portions of the genome. The polymorphism rate between upland and lowland switchgrass ecotypes was found to be similar to rates reported between chloroplast genomes of indica and japonica subspecies of rice which were believed to have diverged 0.2-0.4 million years ago.     

A genome-wide analysis of single nucleotide polymorphism diversity in the world's major cereal crops

Over 3.5 million expressed sequence tags from the major cereal taxa were used to electronically mine over 176 000 putative single nucleotide polymorphisms (SNPs). The density, distribution and degree of linkage between these SNPs were compared among the different taxa. The frequency of sequence polymorphism was lowest in diploid taxa (rice, barley and sorghum), intermediate in tetraploid maize and highest in allohexaploid wheat and octoploid sugarcane. SNPs were further categorized as either intravarietal (differences between gene family members and homoeologues) or varietal (differences between two varieties), and as either co-segregating or non-co-segregating with neighbouring polymorphisms. Varietal co-segregating SNPs represent the best candidates for molecular markers as they show variation between varieties and have a high probability of being validated, as sequencing errors are unlikely to co-segregate with one another. This elite class of SNPs was most abundant in barley and least abundant in wheat and rice. Despite the large number of observed sequence polymorphisms in allohexaploid wheat, only a fraction of those available are likely to make good molecular markers. In addition, we found that rice SNPs up to 10 kb apart were in linkage disequilibrium (LD), but that high levels of LD attributable to population structure confounded the tracking of LD over greater distances.     

Genome-wide searching of single-nucleotide polymorphisms among eight distantly and closely related rice cultivars (Oryza sativa L.) and a wild accession (Oryza rufipogon Griff.).

We searched the genomes of eight rice cultivars (Oryza sativa L. ssp. japonica and ssp. indica) and a wild rice accession (Oryza rufipogon Griffith) for nucleotide polymorphisms, and identified 7805 polymorphic loci, including single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels), in predicted intergenic regions. Polymorphisms are useful as DNA markers for genetic analysis or positional cloning with segregating populations of crosses. Pairwise comparison between cultivars and a neighbor-joining tree calculated from SNPs agreed very well with relationships between rice strains predicted from pedigree data or calculated with other DNA markers such as p-SINE1 and simple sequence repeats (SSRs), suggesting that whole-genome SNP information can be used for analysis of evolutionary relationships. Using multiple SNPs to identify alleles, we drew a map to illustrate the alleles shared among the eight cultivars and the accession. The map revealed that most of the genome is mono- or di-allelic among japonica cultivars, whereas alleles well conserved among modern japonica paddy rice cultivars were often shared with indica cultivars or wild rice, suggesting that the genome structure of modern cultivars is composed of chromosomal segments from various genetic backgrounds. Use of allele-sharing analysis and association analysis were also tested and are discussed.     

Direct Chloroplast Sequencing: Comparison of Sequencing Platforms and Analysis Tools for Whole Chloroplast Barcoding

Direct sequencing of total plant DNA using next generation sequencing technologies generates a whole chloroplast genome sequence that has the potential to provide a barcode for use in plant and food identification. Advances in DNA sequencing platforms may make this an attractive approach for routine plant identification. The HiSeq (Illumina) and Ion Torrent (Life Technology) sequencing platforms were used to sequence total DNA from rice to identify polymorphisms in the whole chloroplast genome sequence of a wild rice plant relative to cultivated rice (cv. Nipponbare). Consensus chloroplast sequences were produced by mapping sequence reads to the reference rice chloroplast genome or by de novo assembly and mapping of the resulting contigs to the reference sequence. A total of 122 polymorphisms (SNPs and indels) between the wild and cultivated rice chloroplasts were predicted by these different sequencing and analysis methods. Of these, a total of 102 polymorphisms including 90 SNPs were predicted by both platforms. Indels were more variable with different sequencing methods, with almost all discrepancies found in homopolymers. The Ion Torrent platform gave no apparent false SNP but was less reliable for indels. The methods should be suitable for routine barcoding using appropriate combinations of sequencing platform and data analysis.     

Microsatellite diversity within Oryza sativa with emphasis on indica-japonica divergence

The molecular evolution of cultivated rice Oryza sativa L. has long been a subject of rice evolutionists. To investigate genetic diversity within and differentiation between the indica and japonica subspecies, 22 accessions of indica and 35 of japonica rice were examined by five microsatellite loci from each chromosome totalling 60 loci. Mean gene diversity value in the indica rice (H=0.678) was 1.18 times larger than in the japonica rice (H=0.574). Taking the sampling effect into consideration, average allele number in the indica rice was 1.40 times higher than that in the japonica rice (14.6 vs 10.4 per variety). Chromosome-based comparisons revealed that nine chromosomes (1, 2, 3, 4, 5, 8, 9, 10 and 11) harboured higher levels of genetic diversity within the indica rice than the japonica rice. An overall estimate of F(ST) was 0.084-0.158, indicating that the differentiation is moderate and 8.4-15.8% of the total genetic variation resided between the indica and japonica groups. Our chromosome-based comparisons further suggested that the extent of the indica-japonica differentiation varied substantially, ranging from 7.62% in chromosome 3 to 28.72% in chromosome 1. Cluster analyses found that most varieties formed merely two clusters for the indica and japonica varieties, in which two japonica varieties and five indica varieties were included in the counterpart clusters, respectively. The 12 chromosome-based trees further showed that 57 rice varieties cannot be clearly clustered together into either the indica or japonica groups, but displayed relatively different clustering patterns. The results suggest that the process of indica japonica differentiation may have proceeded through an extensive contribution by the alleles of the majority in the rice genome.     

Mutational decay and age of chloroplast and mitochondrial genomes transferred recently to angiosperm nuclear chromosomes

Transfers of organelle DNA to the nucleus established several thousand functional genes in eukaryotic chromosomes over evolutionary time. Recent transfers have also contributed nonfunctional plastid (pt)- and mitochondrion (mt)-derived DNA (termed nupts and numts, respectively) to plant nuclear genomes. The two largest transferred organelle genome copies are 131-kb nuptDNA in rice (Oryza sativa) and 262-kb numtDNA in Arabidopsis (Arabidopsis thaliana). These transferred copies were compared in detail with their bona fide organelle counterparts, to which they are 99.77% and 99.91% identical, respectively. No evidence for purifying selection was found in either nuclear integrant, indicating that they are nonfunctional. Mutations attributable to 5-methylcytosine hypermutation have occurred at a 6- to 10-fold higher rate than other point mutations in Arabidopsis numtDNA and rice nuptDNA, respectively, revealing this as a major mechanism of mutational decay for these transferred organelle sequences. Short indels occurred preferentially within homopolymeric stretches but were less frequent than point mutations. The 131-kb nuptDNA is absent in the O. sativa subsp. indica or Oryza rufipogon nuclear genome, suggesting that it was transferred within the O. sativa subsp. japonica lineage and, as revealed by sequence comparisons, after its divergence from the indica chloroplast lineage. The time of the transfer for the rice nupt was estimated as 148,000 (74,000--296,000) years ago and that for the Arabidopsis numtDNA as 88,000 (44,000--176,000) years ago. The results reveal transfer and integration of entire organelle genomes into the nucleus as an ongoing evolutionary process and uncover mutational mechanisms affecting organelle genomes recently transferred into a new mutational environment.     

Genome analysis of the domestic dog (Korean Jindo) by massively parallel sequencing

Although pioneering sequencing projects have shed light on the boxer and poodle genomes, a number of challenges need to be met before the sequencing and annotation of the dog genome can be considered complete. Here, we present the DNA sequence of the Jindo dog genome, sequenced to 45-fold average coverage using Illumina massively parallel sequencing technology. A comparison of the sequence to the reference boxer genome led to the identification of 4 675 437 single nucleotide polymorphisms (SNPs, including 3 346 058 novel SNPs), 71 642 indels and 8131 structural variations. Of these, 339 non-synonymous SNPs and 3 indels are located within coding sequences (CDS). In particular, 3 non-synonymous SNPs and a 26-bp deletion occur in the TCOF1 locus, implying that the difference observed in cranial facial morphology between Jindo and boxer dogs might be influenced by those variations. Through the annotation of the Jindo olfactory receptor gene family, we found 2 unique olfactory receptor genes and 236 olfactory receptor genes harbouring non-synonymous homozygous SNPs that are likely to affect smelling capability. In addition, we determined the DNA sequence of the Jindo dog mitochondrial genome and identified Jindo dog-specific mtDNA genotypes. This Jindo genome data upgrade our understanding of dog genomic architecture and will be a very valuable resource for investigating not only dog genetics and genomics but also human and dog disease genetics and comparative genomics.     

The contribution of intersubspecific hybridization to?the breeding of?super-high-yielding japonica rice in northeast China

Hybridization between indica and japonica rice combined with utilization of ideal plant type has led to the development of high-yielding japonica rice in northern China. However, the contribution at the genomic level of intersubspecific hybridization to the increased yield of northern Chinese japonica rice is uncertain. In this study, we analyzed the genomic pedigree of descendants of hybridization between indica and japonica rice grown in northeastern China between 1963 and 2008. Simple sequence repeat markers indicated that since 1990 the genetic diversity among northern japonica cultivars was enriched. Genome-wide analysis with subspecies-specific indel and intron length polymorphism markers showed indica-allele frequencies were significantly increased in cultivars bred after 1990, and were significantly positively correlated with spikelet number per panicle and significantly negatively correlated with panicle number per plant. Among eight genes controlling agronomic traits, GN1a and GS3 were partially fixed in the genome of northern japonica cultivars. In contrast, Waxy and qSH1 were eliminated, whereas DEP1 and qSW5 were retained. Indica germplasm is an important contributor to the increased yield of northern japonica rice. Breeding for high yield and grain quality in combination is a complicated process and difficult to achieve when relying on only one or several functional genes, thus the selection expertise of the breeder remains critical.     

Characterization of the major fragance gene from an aromatic japonica rice and analysis of its diversity in Asian cultivated rice

In Asian cultivated rice (Oryza sativa L.), aroma is one of the most valuable traits in grain quality and 2-ACP is the main volatile compound contributing to the characteristic popcorn-like odour of aromatic rices. Although the major locus for grain fragrance (frg gene) has been described recently in Basmati rice, this gene has not been characterised in true japonica varieties and molecular information available on the genetic diversity and evolutionary origin of this gene among the different varieties is still limited. Here we report on characterisation of the frg gene in the Azucena variety, one of the few aromatic japonica cultivars. We used a RIL population from a cross between Azucena and IR64, a non-aromatic indica, the reference genomic sequence of Nipponbare (japonica) and 93-11 (indica) as well as an Azucena BAC library, to identify the major fragance gene in Azucena. We thus identified a betaine aldehyde dehydrogenase gene, badh2, as the candidate locus responsible for aroma, which presented exactly the same mutation as that identified in Basmati and Jasmine-like rices. Comparative genomic analyses showed very high sequence conservation between Azucena and Nipponbare BADH2, and a MITE was identified in the promotor region of the BADH2 allele in 93-11. The badh2 mutation and MITE were surveyed in a representative rice collection, including traditional aromatic and non-aromatic rice varieties, and strongly suggested a monophylogenetic origin of this badh2 mutation in Asian cultivated rices. Altogether these new data are discussed here in the light of current hypotheses on the origin of rice genetic diversity.     

Genome Sequence of Herpes Simplex Virus 1 Strain McKrae

The herpes simplex virus 1 (HSV-1) strain McKrae is highly virulent compared to other wild-type strains of HSV-1. To help us better understand the genetic determinants that lead to differences in the pathogenicity of McKrae and other HSV-1 strains, we sequenced its genome. Comparing the sequence of McKrae's genome to that of strain 17 revealed that the genomes differ by at least 752 single nucleotide polymorphisms (SNPs) and 86 insertion/deletion events (indels). Although the majority of these polymorphisms reside in noncoding regions, 241 SNPs and 10 indels alter the protein-coding sequences of 58 open reading frames. Some of these variations are expected to contribute to the pathogenic phenotype of McKrae.     

The extent of linkage disequilibrium in rice (Oryza sativa L.)

Despite its status as one of the world's major crops, linkage disequilibrium (LD) patterns have not been systematically characterized across the genome of Asian rice (Oryza sativa). Such information is critical to fully exploit the genome sequence for mapping complex traits using association techniques. Here we characterize LD in five 500-kb regions of the rice genome in three major cultivated rice varieties (indica, tropical japonica, and temperate japonica) and in the wild ancestor of Asian rice, Oryza rufipogon. Using unlinked SNPs to determine the amount of background linkage disequilibrium in each population, we find that the extent of LD is greatest in temperate japonica (probably >500 kb), followed by tropical japonica (approximately 150 kb) and indica (approximately 75 kb). LD extends over a shorter distance in O. rufipogon (<40 kb) than in any of the O. sativa groups assayed here. The differences in the extent of LD among these groups are consistent with differences in outcrossing and recombination rate estimates. As well as heterogeneity between groups, our results suggest variation in LD patterns among genomic regions. We demonstrate the feasibility of genomewide association mapping in cultivated Asian rice using a modest number of SNPs.     

Genome sequencing of a 239-kb region of rice chromosome 10L reveals a high frequency of gene duplication and a large chloroplast DNA insertion

In this study we describe a 239-kb region on the long arm of rice chromosome 10 that contains a high density (71%) of locally duplicated genes, including 24 copies of a glutathione S-transferase gene. Intriguingly, embedded within this cluster is a large insertion (approximately 33 kb) of rice (Oryza sativa) chloroplast DNA that is derived from two separate regions of the chloroplast genome. We used DNA fiber-based fluorescence in situ hybridization (fiber-FISH) analyses of O. sativa spp. japonica nuclei to confirm that the insertion of organellar DNA was not a cloning artifact. The sequence of the chloroplast insertion is nearly identical (99.7% identity) to the corresponding regions in the published rice chloroplast genome sequence, suggesting that the transfer event occurred recently. PCR amplification and sequence analysis in two subspecies of rice, O. sativa spp. japonica and spp. indica, indicates that the transfer event predated the divergence of these two subspecies. The chloroplast insertion is flanked by a 2.1-kb perfect direct repeat that is unique to this location in the rice genome.     

Genetic diversity analysis of traditional and improved Indonesian rice (Oryza sativa L.) germplasm using microsatellite markers

The archipelago of Indonesia has a long history of rice production across a broad range of rice-growing environments resulting in a diverse array of local Indonesian rice varieties. Although some have been incorporated into modern breeding programs, the vast majority of these landraces remain untapped. To better understand this rich source of genetic diversity we have characterized 330 rice accessions, including 246 Indonesian landraces and 63 Indonesian improved cultivars, using 30 fluorescently-labeled microsatellite markers. The landraces were selected across 21 provinces and include representatives of the classical subpopulations of cere, bulu, and gundil rices. A total of 394 alleles were detected at the 30 simple sequence repeat loci, with an average number of 13 alleles per locus across all accessions, and an average polymorphism information content value of 0.66. Genetic diversity analysis characterized the Indonesian landraces as 68% indica and 32% tropical japonica, with an indica gene diversity of 0.53 and a tropical japonica gene diversity of 0.56, and a Fst of 0.38 between the two groups. All of the improved varieties sampled were indica, and had an average gene diversity of 0.46. A set of high quality Indonesian varieties, including Rojolele, formed a separate cluster within the tropical japonicas. This germplasm presents a valuable source of diversity for future breeding and association mapping efforts.