Comparative transcriptomics method to infer gene coexpression networks and its applications to maize and rice leaf transcriptomes

Multi-omics analyzes of Rosa gigantea illuminate tea scent biosynthesis and release mechanisms

Rose is an important ornamental crop cultivated globally for perfume production. However, our understanding of the mechanisms underlying scent production and molecular breeding for fragrance is hindered by the lack of a reference genome for tea roses. We present the first complete telomere-to-telomere (T2T) genome of Rosa gigantea, with high quality (QV > 60), including detailed characterization of the structural features of repetitive regions. The expansion of genes associated with phenylpropanoid biosynthesis may account for the unique tea scent. We uncover the release rhythm of aromatic volatile organic compounds and their gene regulatory networks through comparative genomics and time-ordered gene co-expression networks. Analyzes of eugenol homologs demonstrate how plants attract pollinators using specialized phenylpropanoids in specific tissues. This study highlights the conservation and utilization of genetic diversity from wild endangered species through multi-omics approaches, providing a scientific foundation for enhancing rose fragrance via de novo domestication.

Deciphering the Role of SVP-Like Genes and Their Key Regulation Networks During Reproductive Cone Development in Pinus tabuliformis

Reproductive development plays an essential role in the perpetuation of genetic material and environmental adaptation. In angiosperms, the Short Vegetative Phase (SVP) serves as a flowering repressor, influencing the development of floral organs. In this study, heterologous transformation of Arabidopsis thaliana with SVP-like genes (PtSVL1 and PtSVL2) derived from Pinus tabuliformis significantly impacted stamen formation and pollen fertility, without altering flowering time. Gene co-expression networks revealed that SVP-like and SOC1-like genes function as key coregulatory transcription factors during the initial stages of cone development in P. tabuliformis. Interestingly, the regulatory module of SOC1 regulated by SVP in angiosperms is absent in conifers and conifer SVP-like exercises its function in a form that is physically bound to SOC1-like. Furthermore, combining the yeast one-hybrid scanning with co-expression network analysis, revealed that SPLs and TPSs were the principal downstream target genes of PtSVL1. Notably, the PtSPL16 promoter is positively regulated by PtSVL1, and overexpression of PtSPL16 results in delayed flowering in Arabidopsis, suggesting that the PtSVL1-PtSPL16 module plays a crucial role in regulating reproductive development in conifers. Collectively, these findings enhance our understanding of the roles of SVP-like genes in conifers and the key regulatory networks centred on PtSVL1 during reproductive cone development.

Identification of key genes involved in lignin and flavonoid accumulation during Tilia tuan seed maturation

KEY MESSAGE

Transcriptomics and phenotypic data analysis identified 24 transcription factors (TFs) that play key roles in regulating the competitive accumulation of lignin and flavonoids. Tilia tuan Szyszyl. (T. tuan) is a timber tree species with important ecological and commercial value. However, its highly lignified pericarp results in a low seed germination rate and a long dormancy period. In addition, it is unknown whether there is an interaction between the biosynthesis of flavonoids and lignin as products of the phenylpropanoid pathway during seed development. To explore the molecular regulatory mechanism of lignin and flavonoid biosynthesis, T. tuan seeds were harvested at five stages (30, 60, 90, 120, and 150 days after pollination) for lignin and flavonoid analyses. The results showed that lignin accumulated rapidly in the early and middle stages (S1, S3, and S4), and rapid accumulation of flavonoids during the early and late stages (S1 and S5). High-throughput RNA sequencing analysis of developing seeds identified 50,553 transcripts, including 223 phenylpropanoid biosynthetic pathway genes involved in lignin accumulation grouped into 3 clusters, and 106 flavonoid biosynthetic pathway genes (FBPGs) grouped into 2 clusters. Subsequent WGCNA and time-ordered gene co-expression network (TO-GCN) analysis revealed that 24 TFs (e.g., TtARF2 and TtWRKY15) were involved in flavonoids and lignin biosynthesis regulation. The transcriptome data were validated by qRT-PCR to analyze the expression profiles of key enzyme-coding genes. This study revealed that there existed a competitive relationship between flavonoid and lignin biosynthesis pathway during the development of T. tuan seeds, that provide a foundation for the further exploration of molecular mechanisms underlying lignin and flavonoid accumulation in T. tuan seeds.

Single-cell RNA-sequencing provides new insights into the cell-specific expression patterns and transcriptional regulation of photosynthetic genes in bermudagrass leaf blades

As an important warm-season turfgrass species, bermudagrass (Cynodon dactylon L.) flourishes in warm areas around the world due to the existence of the C4 photosynthetic pathway. However, how C4 photosynthesis operates in bermudagrass leaves is still poorly understood. In this study, we performed single-cell RNA-sequencing on 5296 cells from bermudagrass leaf blades. Eight cell clusters corresponding to mesophyll, bundle sheath, epidermis and vascular bundle cells were successfully identified using known cell marker genes. Expression profiling indicated that genes encoding NADP-dependent malic enzymes (NADP-MEs) were highly expressed in bundle sheath cells, whereas NAD-ME genes were weakly expressed in all cell types, suggesting C4 photosynthesis of bermudagrass leaf blades might be NADP-ME type rather than NAD-ME type. The results also indicated that starch synthesis-related genes showed preferential expression in bundle sheath cells, whereas starch degradation-related genes were highly expressed in mesophyll cells, which agrees with the observed accumulation of starch-filled chloroplasts in bundle sheath cells. Gene co-expression analysis further revealed that different families of transcription factors were co-expressed with multiple C4 photosynthesis-related genes, suggesting a complex transcription regulatory network of C4 photosynthesis might exist in bermudagrass leaf blades. These findings collectively provided new insights into the cell-specific expression patterns and transcriptional regulation of photosynthetic genes in bermudagrass.

Haplotype-resolved chromosomal-level genome assembly reveals regulatory variations in mulberry fruit anthocyanin content

Understanding the intricate regulatory mechanisms underlying the anthocyanin content (AC) in fruits and vegetables is crucial for advanced biotechnological customization. In this study, we generated high-quality haplotype-resolved genome assemblies for two mulberry cultivars: the high-AC 'Zhongsang5801' (ZS5801) and the low-AC 'Zhenzhubai' (ZZB). Additionally, we conducted a comprehensive analysis of genes associated with AC production. Through genome-wide association studies (GWAS) on 112 mulberry fruits, we identified MaVHAG3, which encodes a vacuolar-type H+-ATPase G3 subunit, as a key gene linked to purple pigmentation. To gain deeper insights into the genetic and molecular processes underlying high AC, we compared the genomes of ZS5801 and ZZB, along with fruit transcriptome data across five developmental stages, and quantified the accumulation of metabolic substances. Compared to ZZB, ZS5801 exhibited significantly more differentially expressed genes (DEGs) related to anthocyanin metabolism and higher levels of anthocyanins and flavonoids. Comparative analyses revealed expansions and contractions in the flavonol synthase (FLS) and dihydroflavonol 4-reductase (DFR) genes, resulting in altered carbon flow. Co-expression analysis demonstrated that ZS5801 displayed more significant alterations in genes involved in late-stage AC regulation compared to ZZB, particularly during the phase stage. In summary, our findings provide valuable insights into the regulation of mulberry fruit AC, offering genetic resources to enhance cultivars with higher AC traits.

Conserved regulatory switches for the transition from natal down to juvenile feather in birds

The transition from natal downs for heat conservation to juvenile feathers for simple flight is a remarkable environmental adaptation process in avian evolution. However, the underlying epigenetic mechanism for this primary feather transition is mostly unknown. Here we conducted time-ordered gene co-expression network construction, epigenetic analysis, and functional perturbations in developing feather follicles to elucidate four downy-juvenile feather transition events. We report that extracellular matrix reorganization leads to peripheral pulp formation, which mediates epithelial-mesenchymal interactions for branching morphogenesis. α-SMA (ACTA2) compartmentalizes dermal papilla stem cells for feather renewal cycling. LEF1 works as a key hub of Wnt signaling to build rachis and converts radial downy to bilateral symmetry. Novel usage of scale keratins strengthens feather sheath with SOX14 as the epigenetic regulator. We show that this primary feather transition is largely conserved in chicken (precocial) and zebra finch (altricial) and discuss the possibility that this evolutionary adaptation process started in feathered dinosaurs.

Transcriptional Regulation of SugarCane Response to Sporisorium scitamineum: Insights from Time-Course Gene Coexpression and Ca2+ Signaling

Sugarcane response to Sporisorium scitamineum is determined by multiple major genes and numerous microeffector genes. Here, time-ordered gene coexpression networks were applied to explore the interaction between sugarcane and S. scitamineum. Totally, 2459 differentially expressed genes were identified and divided into 10 levels, and several stress-related subnetworks were established. Interestingly, the Ca2+ signaling pathway was activated to establish the response to sugarcane smut disease. Accordingly, two CAX genes (ScCAX2 and ScCAX3) were cloned and characterized from sugarcane. They were significantly upregulated under ABA stress but inhibited by MeJA treatment. Furthermore, overexpression of ScCAX2 and ScCAX3 enhanced the susceptibility of transgenic plants to the pathogen infection, suggesting its negative role in disease resistance. A regulatory model for ScCAX genes in disease response was thus depicted. This work helps to clarify the transcriptional regulation of sugarcane response to S. scitamineum stress and the function of the CAX gene in disease response.

Temporal dynamics of genetic architecture governing leaf development in Populus

Leaf development is a multifaceted and dynamic process orchestrated by a myriad of genes to shape the proper size and morphology. The dynamic genetic network underlying leaf development remains largely unknown. Utilizing a synergistic genetic approach encompassing dynamic genome-wide association study (GWAS), time-ordered gene co-expression network (TO-GCN) analyses and gene manipulation, we explored the temporal genetic architecture and regulatory network governing leaf development in Populus. We identified 42 time-specific and 18 consecutive genes that displayed different patterns of expression at various time points. We then constructed eight TO-GCNs that covered the cell proliferation, transition, and cell expansion stages of leaf development. Integrating GWAS and TO-GCN, we postulated the functions of 27 causative genes for GWAS and identified PtoGRF9 as a key player in leaf development. Genetic manipulation via overexpression and suppression of PtoGRF9 revealed its primary influence on leaf development by modulating cell proliferation. Furthermore, we elucidated that PtoGRF9 governs leaf development by activating PtoHB21 during the cell proliferation stage and attenuating PtoLD during the transition stage. Our study provides insights into the dynamic genetic underpinnings of leaf development and understanding the regulatory mechanism of PtoGRF9 in this dynamic process.

Multi-omics identification of a key glycosyl hydrolase gene FtGH1 involved in rutin hydrolysis in Tartary buckwheat (Fagopyrum tataricum)

Rutin, a flavonoid rich in buckwheat, is important for human health and plant resistance to external stresses. The hydrolysis of rutin to quercetin underlies the bitter taste of Tartary buckwheat. In order to identify rutin hydrolysis genes, a 200 genotypes mini-core Tartary buckwheat germplasm resource was re-sequenced with 30-fold coverage depth. By combining the content of the intermediate metabolites of rutin metabolism with genome resequencing data, metabolite genome-wide association analyses (GWAS) eventually identified a glycosyl hydrolase gene FtGH1, which could hydrolyse rutin to quercetin. This function was validated both in Tartary buckwheat overexpression hairy roots and in vitro enzyme activity assays. Mutation of the two key active sites, which were determined by molecular docking and experimentally verified via overexpression in hairy roots and transient expression in tobacco leaves, exhibited abnormal subcellular localization, suggesting functional changes. Sequence analysis revealed that mutation of the FtGH1 promoter in accessions of two haplotypes might be necessary for enzymatic activity. Co-expression analysis and GWAS revealed that FtbHLH165 not only repressed FtGH1 expression, but also increased seed length. This work reveals a potential mechanism behind rutin metabolism, which should provide both theoretical support in the study of flavonoid metabolism and in the molecular breeding of Tartary buckwheat.

Exploration the interaction of cadmium and copper toxic effects in pakchoi (Brassica chinensis L) roots through combinatorial transcriptomic and weighted gene co-expression network analysis

The interaction between cadmium(Cd) and copper(Cu) during combined pollution can lead to more complex toxic effects on humans and plants.However, there is still a lack of sufficient understanding regarding the types of interactions at the plant molecular level and the response strategies of plants to combined pollution. To assess this, we investigated the phenotypic and transcriptomic patterns of pakchoi (Brassica chinensis L) roots in response to individual and combined pollution of Cd and Cu. The results showed that compared to single addition, the translocation factor of heavy metals in roots significantly decreased (p < 0.05) under the combined addition, resulting in higher accumulation of Cd and Cu in the roots. Transcriptomic analysis of pakchoi roots revealed that compared to single pollution, there were 312 and 1926 differentially expressed genes (DEGs) specifically regulated in the Cd2Cu20 and Cd2Cu100 combined treatments, respectively. By comparing the expression of these DEGs among different treatments, we found that the combined pollution of Cd and Cu mainly affected the transcriptome of the roots in an antagonistic manner. Enrichment analysis indicated that pakchoi roots upregulated the expression of genes involved in glucosetransferase activity, phospholipid homeostasis, proton transport, and the biosynthesis of phenylpropanoids and flavonoids to resist Cd and Cu combined pollution. Using weighted gene co-expression network analysis (WGCNA), we identified hub genes related to the accumulation of Cd and Cu in the roots, which mainly belonged to the LBD, thaumatin-like protein, ERF, MYB, WRKY, and TCP transcription factor families. This may reflect a transcription factor-driven trade-off strategy between heavy metal accumulation and growth in pakchoi roots. Additionally, compared to single metal pollution, the expression of genes related to Nramp, cation/H+ antiporters, and some belonging to the ABC transporter family in the pakchoi roots was significantly upregulated under combined pollution. This could lead to increased accumulation of Cd and Cu in the roots. These findings provide new insights into the interactions and toxic mechanisms of multiple metal combined pollution at the molecular level in plants.

Comparative analysis of differential gene expression reveals novel insights into the heteroblastic foliage functional traits of Pinus massoniana seedlings

Pinus massoniana needles, rich in medicinal polysaccharides and flavonoids, undergo heteroblastic foliage, transitioning from primary needles (PN) to secondary needles (SN) during growth, resulting in altered functional traits. Despite its significance, the molecular regulatory mechanisms governing these traits remain unclear. This study employs Iso-Seq and RNA-Seq analyses to explore differentially expressed genes (DEGs) associated with functional traits throughout the main growth season of heteroblastic foliage. Co-expression network analysis identified 34 hub genes and 17 key transcription factors (TFs) influencing light-harvesting antenna, photosystem I and II, crucial in photosynthesis regulation. Additionally, 14 genes involved in polysaccharide metabolism pathways, synthesizing sucrose, glucose, UDP sugars, and xylan, along with four genes in flavonoid biosynthesis pathways, regulating p-coumaroyl-CoA, quercetin, galangin, and myricetin production, exhibited differential expression between PN and SN. Further analysis unveils a highly interconnected network among these genes, forming a pivotal cascade of TFs and DEGs. Therefore, heteroblastic changes significantly impact needle functional traits, potentially affecting the pharmacological properties of PN and SN. Thus, these genomic insights into understanding the molecular-level differences of heteroblastic foliage, thereby establishing a foundation for advancements in the pharmaceutical industry related to needle-derived products.

Drying without dying: A genome database for desiccation-tolerant plants and evolution of desiccation tolerance

Desiccation is typically fatal, but a small number of land plants have evolved vegetative desiccation tolerance (VDT), allowing them to dry without dying through a process called anhydrobiosis. Advances in sequencing technologies have enabled the investigation of genomes for desiccation-tolerant plants over the past decade. However, a dedicated and integrated database for these valuable genomic resources has been lacking. Our prolonged interest in VDT plant genomes motivated us to create the "Drying without Dying" database, which contains a total of 16 VDT-related plant genomes (including 10 mosses) and incorporates 10 genomes that are closely related to VDT plants. The database features bioinformatic tools, such as blast and homologous cluster search, sequence retrieval, Gene Ontology term and metabolic pathway enrichment statistics, expression profiling, co-expression network extraction, and JBrowser exploration for each genome. To demonstrate its utility, we conducted tailored PFAM family statistical analyses, and we discovered that the drought-responsive ABA transporter AWPM-19 family is significantly tandemly duplicated in all bryophytes but rarely so in tracheophytes. Transcriptomic investigations also revealed that response patterns following desiccation diverged between bryophytes and angiosperms. Combined, the analyses provided genomic and transcriptomic evidence supporting a possible divergence and lineage-specific evolution of VDT in plants. The database can be accessed at http://desiccation.novogene.com. We expect this initial release of the "Drying without Dying" plant genome database will facilitate future discovery of VDT genetic resources.

Differential gene expression and potential regulatory network of fatty acid biosynthesis during fruit and leaf development in yellowhorn (Xanthoceras sorbifolium), an oil-producing tree with significant deployment values

Xanthoceras sorbifolium (yellowhorn) is a woody oil plant with super stress resistance and excellent oil characteristics. The yellowhorn oil can be used as biofuel and edible oil with high nutritional and medicinal value. However, genetic studies on yellowhorn are just in the beginning, and fundamental biological questions regarding its very long-chain fatty acid (VLCFA) biosynthesis pathway remain largely unknown. In this study, we reconstructed the VLCFA biosynthesis pathway and annotated 137 genes encoding relevant enzymes. We identified four oleosin genes that package triacylglycerols (TAGs) and are specifically expressed in fruits, likely playing key roles in yellowhorn oil production. Especially, by examining time-ordered gene co-expression network (TO-GCN) constructed from fruit and leaf developments, we identified key enzymatic genes and potential regulatory transcription factors involved in VLCFA synthesis. In fruits, we further inferred a hierarchical regulatory network with MYB-related (XS03G0296800) and B3 (XS02G0057600) transcription factors as top-tier regulators, providing clues into factors controlling carbon flux into fatty acids. Our results offer new insights into key genes and transcriptional regulators governing fatty acid production in yellowhorn, laying the foundation for efforts to optimize oil content and fatty acid composition. Moreover, the gene expression patterns and putative regulatory relationships identified here will inform metabolic engineering and molecular breeding approaches tailored to meet biofuel and bioproduct demands.

How to utilize comparative transcriptomics to dissect morphological diversity in plants

Comparative transcriptomics has emerged as a powerful approach that allows us to unravel the genetic basis of organ morphogenesis and its diversification processes during evolution. However, the application of comparative transcriptomics in studying plant morphological diversity addresses challenges such as identifying homologous gene pairs, selecting appropriate developmental stages for comparison, and extracting biologically meaningful networks. Methods such as phylostratigraphy, clustering, and gene co-expression networks are explored to identify functionally equivalent genes, align developmental stages, and uncover gene regulatory relationships. In the current review, we highlight the importance of these approaches in overcoming the complexity of plant genomes, the impact of heterochrony on stage alignment, and the integration of gene networks with additional data for a comprehensive understanding of morphological evolution.

C4 leaf development and evolution

C4 photosynthesis is more efficient than C3 photosynthesis for two reasons. First, C4 plants have evolved efficient C4 enzymes to suppress wasteful photorespiration and enhance CO2 fixation. Second, C4 leaves have Kranz anatomy in which the veins are surrounded by one layer of bundle sheath (BS) cells and one layer of mesophyll (M) cells. The BS and M cells are functionally well differentiated and also well coordinated for rapid assimilation of atmospheric CO2 and transport of photo-assimilates between the two types of cells. Recent comparative transcriptomics of developing M and BS cells in young maize embryonic leaves revealed not only potential regulators of BS and M cell differentiation but also rapid early BS cell differentiation whereas slower, more prolonged M cell differentiation, contrary to the traditional view of a far simpler process of M cell development. Moreover, new upstream regulators of Kranz anatomy development have been identified and a number of gene co-expression modules for early vascular development have been inferred. Also, a candidate gene regulatory network associated with Kranz anatomy and vascular development has been constructed. Additionally, how whole genome duplication (WGD) may facilitate C4 evolution has been studied and the reasons for why the same WGD event led to successful C4 evolution in Gynandropsis gynandra but not in the sister species Tarenaya hassleriana have been proposed. Finally, new future research directions are suggested.

Epigenetic regulator RNF20 underlies temporal hierarchy of gene expression to regulate postnatal cardiomyocyte polarization

Differentiated cardiomyocytes (CMs) must undergo diverse morphological and functional changes during postnatal development. However, the mechanisms underlying initiation and coordination of these changes remain unclear. Here, we delineate an integrated, time-ordered transcriptional network that begins with expression of genes for cell-cell connections and leads to a sequence of structural, cell-cycle, functional, and metabolic transitions in mouse postnatal hearts. Depletion of histone H2B ubiquitin ligase RNF20 disrupts this gene network and impairs CM polarization. Subsequently, assay for transposase-accessible chromatin using sequencing (ATAC-seq) analysis confirmed that RNF20 contributes to chromatin accessibility in this context. As such, RNF20 is likely to facilitate binding of transcription factors at the promoters of genes involved in cell-cell connections and actin organization, which are crucial for CM polarization and functional integration. These results suggest that CM polarization is one of the earliest events during postnatal heart development and provide insights into how RNF20 regulates CM polarity and the postnatal gene program.

Comparative transcriptome analysis of dioecious floral development in Trachycarpus fortunei using Illumina and PacBio SMRT sequencing

BACKGROUND

Trachycarpus fortunei is a plant with significant economic and ornamental value. Both male and female flowers of T. fortunei originate as bisexual flowers, and selective abortion occurs during floral development. However, the regulatory mechanisms underlying this process remain unclear in T. fortunei. In this study, transcriptome sequencing with Illumina and Pacific BioSciences (PacBio) single-molecule real-time (SMRT) platforms were used to investigate gene expression differences between male and female T. fortunei plants.

RESULTS

A total of 833,137 full-length non-chimeric (FLNC) reads were obtained, and 726,846 high-quality full-length transcripts were identified. A total of 159 genes were differentially expressed between male and female flowers at all development stages. Some of the differentially expressed genes (DEGs) showed male bias, including serine/threonine-protein kinase (STPK), THUMP1 homolog and other genes. Through single-nucleotide polymorphisms(SNPs) identification, 28 genes were considered as potential sex-associated SNPs. Time-Ordered Gene Co-expression Network (TO-GCN) analysis revealed that MADS2 and MADS26 may play important roles in the development of female and male flowers T. fortune plants, respectively.

CONCLUSIONS

These findings provide a genetic basis for flower development and differentiation in T. fortunei, and improve our understanding of the mechanisms underlying sexual differentiation in T. fortunei.

Rice GA3ox1 modulates pollen starch granule accumulation and pollen wall development

The rice GA biosynthetic gene OsGA3ox1 has been proposed to regulate pollen development through the gametophytic manner, but cellular characterization of its mutant pollen is lacking. In this study, three heterozygotic biallelic variants, "-3/-19", "-3/-2" and "-3/-10", each containing one null and one 3bp-deletion allele, were obtained by the CRISPR/Cas9 technique for the functional study of OsGA3ox1. The three homozygotes, "-19/-19", "-2/-2" and "-10/-10", derived from heterozygotic variants, did not affect the development of most vegetative and floral organs but showed a significant reduction in seed-setting rate and in pollen viability. Anatomic characterizations of these mutated osga3ox1 pollens revealed defects in starch granule accumulation and pollen wall development. Additional molecular characterization suggests that abnormal pollen development in the osga3ox1 mutants might be linked to the regulation of transcription factors OsGAMYB, OsTDR and OsbHLH142 during late pollen development. In brief, the rice GA3ox1 is a crucial gene that modulates pollen starch granule accumulation and pollen wall development at the gametophytic phase.

Transition from natal downs to juvenile feathers: conserved regulatory switches in Neoaves

The transition from natal downs for heat conservation to juvenile feathers for simple flight is a remarkable environmental adaptation process in avian evolution. However, the underlying epigenetic mechanism for this primary feather transition is mostly unknown. Here we conducted time-ordered gene co-expression network construction, epigenetic analysis, and functional perturbations in developing feather follicles to elucidate four downy-juvenile feather transition events. We discovered that LEF1 works as a key hub of Wnt signaling to build rachis and converts radial downy to bilateral symmetry. Extracellular matrix reorganization leads to peripheral pulp formation, which mediates epithelial -mesenchymal interactions for branching morphogenesis. ACTA2 compartments dermal papilla stem cells for feather cycling. Novel usage of scale keratins strengthens feather sheath with SOX14 as the epigenetic regulator. We found this primary feather transition largely conserved in chicken (precocious) and zebra finch (altricial) and discussed the possibility that this evolutionary adaptation process started in feathered dinosaurs.

Diel dynamics of multi-omics in elkhorn fern provide new insights into weak CAM photosynthesis

Crassulacean acid metabolism (CAM) has high water-use efficiency (WUE) and is widely recognized to have evolved from C3 photosynthesis. Different plant lineages have convergently evolved CAM, but the molecular mechanism that underlies C3-to-CAM evolution remains to be clarified. Platycerium bifurcatum (elkhorn fern) provides an opportunity to study the molecular changes underlying the transition from C3 to CAM photosynthesis because both modes of photosynthesis occur in this species, with sporotrophophyll leaves (SLs) and cover leaves (CLs) performing C3 and weak CAM photosynthesis, respectively. Here, we report that the physiological and biochemical attributes of CAM in weak CAM-performing CLs differed from those in strong CAM species. We investigated the diel dynamics of the metabolome, proteome, and transcriptome in these dimorphic leaves within the same genetic background and under identical environmental conditions. We found that multi-omic diel dynamics in P. bifurcatum exhibit both tissue and diel effects. Our analysis revealed temporal rewiring of biochemistry relevant to the energy-producing pathway (TCA cycle), CAM pathway, and stomatal movement in CLs compared with SLs. We also confirmed that PHOSPHOENOLPYRUVATE CARBOXYLASE KINASE (PPCK) exhibits convergence in gene expression among highly divergent CAM lineages. Gene regulatory network analysis identified candidate transcription factors regulating the CAM pathway and stomatal movement. Taken together, our results provide new insights into weak CAM photosynthesis and new avenues for CAM bioengineering.

Brassinosteroids affect wood development and properties of Fraxinus mandshurica

Introduction

Xylem development plays a crucial role in wood formation in woody plants. In recent years, there has been growing attention towards the impact of brassinosteroids (BRs) on this xylem development. In the present study, we evaluated the dynamic variation of xylem development in Fraxinus mandshurica (female parent, M8) and a novel interspecific hybrid F. mandshurica × Fraxinus sogdiana (1601) from May to August 2020.

Methods

We obtained RNA-Seq transcriptomes of three tissue types (xylem, phloem, and leaf) to identify the differences in xylem-differentially expressed genes (X-DEGs) and xylem-specifically expressed genes (X-SEGs) in M8 and 1601 variants. We then further evaluated these genes via weighted gene co-expression network analysis (WGCNA) alongside overexpressing FmCPD, a BR biosynthesis enzyme gene, in transient transgenic F. mandshurica.

Results

Our results indicated that the xylem development cycle of 1601 was extended by 2 weeks compared to that of M8. In addition, during the later wood development stages (secondary wall thickening) of 1601, an increased cellulose content (14%) and a reduced lignin content (11%) was observed. Furthermore, vessel length and width increased by 67% and 37%, respectively, in 1601 compared with those of M8. A total of 4589 X-DEGs were identified, including enzymes related to phenylpropane metabolism, galactose metabolism, BR synthesis, and signal transduction pathways. WGCNA identified hub X-SEGs involved in cellulose synthesis and BR signaling in the 1601 wood formation-related module (CESA8, COR1, C3H14, and C3H15); in contrast, genes involved in phenylpropane metabolism were significantly enriched in the M8 wood formation-related module (CCoAOMT and CCR). Moreover, overexpression of FmCPD in transient transgenic F. mandshurica affected the expression of genes associated with lignin and cellulose biosynthesis signal transduction. Finally, BR content was determined to be approximately 20% lower in the M8 xylem than in the 1601 xylem, and the exogenous application of BRs (24-epi brassinolide) significantly increased the number of xylem cell layers and altered the composition of the secondary cell walls in F. mandshurica.

Discussion

Our findings suggest that BR biosynthesis and signaling play a critical role in the differing wood development and properties observed between M8 and 1601 F. mandshurica.

Time-series transcriptome provides insights into the gene regulation network involved in the icariin-flavonoid metabolism during the leaf development of Epimedium pubescens

Herba Epimedii (Epimedium) leaves are rich in prenylated flavonol glycosides (PFGs) with high medicinal value. However, the dynamics and regulatory network of PFG biosynthesis remain largely unclear. Here, we combined metabolite profiling (targeted to PFGs) and a high-temporal-resolution transcriptome to elucidate PFGs' regulatory network in Epimedium pubescens and identified key candidate structural genes and transcription factors (TFs) involved in PFG accumulation. Chemical profile analysis revealed that PFG content was quite different between buds and leaves and displayed a continuous decline with leaf development. The structural genes are the determinant reasons, and they are strictly regulated by TFs under temporal cues. We further constructed seven time-ordered gene co-expression networks (TO-GCNs) of PFG biosynthesis genes (including EpPAL2, EpC4H, EpCHS2, EpCHI2, EpF3H, EpFLS3, and EpPT8), and three flavonol biosynthesis routines were then predicted. The TFs involved in TO-GCNs were further confirmed by WGCNA analysis. Fourteen hub genes, comprising 5 MYBs, 1 bHLH, 1 WD40, 2 bZIPs, 1 BES1, 1 C2H2, 1 Trihelix, 1 HD-ZIP, and 1 GATA were identified as candidate key TFs. The results were further validated by TF binding site (TFBS) analysis and qRT-PCR. Overall, these findings provide valuable information for understanding the molecular regulatory mechanism of PFGs biosynthesis, enriching the gene resources, which will guide further research on PFG accumulation in Epimedium.

Weighted gene co-expression network analysis of nitrogen (N)-responsive genes and the putative role of G-quadruplexes in N use efficiency (NUE) in rice

Rice is an important target to improve crop nitrogen (N) use efficiency (NUE), and the identification and shortlisting of the candidate genes are still in progress. We analyzed data from 16 published N-responsive transcriptomes/microarrays to identify, eight datasets that contained the maximum number of 3020 common genes, referred to as N-responsive genes. These include different classes of transcription factors, transporters, miRNA targets, kinases and events of post-translational modifications. A Weighted gene co-expression network analysis (WGCNA) with all the 3020 N-responsive genes revealed 15 co-expression modules and their annotated biological roles. Protein-protein interaction network analysis of the main module revealed the hub genes and their functional annotation revealed their involvement in the ubiquitin process. Further, the occurrences of G-quadruplex sequences were examined, which are known to play important roles in epigenetic regulation but are hitherto unknown in N-response/NUE. Out of the 3020 N-responsive genes studied, 2298 contained G-quadruplex sequences. We compared these N-responsive genes containing G-quadruplex sequences with the 3601 genes we previously identified as NUE-related (for being both N-responsive and yield-associated). This analysis revealed 389 (17%) NUE-related genes containing G-quadruplex sequences. These genes may be involved in the epigenetic regulation of NUE, while the rest of the 83% (1811) genes may regulate NUE through genetic mechanisms and/or other epigenetic means besides G-quadruplexes. A few potentially important genes/processes identified as associated with NUE were experimentally validated in a pair of rice genotypes contrasting for NUE. The results from the WGCNA and G4 sequence analysis of N-responsive genes helped identify and shortlist six genes as candidates to improve NUE. Further, the hitherto unavailable segregation of genetic and epigenetic gene targets could aid in informed interventions through genetic and epigenetic means of crop improvement.

Transcriptome dynamics provides insights into divergences of the photosynthesis pathway between Saccharum officinarum and Saccharum spontaneum

Saccharum spontaneum and Saccharum officinarum contributed to the genetic background of modern sugarcane cultivars. Saccharum spontaneum has shown a higher net photosynthetic rate and lower soluble sugar than S. officinarum. Here, we analyzed 198 RNA-sequencing samples to investigate the molecular mechanisms for the divergences of photosynthesis and sugar accumulation between the two Saccharum species. We constructed gene co-expression networks based on differentially expressed genes (DEGs) both for leaf developmental gradients and diurnal rhythm. Our results suggested that the divergence of sugar accumulation may be attributed to the enrichment of major carbohydrate metabolism and the oxidative pentose phosphate pathway. Compared with S. officinarum, S. spontaneum DEGs showed a high enrichment of photosynthesis and contained more complex regulation of photosynthesis-related genes. Noticeably, S. spontaneum lacked gene interactions with sulfur assimilation stimulated by photorespiration. In S. spontaneum, core genes related to clock and photorespiration displayed a sensitive regulation by the diurnal rhythm and phase-shift. Small subunit of Rubisco (RBCS) displayed higher expression in the source tissues of S. spontaneum. Additionally, it was more sensitive under a diurnal rhythm, and had more complex gene networks than that in S. officinarum. This indicates that the differential regulation of RBCS Rubisco contributed to photosynthesis capacity divergence in both Saccharum species.

A functional gene module identification algorithm in gene expression data based on genetic algorithm and gene ontology

Since genes do not function individually, the gene module is considered an important tool for interpreting gene expression profiles. In order to consider both functional similarity and expression similarity in module identification, GMIGAGO, a functional Gene Module Identification algorithm based on Genetic Algorithm and Gene Ontology, was proposed in this work. GMIGAGO is an overlapping gene module identification algorithm, which mainly includes two stages: In the first stage (initial identification of gene modules), Improved Partitioning Around Medoids Based on Genetic Algorithm (PAM-GA) is used for the initial clustering on gene expression profiling, and traditional gene co-expression modules can be obtained. Only similarity of expression levels is considered at this stage. In the second stage (optimization of functional similarity within gene modules), Genetic Algorithm for Functional Similarity Optimization (FSO-GA) is used to optimize gene modules based on gene ontology, and functional similarity within gene modules can be improved. Without loss of generality, we compared GMIGAGO with state-of-the-art gene module identification methods on six gene expression datasets, and GMIGAGO identified the gene modules with the highest functional similarity (much higher than state-of-the-art algorithms). GMIGAGO was applied in BRCA, THCA, HNSC, COVID-19, Stem, and Radiation datasets, and it identified some interesting modules which performed important biological functions. The hub genes in these modules could be used as potential targets for diseases or radiation protection. In summary, GMIGAGO has excellent performance in mining molecular mechanisms, and it can also identify potential biomarkers for individual precision therapy.

Transcriptomic mechanisms of reduced PM2.5-Pb retention in the leaves of the low-Pb-accumulation genotype of Chinese cabbage

Atmospheric fine particulate matter (PM_2.5) mainly contributes to Pb accumulation in the edible leaves of Chinese cabbage in North China. It was found that a low-Pb-accumulation (LPA) genotype of Chinese cabbage contained less Pb in leaves than high-Pb-accumulation (HPA) genotype exposed to PM_2.5-Pb. However, there are no data on the transcriptional regulatory mechanisms of foliar PM_2.5-Pb uptake by Chinese cabbage. The present study investigated the retention of PM_2.5-Pb in foliar apoplast and symplasm and the underlying molecular mechanisms of reduced Pb in LPA leaves. It appeared more Pb in apoplast and less Pb in symplasm of LPA leaves, whereas the pattern was opposite in HPA. There were 2646 and 3095 differentially expressed genes (DEGs) in LPA and HPA leaves under PM_2.5-Pb stress with clearly genotype-specific function, respectively. Furthermore, mRNA levels of XTH16 regulating cell wall thickening, PME2 and PME6 involved in cell wall remodification were significantly expressed in LPA, but not in HPA. Meanwhile, foliar PM_2.5-Pb stress downregulated expression of ZIP1, YSL1, and CNGC3 responsible for Pb influx to cell, and upregulated expression of ABCG36 regulated Pb efflux from symplasm in LPA leaves. These results improve our understanding to the mechanisms underlying foliar Pb uptake from PM_2.5-Pb at transcriptomic level.

Upregulation of PD-L1 by SARS-CoV-2 promotes immune evasion

Patients with severe COVID-19 often suffer from lymphopenia, which is linked to T-cell sequestration, cytokine storm, and mortality. However, it remains largely unknown how severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) induces lymphopenia. Here, we studied the transcriptomic profile and epigenomic alterations involved in cytokine production by SARS-CoV-2-infected cells. We adopted a reverse time-order gene coexpression network approach to analyze time-series RNA-sequencing data, revealing epigenetic modifications at the late stage of viral egress. Furthermore, we identified SARS-CoV-2-activated nuclear factor-κB (NF-κB) and interferon regulatory factor 1 (IRF1) pathways contributing to viral infection and COVID-19 severity through epigenetic analysis of H3K4me3 chromatin immunoprecipitation sequencing. Cross-referencing our transcriptomic and epigenomic data sets revealed that coupling NF-κB and IRF1 pathways mediate programmed death ligand-1 (PD-L1) immunosuppressive programs. Interestingly, we observed higher PD-L1 expression in Omicron-infected cells than SARS-CoV-2 infected cells. Blocking PD-L1 at an early stage of virally-infected AAV-hACE2 mice significantly recovered lymphocyte counts and lowered inflammatory cytokine levels. Our findings indicate that targeting the SARS-CoV-2-mediated NF-κB and IRF1-PD-L1 axis may represent an alternative strategy to reduce COVID-19 severity.

Comparative network analysis reveals the dynamics of organic acid diversity during fruit ripening in peach (Prunus persica L. Batsch)

BACKGROUND

Organic acids are important components that determine the fruit flavor of peach (Prunus persica L. Batsch). However, the dynamics of organic acid diversity during fruit ripening and the key genes that modulate the organic acids metabolism remain largely unknown in this kind of fruit tree which yield ranks sixth in the world.

RESULTS

In this study, we used 3D transcriptome data containing three dimensions of information, namely time, phenotype and gene expression, from 5 different varieties of peach to construct gene co-expression networks throughout fruit ripening of peach. With the network inferred, the time-ordered network comparative analysis was performed to select high-acid specific gene co-expression network and then clarify the regulatory factors controlling organic acid accumulation. As a result, network modules related to organic acid synthesis and metabolism under high-acid and low-acid comparison conditions were identified for our following research. In addition, we obtained 20 candidate genes as regulatory factors related to organic acid metabolism in peach.

CONCLUSIONS

The study provides new insights into the dynamics of organic acid accumulation during fruit ripening, complements the results of classical co-expression network analysis and establishes a foundation for key genes discovery from time-series multiple species transcriptome data.

Gapless genome assembly of azalea and multi-omics investigation into divergence between two species with distinct flower color

The genus Rhododendron (Ericaceae), with more than 1000 species highly diverse in flower color, is providing distinct ornamental values and a model system for flower color studies. Here, we investigated the divergence between two parental species with different flower color widely used for azalea breeding. Gapless genome assembly was generated for the yellow-flowered azalea, Rhododendron molle. Comparative genomics found recent proliferation of long terminal repeat retrotransposons (LTR-RTs), especially Gypsy, has resulted in a 125 Mb (19%) genome size increase in species-specific regions, and a significant amount of dispersed gene duplicates (13 402) and pseudogenes (17 437). Metabolomic assessment revealed that yellow flower coloration is attributed to the dynamic changes of carotenoids/flavonols biosynthesis and chlorophyll degradation. Time-ordered gene co-expression networks (TO-GCNs) and the comparison confirmed the metabolome and uncovered the specific gene regulatory changes underpinning the distinct flower pigmentation. B3 and ERF TFs were found dominating the gene regulation of carotenoids/flavonols characterized pigmentation in R. molle, while WRKY, ERF, WD40, C2H2, and NAC TFs collectively regulated the anthocyanins characterized pigmentation in the red-flowered R simsii. This study employed a multi-omics strategy in disentangling the complex divergence between two important azaleas and provided references for further functional genetics and molecular breeding.

Simultaneous detection of miRNA and mRNA at the single-cell level in plant tissues

Detecting the simultaneous presence of a microRNA (miRNA) and a mRNA in a specific tissue can provide support for the prediction that the miRNA regulates the mRNA. Although two such methods have been developed for mammalian tissues, they have a low signal-noise ratio and/or poor resolution at the single-cell level. To overcome these drawbacks, we develop a method that uses sequence-specific miRNA-locked nucleic acid (LNA) and mRNA-LNA probes. Moreover, it augments the detection signal by rolling circle amplification, achieving a high signal-noise ratio at the single-cell level. Dot signals are counted for determining the expression levels of mRNA and miRNA molecules in specific cells. We show a high sequence specificity of our miRNA-LNA probe, revealing that it can discriminate single-base mismatches. Numerical quantification by our method is tested in transgenic rice lines with different gene expression levels. We conduct several applications. First, the spatial expression profiling of osa-miR156 and OsSPL12 in rice leaves reveals their specific expression in mesophyll cells. Second, studying rice and its mutant lines with our method reveals opposite expression patterns of miRNA and its target mRNA in tissues. Third, the dynamic expression profiles of ZmGRF8 and zma-miR396 during maize leaf development provide evidence that zma-miR396 regulates the preferential spatial expression of ZmGRF8 in bundle sheath cells. Finally, our method can be scaled up to simultaneously detect multiple miRNAs and mRNAs in a tissue. Thus, it is a sensitive and versatile technique for studying miRNA regulation of plant tissue development.

E2Fs co-participate in cadmium stress response through activation of MSHs during the cell cycle

Cadmium is one of the most common heavy metal contaminants found in agricultural fields. MutSα, MutSβ, and MutSγ are three different MutS-associated protein heterodimer complexes consisting of MSH2/MSH6, MSH2/MSH3, and MSH2/MSH7, respectively. These complexes have different mismatch recognition properties and abilities to support MMR. However, changes in mismatch repair genes (OsMSH2, OsMSH3, OsMSH6, and OsMSH7) of the MutS system in rice, one of the most important food crops, under cadmium stress and their association with E2Fs, the key transcription factors affecting cell cycles, are poorly evaluated. In this study, we systematically categorized six rice E2Fs and confirmed that OsMSHs were the downstream target genes of E2F using dual-luciferase reporter assays. In addition, we constructed four msh mutant rice varieties (msh2, msh3, msh6, and msh7) using the CRISPR-Cas9 technology, exposed these mutant rice seedlings to different concentrations of cadmium (0, 2, and 4 mg/L) and observed changes in their phenotype and transcriptomic profiles using RNA-Seq and qRT-PCR. We found that the difference in plant height before and after cadmium stress was more significant in mutant rice seedlings than in wild-type rice seedlings. Transcriptomic profiling and qRT-PCR quantification showed that cadmium stress specifically mobilized cell cycle-related genes ATR, CDKB2;1, MAD2, CycD5;2, CDKA;1, and OsRBR1. Furthermore, we expressed OsE2Fs in yeasts and found that heterologous E2F expression in yeast strains regulated cadmium tolerance by regulating MSHs expression. Further exploration of the underlying mechanisms revealed that cadmium stress may activate the CDKA/CYCD complex, which phosphorylates RBR proteins to release E2F, to regulate downstream MSHs expression and subsequent DNA damage repairment, thereby enhancing the response to cadmium stress.

Effects of ocean acidification and tralopyril on bivalve biomineralization and carbon cycling: A study of the Pacific Oyster (Crassostrea gigas)

The combined effects of emerging pollutants and ocean acidification (OA) on marine organisms and marine ecosystems have attracted increasing attention. However, the combined effects of tralopyril and OA on marine organisms and marine ecosystems remain unclear. In this study, Crassostrea gigas (C. gigas) were exposed to tralopyril (1 μg/L) and/or OA (PH = 7.7) for 21 days and a 14-day recovery acclimation. To investigate the stress response and potential molecular mechanisms of C. gigas to OA and tralopyril exposure alone or in combination, as well as the effects of OA and/or tralopyril on bivalve biomineralization and marine carbon cycling. The results showed that the combined toxicity was between that of acidification and tralopyril alone. Single or combined exposure activated the general stress defense responses of C. gigas mantle, affected energy metabolism and biomineralization of the organism and the carbon cycle of the marine ecosystem. Moreover, acidification-induced and tralopyril-induced toxicity showed potential recoverability at molecular and biochemical levels. This study provides a new perspective on the molecular mechanisms of tralopyril toxicity to bivalve shellfish and reveals the potential role of tralopyril and OA on marine carbon cycling.

Effects of exogenous melatonin on expressional differences of immune-related genes in cashmere goats

The interplay between melatonin and immune system is well recognized in humans. The true integration of research on cashmere goat is still far from clear, especially for cashmere goat maintained in wool and cashmere growth. In this study, we applied various approaches to identify the complex regulated network between the immune-related genes and transcription factors (TFs) and to explore the relationship between melatonin and gene expression in cashmere goats. In total, 1,599 and 1756 immune-related genes were found in the blood and skin of cashmere goats, respectively, and 24 differentially expressed immune-related GO terms were highly expressed in blood after melatonin implantation. We studied the melatonin-dependent networks between the TFs and immune-related genes in cashmere goat. The 3 major regulatory networks were interconnected through TFs. The TFs, such as PHF5A, REXO4, STRAP, JUNB, GATAD2A, ZNF710, and VDR, were also expressed in the blood and skin tissue of cashmere goat. In addition, most genes in these networks, such as VDR, JUNB, and Trib3, were involved in WNT pathway, which is related to cashmere wool growth regulation. On the network basis, we developed a knockout mouse model to identify the network interaction. We observed that 8 high-sulfur protein genes, 12 keratin (KRT) genes, and 19 keratin associated protein (KRTAP) genes related to the growth of cashmere wool were almost not expressed in Trib3^−/− rat skin. Our results suggested that the expression of genes related to wool and cashmere growth may be regulated by the interaction network between genes affected by melatonin and immune-related genes. In summary, we outlined some particularly promising ways for future research on immune-related genes of cashmere goats and the role of melatonin in wool and cashmere growth.

Multi-Omics Approaches and Resources for Systems-Level Gene Function Prediction in the Plant Kingdom

In higher plants, the complexity of a system and the components within and among species are rapidly dissected by omics technologies. Multi-omics datasets are integrated to infer and enable a comprehensive understanding of the life processes of organisms of interest. Further, growing open-source datasets coupled with the emergence of high-performance computing and development of computational tools for biological sciences have assisted in silico functional prediction of unknown genes, proteins and metabolites, otherwise known as uncharacterized. The systems biology approach includes data collection and filtration, system modelling, experimentation and the establishment of new hypotheses for experimental validation. Informatics technologies add meaningful sense to the output generated by complex bioinformatics algorithms, which are now freely available in a user-friendly graphical user interface. These resources accentuate gene function prediction at a relatively minimal cost and effort. Herein, we present a comprehensive view of relevant approaches available for system-level gene function prediction in the plant kingdom. Together, the most recent applications and sought-after principles for gene mining are discussed to benefit the plant research community. A realistic tabulation of plant genomic resources is included for a less laborious and accurate candidate gene discovery in basic plant research and improvement strategies.

Comparative study of disease progression for heart failure with different etiologies via time-ordered network analysis

OBJECTIVES

Heart failure (HF), the primary end-stage manifestation of multiple cardiovascular diseases, has become a global epidemic with high morbidity and mortality. However, the mechanisms underlying the pathogenesis of HF with different etiologies have yet to be fully elucidated.

METHODS

In this study, we developed a novel method to determine the dysregulated lncRNA-mRNA regulation pairs (LMRPs) in the different causes that lead to HF. Time-ordered dysregulated lncRNA-mRNA regulation networks were constructed for comparing the HF progression initiated from different causes. Additionally, the random forest and support vector machine classification algorithm were applied to identify HF-related diagnostic biomarkers.

RESULTS

Biological functional analysis indicated that similar functions were detected at the late stage across different causes of HF, whereas different characteristics were revealed during disease progression. Specifically, the disturbance of myocardial energy metabolism might be a cause of dilated cardiomyopathy (DCM) and peripartum cardiomyopathy (PPCM), while immune response appeared earlier in hypertrophic cardiomyopathy (HCM). Inflammatory response during HCM and PPCM progression might be mediated by complement system, whereas ischemic cardiomyopathy (ICM) might be induced by cytokines. Finally, we identified several panels of diagnostic biomarkers for distinguishing HF patients of different etiologies from non-heart failure (NF) controls.

CONCLUSIONS

This study revealed distinct functional characteristics during the progression of HF from different causes and facilitated the discovery of candidate diagnostic biomarkers for HF.

Generating transcriptional regulatory networks from time-ordered stem cell differentiation RNA sequencing data

We describe steps to 1) identify ascending and descending monotonic key genes from time-ordered stem cell differentiation expression data, 2) construct time-ordered transcriptional regulatory networks, and 3) infer the involvement of transcription factors along the differentiation process.

For complete details on the use and execution of this protocol, please refer to Wong et al. (2020).

•

A protocol to identify ascending and descending monotonic pattern genes

•

Temporal gene regulation in development

•

Pipeline can be adapted to other time-series data sets

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

A comparative transcriptomics and eQTL approach identifies SlWD40 as a tomato fruit ripening regulator

Although multiple vital genes with strong effects on the tomato (Solanum lycopersicum) ripening process have been identified via the positional cloning of ripening mutants and cloning of ripening-related transcription factors (TFs), recent studies suggest that it is unlikely that we have fully characterized the gene regulatory networks underpinning this process. Here, combining comparative transcriptomics and expression QTLs, we identified 16 candidate genes involved in tomato fruit ripening and validated them through virus-induced gene silencing analysis. To further confirm the accuracy of the approach, one potential ripening regulator, SlWD40 (WD-40 repeats), was chosen for in-depth analysis. Co-expression network analysis indicated that master regulators such as RIN (ripening inhibitor) and NOR (nonripening) as well as vital TFs including FUL1 (FRUITFUL1), SlNAC4 (NAM, ATAF1,2, and CUC2 4), and AP2a (Activating enhancer binding Protein 2 alpha) strongly co-expressed with SlWD40. Furthermore, SlWD40 overexpression and RNAi lines exhibited substantially accelerated and delayed ripening phenotypes compared with the wild type, respectively. Moreover, transcriptome analysis of these transgenics revealed that expression patterns of ethylene biosynthesis genes, phytoene synthase, pectate lyase, and branched chain amino transferase 2, in SlWD40-RNAi lines were similar to those of rin and nor fruits, which further demonstrated that SlWD40 may act as an important ripening regulator in conjunction with RIN and NOR. These results are discussed in the context of current models of ripening and in terms of the use of comparative genomics and transcriptomics as an effective route for isolating causal genes underlying differences in genotypes.

UV-B and UV-C radiation trigger both common and distinctive signal perceptions and transmissions in Pinus tabuliformis Carr

In plants, ultraviolet (UV)-light is an important driver for growth and natural distribution, and is also a valuable tool for manipulating productivity as well as biotic interactions. Understanding of plant responses to different UV radiation is sparse, especially from a systems biology perspective and particularly for conifers. Here, we evaluated the physiological and transcriptomic responses to the short-term application of high-irradiance UV-B and UV-C waves on Pinus tabuliformis Carr., a major conifer in Northern China. By undertaking time-ordered gene coexpression network analyses and network comparisons incorporating physiological traits and gene expression variation, we uncovered communalities but also differences in P. tabuliformis responses to UV-B and UV-C. Both types of spectral bands caused a significant inhibition of photosynthesis, and conversely, the improvement of antioxidant capacity, flavonoid production and signaling pathways related to stress resistance, indicating a clear switch from predominantly primary metabolism to enhanced defensive metabolism in pine. We isolated distinct subnetworks for photoreceptor-mediated signal transduction, maximum quantum efficiency of photosystem II (Fv/Fm) regulation and flavonoid biosynthesis in response to UV-B and UV-C radiation. From these subnetworks, we further identified phototropins as potentially important elements in both UV-B and UV-C signaling and, for the first time, suggesting peptide hormones to be involved in promoting flavonoid biosynthesis against UV-B, while these hormones seem not to be implicated in the defense against UV-C exposure. The present study employed an effective strategy for disentangling the complex physiological and genetic regulatory mechanisms in a nonmodel plant species, and thus, provides a suitable reference for future functional evaluations and artificial UV-light mediated growing strategies in plant production.

Lilac (Syringa oblata) genome provides insights into its evolution and molecular mechanism of petal color change

Color change during flower opening is common; however, little is understood on the biochemical and molecular basis related. Lilac (Syringa oblata), a well-known woody ornamental plant with obvious petal color changes, is an ideal model. Here, we presented chromosome-scale genome assembly for lilac, resolved the flavonoids metabolism, and identified key genes and potential regulatory networks related to petal color change. The genome assembly is 1.05 Gb anchored onto 23 chromosomes, with a BUSCO score of 96.6%. Whole-genome duplication (WGD) event shared within Oleaceae was revealed. Metabolome quantification identified delphinidin-3-O-rutinoside (Dp3Ru) and cyanidin-3-O-rutinoside (Cy3Ru) as the major pigments; gene co-expression networks indicated WRKY an essential regulation factor at the early flowering stage, ERF more important in the color transition period (from violet to light nearly white), while the MBW complex participated in the entire process. Our results provide a foundation for functional study and molecular breeding in lilac.

The Chromosome-Scale Assembly of the Curcuma alismatifolia Genome Provides Insight Into Anthocyanin and Terpenoid Biosynthesis

Curcuma alismatifolia, a bulbous flower known for its showy bracts, is widely used around the world as a cut flower, potted, and garden plant. Besides its ornamental value, this species is rich in terpenoid metabolites and could serve as a resource for essential oils. Here, we report a chromosome-level genome assembly of C. alismatifolia and describe its biosynthetic pathways for anthocyanins and terpenoids. This high-quality, assembled genome size is 991.3 Mb with a scaffold N50 value of 56.7 Mb. Evolutionary analysis of the genome suggests that C. alismatifolia diverged from Zingiber officinale about 9.7 million years ago, after it underwent a whole-genome duplication. Transcriptome analysis was performed on bracts at five developmental stages. Nine highly expressed genes were identified, encoding for six enzymes downstream of the anthocyanin biosynthetic pathway. Of these, one gene encoding F3'5'H might be a key node in the regulation of bract color formation. Co-expression network analysis showed that MYB, bHLH, NAC, and ERF transcription factors collectively regulated color formation in the bracts. Characterization of terpenoid biosynthesis genes revealed their dispersal and tandem duplications, both of which contributed greatly to the increase in the number of terpene synthase genes in C. alismatifolia, especially to species-specific expansion of sesquiterpene synthase genes. This work facilitates understanding of genetic basis of anthocyanin and terpenoid biosynthesis and could accelerate the selective breeding of C. alismatifolia varieties with higher ornamental and medicinal value.

RSNET: inferring gene regulatory networks by a redundancy silencing and network enhancement technique

Background

Current gene regulatory network (GRN) inference methods are notorious for a great number of indirect interactions hidden in the predictions. Filtering out the indirect interactions from direct ones remains an important challenge in the reconstruction of GRNs. To address this issue, we developed a redundancy silencing and network enhancement technique (RSNET) for inferring GRNs.

Results

To assess the performance of RSNET method, we implemented the experiments on several gold-standard networks by using simulation study, DREAM challenge dataset and Escherichia coli network. The results show that RSNET method performed better than the compared methods in sensitivity and accuracy. As a case of study, we used RSNET to construct functional GRN for apple fruit ripening from gene expression data.

Conclusions

In the proposed method, the redundant interactions including weak and indirect connections are silenced by recursive optimization adaptively, and the highly dependent nodes are constrained in the model to keep the real interactions. This study provides a useful tool for inferring clean networks.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04696-w.

Genomic adaptations for arboreal locomotion in Asian flying treefrogs

SignificanceTo adapt to arboreal lifestyles, treefrogs have evolved a suite of complex traits that support vertical movement and gliding, thus presenting a unique case for studying the genetic basis for traits causally linked to vertical niche expansion. Here, based on two de novo-assembled Asian treefrog genomes, we determined that genes involved in limb development and keratin cytoskeleton likely played a role in the evolution of their climbing systems. Behavioral and morphological evaluation and time-ordered gene coexpression network analysis revealed the developmental patterns and regulatory pathways of the webbed feet used for gliding in Rhacophorus kio.

Transcriptome characterization of Larrea tridentata and identification of genes associated with phenylpropanoid metabolic pathways

Larrea tridentata (Sesse and Moc. ex DC.) Coville (family: Zygophyllaceae) is an aromatic evergreen shrub with resin-covered leaves, known to use in traditional medicine for diverse ailments. It also has immense pharmacological significance due to presence of powerful phenylpropanoids antioxidant, nordihydroguaiaretic acid (NDGA). The RNA sequence/transcriptome analyses connect the genomic information into the discovery of gene function. Hence, the acquaint analysis of L. tridentata is in lieu to characterize the transcriptome, and to identify the candidate genes involved in the phenylpropanoid biosynthetic pathway. To gain molecular insight, the bioinformatics analysis of transcriptome was performed. The total bases covered 48,630 contigs of length greater than 200 bp and above came out to 21,590,549 with an average GC content of 45% and an abundance of mononucleotide, SSR, including C3H, FAR1, and MADS transcription gene families. The best enzyme commission (EC) classification obtained from the assembled sequences represented major abundant enzyme classes e.g., RING-type E3 ubiquitin transferase and non-specific serine/threonine protein kinase. The KEGG pathway analysis mapped into 377 KEGG different metabolic pathways. The enrichment of phenylpropanoid biosynthesis pathways (22 genes i.e., phenylalanine ammonia-lyase, trans-cinnamate 4-monooxygenase, 4-coumarate-CoA ligase, cinnamoyl-CoA reductase, beta-glucosidase, shikimate O-hydroxycinnamoyl transferase, 5-O-(4-coumaroyl)-D-quinate 3'-monooxygenase, cinnamyl-alcohol dehydrogenase, peroxidase, coniferyl-alcohol glucosyltransferase, caffeoyl shikimate esterase, caffeoyl-CoA O-methyltransferase, caffeate O-methyltransferase, coniferyl-aldehyde dehydrogenase, feruloyl-CoA 6-hydroxylase, and ferulate-5-hydroxylase), and expression profile indicated antioxidant, anti-arthritic, and anticancer properties of L. tridentata. The present results could provide an important resource for squeezing biotechnological applications of L. tridentata.

Universal gene co-expression network reveals receptor-like protein genes involved in broad-spectrum resistance in pepper (Capsicum annuum L.)

Receptor-like proteins (RLPs) on plant cells have been implicated in immune responses and developmental processes. Although hundreds of RLP genes have been identified in plants, only a few RLPs have been functionally characterized in a limited number of plant species. Here, we identified RLPs in the pepper (Capsicum annuum) genome and performed comparative transcriptomics coupled with the analysis of conserved gene co-expression networks (GCNs) to reveal the role of core RLP regulators in pepper-pathogen interactions. A total of 102 RNA-seq datasets of pepper plants infected with four pathogens were used to construct CaRLP-targeted GCNs (CaRLP-GCNs). Resistance-responsive CaRLP-GCNs were merged to construct a universal GCN. Fourteen hub CaRLPs, tightly connected with defense-related gene clusters, were identified in eight modules. Based on the CaRLP-GCNs, we evaluated whether hub CaRLPs in the universal GCN are involved in the biotic stress response. Of the nine hub CaRLPs tested by virus-induced gene silencing, three genes (CaRLP264, CaRLP277, and CaRLP351) showed defense suppression with less hypersensitive response-like cell death in race-specific and non-host resistance response to viruses and bacteria, respectively, and consistently enhanced susceptibility to Ralstonia solanacearum and/or Phytophthora capsici. These data suggest that key CaRLPs are involved in the defense response to multiple biotic stresses and can be used to engineer a plant with broad-spectrum resistance. Together, our data show that generating a universal GCN using comprehensive transcriptome datasets can provide important clues to uncover genes involved in various biological processes.

NK cell receptor and ligand composition influences the clearance of SARS-CoV-2

To explore how the immune system controls clearance of SARS-CoV-2, we used a single-cell, mass cytometry-based proteomics platform to profile the immune systems of 21 patients who had recovered from SARS-CoV-2 infection without need for admission to an intensive care unit or for mechanical ventilation. We focused on receptors involved in interactions between immune cells and virus-infected cells. We found that the diversity of receptor repertoires on natural killer (NK) cells was negatively correlated with the viral clearance rate. In addition, NK subsets expressing the receptor DNAM1 were increased in patients who more rapidly recovered from infection. Ex vivo functional studies revealed that NK subpopulations with high DNAM1 expression had cytolytic activities in response to target cell stimulation. We also found that SARS-CoV-2 infection induced the expression of CD155 and nectin-4, ligands of DNAM1 and its paired coinhibitory receptor TIGIT, which counterbalanced the cytolytic activities of NK cells. Collectively, our results link the cytolytic immune responses of NK cells to the clearance of SARS-CoV-2 and show that the DNAM1 pathway modulates host-pathogen interactions during SARS-CoV-2 infection.

Time-ordered dysregulated ceRNA networks reveal disease progression and diagnostic biomarkers in ischemic and dilated cardiomyopathy

Ischemic cardiomyopathy (ICM) and dilated cardiomyopathy (DCM) are the two main causes of heart failure (HF). Despite similar clinical characteristics and common “HF pathways”, ICM and DCM are expected to have different personalized treatment strategies. The underlying mechanisms of ICM and DCM have yet to be fully elucidated. The present study developed a novel computational method for identifying dysregulated long noncoding RNA (lncRNA)–microRNA (miRNA)–mRNA competing endogenous RNA (ceRNA) triplets. Time-ordered dysregulated ceRNA networks were subsequently constructed to reveal the possible disease progression of ICM and DCM based on the method. Biological functional analysis indicated that ICM and DCM had similar features during myocardial remodeling, whereas their characteristics differed during progression. Specifically, disturbance of myocardial energy metabolism may be the main characteristic during DCM progression, whereas early inflammation and response to oxygen are the characteristics that may be specific to ICM. In addition, several panels of diagnostic biomarkers for differentiating non-heart failure (NF) and ICM (NF-ICM), NF-DCM, and ICM-DCM were identified. Our study reveals biological differences during ICM and DCM progression and provides potential diagnostic biomarkers for ICM and DCM.

Analyses of open-access multi-omics data sets reveal genetic and expression characteristics of maize ZmCCT family genes

Flowering in maize (Zea mays) is influenced by photoperiod. The CO, CO-like/COL and TOC1 (CCT) domain protein-encoding genes in maize, ZmCCTs, are particularly important for photoperiod sensitivity. However, little is known about CCT protein-encoding gene number across plant species or among maize inbred lines. Therefore, we analysed CCT protein-encoding gene number across plant species, and characterized ZmCCTs in different inbred lines, including structural variations (SVs), copy number variations (CNVs), expression under stresses, dark-dark (DD) and dark-light (DL) cycles, interaction network and associations with maize quantitative trait loci (QTLs) by referring to the latest v4 genome data of B73. Gene number varied greatly across plant species, more in polyploids than in diploids. The numbers of ZmCCTs identified were 58 in B73, 59 in W22, 48 in Mo17, and 57 in Huangzao4 for temperate maize inbred lines, and 68 in tropical maize inbred line SK. Some ZmCCTs underwent duplications and presented chromosome collinearity. Structural variations and CNVs were found but they had no germplasm specificity. Forty-two ZmCCTs responded to stresses. Expression of 37 ZmCCTs in embryonic leaves during seed germination of maize under DD and DL cycles was roughly divided into five patterns of uphill pattern, downhill-pattern, zigzag-pattern, └-pattern and ⅃-pattern, indicating some of them have a potential to perceive dark and/or dark-light transition. Thirty-three ZmCCTs were co-expressed with 218 other maize genes; and 24 ZmCCTs were associated with known QTLs. The data presented in this study will help inform further functions of ZmCCTs.

A transcriptomic view of cadmium retention in roots of cadmium-safe rice line (Oryza sativa L.)

A better understanding of the mechanisms controlling cadmium (Cd) accumulation in rice will benefit the development of strategies to minimize Cd accumulation in grains. A Cd-safe rice line designated D62B accumulated less than 0.2 mg Cd kg^-1 in brown rice due to its strong capacity for Cd retention in roots. Here transcriptomic was used to clarify the underlying mechanisms of Cd response in roots of D62B compared with a high Cd-accumulating line (Wujin4B). There were 777, 1058 differentially expressed genes (DEGs) in D62B and Wujin4B, respectively, when exposed to Cd. The functions of DEGs were clearly line-specific. Cell wall biosynthesis responded more intensively to Cd stress in D62B, facilitating Cd restriction. Meanwhile, more glutathione (GSH) and phytochelatins synthesized in D62B with the upregulation of sulphur and GSH metabolism. Besides, membrane proteins played critical roles in Cd response in D62B, whereas 18 terms involved in regulation were enriched in Wujin4B. Exogenous GSH further induced the expression of genes related to GSH metabolism and cell wall biosynthesis, leading to the retention of more Cd. Great responsiveness of cell wall biosynthesis and GSH metabolism could be considered the most important specific mechanisms for Cd retention in the roots of Cd-safe rice line.

The Molecular Mechanism of Multiple Organ Dysfunction and Targeted Intervention of COVID-19 Based on Time-Order Transcriptomic Analysis

Coronavirus disease 2019 (COVID-19) pandemic is caused by the novel coronavirus that has spread rapidly around the world, leading to high mortality because of multiple organ dysfunction; however, its underlying molecular mechanism is unknown. To determine the molecular mechanism of multiple organ dysfunction, a bioinformatics analysis method based on a time-order gene co-expression network (TO-GCN) was performed. First, gene expression profiles were downloaded from the gene expression omnibus database (GSE161200), and a TO-GCN was constructed using the breadth-first search (BFS) algorithm to infer the pattern of changes in the different organs over time. Second, Gene Ontology enrichment analysis was used to analyze the main biological processes related to COVID-19. The initial gene modules for the immune response of different organs were defined as the research object. The STRING database was used to construct a protein-protein interaction network of immune genes in different organs. The PageRank algorithm was used to identify five hub genes in each organ. Finally, the Comparative Toxicogenomics Database played an important role in exploring the potential compounds that target the hub genes. The results showed that there were two types of biological processes: the body's stress response and cell-mediated immune response involving the lung, trachea, and olfactory bulb (olf) after being infected by COVID-19. However, a unique biological process related to the stress response is the regulation of neuronal signals in the brain. The stress response was heterogeneous among different organs. In the lung, the regulation of DNA morphology, angiogenesis, and mitochondrial-related energy metabolism are specific biological processes related to the stress response. In particular, an effect on tracheal stress response was made by the regulation of protein metabolism and rRNA metabolism-related biological processes, as biological processes. In the olf, the distinctive stress responses consist of neural signal transmission and brain behavior. In addition, myeloid leukocyte activation and myeloid leukocyte-mediated immunity in response to COVID-19 can lead to a cytokine storm. Immune genes such as SRC, RHOA, CD40LG, CSF1, TNFRSF1A, FCER1G, ICAM1, LAT, LCN2, PLAU, CXCL10, ICAM1, CD40, IRF7, and B2M were predicted to be the hub genes in the cytokine storm. Furthermore, we inferred that resveratrol, acetaminophen, dexamethasone, estradiol, statins, curcumin, and other compounds are potential target drugs in the treatment of COVID-19.

The Performance Comparison of Gene Co-expression Networks of Breast and Prostate Cancer using Different Selection Criteria

Gene co-expression networks (GCN) present undirected relations between genes to understand molecular structures behind the diseases, including cancer. The utilization of various biological datasets and gene network inference (GNI) algorithms can reveal meaningful gene-gene interactions of GCNs. This study applies three GNI algorithms on mRNA gene expression, RNA-Seq, and miRNA-target genes datasets to infer GCNs of breast and prostate cancers. To evaluate the performance of the GCNs, we utilize overlap analysis via literature data, topological assessment, and Gene Ontology-based biological assessment. The results emphasize how the selection of biological datasets and GNI algorithms affect the performance results on different evaluation criteria. GCNs on microarray gene expression data slightly outperform in overlap analysis. Also, GCNs on RNA-Seq and gene expression datasets follow scale-free topology. The biological assessment results are close to each other on all biological datasets. C3NET algorithm-based GCNs did not contain any biological assessment modules; therefore, it is not optimal for biological assessment. GNI algorithms' selection did not change the overlap analysis and topological assessment results. Our primary objective is to compare the performance results of biological datasets and GNI algorithms based on different evaluation criteria. For this purpose, we developed the GNIAP R package that enables users to select different GNI algorithms to infer GCNs. The GNIAP R package also provides literature-based overlap analysis, and topological and biological analyses on GCNs. Users can access the GNIAP R package via https://github.com/ozgurcingiz/GNIAP .