Search

WUSCHEL-related homeobox (WOX) genes are a class of plant-specific transcription factors, regulating the development of multiple tissues. However, the genomic characterizations and expression patterns of WOX genes have not been analyzed in lotus. In this study, 15 NnWOX genes were identified based on the well-annotated reference genome of lotus. According to the phylogenetic analysis, the NnWOX genes were clustered into three clades, i.e., ancient clade, intermediate clade, and WUS clade. Except for the conserved homeobox motif, we further found specific motifs of NnWOX genes in different clades and divergence gene structures, suggesting their distinct functions. In addition, two NnWOX genes in the ancient clade have conserved expression patterns and other NnWOX genes exhibit different expression patterns in lotus tissues, suggesting a low level of functional redundancy in lotus WOX genes. Furthermore, we constructed the gene co-expression networks for each NnWOX gene. Based on weighted gene co-expression network analysis (WGCNA), ten NnWOX genes and their co-expressed genes were assigned to the modules that were significantly related to the cotyledon and seed coat. We further performed RT-qPCR experiments, validating the expression levels of ten NnWOX genes in the co-expression networks. Our study reveals comprehensive genomic features of NnWOX genes in lotus, providing a solid basis for further function studies. [ABSTRACT FROM AUTHOR]

DUSP4 is a member of the DUSP (dual-specificity phosphatase) subfamily that is selective to the mitogen-activated protein kinases (MAPK) and has been implicated in a range of biological processes and functions in Alzheimer's disease (AD). In this study, we utilized the stereotactic delivery of adeno-associated virus (AAV)-DUSP4 to overexpress DUSP4 in the dorsal hippocampus of 5xFAD and wildtype (WT) mice, then used mass spectrometry (MS)-based proteomics along with the label-free quantification to profile the proteome and phosphoproteome in the hippocampus. We identified protein expression and phosphorylation patterns modulated in 5xFAD mice and examined the sex-specific impact of DUSP4 overexpression on the 5xFAD proteome/phosphoproteome. In 5xFAD mice, a substantial number of proteins were up- or down-regulated in both male and female mice in comparison to age and sex-matched WT mice, many of which are involved in AD-related biological processes, such as activated immune response or suppressed synaptic activities. Many proteins in pathways, such as immune response were found to be suppressed in response to DUSP4 overexpression in male 5xFAD mice. In contrast, such a shift was absent in female mice. For the phosphoproteome, we detected an array of phosphorylation sites regulated in 5xFAD compared to WT and modulated via DUSP4 overexpression in each sex. Interestingly, 5xFAD- and DUSP4-associated phosphorylation changes occurred in opposite directions. Strikingly, both the 5xFAD- and DUSP4-associated phosphorylation changes were found to be mostly in neurons and play key roles in neuronal processes and synaptic functions. Site-centric pathway analysis revealed that both the 5xFAD- and DUSP4-associated phosphorylation sites were enriched for a number of kinase sets in females but only a limited number of sets of kinases in male mice. Taken together, our results suggest that male and female 5xFAD mice responded to DUSP4 overexpression via shared and sex-specific molecular mechanisms, which might underly similar reductions in amyloid pathology in both sexes while learning deficits were reduced in only females with DUSP4 overexpression. Finally, we validated our findings with the sex-specific AD-associated proteomes in human cohorts and further developed DUSP4-centric proteomic network models and signaling maps for each sex. [ABSTRACT FROM AUTHOR]

One of the most destructive diseases, Gibberella stalk rot (GSR), caused by Fusarium graminearum, reduces maize yields significantly. An induced resistance response is a potent and cost-effective plant defense against pathogen attack. The functional counterpart of JAs, coronatine (COR), has attracted a lot of interest recently due to its ability to control plant growth and stimulate secondary metabolism. Although several studies have focused on COR as a plant immune elicitor to improve plant resistance to pathogens, the effectiveness and underlying mechanisms of the suppressive ability against COR to F. graminearum in maize have been limited. We investigated the potential physiological and molecular mechanisms of COR in modulating maize resistance to F. graminearum. COR treatment strongly enhanced disease resistance and promoted stomatal closure with H2O2 accumulation, and 10 μg/mL was confirmed as the best concentration. COR treatment increased defense-related enzyme activity and decreased the malondialdehyde content with enhanced antioxidant enzyme activity. To identify candidate resistance genes and gain insight into the molecular mechanism of GSR resistance associated with COR, we integrated transcriptomic and metabolomic data to systemically explore the defense mechanisms of COR, and multiple hub genes were pinpointed using weighted gene correlation network analysis (WGCNA). We discovered 6 significant modules containing 10 candidate genes: WRKY transcription factor (LOC100279570), calcium-binding protein (LOC100382070), NBR1-like protein (LOC100275089), amino acid permease (LOC100382244), glutathione S-transferase (LOC541830), HXXXD-type acyl-transferase (LOC100191608), prolin-rich extensin-like receptor protein kinase (LOC100501564), AP2-like ethylene-responsive transcription factor (LOC100384380), basic leucine zipper (LOC100275351), and glycosyltransferase (LOC606486), which are highly correlated with the jasmonic acid–ethylene signaling pathway and antioxidants. In addition, a core set of metabolites, including alpha-linolenic acid metabolism and flavonoids biosynthesis linked to the hub genes, were identified. Taken together, our research revealed differentially expressed key genes and metabolites, as well as co-expression networks, associated with COR treatment of maize stems after F. graminearum infection. In addition, COR-treated maize had higher JA (JA-Ile and Me-JA) levels. We postulated that COR plays a positive role in maize resistance to F. graminearum by regulating antioxidant levels and the JA signaling pathway, and the flavonoid biosynthesis pathway is also involved in the resistance response against GSR. [ABSTRACT FROM AUTHOR]

Aquaporins play a crucial role in helping water molecules move across cell membranes. While some studies have examined the role of AQPs in model plants like Arabidopsis, their impact on the ability of non-model plants to withstand environmental stress is largely unknown. In this study, we have explored the functions of the PaPIP1-2 gene, which encodes a protein called PIP, in apricot kernels (Prunus armeniaca L.). Our findings reveal how the PaPIP1-2 gene behaves during both dormancy and sprouting phases. Using a network analysis, we identified its interaction with six genes related to cold resistance. The mRNA levels of PaAQP genes, which co-express with cold resistance genes, remain consistent throughout different stages of P. armeniaca flower bud development, including physiological dormancy (PD), ecological dormancy (ED), sprouting period (SP), and germination stage (GS). Furthermore, our investigation of the location of the GFP-tagged PaPIP1-2 protein showed that it is mainly found in the cell membrane. Yeast strains with overexpressed PaPIP1-2 exhibited improved cold resistance and higher protein content. Similarly, when we overexpressed PaPIP1-2 in Arabidopsis, it enhanced the growth of these transgenic plants under cold stress. This improvement was associated with reduced levels of MDA (malondialdehyde); decreased ion leakage; increased proline accumulation; superoxide dismutase (SOD) activity; and the expression of cold resistance genes like AtPUB26, AtBTF3L, AtEBF1-1, and AtRAV1, compared with the wild-type plants. In summary, our results highlight the role of the P. armeniaca PaPIP1-2 gene in enhancing cold resistance and its importance in the dormancy and germination stages. [ABSTRACT FROM AUTHOR]

The BBAA subgenomes of hexaploid common wheat can be 'extracted' to constitute a viable and self-reproducing novel tetraploid wheat, termed extracted tetraploid wheat (ETW). Prior studies have shown ETW manifesting phenotypic abnormalities and alteration in gene expression and epigenetic modifications. No population level investigation has been conducted, leaving the issue unclear regarding whether developmental stability, an essential property evolved in all natural organisms, might have been undermined in ETW. Here, we measured variations in five morphological traits and somatic chromosomal stability in populations of ETW and of its hexaploid donor, a resynthesized hexaploid and a natural tetraploid wheat. We observed phenotypic defects in ETW. Meanwhile, we documented much greater within-population variations in ETW than in the other wheat genotypes, most probably due to disrupted developmental stability in ETW. Also, somatic structural chromosome variations were detected only in ETW. Comparative transcriptome analyses indicated that the disrupted developmental stability of ETW is likely linked to massive dysregulation of genome-wide gene expression rather than to genetic mutations. Population network analysis of gene expression implicated intrinsic connectivity among the variable traits, while gene set enrichment analysis provided possible links between dysregulated gene expression and interlaced trait variation. [ABSTRACT FROM AUTHOR]

Cultivated peanut (Arachis hypogaea L.) is an important economic and oilseed crop worldwide, providing high-quality edible oil and high protein content. Seed size/weight and oil content are two important determinants of yield and quality in peanut breeding. To identify key regulators controlling these two traits, two peanut cultivars with contrasting phenotypes were compared to each other, one having a larger seed size and higher oil content (Zhonghua16, ZH16 for short), while the second cultivar had smaller-sized seeds and lower oil content (Zhonghua6, ZH6). Whole transcriptome analyses were performed on these two cultivars at four stages of seed development. The results showed that ~40% of the expressed genes were stage-specific in each cultivar during seed development, especially at the early stage of development. In addition, we identified a total of 5356 differentially expressed genes (DEGs) between ZH16 and ZH6 across four development stages. Weighted gene co-expression network analysis (WGCNA) based on DEGs revealed multiple hub genes with potential roles in seed size/weight and/or oil content. These hub genes were mainly involved in transcription factors (TFs), phytohormones, the ubiquitin–proteasome pathway, and fatty acid synthesis. Overall, the candidate genes and co-expression networks detected in this study could be a valuable resource for genetic breeding to improve seed yield and quality traits in peanut. [ABSTRACT FROM AUTHOR]

Enhancing the absorption and utilization of phosphorus by crops is an important aim for ensuring food security worldwide. However, the gene regulatory network underlying phosphorus use in foxtail millet remains unclear. In this study, the molecular mechanism underlying low-phosphorus (LP) responsiveness in foxtail millet was evaluated using a comparative transcriptome analysis. LP reduced the chlorophyll content in shoots, increased the anthocyanin content in roots, and up-regulated purple acid phosphatase and phytase activities as well as antioxidant systems (CAT, POD, and SOD). Finally, 13 differentially expressed genes related to LP response were identified and verified using transcriptomic data and qRT-PCR. Two gene co-expression network modules related to phosphorus responsiveness were positively correlated with POD, CAT, and PAPs. Of these, SiPHR1, functionally annotated as PHOSPHATE STARVATION RESPONSE 1, was identified as an MYB transcription factor related to phosphate responsiveness. SiPHR1 overexpression in Arabidopsis significantly modified the root architecture. LP stress caused cellular, physiological, and phenotypic changes in seedlings. SiPHR1 functioned as a positive regulator by activating downstream genes related to LP tolerance. These results improve our understanding of the molecular mechanism underlying responsiveness to LP stress, thereby laying a theoretical foundation for the genetic modification and breeding of new LP-tolerant foxtail millet varieties. [ABSTRACT FROM AUTHOR]

The tea plant, Camellia sinensis (L.) O. Kuntze, is one of the most important beverage crops with significant economic and cultural value. Global climate change and population growth have led to increased salt and drought stress, negatively affecting tea yield and quality. The response mechanism of tea plants to these stresses remains poorly understood due to the lack of reference genome-based transcriptional descriptions. This study presents a high-quality genome-based transcriptome dynamic analysis of C. sinensis' response to salt and drought stress. A total of 2244 upregulated and 2164 downregulated genes were identified under salt and drought stress compared to the control sample. Most of the differentially expression genes (DEGs) were found to involve divergent regulation processes at different time points under stress. Some shared up- and downregulated DEGs related to secondary metabolic and photosynthetic processes, respectively. Weighted gene co-expression network analysis (WGCNA) revealed six co-expression modules significantly positively correlated with C. sinensis' response to salt or drought stress. The MEpurple module indicated crosstalk between the two stresses related to ubiquitination and the phenylpropanoid metabolic regulation process. We identified 1969 salt-responsive and 1887 drought-responsive allele-specific expression (ASE) genes in C. sinensis. Further comparison between these ASE genes and tea plant heterosis-related genes suggests that heterosis likely contributes to the adversity and stress resistance of C. sinensis. This work offers new insight into the underlying mechanisms of C. sinensis' response to salt and drought stress and supports the improved breeding of tea plants with enhanced salt and drought tolerance. [ABSTRACT FROM AUTHOR]

Leaf blight is a fungal disease that mainly affects the growth and development of leaves in plants. To investigate the molecular mechanisms of leaf blight defense in poplar, we performed RNA-Seq and enzyme activity assays on the Populus simonii × Populus nigra leaves inoculated with Alternaria alternate fungus. Through weighted gene co-expression network analysis (WGCNA), we obtained co-expression gene modules significantly associated with SOD and POD activities, containing 183 and 275 genes, respectively. We then constructed a co-expression network of poplar genes related to leaf blight resistance based on weight values. Additionally, we identified hub transcription factors (TFs) and structural genes in the network. The network was dominated by 15 TFs, and four out of them, including ATWRKY75, ANAC062, ATMYB23 and ATEBP, had high connectivity in the network, which might play important functions in leaf blight defense. In addition, GO enrichment analysis revealed a total of 44 structural genes involved in biotic stress, resistance, cell wall and immune-related biological processes in the network. Among them, there were 16 highly linked structural genes in the central part, which may be directly involved in poplar resistance to leaf blight. The study explores key genes associated with leaf blight defense in poplar, which further gains an understanding of the molecular mechanisms of biotic stress response in plants. [ABSTRACT FROM AUTHOR]

Aroma is a crucial attribute affecting the quality of pepper and its processed products, which has significant commercial value. However, little is known about the composition of volatile aroma compounds (VACs) in pepper fruits and their potential molecular regulatory mechanisms. In this study, HS-SPME-GC-MS combined with transcriptome sequencing is used to analyze the composition and formation mechanism of VACs in different kinds and development stages of pepper fruits. The results showed that 149 VACs, such as esters, alcohols, aldehydes, and terpenoids, were identified from 4 varieties and 3 development stages, and there were significant quantitative differences among different samples. Volatile esters were the most important aroma components in pepper fruits. PCA analysis showed that pepper fruits of different developmental stages had significantly different marker aroma compounds, which may be an important provider of pepper's characteristic aroma. Transcriptome analysis showed that many differential genes (DEGs) were enriched in the metabolic pathways related to the synthesis of VACs, such as fatty acids, amino acids, MVA, and MEP in pepper fruits. In addition, we identified a large number of differential transcription factors (TFs) that may regulate the synthesis of VACs. Combined analysis of differential aroma metabolites and DEGs identified two co-expression network modules highly correlated with the relative content of VACs in pepper fruit. This study confirmed the basic information on the changes of VACs in the fruits of several Chinese spicy peppers at different stages of development, screened out the characteristic aroma components of different varieties, and revealed the molecular mechanism of aroma formation, providing a valuable reference for the quality breeding of pepper. [ABSTRACT FROM AUTHOR]

Ultraviolet-B (UV-B) radiation is a major abiotic stress that dragon trees are exposed to during their growth and development; however, it is also an environmental signal perceived by plants that affects the flavonoid pathway. Previous studies have demonstrated that amounts of flavonoids are contained in dragon tree resin, otherwise known as dragon's blood. However, the traits and mechanism involved in the UV-B-mediated increase in flavonoids in dragon trees are still unknown. Here, we studied the response of Dracaena cambodiana under full solar UV-B radiation. The results showed that the contents of total flavonoids in D. cambodiana significantly increased after UV-B radiation exposure. Then, the transcriptome was used for determining the interactive mechanism of flavonoid accumulation and UV-B stress. Differential expression analyses identified 34 differentially expressed genes (DEGs) involved in flavonoid synthesis; specifically, 24 of the identified DEGs were significantly up-regulated after UV-B radiation exposure. In addition, 57 DEGs involved in Ca2+/kinase sensors, 58 DEGs involved in ROS scavenging and the plant hormone pathway, and 116 DEGs transcription factors in 5 families were further identified and analyzed. Finally, we deduced the potential mechanism of UV-B-promoting flavonoid formation to neutralize ROS damage derived from UV-B radiation in D. cambodiana based on the gene co-expression network and previous studies from other plants. Considering that wild dragon tree populations are currently highly threatened by anthropogenic and natural stressors, the interactive studies between D. cambodiana plants and UV-B radiation provide valuable information toward understanding the mechanism of dragon's blood formation and help us reveal the evolution of D. cambodiana, with the eventual goal of aiding in the global conservation of this precious biological resource. [ABSTRACT FROM AUTHOR]

Drought is one of the most serious abiotic stressors in the environment, restricting agricultural production by reducing plant growth, development, and productivity. To investigate such a complex and multifaceted stressor and its effects on plants, a systems biology-based approach is necessitated, entailing the generation of co-expression networks, identification of high-priority transcription factors (TFs), dynamic mathematical modeling, and computational simulations. Here, we studied a high-resolution drought transcriptome of Arabidopsis. We identified distinct temporal transcriptional signatures and demonstrated the involvement of specific biological pathways. Generation of a large-scale co-expression network followed by network centrality analyses identified 117 TFs that possess critical properties of hubs, bottlenecks, and high clustering coefficient nodes. Dynamic transcriptional regulatory modeling of integrated TF targets and transcriptome datasets uncovered major transcriptional events during the course of drought stress. Mathematical transcriptional simulations allowed us to ascertain the activation status of major TFs, as well as the transcriptional intensity and amplitude of their target genes. Finally, we validated our predictions by providing experimental evidence of gene expression under drought stress for a set of four TFs and their major target genes using qRT-PCR. Taken together, we provided a systems-level perspective on the dynamic transcriptional regulation during drought stress in Arabidopsis and uncovered numerous novel TFs that could potentially be used in future genetic crop engineering programs. [ABSTRACT FROM AUTHOR]

Salinity and alkalinity are among the important factors affecting the distribution, survival, growth and physiology of aquatic animals. Chinese sea bass (Lateolabrax maculatus) is an important aquaculture fish species in China that can widely adapt to diverse salinities from freshwater (FW) to seawater (SW) but moderately adapt to highly alkaline water (AW). In this study, juvenile L. maculatus were exposed to salinity change (SW to FW) and alkalinity stress (FW to AW). Coordinated transcriptomic responses in L. maculatus gills were investigated and based on the weighted gene co-expression network analysis (WGCNA), 8 and 11 stress-responsive modules (SRMs) were identified for salinity change and alkalinity stress, respectively, which revealed a cascade of cellular responses to oxidative and osmotic stress in L. maculatus gills. Specifically, four upregulated SRMs were enriched with induced differentially expressed genes (DEGs) for alkalinity stress, mainly corresponding to the functions of "extracellular matrix" and "anatomical structure", indicating a strong cellular response to alkaline water. Both "antioxidative activity" and "immune response" functions were enriched in the downregulated alkaline SRMs, which comprised inhibited alkaline specific DEGs, revealing the severely disrupted immune and antioxidative functions under alkalinity stress. These alkaline-specific responses were not revealed in the salinity change groups with only moderately inhibited osmoregulation and induced antioxidative response in L. maculatus gills. Therefore, the results revealed the diverse and correlated regulation of the cellular process and stress response in saline-alkaline water, which may have arisen through the functional divergence and adaptive recruitment of the co-expression genes and will provide vital insights for the development of L. maculatus cultivation in alkaline water. [ABSTRACT FROM AUTHOR]

Subtilases (SBTs), which belong to the serine peptidases, control plant development by regulating cell wall properties and the activity of extracellular signaling molecules, and affect all stages of the life cycle, such as seed development and germination, and responses to biotic and abiotic environments. In this study, 146 Gossypium hirsutum, 138 Gossypium barbadense, 89 Gossypium arboreum and 84 Gossypium raimondii SBTs were identified and divided into six subfamilies. Cotton SBTs are unevenly distributed on chromosomes. Synteny analysis showed that the members of SBT1 and SBT4 were expanded in cotton compared to Arabidopsis thaliana. Co-expression network analysis showed that six Gossypium arboreum SBT gene family members were in a network, among which five SBT1 genes and their Gossypium hirsutum and Arabidopsis thaliana direct homologues were down-regulated by salt treatment, indicating that the co-expression network might share conserved functions. Through co-expression network and annotation analysis, these SBTs may be involved in the biological processes of auxin transport, ABA signal transduction, cell wall repair and root tissue development. In summary, this study provides valuable information for the study of SBT genes in cotton and excavates SBT genes in response to salt stress, which provides ideas for cotton breeding for salinity resistance. [ABSTRACT FROM AUTHOR]

Eggshell plays an essential role in preventing physical damage and microbial invasions. Therefore, the analysis of genetic regulatory mechanisms of eggshell quality deterioration during aging in laying hens is important for the biosecurity and economic performance of poultry egg production worldwide. This study aimed to compare the differences in the expression profiles of long non-coding RNAs (lncRNAs) and mRNAs between old and young laying hens by the method of high-throughput RNA sequencing to identify candidate genes associated with aging in the uterus of laying hens. Overall, we detected 176 and 383 differentially expressed (DE) lncRNAs and mRNAs, respectively. Moreover, functional annotation analysis based on the Gene Ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) databases revealed that DE-lncRNAs and DE-mRNAs were significantly enriched in "phosphate-containing compound metabolic process", "mitochondrial proton-transporting ATP synthase complex", "inorganic anion transport", and other terms related to eggshell calcification and cuticularization. Through integrated analysis, we found that some important genes such as FGF14, COL25A1, GPX8, and GRXCR1 and their corresponding lncRNAs were expressed differentially between two groups, and the results of quantitative real-time polymerase chain reaction (qPCR) among these genes were also in excellent agreement with the sequencing data. In addition, our study found that TCONS_00181492, TCONS_03234147, and TCONS_03123639 in the uterus of laying hens caused deterioration of eggshell quality in the late laying period by up-regulating their corresponding target genes FGF14, COL25A1, and GRXCR1 as well as down-regulating the target gene GPX8 by TCONS_01464392. Our findings will provide a valuable reference for the development of breeding programs aimed at breeding excellent poultry with high eggshell quality or regulating dietary nutrient levels to improve eggshell quality. [ABSTRACT FROM AUTHOR]

Glioblastoma (GBM) is a type of brain cancer that is typically very aggressive and difficult to treat. Glioblastoma cases have been reported to have increased during COVID-19. The mechanisms underlying this comorbidity, including genomic interactions, tumor differentiation, immune responses, and host defense, are not completely explained. Therefore, we intended to investigate the differentially expressed shared genes and therapeutic agents which are significant for these conditions by using in silico approaches. Gene expression datasets of GSE68848, GSE169158, and GSE4290 studies were collected and analyzed to identify the DEGs between the diseased and the control samples. Then, the ontology of the genes and the metabolic pathway enrichment analysis were carried out for the classified samples based on expression values. Protein–protein interactions (PPI) map were performed by STRING and fine-tuned by Cytoscape to screen the enriched gene module. In addition, the connectivity map was used for the prediction of potential drugs. As a result, 154 overexpressed and 234 under-expressed genes were identified as common DEGs. These genes were found to be significantly enriched in the pathways involved in viral diseases, NOD-like receptor signaling pathway, the cGMP-PKG signaling pathway, growth hormone synthesis, secretion, and action, the immune system, interferon signaling, and the neuronal system. STAT1, CXCL10, and SAMDL were screened out as the top 03 out of the top 10 most critical genes among the DEGs from the PPI network. AZD-8055, methotrexate, and ruxolitinib were predicted to be the possible agents for the treatment. The current study identified significant key genes, common metabolic signaling networks, and therapeutic agents to improve our perception of the common mechanisms of GBM–COVID-19. [ABSTRACT FROM AUTHOR]

Age is an important physiological factor that affects the metabolism and immune function of beef cattle. While there have been many studies using the blood transcriptome to study the effects of age on gene expression, few have been reported on beef cattle. To this end, we used the blood transcriptomes of Japanese black cattle at different ages as the study subjects and screened 1055, 345, and 1058 differential expressed genes (DEGs) in the calf vs. adult, adult vs. old, and calf vs. old comparison groups, respectively. The weighted co-expression network consisted of 1731 genes. Finally, blue, brown, and yellow age-specific modules were obtained, in which genes were enriched in signaling pathways related to growth and development and immune metabolic dysfunction, respectively. Protein-protein interaction (PPI) analysis showed gene interactions in each specific module, and 20 of the highest connectivity genes were chosen as potential hub genes. Finally, we identified 495, 244, and 1007 genes by exon-wide selection signature (EWSS) analysis of different comparison groups. Combining the results of hub genes, we found that VWF, PARVB, PRKCA, and TGFB1I1 could be used as candidate genes for growth and development stages of beef cattle. CORO2B and SDK1 could be used as candidate marker genes associated with aging. In conclusion, by comparing the blood transcriptome of calves, adult cattle, and old cattle, the candidate genes related to immunity and metabolism affected by age were identified, and the gene co-expression network of different age stages was constructed. It provides a data basis for exploring the growth, development, and aging of beef cattle. [ABSTRACT FROM AUTHOR]

Nitrogen metabolism (NM) plays an essential role in response to abiotic stresses for plants. Enzyme activities have been extensively studied for nitrogen metabolism-associated pathways, but the knowledge of nitrogen metabolism-associated genes involved in stress response is still limited, especially for rice. In this study, we performed the genome-wide characterization of the genes putatively involved in nitrogen metabolism. A total of 1110 potential genes were obtained to be involved in nitrogen metabolism from eight species (Arabidopsis thaliana (L.) Heynh., Glycine max (L.) Merr., Brassica napus L., Triticum aestivum L., Sorghum bicolor L., Zea mays L., Oryza sativa L. and Amborella trichopoda Baill.), especially 104 genes in rice. The comparative phylogenetic analysis of the superfamily revealed the complicated divergence of different NM genes. The expression analysis among different tissues in rice indicates the NM genes showed diverse functions in the pathway of nitrogen absorption and assimilation. Distinct expression patterns of NM genes were observed in rice under drought stress, heat stress, and salt stress, indicating that the NM genes play a curial role in response to abiotic stress. Most NM genes showed a down-regulated pattern under heat stress, while complicated expression patterns were observed for different genes under salt stress and drought stress. The function of four representative NM genes (OsGS2, OsGLU, OsGDH2, and OsAMT1;1) was further validated by using qRT-PCR analysis to confirm their responses to these abiotic stresses. Based on the predicted transcription factor binding sites (TFBSs), we built a co-expression regulatory network containing transcription factors (TFs) and NM genes, of which the constructed ERF and Dof genes may act as the core genes to respond to abiotic stresses. This study provides novel sights to the interaction between nitrogen metabolism and the response to abiotic stresses. [ABSTRACT FROM AUTHOR]

The lignified tissue in the secondary stem is the main source of wood. In this study, we applied RNA-Seq analysis to the poplar stems in three developmental stages, including primary stem (PS), transitional stem (TS), and secondary stem (SS), to identify a total of 2028 genes that were highly expressed in the SS. Gene annotation indicated that the functions of these genes are mainly involved in cell wall biosynthesis, xylem development, and programmed cell death (PCD) processes. Subsequently, we explored the expression pattern of these genes at various developmental stages in the horizontal direction of the wood by ASPwood. The expression of these genes was modularized and correlated with the percentage of lignified xylem, using weighted gene co-expression network analysis (WGCNA). Among the genes, as many as 690 were identified as directly associated with lignification in the SS. In addition, the gene promoter cis-elements and protein interactions were predicted by PlantRegMap and STRING, respectively. The results were introduced into a co-expression network to confirm their relationship. We eventually found 54 TFs dominating this network, of which ADOF1, ATMYB3, AtbZIP44 (Potri.005G231300), ANAC043, ATWRKY40, ATEBP (Potri.010G006800), ARF5, anac075, RAP2.1, ARF16, AT- HSFB3, Potri.014G050000 (from WRKY family), HAT22, AT-HSFB2B, and AtWRKY20 had extremely high connectivity, which may play an important role in the lignification of wood formation at secondary stages. [ABSTRACT FROM AUTHOR]

As one of the most important crops for oil, olive (Olea europaea L.) is well-known worldwide for its commercial product "virgin olive oil" containing high-content fatty acids and many secondary metabolites. The molecular mechanisms underlying the enhanced oil content in olive remain unclear. To further investigate the molecular mechanisms of olive oil biosynthesis, we selected two olive cultivars, i.e., Kalinjot (JZ) and Coratina (KLD), at three maturity stages (MI-1, MI-3, and MI-6) for transcriptomic analysis based on Nanopore sequencing. Significant differences were observed in oil content between JZ and KLD during three maturity stages. Enrichment analysis revealed significant enrichment of differentially expressed genes (DEGs) in metabolic pathways of photosynthesis, amino acid biosynthesis, response to stress, and energy metabolism, in particular, fatty acid metabolism. A total of 170 (31.54% of 539 genes involved in oil synthesis) DEGs were further investigated based on expression analysis to identify their molecular functions in oil biosynthesis in olive. A co-expression network based on 714 transcription factors and their targeted genes in oil biosynthesis was constructed. Our study provided novel experimental evidence to investigate the molecular mechanisms of olive oil biosynthesis and to improve the breeding of olive varieties with enhanced oil contents. [ABSTRACT FROM AUTHOR]

Bud dormancy is a genotype-dependent mechanism observed in Prunus species in which bud growth is inhibited, and the accumulation of a specific amount of chilling (endodormancy) and heat (ecodormancy) is necessary to resume growth and reach flowering. We analyzed publicly available transcriptome data from fifteen cultivars of four Prunus species (almond, apricot, peach, and sweet cherry) sampled at endo- and ecodormancy points to identify conserved genes and pathways associated with dormancy control in the genus. A total of 13,018 genes were differentially expressed during dormancy transitions, of which 139 and 223 were of interest because their expression profiles correlated with endo- and ecodormancy, respectively, in at least one cultivar of each species. The endodormancy-related genes comprised transcripts mainly overexpressed during chilling accumulation and were associated with abiotic stresses, cell wall modifications, and hormone regulation. The ecodormancy-related genes, upregulated after chilling fulfillment, were primarily involved in the genetic control of carbohydrate regulation, hormone biosynthesis, and pollen development. Additionally, the integrated co-expression network of differentially expressed genes in the four species showed clusters of co-expressed genes correlated to dormancy stages and genes of breeding interest overlapping with quantitative trait loci for bloom time and chilling and heat requirements. [ABSTRACT FROM AUTHOR]

The F-box family exists in a wide variety of plants and plays an extremely important role in plant growth, development and stress responses. However, systematic studies of F-box family have not been reported in populus trichocarpa. In the present study, 245 PtrFBX proteins in total were identified, and a phylogenetic tree was constructed on the basis of their C-terminal conserved domains, which was divided into 16 groups (A–P). F-box proteins were located in 19 chromosomes and six scaffolds, and segmental duplication was main force for the evolution of the F-box family in poplar. Collinearity analysis was conducted between poplar and other species including Arabidopsis thaliana, Glycine max, Anemone vitifolia Buch, Oryza sativa and Zea mays, which indicated that poplar has a relatively close relationship with G. max. The promoter regions of PtrFBX genes mainly contain two kinds of cis-elements, including hormone-responsive elements and stress-related elements. Transcriptome analysis indicated that there were 82 differentially expressed PtrFBX genes (DEGs), among which 64 DEGs were in the roots, 17 in the leaves and 26 in the stems. In addition, a co-expression network analysis of four representative PtrFBX genes indicated that their co-expression gene sets were mainly involved in abiotic stress responses and complex physiological processes. Using bioinformatic methods, we explored the structure, evolution and expression pattern of F-box genes in poplar, which provided clues to the molecular function of F-box family members and the screening of salt-tolerant PtrFBX genes. [ABSTRACT FROM AUTHOR]

Auxin receptor plays a significant role in the plant auxin signalling pathway in response to abiotic stress. Recently, we found that transgenic rice overexpressing ABP57 had higher drought tolerance than the wild-type cultivar, MR219, due to the fact of its enhanced leaf photosynthetic rate and yields under drought stress. We performed a microarray study on this line to investigate the underlying mechanisms contributing to the observed phenotype. After microarray data filtering, 3596 genes were subjected to modular gene co-expression network (mGCN) development using CEMiTool, an R package. We identified highly related genes in 12 modules that could act to specific responses towards drought or any of the abiotic stress types. Gene set enrichment and overrepresentation analyses for modules extracted two highly upregulated modules that are involved in drought-related biological processes such as transmembrane transport of metal ions and response to oxidative stress. Finally, 123 hub genes were identified in all modules after integrating co-expression information with physical interaction data. In addition, the interplay of significant pathways between the metabolism of chlorophyll and flavonoid and the signalling pathways of MAPK, IAA, and SA inferred the concurrent involvement of stress tolerance response. Collectively, our findings seek new future directions for breeding strategies in rice tolerant improvements. [ABSTRACT FROM AUTHOR]

Poplars have a strong response to light, and shade is one of the main environmental factors that limits the growth and development of poplars. Exploring the response mechanism of the developing xylem of poplar under shady conditions is of great reference significance for improving wood yields. In this study, three excellent hybrids of poplar (Populus euramericana 'Zhonglin46' (Pe), Populus deltoides '27-17' (Pd), Populus × 'Wq156' (Pw) were studied under shady conditions. Based on the phenotypic data and developing a xylem transcriptome analysis, the molecular mechanism of poplars' response to shade was preliminarily revealed, and the core regulatory genes responding to shade were identified by a weighted co-expression network analysis (WGCNA). The results showed that Pw growth was significantly affected by shade, while Pe growth was slightly affected by shade. An enrichment analysis of 13,675 differentially expressed genes (DEGs) found that shade affected the expression of genes related to the glutathione metabolic pathway. The WGCNA analysis identified two modules ("Brown" and "Purple") related to the shade response and discovered seven hub genes. These hub genes were related to xylem development, vascular cambium division, stomatal development, and phytochrome A signal transduction. These results provide important basic information for gaining insight into the molecular response to shade in different poplar hybrids. [ABSTRACT FROM AUTHOR]

Salt stress negatively impacts crop production worldwide. Genetic diversity among barley (Hordeum vulgare) landraces adapted to adverse conditions should provide a valuable reservoir of tolerance genes for breeding programs. To identify molecular and biochemical differences between barley genotypes, transcriptomic and antioxidant enzyme profiles along with several morpho-physiological features were compared between salt-tolerant (Boulifa) and salt-sensitive (Testour) genotypes subjected to salt stress. Decreases in biomass, photosynthetic parameters, and relative water content were low in Boulifa compared to Testour. Boulifa had better antioxidant protection against salt stress than Testour, with greater antioxidant enzymes activities including catalase, superoxide dismutase, and guaiacol peroxidase. Transcriptome assembly for both genotypes revealed greater accumulation of differentially expressed transcripts in Testour compared to Boulifa, emphasizing the elevated transcriptional response in Testour following salt exposure. Various salt-responsive genes, including the antioxidant catalase 3, the osmoprotectant betaine aldehyde dehydrogenase 2, and the transcription factors MYB20 and MYB41, were induced only in Boulifa. By contrast, several genes associated with photosystems I and II, and light receptor chlorophylls A and B, were more repressed in Testour. Co-expression network analysis identified specific gene modules correlating with differences in genotypes and morpho-physiological traits. Overall, salinity-induced differential transcript accumulation underlies the differential morpho-physiological response in both genotypes and could be important for breeding salt tolerance in barley. [ABSTRACT FROM AUTHOR]

NAC (NAM, ATAF1/2, and CUC2) transcription factors (TFs) are one of the largest plant-specific TF families and play a pivotal role in adaptation to abiotic stresses. The genome-wide analysis of NAC TFs is still absent in Zanthoxylum bungeanum. Here, 109 ZbNAC proteins were identified from the Z. bungeanum genome and were classified into four groups with Arabidopsis NAC proteins. The 109 ZbNAC genes were unevenly distributed on 46 chromosomes and included 4 tandem duplication events and 17 segmental duplication events. Synteny analysis of six species pairs revealed the closely phylogenetic relationship between Z. bungeanum and C. sinensis. Twenty-four types of cis-elements were identified in the ZbNAC promoters and were classified into three types: abiotic stress, plant growth and development, and response to phytohormones. Co-expression network analysis of the ZbNACs revealed 10 hub genes, and their expression levels were validated by real-time quantitative polymerase chain reaction (qRT-PCR). Finally, ZbNAC007, ZbNAC018, ZbNAC047, ZbNAC072, and ZbNAC079 were considered the pivotal NAC genes for drought tolerance in Z. bungeanum. This study represented the first genome-wide analysis of the NAC family in Z. bungeanum, improving our understanding of NAC proteins and providing useful information for molecular breeding of Z. bungeanum. [ABSTRACT FROM AUTHOR]

Simple Summary: Hepatocellular carcinoma (HCC) is the most common malignancy of liver cancer. However, treatment of HCC is still severely limited due to limitation of drug therapy. We aimed to screen more possible target genes and candidate drugs for HCC, exploring the possibility of drug treatments from systems biological perspective. We identified ten candidate target genes, which are hub genes in HCC co-expression networks, which also possess significant prognostic value in two independent HCC cohorts. The rationality of these target genes was well demonstrated through variety analyses of patient expression profiles. We then screened candidate drugs for target genes and finally identified withaferin-a and mitoxantrone as the candidate drug for HCC treatment. The drug effectiveness was validated in in vitro model and computational analysis, providing more evidence for our drug repositioning method and results. Hepatocellular carcinoma (HCC) is a malignant liver cancer that continues to increase deaths worldwide owing to limited therapies and treatments. Computational drug repurposing is a promising strategy to discover potential indications of existing drugs. In this study, we present a systematic drug repositioning method based on comprehensive integration of molecular signatures in liver cancer tissue and cell lines. First, we identify robust prognostic genes and two gene co-expression modules enriched in unfavorable prognostic genes based on two independent HCC cohorts, which showed great consistency in functional and network topology. Then, we screen 10 genes as potential target genes for HCC on the bias of network topology analysis in these two modules. Further, we perform a drug repositioning method by integrating the shRNA and drug perturbation of liver cancer cell lines and identifying potential drugs for every target gene. Finally, we evaluate the effects of the candidate drugs through an in vitro model and observe that two identified drugs inhibited the protein levels of their corresponding target genes and cell migration, also showing great binding affinity in protein docking analysis. Our study demonstrates the usefulness and efficiency of network-based drug repositioning approach to discover potential drugs for cancer treatment and precision medicine approach. [ABSTRACT FROM AUTHOR]

Several species in Brassicaceae produce glucosinolates (GSLs) to protect themselves against pests. As demonstrated in A. thaliana, the reallocation of defence compounds, of which GSLs are a major part, is highly dependent on transport processes and serves to protect high-value tissues such as reproductive tissues. This study aimed to identify potential GSL-transporter proteins (TPs) using a network-biology approach. The known A. thaliana GSL genes were retrieved from the literature and pathway databases and searched against several co-expression databases to generate a gene network consisting of 1267 nodes and 14,308 edges. In addition, 1151 co-expressed genes were annotated, integrated, and visualised using relevant bioinformatic tools. Based on three criteria, 21 potential GSL genes encoding TPs were selected. The AST68 and ABCG40 potential GSL TPs were chosen for further investigation because their subcellular localisation is similar to that of known GSL TPs (SULTR1;1 and SULTR1;2) and ABCG36, respectively. However, AST68 was selected for a molecular-docking analysis using AutoDOCK Vina and AutoDOCK 4.2 with the generated 3D model, showing that both domains were well superimposed on the homologs. Both molecular-docking tools calculated good binding-energy values between the sulphate ion and Ser419 and Val172, with the formation of hydrogen bonds and van der Waals interactions, respectively, suggesting that AST68 was one of the sulphate transporters involved in GSL biosynthesis. This finding illustrates the ability to use computational analysis on gene co-expression data to screen and characterise plant TPs on a large scale to comprehensively elucidate GSL metabolism in A. thaliana. Most importantly, newly identified potential GSL transporters can serve as molecular tools in improving the nutritional value of crops. [ABSTRACT FROM AUTHOR]

Quinoa is a gluten-free food crop that contains all the essential amino acids and vitamins. The selection of proper housekeeping and tissue-specific genes is the crucial prerequisite for gene expression analysis using the common approach, real-time quantitative PCR (RT-qPCR). In this study, we identified 40 novel candidate housekeeping genes by the minimum transcript per million (TPM), coefficient of variation (CV) and maximum fold change (MFC) methods and 19 candidate tissuespecific genes by the co-expression network method based on an RNA-seq dataset that included 53 stem, leaf, flower and seed samples, as well as additional shoot and root samples under different stresses. The expression stability of 12 housekeeping and tissue-specific genes, as well as that of another two traditionally used housekeeping genes, was further evaluated using qPCR and ranked using NormFinder, BestKeeper and the comparative delta-Ct method. The results demonstrated that MIF, RGGA, VATE and UBA2B were ranked as the top four most stable candidate housekeeping genes. qPCR analysis also revealed three leaf-specific genes and five root-specific genes, but no stem-specific gene was identified. Gene Ontology (GO) enrichment analysis identified that housekeeping genes were mainly enriched in the small molecule metabolic process, organonitrogen compound metabolic process, NAD binding and ligase activity. In addition, tissue-specific genes are closely associated with the major functions of a specific tissue. Specifically, GO terms "photosynthesis" and "thylakoid" were most significantly overrepresented in candidate leaf-specific genes. The novel housekeeping and tissue-specific genes in our study will enable better normalization and quantification of transcript levels in quinoa. [ABSTRACT FROM AUTHOR]

Single-cell RNA sequencing is a powerful technology for obtaining transcriptomes at single-cell resolutions. However, it suffers from dropout events (i.e., excess zero counts) since only a small fraction of transcripts get sequenced in each cell during the sequencing process. This inherent sparsity of expression profiles hinders further characterizations at cell/gene-level such as cell type identification and downstream analysis. To alleviate this dropout issue we introduce a network-based method, netImpute, by leveraging the hidden information in gene co-expression networks to recover real signals. netImpute employs Random Walk with Restart (RWR) to adjust the gene expression level in a given cell by borrowing information from its neighbors in a gene co-expression network. Performance evaluation and comparison with existing tools on simulated data and seven real datasets show that netImpute substantially enhances clustering accuracy and data visualization clarity, thanks to its effective treatment of dropouts. While the idea of netImpute is general and can be applied with other types of networks such as cell co-expression network or protein–protein interaction (PPI) network, evaluation results show that gene co-expression network is consistently more beneficial, presumably because PPI network usually lacks cell type context, while cell co-expression network can cause information loss for rare cell types. Evaluation results on several biological datasets show that netImpute can more effectively recover missing transcripts in scRNA-seq data and enhance the identification and visualization of heterogeneous cell types than existing methods. [ABSTRACT FROM AUTHOR]

The mitogen-activated protein kinase (MAPK) cascade is a universal signal transduction module that plays a vital role in regulating growth and development, as well as environmental stress responses in plants. Wheat is one of the most important crops worldwide. Although the MAPK kinase kinase (MAP3K) family in wheat has been investigated, the MAPK and MAPK kinase (MAP2K) gene families remain unknown at present. Here, 54 MAPK and 18 MAPKK genes were identified in wheat using recent genomic information. Phylogenetic analysis of Triticum aestivum L. MAPKs and MAPKKs (TaMAPKs and TaMAPKKs) together with homologous genes from other species classified them into four groups, and the clustering was consistent with the genomic exon/intron structures. Conserved motif analysis found that MAPK proteins contained a typical TXY phosphorylation site and MAPKK proteins contained an S/T-X5-S/T motif. RNA-seq data mapping analysis showed that MAPK and MAPKK genes in group IV had tissue-specific expression profiles, whereas each group I member showed relatively high expression in all organs. Expression patterns of TaMAPK and TaMAPKK genes under stress conditions were also investigated and stress-responsive candidates were identified. Co-expression network analysis identified 11 TaMAPK genes and 6 TaMAPKK genes involved in the interaction network pathway. Overall, this study provided useful information for evolutionary and functional surveys of MAPK and MAPKK gene families in wheat and beyond. [ABSTRACT FROM AUTHOR]

Co-expression network analyses provide insights into the molecular interactions underlying complex traits and diseases. In this study, co-expression network analysis was performed to detect expression patterns (modules or clusters) of microRNAs (miRNAs) during lactation, and to identify miRNA regulatory mechanisms for milk yield and component traits (fat, protein, somatic cell count (SCC), lactose, and milk urea nitrogen (MUN)) via miRNA target gene enrichment analysis. miRNA expression (713 miRNAs), and milk yield and components (Fat%, Protein%, lactose, SCC, MUN) data of nine cows at each of six different time points (day 30 (D30), D70, D130, D170, D230 and D290) of an entire lactation curve were used. Four modules or clusters (GREEN, BLUE, RED and TURQUOISE) of miRNAs were identified as important for milk yield and component traits. The GREEN and BLUE modules were significantly correlated (|r| > 0.5) with milk yield and lactose, respectively. The RED and TURQUOISE modules were significantly correlated (|r| > 0.5) with both SCC and lactose. In the GREEN module, three abundantly expressed miRNAs (miR-148a, miR-186 and miR-200a) were most significantly correlated to milk yield, and are probably the most important miRNAs for this trait. DDR1 and DDHX1 are hub genes for miRNA regulatory networks controlling milk yield, while HHEX is an important transcription regulator for these networks. miR-18a, miR-221/222 cluster, and transcription factors HOXA7, and NOTCH 3 and 4, are important for the regulation of lactose. miR-142, miR-146a, and miR-EIA17-14144 (a novel miRNA), and transcription factors in the SMAD family and MYB, are important for the regulation of SCC. Important signaling pathways enriched for target genes of miRNAs of significant modules, included protein kinase A and PTEN signaling for milk yield, eNOS and Noth signaling for lactose, and TGF β, HIPPO, Wnt/β-catenin and cell cycle signaling for SCC. Relevant enriched gene ontology (GO)-terms related to milk and mammary gland traits included cell differentiation, G-protein coupled receptor activity, and intracellular signaling transduction. Overall, this study uncovered regulatory networks in which miRNAs interacted with each other to regulate lactation traits. [ABSTRACT FROM AUTHOR]

Nitrogen metabolism (NM) plays an essential role in response to abiotic stresses for plants. Enzyme activities have been extensively studied for nitrogen metabolism-associated pathways, but the knowledge of nitrogen metabolism-associated genes involved in stress response is still limited, especially for rice. In this study, we performed the genome-wide characterization of the genes putatively involved in nitrogen metabolism. A total of 1110 potential genes were obtained to be involved in nitrogen metabolism from eight species ( Arabidopsis thaliana (L.) Heynh., Glycine max (L.) Merr., Brassica napus L., Triticum aestivum L., Sorghum bicolor L., Zea mays L., Oryza sativa L. and Amborella trichopoda Baill.), especially 104 genes in rice. The comparative phylogenetic analysis of the superfamily revealed the complicated divergence of different NM genes. The expression analysis among different tissues in rice indicates the NM genes showed diverse functions in the pathway of nitrogen absorption and assimilation. Distinct expression patterns of NM genes were observed in rice under drought stress, heat stress, and salt stress, indicating that the NM genes play a curial role in response to abiotic stress. Most NM genes showed a down-regulated pattern under heat stress, while complicated expression patterns were observed for different genes under salt stress and drought stress. The function of four representative NM genes ( OsGS2 , OsGLU , OsGDH2 , and OsAMT1;1 ) was further validated by using qRT-PCR analysis to confirm their responses to these abiotic stresses. Based on the predicted transcription factor binding sites (TFBSs), we built a co-expression regulatory network containing transcription factors (TFs) and NM genes, of which the constructed ERF and Dof genes may act as the core genes to respond to abiotic stresses. This study provides novel sights to the interaction between nitrogen metabolism and the response to abiotic stresses.

Salinity is increasingly becoming a significant problem for the most important yet intrinsically salt-sensitive grain legume chickpea. Chickpea is extremely sensitive to salinity during the reproductive phase. Therefore, it is essential to understand the molecular mechanisms by comparing the transcriptomic dynamics between the two contrasting genotypes in response to salt stress. Chickpea exhibits considerable genetic variation amongst improved cultivars, which show better yields in saline conditions but still need to be enhanced for sustainable crop production. Based on previous extensive multi-location physiological screening, two identified genotypes, JG11 (salt-tolerant) and ICCV2 (salt-sensitive), were subjected to salt stress to evaluate their phenological and transcriptional responses. RNA-Sequencing is a revolutionary tool that allows for comprehensive transcriptome profiling to identify genes and alleles associated with stress tolerance and sensitivity. After the first flowering, the whole flower from stress-tolerant and sensitive genotypes was collected. A total of ~300 million RNA-Seq reads were sequenced, resulting in 2022 differentially expressed genes (DEGs) in response to salt stress. Genes involved in flowering time such as FLOWERING LOCUS T (FT) and pollen development such as ABORTED MICROSPORES (AMS), rho-GTPase, and pollen-receptor kinase were significantly differentially regulated, suggesting their role in salt tolerance. In addition to this, we identify a suite of essential genes such as MYB proteins, MADS-box, and chloride ion channel genes, which are crucial regulators of transcriptional responses to salinity tolerance. The gene set enrichment analysis and functional annotation of these genes in flower development suggest that they can be potential candidates for chickpea crop improvement for salt tolerance. [ABSTRACT FROM AUTHOR]

Black pepper (Piper nigrum L.), is dubbed "the King of Spices". However, the lack of genic knowledge has limited the understanding of its physiological processes and hindered the development of its molecular breeding. The SBP-box gene family is an important family in plant development and integrates multiple physiological processes. Here, we made a genome-wide identification of the pepper SBP-box gene family to provide evolutionary and functional information about this conserved transcription factor. In total, 34 SBP genes were identified in pepper. All these pepper SBP genes were clustered into eight groups, and one pepper group was not found in Arabidopsis thaliana. Segment duplications played the most important role in the expansion process of pepper SBP genes, and all these duplications were subjected to purifying selection. Half of pepper SBP genes were found miR156 target sites, and 17 miR156s were predicted. The tissue expression analysis revealed the differential expression of pepper SBP genes. Eleven SBP genes were found in four co-expression networks, and the GO enrichment further provides a functional prediction for pepper SBP genes. This study lays a foundation for further studies of pepper and provides a valuable reference for functional mining of pepper SBP genes. [ABSTRACT FROM AUTHOR]

Gallbladder cancer (GBC) has a lower incidence rate among the population relative to other cancer types but is a major contributor to the total number of biliary tract system cancer cases. GBC is distinguished from other malignancies by its high mortality, marked geographical variation and poor prognosis. To date no systemic targeted therapy is available for GBC. The main objective of this study is to determine the molecular signatures correlated with GBC development using integrative systems level approaches. We performed analysis of publicly available transcriptomic data to identify differentially regulated genes and pathways. Differential co-expression network analysis and transcriptional regulatory network analysis was performed to identify hub genes and hub transcription factors (TFs) associated with GBC pathogenesis and progression. Subsequently, we assessed the epithelial-mesenchymal transition (EMT) status of the hub genes using a combination of three scoring methods. The identified hub genes including, CDC6, MAPK15, CCNB2, BIRC7, L3MBTL1 were found to be regulators of cell cycle components which suggested their potential role in GBC pathogenesis and progression. [ABSTRACT FROM AUTHOR]

Once considered nonfunctional, pseudogene transcripts are now known to provide valuable information for cancer susceptibility, including head and neck cancer (HNC), a serious health problem worldwide, with about 50% unimproved overall survival over the last decades. The present review focuses on the role of pseudogene transcripts involved in HNC risk and prognosis. We combined current literature and in silico analyses from The Cancer Genome Atlas (TCGA) database to identify the most deregulated pseudogene transcripts in HNC and their genetic variations. We then built a co-expression network and performed gene ontology enrichment analysis to better understand the pseudogenes' interactions and pathways in HNC. In the literature, few pseudogenes have been studied in HNC. Our in silico analysis identified 370 pseudogene transcripts associated with HNC, where SPATA31D5P, HERC2P3, SPATA31C2, MAGEB6P1, SLC25A51P1, BAGE2, DNM1P47, SPATA31C1, ZNF733P and OR2W5 were found to be the most deregulated and presented several genetic alterations. NBPF25P, HSP90AB2P, ZNF658B and DPY19L2P3 pseudogenes were predicted to interact with 12 genes known to participate in HNC, DNM1P47 was predicted to interact with the TP53 gene, and HLA-H pseudogene was predicted to interact with HLA-A and HLA-B genes. The identified pseudogenes were associated with cancer biology pathways involving cell communication, response to stress, cell death, regulation of the immune system, regulation of gene expression, and Wnt signaling. Finally, we assessed the prognostic values of the pseudogenes with the Kaplan–Meier Plotter database, and found that expression of SPATA31D5P, SPATA31C2, BAGE2, SPATA31C1, ZNF733P and OR2W5 pseudogenes were associated with patients' survival. Due to pseudogene transcripts' potential for cancer diagnosis, progression, and as therapeutic targets, our study can guide new research to HNC understanding and development of new target therapies. [ABSTRACT FROM AUTHOR]

Systemic sclerosis (SSc) is a chronic autoimmune disease mainly affecting the connective tissue. In SSc patients, monocytes are increased in circulation, infiltrate affected tissues, and show a pro-inflammatory activation status, including the so-called interferon (IFN) signature. We previously demonstrated that the dysregulation of the IFN response in SSc monocytes is sustained by altered epigenetic factors as well as by upregulation of the long non-coding RNA (lncRNA) NRIR. Considering the enormously diverse molecular functions of lncRNAs in immune regulation, the present study investigated the genome-wide profile of lncRNAs in SSc monocytes, with the aim to further unravel their possible role in monocyte dysregulation and disease pathogenesis. Transcriptomic data from two independent cohorts of SSc patients identified 886 lncRNAs with an altered expression in SSc monocytes. Differentially expressed lncRNAs were correlated with neighboring protein coding genes implicated in the regulation of IFN responses and apoptotic signaling in SSc monocytes. In parallel, gene co-expression network analysis identified the lncRNA PSMB8-AS1 as a top-ranking hub gene in co-expression modules implicated in cell activation and response to viral and external stimuli. Functional characterization of PSMB8-AS1 in monocytes demonstrated that this lncRNA is involved in the secretion of IL-6 and TNFα, two pivotal pro-inflammatory cytokines altered in the circulation of SSc patients and associated with fibrosis and disease severity. Collectively, our data showed that lncRNAs are linked to monocyte dysregulation in SSc, and highlight their potential contribution to disease pathogenesis. [ABSTRACT FROM AUTHOR]

Hematopoietic differentiation is a well-orchestrated process by many regulators such as transcription factor and long non-coding RNAs (lncRNAs). However, due to the large number of lncRNAs and the difficulty in determining their roles, the study of lncRNAs is a considerable challenge in hematopoietic differentiation. Here, through gene co-expression network analysis over RNA-seq data generated from representative types of mouse myeloid cells, we obtained a catalog of potential key lncRNAs in the context of mouse myeloid differentiation. Then, employing a widely used in vitro cell model, we screened a novel lncRNA, named Gdal1 (Granulocytic differentiation associated lncRNA 1), from this list and demonstrated that Gdal1 was required for granulocytic differentiation. Furthermore, knockdown of Cebpe, a principal transcription factor of granulocytic differentiation regulation, led to down-regulation of Gdal1, but not vice versa. In addition, expression of genes involved in myeloid differentiation and its regulation, such as Cebpa, were influenced in Gdal1 knockdown cells with differentiation blockage. We thus systematically identified myeloid differentiation associated lncRNAs and substantiated the identification by investigation of one of these lncRNAs on cellular phenotype and gene regulation levels. This study promotes our understanding of the regulation of myeloid differentiation and the characterization of roles of lncRNAs in hematopoietic system. [ABSTRACT FROM AUTHOR]

Steroid metabolism is a fundamental process in the porcine testis to provide testosterone but also estrogens and androstenone, which are essential for the physiology of the boar. This study concerns boars at an early stage of puberty. Using a RT-qPCR approach, we showed that the transcriptional activities of several genes providing key enzymes involved in this metabolism (such as CYP11A1) are correlated. Surprisingly, HSD17B3, a key gene for testosterone production, was absent from this group. An additional weighted gene co-expression network analysis was performed on two large sets of mRNA-seq to identify co-expression modules. Of these modules, two containing either CYP11A1 or HSD17B3 were further analyzed. This comprehensive correlation meta-analysis identified a group of 85 genes with CYP11A1 as hub gene, but did not allow the characterization of a robust correlation network around HSD17B3. As the CYP11A1-group includes most of the genes involved in steroid synthesis pathways (including LHCGR encoding for the LH receptor), it may control the synthesis of most of the testicular steroids. The independent expression of HSD17B3 probably allows part of the production of testosterone to escape this control. This CYP11A1-group contained also INSL3 and AGT genes encoding a peptide hormone and an angiotensin peptide precursor, respectively. [ABSTRACT FROM AUTHOR]

The recently emerged SARS-CoV-2 virus is responsible for the ongoing COVID-19 pandemic that has rapidly developed into a global public health threat. Patients severely affected with COVID-19 present distinct clinical features, including acute respiratory disorder, neutrophilia, cytokine storm, and sepsis. In addition, multiple pro-inflammatory cytokines are found in the plasma of such patients. Transcriptome sequencing of different specimens obtained from patients suffering from severe episodes of COVID-19 shows dynamics in terms of their immune responses. However, those host factors required for SARS-CoV-2 propagation and the underlying molecular mechanisms responsible for dysfunctional immune responses during COVID-19 infection remain elusive. In the present study, we analyzed the mRNA-long non-coding RNA (lncRNA) co-expression network derived from publicly available SARS-CoV-2-infected transcriptome data of human lung epithelial cell lines and bronchoalveolar lavage fluid (BALF) from COVID-19 patients. Through co-expression network analysis, we identified four differentially expressed lncRNAs strongly correlated with genes involved in various immune-related pathways crucial for cytokine signaling. Our findings suggest that the aberrant expression of these four lncRNAs can be associated with cytokine storms and anti-viral responses during severe SARS-CoV-2 infection of the lungs. Thus, the present study uncovers molecular interactions behind the cytokine storm activation potentially responsible for hyper-inflammatory responses in critical COVID-19 patients. [ABSTRACT FROM AUTHOR]

Simple Summary: The amount and fatty acid composition of intramuscular fat are important features for the qualitative characteristics of processed and fresh meat products, but the knowledge of the key molecular drivers controlling these traits is still scant. To this aim, the present study investigated the co-expression networks of genes related to variations in the major fatty acids deposited in pig Semimembranosus muscle. Palmitic and palmitoleic acid contents were associated with a downregulation of genes involved in autophagy, mitochondrial fusion, and mitochondrial activity, suggesting that the deposition of these fatty acids may be enhanced in muscles with a reduced mitochondrial function. A higher proportion of oleic acid and a reduction in the percentages of n-6 and n-3 polyunsaturated fatty acids were related to changes in the mRNA levels of genes involved in Mitogen-Activated Protein Kinase (MAPK) signaling. The obtained results indicated gene expression networks and new candidate genes associated with the studied traits. Further studies are needed to confirm our results and identify in the discussed genes molecular markers for future selection schemes aimed at improving pork nutritional and technological quality. Furthermore, as pigs are considered reliable animal models for several human conditions, the obtained results may also be of interest for improving the knowledge of the molecular pathways associated with obesity and diabetes. To date, high-throughput technology such as RNA-sequencing has been successfully applied in livestock sciences to investigate molecular networks involved in complex traits, such as meat quality. Pork quality depends on several organoleptic, technological, and nutritional characteristics, and it is also influenced by the fatty acid (FA) composition of intramuscular fat (IMF). To explore the molecular networks associated with different IMF FA compositions, the Semimembranosus muscle (SM) from two groups of Italian Large White (ILW) heavy pigs divergent for SM IMF content was investigated using transcriptome analysis. After alignment and normalization, the obtained gene counts were used to perform the Weighted Correlation Network Analysis (WGCNA package in R environment). Palmitic and palmitoleic contents showed association with the same gene modules, comprising genes significantly enriched in autophagy, mitochondrial fusion, and mitochondrial activity. Among the key genes related to these FAs, we found TEAD4, a gene regulating mitochondrial activity that seems to be a promising candidate for further studies. On the other hand, the genes comprised in the modules associated with the IMF contents of oleic, n-6, and n-3 polyunsaturated FAs (PUFAs) were significantly enriched in Mitogen-Activated Protein Kinase (MAPK) signaling, in agreement with previous studies suggesting that several MAPK players may have a primary role in regulating lipid deposition. These results give an insight into the molecular cascade associated with different IMF FA composition in ILW heavy pigs. Further studies are needed to validate the results and confirm whether some of the identified key genes may be effective candidates for pork quality. [ABSTRACT FROM AUTHOR]

The final size of plant leaves is strictly controlled by environmental and genetic factors, which coordinate cell expansion and cell cycle activity in space and time; however, the regulatory mechanisms of leaf growth are still poorly understood. Ginkgo biloba is a dioecious species native to China with medicinally and phylogenetically important characteristics, and its fan-shaped leaves are unique in gymnosperms, while the mechanism of G. biloba leaf development remains unclear. In this study we studied the transcriptome of G. biloba leaves at three developmental stages using high-throughput RNA-seq technology. Approximately 4167 differentially expressed genes (DEGs) were obtained, and a total of 12,137 genes were structure optimized together with 732 new genes identified. More than 50 growth-related factors and gene modules were identified based on DEG and Weighted Gene Co-expression Network Analysis. These results could remarkably expand the existing transcriptome resources of G. biloba, and provide references for subsequent analysis of ginkgo leaf development. [ABSTRACT FROM AUTHOR]

Approximately 10% of agricultural land is subject to periodic flooding, which reduces the growth, survivorship, and yield of most crops, reinforcing the need to understand and enhance flooding resistance in our crops. Here, we generated RNA-Seq data from leaf and root tissue of domesticated sunflower to explore differences in gene expression and alternative splicing (AS) between a resistant and susceptible cultivar under both flooding and control conditions and at three time points. Using a combination of mixed model and gene co-expression analyses, we were able to separate general responses of sunflower to flooding stress from those that contribute to the greater tolerance of the resistant line. Both cultivars responded to flooding stress by upregulating expression levels of known submergence responsive genes, such as alcohol dehydrogenases, and slowing metabolism-related activities. Differential AS reinforced expression differences, with reduced AS frequencies typically observed for genes with upregulated expression. Significant differences were found between the genotypes, including earlier and stronger upregulation of the alcohol fermentation pathway and a more rapid return to pre-flooding gene expression levels in the resistant genotype. Our results show how changes in the timing of gene expression following both the induction of flooding and release from flooding stress contribute to increased flooding tolerance. [ABSTRACT FROM AUTHOR]

The ability to expand crop plantations without irrigation is a major goal to increase agriculture sustainability. To achieve this end, we need to understand the mechanisms that govern plant growth responses under drought conditions. In this study, we combined physiological, transcriptomic, and genomic data to provide a comprehensive picture of drought and recovery responses in the leaves and roots of sugarcane. Transcriptomic profiling using oligoarrays and RNA-seq identified 2898 (out of 21,902) and 46,062 (out of 373,869) transcripts as differentially expressed, respectively. Co-expression analysis revealed modules enriched in photosynthesis, small molecule metabolism, alpha-amino acid metabolism, trehalose biosynthesis, serine family amino acid metabolism, and carbohydrate transport. Together, our findings reveal that carbohydrate metabolism is coordinated with the degradation of amino acids to provide carbon skeletons to the tricarboxylic acid cycle. This coordination may help to maintain energetic balance during drought stress adaptation, facilitating recovery after the stress is alleviated. Our results shed light on candidate regulatory elements and pave the way to biotechnology strategies towards the development of drought-tolerant sugarcane plants. [ABSTRACT FROM AUTHOR]

Drought stress limits plant growth and productivity. It triggers many responses by inducing changes in plant morphology and physiology. KDML105 rice is a key rice variety in Thailand and is normally grown in the northeastern part of the country. The chromosome segment substitution lines (CSSLs) were developed by transferring putative drought tolerance loci (QTLs) on chromosome 1, 3, 4, 8, or 9 into the KDML105 rice genome. CSSL104 is a drought-tolerant line with higher net photosynthesis and leaf water potential than KDML105 rice. The analysis of CSSL104 gene regulation identified the loci associated with these traits via gene co-expression network analysis. Most of the predicted genes are involved in the photosynthesis process. These genes are also conserved in Arabidopsis thaliana. Seven genes encoding chloroplast proteins were selected for further analysis through characterization of Arabidopsis tagged mutants. The response of these mutants to drought stress was analyzed daily for seven days after treatment by scoring green tissue areas via the PlantScreen™ XYZ system. Mutation of these genes affected green areas of the plant and stability index under drought stress, suggesting their involvement in drought tolerance. [ABSTRACT FROM AUTHOR]

Rice (Oryza sativa L.) is more sensitive to drought stress than other cereals. To dissect molecular mechanisms underlying drought-tolerant yield in rice, we applied differential expression and co-expression network approaches to transcriptomes from flag-leaf and emerging panicle tissues of a drought-tolerant yield introgression line, DTY-IL, and the recurrent parent Swarna, under moderate reproductive-stage drought stress. Protein turnover and efficient reactive oxygen species scavenging were found to be the driving factors in both tissues. In the flag-leaf, the responses further included maintenance of photosynthesis and cell wall reorganization, while in the panicle biosynthesis of secondary metabolites was found to play additional roles. Hub genes of importance in differential drought responses included an expansin in the flag-leaf and two peroxidases in the panicle. Overlaying differential expression data with allelic variation in DTY-IL quantitative trait loci allowed for the prioritization of candidate genes. They included a differentially regulated auxin-responsive protein, with DTY-IL-specific amino acid changes in conserved domains, as well as a protein kinase with a DTY-IL-specific frameshift in the C-terminal region. The approach highlights how the integration of differential expression and allelic variation can aid in the discovery of mechanism and putative causal contribution underlying quantitative trait loci for drought-tolerant yield. [ABSTRACT FROM AUTHOR]

Long non-coding RNA (lncRNA) can regulate several aspects of gene expression, being associated with complex phenotypes in humans and livestock species. In taurine beef cattle, recent evidence points to the involvement of lncRNA in feed efficiency (FE), a proxy for increased productivity and sustainability. Here, we hypothesized specific regulatory roles of lncRNA in FE of indicine cattle. Using RNA-Seq data from the liver, muscle, hypothalamus, pituitary gland and adrenal gland from Nellore bulls with divergent FE, we submitted new transcripts to a series of filters to confidently predict lncRNA. Then, we identified lncRNA that were differentially expressed (DE) and/or key regulators of FE. Finally, we explored lncRNA genomic location and interactions with miRNA and mRNA to infer potential function. We were able to identify 126 relevant lncRNA for FE in Bos indicus, some with high homology to previously identified lncRNA in Bos taurus and some possible specific regulators of FE in indicine cattle. Moreover, lncRNA identified here were linked to previously described mechanisms related to FE in hypothalamus-pituitary-adrenal axis and are expected to help elucidate this complex phenotype. This study contributes to expanding the catalogue of lncRNA, particularly in indicine cattle, and identifies candidates for further studies in animal selection and management. [ABSTRACT FROM AUTHOR]

RWP-RK is a plant-specific family of transcription factors, involved in nitrate response, gametogenesis, and nodulation. However, genome-wide characterization, phylogeny, and the regulation of RWP-RK genes in the nodulating and non-nodulating plant species of nitrogen-fixing clade (NFC) are widely unknown. Therefore, we identified a total of 292 RWP-RKs, including 278 RWP-RKs from 25 NFC species and 14 RWP-RKs from the outgroup, Arabidopsis thaliana. We classified the 292 RWP-RKs in two subfamilies: the NIN-like proteins (NLPs) and the RWP-RK domain proteins (RKDs). The transcriptome and phylogenetic analysis of RWP-RKs suggested that, compared to RKD genes, the NLP genes were just upregulated in nitrate response and nodulation. Moreover, nodule-specific NLP genes of some nodulating NFC species may have a common ancestor (OG0002084) with AtNLP genes in A. thaliana. Further, co-expression networks of A.thaliana under N-starvation and N-supplementation conditions revealed that there is a higher correlation between expression of AtNLP genes and symbiotic genes during N-starvation. In P. vulgaris, we confirmed that N-starvation stimulated nodulation by regulating expression of PvNLP2, closely related to AtNLP6 and AtNLP7 with another common origin (OG0004041). Taken together, we concluded that different origins of the NLP genes involved in both N-starvation response and specific expression of nodulation would contribute to the evolution of nodulation in NFC plant species. Our results shed light on the phylogenetic relationships of NLP genes and their differential regulation in nitrate response of A. thaliana and nodulation of NFC. [ABSTRACT FROM AUTHOR]