Genome size variation among and within Camellia species by using flow cytometric analysis

Insights into the Genomic Background of Nine Common Chinese Medicinal Plants by Flow Cytometry and Genome Survey

Medicinal plants have long played a crucial role in healthcare systems, but limited genomic information on these species has impeded the integration of modern biological technologies into medicinal plant research. In this study, we selected nine common medicinal plants, each belonging to a different plant family, including Sarcandra glabra (Chloranthaceae), Nekemias grossedentata (Vitaceae), Uraria crinita (Fabaceae), Gynostemma pentaphyllum (Cucurbitaceae), Reynoutria japonica (Polygonaceae), Pseudostellaria heterophylla (Caryophyllaceae), Morinda officinalis (Rubiaceae), Vitex rotundifolia (Lamiaceae), and Gynura formosana (Asteraceae), to estimate their genome sizes and conduct preliminary genomic surveys. The estimated genome sizes by flow cytometry were 3.66 Gb, 0.65 Gb, 0.58 Gb, 1.02 Gb, 3.96 Gb, 2.99 Gb, 0.43 Gb, 0.78 Gb, and 7.27 Gb, respectively. The genome sizes of M. officinalis, R. japonica, and G. pentaphyllum have been previously reported. Comparative analyses suggest that variations in genome size may arise due to differences in measurement methods and sample sources. Therefore, employing multiple approaches to assess genome size is necessary to provide more reliable information for further genomic research. Based on the genome survey, species with considerable genome size variation or polyploidy, such as G. pentaphyllum, should undergo a ploidy analysis in conjunction with population genomics studies to elucidate the development of the diversified genome size. Additionally, a genome survey of U. crinita, a medicinal plant with a relatively small genome size (509.08 Mb) and of considerable interest in southern China, revealed a low heterozygosity rate (0.382%) and moderate repeat content (51.24%). Given the limited research costs, this species represents a suitable candidate for further genomic studies on Leguminous medicinal plants characteristic of southern China. This foundational genomic information will serve as a critical reference for the sustainable development and utilization of these medicinal plants.

Evolutionary Histories of Camellia japonica and Camellia rusticana

The genus Camellia is widely distributed, primarily in East Asia. Camellia japonica is located at the northern limit of this genus distribution, and understanding changes in its distribution is crucial for understanding the evolution of plants in this region, as well as their relationship with geological history and climate change. Moreover, the classification of sect. Camellia in Japan has not been clarified. Therefore, this study aims to understand the evolutionary history of the Japanese sect. Camellia. The genetic population structure was analysed using SNP data and MIG-seq. The relationship between the Japanese sect. Camellia, including the related species in China, was further inferred from the phylogeny generated by RA x ML, SplitsTree and PCA. Population genetic structure was inferred using a Bayesian clustering method (ADMIXTURE). We subsequently employed approximate Bayesian computation, which was further supported by the coalescent simulations (DIYABC, fastsimcoal and Bayesian Skyline Plots) to explore the changes in population, determining which events appropriately explain the phylogeographical signature. Ecological niche modelling was combined with genetic analyses to compare current and past distributions. The analyses consistently showed that C. japonica and C. rusticana are distinct, having diverged from each other during the Middle to Late Miocene period. Furthermore, C. japonica differentiated into four major populations (North, South, Ryukyu-Taiwan and Continent). The Japanese sect. Camellia underwent speciation during archipelago formation, reflecting its ancient evolutionary history compared with other native Japanese plants. C. rusticana did not diverge from C. japonica in snow-rich environments during the Quaternary period. Our results suggest that both species have been independent since ancient times and that ancestral populations of C. japonica have persisted in northern regions. Furthermore, the C. japonica population on the continent is hypothesised to have experienced a reverse-colonisation event from southern Japan during the late Pleistocene glaciation.

OHDLF: A Method for Selecting Orthologous Genes for Phylogenetic Construction and Its Application in the Genus Camellia

Accurate phylogenetic tree construction for species without reference genomes often relies on de novo transcriptome assembly to identify single-copy orthologous genes. However, challenges such as whole-genome duplication (WGD), heterozygosity, gene duplication, and loss can hinder the selection of these genes, leading to limited data for constructing reliable species trees. To address these issues, we developed a new analytical pipeline, OHDLF (Orthologous Haploid Duplication and Loss Filter), which filters orthologous genes from transcript data and adapts parameter settings based on genomic characteristics for further phylogenetic tree construction. In this study, we applied OHDLF to the genus Camellia and evaluated its effectiveness in constructing phylogenetic trees. The results highlighted the pipeline's ability to handle challenges like high heterozygosity and recent gene duplications by selectively retaining genes with a missing rate and merging duplicates with high similarity. This approach ensured the preservation of informative sites and produced a highly supported consensus tree for Camellia. Additionally, we evaluate the accuracy of the OHDLF phylogenetic trees for different species, demonstrating that the OHDLF pipeline provides a flexible and effective method for selecting orthologous genes and constructing accurate phylogenetic trees, adapting to the genomic characteristics of various plant groups.

The complex hexaploid oil-Camellia genome traces back its phylogenomic history and multi-omics analysis of Camellia oil biosynthesis

Oil-Camellia (Camellia oleifera), belonging to the Theaceae family Camellia, is an important woody edible oil tree species. The Camellia oil in its mature seed kernels, mainly consists of more than 90% unsaturated fatty acids, tea polyphenols, flavonoids, squalene and other active substances, which is one of the best quality edible vegetable oils in the world. However, genetic research and molecular breeding on oil-Camellia are challenging due to its complex genetic background. Here, we successfully report a chromosome-scale genome assembly for a hexaploid oil-Camellia cultivar Changlin40. This assembly contains 8.80 Gb genomic sequences with scaffold N50 of 180.0 Mb and 45 pseudochromosomes comprising 15 homologous groups with three members each, which contain 135 868 genes with an average length of 3936 bp. Referring to the diploid genome, intragenomic and intergenomic comparisons of synteny indicate homologous chromosomal similarity and changes. Moreover, comparative and evolutionary analyses reveal three rounds of whole-genome duplication (WGD) events, as well as the possible diversification of hexaploid Changlin40 with diploid occurred approximately 9.06 million years ago (MYA). Furthermore, through the combination of genomics, transcriptomics and metabolomics approaches, a complex regulatory network was constructed and allows to identify potential key structural genes (SAD, FAD2 and FAD3) and transcription factors (AP2 and C2H2) that regulate the metabolism of Camellia oil, especially for unsaturated fatty acids biosynthesis. Overall, the genomic resource generated from this study has great potential to accelerate the research for the molecular biology and genetic improvement of hexaploid oil-Camellia, as well as to understand polyploid genome evolution.

Genome assembly of the rare and endangered Grantham's camellia, Camellia granthamiana

Grantham's camellia (Camellia granthamiana Sealy) is a rare and endangered tea species discovered in Hong Kong in 1955 and endemic to southern China. Despite its high conservation value, the genomic resources of C. granthamiana are limited. Here, we present a chromosome-scale draft genome of the tetraploid C. granthamiana (2n = 4x = 60), combining PacBio long-read sequencing and Omni-C data. The assembled genome size is ∼2.4 Gb, with most sequences anchored to 15 pseudochromosomes resembling a monoploid genome. The genome has high contiguity, with a scaffold N50 of 139.7 Mb, and high completeness (97.8% BUSCO score). Our gene model prediction resulted in 68,032 protein-coding genes (BUSCO score of 90.9%). We annotated 1.65 Gb of repeat content (68.48% of the genome). Our Grantham's camellia genome assembly is a valuable resource for investigating Grantham's camellia's biology, ecology, and phylogenomic relationships with other Camellia species, and provides a foundation for further conservation measures.

Genome size, chromosome number variation and its correlation with stomatal characters for assessment of ploidy levels in a core subset of mulberry (Morus spp.) germplasm

The large size of the germplasm collection along with scanty information on their cytological and genome constitution have hindered well-planned breeding schemes in mulberry. To address the issue, a study was undertaken to investigate the variability in DNA content and genome size, chromosome number, ploidy and its relation with important stomatal characteristics among 162 mulberry germplasm. These germplasm comprise a core subset of 150 collections along with a representative collection of different mulberry species including the wild. Among the germplasm belonging to 16 species, we identified 122 diploids (2n = 28), 4 aneuploids (2n = 30), 13 triploids (2n = 42), 15 tetraploids (2n = 56), 7 hexaploids (2n = 84) and 1 dodecosaploid (2n = 308) based on the chromosome count. Most of the cultivated mulberries are found to be diploids. The mean nuclear 2C DNA content estimated by Flow cytometry, varied from 0.723±0.006 pg (M. australis, 2n = 2x) to 7.732 pg (M. nigra, 2n = 22x). The 2C DNA content positively correlated with the ploidy status and stomatal length (r = 0.814, p<0.001). Based on the 1Cx value, the study also suggests that the majority of the polyploid species have experienced genome downsizing in relation to their diploid progenitors. This study provides the most essential information on chromosome number, ploidy and DNA content to facilitate the utilization of a core subset of germplasm in the mulberry breeding program.

Genetic differentiation and genetic structure of mixed-ploidy Camellia hainanica populations

Camellia hainanica, which is common in China's Hainan Province, is an important woody olive tree species. Due to many years of geographic isolation, C. hainanica has not received the attention it deserves, which limits the exploitation of germplasm resources. Therefore, it is necessary to study population genetic characteristics for further utilization and conservation of C. hainanica. In this study, 96 individuals in six wild Camellia hainanica populations were used for ploidy analysis of the chromosome number, and the genetic diversity and population structure were investigated using 12 pairs of SSR primers. The results show complex ploidy differentiation in C. hainanica species. The ploidy of wild C. hainanica includes tetraploid, pentaploid, hexaploid, heptaploid, octoploid and decaploid species. Genetic analysis shows that genetic diversity and genetic differentiation among populations are low. Populations can be divided into two clusters based on their genetic structure, which matches their geographic location. Finally, to further maintain the genetic diversity of C. hainanica, ex-situ cultivation and in-situ management measures should be considered to protect it in the future.

Genomic and genetic advances of oiltea-camellia (Camellia oleifera)

Oiltea-camellia (C. oleifera) is a widely cultivated woody oil crop in Southern China and Southeast Asia. The genome of oiltea-camellia was very complex and not well explored. Recently, genomes of three oiltea-camellia species were sequenced and assembled, multi-omic studies of oiltea-camellia were carried out and provided a better understanding of this important woody oil crop. In this review, we summarized the recent assembly of the reference genomes of oiltea-camellia, genes related to economic traits (flowering, photosynthesis, yield and oil component), disease resistance (anthracnose) and environmental stress tolerances (drought, cold, heat and nutrient deficiency). We also discussed future directions of integrating multiple omics for evaluating genetic resources and mining key genes of important traits, and the application of new molecular breeding and gene editing technologies to accelerate the breeding process of oiltea-camellia.

Comparative Genomic Analysis Uncovers the Chloroplast Genome Variation and Phylogenetic Relationships of Camellia Species

Camellia is the largest genus in the family Theaceae. Due to phenotypic diversity, frequent hybridization, and polyploidization, an understanding of the phylogenetic relationships between Camellia species remains challenging. Comparative chloroplast (cp) genomics provides an informative resource for phylogenetic analyses of Camellia. In this study, 12 chloroplast genome sequences from nine Camellia species were determined using Illumina sequencing technology via de novo assembly. The cp genome sizes ranged from 156,545 to 157,021 bp and were organized into quadripartite regions with the typical angiosperm cp genomes. Each genome harbored 87 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes in the same order and orientation. Differences in long and short sequence repeats, SNPs, and InDels were detected across the 12 cp genomes. Combining with the complete cp sequences of seven other species in the genus Camellia, a total of nine intergenic sequence divergent hotspots and 14 protein-coding genes with high sequence polymorphism were identified. These hotspots, especially the InDel (~400 bp) located in atpH-atpI region, had sufficient potential to be used as barcode markers for further phylogenetic analysis and species identification. Principal component and phylogenetic analysis suggested that regional constraints, rather than functional constraints, strongly affected the sequence evolution of the cp genomes in this study. These cp genomes could facilitate the development of new molecular markers, accurate species identification, and investigations of the phylogenomic relationships of the genus Camellia.

Altitudinal Patterns in Adaptive Evolution of Genome Size and Inter-Genome Hybridization Between Three Elymus Species From the Qinghai–Tibetan Plateau

Genome size variation and hybridization occur frequently within or between plant species under diverse environmental conditions, which enrich species diversification and drive the evolutionary process. Elymus L. is the largest genus in Triticeae with five recognized basic genomes (St, H, P, W, and Y). However, the data on population cytogenetics of Elymus species are sparse, especially whether genome hybridization and chromosomal structure can be affected by altitude are still unknown. In order to explore the relationship between genome sizes, we studied interspecific hybridization and altitude of Elymus species at population genetic and cytological levels. Twenty-seven populations at nine different altitudes (2,800–4,300 m) of three Elymus species, namely, hexaploid E. nutans (StHY, 2n = 6x = 42), tetraploid E. burchan-buddae (StY, 2n = 4x = 28), and E. sibiricus (StH, 2n = 4x = 28), were sampled from the Qinghai–Tibetan Plateau (QTP) to estimate whether intraspecific variation could affect the genomic relationships by genomic in situ hybridization (GISH), and quantify the genome size of Elymus among different altitude ecological groups by flow cytometry. The genome size of E. nutans, E. burchan-buddae, and E. sibiricus varied from 12.38 to 22.33, 8.81 to 18.93, and 11.46 to 20.96 pg/2C with the averages of 19.59, 12.39, and 16.85 pg/2C, respectively. The curve regression analysis revealed a strong correlation between altitude and nuclear DNA content in three Elymus species. In addition, the chromosomes of the St and Y genomes demonstrated higher polymorphism than that of the H genome. Larger genome size variations occurred in the mid-altitude populations (3,900–4,300 m) compared with other-altitude populations, suggesting a notable altitudinal pattern in genome size variation, which shaped genome evolution by altitude. This result supports our former hypothesis that genetic richness center at medium altitude is useful and valuable for species adaptation to highland environmental conditions, germplasm utilization, and conservation.

Admixture of divergent genomes facilitates hybridization across species in the family Brassicaceae

Hybridization and polyploidization are pivotal to plant evolution. Genetic crosses between distantly related species are rare in nature due to reproductive barriers but how such hurdles can be overcome is largely unknown. Here we report the hybrid genome structure of xBrassicoraphanus, a synthetic allotetraploid of Brassica rapa and Raphanus sativus. We performed cytogenetic analysis and de novo genome assembly to examine chromosome behaviors and genome integrity in the hybrid. Transcriptome analysis was conducted to investigate expression of duplicated genes in conjunction with epigenome analysis to address whether genome admixture entails epigenetic reconfiguration. Allotetraploid xBrassicoraphanus retains both parental chromosomes without genome rearrangement. Meiotic synapsis formation and chromosome exchange are avoided between nonhomologous progenitor chromosomes. Reconfiguration of transcription network occurs, and less divergent cis-elements of duplicated genes are associated with convergent expression. Genome-wide DNA methylation asymmetry between progenitors is largely maintained but, notably, B. rapa-originated transposable elements are transcriptionally silenced in xBrassicoraphanus through gain of DNA methylation. Our results demonstrate that hybrid genome stabilization and transcription compatibility necessitate epigenome landscape adjustment and rewiring of cis-trans interactions. Overall, this study suggests that a certain extent of genome divergence facilitates hybridization across species, which may explain the great diversification and expansion of angiosperms during evolution.

Chromosome-level genome of Camellia lanceoleosa provides a valuable resource for understanding genome evolution and self-incompatibility

The section Oleifera (Theaceae) has attracted attention for the high levels of unsaturated fatty acids found in its seeds. Here, we report the chromosome-scale genome of the sect. Oleifera using diploid wild Camellia lanceoleosa with a final size of 3.00 Gb and an N50 scaffold size of 186.43 Mb. Repetitive sequences accounted for 80.63% and were distributed unevenly across the genome. Camellia lanceoleosa underwent a whole-genome duplication event approximately 65 million years ago (65 Mya), prior to the divergence of C. lanceoleosa and Camellia sinensis (approx. 6-7 Mya). Syntenic comparisons of these two species elucidated the genomic rearrangement, appearing to be driven in part by the activity of transposable elements. The expanded and positively selected genes in C. lanceoleosa were significantly enriched in oil biosynthesis, and the expansion of homomeric acetyl-coenzyme A carboxylase (ACCase) genes and the seed-biased expression of genes encoding heteromeric ACCase, diacylglycerol acyltransferase, glyceraldehyde-3-phosphate dehydrogenase and stearoyl-ACP desaturase could be of primary importance for the high oil and oleic acid content found in C. lanceoleosa. Theanine and catechins were present in the leaves of C. lanceoleosa. However, caffeine can not be dectected in the leaves but was abundant in the seeds and roots. The functional and transcriptional divergence of genes encoding SAM-dependent N-methyltransferases may be associated with caffeine accumulation and distribution. Gene expression profiles, structural composition and chromosomal location suggest that the late-acting self-incompatibility of C. lanceoleosa is likely to have favoured a novel mechanism co-occurring with gametophytic self-incompatibility. This study provides valuable resources for quantitative and qualitative improvements and genome assembly of polyploid plants in sect. Oleifera.

Genetic Diversity and Genome Size Variability in the Russian Genebank Collection of Tea Plant [Camellia sinensis (L). O. Kuntze]

The tea collection of the FRC SSC RAS (Sochi, Maykop in Russia) represents one of the northernmost germplasm comprising a number of locally derived cultivars and ɣ-irradiation mutants. The latter are often characterized by larger genome size, which may lead to better adaptation to biotic and abiotic stress. Such genotypes may be a valuable genetic resource for better adaptability to extreme environmental conditions, which could enable tea cultivation outside global growing regions. Microsatellite markers are often the best choice for genetic diversity analysis in genebank collections. However, their use in polyploid species is questionable because simple sequence repeat (SSR) allele dosage cannot be readily determined. Therefore, the efficiency of SSR and start codon targeted (SCoT) markers was investigated using 43 selected cultivars from the Russian genebank collection derived from mutant breeding and clonal selection. Previously, the increase in genome size was confirmed in 18 mutants within this collection. Despite the presence of polyploid tea genotypes, our study revealed higher efficiency of SSR markers than SCoT markers. Subsequent SSR analysis of the 106 genotypes in the Russian genebank collection revealed three distinct genetic clusters after STRUCTURE analysis. Greater genetic variation was observed within genetic clusters than between clusters, indicating low genetic variation between collections. Nevertheless, the northernmost tea collection exhibited a greater genetic distance from the other two clusters than they did from each other. Close genetic relationships were found between many cultivars with particularly large leaves and mutant forms. Pearson's correlation analysis revealed a significant, moderate correlation between genome size and leaf area size. Our study shows that microsatellite fingerprinting is useful to estimate the genetic diversity and genetic background of tea germplasm in Russia despite polyploid tea accessions. Thus, the results of our study contribute to the development of future tea germplasm conservation strategies and modern tea breeding programs.

Comparative transcriptomic analysis reveals genes related to the rapid accumulation of oleic acid in Camellia chekiangoleosa, an oil tea plant with early maturity and large fruit

Camellia chekiangoleosa has a higher oleic acid content and a shorter reproductive cycle than typical oil tea plants. It was intensively sampled over six C. chekiangoleosa seed development stages. The content of fatty acids determined by GC showed that the accumulation of fatty acids gradually increased from the S1 to S5 stages, and the maximum concentration was reached in S5. Then, fatty acids declined slightly in S6. The main fatty acid component showed the same accumulation trend as the total fatty acids, except linolenic acid, which remained at a low level throughout seed developmental stages. Changes in the expression of fatty acid accumulation-related genes were monitored using second-generation and SMRT full-length transcriptome sequencing. Finally, 18.92 G accurate and reliable data were obtained. Differential expression analysis and weighted coexpression analysis revealed two "gene modules" significantly associated with oleic acid and linoleic acid contents, and the high expression of ENR, KAS I, and KAS II, which accumulate substrates for oleic acid synthesis, was thought to be responsible for the rapid accumulation of fatty acids in the early stage. The rapid increase in fatty acids in the second stage may be closely related to the synergy between the high expression of SAD and low expression of FAD2. In addition, many transcription factors, such as ERF, GRAS, GRF, MADS, MYB and WRKY, may be involved in the fatty acid synthesis. Our data provide a rich resource for further studies on the regulation of fatty acid synthesis in C. chekiangoleosa.

Genome Size Variation and Evolution Driven by Transposable Elements in the Genus Oryza

Genome size variation and evolutionary forces behind have been long pursued in flowering plants. The genus Oryza, consisting of approximately 25 wild species and two cultivated rice, harbors eleven extant genome types, six of which are diploid (AA, BB, CC, EE, FF, and GG) and five of which are tetraploid (BBCC, CCDD, HHJJ, HHKK, and KKLL). To obtain the most comprehensive knowledge of genome size variation in the genus Oryza, we performed flow cytometry experiments and estimated genome sizes of 166 accessions belonging to 16 non-AA genome Oryza species. k-mer analyses were followed to verify the experimental results of the two accessions for each species. Our results showed that genome sizes largely varied fourfold in the genus Oryza, ranging from 279 Mb in Oryza brachyantha (FF) to 1,203 Mb in Oryza ridleyi (HHJJ). There was a 2-fold variation (ranging from 570 to 1,203 Mb) in genome size among the tetraploid species, while the diploid species had 3-fold variation, ranging from 279 Mb in Oryza brachyantha (FF) to 905 Mb in Oryza australiensis (EE). The genome sizes of the tetraploid species were not always two times larger than those of the diploid species, and some diploid species even had larger genome sizes than those of tetraploids. Nevertheless, we found that genome sizes of newly formed allotetraploids (BBCC-) were almost equal to totaling genome sizes of their parental progenitors. Our results showed that the species belonging to the same genome types had similar genome sizes, while genome sizes exhibited a gradually decreased trend during the evolutionary process in the clade with AA, BB, CC, and EE genome types. Comparative genomic analyses further showed that the species with different rice genome types may had experienced dissimilar amplification histories of retrotransposons, resulting in remarkably different genome sizes. On the other hand, the closely related rice species may have experienced similar amplification history. We observed that the contents of transposable elements, long terminal repeats (LTR) retrotransposons, and particularly LTR/Gypsy retrotransposons varied largely but were significantly correlated with genome sizes. Therefore, this study demonstrated that LTR retrotransposons act as an active driver of genome size variation in the genus Oryza.

A phylogenetic analysis of Momordica (Cucurbitaceae) in India based on karyo-morphology, nuclear DNA content and rDNA ITS1-5.8S-ITS2 sequences

The infrageneric delimitation of Momordica, a medicinally important genus of Cucurbitaceae, is ill-defined until date. Momordica chromosomes are extremely small and are difficult to stain and visualize because of the dense cytoplasmic background. We have conducted karyomorphometric analysis by EMA method in five Indian Momordica species, and the nuclear genome sizes were estimated by flow cytometry for the first time. The somatic chromosome numbers ranged from 2n = 18 to 56 in the species. We have resolved previously disputed chromosome numbers in M. cymbalaria and M. dioica as 2n = 18 (lowest) and 2n = 56, respectively. Chromosome counts in the other species were re-confirmed as 2n = 22 in M. charantia, 2n = 28 in M. cochinchinensis and 2n = 56 in M. subangulata. The largest genome size was recorded in M. cymbalaria (3.74 pg 2C^-1), while the smallest size (0.72 pg 2C^-1) was detected in M. charantia var. charantia. The nuclear genome sizes were analysed in comparison to chromosome numbers and total chromosome lengths of the species. Karyomorphometric indices showed comparable symmetric karyotypes in the species except in M. cymbalaria having tendency towards asymmetry. The UPGMA phenogram and principle component analysis based on nuclear DNA contents and karyomorphometric parameters demonstrated interspecies differences, intraspecific distinction within M. charantia varieties and highlighted distinction of M. cymbalaria. This study was further supported by the rDNA ITS sequence-based phylogenetic analysis which revealed the monophyletic origin of the Indian members of Momordica and clarified the intraspecies relationship among the studied members. As a whole, the study brought out new insights on species diversification within the genus Momordica in India and would benefit further studies on biosystematics and plant breeding programmes.

Insights into triterpene synthesis and unsaturated fatty-acid accumulation provided by chromosomal-level genome analysis of Akebia trifoliata subsp. australis

Akebia trifoliata subsp. australis is a well-known medicinal and potential woody oil plant in China. The limited genetic information available for A. trifoliata subsp. australis has hindered its exploitation. Here, a high-quality chromosome-level genome sequence of A. trifoliata subsp. australis is reported. The de novo genome assembly of 682.14 Mb was generated with a scaffold N50 of 43.11 Mb. The genome includes 25,598 protein-coding genes, and 71.18% (485.55 Mb) of the assembled sequences were identified as repetitive sequences. An ongoing massive burst of long terminal repeat (LTR) insertions, which occurred ~1.0 million years ago, has contributed a large proportion of LTRs in the genome of A. trifoliata subsp. australis. Phylogenetic analysis shows that A. trifoliata subsp. australis is closely related to Aquilegia coerulea and forms a clade with Papaver somniferum and Nelumbo nucifera, which supports the well-established hypothesis of a close relationship between basal eudicot species. The expansion of UDP-glucoronosyl and UDP-glucosyl transferase gene families and β-amyrin synthase-like genes and the exclusive contraction of terpene synthase gene families may be responsible for the abundant oleanane-type triterpenoids in A. trifoliata subsp. australis. Furthermore, the acyl-ACP desaturase gene family, including 12 stearoyl-acyl-carrier protein desaturase (SAD) genes, has expanded exclusively. A combined transcriptome and fatty-acid analysis of seeds at five developmental stages revealed that homologs of SADs, acyl-lipid desaturase omega fatty acid desaturases (FADs), and oleosins were highly expressed, consistent with the rapid increase in the content of fatty acids, especially unsaturated fatty acids. The genomic sequences of A. trifoliata subsp. australis will be a valuable resource for comparative genomic analyses and molecular breeding.

Concerted and birth-and-death evolution of 26S ribosomal DNA in Camellia L

BACKGROUND AND AIMS

The ribosomal DNA (rDNA) gene family, encoding ribosomal RNA (rRNA), has long been regarded as an archetypal example illustrating the model of concerted evolution. However, controversy is arising, as rDNA in many eukaryotic species has been proved to be polymorphic. Here, a metagenomic strategy was applied to detect the intragenomic polymorphism as well as the evolutionary patterns of 26S rDNA across the genus Camellia.

METHODS

Degenerate primer pairs were designed to amplify the 26S rDNA fragments from different Camellia species. The amplicons were then paired-end sequenced on the Illumina MiSeq platform.

KEY RESULTS

An extremely high level of rDNA polymorphism existed universally in Camellia. However, functional rDNA was still the major component of the family, and was relatively conserved among different Camellia species. Sequence variations mainly came from rRNA pseudogenes and favoured regions that are rich in GC. Specifically, some rRNA pseudogenes have existed in the genome for a long time, and have even experienced several expansion events, which has greatly enriched the abundance of rDNA polymorphism.

CONCLUSIONS

Camellia represents a group in which rDNA is subjected to a mixture of concerted and birth-and-death evolution. Some rRNA pseudogenes may still have potential functions. Conversely, when released from selection constraint, they can evolve in the direction of decreasing GC content and structural stability through a methylation-induced process, and finally be eliminated from the genome.

SMRT- and Illumina-based RNA-seq analyses unveil the ginsinoside biosynthesis and transcriptomic complexity in Panax notoginseng

Panax notoginseng is one of the most widely used traditional Chinese herbs with particularly valued roots. Triterpenoid saponins are mainly specialized secondary metabolites, which medically act as bioactive components. Knowledge of the ginsenoside biosynthesis in P. notoginseng, which is of great importance in the industrial biosynthesis and genetic breeding program, remains largely undetermined. Here we combined single molecular real time (SMRT) and Second-Generation Sequencing (SGS) technologies to generate a widespread transcriptome atlas of P. notoginseng. We mapped 2,383 full-length non-chimeric (FLNC) reads to adjacently annotated genes, corrected 1,925 mis-annotated genes and merged into 927 new genes. We identified 8,111 novel transcript isoforms that have improved the annotation of the current genome assembly, of which we found 2,664 novel lncRNAs. We characterized more alternative splicing (AS) events from SMRT reads (20,015 AS in 6,324 genes) than Illumina reads (18,498 AS in 9,550 genes), which contained a number of AS events associated with the ginsenoside biosynthesis. The comprehensive transcriptome landscape reveals that the ginsenoside biosynthesis predominantly occurs in flowers compared to leaves and roots, substantiated by levels of gene expression, which is supported by tissue-specific abundance of isoforms in flowers compared to roots and rhizomes. Comparative metabolic analyses further show that a total of 17 characteristic ginsenosides increasingly accumulated, and roots contained the most ginsenosides with variable contents, which are extraordinarily abundant in roots of the three-year old plants. We observed that roots were rich in protopanaxatriol- and protopanaxadiol-type saponins, whereas protopanaxadiol-type saponins predominated in aerial parts (leaves, stems and flowers). The obtained results will greatly enhance our understanding about the ginsenoside biosynthetic machinery in the genus Panax.

The Reference Genome of Tea Plant and Resequencing of 81 Diverse Accessions Provide Insights into Its Genome Evolution and Adaptation

Tea plant is an important economic crop, which is used to produce the world's oldest and most widely consumed tea beverages. Here, we present a high-quality reference genome assembly of the tea plant (Camellia sinensis var. sinensis) consisting of 15 pseudo-chromosomes. LTR retrotransposons (LTR-RTs) account for 70.38% of the genome, and we present evidence that LTR-RTs play critical roles in genome size expansion and the transcriptional diversification of tea plant genes through preferential insertion in promoter regions and introns. Genes, particularly those coding for terpene biosynthesis proteins, associated with tea aroma and stress resistance were significantly amplified through recent tandem duplications and exist as gene clusters in tea plant genome. Phylogenetic analysis of the sequences of 81 tea plant accessions with diverse origins revealed three well-differentiated tea plant populations, supporting the proposition for the southwest origin of the Chinese cultivated tea plant and its later spread to western Asia through introduction. Domestication and modern breeding left significant signatures on hundreds of genes in the tea plant genome, particularly those associated with tea quality and stress resistance. The genomic sequences of the reported reference and resequenced tea plant accessions provide valuable resources for future functional genomics study and molecular breeding of improved cultivars of tea plants.

The draft genome sequence of an upland wild rice species, Oryza granulata

Exploiting novel gene sources from wild relatives has proven to be an efficient approach to advance crop genetic breeding efforts. Oryza granulata, with the GG genome type, occupies the basal position of the Oryza phylogeny and has the second largest genome (~882 Mb). As an upland wild rice species, it possesses renowned traits that distinguish it from other Oryza species, such as tolerance to shade and drought, immunity to bacterial blight and resistance to the brown planthopper. Here, we generated a 736.66-Mb genome assembly of O. granulata with 40,131 predicted protein-coding genes. With Hi-C data, for the first time, we anchored ~98.2% of the genome assembly to the twelve pseudo-chromosomes. This chromosome-length genome assembly of O. granulata will provide novel insights into rice genome evolution, enhance our efforts to search for new genes for future rice breeding programmes and facilitate the conservation of germplasm of this endangered wild rice species.

Improved hybrid de novo genome assembly and annotation of African wild rice, Oryza longistaminata, from Illumina and PacBio sequencing reads

African wild rice Oryza longistaminata, one of the eight AA- genome species in the genus Oryza, possesses highly valued traits, such as the rhizomatousness for perennial rice breeding, strong tolerance to biotic and abiotic stresses, and high biomass production on poor soils. To obtain the high-quality reference genome for O. longistaminata we employed a hybrid assembly approach through incorporating Illumina and PacBio sequencing datasets. The final genome assembly comprised only 107 scaffolds and was approximately ∼363.5 Mb, representing ∼92.7% of the estimated African wild rice genome (∼392 Mb). The N50 lengths of the assembled contigs and scaffolds were ∼46.49 Kb and ∼6.83 Mb, indicating ∼3.72-fold and ∼18.8-fold improvement in length compared to the earlier released assembly (∼12.5 Kb and 364 Kb, respectively). Aided with Hi-C data and syntenic relationship with O. sativa, these assembled scaffolds were anchored into 12 pseudo-chromosomes. Genome annotation and comparative genomic analysis reveal that lineage-specific expansion of gene families that respond to biotic- and abiotic stresses are of great potential for mining novel alleles to overcome major diseases and abiotic adaptation in rice breeding programs. This reference genome of African wild rice will greatly enlarge the existing database of rice genome resources and unquestionably form a solid base to understand genomic basis underlying highly valued phenotypic traits and search for novel gene sources in O. longistaminata for the future rice breeding programs.

Rooibos (Aspalathus linearis) Genome Size Estimation Using Flow Cytometry and K-Mer Analyses

Plant genomes provide information on biosynthetic pathways involved in the production of industrially relevant compounds. Genome size estimates are essential for the initiation of genome projects. The genome size of rooibos (Aspalathus linearis species complex) was estimated using DAPI flow cytometry and k-mer analyses. For flow cytometry, a suitable nuclei isolation buffer, plant tissue and a transport medium for rooibos ecotype samples collected from distant locations were identified. When using radicles from commercial rooibos seedlings, Woody Plant Buffer and Vicia faba as an internal standard, the flow cytometry-estimated genome size of rooibos was 1.24 ± 0.01 Gbp. The estimates for eight wild rooibos growth types did not deviate significantly from this value. K-mer analysis was performed using Illumina paired-end sequencing data from one commercial rooibos genotype. For biocomputational estimation of the genome size, four k-mer analysis methods were investigated: A standard formula and three popular programs (BBNorm, GenomeScope, and FindGSE). GenomeScope estimates were strongly affected by parameter settings, specifically CovMax. When using the complete k-mer frequency histogram (up to 9 × 105), the programs did not deviate significantly, estimating an average rooibos genome size of 1.03 ± 0.04 Gbp. Differences between the flow cytometry and biocomputational estimates are discussed.

A draft genome of a field-collected Steinernema feltiae strain NW

Advances in sequencing technologies have accelerated our understanding of the complex genetic network of organisms and genomic divergences that are linked to evolutionary processes. While many model organisms and laboratory strains have been sequenced, wild populations are underrepresented in the growing list of sequenced genomes. Here, we present a de novo assembly of Steinernema feltiae, strain NW, collected from a working agricultural field in south central Washington, USA. Leveraging Pacific Biosciences (PacBio) long reads, we sequenced strain NW to a high depth (99×). The resulting de novo assembly is significantly larger than the previous assembly generated from the laboratory strain SN, with a noticeable improvement in continuity and completeness. Comparative analysis of two assemblies revealed numerous single nucleotide polymorphisms (SNPs), breakpoints, and indels present between the two genomes. This alternative genome resource and annotation could benefit the research community to examine the genetic foundation of evolutionary processes as well as genomic variation among conspecific populations.

Artificial Chromosome Doubling in Allotetraploid Calendula officinalis

Calendula officinalis L. is known as an ornamental plant as well as a source of biochemical compounds used in cosmetics and industry. C. officinalis has a complex karyotype. Published chromosome numbers differ between 2n = 4x = 28 or 32. We have estimated genome sizes in nine commercial cultivars and evaluated the ploidy level by karyotyping and fluorescent in situ hybridization (FISH) using 5S and 45S rDNA loci. The detection of chromosome sets of two rather than four homologues would suggest that C. officinalis has an allotetraploid background. In addition, four signals for 45S but only two for 5S were found by using FISH. Artificial chromosome doubling is a common technique in plant breeding, as polyploidization results in several consequences for plant growth and development. Especially the suggested allotetraploid background in C. officinalis is interesting when examining the effect of chromosome doubling on the plant phenotype. Here we describe chromosome doubling of three allotetraploid cultivars of C. officinalis, 'Nova,' 'WUR 1553-7' and 'Orange Beauty'. Three antimitotic agents - colchicine, oryzalin and trifluralin - were used in different concentrations to find the combination of the best agent and the best dosage to obtain octaploids. For all three cultivars a few octaploids were obtained. A concentration of 200 and 400 ppm of colchicine was most efficient for chromosome doubling in 'Nova' and 'Orange Beauty,' respectively. For 'WUR 1553-7' the treatment with 20 ppm oryzalin was also effective. Cell numbers and first observations of the phenotype in the chromosome doubled plants show thicker leaves and bigger cells, as commonly observed after ploidy doubling. Due to the low number of chromosome doubled plants obtained more elaborate phenotyping will be performed on following generations cultivated under field conditions.

Tea plant genomics: achievements, challenges and perspectives

Tea is among the world's most widely consumed non-alcoholic beverages and possesses enormous economic, health, and cultural values. It is produced from the cured leaves of tea plants, which are important evergreen crops globally cultivated in over 50 countries. Along with recent innovations and advances in biotechnologies, great progress in tea plant genomics and genetics has been achieved, which has facilitated our understanding of the molecular mechanisms of tea quality and the evolution of the tea plant genome. In this review, we briefly summarize the achievements of the past two decades, which primarily include diverse genome and transcriptome sequencing projects, gene discovery and regulation studies, investigation of the epigenetics and noncoding RNAs, origin and domestication, phylogenetics and germplasm utilization of tea plant as well as newly developed tools/platforms. We also present perspectives and possible challenges for future functional genomic studies that will contribute to the acceleration of breeding programs in tea plants.

Estimation of nuclear DNA content and its variation among Indian Tea accessions by flow cytometry

Nuclear DNA content and genome size variation among 36 Indian tea accessions were analyzed by flow cytometry. Initial standardization of protocols for isolation of nuclei, DNA staining and selection of an internal standard for tea accessions which have significantly high amount of phenolic secondary metabolites in their cytosol was carried out. Results obtained revealed that 2C DNA content of Indian tea is 7.46 pg which corresponds to 1C genome size of 3673 Mb. Inter accession variation in 2C DNA content was also observed among 35 diploid taxa ranging from 7.23 to 7.73 pg which was significant at 1% probability level. The 2C DNA content of triploid (UPASI 3) was observed to be 11.47 pg which is concurrent with the expected value. Results obtained showed that Assam and Cambod type tea accession have higher 2C DNA content of 7.73 pg whereas Assam Cambod hybrids and Assam China hybrids have reduction in DNA content with 2C amounts, 7.23 and 7.32 pg DNA respectively. The present study suggests that the species involved in origin of Indian tea must have differed in their genome sizes owing to significant inter accession variation in nuclear DNA content.

Circular RNA architecture and differentiation during leaf bud to young leaf development in tea (Camellia sinensis)

Circular RNA (circRNA) discovery, expression patterns and experimental validation in developing tea leaves indicates its correlation with circRNA-parental genes and potential roles in ceRNA interaction network. Circular RNAs (circRNAs) have recently emerged as a novel class of abundant endogenous stable RNAs produced by circularization with regulatory potential. However, identification of circRNAs in plants, especially in non-model plants with large genomes, is challenging. In this study, we undertook a systematic identification of circRNAs from different stage tissues of tea plant (Camellia sinensis) leaf development using rRNA-depleted circular RNA-seq. By combining two state-of-the-art detecting tools, we characterized 3174 circRNAs, of which 342 were shared by each approach, and thus considered high-confidence circRNAs. A few predicted circRNAs were randomly chosen, and 20 out of 24 were experimental confirmed by PCR and Sanger sequencing. Similar in other plants, tissue-specific expression was also observed for many C. sinensis circRNAs. In addition, we found that circRNA abundances were positively correlated with the mRNA transcript abundances of their parental genes. qRT-PCR validated the differential expression patterns of circRNAs between leaf bud and young leaf, which also indicated the low expression abundance of circRNAs compared to the standard mRNAs from the parental genes. We predicted the circRNA-microRNA interaction networks, and 54 of the differentially expressed circRNAs were found to have potential tea plant miRNA binding sites. The gene sets encoding circRNAs were significantly enriched in chloroplasts related GO terms and photosynthesis/metabolites biosynthesis related KEGG pathways, suggesting the candidate roles of circRNAs in photosynthetic machinery and metabolites biosynthesis during leaf development.

Genome survey sequencing of red swamp crayfish Procambarus clarkii

Red swamp crayfish, Procambarus clarkii, presently is an important aquatic commercial species in China. The crayfish is a hot area of research focus, and its genetic improvement is quite urgent for the crayfish aquaculture in China. However, the knowledge of its genomic landscape is limited. In this study, a survey of P. clarkii genome was investigated based on Illumina's Solexa sequencing platform. Meanwhile, its genome size was estimated using flow cytometry. Interestingly, the genome size estimated is about 8.50 Gb by flow cytometry and 1.86 Gb with genome survey sequencing. Based on the assembled genome sequences, total of 136,962 genes and 152,268 exons were predicted, and the predicted genes ranged from 150 to 12,807 bp in length. The survey sequences could help accelerate the progress of gene discovery involved in genetic diversity and evolutionary analysis, even though it could not successfully applied for estimation of P. clarkii genome size.

Comparative genomics analysis reveals gene family expansion and changes of expression patterns associated with natural adaptations of flowering time and secondary metabolism in yellow Camellia

Yellow-flowering species are unique in the genus Camellia not only for their bright yellow pigments but also the health-improving substances in petals. However, little is known regarding the biosynthesis pathways of pigments and secondary metabolites. Here, we performed comparative genomics studies in two yellow-flowered species of the genus Camellia with distinctive flowering periods. We obtained 112,190 and 89,609 unigenes from Camellia nitidissima and Camellia chuongtsoensis, respectively, and identified 9547 gene family clusters shared with various plant species and 3414 single-copy gene families. Global gene expression analysis revealed six comparisons of differentially expressed gene sets in different developmental stages of floral bud. Through the identification of orthologous pairs, conserved and specific differentially expressed genes (DEGs) between species were compared. Functional enrichment analysis suggested that the gibberellin (GA) biosynthesis pathway might be related to the alteration of flowering responses. Furthermore, the expression patterns of secondary metabolism pathway genes were analyzed between yellow- and red-flowered Camellias. We showed that the key enzymes involved in glycosylation of flavonoids displayed differential expression patterns, indicating that the direct glycosylation of flavonols rather than anthocyanins was pivotal to coloration and health-improving metabolites in the yellow Camellia petals. Finally, the gene family analysis of UDP-glycosyltransferases revealed an expansion of group C members in C. nitidissima. Through comparative genomics analysis, we demonstrate that changes of gene expression and gene family members are critical to the variation of natural traits. This work provides valuable insights into the molecular regulation of trait adaptations of floral pigmentation and flowering timing.

Intraspecific variations in cardamom (Elettaria cardamomum Maton): assessment of genomic diversity by flow cytometry, cytological studies and ISSR analysis

BACKGROUND

The main goal of the work was to analyse intraspecific variation in Elettaria cardamomum Maton (cardamom) using genome size, cytological studies and molecular marker data. Nuclear DNA content and molecular marker details furnish data on genome size and genetic diversity respectively among the studied accessions and both complement each other for evolutionary and taxonomic studies.

RESULTS

The relative 2C genome size and total number of base pairs of cardamom was determined through flow cytometric analysis using propidium iodide staining. The nuclear DNA content was estimated in various sections of the species representing individuals from wild and cultivar genotypes following Zea mays L. CE-777 (2C = 5.43 pg) as internal reference standard. Chromosome number from growing root tip was examined following standard protocols. Twenty-six ISSR primers that generated polymorphic bands were used for genetic diversity analysis of the thirty accessions of cardamom. Estimated nuclear 2C DNA content ranged from 2.57 to 3.22 pg demonstrating 1.25-fold variation. The mean amount of 2C nuclear DNA of the cardamom was calculated as 2.87 pg which is equivalent of 2806 Mbp as the diploid genome size. The chromosome number was found to be 2n = 48. Among the thirty accessions of cardamom studied using ISSR markers, C53 (feral from Bonacaud) showed a very prominent level of genetic diversity and was lowest for C96 (Avinash-I, a released variety from Indian Institute of Spices Research, Kozhikode).

CONCLUSION

These analyses revealed the existence of genetic variability within the studied cardamom accessions. The plant specimens also differed significantly in their genome size. However, the genetic variability parameters did not show any correlation with genome size.

Testing the link between genome size and growth rate in maize

Little is known about the factors driving within species Genome Size (GS) variation. GS may be shaped indirectly by natural selection on development and adaptative traits. Because GS variation is particularly pronounced in maize, we have sampled 83 maize inbred lines from three well described genetic groups adapted to contrasted climate conditions: inbreds of tropical origin, Flint inbreds grown in temperate climates, and Dent inbreds distributed in the Corn Belt. As a proxy for growth rate, we measured the Leaf Elongation Rate maximum during nighttime (LER_max) as well as GS in all inbred lines. In addition we combined available and new nucleotide polymorphism data at 29,090 sites to characterize the genetic structure of our panel. We found significant variation for both LER_max and GS among groups defined by our genetic structuring. Tropicals displayed larger GS than Flints while Dents exhibited intermediate values. LER_max followed the opposite trend with greater growth rate in Flints than in Tropicals. In other words, LER_max and GS exhibited a significantly negative correlation (r = − 0.27). However, this correlation was driven by among-group variation rather than within-group variation—it was no longer significant after controlling for structure and kinship among inbreds. Our results indicate that selection on GS may have accompanied ancient maize diffusion from its center of origin, with large DNA content excluded from temperate areas. Whether GS has been targeted by more intense selection during modern breeding within groups remains an open question.

Genetic Divergence between Camellia sinensis and Its Wild Relatives Revealed via Genome-Wide SNPs from RAD Sequencing

Tea is one of the most popular beverages across the world and is made exclusively from cultivars of Camellia sinensis. Many wild relatives of the genus Camellia that are closely related to C. sinensis are native to Southwest China. In this study, we first identified the distinct genetic divergence between C. sinensis and its wild relatives and provided a glimpse into the artificial selection of tea plants at a genome-wide level by analyzing 15,444 genomic SNPs that were identified from 18 cultivated and wild tea accessions using a high-throughput genome-wide restriction site-associated DNA sequencing (RAD-Seq) approach. Six distinct clusters were detected by phylogeny inferrence and principal component and genetic structural analyses, and these clusters corresponded to six Camellia species/varieties. Genetic divergence apparently indicated that C. taliensis var. bangwei is a semi-wild or transient landrace occupying a phylogenetic position between those wild and cultivated tea plants. Cultivated accessions exhibited greater heterozygosity than wild accessions, with the exception of C. taliensis var. bangwei. Thirteen genes with non-synonymous SNPs exhibited strong selective signals that were suggestive of putative artificial selective footprints for tea plants during domestication. The genome-wide SNPs provide a fundamental data resource for assessing genetic relationships, characterizing complex traits, comparing heterozygosity and analyzing putatitve artificial selection in tea plants.

Transcriptome Analysis Identifies Candidate Genes Related to Triacylglycerol and Pigment Biosynthesis and Photoperiodic Flowering in the Ornamental and Oil-Producing Plant, Camellia reticulata (Theaceae)

Camellia reticulata, which is native to Southwest China, is famous for its ornamental flowers and high-quality seed oil. However, the lack of genomic information for this species has largely hampered our understanding of its key pathways related to oil production, photoperiodic flowering process and pigment biosynthesis. Here, we first sequenced and characterized the transcriptome of a diploid C. reticulata in an attempt to identify genes potentially involved in triacylglycerol biosynthesis (TAGBS), photoperiodic flowering, flavonoid biosynthesis (FlaBS), carotenoid biosynthesis (CrtBS) pathways. De novo assembly of the transcriptome provided a catalog of 141,460 unigenes with a total length of ~96.1 million nucleotides (Mnt) and an N50 of 1080 nt. Of them, 22,229 unigenes were defined as differentially expressed genes (DEGs) across five sequenced tissues. A large number of annotated genes in C. reticulata were found to have been duplicated, and differential expression patterns of these duplicated genes were commonly observed across tissues, such as the differential expression of SOC1_a, SOC1_b, and SOC1_c in the photoperiodic flowering pathway. Up-regulation of SAD_a and FATA genes and down-regulation of FAD2_a gene in the TAGBS pathway in seeds may be relevant to the ratio of monounsaturated fatty acid (MUFAs) to polyunsaturated fatty acid (PUFAs) in seed oil. MYBF1, a transcription regulator gene of the FlaBS pathway, was found with great sequence variation and alteration of expression patterns, probably resulting in functionally evolutionary differentiation in C. reticulata. MYBA1_a and some anthocyanin-specific biosynthetic genes in the FlaBS pathway were highly expressed in both flower buds and flowers, suggesting important roles of anthocyanin biosynthesis in flower development. Besides, a total of 40,823 expressed sequence tag simple sequence repeats (EST-SSRs) were identified in the C. reticulata transcriptome, providing valuable marker resources for further basic and applied researches on this economically important Camellia plant.

Quantitative testing of the methodology for genome size estimation in plants using flow cytometry: a case study of the Primulina genus

Flow cytometry (FCM) is a commonly used method for estimating genome size in many organisms. The use of FCM in plants is influenced by endogenous fluorescence inhibitors and may cause an inaccurate estimation of genome size; thus, falsifying the relationship between genome size and phenotypic traits/ecological performance. Quantitative optimization of FCM methodology minimizes such errors, yet there are few studies detailing this methodology. We selected the genus Primulina, one of the most representative and diverse genera of the Old World Gesneriaceae, to evaluate the methodology effect on determining genome size. Our results showed that buffer choice significantly affected genome size estimation in six out of the eight species examined and altered the 2C-value (DNA content) by as much as 21.4%. The staining duration and propidium iodide (PI) concentration slightly affected the 2C-value. Our experiments showed better histogram quality when the samples were stained for 40 min at a PI concentration of 100 μg ml(-1). The quality of the estimates was not improved by 1-day incubation in the dark at 4°C or by centrifugation. Thus, our study determined an optimum protocol for genome size measurement in Primulina: LB01 buffer supplemented with 100 μg ml(-1) PI and stained for 40 min. This protocol also demonstrated a high universality in other Gesneriaceae genera. We report the genome size of nine Gesneriaceae species for the first time. The results showed substantial genome size variation both within and among the species, with the 2C-value ranging between 1.62 and 2.71 pg. Our study highlights the necessity of optimizing the FCM methodology prior to obtaining reliable genome size estimates in a given taxon.

De novo transcriptome assembly of the wild relative of tea tree (Camellia taliensis) and comparative analysis with tea transcriptome identified putative genes associated with tea quality and stress response

Background

Camellia taliensis is one of the most important wild relatives of cultivated tea tree, C. sinensis. The species extensively occupies mountainous habitats representing a wide-range abiotic tolerance and biotic resistance and thus harbors valuable gene resources that may greatly benefit genetic improvement of cultivated tea tree. However, owning to a large genome size of ~3 Gb and structurally complex genome, there are fairly limited genetic information and particularly few genomic resources publicly available for this species. To better understand the key pathways determining tea flavor and enhance tea tree breeding programs, we performed a high-throughput transcriptome sequencing for C. taliensis.

Results

In this study, approximate 241.5 million high-quality paired-end reads, accounting for ~24 Gb of sequence data, were generated from tender shoots, young leaves, flower buds and flowers using Illumina HiSeq 2000 platform. De novo assembly with further processing and filtering yielded a set of 67,923 transcripts with an average length of 685 bp and an N50 of 995 bp. Based on sequence similarity searches against public databases, a total of 39,475 transcripts were annotated with gene descriptions, conserved protein domains or gene ontology (GO) terms. Candidate genes for major metabolic pathways involved in tea quality were identified and experimentally validated using RT-qPCR. Further gene expression profiles showed that they are differentially regulated at different developmental stages. To gain insights into the evolution of these genes, we aligned them to the previously cloned orthologous genes in C. sinensis, and found that considerable nucleotide variation within several genes involved in important secondary metabolic biosynthesis pathways, of which flavone synthase II gene (FNSII) is the most variable between these two species. Moreover, comparative analyses revealed that C. taliensis shows a remarkable expansion of LEA genes, compared to C. sinensis, which might contribute to the observed stronger stress resistance of C. taliensis.

Conclusion

We reported the first large-coverage transcriptome datasets for C. taliensis using the next-generation sequencing technology. Such comprehensive EST datasets provide an unprecedented opportunity for identifying genes involved in several major metabolic pathways and will accelerate functional genomic studies and genetic improvement efforts of tea trees in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1494-4) contains supplementary material, which is available to authorized users.

Thirteen Camellia chloroplast genome sequences determined by high-throughput sequencing: genome structure and phylogenetic relationships

BACKGROUND

Camellia is an economically and phylogenetically important genus in the family Theaceae. Owing to numerous hybridization and polyploidization, it is taxonomically and phylogenetically ranked as one of the most challengingly difficult taxa in plants. Sequence comparisons of chloroplast (cp) genomes are of great interest to provide a robust evidence for taxonomic studies, species identification and understanding mechanisms that underlie the evolution of the Camellia species.

RESULTS

The eight complete cp genomes and five draft cp genome sequences of Camellia species were determined using Illumina sequencing technology via a combined strategy of de novo and reference-guided assembly. The Camellia cp genomes exhibited typical circular structure that was rather conserved in genomic structure and the synteny of gene order. Differences of repeat sequences, simple sequence repeats, indels and substitutions were further examined among five complete cp genomes, representing a wide phylogenetic diversity in the genus. A total of fifteen molecular markers were identified with more than 1.5% sequence divergence that may be useful for further phylogenetic analysis and species identification of Camellia. Our results showed that, rather than functional constrains, it is the regional constraints that strongly affect sequence evolution of the cp genomes. In a substantial improvement over prior studies, evolutionary relationships of the section Thea were determined on basis of phylogenomic analyses of cp genome sequences.

CONCLUSIONS

Despite a high degree of conservation between the Camellia cp genomes, sequence variation among species could still be detected, representing a wide phylogenetic diversity in the genus. Furthermore, phylogenomic analysis was conducted using 18 complete cp genomes and 5 draft cp genome sequences of Camellia species. Our results support Chang's taxonomical treatment that C. pubicosta may be classified into sect. Thea, and indicate that taxonomical value of the number of ovaries should be reconsidered when classifying the Camellia species. The availability of these cp genomes provides valuable genetic information for accurately identifying species, clarifying taxonomy and reconstructing the phylogeny of the genus Camellia.