Telomere-to-telomere genome assembly of Eleocharis dulcis and expression profiles during corm development.

Eleocharis dulcis (Burm. f.) Trin. ex Hensch., commonly known as Chinese water chestnut, is a traditional aquatic vegetable in China, and now is widely cultivated throughout the world because of its high nutritional value and unique tastes. Here, we report the assembly of a 493.24 Mb telomere-to-telomere (T2T) genome of E. dulcis accomplished by integrating ONT ultra-long reads, PacBio long reads and Hi-C data. The reference genome was anchored onto 111 gap-free chromosomes, containing 48.31% repeat elements and 33,493 predicted protein-coding genes. Whole genome duplication (WGD) and inter-genomic synteny analyses indicated that chromosome breakage and genome duplication in E. dulcis possibly occurred multiple times during genome evolution after its divergence from a common ancestor with Rhynchospora breviuscula at ca. 35.6 Mya. A comparative time-course transcriptome analysis of corm development revealed different patterns of gene expression between cultivated and wild accessions with the highest number of differentially expressed genes (DEGs, 15,870) at the middle swelling stage and some of the DEGs were significantly enriched for starch metabolic process.

Chromosome-level assemblies of cultivated water chestnut Trapa bicornis and its wild relative Trapa incisa.

Water chestnut (Trapa L.) is a floating-leaved aquatic plant with high edible and medicinal value. In this study, we presented chromosome-level genome assemblies of cultivated large-seed species Trapa bicornis and its wild small-seed relative Trapa incisa by using PacBio HiFi long reads and Hi-C technology. The T. bicornis and T. incisa assemblies consisted of 479.90 Mb and 463.97 Mb contigs with N50 values of 13.52 Mb and 13.77 Mb, respectively, and repeat contents of 62.88% and 62.49%, respectively. A total of 33,306 and 33,315 protein-coding genes were predicted in T. bicornis and T. incisa assemblies, respectively. There were 159,232 structural variants affecting more than 11 thousand genes detected between the two genomes. The phylogenetic analysis indicated that the lineage leading to Trapa was diverged from the lineage to Sonneratia approximately 23 million years ago. These two assemblies provide valuable resources for future evolutionary and functional genomic research and molecular breeding of water chestnut.

Unprecedented Intraindividual Structural Heteroplasmy in Eleocharis (Cyperaceae, Poales) Plastomes.

Plastid genomes (plastomes) of land plants have a conserved quadripartite structure in a gene-dense unit genome consisting of a large inverted repeat that separates two single copy regions. Recently, alternative plastome structures were suggested in Geraniaceae and in some conifers and Medicago the coexistence of inversion isomers has been noted. In this study, plastome sequences of two Cyperaceae, Eleocharis dulcis (water chestnut) and Eleocharis cellulosa (gulf coast spikerush), were completed. Unlike the conserved plastomes in basal groups of Poales, these Eleocharis plastomes have remarkably divergent features, including large plastome sizes, high rates of sequence rearrangements, low GC content and gene density, gene duplications and losses, and increased repetitive DNA sequences. A novel finding among these features was the unprecedented level of heteroplasmy with the presence of multiple plastome structural types within a single individual. Illumina paired-end assemblies combined with PacBio single-molecule real-time sequencing, long-range polymerase chain reaction, and Sanger sequencing data identified at least four different plastome structural types in both Eleocharis species. PacBio long read data suggested that one of the four E. dulcis plastome types predominates.

Identification and characterization of water chestnut Soymovirus-1 (WCSV-1), a novel Soymovirus in water chestnuts (Eleocharis dulcis).

BACKGROUND: A disease of unknown etiology in water chestnut plants (Eleocharis dulcis) was reported in China between 2012 and 2014. High throughput sequencing of small RNA (sRNA) combined with bioinformatics, and molecular identification based on PCR detection with virus-specific primers and DNA sequencing is a desirable approach to identify an unknown infectious agent. In this study, we employed this approach to identify viral sequences in water chestnut plants and to explore the molecular interaction of the identified viral pathogen and its natural plant host. RESULTS: Based on high throughput sequencing of virus-derived small RNAs (vsRNA), we identified the sequence a new-to-science double-strand DNA virus isolated from water chestnut cv. 'Tuanfeng' samples, a widely grown cultivar in Hubei province, China, and analyzed its genomic organization. The complete genomic sequence is 7535 base-pairs in length, and shares 42-52% nucleotide sequence identity with viruses in the Caulimoviridae family. The virus contains nine predicated open reading frames (ORFs) encoding nine hypothetical proteins, with conserved domains characteristic of caulimoviruses. Phylogenetic analyses at the nucleotide and amino acid levels indicated that the virus belongs to the genus Soymovirus. The virus is tentatively named Water chestnut soymovirus-1 (WCSV-1). Phylogenetic analysis of the putative viral polymerase protein suggested that WCSV-1 is distinct to other well established species in the Soymovirus genus. This conclusion was supported by phylogenetic analyses of the amino acid sequences encoded by ORFs I, IV, VI, or VII. The sRNA bioinformatics showed that the majority of the vsRNAs are 22-nt in length with a preference for U at the 5'-terminal nucleotide. The vsRNAs are unevenly distributed over both strands of the entire WCSV-1 circular genome, and are clustered into small defined regions. In addition, we detected WCSV-1 in asymptomatic and symptomatic water chestnut samples collected from different regions of China by using PCR. RNA-seq assays further confirmed the presence of WCSV-1-derived viral RNA in infected plants. CONCLUSIONS: This is the first discovery of a dsDNA virus in the genus Soymovirus infecting water chestnuts. Data presented also add new information towards a better understanding of the co-evolutionary mechanisms between the virus and its natural plant host.

Analysis of retrotransposon abundance, diversity and distribution in holocentric Eleocharis (Cyperaceae) genomes.

Background and Aims: Long terminal repeat-retrotransposons (LTR-RTs) comprise a large portion of plant genomes, with massive repeat blocks distributed across the chromosomes. Eleocharis species have holocentric chromosomes, and show a positive correlation between chromosome numbers and the amount of nuclear DNA. To evaluate the role of LTR-RTs in karyotype diversity in members of Eleocharis (subgenus Eleocharis), the occurrence and location of different members of the Copia and Gypsy superfamilies were compared, covering interspecific variations in ploidy levels (considering chromosome numbers), DNA C-values and chromosomal arrangements. Methods: The DNA C-value was estimated by flow cytometry. Genomes of Eleocharis elegans and E. geniculata were partially sequenced using Illumina MiSeq assemblies, which were a source for searching for conserved proteins of LTR-RTs. POL domains were used for recognition, comparing families and for probe production, considering different families of Copia and Gypsy superfamilies. Probes were obtained by PCR and used in fluorescence in situ hybridization (FISH) against chromosomes of seven Eleocharis species. Key Results: A positive correlation between ploidy levels and the amount of nuclear DNA was observed, but with significant variations between samples with the same ploidy levels, associated with repetitive DNA fractions. LTR-RTs were abundant in E. elegans and E. geniculata genomes, with a predominance of Copia Sirevirus and Gypsy Athila/Tat clades. FISH using LTR-RT probes exhibited scattered and clustered signals, but with differences in the chromosomal locations of Copia and Gypsy. The diversity in LTR-RT locations suggests that there is no typical chromosomal distribution pattern for retrotransposons in holocentric chromosomes, except the CRM family with signals distributed along chromatids. Conclusions: These data indicate independent fates for each LTR-RT family, including accumulation between and within chromosomes and genomes. Differential activity and small changes in LTR-RTs suggest a secondary role in nuclear DNA variation, when compared with ploidy changes.

Transcriptome Analysis of Gene Expression during Chinese Water Chestnut Storage Organ Formation.

The product organ (storage organ; corm) of the Chinese water chestnut has become a very popular food in Asian countries because of its unique nutritional value. Corm formation is a complex biological process, and extensive whole genome analysis of transcripts during corm development has not been carried out. In this study, four corm libraries at different developmental stages were constructed, and gene expression was identified using a high-throughput tag sequencing technique. Approximately 4.9 million tags were sequenced, and 4,371,386, 4,372,602, 4,782,494, and 5,276,540 clean tags, including 119,676, 110,701, 100,089, and 101,239 distinct tags, respectively, were obtained after removal of low-quality tags from each library. More than 39% of the distinct tags were unambiguous and could be mapped to reference genes, while 40% were unambiguous tag-mapped genes. After mapping their functions in existing databases, a total of 11,592, 10,949, 10,585, and 7,111 genes were annotated from the B1, B2, B3, and B4 libraries, respectively. Analysis of the differentially expressed genes (DEGs) in B1/B2, B2/B3, and B3/B4 libraries showed that most of the DEGs at the B1/B2 stages were involved in carbohydrate and hormone metabolism, while the majority of DEGs were involved in energy metabolism and carbohydrate metabolism at the B2/B3 and B3/B4 stages. All of the upregulated transcription factors and 9 important genes related to product organ formation in the above four stages were also identified. The expression changes of nine of the identified DEGs were validated using a quantitative PCR approach. This study provides a comprehensive understanding of gene expression during corm formation in the Chinese water chestnut.

Leaf transcriptome analysis and development of SSR markers in water chestnut (Eleocharis dulcis).

Water chestnut (Eleocharis dulcis) is an important aquatic crop in China; however, transcriptomic and genomic data in public databases are limited. To identify genes and development molecular markers, high-throughput transcriptome sequencing was applied to generate transcript sequences from water chestnut leaf. More than 24 million reads were obtained, trimmed, and assembled into 40,796 contigs with an average length of 616.6 bp. Sequence similarity analyses against 4 public databases (NR, GO, KEGG, KOG) revealed 17,628 contigs that could be annotated with gene descriptions, conserved protein domains, or gene ontology terms. Among the important metabolic pathways, 27 genes were related to starch synthesis and 13 genes were in the steroid synthetic pathway. In addition, 2570 cDNA simple sequence repeats were identified as potential molecular markers in our contigs. One hundred pairs of polymerase chain reaction primers were designed and used for validation of the amplification. The results revealed that 87 primer pairs were successfully amplified in initial screening tests. Overall, this transcriptome dataset and these markers can serve as a platform for further gene expression studies, functional genomic studies, and marker-assisted selection in E. dulcis.

Transcriptomic identification and expression of starch and sucrose metabolism genes in the seeds of Chinese chestnut (Castanea mollissima).

The Chinese chestnut (Castanea mollissima) seed provides a rich source of carbohydrates as food and feed. However, little is known about starch biosynthesis in the seeds. The objectives of this study were to determine seed composition profiles and identify genes involved in starch and sucrose metabolism. Metabolite analysis showed that starch was the major component and rapidly accumulated during seed endosperm development. Amylopectin was approximately 3-fold of amylose content in chestnut starch. Illumina platform-based transcriptome sequencing generated 56671 unigenes in two cDNA libraries from seed endosperms collected at 45 and 75 days after flowering (DAF). A total of 1537 unigenes showed expression differences ≥2-fold in the two stages of seeds including 570 up-regulated and 967 down-regulated unigenes. One hundred and fifty-two unigenes were identified as involved in starch and sucrose metabolism, including 1 for glycogenin glucosyltransferase, 4 for adenylate transporter (brittle1-type), 3 for ADP-glucose pyrophosphorylase (AGP, not brittle2- or shrunken2-type), 3 for starch synthase (SS), 2 for starch branching enzyme, 5 for starch debranching enzyme, 11 for sucrose synthase, and 3 for sucrose-phosphate synthase. Among them, 58 unigenes showed a ≥2-fold expression difference between the 45 and 75 DAF seeds including 11 up- and 47 down-regulated unigenes. The expression of 21 unigenes putatively coding for major enzymes in starch and sucrose metabolism was validated by qPCR using RNA from five seed stages. Expression profiles and correlation analysis indicated that the mRNA levels of AGP (large and small subunits), granule-bound SS2, and soluble SS1 and SS4 were well-correlated with starch accumulation in the seeds. This study suggests that the starch biosynthesis pathway in Chinese chestnut is similar to that of potato tuber/Arabidopsis leaf and differs from that of maize endosperm. The information provides valuable metabolite and genetic resources for future research in starch and sucrose metabolism in Chinese chestnut tree.

Major alterations in transcript profiles between C3-C4 and C4 photosynthesis of an amphibious species Eleocharis baldwinii.

Engineering C4 photosynthetic metabolism into C3 crops is regarded as a major strategy to increase crop productivity, and clarification of the evolutionary processes of C4 photosynthesis can help the better use of this strategy. Here, Eleocharis baldwinii, a species in which C4 photosynthesis can be induced from a C3-C4 state under either environmental or ABA treatments, was used to identify the major transcriptional modifications during the process from C3-C4 to C4. The transcriptomic comparison suggested that in addition to the major differences in C4 core pathway, the pathways of glycolysis, citrate acid metabolism and protein synthesis were dramatically modified during the inducement of C4 photosynthetic states. Transcripts of many transporters, including not only metabolite transporters but also ion transporters, were dramatically increased in C4 photosynthetic state. Many candidate regulatory genes with unidentified functions were differentially expressed in C3-C4 and C4 photosynthetic states. Finally, it was indicated that ABA, auxin signaling and DNA methylation play critical roles in the regulation of C4 photosynthesis. In summary, by studying the different photosynthetic states of the same species, this work provides the major transcriptional differences between C3-C4 and C4 photosynthesis, and many of the transcriptional differences are potentially related to C4 development and therefore are the potential targets for reverse genetics studies.

Identification of C4 photosynthesis metabolism and regulatory-associated genes in Eleocharis vivipara by SSH.

This is the first effort to investigate the candidate genes involved in kranz developmental regulation and C(4) metabolic fluxes in Eleocharis vivipara, which is a leafless freshwater amphibious plant and possesses a distinct culms anatomy structure and photosynthetic pattern in contrasting environments. A terrestrial specific SSH library was constructed to investigate the genes involved in kranz anatomy developmental regulation and C(4) metabolic fluxes. A total of 73 ESTs and 56 unigenes in 384 clones were identified by array hybridization and sequencing. In total, 50 unigenes had homologous genes in the databases of rice and Arabidopsis. The real-time quantitative PCR results showed that most of the genes were accumulated in terrestrial culms and ABA-induced culms. The C(4) marker genes were stably accumulated during the culms development process in terrestrial culms. With respect to C(3) culms, C(4) photosynthesis metabolism consumed much more transporters and translocators related to ion metabolism, organic acids and carbohydrate metabolism, phosphate metabolism, amino acids metabolism, and lipids metabolism. Additionally, ten regulatory genes including five transcription factors, four receptor-like proteins, and one BURP protein were identified. These regulatory genes, which co-accumulated with the culms developmental stages, may play important roles in culms structure developmental regulation, bundle sheath chloroplast maturation, and environmental response. These results shed new light on the C(4) metabolic fluxes, environmental response, and anatomy structure developmental regulation in E. vivipara.

Correlated evolution of LTR retrotransposons and genome size in the genus Eleocharis.

BACKGROUND: Transposable elements (TEs) are considered to be an important source of genome size variation and genetic and phenotypic plasticity in eukaryotes. Most of our knowledge about TEs comes from large genomic projects and studies focused on model organisms. However, TE dynamics among related taxa from natural populations and the role of TEs at the species or supra-species level, where genome size and karyotype evolution are modulated in concert with polyploidy and chromosomal rearrangements, remain poorly understood. We focused on the holokinetic genus Eleocharis (Cyperaceae), which displays large variation in genome size and the occurrence of polyploidy and agmatoploidy/symploidy. We analyzed and quantified the long terminal repeat (LTR) retrotransposons Ty1-copia and Ty3-gypsy in relation to changes in both genome size and karyotype in Eleocharis. We also examined how this relationship is reflected in the phylogeny of Eleocharis. RESULTS: Using flow cytometry, we measured the genome sizes of members of the genus Eleocharis (Cyperaceae). We found positive correlation between the independent phylogenetic contrasts of genome size and chromosome number in Eleocharis. We analyzed PCR-amplified sequences of various reverse transcriptases of the LTR retrotransposons Ty1-copia and Ty3-gypsy (762 sequences in total). Using real-time PCR and dot blot approaches, we quantified the densities of Ty1-copia and Ty3-gypsy within the genomes of the analyzed species. We detected an increasing density of Ty1-copia elements in evolutionarily younger Eleocharis species and found a positive correlation between Ty1-copia densities and C/n-values (an alternative measure of monoploid genome size) in the genus phylogeny. In addition, our analysis of Ty1-copia sequences identified a novel retrotransposon family named Helos1, which is responsible for the increasing density of Ty1-copia. The transition:transversion ratio of Helos1 sequences suggests that Helos1 recently transposed in later-diverging Eleocharis species. CONCLUSIONS: Using several different approaches, we were able to distinguish between the roles of LTR retrotransposons, polyploidy and agmatoploidy/symploidy in shaping Eleocharis genomes and karyotypes. Our results confirm the occurrence of both polyploidy and agmatoploidy/symploidy in Eleocharis. Additionally, we introduce a new player in the process of genome evolution in holokinetic plants: LTR retrotransposons.

Distribution of 45S and 5S rDNA sites in 23 species of Eleocharis (Cyperaceae).

Studies of rDNA location in holocentric chromosomes of the Cyperaceae are scarce, but a few reports have indicated the occurrence of multiple 45S rDNA sites at terminal positions, and in the decondensed state of these regions in prometaphase/metaphase. To extend our knowledge of the number 45S and 5S rDNA sites and distribution in holocentric chromosomes of the Cyperaceae, 23 Brazilian species of Eleocharis were studied. FISH showed 45S rDNA signals always located in terminal regions, which varied from two (E. bonariensis with 2n = 20) to ten (E. flavescens with 2n = 10 and E. laeviglumis with 2n = 60). 5S rDNA showed less variation, with 16 species exhibiting two sites and 7 species four sites, preferentially at terminal positions, except for four species (E. subarticulata, E. flavescens, E. sellowiana and E. geniculata) that showed interstitial sites. The results are discussed in order to understand the predominance of terminal rDNA sites, the mechanisms involved in the interstitial positioning of 5S rDNA sites in some species, and the events of amplification and dispersion of 45S rDNA terminal sites.

Chromosome reduction in Eleocharis maculosa (Cyperaceae).

Chromosome numbers in Cyperaceae lower than the typical basic number x = 5 have been described for only three species: Rhynchospora tenuis (n = 2), Fimbristylis umbellaris (n = 3) and Eleocharis subarticulata (n = 3). Eleocharis maculosa is recorded here as the fourth species of Cyperaceae that has a chromosome number lower than 2n = 10, with 2n = 8, 7 and 6. The karyotype differentiation in E. maculosa was studied using conventional staining (mitosis and meiosis), FISH with 45S and 5S rDNA and telomere probes. The results allow us to determine which chromosomes of the chromosome race with 2n = 10 fused to form the remaining reduced numbers, as well as to understand how the symploidy and translocation mechanisms were important in karyotype differentiation and the formation of chromosome races in Eleocharis.

Molecular phylogeny of Japanese Eleocharis (Cyperaceae) based on ITS sequence data, and chromosomal evolution.

ITS sequence data were used to estimate the phylogeny of 24 Japanese Eleocharis species and to make karyomorphological observations on 19 of these taxa. Two major clades were identified in Japanese Eleocharis molecular phylogenetic trees: (1) one including all species of section Limnochloa, and (2) another comprising two sections, Pauciflorae and Eleocharis. Phylogenetic analysis including both Japanese and North American species also shows strong support for monophyly of the Mutatae/ Limnochloa clade. The width of the spikelets in the species Mutatae/ Limnochloa is the same as that of the culms, indicating that the relative widths of spikelets and culms are useful characteristics for classification. Two major clades were supported by karyomorphological data. All taxa of section Limnochloa had very small chromosomes, while sections Pauciflorae and Eleocharis had large chromosomes. The basic chromosome number of sections Eleocharis and Pauciflorae is thought to be x=5. Chromosomal evolution in the genus Eleocharis with diffuse centromeric chromosomes may be caused by both aneuploidization and polyploidization. Our data suggest that a 3-bp insertion near the 3' end of the 5.8S gene is useful for intrageneric delimitations of the genus Eleocharis.