Sampath Perumal

Complete chloroplast and ribosomal sequences for 30 accessions elucidate evolution of Oryza AA genome species.

Cytoplasmic chloroplast (cp) genomes and nuclear ribosomal DNA (nR) are the primary sequences used to understand plant diversity and evolution. We introduce a high-throughput method to simultaneously obtain complete cp and nR sequences using Illumina platform whole-genome sequence. We applied the method to 30 rice specimens belonging to nine Oryza species. Concurrent phylogenomic analysis using cp and nR of several of specimens of the same Oryza AA genome species provides insight into the evolution and domestication of cultivated rice, clarifying three ambiguous but important issues in the evolution of wild Oryza species. First, cp-based trees clearly classify each lineage but can be biased by inter-subspecies cross-hybridization events during speciation. Second, O. glumaepatula, a South American wild rice, includes two cytoplasm types, one of which is derived from a recent interspecies hybridization with O. longistminata. Third, the Australian O. rufipogan-type rice is a perennial form of O. meridionalis.

Genotyping-by-sequencing map permits identification of clubroot resistance QTLs and revision of the reference genome assembly in cabbage (Brassica oleracea L.)

Clubroot is a devastating disease caused by Plasmodiophora brassicae and results in severe losses of yield and quality in Brassica crops. Many clubroot resistance genes and markers are available in Brassica rapa but less is known in Brassica oleracea. Here, we applied the genotyping-by-sequencing (GBS) technique to construct a high-resolution genetic map and identify clubroot resistance (CR) genes. A total of 43,821 SNPs were identified using GBS data for two parental lines, one resistant and one susceptible lines to clubroot, and 18,187 of them showed >5× coverage in the GBS data. Among those, 4,103 were credibly genotyped for all 78 F2 individual plants. These markers were clustered into nine linkage groups spanning 879.9 cM with an average interval of 1.15 cM. Quantitative trait loci (QTLs) survey based on three rounds of clubroot resistance tests using F2 : 3 progenies revealed two and single major QTLs for Race 2 and Race 9 of P. brassicae, respectively. The QTLs show similar locations to the previously reported CR loci for Race 4 in B. oleracea but are in different positions from any of the CR loci found in B. rapa. We utilized two reference genome sequences in this study. The high-resolution genetic map developed herein allowed us to reposition 37 and 2 misanchored scaffolds in the 02–12 and TO1000DH genome sequences, respectively. Our data also support additional positioning of two unanchored 3.3 Mb scaffolds into the 02–12 genome sequence.

Microsatellite-based analysis of genetic diversity in 91 commercial Brassica oleracea L. cultivars belonging to six varietal groups

Brassica oleracea L. includes various types of important vegetables that show extremely diverse phenotypes. To elucidate the genetic diversity and relationships among commercial cultivars derived by different companies throughout the world, we characterized the diversity and genetic structure of 91 commercial B. oleracea cultivars belonging to six varietal groups, including cabbage, broccoli, cauliflower, kohlrabi, kale and kai-lan. We used 69 polymorphic microsatellite markers showing a total of 359 alleles with an average number of 5.20 alleles per locus. Polymorphism information content (PIC) values ranged from 0.06 to 0.73, with an average of 0.40. Among the six varietal groups, kohlrabi cultivars exhibited the highest heterozygosity level, whereas kale cultivars showed the lowest. Based on genetic similarity values, an UPGMA clustering dendrogram and a two-dimensional scale diagram (PCoA) were generated to analyze genetic diversity. The cultivars were clearly separated into six different clusters with a tendency to cluster into varietal groups. Model-based structure analysis revealed six genetic groups, in which cabbage cultivars were divided into two subgroups that were differentiated by their head shape, whereas cauliflower and kai-lan cultivars clustered together into a single group. Furthermore, we identified 18 SSR markers showing 27 unique alleles specific to only one cultivar that can be used to discriminate 22 cultivars from the others. Our phylogenetic and population structure analysis presents new insights into the genetic structure and relationships among 91 B. oleracea cultivars and provides valuable information for breeding of B. oleracea species. In addition, we demonstrate the utility of SSR markers as a powerful tool for discriminating between the cultivars. The SSR markers described herein will also be helpful for Distinctness, Uniformity and Stability (DUS) test of new cultivars.

Genome and evolution of the shade‐requiring medicinal herb Panax ginseng

Summary Panax ginseng C. A. Meyer, reputed as the king of medicinal herbs, has slow growth, long generation time, low seed production and complicated genome structure that hamper its study. Here, we unveil the genomic architecture of tetraploid P. ginseng by de novo genome assembly, representing 2.98 Gbp with 59 352 annotated genes. Resequencing data indicated that diploid Panax species diverged in association with global warming in Southern Asia, and two North American species evolved via two intercontinental migrations. Two whole genome duplications (WGD) occurred in the family Araliaceae (including Panax) after divergence with the Apiaceae, the more recent one contributing to the ability of P. ginseng to overwinter, enabling it to spread broadly through the Northern Hemisphere. Functional and evolutionary analyses suggest that production of pharmacologically important dammarane‐type ginsenosides originated in Panax and are produced largely in shoot tissues and transported to roots; that newly evolved P. ginseng fatty acid desaturases increase freezing tolerance; and that unprecedented retention of chlorophyll a/b binding protein genes enables efficient photosynthesis under low light. A genome‐scale metabolic network provides a holistic view of Panax ginsenoside biosynthesis. This study provides valuable resources for improving medicinal values of ginseng either through genomics‐assisted breeding or metabolic engineering.

The complete chloroplast genome of two Brassica species, Brassica nigra and B. Oleracea

Abstract The two Brassica species, Brassica nigra and Brassica oleracea, are important agronomic crops. The chloroplast genome sequences were generated by de novo assembly using whole genome next-generation sequences. The chloroplast genomes of B. nigra and B. oleracea were 153 633 bp and 153 366 bp in size, respectively, and showed conserved typical chloroplast structure. The both chloroplast genomes contained a total of 114 genes including 80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Phylogenetic analysis revealed that B. oleracea is closely related to B. rapa and B. napus but B. nigra is more diverse than the neighbor species Raphanus sativus.

Repeat Evolution in Brassica rapa (AA), B. oleracea (CC), and B. napus (AACC) Genomes

The genus Brassica is an important resource for major agricultural products such as oils, vegetable and fodder. The Brassiceae tribe-specific whole-genome triplication that occurred ~15.9 million years ago influenced the speciation and morphological diversification that has been exploited in agriculture, making Brassica an excellent model system for studying polyploidization- mediated evolution. Genome sequencing and comparative genome analysis have revealed conserved structures and uncovered the genome evolution of Brassica species. While chromosome shuffling and asymmetric subgenome gene retention are widely reported in Brassica species, limited information is available about the dynamics of repetitive elements (REs), which are central to epigenetic mechanisms and thus play a pivotal role in plant genome adaptation and evolution. The assembled reference genome sequences of B. rapa (AA) and B. oleracea (CC), and their derived allotetraploid, B. napus (AACC), cover 58%, 86%, and 75% of their respective estimated genome sizes. The remaining non-assembled genome portions vary between these three genome sequences, and the major components remain hidden in each genome. Here, we review the dynamics of the major Brassica repeats that have played roles in speciation of the AA, CC, and AACC genomes. We show that 10 major Brassica repeats appear to occupy more than 50% of each respective unassembled genome sequence, yet represent less than 1% of assembled reference genome sequences. We have estimated their genome proportions using whole-genome Illumina reads and cytogenetic analyses in an attempt to understand the role of these repeats in genome evolution.

Elucidating the major hidden genomic components of the A, C, and AC genomes and their influence on Brassica evolution

Decoding complete genome sequences is prerequisite for comprehensive genomics studies. However, the currently available reference genome sequences of Brassica rapa (A genome), B. oleracea (C) and B. napus (AC) cover 391, 540, and 850 Mbp and represent 80.6, 85.7, and 75.2% of the estimated genome size, respectively, while remained are hidden or unassembled due to highly repetitive nature of these genome components. Here, we performed the first comprehensive genome-wide analysis using low-coverage whole-genome sequences to explore the hidden genome components based on characterization of major repeat families in the B. rapa and B. oleracea genomes. Our analysis revealed 10 major repeats (MRs) including a new family comprising about 18.8, 10.8, and 11.5% of the A, C and AC genomes, respectively. Nevertheless, these 10 MRs represented less than 0.7% of each assembled reference genome. Genomic survey and molecular cytogenetic analyses validates our insilico analysis and also pointed to diversity, differential distribution, and evolutionary dynamics in the three Brassica species. Overall, our work elucidates hidden portions of three Brassica genomes, thus providing a resource for understanding the complete genome structures. Furthermore, we observed that asymmetrical accumulation of the major repeats might be a cause of diversification between the A and C genomes.

Rapid amplification of four retrotransposon families promoted speciation and genome size expansion in the genus Panax

Genome duplication and repeat multiplication contribute to genome evolution in plants. Our previous work identified a recent allotetraploidization event and five high-copy LTR retrotransposon (LTR-RT) families PgDel, PgTat, PgAthila, PgTork, and PgOryco in Panax ginseng. Here, using whole-genome sequences, we quantified major repeats in five Panax species and investigated their role in genome evolution. The diploids P. japonicus, P. vietnamensis, and P. notoginseng and the tetraploids P. ginseng and P. quinquefolius were analyzed alongside their relative Aralia elata. These species possess 0.8–4.9 Gb haploid genomes. The PgDel, PgTat, PgAthila, and PgTork LTR-RT superfamilies accounted for 39–52% of the Panax species genomes and 17% of the A. elata genome. PgDel included six subfamily members, each with a distinct genome distribution. In particular, the PgDel1 subfamily occupied 23–35% of the Panax genomes and accounted for much of their genome size variation. PgDel1 occupied 22.6% (0.8 Gb of 3.6 Gb) and 34.5% (1.7 Gb of 4.9 Gb) of the P. ginseng and P. quinquefolius genomes, respectively. Our findings indicate that the P. quinquefolius genome may have expanded due to rapid PgDel1 amplification over the last million years as a result of environmental adaptation following migration from Asia to North America.

The complete chloroplast genome sequence of Indian barnyard millet, Echinochloa frumentacea (Poaceae)

Abstract We report the complete chloroplast genome sequence of the important nutritious millet crop Indian barnyard millet, Echinochloa frumentacea Link (Poaceae). The size of the circular chloroplast genome is 139 593 bp in length with a typical quadripartite structure, containing pair of inverted repeats of 22 618 bp, flanked by large single copy and small single copy regions of 81 839 bp, 12 518 bp, respectively. Overall GC content of the genome was 38.6% and consists of 112 individual genes, including 77 protein coding genes, 30 tRNA genes, four rRNA genes and one conserved open reading frame. In addition, phylogenetic analysis with grass species has revealed that E. frumentacea is diverged around 1.9–2.7 million years with its close relatives, E. oryzicola and E. crus-galli.

Re-exploration of U's Triangle Brassica Species Based on Chloroplast Genomes and 45S nrDNA Sequences

The concept of U’s triangle, which revealed the importance of polyploidization in plant genome evolution, described natural allopolyploidization events in Brassica using three diploids [B. rapa (A genome), B. nigra (B), and B. oleracea (C)] and derived allotetraploids [B. juncea (AB genome), B. napus (AC), and B. carinata (BC)]. However, comprehensive understanding of Brassica genome evolution has not been fully achieved. Here, we performed low-coverage (2–6×) whole-genome sequencing of 28 accessions of Brassica as well as of Raphanus sativus [R genome] to explore the evolution of six Brassica species based on chloroplast genome and ribosomal DNA variations. Our phylogenomic analyses led to two main conclusions. (1) Intra-species-level chloroplast genome variations are low in the three allotetraploids (2~7 SNPs), but rich and variable in each diploid species (7~193 SNPs). (2) Three allotetraploids maintain two 45SnrDNA types derived from both ancestral species with maternal dominance. Furthermore, this study sheds light on the maternal origin of the AC chloroplast genome. Overall, this study clarifies the genetic relationships of U’s triangle species based on a comprehensive genomics approach and provides important genomic resources for correlative and evolutionary studies.