The Phosphatome of Medicinal and Edible Fungus Wolfiporia cocos.

Wolfiporia cocos is an important medicinal and edible fungus that grows in association with pine trees, and its dried sclerotium has been used as a traditional medicine in China for centuries. However, the commercial production of W. cocos sclerotia is currently limited by shortages in pine wood resources. Since protein phosphatases (PPs) play significant roles in growth, signal transduction, development, metabolism, sexual reproduction, cell cycle, and environmental stress responses in fungi, the phosphatome of W. cocos was analyzed in this study by identifying PP genes, studying transcript profiles and assigning PPs to orthologous groups. Fifty-four putative PP genes were putatively identified in W. cocos genome based on homologous sequences searching using BLASTx program against the Saccharomyces cerevisiae, Fusarium graminearum, and Sclerotinia sclerotiorum databases. Based on known and presumed functions of orthologues of these PP genes found in other fungi, the putative roles of these W. cocos PPs in colonization, hyphal growth, sclerotial formation, secondary metabolism, and stress tolerance to environment were discussed in this study. And the level of transcripts from PP genes in the mycelium and sclerotium stages was also analyzed by qRT-PCR. Our study firstly identified and functional discussed the phosphatome in the medicinal and edible fungus W. cocos. The data from our study contribute to a better understanding of PPs potential roles in various cellar processes of W. cocos, and systematically provide comprehensive and novel insights into W. cocos economically important traits that could be extended to other fungi.

Functional divergence and origin of the DAG-like gene family in plants.

The nuclear-encoded DAG-like (DAL) gene family plays critical roles in organelle C-to-U RNA editing in Arabidopsis thaliana. However, the origin, diversification and functional divergence of DAL genes remain unclear. Here, we analyzed the genomes of diverse plant species and found that: DAL genes are specific to spermatophytes, all DAL genes share a conserved gene structure and protein similarity with the inhibitor I9 domain of subtilisin genes found in ferns and mosses, suggesting that DAL genes likely arose from I9-containing proproteases via exon shuffling. Based on phylogenetic inference, DAL genes can be divided into five subfamilies, each composed of putatively orthologous and paralogous genes from different species, suggesting that all DAL genes originated from a common ancestor in early seed plants. Significant type I functional divergence was observed in 6 of 10 pairwise comparisons, indicating that shifting functional constraints have contributed to the evolution of DAL genes. This inference is supported by the finding that functionally divergent amino acids between subfamilies are predominantly located in the DAL domain, a critical part of the RNA editosome. Overall, these findings shed light on the origin of DAL genes in spermatophytes and outline functionally important residues involved in the complexity of the RNA editosome.

De novo assembly and transcriptome analysis of sclerotial development in Wolfiporia cocos.

Wolfiporia cocos Ryvarden et Gilbertson, a well-known medicinal fungus in the Basidiomycetes, is widely distributed in East Asia. Its dried sclerotium, which is known as Fuling in China, has been used as a traditional crude drug in Chinese traditional medicine for thousand years. However, little is known about how the sclerotium is developed at the genetic level. In this study, the de novo sequencing of sclerotia of W. cocos (S1_initial stage; S2_developmental stage and S3_mature stage) was carried out by illumina HiSeq 2000 technology. 27,438 unigenes were assembled from ~30Gbp raw data, and 12,093 unigenes were significantly annotated. The analysis of expression profiles during development returned 304 differentially expressed genes (DEGs), which were clustered into four different groups according to their expression trends. Especially for the maturation stage (S3), the sclerotium exhibited a markedly different expression profile from other stages. We further showed that peroxisome, unsaturation of fatty acids and degradation pathway were respectively prevalent in S1, S2 and S3 stages as evidenced by enrichment analysis. To our knowledge, this study represents the first report of sclerotial development transcriptomics in W. cocos. The obtained results provide novel insights into the developmental biology of the sclerotia, which is helpful for future studies about cultivation and breeding of W. cocos.

De Novo Analysis of Wolfiporia cocos Transcriptome to Reveal the Differentially Expressed Carbohydrate-Active Enzymes (CAZymes) Genes During the Early Stage of Sclerotial Growth.

The sclerotium of Wolfiporia cocos has been used as an edible mushroom and/or a traditional herbal medicine for centuries. W. cocos sclerotial formation is dependent on parasitism of the wood of Pinus species. Currently, the sclerotial development mechanisms of W. cocos remain largely unknown and the lack of pine resources limit the commercial production. The CAZymes (carbohydrate-active enzymes) play important roles in degradation of the plant cell wall to provide carbohydrates for fungal growth, development, and reproduction. In this study, the transcript profiles from W. cocos mycelium and 2-months-old sclerotium, the early stage of sclerotial growth, were specially analyzed using de novo sequencing technology. A total of 142,428,180 high-quality reads of mycelium and 70,594,319 high-quality reads of 2-months-old sclerotium were obtained. Additionally, differentially expressed genes from the W. cocos mycelium and 2-months-old sclerotium stages were analyzed, resulting in identification of 69 CAZymes genes which were significantly up-regulated during the early stage of sclerotial growth compared to that of in mycelium stage, and more than half of them belonged to glycosyl hydrolases (GHs) family, indicating the importance of W. cocos GHs family for degrading the pine woods. And qRT-PCR was further used to confirm the expression pattern of these up-regulated CAZymes genes. Our results will provide comprehensive CAZymes genes expression information during W. cocos sclerotial growth at the transcriptional level and will lay a foundation for functional genes studies in this fungus. In addition, our study will also facilitate the efficient use of limited pine resources, which is significant for promoting steady development of Chinese W. cocos industry.

Polysaccharides from Rhizoma Panacis Majoris and its anti-oxidant activity.

Response surface method (RSM) was employed to optimize the extraction conditions of polysaccharides from Rhizoma Panacis Majoris (the rhizomes of Panax japonicus C. A. Mey. var. major (Burk.) C. Y. Wu et K. M. Feng) (RPMP), a well-known Chinese traditional medicine. In order to obtain the optimal processing parameters, a three-variable Box-Behnken designs (BBD) were applied for experimental designs. RSM analysis indicated the good correspondence between experimental and predicted values, the optimal conditions for the yield of polysaccharides were as follows: the ultrasound time is 31.15 min, extraction temperature is 92.50 °C, and the ratio of water to raw material is 40 mL/g. The maximum value (13.87 ± 0.16%) of the yield of polysaccharides was obtained under these optimal conditions. The molecular weight (MW) was determined to be 1.48 × 10(5)(±0.39%)Da by HPSEC-MALLS-RID chromatography system. FT-IR spectra demonstrated obvious characteristic peaks of polysaccharides. The antioxidant activities of RPMP were investigated including scavenging activity of hydrogen radicals, ABTS radicals, and free radicals of superoxide anion in vitro, and the results exhibited that RPMP had a good potential for antioxidant.

[Transcriptome profiling and analysis of Panax japonicus var. major].

The rhizome of Panax japonicus var. major have been used as the natural medicinal agent by Chinese traditional doctors for more than thousand years. Most of the therapeutic effects of P. japonicus var. major had been reported due to the presence of tetracyclic or pentacyclic triterpene saponins. In this study, Illumina pair-end RNA-sequencing and de novo splicing were done in order to understand the pathway of triterpenoid saponins in this species. The valid reads data of 15. 6 Gb were obtained. The 62 240 unigenes were finally obtained by de novo splicing. After annotation, we discovered 19 unigenes involved in ginsenoside backbone biosynthesis. Additionally, 69 unigenes and 18 unigenes were predicted to have potential function of cytochrome P450 and UDP-glycosyltransferase based on the annotation results, which may encode enzymes responsible for ginsenoside backbone modification. This study provides global expressed datas for P. japonicus var. major, which will contribute significantly to further genome-wide research and analysis for this species.

KRN4 Controls Quantitative Variation in Maize Kernel Row Number.

Kernel row number (KRN) is an important component of yield during the domestication and improvement of maize and controlled by quantitative trait loci (QTL). Here, we fine-mapped a major KRN QTL, KRN4, which can enhance grain productivity by increasing KRN per ear. We found that a ~3-Kb intergenic region about 60 Kb downstream from the SBP-box gene Unbranched3 (UB3) was responsible for quantitative variation in KRN by regulating the level of UB3 expression. Within the 3-Kb region, the 1.2-Kb Presence-Absence variant was found to be strongly associated with quantitative variation in KRN in diverse maize inbred lines, and our results suggest that this 1.2-Kb transposon-containing insertion is likely responsible for increased KRN. A previously identified A/G SNP (S35, also known as Ser220Asn) in UB3 was also found to be significantly associated with KRN in our association-mapping panel. Although no visible genetic effect of S35 alone could be detected in our linkage mapping population, it was found to genetically interact with the 1.2-Kb PAV to modulate KRN. The KRN4 was under strong selection during maize domestication and the favorable allele for the 1.2-Kb PAV and S35 has been significantly enriched in modern maize improvement process. The favorable haplotype (Hap1) of 1.2-Kb-PAV-S35 was selected during temperate maize improvement, but is still rare in tropical and subtropical maize germplasm. The dissection of the KRN4 locus improves our understanding of the genetic basis of quantitative variation in complex traits in maize.

De novo characterization of Panax japonicus C. A. Mey transcriptome and genes related to triterpenoid saponin biosynthesis.

Panax japonicus C.A.Mey, the traditional medicinal herb in the Araliaceae family, has been used as a tonic, anti-inflammatory and haemostatic agent in China for more than thousand years. Its clinical effects are mainly due to the presence of triterpenoid saponins. However, little is known at the genetic level about how saponins are biosynthesized in this plant. We have therefore performed the de novo transcriptome assembly and high throughput RNA-seq analysis for P. japonicus. 66,403 unigenes were assembled from 19.6 Gbp raw data, and 34,639 unigenes were annotated. After annotation, 29 unigenes involved in putative backbone biosynthesis of triterpenoid saponin were selected. Additionally, 34 Cytochrome P450 and 18 UDP-glycosyltransferase unigenes were predicted based on the annotation, which were related to the saponin backbone modification. The expression level of related key genes were further verified by qPCR analysis. The results of this study provide the most comprehensive expressed sequence resources for P. japonicus, which will enlarge the available P. japonicus gene pool and facilitate further genome-wide research and analyses in this species.

Multiple loci not only Rf3 involved in the restoration ability of pollen fertility, anther exsertion and pollen shedding to S type cytoplasmic male sterile in maize.

KEY MESSAGE: Thirty loci for fertility restoration of pollen fertility, anther exsertion and pollen shedding to maize CMS-S were identified by GWAS. S type cytoplasmic male sterile (CMS-S) is the main type of CMS in maize; poor understanding of the genetic architecture of fertility restoration to CMS-S is one of the reasons to impede its utility in hybrid breeding. In this study, genome-wide identification of genetic loci for fertility restoration ability to CMS-S was firstly conducted with a set of testcrossing association mapping panel in three environments. A total of 19, 3 and 8 significant loci (P < 1.8 × 10(-6), α = 1) for pollen fertility, anther exsertion and pollen shedding were identified, respectively, and individual locus explained up to 28.26% of phenotypic variation. Of them, only Rf3, the main restorer-fertility gene of CMS-S, was identified for the three traits simultaneously. In addition, 83 candidate genes within the 100 kb extension regions of these loci were predicted. These results revealed that besides Rf3 multiple genetic loci and mechanisms are involved in the fertility restoration ability to CMS-S. Results in this study would provide important information for understanding the genetic architecture of fertility restoration to CMS-S in maize.

Genetic architecture of maize kernel row number and whole genome prediction.

KEY MESSAGE: Maize kernel row number might be dominated by a set of large additive or partially dominant loci and several small dominant loci and can be accurately predicted by fewer than 300 top KRN-associated SNPs. Kernel row number (KRN) is an important yield component in maize and directly affects grain yield. In this study, we combined linkage and association mapping to uncover the genetic architecture of maize KRN and to evaluate the phenotypic predictability using these detected loci. A genome-wide association study revealed 31 associated single nucleotide polymorphisms (SNPs) representing 17 genomic loci with an effect in at least one of five individual environments and the best linear unbiased prediction (BLUP) over all environments. Linkage mapping in three F2:3 populations identified 33 KRN quantitative trait loci (QTLs) representing 21 QTLs common to several population/environments. The majority of these common QTLs that displayed a large effect were additive or partially dominant. We found 70% KRN-associated genomic loci were mapped in KRN QTLs identified in this study, KRN-associated SNP hotspots detected in NAM population and/or previous identified KRN QTL hotspots. Furthermore, the KRN of inbred lines and hybrids could be predicted by the additive effect of the SNPs, which was estimated using inbred lines as a training set. The prediction accuracy using the top KRN-associated tag SNPs was obviously higher than that of the randomly selected SNPs, and approximately 300 top KRN-associated tag SNPs were sufficient for predicting the KRN of the inbred lines and hybrids. The results suggest that the KRN-associated loci and QTLs that were detected in this study show great potential for improving the KRN with genomic selection in maize breeding.

An Efficient PEG/CaCl2-Mediated Transformation Approach for the Medicinal Fungus Wolfiporia cocos.

Sclerotia of Wolfiporia cocos are of medicinal and culinary value. The genes and molecular mechanisms involved in W. cocos sclerotial formation are poorly investigated because of the lack of a suitable and reproducible transformation system for W. cocos. In this study, a PEG/ CaCl2-mediated genetic transformation system for W. cocos was developed. The promoter Pgpd from Ganoderma lucidum effectively drove expression of the hygromycin B phosphotransferase gene in W. cocos, and approximately 30 transformants were obtained per 10 µg DNA when the protoplast suspension density was 10(6) protoplasts/ml. However, no transformants were obtained under the regulation of the PtrpC promoter from Aspergillus nidulans.

Resource investigation of traditional medicinal plant Panax japonicus (T.Nees) C.A. Mey and its varieties in China.

ETHNOPHARMACOLOGICAL RELEVANCE: Panax japonicus, the perennial herb in the Araliaceae family, was used as the natural medicinal herb by Chinese traditional doctors for more than thousand years. Its rhizome was mainly used as a tonic, anti-inflammatory and hemostatic agent in China. Most of the therapeutic effects of P. japonicus had been reported due to the presence of tetracyclic or pentacyclic triterpene saponins. Volatile oil, polysaccharides and amino acids had also been found in P. japonicus species and reported in the pharmacological functions. AIM OF THE STUDY: A three-year survey was conducted to determine the current resource status of P. japonicus (T.Nees) C. A. Mey and its varieties (P. japonicus var. major (Burkill) C.Y.Wu & Feng and P. japonicus var. bipinnatifidus (Seem.) C.Y.Wu & Feng) in 10 provinces of southern and southwestern China. METHODS AND RESULTS: Whole plants were sampled at 64 sites. Resource distribution, habitat type, morphological variation and market trend of them were studied and discussed. The natural resource in China is rarely available due to extensive exploitation and continual environment deterioration in recent decades, Abundance of P. japonicus was much lower than previous records, mainly found in Hubei, Sichuan, Guizhou and Yunnan province. Wild resources of P.japonicus var. major and P.japonicus var. bipinnatifidus were even scarcer, only found in Guizhou and Yunan province. Despite their dramatic rise of market trend, the artificial cultivation of them was still not fully developed in China, but progressed rapidly in Hubei province. CONCLUSION: In this study, we synthesized our understandings of the current resource state of P. japonicus׳s existence, variation and cultivation in China. This study will aid further investigations and increased protection of these plants, which are very valuable to traditional herbal medicine.

Biomass enzymatic saccharification is determined by the non-KOH-extractable wall polymer features that predominately affect cellulose crystallinity in corn.

Corn is a major food crop with enormous biomass residues for biofuel production. Due to cell wall recalcitrance, it becomes essential to identify the key factors of lignocellulose on biomass saccharification. In this study, we examined total 40 corn accessions that displayed a diverse cell wall composition. Correlation analysis showed that cellulose and lignin levels negatively affected biomass digestibility after NaOH pretreatments at p<0.05 & 0.01, but hemicelluloses did not show any significant impact on hexoses yields. Comparative analysis of five standard pairs of corn samples indicated that cellulose and lignin should not be the major factors on biomass saccharification after pretreatments with NaOH and H2SO4 at three concentrations. Notably, despite that the non-KOH-extractable residues covered 12%-23% hemicelluloses and lignin of total biomass, their wall polymer features exhibited the predominant effects on biomass enzymatic hydrolysis including Ara substitution degree of xylan (reverse Xyl/Ara) and S/G ratio of lignin. Furthermore, the non-KOH-extractable polymer features could significantly affect lignocellulose crystallinity at p<0.05, leading to a high biomass digestibility. Hence, this study could suggest an optimal approach for genetic modification of plant cell walls in bioenergy corn.

Genome-wide identification, evolution, and expression analysis of RNA-binding glycine-rich protein family in maize.

The RNA-binding glycine-rich protein (RB-GRP) family is characterized by the presence of a glycine-rich domain arranged in (Gly)n-X repeats and an RNA-recognition motif (RRM). RB-GRPs participate in varied physiological and biochemical processes especially in the stress response of plants. In this study, a total of 23 RB-GRPs distributed on 10 chromosomes were identified in maize (Zea mays L.), and they were divided into four subgroups according to their conserved domain architecture. Five pairs of paralogs were identified, while none of them was located on the same chromosomal region, suggesting that segmental duplication is predominant in the duplication events of the RB-GRPs in maize. Comparative analysis of RB-GRPs in maize, Arabidopsis (Arabidopsis thaliana L.), rice (Oryza sativa L.), and wheat (Triticum aestivum) revealed that two exclusive subgroups were only identified in maize. Expression of eight ZmRB-GRPs was significantly regulated by at least two kinds of stresses. In addition, cis-elements predicted in the promoter regions of the ZmRB-GRPs also indicated that these ZmRB-GRPs would be involved in stress response of maize. The preliminary genome-wide analysis of the RB-GRPs in maize would provide useful information for further study on the function of the ZmRB-GRPs.

Characterization and antioxidant activities of degraded polysaccharides from Poria cocos sclerotium.

Poria cocos F.A.Wolf is a Chinese traditional medicine used to treat chronic gastritis, edema, nephrosis, gastric atony, and acute gastroenteric catarrh. Polysaccharides are the main active component of P. cocos. We obtained polysaccharides PCP-1, PCP-2, and PCP-3 from the degradation of P. cocos polysaccharides (PCP) with different concentrations of H2O2 solution. Molecular weights were determined by high performance size exclusion chromatography. HPLC analysis of monosaccharide composition confirmed that PCP-1, PCP-2, and PCP-3 are heteropolysaccharides composed of glucose and arabinose. IR spectra indicated obvious characteristic peaks of polysaccharides. The antioxidant activities of these polysaccharides were evaluated by established in vitro systems, including scavenging activity of hydroxyl radicals, ABTS radicals, and ferrous ions. The degradation polysaccharides exhibited obvious and concentration-dependent antioxidant properties. In addition, DNA binding analysis showed that PCP-1 had a stronger capacity than other polysaccharides to interact with DNA. However, each polysaccharide had a certain capacity for DNA damage protection.

Genome-wide identification, evolution and expression analysis of mTERF gene family in maize.

Plant mitochondrial transcription termination factor (mTERF) genes comprise a large family with important roles in regulating organelle gene expression. In this study, a comprehensive database search yielded 31 potential mTERF genes in maize (Zea mays L.) and most of them were targeted to mitochondria or chloroplasts. Maize mTERF were divided into nine main groups based on phylogenetic analysis, and group IX represented the mitochondria and species-specific clade that diverged from other groups. Tandem and segmental duplication both contributed to the expansion of the mTERF gene family in the maize genome. Comprehensive expression analysis of these genes, using microarray data and RNA-seq data, revealed that these genes exhibit a variety of expression patterns. Environmental stimulus experiments revealed differential up or down-regulation expression of maize mTERF genes in seedlings exposed to light/dark, salts and plant hormones, respectively, suggesting various important roles of maize mTERF genes in light acclimation and stress-related responses. These results will be useful for elucidating the roles of mTERF genes in the growth, development and stress response of maize.

Genetic analysis and major QTL detection for maize kernel size and weight in multi-environments.

KEY MESSAGE: Twelve major QTL in five optimal clusters and several epistatic QTL are identified for maize kernel size and weight, some with pleiotropic will be promising for fine-mapping and yield improvement. Kernel size and weight are important target traits in maize (Zea mays L.) breeding programs. Here, we report a set of quantitative trait loci (QTL) scattered through the genome and significantly controlled the performance of four kernel traits including length, width, thickness and weight. From the cross V671 (large kernel) × Mc (small kernel), 270 derived F2:3 families were used to identify QTL of maize kernel-size traits and kernel weight in five environments, using composite interval mapping (CIM) for single-environment analysis along with mixed linear model-based CIM for joint analysis. These two mapping strategies identified 55 and 28 QTL, respectively. Among them, 6 of 23 coincident were detected as interacting with environment. Single-environment analysis showed that 8 genetic regions on chromosomes 1, 2, 4, 5 and 9 clustered more than 60 % of the identified QTL. Twelve stable major QTLs accounting for over 10 % of phenotypic variation were included in five optimal clusters on the genetic region of bins 1.02-1.03, 1.04-1.06, 2.05-2.07, 4.07-4.08 and 9.03-9.04; the addition and partial dominance effects of significant QTL play an important role in controlling the development of maize kernel. These putative QTL may have great promising for further fine-mapping with more markers, and genetic improvement of maize kernel size and weight through marker-assisted breeding.

Polysaccharides from Panax japonicus C.A. Meyer and their antioxidant activities.

Polysaccharides named as PP1, PP2, PP3, PP4 and PP5 were extracted and isolated from the rhizomes of Panax japonicus C.A. Meyer, a well-known Chinese traditional medicine, by controlling the final concentration of solution to precipitate the polysaccharides. The molecular weight of polysaccharides was determined by HGPLC chromatography system. The monosaccharide composition was analyzed by Gas chromatography on an Agilent 6890A instrument using a DB-35MS column and flame-ionization detector. All of the polysaccharides were heteropolysaccharides and consisted of arabinose, glucose and galactose. The content of arabinose increased with the increasing of ethanol concentration and PP5 had the most arabinose content in these samples. IR spectra indicated obvious characteristic peaks of polysaccharide, the presence of uronic acids. Their antioxidant activities were evaluated by various established in vitro systems, including scavenging activity of superoxide anion, hydroxyl radical, ABTS and DPPH radicals. Both samples showed inhibitory effects on superoxide, hydroxyl, ABTS and DPPH radical. And PP5 shows more clearly and relatively stronger capacity than other polysaccharides on the protective effect of DNA damage.

Dynamic QTL analysis and candidate gene mapping for waterlogging tolerance at maize seedling stage.

Soil waterlogging is one of the major abiotic stresses adversely affecting maize growth and yield. To identify dynamic expression of genes or quantitative trait loci (QTL), QTL associated with plant height, root length, root dry weight, shoot dry weight and total dry weight were identified via conditional analysis in a mixed linear model and inclusive composite interval mapping method at three respective periods under waterlogging and control conditions. A total of 13, 19 and 23 QTL were detected at stages 3D|0D (the period during 0-3 d of waterlogging), 6D|3D and 9D|6D, respectively. The effects of each QTL were moderate and distributed over nine chromosomes, singly explaining 4.14-18.88% of the phenotypic variation. Six QTL (ph6-1, rl1-2, sdw4-1, sdw7-1, tdw4-1 and tdw7-1) were identified at two consistent stages of seedling development, which could reflect a continuous expression of genes; the remaining QTL were detected at only one stage. Thus, expression of most QTL was influenced by the developmental status. In order to provide additional evidence regarding the role of corresponding genes in waterlogging tolerance, mapping of Expressed Sequence Tags markers and microRNAs were conducted. Seven candidate genes were observed to co-localize with the identified QTL on chromosomes 1, 4, 6, 7 and 9, and may be important candidate genes for waterlogging tolerance. These results are a good starting point for understanding the genetic basis for selectively expressing of QTL in different stress periods and the common genetic control mechanism of the co-localized traits.