Search

Genome‐wide 5hmC profiling in cell‐free DNA acquired from CRC patients, adenoma patients, and healthy individuals revealed that the differential 5hmC‐modified regions were gathered into four clusters with no overlap, although there are a few overlapped genes shared between the different clusters. CRC patients with adenoma history showed exclusive 5hmC‐gain characteristics, which was consistent with the ‘parallel’ evolution hypothesis in adenoma., Approximately 85% colorectal cancers (CRCs) are thought to evolve through the adenoma‐to‐carcinoma sequence associated with specific molecular alterations, including the 5‐hydroxymethylcytosine (5hmC) signature in circulating cell‐free DNA (cfDNA). To explore colorectal disease progression and evaluate the use of cfDNA as a potential diagnostic factor for CRC screening, here, we performed genome‐wide 5hmC profiling in plasma cfDNA and tissue genomic DNA (gDNA) acquired from 101 samples (63 plasma and 38 tissues), collected from 21 early‐stage CRC patients, 21 AD patients, and 21 healthy controls (HC). The gDNA and cfDNA 5hmC signatures identified in gene bodies and promoter regions in CRC and AD groups were compared with those in HC group. All the differential 5hmC‐modified regions (DhMRs) were gathered into four clusters: Disease‐enriched, AD‐enriched, Disease‐lost, and AD‐lost, with no overlap. AD‐related clusters, AD‐enriched and AD‐lost, displayed the unique 5hmC signals in AD patients. Disease‐enriched and Disease‐lost clusters indicated the general 5hmC changes when colorectal lesions occurred. Cancer patients with a confirmable adenoma history segmentally gathered in AD‐enriched clusters. KEGG functional enrichment and GO analyses determined distinct differential 5hmC‐modified profiles in cfDNA of HC individuals, AD, and CRC patients. All patients had comprehensive 5hmC signatures where Disease‐enriched and Disease‐lost DhMR clusters demonstrated similar epigenetic modifications, while AD‐enriched and AD‐lost DhMR clusters indicated complicated subpopulations in adenoma. Analysis of CRC patients with adenoma history showed exclusive 5hmC‐gain characteristics, consistent with the ‘parallel’ evolution hypothesis in adenoma, either developed through the adenoma‐to‐carcinoma sequence or not. These findings deepen our understanding of colorectal disease and suggest that the 5hmC modifications of different pathological subtypes (cancer patients with or without adenoma history) could be used to screen early‐stage CRC and assess adenoma malignancy with large‐scale follow‐up studies in the future.

Background Panax ginseng C. A. Mey is one of famous medicinal herb plant species. Its major bioactive compounds are various ginsenosides in roots and rhizomes. It is commonly accepted that ginsenosides are synthesized from terpene precursors, IPP and DMAPP, through the cytoplasmic mevalonate (MVA) pathway. Another plastic 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway was proved also contributing to ginsenoside generation in the roots of P. ginseng by using specific chemical inhibitors recently. But their gene expression characteristics are still under reveal in P. ginseng. With the development of the high-throughput next generation sequencing (NGS) technologies, we have opportunities to discover more about the complex ginsenoside biosynthesis pathways in P. ginseng. Results We carried out deep RNA sequencing and comprehensive analyses on the ginseng root samples of 1–5 years old and five different tissues of 5 years old ginseng plants. The de novo assembly totally generated 48,165 unigenes, including 380 genes related to ginsenoside biosynthesis and all the genes encoding the enzymes of the MEP pathway and the MVA pathway. We further illustrated the gene expression profiles related to ginsenoside biosynthesis among 1–5 year-old roots and different tissues of 5 year-old ginseng plants. Particularly for the first time, we revealed that the gene transcript abundances of the MEP pathway were similar to those of the MVA pathway in ginseng roots but higher in ginseng leaves. The IspD was predicated to be the rate-limiting enzyme in the MEP pathway through both co-expression network and gene expression profile analyses. Conclusions At the transcriptional level, the MEP pathway has similar contribution to ginsenoside biosynthesis in ginseng roots, but much higher in ginseng leaves, compared with the MVA pathway. The IspD might be the key enzyme for ginsenoside generation through the MEP pathway. These results provide new information for further synthetic biology study on ginsenoside metabolic regulation. Electronic supplementary material The online version of this article (10.1186/s12864-019-5718-x) contains supplementary material, which is available to authorized users.

The carbohydrate-rich cell walls of land plants and algae have been the focus of much interest given the value of cell wall-based products to our current and future economies. Hydroxyproline-rich glycoproteins (HRGPs), a major group of wall glycoproteins, play important roles in plant growth and development, yet little is known about how they have evolved in parallel with the polysaccharide components of walls. We investigate the origins and evolution of the HRGP superfamily, which is commonly divided into three major multigene families: the arabinogalactan proteins (AGPs), extensins (EXTs), and proline-rich proteins. Using motif and amino acid bias, a newly developed bioinformatics pipeline, we identified HRGPs in sequences from the 1000 Plants transcriptome project (www.onekp.com). Our analyses provide new insights into the evolution of HRGPs across major evolutionary milestones, including the transition to land and the early radiation of angiosperms. Significantly, data mining reveals the origin of glycosylphosphatidylinositol (GPI)-anchored AGPs in green algae and a 3- to 4-fold increase in GPI-AGPs in liverworts and mosses. The first detection of cross-linking (CL)-EXTs is observed in bryophytes, which suggests that CL-EXTs arose though the juxtaposition of preexisting SPn EXT glycomotifs with refined Y-based motifs. We also detected the loss of CL-EXT in a few lineages, including the grass family (Poaceae), that have a cell wall composition distinct from other monocots and eudicots. A key challenge in HRGP research is tracking individual HRGPs throughout evolution. Using the 1000 Plants output, we were able to find putative orthologs of Arabidopsis pollen-specific GPI-AGPs in basal eudicots.

The common carp, Cyprinus carpio, is one of the most important cyprinid species and globally accounts for 10% of freshwater aquaculture production. Here we present a draft genome of domesticated C. carpio (strain Songpu), whose current assembly contains 52,610 protein-coding genes and approximately 92.3% coverage of its paleotetraploidized genome (2n = 100). The latest round of whole-genome duplication has been estimated to have occurred approximately 8.2 million years ago. Genome resequencing of 33 representative individuals from worldwide populations demonstrates a single origin for C. carpio in 2 subspecies (C. carpio Haematopterus and C. carpio carpio). Integrative genomic and transcriptomic analyses were used to identify loci potentially associated with traits including scaling patterns and skin color. In combination with the high-resolution genetic map, the draft genome paves the way for better molecular studies and improved genome-assisted breeding of C. carpio and other closely related species.

Three Gracilaria species, G. chouae, G. blodgettii, G. vermiculophylla and a close relative species, Gracilariopsis lemaneiformis which is now nominated as Gracilaria lemaneiformis, are the typically indigenous species which are important resources for the production of special proteins, phycobilisomes, special carbohydrates, and agar in China. In this study, de novo transcriptome sequencing on these four species using the next generation sequencing technology was performed for the first time. Functional annotations on assembled sequencing reads showed that the transcriptomic profiles were quite different between G. lemaneiformis and other three Gracilaria species. Comparative analysis of differential gene expression related to carbohydrate and phycobiliprotein metabolisms also showed that the expression profiles of these essential genes were different in four species. The genes encoding allophycocyanin, phycocyanin and phycoerythrin were further examined in four species and their deduced amino acid sequences were used for phylogenetic analysis to confirm that G. lemaneiformis had close relationship to genus Gracilaria, as well as that within genus Gracilaria, G. chouae had closer relationship to G. vermiculophylla rather than to G. blodgettii. The de novo transcriptome study on four species provided a valuable genomic resource for further understanding and analysis on biological and evolutionary study among marine algae.

Species of Sargassum are distributed worldwide, and are of great ecological and economic importance in marine ecosystems and bioresources. In this study, transcriptome sequencings of six Sargassum species were performed for the first time using an Illumina platform. For each sample, a total of 2.1–2.5 Gb of nucleotides are collected and assembled into 69 871–116 790 scaffolds, with an average length of 410–550 bp and N50 length of 756–1 462 bp. A total of 20 512–28 684 unigenes of each sample were annotated and compared well with known gene sequences from nr database. Clusters of Orthologous Groups (COG), gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were also performed for further understanding of gene functions and regulation pathways. Gene expression levels were calculated based on RPKM values and compared among these species, especially for those genes related to carbohydrate metabolism. Cluster analyses indicated that the differences of global gene expression between S. fusiforme, which was nominated as Hizikia fusiformis before, and other five species were not significant. Further phylogenetic analysis of 108 orthologous genes confirmed that S. fusiforme had closer relationship with S. hemiphyllum rather than S. horneri. These transcriptome data provided valuable information for better understanding of genome and gene characteristics of Sargassum algae and benefiting comparative and phylogenetic studies of Phaeophyceae species in future studies.

In eukaryotes, protein deacetylation is carried out by two well-conserved his tone deacetylase (HDAC) families: RPD3/HDA1 and SIR2. Intriguingly, model plants such as Arabidopsis express an additional plant-specific HDAC family, termed type-2 HDACs (HD2s). Transcriptomic analyses from more than 1300 green plants generated by the 1000 plants (1KP) consortium showed that HD2s appeared early in green plant evolution, the first members being detected in several streptophyte green alga. The HD2 family has expanded via several rounds of successive duplication; members are expressed in all major green plant clades. Interestingly, angiosperm species express new HD2 genes devoid of a zinc-finger domain, one of the main structural features of HD2s. These variants may have been associated with the origin and/or the biology of the ovule/seed.

Both Calvin-Benson-Bassham (C3) and Hatch-Slack (C4) cycles are most important autotrophic CO2 fixation pathways on today’s Earth. C3 cycle is believed to be originated from cyanobacterial endosymbiosis. However, studies on evolution of different biochemical variants of C4 photosynthesis are limited to tracheophytes and origins of C4-cycle genes are not clear till now. Our comprehensive analyses on bioinformatics and phylogenetics of novel transcriptomic sequencing data of 21 rhodophytes and 19 Phaeophyceae marine species and public genomic data of more algae, tracheophytes, cyanobacteria, proteobacteria and archaea revealed the origin and evolution of C4 cycle-related genes. Almost all of C4-related genes were annotated in extensive algal lineages with proteobacterial or archaeal origins, except for phosphoenolpyruvate carboxykinase (PCK) and aspartate aminotransferase (AST) with both cyanobacterial and archaeal/proteobacterial origin. Notably, cyanobacteria may not possess complete C4 pathway because of the flawed annotation of pyruvate orthophosphate dikinase (PPDK) genes in public data. Most C4 cycle-related genes endured duplication and gave rise to functional differentiation and adaptation in different algal lineages. C4-related genes of NAD-ME (NAD-malic enzyme) and PCK subtypes exist in most algae and may be primitive ones, while NADP-ME (NADP-malic enzyme) subtype genes might evolve from NAD-ME subtype by gene duplication in chlorophytes and tracheophytes.

We generated hairy root cultures from two closely related species, Adenophora lobophylla and A. potaninii (Campanulaeae) and carried out a comparative study on their cadmium (Cd) response. A. lobophylla is an endangered species while A. potaninii is widely distributed in the same habitat. Upon exposure to Cd concentrations higher than 50 μM, more extensive growth inhibition and higher Cd accumulation were detected in the hairy root of A. lobophylla. Cd treatment affected the protein content in both the species. Phytochelatins (PCs) have been isolated and characterized from the hairy roots for both species. They shared structure similarities but showed different accumulation kinetics. The content of reduced glutathione (GSH) and cysteine (Cys) differs in both the species and they show different changes upon Cd challenge. The results suggested that these two species might employ different strategy for Cd detoxification. A. lobophylla is capable of synthesizing high level of PCs while a Cd exclusion system and a tighter homeostasis mechanism(s) to maintain the cellular GSH level could have been evolved in A. potanini in additon to its capability of synthesizing PCs.