SIRT1 and ZNF350 as novel biomarkers for osteoporosis: a bioinformatics analysis and experimental validation

BACKGROUND

Osteoporosis (OP) is characterized by bone mass decrease and bone tissue microarchitectural deterioration in bone tissue. This study identified potential biomarkers for early diagnosis of OP and elucidated the mechanism of OP.

METHODS

Gene expression profiles were downloaded from Gene Expression Omnibus (GEO) for the GSE56814 dataset. A gene co-expression network was constructed using weighted gene co-expression network analysis (WGCNA) to identify key modules associated with healthy and OP samples. Functional enrichment analysis was conducted using the R clusterProfiler package for modules to construct the transcriptional regulatory factor networks. We used the "ggpubr" package in R to screen for differentially expressed genes between the two samples. Gene set variation analysis (GSVA) was employed to further validate hub gene expression levels between normal and OP samples using RT-PCR and immunofluorescence to evaluate the potential biological changes in various samples.

RESULTS

There was a distinction between the normal and OP conditions based on the preserved significant module. A total of 100 genes with the highest MM scores were considered key genes. Functional enrichment analysis suggested that the top 10 biological processes, cellular component and molecular functions were enriched. The Toll-like receptor signaling pathway, TNF signaling pathway, PI3K-Akt signaling pathway, osteoclast differentiation, JAK-STAT signaling pathway, and chemokine signaling pathway were identified by Kyoto Encyclopedia of Genes and Genomes pathway analysis. SIRT1 and ZNF350 were identified by Wilcoxon algorithm as hub differentially expressed transcriptional regulatory factors that promote OP progression by affecting oxidative phosphorylation, apoptosis, PI3K-Akt-mTOR signaling, and p53 pathway. According to RT-PCR and immunostaining results, SIRT1 and ZNF350 levels were significantly higher in OP samples than in normal samples.

CONCLUSION

SIRT1 and ZNF350 are important transcriptional regulatory factors for the pathogenesis of OP and may be novel biomarkers for OP treatment.

NnSUS1 encodes a sucrose synthase involved in sugar accumulation in lotus seed cotyledons

Fresh lotus seeds are gaining favor with consumers for their crunchy texture and natural sweetness. However, the intricacies of sugar accumulation in lotus seeds remain elusive, which greatly hinders the quality improvement of fresh lotus seeds. This study endeavors to elucidate this mechanism by identifying and characterizing the sucrose synthase (SUS) gene family in lotus. Comprising five distinct members, namely NnSUS1 to NnSUS5, each gene within this family features a C-terminal glycosyl transferase1 (GT1) domain. Among them, NnSUS1 is the predominately expressed gene, showing high transcript abundance in the floral organs and cotyledons. NnSUS1 was continuously up-regulated from 6 to 18 days after pollination (DAP) in lotus cotyledons. Furthermore, NnSUS1 demonstrates co-expression relationships with numerous genes involved in starch and sucrose metabolism. To investigate the function of NnSUS1, a transient overexpression system was established in lotus cotyledons, which confirmed the gene's contribution to sugar accumulation. Specifically, transient overexpression of NnSUS1 in seed cotyledons leads to a significant increase in the levels of total soluble sugar, including sucrose and fructose. These findings provide valuable theoretical insights for improving sugar content in lotus seeds through molecular breeding methods.

Functional identification of specialized diterpene synthases from Chamaecyparis obtusa and C. obtusa var. formosana to illustrate the putative evolution of diterpene synthases in Cupressaceae

Chamaecyparis obtusa and C. obtusa var. formosana of the Cupressaceae family are well known for their fragrance and excellent physical properties. To investigate the biosynthesis of unique diterpenoid compounds, diterpene synthase genes for specialized metabolite synthesis were cloned from C. obtusa and C. obtusa var. formosana. Using an Escherichia coli co-expression system, eight diterpene synthases (diTPSs) were characterized. CoCPS and CovfCPS are class II monofunctional (+)-copalyl diphosphate synthases [(+)-CPSs]. Class I monofunctional CoLS and CovfLS convert (+)-copalyl diphosphate [(+)-CPP] to levopimaradiene, CoBRS, CovfBRS1, and CovfBRS3 convert (+)-CPP to (-)-beyerene, and CovfSDS converts (+)-CPP to (-)-sandaracopimaradiene. These enzymes are all monofunctional diterpene syntheses in Cupressaceae family of gymnosperm, and differ from those in Pinaceae. The discovery of the enzyme responsible for the biosynthesis of tetracyclic diterpene (-)-beyerene was characterized for the first time. Diterpene synthases with different catalytic functions exist in closely related species within the Cupressaceae family, indicating that this group of monofunctional diterpene synthases is particularly prone to the evolution of new functions and development of species-specific specialized diterpenoid constituents.

Exploiting the biocatalytic potential of co-expressed l-fucose isomerase and d-tagatose 3-epimerase for the biosynthesis of 6-deoxy-l-sorbose

6-Deoxy-l-sorbose (6-DLS) is an imperative rare sugar employed in food, agriculture, pharmaceutical and cosmetic industeries. However, it is a synthetic and very expensive rare sugars, previously synthesized by chemo-enzymatic methods through a long chain of chemical processes. Recently, enzymatic synthesis of rare sugars has attracted a lot of attention due to its advantages over synthetic methods. In this work, a promising approach for the synthesis of 6-DLS from an inexpensive sugar l-fucose was identified. The genes for l-fucose isomerase from Paenibacillus rhizosphaerae (Pr-LFI) and genes for d-tagatose-3-epimerase from Caballeronia fortuita (Cf-DTE) have been used for cloning and co-expression in Escherichia coli, developed a recombinant plasmid harboring pANY1-Pr-LFI/Cf-DTE vector. The recombinant co-expression system exhibited an optimum activity at 50 °C of temperature and pH 6.5 in the presence of Co2+ metal ion which inflated the catalytic activity by 6.8 folds as compared to control group with no metal ions. The recombinant co-expressed system was stable up to more than 50 % relative activity after 12 h and revealed a melting temperature (Tm) of 63.38 °C exhibiting half-life of 13.17 h at 50 °C. The co-expression system exhibited, 4.93, 11.41 and 16.21 g/L of 6-DLS production from initial l-fucose concentration of 30, 70 and 100 g/L, which equates to conversion yield of 16.44 %, 16.30 % and 16.21 % respectively. Generally, this study offers a promising strategy for the biological production of 6-DLS from an inexpensive substrate l-fucose in slightly acidic conditions with the aid of co-expression system harboring Pr-LFI and CF-DTE genes.

Enriched atlas of lncRNA and protein-coding genes for the GRCg7b chicken assembly and its functional annotation across 47 tissues

Gene atlases for livestock are steadily improving thanks to new genome assemblies and new expression data improving the gene annotation. However, gene content varies across databases due to differences in RNA sequencing data and bioinformatics pipelines, especially for long non-coding RNAs (lncRNAs) which have higher tissue and developmental specificity and are harder to consistently identify compared to protein coding genes (PCGs). As done previously in 2020 for chicken assemblies galgal5 and GRCg6a, we provide a new gene atlas, lncRNA-enriched, for the latest GRCg7b chicken assembly, integrating "NCBI RefSeq", "EMBL-EBI Ensembl/GENCODE" reference annotations and other resources such as FAANG and NONCODE. As a result, the number of PCGs increases from 18,022 (RefSeq) and 17,007 (Ensembl) to 24,102, and that of lncRNAs from 5789 (RefSeq) and 11,944 (Ensembl) to 44,428. Using 1400 public RNA-seq transcriptome representing 47 tissues, we provided expression evidence for 35,257 (79%) lncRNAs and 22,468 (93%) PCGs, supporting the relevance of this atlas. Further characterization including tissue-specificity, sex-differential expression and gene configurations are provided. We also identified conserved miRNA-hosting genes with human counterparts, suggesting common function. The annotated atlas is available at gega.sigenae.org.

Inference of gene networks using gene expression data with applications

Gene networks (GNs) use graphs to represent the interaction relationships between genes. Large-scale GNs are often sparse and contain hub genes that interact with many other genes. In this paper, we propose a novel method called NetARD, which utilizes Automatic Relevance Determination (ARD) to estimate partial correlations, to infer GNs with the hub genes from gene expression data. We test NetARD on simulated GNs and in silico GNs, and it outperforms existing methods. In our high-throughput gene expression data analysis, we integrate the NetARD into a method called GN Co-expression Extension (GNCE). This approach infers the GNs of co-expressed genes, with genes from a predefined GN serving as hub genes. We validate this approach by extending the core GN of transcription factor genes of E. coli using microarray data. In an application example, we identify biological process (BP) Gene Ontology (GO) terms that are significantly involved in cancer progression. This task is accomplished by analyzing the GN inferred through GNCE using the core GN associated with the colorectal cancer pathway and RNA-seq data.

In Silico Comparison of WRKY Transcription Factors in Wild and Cultivated Soybean and Their Co-expression Network Arbitrating Disease Resistance

WRKY Transcription factors (TFs) play critical roles in plant defence mechanisms that are activated in response to biotic and abiotic stresses. However, information on the Glycine soja WRKYs (GsoWRKYs) is scarce. Owing to its importance in soybean breeding, here we identified putative WRKY TFs in wild soybean, and compared the results with Glycine max WRKYs (GmaWRKYs) by phylogenetic, conserved motif, and duplication analyses. Moreover, we explored the expression trends of WRKYs in G. max (oomycete, fungi, virus, bacteria, and soybean cyst nematode) and G. soja (soybean cyst nematode), and identified commonly expressed WRKYs and their co-expressed genes. We identified, 181 and 180 putative WRKYs in G. max and G. soja, respectively. Though the number of WRKYs in both studied species is almost the same, they differ in many ways, i.e., the number of WRKYs on corresponding chromosomes, conserved domain structures, WRKYGQK motif variants, and zinc-finger motifs. WRKYs in both species grouped in three major clads, i.e., I-III, where group-II had sub-clads IIa-IIe. We found that GsoWRKYs expanded mostly through segmental duplication. A large number of WRKYs were expressed in response to biotic stresses, i.e., Phakospora pachyrhizi, Phytoplasma, Heterodera glycines, Macrophomina phaseolina, and Soybean mosaic virus; 56 GmaWRKYs were commonly expressed in soybean plants infected with these diseases. Finally, 30 and 63 GmaWRKYs and GsoWRKYs co-expressed with 205 and 123 non-WRKY genes, respectively, indicating that WRKYs play essential roles in biotic stress tolerance in Glycine species.

Prioritizing Metabolic Gene Regulators through Multi-Omic Network Integration in Maize

Elucidating gene regulatory networks (GRNs) is a major area of study within plant systems biology. Phenotypic traits are intricately linked to specific gene expression profiles. These expression patterns arise primarily from regulatory connections between sets of transcription factors (TFs) and their target genes. In this study, we integrated publicly available co-expression networks derived from more than 6,000 RNA-seq samples, 283 protein-DNA interaction assays, and 16 million of SNPs used to identify expression quantitative loci (eQTL), to construct TF-target networks. In total, we analyzed ~4.6M interactions to generate four distinct types of TF-target networks: co-expression, protein-DNA interaction (PDI), trans-expression quantitative loci (trans-eQTL), and cis-eQTL combined with PDIs. To improve the functional annotation of TFs based on its target genes, we implemented three different strategies to integrate these four types of networks. We subsequently evaluated the effectiveness of our method through loss-of function mutant and random networks. The multi-network integration allowed us to identify transcriptional regulators of hormone-, metabolic- and development-related processes. Finally, using the topological properties of the fully integrated network, we identified potentially functional redundant TF paralogs. Our findings retrieved functions previously documented for numerous TFs and revealed novel functions that are crucial for informing the design of future experiments. The approach here-described lays the foundation for the integration of multi-omic datasets in maize and other plant systems.

Identification of environment-insensitive genes for oil content by combination of transcriptome and genome-wide association analysis in rapeseed

BACKGROUND

The primary objective of rapeseed breeding is to enhance oil content, which is predominantly influenced by environmental factors. However, the molecular mechanisms underlying the impact of these environmental factors on oil accumulation remain inadequately elucidated. In this study, we used transcriptome data from two higher (HOC) and two lower oil content (LOC) inbred lines at 35 days after pollination (DAP) to investigate genes exhibiting stable expression across three different environments. Meanwhile, a genome-wide association study (GWAS) was utilized to detect candidate genes exhibiting significant associations with seed oil content across three distinct environments.

RESULTS

The study found a total of 405 stable differentially expressed genes (DEGs), including 25 involved in lipid/fatty acid metabolism and 14 classified as transcription factors. Among these genes, BnBZIP10-A09, BnMYB61-A06, BnAPA1-A08, BnPAS2-A10, BnLCAT3-C05 and BnKASIII-C09 were also found to exhibit significant associations with oil content across multiple different environments based on GWAS of 50 re-sequenced semi-winter rapeseed inbred lines and previously reported intervals. Otherwise, we revealed the presence of additive effects among BnBZIP10-A09, BnKASIII-C09, BnPAS2-A10 and BnAPA1-A08, resulting in a significant increase in seed oil content. Meanwhile, the majority of these stable DEGs are interconnected either directly or indirectly through co-expression network analysis, thereby giving rise to an elaborate molecular network implicated in the potential regulation of seed oil accumulation and stability.

CONCLUSIONS

The combination of transcription and GWAS revealed that natural variation in six environment-insensitive gene regions exhibited significant correlations with seed oil content phenotypes. These results provide important molecular marker information for us to further improve oil content accumulation and stability in rapeseed.

Genome-wide analysis of long noncoding RNAs in response to salt stress in Nicotiana tabacum

BACKGROUND

Long noncoding RNAs (lncRNAs) have been shown to play important roles in the response of plants to various abiotic stresses, including drought, heat and salt stress. However, the identification and characterization of genome-wide salt-responsive lncRNAs in tobacco (Nicotiana tabacum L.) have been limited. Therefore, this study aimed to identify tobacco lncRNAs in roots and leaves in response to different durations of salt stress treatment.

RESULTS

A total of 5,831 lncRNAs were discovered, with 2,428 classified as differentially expressed lncRNAs (DElncRNAs) in response to salt stress. Among these, only 214 DElncRNAs were shared between the 2,147 DElncRNAs in roots and the 495 DElncRNAs in leaves. KEGG pathway enrichment analysis revealed that these DElncRNAs were primarily associated with pathways involved in starch and sucrose metabolism in roots and cysteine and methionine metabolism pathway in leaves. Furthermore, weighted gene co-expression network analysis (WGCNA) identified 15 co-expression modules, with four modules strongly linked to salt stress across different treatment durations (MEsalmon, MElightgreen, MEgreenyellow and MEdarkred). Additionally, an lncRNA-miRNA-mRNA network was constructed, incorporating several known salt-associated miRNAs such as miR156, miR169 and miR396.

CONCLUSIONS

This study enhances our understanding of the role of lncRNAs in the response of tobacco to salt stress. It provides valuable information on co-expression networks of lncRNA and mRNAs, as well as networks of lncRNAs-miRNAs-mRNAs. These findings identify important candidate lncRNAs that warrant further investigation in the study of plant-environment interactions.

Modular bioengineering of whole-cell catalysis for sialo-oligosaccharide production: coordinated co-expression of CMP-sialic acid synthetase and sialyltransferase

BACKGROUND

In whole-cell bio-catalysis, the biosystems engineering paradigm shifts from the global reconfiguration of cellular metabolism as in fermentation to a more focused, and more easily modularized, optimization of comparably short cascade reactions. Human milk oligosaccharides (HMO) constitute an important field for the synthetic application of cascade bio-catalysis in resting or non-living cells. Here, we analyzed the central catalytic module for synthesis of HMO-type sialo-oligosaccharides, comprised of CMP-sialic acid synthetase (CSS) and sialyltransferase (SiaT), with the specific aim of coordinated enzyme co-expression in E. coli for reaction flux optimization in whole cell conversions producing 3'-sialyllactose (3SL).

RESULTS

Difference in enzyme specific activity (CSS from Neisseria meningitidis: 36 U/mg; α2,3-SiaT from Pasteurella dagmatis: 5.7 U/mg) was compensated by differential protein co-expression from tailored plasmid constructs, giving balance between the individual activities at a high level of both (α2,3-SiaT: 9.4 × 102 U/g cell dry mass; CSS: 3.4 × 102 U/g cell dry mass). Finally, plasmid selection was guided by kinetic modeling of the coupled CSS-SiaT reactions in combination with comprehensive analytical tracking of the multistep conversion (lactose, N-acetyl neuraminic acid (Neu5Ac), cytidine 5'-triphosphate; each up to 100 mM). The half-life of SiaT in permeabilized cells (≤ 4 h) determined the efficiency of 3SL production at 37 °C. Reaction at 25 °C gave 3SL (40 ± 4 g/L) in ∼ 70% yield within 3 h, reaching a cell dry mass-specific productivity of ∼ 3 g/(g h) and avoiding intermediary CMP-Neu5Ac accumulation.

CONCLUSIONS

Collectively, balanced co-expression of CSS and SiaT yields an efficient (high-flux) sialylation module to support flexible development of E. coli whole-cell catalysts for sialo-oligosaccharide production.

Bayesian Optimized sample-specific Networks Obtained By Omics data (BONOBO)

Gene regulatory networks (GRNs) are effective tools for inferring complex interactions between molecules that regulate biological processes and hence can provide insights into drivers of biological systems. Inferring co-expression networks is a critical element of GRN inference as the correlation between expression patterns may indicate that genes are coregulated by common factors. However, methods that estimate co-expression networks generally derive an aggregate network representing the mean regulatory properties of the population and so fail to fully capture population heterogeneity. To address these concerns, we introduce BONOBO (Bayesian Optimized Networks Obtained By assimilating Omics data), a scalable Bayesian model for deriving individual sample-specific co-expression networks by recognizing variations in molecular interactions across individuals. For every sample, BONOBO assumes a Gaussian distribution on the log-transformed centered gene expression and a conjugate prior distribution on the sample-specific co-expression matrix constructed from all other samples in the data. Combining the sample-specific gene expression with the prior distribution, BONOBO yields a closed-form solution for the posterior distribution of the sample-specific co-expression matrices, thus making the method extremely scalable. We demonstrate the utility of BONOBO in several contexts, including analyzing gene regulation in yeast transcription factor knockout studies, prognostic significance of miRNA-mRNA interaction in human breast cancer subtypes, and sex differences in gene regulation within human thyroid tissue. We find that BONOBO outperforms other sample-specific co-expression network inference methods and provides insight into individual differences in the drivers of biological processes.

Long noncoding RNAs emerge from transposon-derived antisense sequences and may contribute to infection stage-specific transposon regulation in a fungal phytopathogen

BACKGROUND

The genome of the obligate biotrophic phytopathogenic barley powdery mildew fungus Blumeria hordei is inflated due to highly abundant and possibly active transposable elements (TEs). In the absence of the otherwise common repeat-induced point mutation transposon defense mechanism, noncoding RNAs could be key for regulating the activity of TEs and coding genes during the pathogenic life cycle.

RESULTS

We performed time-course whole-transcriptome shotgun sequencing (RNA-seq) of total RNA derived from infected barley leaf epidermis at various stages of fungal pathogenesis and observed significant transcript accumulation and time point-dependent regulation of TEs in B. hordei. Using a manually curated consensus database of 344 TEs, we discovered phased small RNAs mapping to 104 consensus transposons, suggesting that RNA interference contributes significantly to their regulation. Further, we identified 5,127 long noncoding RNAs (lncRNAs) genome-wide in B. hordei, of which 823 originated from the antisense strand of a TE. Co-expression network analysis of lncRNAs, TEs, and coding genes throughout the asexual life cycle of B. hordei points at extensive positive and negative co-regulation of lncRNAs, subsets of TEs and coding genes.

CONCLUSIONS

Our work suggests that similar to mammals and plants, fungal lncRNAs support the dynamic modulation of transcript levels, including TEs, during pivotal stages of host infection. The lncRNAs may support transcriptional diversity and plasticity amid loss of coding genes in powdery mildew fungi and may give rise to novel regulatory elements and virulence peptides, thus representing key drivers of rapid evolutionary adaptation to promote pathogenicity and overcome host defense.

Identification of hub genes associated with somatic cell score in dairy cow

CONTEXT

Somatic cell count (SCC) is used as an indicator of udder health. The log transformation of SCC is called somatic cell score (SCS).

AIM

Several QTL and genes have been identified that are associated with SCS. This study aimed to identify the most important genes associated with SCS.

METHODS

This study compiled 168 genes that were reported to be significantly linked to SCS. Pathway analysis and network analysis were used to identify hub genes.

KEY RESULTS

Pathway analysis of these genes identified 73 gene ontology (GO) terms associated with SCS. These GO terms are associated with molecular function, biological processes, and cellular components, and the identified pathways are directly or indirectly linked with the immune system. In this study, a gene network was constructed, and from this network, the 17 hub genes (CD4, CXCL8, TLR4, STAT1, TLR2, CXCL9, CCR2, IGF1, LEP, SPP1, GH1, GHR, VWF, TNFSF11, IL10RA, NOD2, and PDGFRB) associated to SCS were identified. The subnetwork analysis yielded 10 clusters, with cluster 1 containing all identified hub genes (except for the VWF gene).

CONCLUSION

Most hub genes and pathways identified in our study were mainly involved in inflammatory and cytokine responses.

IMPLICATIONS

Result obtained in current study provides knowledge of the genetic basis and biological mechanisms controlling SCS. Therefore, the identified hub genes may be regarded as the main gene for the genomic selection of mastitis resistance.

Interaction between phenylpropane metabolism and oil accumulation in the developing seed of Brassica napus revealed by high temporal-resolution transcriptomes

BACKGROUND

Brassica napus is an important oilseed crop providing high-quality vegetable oils for human consumption and non-food applications. However, the regulation between embryo and seed coat for the synthesis of oil and phenylpropanoid compounds remains largely unclear.

RESULTS

Here, we analyzed the transcriptomes in developing seeds at 2-day intervals from 14 days after flowering (DAF) to 64 DAF. The 26 high-resolution time-course transcriptomes are clearly clustered into five distinct groups from stage I to stage V. A total of 2217 genes including 136 transcription factors, are specifically expressed in the seed and show high temporal specificity by being expressed only at certain stages of seed development. Furthermore, we analyzed the co-expression networks during seed development, which mainly included master regulatory transcription factors, lipid, and phenylpropane metabolism genes. The results show that the phenylpropane pathway is prominent during seed development, and the key enzymes in the phenylpropane metabolic pathway, including TT5, BAN, and the transporter TT19, were directly or indirectly related to many key enzymes and transcription factors involved in oil accumulation. We identified candidate genes that may regulate seed oil content based on the co-expression network analysis combined with correlation analysis of the gene expression with seed oil content and seed coat content.

CONCLUSIONS

Overall, these results reveal the transcriptional regulation between lipid and phenylpropane accumulation during B. napus seed development. The established co-expression networks and predicted key factors provide important resources for future studies to reveal the genetic control of oil accumulation in B. napus seeds.

JustRNA: a database of plant long noncoding RNA expression profiles and functional network

Long noncoding RNAs (lncRNAs) are regulatory RNAs involved in numerous biological processes. Many plant lncRNAs have been identified, but their regulatory mechanisms remain largely unknown. A resource that enables the investigation of lncRNA activity under various conditions is required because the co-expression between lncRNAs and protein-coding genes may reveal the effects of lncRNAs. This study developed JustRNA, an expression profiling resource for plant lncRNAs. The platform currently contains 1 088 565 lncRNA annotations for 80 plant species. In addition, it includes 3692 RNA-seq samples derived from 825 conditions in six model plants. Functional network reconstruction provides insight into the regulatory roles of lncRNAs. Genomic association analysis and microRNA target prediction can be employed to depict potential interactions with nearby genes and microRNAs, respectively. Subsequent co-expression analysis can be employed to strengthen confidence in the interactions among genes. Chromatin immunoprecipitation sequencing data of transcription factors and histone modifications were integrated into the JustRNA platform to identify the transcriptional regulation of lncRNAs in several plant species. The JustRNA platform provides researchers with valuable insight into the regulatory mechanisms of plant lncRNAs. JustRNA is a free platform that can be accessed at http://JustRNA.itps.ncku.edu.tw.

A DNA vaccine (EG95-PT1/2/3-IL2) encoding multi-epitope antigen and IL-2 provokes efficient and long-term immunity to echinococcosis

Echinococcosis is a highly prevalent global zoonosis, and vaccines are required. The commercial vaccine based on a protein-based subunit (EG95), however, is limited by its insufficient cellular immunity, a short protection period, and limited prevention against novel mutant strains. Herein, we applied bioinformatics to develop a DNA vaccine (pEG95-IL2) expressing both multi-epitope-based antigens (EG95-PT1/2/3) and an IL-2 adjuvant to regulate T cell differentiation and memory cell response. EG95-PT1/2/3 was screened with hierarchical structure prediction from the epitope conformation of B cells with high confidence across various species to guarantee immunogenicity. Importantly, cationic arginine-rich lipid nanoparticles (RNP) were utilized as a delivery vehicle to form lipoplexes that had a transfection efficiency of nearly two orders of magnitude greater than that of commercial reagents (Lipofectamine 2000 and polyethyleneimine) with both immune and nonimmune cells (DC2.4 and L929 cells, respectively). RNP/pEG95-IL2 lipoplexes displayed a robust and long-term antigen expression, as well as adjuvant effects during the immunization. Consequently, intramuscular injection of RNP/pEG95-IL2 elicited similar humoral immune responses and significantly greater cellular responses in mice when compared with those of the commercial vaccine. In addition, the inoculation protocol of RNP/pEG95-IL2 with sequential booster further strengthens cellular immunity in comparison with the homologous booster. Those findings provide a promising strategy for improving plasmid vaccine efficacy.

CryptoCEN: A Co-Expression Network for Cryptococcus neoformans reveals novel proteins involved in DNA damage repair

Elucidating gene function is a major goal in biology, especially among non-model organisms. However, doing so is complicated by the fact that molecular conservation does not always mirror functional conservation, and that complex relationships among genes are responsible for encoding pathways and higher-order biological processes. Co-expression, a promising approach for predicting gene function, relies on the general principal that genes with similar expression patterns across multiple conditions will likely be involved in the same biological process. For Cryptococcus neoformans, a prevalent human fungal pathogen greatly diverged from model yeasts, approximately 60% of the predicted genes in the genome lack functional annotations. Here, we leveraged a large amount of publicly available transcriptomic data to generate a C. neoformans Co-Expression Network (CryptoCEN), successfully recapitulating known protein networks, predicting gene function, and enabling insights into the principles influencing co-expression. With 100% predictive accuracy, we used CryptoCEN to identify 13 new DNA damage response genes, underscoring the utility of guilt-by-association for determining gene function. Overall, co-expression is a powerful tool for uncovering gene function, and decreases the experimental tests needed to identify functions for currently under-annotated genes.

Efficient expression of a cytokine combination in Saccharomyces cerevisiae for cultured meat production

Cultured meat technology is a novel and promising alternative strategy for meat production, and it provides an efficient, safe, and sustainable way to supply animal protein. Cytokines play an important role in promoting the rapid proliferation of cells, but the high cost and potential food safety concerns of commercial cytokines have hindered their application in large-scale cultured meat production. Herein, Saccharomyces cerevisiae C800 was used as a starting strain in which four cytokines were exogenously expressed simultaneously using the Cre-loxP system, including long-chain human insulin growth factor-1, platelet-derived growth factor-BB, basic fibroblast growth factor, and epidermal growth factor. Through promoter optimization, endogenous protease knockout, genomic co-expression, expression frame gene order optimization, and fermentation optimization, a recombinant strain CPK2B2 co-expressing four cytokines was obtained with a yield of 18.35 mg/L. After cell lysis and filter sterilization, the CPK2B2 lysate was directly added to the culture medium of porcine muscle satellite cells (MuSCs). CPK2B2 lysate promoted the growth of MuSCs and increased the proportion of G2/S cells and EdU⁺ cells significantly, indicating its efficacy in promoting cell proliferation. This study provides a simple and cost-saving strategy by using S. cerevisiae to produce a recombinant cytokine combination for cultured meat production.

The landscape of the long non-coding RNAs in developing mouse retinas

BACKGROUND

The long non-coding RNAs (lncRNAs) are critical regulators of diverse biological processes. Nevertheless, a global view of its expression and function in the mouse retina, a crucial model for neurogenesis study, still needs to be made available.

RESULTS

Herein, by integrating the established gene models and the result from ab initio prediction using short- and long-read sequencing, we characterized 4,523 lncRNA genes (MRLGs) in developing mouse retinas (from the embryonic day of 12.5 to the neonatal day of P28), which was so far the most comprehensive collection of retinal lncRNAs. Next, derived from transcriptomics analyses of different tissues and developing retinas, we found that the MRLGs were highly spatiotemporal specific in expression and played essential roles in regulating the genesis and function of mouse retinas. In addition, we investigated the expression of MRLGs in some mouse mutants and revealed that 97 intergenic MRLGs might be involved in regulating differentiation and development of retinal neurons through Math5, Isl1, Brn3b, NRL, Onecut1, or Onecut2 mediated pathways.

CONCLUSIONS

In summary, this work significantly enhanced our knowledge of lncRNA genes in mouse retina development and provided valuable clues for future exploration of their biological roles.

Genome-wide identification of the NAC transcription factor family and regulation of metabolites under salt stress in Isatis indigotica

NAC (NAM, ATAF1/2 and CUC2) transcription factors (TFs) are a class of TFs families unique to plants, which not only play an important role in the growth and developmental stages of plants but also function in response to stress and regulation of secondary metabolite biosynthesis. However, there are few studies on NAC genes in the medicinal plant Isatis indigotica. In this study, 96 IiNAC genes were identified based on the whole-genome data of I. indigotica, distributed in seven chromosomes and three contigs. IiNAC genes were structurally conserved and divided into 15 subgroups. Cis-elements were identified in the promoter region of the IiNAC gene in response to plant growth and development, abiotic stresses and hormones. In addition, transcriptome and metabolome data of I. indigotica leaves under salt stress were analyzed to construct a network of IiNAC gene co-expression and metabolite association. Ten differentially expressed IiNAC genes were co-expressed with 109 TFs, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that most of these genes were associated with plant growth and development and abiotic stress responses. Eleven IiNAC genes were positively associated with 72 metabolites. Eleven IiNAC genes were positively or negatively associated with 47 metabolites through 37 TFs. Commonly associated secondary metabolites include two terpenoids, abscisic acid and bilobalide, two flavonoids, dihydrokaempferol and syringaldehyde, a coumarin, 7-methoxycoumarin, an alkaloid, lupinine, and quinone dihydrotanshinone I. This study provides important data to support the identification of the NAC gene family in I. indigotica and the regulatory functions of IiNAC genes in metabolites under salt stress.

Diagnostic and therapeutic potential of LINC01929 as an oncogenic LncRNA in human cancers

Recent findings have shown the significant role of long non-coding RNAs in the pathogenesis of various cancers. In this regard, the variation in the expression of LINC01929 was explored in various cancers to explore its impact on the development of diverse malignancies and cancers. The data of the cancer genome atlas (TCGA) were utilized to evaluate the changes in the expression of LINC01929 in various cancers, as well as its relationship with the patients' survival rate. The co-expression of the genes and data merging of TCGA were utilized to identify the LINC01929-associated pathways. The samples of colorectal, gastric, and breast cancers were also examined by the RT-qPCR to confirm the results and evaluate the expression of LINC01929 in the mentioned cancers. In silico investigations indicated a remarkable enhancement in the expression of LINC01929 within the tumor tissues compared to normal samples in 10 types of cancer. Based on the survival results, the increase in the LINC01929 expression is linked to poor prognosis of bladder, breast, colorectal, kidney, and liver cancers. The gene co-expression network showed the strong co-expression of LINC01929 with genes involved in the metastatic pathways including COL5A1. RT-qPCR findings showed a remarkable increment in the expression level of LINC01929 in the colorectal, gastric, and breast tumor tissues versus the adjacent normal tissues. A significant and strong relationship was also found between the expression of LINC01929 and COL5A1. This study indicated a significant enhancement in the expression level of LINC01929 in various cancer types, accompanied by the mortality rate. Moreover, LINC01929 exhibited a strong co-expression with the metastatic genes such as COL5A1. As an oncogene and regulator of the metastatic pathways, LINC01929 can be a proper candidate for diagnostic and therapeutic purposes.

Evaluation of differentially expressed genes during replication using gene expression landscape of monkeypox-infected MK2 cells: A bioinformatics and systems biology approach to understanding the genomic pattern of viral replication

PURPOSE

The current outbreak of monkeypox (MPX) has created colossal concerns. However, immense research gaps have been noted in our understanding of the replication process, machinery, and genomic landscape during host cell infection. To fill this gap, differentially expressed genes (DEGs) were comprehensively analyzed during viral replication in host (MK2) cells.

METHODS

We used a microarray GEO dataset which was divided into three groups: control, MPXV-infected MK2 cells at 3 h, and MPXV-infected MK2 cells at 7 h. Using the dataset, DEG analysis, PPI network analysis, co-expression, and pathway analysis were conducted using bioinformatics, systems biology, and statistical approaches.

RESULTS

We identified 250 DEGs and 24 top-ranked genes. During the DEG analysis, we identified eight up-regulated genes (LOC695323, TMEM107, LOC695427, HIST1H2AD, LOC705469, PMAIP1, HIST1H2BJ, and HIST1H3D) and 16 down-regulated genes (HOXA9, BAMBI, LMO4, PAX6, AJUBA, CREBRF, CD24, JADE1, SLC7A11, EID2, SOX4, B4GALT5, PPARGC1A, BUB3, SOS2, and CDK19). We also developed PPI networks and performed co-expression analyses using the top-ranked genes. Furthermore, five genes were listed for co-expression pattern analysis.

CONCLUSIONS

This study will help in better understanding the replication process, machinery, and genomic landscape of the virus. This will further aid the discovery and development of therapeutics against viruses.

A novel process using immobilized engineered strain HC01 cells with co-expressed D-hydantoinase and D-N-carbamoylase as biocatalyst for the production of D-Valine

The demand for D-Valine increases because of its wide range of use. A whole-cell biocatalyst for the production of D-Valine from 5'-isopropyl hydantoin by co-expression of the D-hydantoinase (hyd) gene from Pseudomonas putida YZ-26 and D-N-carbamoylase (cab) gene from Sinorhizobium sp. SS-ori in Escherichia coli BL 21 (DE3) was developed. The expression condition of the engineered strain HC01 co-expressing D-hydantoinase (HYD) and D-N-carbamoylase (CAB) was optimized. HYD and CAB reached the highest activities (4.65 and 0.75 U/ml-broth) after inducing for 8-12 h. Subsequently, the cells of HC01 were immobilized in the form of Ca²⁺-alginate beads, and the optimal conditions for immobilizing were obtained as 2.5% gel concentration and 0.029 g/mL cell concentration in the presence of 3% CaCl₂. The thermostability of immobilized cells was 5 ℃ higher than that of free cells in the same condition. And the divalent metal ions such as Mn²⁺, Mg²⁺, Cu²⁺, Co²⁺, and Ni²⁺ did not significantly affect the enzymatic activity of HYD and CAB in immobilized cells. Bioconversion rate reached to 91% after a 42-h reaction when the substrate concentration was 50 mmol/L with the initial pH of 9.0 under the nitrogen protection. This method provides D-Valine with optical purity of 97% and an overall yield of 72%. Furthermore, the immobilized cells can be reused for more than 7 cycles and maintain their capacity of over 70%. Hence the immobilized cells of engineered strain HC01 could potentially be used to prepare D-Valine.

Co-expression of different proteins in Escherichia coli using plasmids with identical origins of replication

It is generally accepted that the use of two different plasmids with the identical origins of replication in bacteria is not desirable due to their "incompatibility". The utilization of the same bacterial enzymatic apparatus for replication of different plasmids is thought to cause a significant redistribution in favor of one of them. In the present work, examining co-expression of two different fluorescent proteins in Escherichia coli, we have shown that the use of highly homologous plasmids with identical origins of replication and providing resistance to different antibiotics results in high representation of both plasmids in bacteria. Meanwhile, the level of gene expression and the amount of proteins produced may differ and is determined mostly by their sequence rather than by the "incompatibility" of the plasmids.

Promoting soluble expression of hybrid mussel foot proteins by SUMO-TrxA tags for production of mussel glue

Mussel foot proteins (Mfps) display application potential with strong adhesion, enabling mussels to adhere firmly to various surfaces. Mytilus galloprovincialis foot protein 3B (Mgfp-3B) exhibits this characteristic remarkably. However, it remains a challenge for further research due to the low soluble expression of heterologous production. In this study, a small ubiquitin-related modifier (SUMO) and thioredoxin A (TrxA), which catalyzed the proper folding of disulfide bridges, were selected to increase the soluble expression of mfps. An additional ribosome binding site was introduced between the molecular chaperones and Mgfp-3B (fp-3) to form a bicistronic translation-coupled expression vector for co-expression. The results revealed that the combination of SUMO-TrxA increased the soluble expression of fp-3 by 18.07 %. Furthermore, the SUMO-TrxA also boosted the soluble expression of hybrid mfps Mgfp-3B-Mfp-1 (fp-3-1) by 11.29 %, Mgfp-3B-Mgfp-3B (fp-3-3) by 19.91 %, and Mgfp-3B-Mgfp-5 (fp-3-5) by 14.03 %. Ultimately, by high cell density cultivation in a 5 L bioreactor, the yields of fp-3, fp-3-3, and fp-3-5 co-expressed with SUMO-TrxA reached 217.75 mg/L, 127.2 mg/L, and 97.28 mg/L, respectively. Consequently, soluble production of mfps holds great potential for the sustainable supply of protein adhesive materials.

Bidirectional hybrid erythritol-inducible promoter for synthetic biology in Yarrowia lipolytica

BACKGROUND

The oleaginous yeast Yarrowia lipolytica is increasingly used as a chassis strain for generating bioproducts. Several hybrid promoters with different strengths have been developed by combining multiple copies of an upstream activating sequence (UAS) associated with a TATA box and a core promoter. These promoters display either constitutive, phase-dependent, or inducible strong expression. However, there remains a lack of bidirectional inducible promoters for co-expressing genes in Y. lipolytica.

RESULTS

This study built on our previous work isolating and characterizing the UAS of the erythritol-induced genes EYK1 and EYD1 (UAS-eyk1). We found an erythritol-inducible bidirectional promoter (BDP) located in the EYK1-EYL1 intergenic region. We used the BDP to co-produce YFP and RedStarII fluorescent proteins and demonstrated that the promoter's strength was 2.7 to 3.5-fold stronger in the EYL1 orientation compared to the EYK1 orientation. We developed a hybrid erythritol-inducible bidirectional promoter (HBDP) containing five copies of UAS-eyk1 in both orientations. It led to expression levels 8.6 to 19.2-fold higher than the native bidirectional promoter. While the BDP had a twofold-lower expression level than the strong constitutive TEF promoter, the HBDP had a 5.0-fold higher expression level when oriented toward EYL1 and a 2.4-fold higher expression level when oriented toward EYK1. We identified the optimal media for BDP usage by exploring yeast growth under microbioreactor conditions. Additionally, we constructed novel Golden Gate biobricks and a destination vector for general use.

CONCLUSIONS

In this research, we developed novel bidirectional and hybrid bidirectional promoters of which expression can be fine-tuned, responding to the need for versatile promoters in the yeast Y. lipolytica. This study provides effective tools that can be employed to smoothly adjust the erythritol-inducible co-expression of two target genes in biotechnology applications. BDPs developed in this study have potential applications in the fields of heterologous protein production, metabolic engineering, and synthetic biology.

Differential Expression Feature Extraction (DEFE): A Case Study in Wheat FHB RNA-Seq Data Analysis

In differential gene expression data analysis, one objective is to identify groups of co-expressed genes from a large dataset in order to detect the association between such a group of genes and an experimental condition. This is often done through a clustering approach, such as k-means or bipartition hierarchical clustering, based on particular similarity measures in the grouping process. In such a dataset, the gene differential expression itself is an innate attribute that can be used in the feature extraction process. For example, in a dataset consisting of multiple treatments versus their controls, the expression of a gene in each treatment would have three possible behaviors, upregulated, downregulated, or unchanged. We present in this chapter, a differential expression feature extraction (DEFE) method by using a string consisting of three numerical values at each character to denote such behavior, i.e., 1 = up, 2 = down, and 0 = unchanged, which results in up to 3^B differential expression patterns across all B comparisons. This approach has been successfully applied in many research projects, and among these, we demonstrate the strength of DEFE in a case study on RNA-sequencing (RNA-seq) data analysis of wheat challenged with the phytopathogenic fungus, Fusarium graminearum. Combinations of multiple schemes of DEFE patterns revealed groups of genes putatively associated with resistance or susceptibility to FHB.

Comprehensive expression analyses of the ABCG subfamily reveal SvABCG17 as a potential transporter of lignin monomers in the model C4 grass Setaria viridis

Although several aspects of lignin metabolism have been extensively characterized, the mechanism(s) by which lignin monomers are transported across the plasma membrane remains largely unknown. Biochemical, proteomic, expression and co-expression analyses from several plant species support the involvement of active transporters, mainly those belonging to the ABC superfamily. Here, we report on the genome-wide characterization of the ABCG gene subfamily in the model C4 grass Setaria viridis and further identification of the members potentially involved in monolignol transport. A total of 48 genes encoding SvABCGs were found in the S. viridis genome, from which 21 SvABCGs were classified as full-size transporters and 27 as half-size transporters. Comprehensive analysis of the ABCG subfamily in S. viridis based on expression and co-expression analyses support a role for SvABCG17 in monolignol transport: (i) SvABCG17 is orthologous to AtABCG29, a monolignol transporter in Arabidopsis thaliana; (ii) SvABCG17 displays a similar expression profile to that of lignin biosynthetic genes in a set of different S. viridis tissues and along the elongating internode; (iii) SvABCG17 is highly co-expressed with lignin-related genes in a public transcriptomic database; (iv) SvABCG17displays particularly high expression in the top of the S. viridis elongating internode, a tissue undergoing active lignification; (v) SvABCG17 mRNA localization coincides with the histochemical pattern of lignin deposition; and (vi) the promoter of SvABCG17 is activated by secondary cell wall-associated transcription factors, especially by lignin-specific activators of the MYB family. Further studies might reveal further aspects of this potential monolignol transporter, including its real substrate specificity and whether it works redundantly with other ABC members during S. viridis lignification.

Co-Expression of Chaperones for Improvement of Soluble Expression and Purification of An Anti-HER2 scFv in Escherichia Coli

Background

Single-chain fragment variable (scFv) is one of the most commonly used antibody fragments. They offer some advantages over full-length antibodies, including better penetration to target tissues. However, their functional production has been a challenge for manufacturers due to the potential misfolding and formation of inclusion bodies. Here we evaluated the soluble expression and purification of molecular chaperone co-expression.

Materials and Methods

E. coli BL21(DE3) cells were co-transformed with the mixture of plasmids pKJE7 and pET22b-scFv by the electroporation method. First, L-arabinose was added to induce the expression of molecular chaperones, and then IPTG was used as an inducer to start the expression of anti-HER2 scFv. The effect of cultivation temperature and IPTG concentration on soluble expression of the protein with or without chaperones was evaluated. The soluble expressed protein was subjected to native purification using the Ni-NTA affinity column.

Results

SDS-PAGE analysis confirmed the successful co-expression of anti-HER2-scFv and DnaK/DnaJ/GrpE chaperones. Co-expression with chaperones and low-temperature cultivation synergistically improved the soluble expression of anti-HER2 scFv. Co-expression with chaperone also exhibited an approximately four-fold increase in the final yield of purified soluble protein.

Conclusion

The combination of co-expression with chaperones and low temperature presented in this work may be useful for the improvement of commercial production of other scFvs in E. coli as functionally bioactive and soluble form.

Expression and Clinical Significance of Various Checkpoint Molecules in Advanced Osteosarcoma: Possibilities for Novel Immunotherapy

OBJECTIVES

The fact that studies on anti-programmed cell death 1 (PD-1) or its relevant ligand 1 (PD-L1) have yielded such few responses greatly decreases the confidence in immunotherapy with checkpoint inhibitors for advanced osteosarcoma. We intended to characterize the expression of various checkpoint molecules with immunohistochemistry in osteosarcoma specimens and analyzed the relationship of the expression of these checkpoint molecules with patients' clinical courses.

METHODS

This study was a retrospective non-intervention study from August 1st 2017 to March 1st 2020. Immunohistochemistry for B7-H3 (CD276, Cluster of Differentiation 276), CD47 (Cluster of Differentiation 47), PD-L1 (programmed cell death ligand 1), TIM3 (mucin-domain containing-3), TGF-β (TransformingGrowth Factor β), CXCR 4 (Chemokine Receptor 4), CD27 (Cluster of Differentiation 27), IDO1 (Indoleamine 2,3-dioxygenase 1), KIRs (Killer cell Immunoglobulin-like Receptors), and SDF-1 (Stromal cell-Derived Factor-1) was performed on 35 resected osteosarcoma specimens. Patients progressed upon first-line chemotherapy with evaluable lesions were qualified for this study, and their specimens previously stored in the pathological department repository would be retrieved for analysis. Associations between the immunohischemistry markers and clinicopathological variables and survival were evaluated by the χ² displayed by cross-table, Cox proportional hazards regression model, and Kaplan-Meier plots.

RESULTS

The positive rates of B7-H3, CD47, PD-L1, TIM3, and TGF-β expression in this sample of 35 heavily treated osteosarcomas were 29% (10/35), 15% (5/35), 9% (3/35), 6% (2/35), and 6% (2/35), respectively, and diverse staining intensities were observed. Among these advanced patients, 15/35 (43%) had positive checkpoint expression, of which 33% (5/15) showed evidence of the co-expression of more than one checkpoint molecule. We did not find any obvious correlation with clinicopathological characteristics and the positive expression of these molecules.

CONCLUSIONS

The present study highlights that only a small subset of progressive osteosarcomas, which had been heavily-treated, expressed tumor immune-associated checkpoint molecules, of which B7-H3 was the most positively expressed checkpoint and might be a promising target for further osteosarcoma investigation.

Coenzyme self-sufficiency system-recent advances in microbial production of high-value chemical phenyllactic acid

Phenyllactic acid (PLA), a natural antimicrobial substance, has many potential applications in the food, animal feed, pharmaceutical and cosmetic industries. However, its production is limited by the complex reaction steps involved in its chemical synthesis. Through advances in metabolic engineering and synthetic biology strategies, enzymatic or whole-cell catalysis was developed as an alternative method for PLA production. Herein, we review recent developments in metabolic engineering and synthetic biology strategies that promote the microbial production of high-value PLA. Specially, the advantages and disadvantages of the using of the three kinds of substrates, which includes phenylpyruvate, phenylalanine and glucose as starting materials by natural or engineered microbes is summarized. Notably, the bio-conversion of PLA often requires the consumption of expensive coenzyme NADH. To overcome the issues of NADH regeneration, efficiently internal cofactor regeneration systems constructed by co-expressing different enzyme combinations composed of lactate dehydrogenase with others for enhancing the PLA production, as well as their possible improvements, are discussed. In particular, the construction of fusion proteins with different linkers can achieve higher PLA yield and more efficient cofactor regeneration than that of multi-enzyme co-expression. Overall, this review provides a comprehensive overview of PLA biosynthesis pathways and strategies for increasing PLA yield through biotechnology, providing future directions for the large-scale commercial production of PLA and the expansion of downstream applications.

Introducing Protein Homeostasis and Glycogen Synthesis as Two Targets of Blue Light Radiation in Lentinula edodes

Introduction: There are documents about the biological effects of blue light radiation on different organisms. An understanding of the molecular mechanism of radiation effects on biological samples is an important event which has attracted researchers' attention. Determining the critical dysregulated proteins of Lentinula edodes following blue light radiation is the aim of this study. Methods: 22 differentially expressed proteins of L. edodes in response to 300 lux of blue light were extracted from the related literature. Experimental, text mining and co-expression connections between the queried proteins were assessed via the STRING database. The maps were compared and the critical proteins were identified. Results: Among the 21 queried proteins, six individuals including heat shock HSP70 protein, 20S proteasome subunit, 26S proteasome subunit P45, Aspartate aminotransferase, phosphopyruvate hydratase, and phosphoglucomutase were highlighted as the critical proteins in response to blue light radiation. Conclusion: The finding indicates that protein homeostasis and glycogen synthesis are affected by blue light radiation. Due to the critical roles of proteins as enzymes and structural elements in life maintenance and involvement of glycogen synthesis in energy consumption, blue light radiation can be considered as a life promotional agent in future investigations.

RiceNCexp: a rice non-coding RNA co-expression atlas based on massive RNA-seq and small-RNA seq data

Non-coding RNAs (ncRNAs) play important roles in regulating expression of protein-coding genes. Although gene expression databases have emerged in a timely manner, a comprehensive expression database for ncRNAs is still lacking. Herein, we constructed a rice ncRNA co-expression atlas (RiceNCexp), based on 491 RNA-seq and 274 small RNA (sRNA)-seq datasets. RiceNCexp hosts four types of ncRNAs, namely lncRNAs, PHAS genes, miRNAs, and phasiRNAs. RiceNCexp provides comprehensive expression information for rice ncRNAs in 22 tissues/organs, an efficient tau-based mining tool for tissue-specific ncRNAs, and the robust co-expression analysis among ncRNAs or between ncRNAs and protein-coding genes, based on 116 pairs of RNA-seq and sRNA-seq libraries from the same experiments. In summary, RiceNCexp is a user-friendly and comprehensive rice ncRNA co-expression atlas and can be freely accessed at https://cbi.njau.edu.cn/RiceNCexp/.

Improving AOX1 promoter efficiency by overexpression of Mit1 transcription factor

BACKGROUND

Reprogramming in transcriptional regulation provides an effective tool for adjusting cellular metabolic activities. The strong methanol-inducible alcohol oxidase-1 promoter (pAOX1) is commonly used for heterologous gene expression in the yeast Pichia pastoris. Here, we present a novel Pichia pastoris strain engineered to co-express methanol-induced transcription factor 1 (Mit1) and the target protein. Mit1 upregulates pAOX1 in response to methanol.

METHODS AND RESULTS

Two model proteins (VEGF and eGFP) have been used as the target proteins under the control of pAOX1. The sequence of Mit1 had obtained from the yeast genome and likewise cloned under the control of pAOX1. The results indicated a 1.9 and 2.2 fold increase in the detected VEGF and eGFP, respectively, when co-expressed with Mit1. Furthermore, the double-recombinant cells, containing Mit-1 and eGFP, produced 1.3 fold more eGFP when the methanol feeding concentration was doubled. The real-time PCR indicated a slight increase in the Mit1 expression, probably due to the negative regulatory feedback loop that exists for the intrinsic yeast Mit1. Overexpression of Mit1 also led to duplication of AOX1 enzyme activity, which may enhance the yeast cells' capacity for methanol detoxification.

CONCLUSION

Overexpression of Mit1 could be considered a promising strategy for upregulation of target recombinant proteins in Pichia pastoris. Intracellular overexpression of Mit1 upregulates the heterologous target gene (eGFP) production, which is expressed under the control of pAOX1.

Characterization of divergent promoters PmaiA and Phyd from Gordonia: Co-expression and regulation by CRP

Divergent promoters are often responsible for controlling gene expression of related genes of the same pathway or for coordinating regulation at different time points. There are relatively few reports on characterization of divergent promoters in bacteria. In the present study, microarray profiling was carried out to analyze gene expression during growth of Gordonia sp. IITR100, which led to the identification of 35 % of adjacent gene candidates that are divergently transcribed. We focus here on the in-depth characterization of one such pair of genes. Two divergent promoters, PmaiA and Phyd, drive the expression of genes encoding maleate cis-trans isomerase (maiA) and hydantoinase (hyd), respectively. Our findings reveal asymmetric promoter activity with higher activity in the reverse orientation (Phyd) as compared to the forward orientation (PmaiA). Minimal promoter region for each orientation was identified by deletion mapping. Deletion of a 5'-untranslated region of each gene resulted in an increase in promoter activity. A putative binding site for CRP (Catabolite Repressor Protein) transcription regulator was also identified in the 80 bp common regulatory region between the -35 hexamers of the two promoters. The results of this study suggest that CRP-mediated repression of PmaiA occurs only in the cells grown in glucose. Phyd, on the other hand, is not repressed by CRP. However, deletion of the CRP binding site located between -95 to -110 upstream to the transcription start site of the maiA gene resulted in increased activity of PmaiA and decreased activity of Phyd. A single CRP binding site, therefore, affects the two promoters differently.

Genome-wide identification of long non-coding (lncRNA) in Nilaparvata lugens's adaptability to resistant rice

Background

The brown planthopper (BPH), Nilaparvata lugens (Stål), is a very destructive pest that poses a major threat to rice plants worldwide. BPH and rice have developed complex feeding and defense strategies in the long-term co-evolution.

Methods

To explore the molecular mechanism of BPH's adaptation to resistant rice varieties, the lncRNA expression profiles of two virulent BPH populations were analyzed. The RNA-seq method was used to obtain the lncRNA expression data in TN1 and YHY15.

Results

In total, 3,112 highly reliable lncRNAs in TN1 and YHY15 were identified. Compared to the expression profiles between TN1 and YHY15, 157 differentially expressed lncRNAs, and 675 differentially expressed mRNAs were identified. Further analysis of the possible regulation relationships between differentially expressed lncRNAs and differentially expressed mRNAs, identified three pair antisense targets, nine pair cis-regulation targets, and 3,972 pair co-expressed targets. Function enriched found arginine and proline metabolism, glutathione metabolism, and carbon metabolism categories may significantly affect the adaptability in BPH when it is exposed to susceptible and resistant rice varieties. Altogether, it provided scientific data for the study of lncRNA regulation of brown planthopper resistance to rice. These results are helpful in the development of new control strategies for host defense against BPH and breeding rice for high yield.

Transcriptome profiling and co-expression network analysis of lncRNAs and mRNAs in colorectal cancer by RNA sequencing

Background

Long non-coding RNAs (lncRNAs) are widely involved in the pathogenesis of cancers. However, biological roles of lncRNAs in occurrence and progression of colorectal cancer (CRC) remain unclear. The current study aimed to evaluate the expression pattern of lncRNAs and messenger RNAs (mRNAs).

Methods

RNA sequencing (RNA-Seq) in CRC tissues and adjacent normal tissues from 6 CRC patients was performed and functional lncRNA-mRNA co-expression network was constructed afterwards. Gene enrichment analysis was demonstrated using DAVID 6.8 tool. Reverse transcription quantitative polymerase chain reaction (RT-qPCR) was used to validate the expression pattern of differentially expressed lncRNAs. Pearson correlation analysis was applied to evaluate the relationships between selected lncRNAs and mRNAs.

Results

One thousand seven hundred and sixteenth differentially expressed mRNAs and 311 differentially expressed lncRNAs were screened out. Among these, 568 mRNAs were up-regulated while 1148 mRNAs down-regulated, similarly 125 lncRNAs were up-regulated and 186 lncRNAs down-regulated. In addition, 1448 lncRNA–mRNA co-expression pairs were screened out from 940,905 candidate lncRNA-mRNA pairs. Gene enrichment analysis revealed that these lncRNA-related mRNAs are associated with cell adhesion, collagen adhesion, cell differentiation, and mainly enriched in ECM-receptor interaction and PI3K-Akt signaling pathways. Finally, RT-qPCR results verified the expression pattern of lncRNAs, as well as the relationships between lncRNAs and mRNAs in 60 pairs of CRC tissues.

Conclusions

In conclusion, these results of the RNA-seq and bioinformatic analysis strongly suggested that the dysregulation of lncRNA is involved in the complicated process of CRC development, and providing important insight regarding the lncRNAs involved in CRC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09878-6.

Expression analysis of peptidergic enteroendocrine cells in the silkworm Bombyx mori

Enteroendocrine cells (ECs) in the insect midgut respond to physiological changes in the intestine by releasing multiple peptides to control food intake, gastrointestinal activity and systemic metabolism. Here, we performed a comprehensive mapping of ECs producing different regulatory peptides in the larval midgut of Bombyx mori. In total, we identified 20 peptide genes expressed in different ECs in specific regions of the midgut. Transcript-specific in situ hybridisation combined with antibody staining revealed approximately 30 subsets of ECs, each producing a unique peptide or a combination of several different peptides. Functional significance of this diversity and specific roles of different enteroendocrine peptides are largely unknown. Results of this study highlight the importance of the midgut as a major endocrine/paracrine source of regulatory molecules in insects and provide important information to clarify functions of ECs during larval feeding and development.

Evidence from Co-expression Analysis for the Involvement of Amidase and INS in the Tryptophan-Independent Pathway of IAA Synthesis in Arabidopsis

The reverse genetic approach has uncovered indole synthase (INS) as the first enzyme in the tryptophan (trp)-independent pathway of IAA synthesis. The importance of INS was reevaluated suggesting it may interact with tryptophan synthase B (TSB) and therefore involved in the trp-dependent pathway. Thus, the main aim of this study was to clarify the route of INS through the analysis of Arabidopsis genome. Analysis of the top 2000 co-expression gene lists in general and specific conditions shows that TSA is strongly positively co-expressed with TSB in general, hormone, and abiotic conditions with mutual ranks of 89, 38, and 180 respectively. Moreover, TSA is positively correlated with TSB (0.291). However, INS was not found in any of these coexpressed gene lists and negatively correlated with TSB (- 0.046) suggesting unambiguously that these two routes are separately and independently operated. So far, the remaining steps in the INS pathway have remained elusive. Among all enzymes reported to have a role in IAA synthesis, amidase was found to strongly positively co-expressed with INS in general and light conditions with mutual ranks of 116 and 141 respectively. Additionally, amidase1 was found to positively correlate with INS (0.297) and negatively coexpressed with TSB concluding that amidase may exclusively involve in the trp-independent pathway.

RNA-sequencing of human aortic valves identifies that miR-629-3p and TAGLN miRNA-mRNA pair involving in calcified aortic valve disease

This study aimed to uncover the microRNA and messenger RNA (miRNA/mRNA) interactions in the pathophysiological process of calcified aortic valve disease (CAVD) of the human aortic valve. RNA sequencing of six selected samples (3 healthy control samples vs. 3 CAVD samples) was performed to obtain mRNA and miRNA sequences, and differential expression (DE) analysis of miRNA and mRNAs was performed. To build a CAVD-specific miRNA-mRNA interactome, the upregulated mRNAs and downregulated miRNAs were selected, followed by the establishment of inverse DE of mRNA-miRNA co-expression network based on Pearson's correlation coefficient using miRanda in the R language software. Subsequently, pathway enrichment analysis was performed to elucidate CAVD-related pathways that were likely mediated by miRNA regulatory mechanisms. In addition, miRNAs with an mRNA correlation greater than 0.9 in the co-expression network were selected for anti-calcification verification in a CAVD cellular model. We identified 216 mRNAs (99 downregulated and 117 upregulated) and 602 miRNAs (371 downregulated and 231 upregulated) that were differentially expressed between CAVD and healthy aortic valves. After applying Pearson's correlation toward miRNA-mRNA targets, a regulatory network of 67 miRNAs targeting 76 mRNAs was created. The subsequent pathway enrichment analysis of these targeted mRNAs elucidated that genes within the focal adhesion pathway are likely mediated by miRNA regulatory mechanisms. The selected hsa-miR-629-3p and TAGLN pair exhibited anti-calcification effects on osteogenic differentiation-induced human aortic valve interstitial cells (hVICs). On integrating the miRNA and mRNA sequencing data for healthy aortic valves and those with CAVD, the CAVD-associated miRNA-mRNA interactome and related pathways were elucidated. Additional cell function data demonstrated anti-calcification effects of the selected hsa-miR-629-3p targeting TAGLN, validating that it is a potential therapeutic target for inhibiting CAVD.

Integrated analysis of the clinical consequence and associated gene expression of ALK in ALK-positive human cancers

Anaplastic lymphoma kinase (ALK) is a tyrosine kinase receptor that is genetically altered in several cancers, including NSCLC, melanoma, lymphoma, and other tumors. Although ALK is associated with various cancers, the relationship between ALK expression and patient prognosis in different cancers is poorly understood. Here, using multidimensional approaches, we revealed the correlation between ALK expression and the clinical outcomes of patients with LUAD, melanoma, OV, DLBC, AML, and BC. We analyzed ALK transcriptional expression, patient survival rate, genetic alteration, protein network, and gene and microRNA (miRNA) co-expression. Compared to that in normal tissues, higher ALK expression was found in LUAD, melanoma, and OV, which are associated with poor patient survival rates. In contrast, lower transcriptional expression was found to decrease the survival rate of patients with DLBC, AML, and BC. A total of 202 missense mutations, 17 truncating mutations, 7 fusions, and 3 in-frame mutations were identified. Further, 17 genes and 19 miRNAs were found to be exclusively co-expressed and echinoderm microtubule-associated protein-like 4 (EML4) was identified as the most positively correlated gene (log odds ratio >3). The gene ontology and signaling pathways of the genes co-expressed with ALK in these six cancers were also identified. Our findings offer a basis for ALK as a prognostic biomarker and therapeutic target in cancers, which will potentially contribute to precision oncology and assist clinicians in identifying suitable treatment options.

Genome-wide identification of R2R3-MYB gene family and association with anthocyanin biosynthesis in Brassica species

Brassica species include important oil crops and vegetables in the world. The R2R3-MYB gene participates in a variety of plant functions, including the activation or inhibition of anthocyanin biosynthesis. Although previous studies have reported its phylogenetic relationships, gene structures, and expression patterns in Arabidopsis, the number and sequence variation of this gene family in Brassica crops and its involvement in the natural quantitative variation in anthocyanin biosynthesis regulation are still largely unknown. In this study, by using whole genome sequences and comprehensive genome-wide comparative analysis among the six cultivated Brassica species, 2120 R2R3-MYB genes were identified in six Brassica species, in total These R2R3-MYB genes were phylogenetically clustered into 12 groups. The R2R3-MYB family between A and C subgenomes showed better collinearity than between B and C and between A and B. From comparing transcriptional changes of five Brassica species with the purple and green leaves for the detection of the R2R3-MYB genes associated with anthocyanin biosynthesis, 7 R2R3-MYB genes were co-differentially expressed. The promoter and structure analysis of these genes showed that some variations between non-coding region, but they were highly conserved at the protein level and spatial structure. Co-expression analysis of anthocyanin-related genes and R2R3-MYBs indicated that MYB90 was strongly co-expressed with TT8, and they were co-expressed with structural genes F3H, LDOX, ANS and UF3GT at the same time. These results further clarified the roles of the R2R3-MYBs for leaf coloration in Brasica species, which provided new insights into the functions of the R2R3-MYB gene family in Brasica species.

Transcriptome analysis reveals multiple effects of nitrogen accumulation and metabolism in the roots, shoots, and leaves of potato (Solanum tuberosum L.)

BACKGROUND

Nitrogen (N) is a major element and fundamental constituent of grain yield. N fertilizer plays an essential role in the roots, shoots, and leaves of crop plants. Here, we obtained two N-sensitive potato cultivars.

RESULTS

The plants were cultivated in the pots using N-deficient and N-sufficient conditions. Crop height, leaf chlorophyll content, dry matter, and N-accumulation significantly decreased under N-deficient conditions. Furthermore, we performed a comprehensive analysis of the phenotype and transcriptome, GO terms, and KEGG pathways. We used WGCNA of co-expressed genes, and 116 differentially expressed hub genes involved in photosynthesis, nitrogen metabolism, and secondary metabolites to generate 23 modules. Among those modules, six NRT gene families, four pigment genes, two auxin-related genes, and two energy-related genes were selected for qRT-PCR validation.

CONCLUSIONS

Overall, our study demonstrates the co-expressed genes and potential pathways associated with N transport and accumulation in potato cultivars' roots, shoots, and leaves under N-deficient conditions. Therefore, this study provides new ideas to conduct further research on improving nitrogen use efficiency in potatoes.

Long Non-coding RNA and mRNA Co-expression Network Reveals Novel Players in Pleomorphic Xanthoastrocytoma

Histological interpretation of the rare pleomorphic xanthoastrocytoma (PXA) has been the holy grail for treatment options. However, no stand-alone clinical interventions have been developed owing to the lack of gene expression profiling data in PXA/APXA patients. We first time report the comprehensive analyses of the coding as well as long non-coding RNA (lncRNA) signatures of PXA/APXA patients. Several genes such as IGFBP2, NF1, FOS, ERBB2, and lncRNAs such as NEAT1, HOTAIRM1, and GAS5 known to play crucial roles in glioma patients were also deregulated in PXA patients suggesting the commonality in the molecular signatures. PPI network, co-expression, and lncRNA-mRNA interaction studies unraveled hub genes (such as ERBB2, FOS, RPA1) and networks that may play a critical role in PXA biology. The most enriched pathways based on gene profiles were related to TLR, chemokine, MAPK, Rb, and PI3K-Akt signaling pathways. The lncRNA targets were enriched in glucuronidation, adipogenesis, TGF-beta signaling, EGF/EGFR signaling, and cell cycle pathways. Interestingly, several mRNAs like PARVG, and ABI2 were found to be targeted by multiple lncRNAs suggesting a tight control of their levels. Some of the most prominent lncRNA-mRNA pairs were LOC728730: MRPL9, XLOC_l2_011987: ASIC2, lnc-C1QTNF5-1: RNF26. Notably, several lncRNAs such as lnc-CETP-1, lnc-XRCC3-1, lnc-RPL31-1, lnc-USP13-1, and MAPKAPK5-AS1, and genes such as RPA1, NTRK3, and CNRP1 showed strong correlation to the progression-free survival of PXA patients suggesting their potential as novel biomarkers. Overall, the findings of this study may facilitate the development of a new realm of RNA biology in PXA that may have clinical significance in the future.

Characterization of virus-like particles assembled by co-expression of BmCPV capsid shell protein and large protrusion protein

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) is a typical single-layer capsid dsRNA virus belonging to the Reoviridae family of the Cypovirus genus. Previous studies have shown that the BmCPV major capsid shell protein (CSP) has the ability to self-assemble into virus-like particles (VLPs), and cryo-electron microscopy of the BmCPV virions has revealed a tight mutual binding region between CSP and another capsid protein known as the Large Protrusion Protein (LPP), which further stabilizes the capsid shell. In this study, the multi-gene baculovirus expression system, Ac-MultiBac, was used to produce both solely CSP-based and CSP-LPP co-assembled VLPs. Transmission electron microscopy (TEM) results showed that addition of LPP did not affect the assembly of VLPs resulting in almost identical structure in both cases. However, ex vivo administration of VLPs to silkworm midgut tissue showed that CSP-based VLPs did not induce a significant transcriptional response in the innate immunity and RNAi gene cascades, compared to the co-assembled CSP-LPP based VLPs and the natural BmCPV virions isolated from polyhedra. The experimental results indicate that CSP and LPP attach tightly ("Plug and Display" model with CSP acting as "catcher" and LPP as "tag") to form VLPs that have a structure similar to that of the native CPV virions. Moreover, our results showed that the formation of VLPs with the two BmCPV capsid proteins is feasible, which can form the basis for the production of BmCPV-based VLPs as a new type of biological material to display exogenous proteins.

Gene Co-Expression Network Analysis Reveals the Correlation Patterns Among Genes in Different Temperature Stress Adaptation of Manila Clam

The Manila clam (Ruditapes philippinarum) is one of the most important aquaculture species and widely distributed along the coasts of China, Japan, and Korea. Due to its wide distribution, it can tolerate a wide range of temperature. Studying the gene expression profiles of clam gills had found differentially expressed genes (DEGs) and pathway involved in temperature stress tolerance. A systematic study of cellular response to temperature stress may provide insights into the mechanism of acquired tolerance. Here, weighted gene co-expression network analysis (WGCNA) was carried out using RNA-seq data from gill transcriptome in response to high and low temperature stress. There are a total 32 gene modules, of which 18 gene modules were identified as temperature-related modules. Blue module was one significantly correlated with temperature which was associated with cellular metabolism, apoptosis pathway, ER stress, and others.

Comparative transcriptomic analysis reveals region-specific expression patterns in different beef cuts

BACKGROUND

Beef cuts in different regions of the carcass have different meat quality due to their distinct physiological function. The objective of this study was to characterize the region-specific expression differences using comparative transcriptomics analysis among five representative beef cuts (tenderloin, longissimus lumborum, rump, neck, chuck).

RESULTS

We obtained 15,701 expressed genes in 30 muscle samples across five regions from carcass meat. We identified a total of 80 region-specific genes (RSGs), ranging from three (identified in the rump cut) to thirty (identified in the longissimus lumborum cut), and detected 25 transcription factors (TFs) for RSGs. Using a co-expression network analysis, we detected seven region-specific modules, including three positively correlated modules and four negatively correlated modules. We finally obtained 91 candidate genes related to meat quality, and the functional enrichment analyses showed that these genes were mainly involved in muscle fiber structure (e.g., TNNI1, TNNT1), fatty acids (e.g., SCD, LPL), amino acids (ALDH2, IVD, ACADS), ion channel binding (PHPT1, SNTA1, SUMO1, CNBP), protein processing (e.g., CDC37, GAPDH, NRBP1), as well as energy production and conversion (e.g., ATP8, COX8B, NDUFB6). Moreover, four candidate genes (ALDH2, CANX, IVD, PHPT1) were validated using RT-qPCR analyses which further supported our RNA-seq results.

CONCLUSIONS

Our results provide valuable insights into understanding the transcriptome regulation of meat quality in different beef cuts, and these findings may further help to improve the selection for health-beneficial meat in beef cattle.

Engineering cascade biocatalysis in whole cells for bottom-up synthesis of cello-oligosaccharides: flux control over three enzymatic steps enables soluble production

BACKGROUND

Soluble cello-oligosaccharides (COS, β-1,4-D-gluco-oligosaccharides with degree of polymerization DP 2-6) have been receiving increased attention in different industrial sectors, from food and feed to cosmetics. Development of large-scale COS applications requires cost-effective technologies for their production. Cascade biocatalysis by the three enzymes sucrose-, cellobiose- and cellodextrin phosphorylase is promising because it enables bottom-up synthesis of COS from expedient substrates such as sucrose and glucose. A whole-cell-derived catalyst that incorporates the required enzyme activities from suitable co-expression would represent an important step towards making the cascade reaction fit for production. Multi-enzyme co-expression to reach distinct activity ratios is challenging in general, but it requires special emphasis for the synthesis of COS. Only a finely tuned balance between formation and elongation of the oligosaccharide precursor cellobiose results in the desired COS.

RESULTS

Here, we show the integration of cellodextrin phosphorylase into a cellobiose-producing whole-cell catalyst. We arranged the co-expression cassettes such that their expression levels were upregulated. The most effective strategy involved a custom vector design that placed the coding sequences for cellobiose phosphorylase (CbP), cellodextrin phosphorylase (CdP) and sucrose phosphorylase (ScP) in a tricistron in the given order. The expression of the tricistron was controlled by the strong T7_lacO promoter and strong ribosome binding sites (RBS) for each open reading frame. The resulting whole-cell catalyst achieved a recombinant protein yield of 46% of total intracellular protein in an optimal ScP:CbP:CdP activity ratio of 10:2.9:0.6, yielding an overall activity of 315 U/g dry cell mass. We demonstrated that bioconversion catalyzed by a semi-permeabilized whole-cell catalyst achieved an industrial relevant COS product titer of 125 g/L and a space-time yield of 20 g/L/h. With CbP as the cellobiose providing enzyme, flux into higher oligosaccharides (DP ≥ 6) was prevented and no insoluble products were formed after 6 h of conversion.

CONCLUSIONS

A whole-cell catalyst for COS biosynthesis was developed. The coordinated co-expression of the three biosynthesis enzymes balanced the activities of the individual enzymes such that COS production was maximized. With the flux control set to minimize the share of insolubles in the product, the whole-cell synthesis shows a performance with respect to yield, productivity, product concentration and quality that is promising for industrial production.

SNMRS: An advanced measure for Co-expression network analysis

The challenge of identifying modules in a gene interaction network is important for a better understanding of the overall network architecture. In this work, we develop a novel similarity measure called Scaling-and-Shifting Normalized Mean Residue Similarity (SNMRS), based on the existing NMRS technique [1]. SNMRS yields correlation values in the range of 0 to +1 corresponding to negative and positive dependency. To study the performance of our measure, internal validation of extracted clusters resulting from different methods is carried out. Based on the performance, we choose hierarchical clustering and apply the same using the corresponding dissimilarity (distance) values of SNMRS scores, and utilize a dynamic tree cut method for extracting dense modules. The modules are validated using a literature search, KEGG pathway analysis, and gene-ontology analyses on the genes that make up the modules. Moreover, our measure can handle absolute, shifting, scaling, and shifting-and-scaling correlations and provides better performance than several other measures in terms of cluster-validity indices. Also, SNMRS based module detection method results in interesting biologically relevant patterns from gene microarray and RNA-seq dataset. A set of crucial genes having high relevance with the ESCC are also identified.