Transcriptome Profile Analysis from Different Sex Types of Ginkgo biloba L

Expression profiling of the phenylalanine ammonia-lyase (PAL) gene family in ginkgo biloba L

The PAL gene family plays an important role in plant growth, development, and response to abiotic stresses and has been identified in a variety of plants. However, a systematic characterization is still lacking in Ginkgo biloba. Using a bioinformatics approach, 11 GbPAL members of the PAL gene family identified in ginkgo were identified in this study. The protein structure and physicochemical properties indicated that the GbPAL genes were highly similar. Based on their exon-intron structures, they can be classified into three groups. A total of 62 cis-elements for hormone, light, and abiotic stress responses were identified in the promoters of GbPAL genes, indicating that PAL is a multifunctional gene family. GbPAL genes were specifically expressed in different tissues and ploidy of ginkgo. These results provide a theoretical basis for further studies on the functional expression of the GbPAL genes.

Effects of High Temperature-Triggered Transcriptomics on the Physiological Adaptability of Cenococcum geophilum, an Ectomycorrhizal Fungus

High temperature stress caused by global warming presents a challenge to the healthy development of forestry. Cenococcum geophilum is a common ectomycorrhizal fungus (ECMF) in the forest system and has become an important fungus resource with application potential in forest vegetation restoration. In this study, three sensitive isolates of C. geophilum (ChCg01, JaCg144 and JaCg202) and three tolerant isolates of C. geophilum (ACg07, ChCg28 and ChCg100) were used to analyze the physiological and molecular responses to high temperature. The results showed that high temperature had a significant negative effect on the growth of sensitive isolates while promoting the growth of tolerant isolates. The antioxidative enzymes activity of C. geophilum isolates increased under high temperature stress, and the SOD activity of tolerant isolates (A07Cg and ChCg100) was higher than that of sensitive isolates (ChCg01 and JaCg202) significantly. The tolerant isolates secreted more succinate, while the sensitive isolates secreted more oxalic acid under high temperature stress. Comparative transcriptomic analysis showed that differentially expressed genes (DEGs) of six C. geophilum isolates were significantly enriched in "antioxidant" GO entry in the molecular. In addition, the "ABC transporters" pathway and the "glyoxylate and dicarboxylic acid metabolic" were shared in the three tolerant isolates and the three sensitive isolates, respectively. These results were further verified by RT-qPCR analysis. In conclusion, our findings suggest that C. geophilum can affect the organic acid secretion and increase antioxidant enzyme activity in response to high temperature by upregulating related genes.

NBLAST: a graphical user interface-based two-way BLAST software with a dot plot viewer

BLAST, a basic bioinformatics tool for searching local sequence similarity, has been one of the most widely used bioinformatics programs since its introduction in 1990. Users generally use the web-based NCBI-BLAST program for BLAST analysis. However, users with large sequence data are often faced with a problem of upload size limitation while using the web-based BLAST program. This proves inconvenient as scientists often want to run BLAST on their own data, such as transcriptome or whole genome sequences. To overcome this issue, we developed NBLAST, a graphical user interface-based BLAST program that employs a two-way system, allowing the use of input sequences either as "query" or "target" in the BLAST analysis. NBLAST is also equipped with a dot plot viewer, thus allowing researchers to create custom database for BLAST and run a dot plot similarity analysis within a single program. It is available to access to the NBLAST with http://nbitglobal.com/nblast.

Co-expression network analysis of genes and networks associated with wheat pistillody

Crop male sterility has great value in theoretical research and breeding application. HTS-1, whose stamens transformed into pistils or pistil-like structures, is an important male sterility material selecting from Chinese Spring three-pistil (CSTP) wheat. However the molecular mechanism of pistillody development in HTS-1 remains a mystery. RNA-seq data of 11 wheat tissues were obtained from the National Center for Biotechnology Information (NCBI), including the stamens of CSTP and the pistils and pistillodic stamen of HTS-1. The Salmon program was utilized to quantify the gene expression levels of the 11 wheat tissues; and gene quantification results were normalized by transcripts per million (TPM). In total, 58,576 genes were used to construct block-wise network by co-expression networks analysis (WGCNA) R package. We obtained all of modules significantly associated with the 11 wheat tissues. AgriGO V2.0 was used to do Gene Ontology (GO) enrichment analysis; and genes and transcription factors (TFs) in these significant modules about wheat pistillody development were identified from GO enrichment results. Basic local alignment search tool (BLAST) was used to align HTS-1 proteins with the published pistillody-related proteins and TFs. Genes about wheat pistillody development were analyzed and validated by qRT-PCR. The MEturquoise, MEsaddlebrown, MEplum, MEcoral1, MElightsteelblue1, and MEdarkslateblue modules were significantly corelated to pistillodic stamen (correlation p < 0.05). Moreover, 206 genes related to carpel development (GO:0048440) or gynoecium development (GO:0048467) were identified only in the MEturquoise module by Gene Ontology (GO) analysis, and 42 of 206 genes were hub genes in MEturquoise module. qRT-PCR results showed that 38 of the 42 hub genes had highly expressed in pistils and pistillodic stamens than in stamens. A total of 15 pistillody development-related proteins were validated by BLAST. Transcription factors (TFs) were also analyzed in the MEturquoise module, and 618 TFs were identified. In total, 56 TFs from 11 families were considered to regulate the development of pistillodic stamen. The co-expression network showed that six of HB and three of BES1 genes were identified in 42 hub genes. This indicated that TFs played important roles in wheat pistillody development. In addition, there were 11 of ethylene-related genes connected with TFs or hub genes, suggesting the important roles of ethylene-related genes in pistillody development. These results provide important insights into the molecular interactions underlying pistillody development.

Systematic Identification and Expression Analysis of the Auxin Response Factor (ARF) Gene Family in Ginkgo biloba L

Auxin participates in various physiological and molecular response-related developmental processes and is a pivotal hormone that regulates phenotypic formation in plants. Auxin response factors (ARFs) are vital transcription factors that mediate downstream auxin signaling by explicitly binding to auxin-responsive genes' promoters. Here, to investigate the possible developmental regulatory functions of ARFs in Ginkgo biloba, through employing comprehensive bioinformatics, we recognized 15 putative GbARF members. Conserved domains and motifs, gene and protein structure, gene duplication, GO enrichment, transcriptome expression profiles, and qRT-PCR all showed that Group I and III members were highly conserved. Among them, GbARF10b and GbARF10a were revealed as transcriptional activators in the auxin response for the development of Ginkgo male flowers through sequences alignment, cis-elements analysis and GO annotation; the results were corroborated for the treatment of exogenous SA. Moreover, the GbARFs expansion occurred predominantly by segmental duplication, and most GbARFs have undergone purifying selection. The Ka/Ks ratio test identified the functional consistence of GbARF2a and GbARF2c, GbARF10b, and GbARF10a in tissue expression profiles and male flower development. In summary, our study established a new research basis for exploring Ginkgo GbARF members' roles in floral organ development and hormone response.

Integrative Analyses of Transcriptomics and Metabolomics in Sex Differentiation of Mulberry Flowers

Sex determination and sex differentiation of plants are important physiological processes of plant development. Mulberry (Morus indica L.) is an important economic tree being cultivated in sericulture countries, and mulberry leaf is commonly used for sericulture. The transcriptomic and metabolomic differences between the staminate flowers (SFs) and pistillate flowers (PFs) of mulberry were investigated by RNA sequencing and ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). Overall, we uncovered 4,230 genes and 209 metabolites are significantly differentially expressed between the SFs and PFs of mulberry. The combined transcriptomic and metabolomic analysis revealed these differentially expressed genes (DEGs) and differentially expressed metabolites (DEMs) are involved in flavonoid biosynthesis, galactose metabolism, plant–pathogen interaction, and starch and sucrose metabolism, and these detected DEGs and DEMs may be associated with sex differentiation of mulberry through the regulation of the enrichment pathways, such as the MAPK pathway, flavonoid biosynthesis, galactose metabolism, plant–pathogen interaction, and starch and sucrose metabolism. This study will provide a rich source for the analysis of the molecular mechanism of mulberry sex differentiation processes.

Construction of a reference transcriptome for the analysis of male sterility in sugi (Cryptomeria japonica D. Don) focusing on MALE STERILITY 1 (MS1)

Sugi (Cryptomeria japonica D. Don) is an important conifer used for afforestation in Japan. As the genome of this species is 11 Gbps, it is too large to assemble within a short timeframe. Transcriptomics is one approach that can address this deficiency. Here we designed a workflow consisting of three stages to de novo assemble transcriptome using Oases and Trinity. The three transcriptomic stage used were independent assembly, automatic and semi-manual integration, and refinement by filtering out potential contamination. We identified a set of 49,795 cDNA and an equal number of translated proteins. According to the benchmark set by BUSCO, 87.01% of cDNAs identified were complete genes, and 78.47% were complete and single-copy genes. Compared to other full-length cDNA resources collected by Sanger and PacBio sequencers, the extent of the coverage in our dataset was the highest, indicating that these data can be safely used for further studies. When two tissue-specific libraries were compared, there were significant expression differences between male strobili and leaf and bark sets. Moreover, subtle expression difference between male-fertile and sterile libraries were detected. Orthologous genes from other model plants and conifer species were identified. We demonstrated that our transcriptome assembly output (CJ3006NRE) can serve as a reference transcriptome for future functional genomics and evolutionary biology studies.

Genome-Wide Identification and Coexpression Network Analysis of DNA Methylation Pathway Genes and Their Differentiated Functions in Ginkgo biloba L

DNA methylation plays a vital role in diverse biological processes. DNA methyltransferases (DNMTs) genes and RNA-directed DNA methylation (RdDM)-related genes are key genes responsible for establishing and maintaining genome DNA methylation in plants. In the present study, we systematically identified nine GbDNMTs in Ginkgo biloba, including the three common families of GbMET1a/1b, GbCMT2, and GbDRMa/b/2a/2b/2c, and a fourth family—GbDNMT3—which is absent in most angiosperms. We also identified twenty RdDM-related genes, including four GbDCLs, six GbAGOs, and ten GbRDRs. Expression analysis of the genes showed the different patterns of individual genes, and 15 of 29 genes displayed expression change under five types of abiotic stress. Gene coexpression analysis and weighted gene co-expression network analysis (WGCNA) using 126 public transcriptomic datasets revealed that these genes were clustered into two groups. In group I, genes covered members from all six families which were preferentially expressed in the ovulate strobile and fruit. A gene ontology (GO) enrichment analysis of WGCNA modules indicated that group I genes were most correlated with the biological process of cell proliferation. Group II only consisted of RdDM-related genes, including GbDRMs, GbAGOs, and GbRDRs, but no GbDCLs, and these genes were specifically expressed in the cambium, suggesting that they may function in a dicer-like (DCL)-independent RdDM pathway in specific tissues. The gene module related to group II was most enriched in signal transduction, cell communication, and the response to the stimulus. These results demonstrate that gene family members could be conserved or diverged across species, and multi-member families in the same pathway may cluster into different modules to function differentially. The study provides insight into the DNA methylation genes and their possible functions in G. biloba, laying a foundation for the further study of DNA methylation in gymnosperms.

Screening and identification of miRNAs related to sexual differentiation of strobili in Ginkgo biloba by integration analysis of small RNA, RNA, and degradome sequencing

BACKGROUND

Ginkgo biloba, a typical dioecious plant, is a traditional medicinal plant widely planted. However, it has a long juvenile period, which severely affected the breeding and cultivation of superior ginkgo varieties.

RESULTS

In order to clarify the complex mechanism of sexual differentiation in G. biloba strobili. Here, a total of 3293 miRNAs were identified in buds and strobili of G. biloba, including 1085 known miRNAs and 2208 novel miRNAs using the three sequencing approaches of transcriptome, small RNA, and degradome. Comparative transcriptome analysis screened 4346 and 7087 differentially expressed genes (DEGs) in male buds (MB) _vs_ female buds (FB) and microstrobilus (MS) _vs_ ovulate strobilus (OS), respectively. A total of 6032 target genes were predicted for differentially expressed miRNA. The combined analysis of both small RNA and transcriptome datasets identified 51 miRNA-mRNA interaction pairs that may be involved in the process of G. biloba strobili sexual differentiation, of which 15 pairs were verified in the analysis of degradome sequencing.

CONCLUSIONS

The comprehensive analysis of the small RNA, RNA and degradome sequencing data in this study provided candidate genes and clarified the regulatory mechanism of sexual differentiation of G. biloba strobili from multiple perspectives.

Full-length sequencing of ginkgo transcriptomes for an in-depth understanding of flavonoid and terpenoid trilactone biosynthesis

Ginkgo biloba L. is regarded as the most ancient living tree, and its kernel has been used as a traditional Chinese medicine for more than 2,000 years. The leaf extracts of this tree have been among the bestselling herbal remedies in Western countries since the last century. To understand the biosynthesis of the pharmacologically active ingredients in G. biloba, flavonoids and terpenoid trilactones (TTLs), we sequenced the transcriptomes of G. biloba leaves, kernels and testae with Iso-Seq and RNA-Seq technologies and obtained 152,524 clean consensus reads. When these reads were used to improve the annotation of the G. biloba genome, 4,856 novel genes, 25,583 new isoforms of previously annotated genes and 4,363 lncRNAs were discovered. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses indicated that genes involved in growth, regulation and response to stress were more likely to be regulated by alternative splicing (AS) or alternative polyadenylation (APA), which represent the two most important posttranscriptional regulation mechanisms. It was found that some of the characterized genes involved in the biosynthesis of flavonoids and TTLs were also possibly regulated by AS and APA. Using phylogenetic and gene expression pattern analyses, some candidate genes for the biosynthesis of flavonoids and TTLs were screened. After qRT-PCR validation, the final candidate genes for flavonoid biosynthesis included three UDP-glycosyltransferases and one MYB transcription factor, while the candidate genes for TTL biosynthesis included two cytochrome P450 and one WRKY transcription factor. Our study suggested that Iso-Seq may play an important role in improving genome annotation, elucidating AS and APA mechanisms and discovering candidate genes involved in the biosynthesis of some secondary metabolites.

Cell Wall Bulking by Maleic Anhydride for Wood Durability Improvement

Wood is susceptible to swelling deformation and decay fungi due to moisture adsorption that originates from the dynamic nanopores of the cell wall and the abundant hydroxyl groups in wood components. This study employed as a modifier maleic anhydride (MAn), with the help of acetone as solvent, to diffuse into the wood cell wall, bulk nanopores, and further chemically bond to the hydroxyl groups of wood components, reducing the numbers of free hydroxyl groups and weakening the diffusion of water molecules into the wood cell wall. The derived MAn-bulked wood, compared to the control wood, presented a reduction in water absorptivity (RWA) of ~23% as well as an anti-swelling efficiency (ASE) of ~39% after immersion in water for 228 h, and showed an improvement in decay resistance of 81.42% against white-rot fungus and 69.79% against brown-rot fungus, respectively. The method of combined cell wall bulking and hydroxyl group bonding could effectively improve the dimensional stability and decay resistance with lower doses of modifier, providing a new strategy for wood durability improvement.

Comparative RNA-Seq Analysis of Differentially Expressed Genes in the Epididymides of Yak and Cattleyak

Background

Cattleyak are the Fl hybrids between (♀) yak (Bos grunniens) and (♂) cattle (Bos taurus). Cattleyak exhibit higher capability in adaptability to a harsh environment and display much higher performances in production than the yak and cattle. The cattleyak, however, are females fertile but males sterile. All previous studies greatly focused on testes tissues to study the mechanism of male infer-tility in cattleyak. However, so far, no transcriptomic study has been conducted on the epididymides of yak and cattleyak.

Objective

Our objective was to perform comparative transcriptome analysis between the epididymides of yak and cattleyak and predict the etiology of male infertility in cattleyak.Methods: We performed comparative transcriptome profiles analysis by mRNA sequencing in the epidi-dymides of yak and cattleyak.

Results

In total 3008 differentially expressed genes (DEGs) were identified in cattleyak, out of which 1645 DEGs were up-regulated and 1363 DEGs were down-regulated. Thirteen DEGs were validated by quantitative real-time PCR. DEGs included certain genes that were associated with spermatozoal matura-tion, motility, male fertility, water and ion channels, and beta-defensins. LCN9, SPINT4, CES5A, CD52, CST11, SERPINA1, CTSK, FABP4, CCR5, GRIA2, ENTPD3, LOC523530 and DEFB129, DEFB128, DEFB127, DEFB126, DEFB124, DEFB122A, DEFB122, DEFB119 were all downregu-lated, whereas NRIP1 and TMEM212 among top 30 DEGs were upregulated. Furthermore, protein processing in endoplasmic reticulum pathway was ranked at top-listed three significantly enriched KEGG pathways that as a consequence of abnormal expression of ER-associated genes in the entire ER protein processing pathway might have been disrupted in male cattleyak which resulted in the down-regulation of several important genes. All the DEGs enriched in this pathway were downregulated ex-cept NEF.

Conclusion

Taken together, our findings revealed that there were marked differences in the epididymal transcriptomic profiles of yak and cattleyak. The DEGs were involved in spermatozoal maturation, mo-tility, male fertility, water and ion channels, and beta-defensins. Abnormal expression of ER-associated genes in the entire ER protein processing pathway may have disrupted protein processing pathway in male cattleyak resulting in the downregulation of several important genes involved in sperm maturation, motility and defense.

Metabolomics Integrated with Transcriptomics Reveals Redirection of the Phenylpropanoids Metabolic Flux in Ginkgo biloba

Ginkgo biloba is a monotypic species native to China with great economic and ecological values. Leaves extract of this tree contains about 24% flavonoids, which are widely used in the pharmaceutical industry. However, the flavonoids biosynthesis pathway is poorly understood in Ginkgo. In this study, we comprehensively compared the transcriptome and metabolite profiles of Ginkgo high-flavonoids mutant (ZY1) and Anlu1 (control) leaves. A total of 122 significantly changed metabolites and 1683 differentially expressed genes (DEGs), including 45 transcription factors, were identified in ZY1 compared to those in Anlu1. An integrated analysis of metabolic and transcriptomic data revealed that the abundances of some major flavonoids (especially flavone and flavonol) were most significantly increased, while other phenylpropanoid-derived products and lipids showed the most largely reduced abundances in ZY1 compared to those in Anlu1. Quantitative real-time polymerase chain reaction results confirmed the alterations in the expression levels of genes encoding components of pathways involved in phenylpropanoids and lipids. The redirection of metabolic flux may contribute to increased accumulation of flavonoid levels in ZY1 leaves. Our results provide valuable information for metabolic engineering of Ginkgo flavonoids biosynthesis.

Comprehensive transcriptome analysis and flavonoid profiling of Ginkgo leaves reveals flavonoid content alterations in day-night cycles

Ginkgo leaves are raw materials for flavonoid extraction. Thus, the timing of their harvest is important to optimize the extraction efficiency, which benefits the pharmaceutical industry. In this research, we compared the transcriptomes of Ginkgo leaves harvested at midday and midnight. The differentially expressed genes with the highest probabilities in each step of flavonoid biosynthesis were down-regulated at midnight. Furthermore, real-time PCR corroborated the transcriptome results, indicating the decrease in flavonoid biosynthesis at midnight. The flavonoid profiles of Ginkgo leaves harvested at midday and midnight were compared, and the total flavonoid content decreased at midnight. A detailed analysis of individual flavonoids showed that most of their contents were decreased by various degrees. Our results indicated that circadian rhythms affected the flavonoid contents in Ginkgo leaves, which provides valuable information for optimizing their harvesting times to benefit the pharmaceutical industry.