Phylogeography of Asian wild rice, Oryza rufipogon: a genome-wide view

Asian wild rice (Oryza rufipogon) that ranges widely across the eastern and southern part of Asia is recognized as the direct ancestor of cultivated Asian rice (O. sativa). Studies of the geographic structure of O. rufipogon, based on chloroplast and low‐copy nuclear markers, reveal a possible phylogeographic signal of subdivision in O. rufipogon. However, this signal of geographic differentiation is not consistently observed among different markers and studies, with often conflicting results. To more precisely characterize the phylogeography of O. rufipogon populations, a genome‐wide survey of unlinked markers, intensively sampled from across the entire range of O. rufipogon is critical. In this study, we surveyed sequence variation at 42 genome‐wide sequence tagged sites (STS) in 108 O. rufipogon accessions from throughout the native range of the species. Using Bayesian clustering, principal component analysis and amova , we conclude that there are two genetically distinct O. rufipogon groups, Ruf‐I and Ruf‐II. The two groups exhibit a clinal variation pattern generally from north‐east to south‐west. Different from many earlier studies, Ruf‐I, which is found mainly in China and the Indochinese Peninsula, shows genetic similarity with one major cultivated rice variety, O. satvia indica, whereas Ruf‐II, mainly from South Asia and the Indochinese Peninsula, is not found to be closely related to cultivated rice varieties. The other major cultivated rice variety, O. sativa japonica, is not found to be similar to either O. rufipogon groups. Our results support the hypothesis of a single origin of the domesticated O. sativa in China. The possible role of palaeoclimate, introgression and migration–drift balance in creating this clinal variation pattern is also discussed.     

A chromosome-level genome assembly of the beet armyworm Spodoptera exigua

BackgroundThe beet armyworm Spodoptera exigua is a polyphagous caterpillar that causes serious damage to many species of crops and vegetables. To gain insight into how this polyphagous insect differs from less harmful oligophagous species, we generated a chromosome-level assembly and compared it to closely related species with the same or different feeding habits.ResultsBased on Illumina and Pacific Biosciences data and Hi-C technology, 425.6 Mb of genome sequences were anchored and oriented into 31 linkage groups, with an N50 length of 14.8 Mb. A total of 24,649 gene models were predicted, of which 97.4% were identified in the genome assembly. Chemosensory genes are vital for locating food: of the four main families, odorant-binding proteins, chemosensory proteins and olfactory receptors showed little difference, whereas gustatory receptors are greatly expanded in S. exigua. Examination of other polyphagous insects confirmed this difference from oligophagous congeners and further identified the bitter receptor subfamily as being particularly affected.ConclusionOur high-quality genome sequence for beet armyworm identified a key expansion of the bitter gustatory receptor subfamily in this and other pests that differs crucially from more benign relatives and offers insight into the biology and possible future means of control for these economically important insects.     

Nucleotide polymorphism in the Adh2 region of the wild rice Oryza rufipogon

DNA variation in the alcohol dehydrogenase (Adh2) region of the wild rice Oryza rufipogon and its related species was analyzed to clarify maintenance mechanisms of the DNA variation in these species. A dimorphic pattern was detected in the Adh2 region of O. rufipogon. The silent nucleotide diversity (π) in the Adh2 region in O. rufipogon was 0.011, which was higher than that of the Adh1 region in O. rufipogon. Especially, a high nucleotide diversity was detected at synonymous sites of the catalytic domain 1. Average nucleotide diversity at silent sites within each of the dimorphic sequence types of the Adh2 region was similar to that in the Adh1 region, indicating that the high level of silent polymorphism in the Adh2 region was caused by the difference between the dimorphic sequence types. On the other hand, the level of replacement polymorphism in the Adh2 region was as low as that in the Adh1 region. The neutrality test of Fu and Li indicated significantly negative deviation from the neutral mutation model for the replacement sites of the Adh2 region. This result suggests purifying selection on the replacement sites of the Adh2 region, as detected for the Adh1 region. Significant linkage disequilibria (16.4% of the tests) were detected between the Adh1 and Adh2 regions. Even when nonrandom association was tested for the strains belonging to one of the divergent sequence types of the Adh2 region, significant interlocus linkage disequilibria were detected. The close physical distance and/or epistasis between the two Adh regions could be invoked to explain these nonrandom associations.     

The complete nucleotide sequence of wild rice (Oryza nivara) chloroplast genome: first genome wide comparative sequence analysis of wild and cultivated rice

We determined the complete nucleotide sequence of the chloroplast genome of wild rice, Oryza nivara and compared it with the corresponding published sequence of relative cultivated rice, Oryza sativa. The genome was 134,494 bp long with a large single-copy region of 80,544 bp, a small single-copy region of 12,346 bp and two inverted repeats of 20,802 bp each. The overall A+T content was 61.0%. The O. nivara chloroplast genome encoded identical functional genes to O. sativa in the same order along the genome. On the other hand, detailed analysis revealed 57 insertion, 61 deletion and 159 base substitution events in the entire chloroplast genome of O. nivara. Among substitutions, transversions were much higher than transitions with the former even more frequent than the latter in the coding region. Most of the insertions/deletions were single-base but a few large length mutations were also detected. The frequency of insertion/deletion events was more in the coding region within inverted repeats. In contrast, a very few substitution events were identified in the coding region. Polymorphism was observed among rice cultivars at loci of large insertion/deletion events. This is the first report describing comparative and genome wide chloroplast analysis between a wild and cultivated crop.     

Plant regeneration from protoplasts of wild rice (Oryza rufipogon Griff.)

Embryogenic callus initiated from basal segments of micropropagated shoots of Oryza rufipogon were used to initiate cell suspension cultures. After approximately 3 months these cultures were capable of yielding large numbers of protoplasts which underwent sustained division in agarose-solidified medium at a frequency comparable to that observed with Japonica rice protoplasts in previous studies. O. rufipogon plants were reproducibly regenerated from the protoplast-derived callus and are currently being grown to maturity. This is the first report of plant regeneration from protoplasts of a wild species of Oryza.Abbreviations BAP 6-benzyalamino purine - 2,4-D 2,4-dichlorophenoxyacetic acid - MES 2[N-morpholino] ethanesulphonic acid - NAA -naphthalene acetic acid - PAR photosynthetically active radiation - SCV settled-cell volume     

Genetic variation in the chloroplast genome suggests multiple domestication of cultivated Asian rice (Oryza sativa L.).

Two hundred and seventy-five accessions of cultivated Asian rice and 44 accessions of AA genome Oryza species were classified into 8 chloroplast (cp) genome types (A-H) based on insertion-deletion events at 3 regions (8K, 57K, and 76K) of the cp genome. The ancestral cp genome type was determined according to the frequency of occurrence in Oryza species and the likely evolution of the variable 57K region of the cp genome. When 2 nucleotide substitutions (AA or TT) were taken into account, these 8 cp types were subdivided into 11 cp types. Most indica cultivars had 1 of 3 cp genome types that were also identified in the wild relatives of rice, O. nivara and O. rufipogon, suggesting that the 3 indica cp types had evolved from distinct gene pools of the O. rufipogon - O. nivara complex. The majority of japonica cultivars had 1 of 3 different cp genome types. One of these 3 was identified in O. rufipogon, suggesting that at least 1 japonica type is derived from O. rufipogon with the same cp genome type. These results provide evidence to support a polyphyletic origin of cultivated Asian rice from at least 4 principal lineages in the O. rufipogon - O. nivara complex.     

De novo chromosome-level assembly of the Centella asiatica genome

Centella asiatica is a herbaceous, perennial species indigenous to India and Southeast Asia. C. asiatica possesses several medicinal properties: anti-aging, anti-inflammatory, wound healing and memory enhancing. The lack of available genomics resources significantly impedes the improvement of C. asiatica varieties through molecular breeding. Here, we combined the 10× Genomics linked-read technology and the long-range HiC technique to obtain the genome assembly. The final assembly contained nine pseudomolecules, corresponding to the haploid chromosome number in C. asiatica. These nine chromosomes covered 402,536,584 bases or 93.6% of the 430-Mb assembly. Comparative genomics analyses based on single-copy orthologous genes showed that C. asiatica and the common ancestor of Coriandrum sativum (coriander) and Daucus carota (carrot) diverged about 48 million years ago. This assembly provides a valuable reference genome for future molecular studies, varietal development through marker-assisted breeding and comparative genomics studies in C. asiatica.     

A chromosome-level genome assembly of the striped catfish (Pangasianodon hypophthalmus)

Striped catfish (Pangasianodon hypophthalmus), belonging to the Pangasiidae family, has become an economically important fish with wide cultivation in Southeast Asia. Owing to the high-fat trait, it is always considered as an oily fish. In our present study, a high-quality genome assembly of the striped catfish was generated by integration of Illumina short reads, Nanopore long reads and Hi-C data. A 731.7-Mb genome assembly was finally obtained, with a contig N50 of 3.5 Mb, a scaffold N50 of 29.5 Mb, and anchoring of 98.46% of the assembly onto 30 pseudochromosomes. The genome contained 36.9% repeat sequences, and a total 18,895 protein-coding genes were predicted. Interestingly, we identified a tandem triplication of fatty acid binding protein 1 gene (fabp1; thereby named as fabp1-1, fabp1-2 and fabp1-3 respectively), which may be related to the high fat content in striped catfish. Meanwhile, the FABP1-2 and -3 isoforms differed from FABP1-1 by several missense mutations including R126T, which may affect the fatty acid binding properties. In summary, we report a high-quality chromosome-level genome assembly of the striped catfish, which provides a valuable genetic resource for biomedical studies on the high-fat trait, and lays a solid foundation for practical aquaculture and molecular breeding of this international teleost species.     

A chromosome-level genome assembly of the potato grouper (Epinephelus tukula)

The potato grouper, Epinephelus tukula, is one of the largest coral reef teleost, and it is an important germplasm resource for selection and cross breeding. Here we report a potato grouper genome assembly generated using PacBio long-read sequencing, Illumina sequencing and high-throughput chromatin conformation capture (Hi-C) technology. The genome size was 1.13 Gb, with a total of 508 contigs anchored into 24 chromosomes. The scaffold N50 was 42.65 Mb. For the genome models, our assembled genome contained 98.11% complete BUSCO with the vertebrata_odb9 database. One more copies of Gh and Hsp90b1 were identified in the E. tukula genome, which might contribute to its fast growth and high resistance to stress. In addition, 435 putative antimicrobial peptide (AMP) genes were identified in the potato grouper. This study provides a good reference for whole genome selective breeding of the potato grouper and for future development of novel marine drugs.     

Molecular evolution of the Pi-ta gene resistant to rice blast in wild rice (Oryza rufipogon)

Rice blast disease resistance to the fungal pathogen Magnaporthe grisea is triggered by a physical interaction between the protein products of the host R (resistance) gene, Pi-ta, and the pathogen Avr (avirulence) gene, AVR-pita. The genotype variation and resistant/susceptible phenotype at the Pi-ta locus of wild rice (Oryza rufipogon), the ancestor of cultivated rice (O. sativa), was surveyed in 36 locations worldwide to study the molecular evolution and functional adaptation of the Pi-ta gene. The low nucleotide polymorphism of the Pi-ta gene of O. rufipogon was similar to that of O. sativa, but greatly differed from what has been reported for other O. rufipogon genes. The haplotypes can be subdivided into two divergent haplogroups named H1 and H2. H1 is derived from H2, with nearly no variation and at a low frequency. H2 is common and is the ancestral form. The leucine-rich repeat (LRR) domain has a high pi(non)/pi(syn) ratio, and the low polymorphism of the Pi-ta gene might have primarily been caused by recurrent selective sweep and constraint by other putative physiological functions. Meanwhile, we provide data to show that the amino acid Ala-918 of H1 in the LRR domain has a close relationship with the resistant phenotype. H1 might have recently arisen during rice domestication and may be associated with the scenario of a blast pathogen-host shift from Italian millet to rice.     

Genetic diversity among perennial wild rice Oryza rufipogon Griff., in the Mekong Delta

Oryza rufipogon Griff. is a perennial species of wild rice widely distributed along the channels and rivers of the Mekong Delta, Vietnam. This study attempted to find centers of diversity among wild rice populations in this area and their inter‐relationships. The highest genetic diversity was found in the Dong Thap population and the lowest in the Can Tho population. Maternal diversity evaluated using chloroplast INDELs detected ten plastid types, five of which were novel relative to other Asian countries. The mitochondrial genome suggested two unique deletions. One 699‐bp deletion via short tandem repeats was accompanied by another deletion including orf153. All accessions carrying the mitochondrial type were found in a particular plastid type. This unique maternal lineage was confined to specific channels where it showed vigorous vegetative growth in comparison to upstream areas where various maternal lineages and maximum genetic diversity occurred. This area along the Mekong Delta is a center of not only nuclear but also maternal diversity.     

Comparison of the genetic diversity of common wild rice (Oryza rufipogon Griff.) and cultivated rice (O. sativa L.) using RFLP markers

Forty fourth single-copy RFLP markers were used to evaluate the genetic diversity of 122 accessions of common wild rice (CWR, Oryza rufipogon Griff.) and 75 entries of cultivated rice (Oryza sativa L. ) from more than ten Asian countries. A comparison of the parameters showing genetic diversity, including the percentage of polymorphic loci (P), the average number of alleles per locus (A), the number of genotypes (Ng), the average heterozygosity (Ho) and the average genetic multiplicity (Hs) of CWR and indica and japonica subspecies of cultivated rice from different countries and regions, indicated that CWR from China possesses the highest genetic diversity, followed by CWR from South Asia and Southeast Asia. The genetic diversity of CWR from India is the second highest. Although the average gene diversity (Hs)of the South Asian CWR is higher than that of the Southeast Asian CWR, its percentage of polymorphic loci (P), number of alleles (Na) and number of genotypes (Ng) are all smaller. It was also found that the genetic diversity of cultivated rice is obviously lower than that of CWR. At the 44 loci investigated, the number of polymorphic loci of cultivated rice is only 3/4 that of CWR, while the number of alleles, 60%, and the number of genotypes is about 1/2 that of CWR. Of the two subspecies studied, the genetic diversity of indica is higher than that of japonica. The average heterozygosity of the Chinese CWR is the highest among all the entries studied. The average heterozygosity of CWR is about two-times that of cultivated rice. It is suggested that during the course of evolution from wild rice to cultivated rice, many alleles were lost through natural and human selection, leading to the lower heterozygosity and genetic diversity of the cultivated rice. Received: 19 May 1999 / Accepted: 26 April 2000     

野生稻Oryza nivara和O. rufipogon DNA甲基化多样性

本文调查研究了野生稻群体内及群体间的DNA甲基化多样性。选取与亚洲栽培稻近缘的两个野生种Oryza nivara和O. rufipogon作为研究对象, 采用改进的MSAP (methylation-sensitive amplification polymorphism)技术对其基因组CCGG位点的甲基化多样性进行了分析。结果表明: 在同一个IRGC(the International Rice Germplasm Center)编号群体内的不同个体间, 基因组甲基化条带高度一致; 而在不同编号群体间, 甲基化条带表现为多态。其中后者又可以分为两类: 条带模式高度一致的Class I和条带模式呈多态性的Class II。将上述两类甲基化片段的编码基因与栽培稻粳稻(O. sativaL. subsp. japonica)和籼稻(O. sativa L. subsp. indica)两个亚种的同源基因进行序列比对发现, 在进化趋势上Class I表现得比较保守, 而Class II较为活跃。DNA甲基化多样性作为标志遗传多样性的一种信息来源, 其在群体分化及物种进化过程中的作用还需要进一步探讨。     

Identification and characterization of Australian wild rice strains of Oryza meridionalis and Oryza rufipogon by SINE insertion polymorphism

Of the rice species with an AA genome, Oryza meridionalis has been identified in northern Australia as a species of the annual type, among those previously classified as Oryza perennis, Oryza rufipogon or Oryza nivara. This notion has, however, led to some confusion to determine which strains belong to O. meridionalis and how different these strains are from the O. rufipogon strains of the annual type. In this paper, we examined Australian wild rice strains for the presence or absence of p-SINE1 members, which have been used for identification of the strains of species with the AA genome, by PCR using primers that hybridize to the sequences flanking each p-SINE1 member. The rice strains examined include perennial and annual strains, which have previously been described as O. rufipogon. We found that all the annual strains and other strains, whose types have not been determined, have p-SINE1 members that are specifically present at the corresponding loci in the standard strains of O. meridionalis, but do not have those which are specifically present at the corresponding loci in the strains of the other species with the AA genome. The perennial strains, however, have p-SINE1 members that are specifically present at the corresponding loci in the standard O. rufipogon strains of either the annual or the perennial type, but do not have those which are specifically present at the corresponding loci in the strains of the other species with the AA genome, including O. meridionalis. These findings support the previous notion that O. meridionalis consists of the annual strains and is a distinct species from O. rufipogon. The p-SINE1 members used in this study appear to be very useful for classification of any wild rice strains of the AA-genome species, even when one has limited knowledge of morphology, taxonomy, physiology, and biochemistry of rice strains.