Genome-wide identification and co-expression network analysis of nuclear factor-Y in barley revealed potential functions in salt stress

Harnessing systems biology approach for characterization of carotenoid biosynthesis pathways in microalgae

Systems biology is an interdisciplinary field that aims to understand complex biological processes at the system level. The data, driven by high-throughput omics technologies, can be used to study the underpinning mechanisms of metabolite production under different conditions to harness this knowledge for the construction of regulatory networks, protein networks, metabolic models, and engineering of strains with enhanced target metabolite production in microalgae. In the current study, we comprehensively reviewed the recent progress in the application of these technologies for the characterization of carotenoid biosynthesis pathways in microalgae. Moreover, harnessing integrated approaches such as network analysis, meta-analysis, and machine learning models for deciphering the complexity of carotenoid biosynthesis pathways were comprehensively discussed.

Global transcriptome analysis identifies critical functional modules associated with multiple abiotic stress responses in microalgae Chromochloris zofingiensis

In the current study, systems biology approach was applied to get a deep insight regarding the regulatory mechanisms of Chromochloris zofingiensis under overall stress conditions. Meta-analysis was performed using p-values combination of differentially expressed genes. To identify the informative models related to stress conditions, two distinct weighted gene co-expression networks were constructed and preservation analyses were performed using medianRankand Zsummary algorithms. Moreover, functional enrichment analysis of non-preserved modules was performed to shed light on the biological performance of underlying genes in the non-preserved modules. In the next step, the gene regulatory networks between top hub genes of non-preserved modules and transcription factors were inferred using ensemble of trees algorithm. Results showed that the power of beta = 7 was the best soft-thresholding value to ensure a scale-free network, leading to the determination of 12 co-expression modules with an average size of 128 genes. Preservation analysis showed that the connectivity pattern of the six modules including the blue, black, yellow, pink, greenyellow, and turquoise changed during stress condition which defined as non-preserved modules. Examples of enriched pathways in non-preserved modules were Oxidative phosphorylation", "Vitamin B6 metabolism", and "Arachidonic acid metabolism". Constructed regulatory network between identified TFs and top hub genes of non-preserved module such as Cz06g10250, Cz03g12130 showed that some specific TFs such as C3H and SQUAMOSA promoter binding protein (SBP) specifically regulates the specific hubs. The current findings add substantially to our understanding of the stress responsive underlying mechanism of C. zofingiensis for future studies and metabolite production programs.

Expression profiling of Nuclear Factor-Y (NF-Y) transcription factors during dehydration and salt stress in finger millet reveals potential candidate genes for multiple stress tolerance

MAIN CONCLUSION

Expression profiling of NF-Y transcription factors during dehydration and salt stress in finger millet genotypes contrastingly differing in tolerance levels identifies candidate genes for further characterization and functional studies. The Nuclear Factor-Y (NF-Y) transcription factors are known for imparting abiotic stress tolerance in different plant species. However, there is no information on the role of this transcription factor family in naturally drought-tolerant crop finger millet (Eleusine coracana L.). Therefore, interpretation of expression profiles against drought and salinity stress may provide valuable insights into specific and/or overlapping expression patterns of Eleusine coracana Nuclear Factor-Y (EcNF-Y) genes. Given this, we identified 59 NF-Y (18 NF-YA, 23 NF-YB, and 18 NF-YC) encoding genes and designated them EcNF-Y genes. Expression profiling of these genes was performed in two finger millet genotypes, PES400 (dehydration and salt stress tolerant) and VR708 (dehydration and salt stress sensitive), subjected to PEG-induced dehydration and salt (NaCl) stresses at different time intervals (0, 6, and 12 h). The qRT-PCR expression analysis reveals that the six EcNF-Y genes namely EcNF-YA1, EcNF-YA5, EcNF-YA16, EcNF-YB6, EcNF-YB10, and EcNF-YC2 might be associated with tolerance to both dehydration and salinity stress in early stress condition (6 h), suggesting the involvement of these genes in multiple stress responses in tolerant genotype. In contrast, the transcript abundance of finger millet EcNF-YA5 genes was also observed in the sensitive genotype VR708 under late stress conditions (12 h) of both dehydration and salinity stress. Therefore, the EcNF-YA5 gene might be important for adaptation to salinity and dehydration stress in sensitive finger millet genotypes. Therefore, this gene could be considered as a susceptibility determinant, which can be edited to impart tolerance. The phylogenetic analyses revealed that finger millet NF-Y genes share strong evolutionary and functional relationship to NF-Ys governing response to abiotic stresses in rice, sorghum, maize, and wheat. This is the first report of expression profiling of EcNF-Ys genes identified from the finger millet genome and reveals potential candidate for enhancing dehydration and salt tolerance.

Trancriptome data mining in combination with co-expression network analysis identifies the functional modules and critical regulators in Hordeum vulgare L. in response to cold stress

Cold stress, as an abiotic stress, is one of the most limiting factors which pose a great threat to the plant's productivity. To understand the transcriptional regulation and connectivity pattern of genes involved in barley cold stress responses, co-expression network analysis was performed based on the global transcriptome profiling. The microarray datasets related to cold stress treatments were retrieved from the Gene Expression Omnibus (GEO) and Array express databases. Four microarray datasets related to cold stress-responsive transcriptome in barley were included in our study. Gene co-expression analysis was constructed using WGCNA method. Module-Trait Relationships (MTR) analysis and hub genes determination and validation were carried out. Finally, transcription factor and kinase regulatory networks were Inferred using machine learning algorithm. The co-expression modules were determined using beta index = 10. In total 13 co-expressed modules were identified with an average size of 153 genes. Functional enrichment based on gene ontology (GO) showed that each of the stress related significant modules were enriched in different biological processes. Annotation of significant modules identifies some TFs and Kinases such as ethylene-responsive transcription factor 1-like, transcription factor PCL1-like, transcription factor MYC2, WRKY, serine/threonine-protein kinase PBL7, and receptor-like protein kinase At2g42960 were contributed in barley cold stress response. Our analysis highlighted the functional importance of ABA signaling pathway, ROS signaling, defensive and protective proteins, degrading protein, Ca2+ related signaling, ribosome-mediated translation and etc. in responding of barley to cold stress condition. The current findings add substantially to our understanding of the cold responsive underlying mechanism of barley which can serve in future studies and breeding programs.

Meta-analysis of transcriptomic profiles in Dunaliella tertiolecta reveals molecular pathway responses to different abiotic stresses

Microalgae are photosynthetic organisms and a potential source of sustainable metabolite production. However, different stress conditions might affect the production of various metabolites. In this study, a meta-analysis of RNA-seq experiments in Dunaliella tertiolecta was evaluated to compare metabolite biosynthesis pathways in response to abiotic stress conditions such as high light, nitrogen deficiency and high salinity. Results showed downregulation of light reaction, photorespiration, tetrapyrrole and lipid-related pathways occurred under salt stress. Nitrogen deficiency mostly induced the microalgal responses of light reaction and photorespiration metabolism. Phosphoenol pyruvate carboxylase, phosphoglucose isomerase, bisphosphoglycerate mutase and glucose-6-phosphate-1-dehydrogenase (involved in central carbon metabolism) were commonly upregulated under salt, light and nitrogen stresses. Interestingly, the results indicated that the meta-genes (modules of genes strongly correlated) were located in a hub of stress-specific protein-protein interaction (PPI) network. Module enrichment of meta-genes PPI networks highlighted the cross-talk between photosynthesis, fatty acids, starch and sucrose metabolism under multiple stress conditions. Moreover, it was observed that the coordinated expression of the tetrapyrrole intermediated with meta-genes was involved in starch biosynthesis. Our results also showed that the pathways of vitamin B6 metabolism, methane metabolism, ribosome biogenesis and folate biosynthesis responded specifically to different stress factors. Since the results of this study revealed the main pathways underlying the abiotic stress, they might be applied in optimised metabolite production by the microalga Dunaliella in future studies. PRISMA check list was also included in the study.

CRISPR-Cas systems and diversity of targeting phages in Lactobacillus johnsonii strains; insights from genome mining approach

Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (CAS) genes make up bacteria's adaptive immune system. These genes protect bacteria from being eaten by bacteriophages. In this study, CRISPR-Cas systems were characterized using a genomic approach. For this purpose, genome sequences of Lactobacillus johnsonii strains were retrieved, and the diversity, occurrence, and evolution of the CRISPR-Cas systems were analyzed. Then, homology analyses of spacer sequences in identified CRISPR arrays were performed to analyze and characterize the diversity of target phages and plasmids. Finally, the evolutionary paths of spaceromes in each subtype of CRISPR arrays were performed using acquisition and deletion events surveyed under the selective pressure of foreign plasmids and phages. Results showed that 138 strains contain valid CRISPR-Cas structures (CRISPR loci together with the Cas genes) in their genomes, which accounted for about 17% of the L. johnsonii studied strains belonging to subtypes II-A, I-E, and I-C. Moreover, results indicated that some specific groups of plasmids were targeted with specific CRISPR array systems. Homology analysis of spacer sequences with phage genomes also revealed that spacers of strains that harbored only CRISPR-Cas subtype-II targeted a greater diversity of foreign phages. In conclusion, the current study indicates that there is great diversity between the CRISPR-Cas systems identified in L. johnsonii strains. Such diverse CRISPR-Cas systems indicate that these systems are naturally active and important in terms of adaptive immunity and evolutionary relationships.

"Millet Models" for harnessing nuclear factor-Y transcription factors to engineer stress tolerance in plants: current knowledge and emerging paradigms

MAIN CONCLUSION

The main purpose of this review is to shed light on the role of millet models in imparting climate resilience and nutritional security and to give a concrete perspective on how NF-Y transcription factors can be harnessed for making cereals more stress tolerant. Agriculture faces significant challenges from climate change, bargaining, population, elevated food prices, and compromises with nutritional value. These factors have globally compelled scientists, breeders, and nutritionists to think of some options that can combat the food security crisis and malnutrition. To address these challenges, mainstreaming the climate-resilient and nutritionally unparalleled alternative crops like millet is a key strategy. The C4 photosynthetic pathway and adaptation to low-input marginal agricultural systems make millets a powerhouse of important gene and transcription factor families imparting tolerance to various kinds of biotic and abiotic stresses. Among these, the nuclear factor-Y (NF-Y) is one of the prominent transcription factor families that regulate diverse genes imparting stress tolerance. The primary purpose of this article is to shed light on the role of millet models in imparting climate resilience and nutritional security and to give a concrete perspective on how NF-Y transcription factors can be harnessed for making cereals more stress tolerant. Future cropping systems could be more resilient to climate change and nutritional quality if these practices were implemented.

Characterization of HSP70 family in watermelon (Citrullus lanatus): identification, structure, evolution, and potential function in response to ABA, cold and drought stress

Watermelon (Citrullus lanatus) as a crop with important economic value, is widely cultivated around the world. The heat shock protein 70 (HSP70) family in plant is indispensable under stress conditions. However, no comprehensive analysis of watermelon HSP70 family is reported to date. In this study, 12 ClHSP70 genes were identified from watermelon, which were unevenly located in 7 out of 11 chromosomes and divided into three subfamilies. ClHSP70 proteins were predicted to be localized primarily in cytoplasm, chloroplast, and endoplasmic reticulum. Two pairs of segmental repeats and 1 pair of tandem repeats existed in ClHSP70 genes, and ClHSP70s underwent strong purification selection. There were many abscisic acid (ABA) and abiotic stress response elements in ClHSP70 promoters. Additionally, the transcriptional levels of ClHSP70s in roots, stems, true leaves, and cotyledons were also analyzed. Some of ClHSP70 genes were also strongly induced by ABA. Furthermore, ClHSP70s also had different degrees of response to drought and cold stress. The above data indicate that ClHSP70s may be participated in growth and development, signal transduction and abiotic stress response, laying a foundation for further analysis of the function of ClHSP70s in biological processes.

Comprehensive Analysis of the NF-YB Gene Family and Expression under Abiotic Stress and Hormone Treatment in Larix kaempferi

NF-YB, a subfamily of Nuclear Factor Y (NF-Y) transcription factor, play crucial role in many biological processes of plant growth and development and abiotic stress responses, and they can therefore be good candidate factors for breeding stress-resistant plants. However, the NF-YB proteins have not yet been explored in Larix kaempferi, a tree species with high economic and ecological values in northeast China and other regions, limiting the breeding of anti-stress L. kaempferi. In order to explore the roles of NF-YB transcription factors in L. kaempferi, we identified 20 LkNF-YB family genes from L. kaempferi full-length transcriptome data and carried out preliminary characterization of them through series of analyses on their phylogenetic relationships, conserved motif structure, subcellular localization prediction, GO annotation, promoter cis-acting elements as well as expression profiles under treatment of phytohormones (ABA, SA, MeJA) and abiotic stresses (salt and drought). The LkNF-YB genes were classified into three clades through phylogenetic analysis and belong to non-LEC1 type NF-YB transcription factors. They have 10 conserved motifs; all genes contain a common motif, and their promoters have various phytohormones and abiotic stress related cis-acting elements. Quantitative real time reverse transcription PCR (RT-qPCR) analysis showed that the sensitivity of the LkNF-YB genes to drought and salt stresses was higher in leaves than roots. The sensitivity of LKNF-YB genes to ABA, MeJA, SA stresses was much lower than that to abiotic stress. Among the LkNF-YBs, LkNF-YB3 showed the strongest responses to drought and ABA treatments. Further protein interaction prediction analysis for LkNF-YB3 revealed that LkNF-YB3 interacts with various factors associated with stress responses and epigenetic regulation as well as NF-YA/NF-YC factors. Taken together, these results unveiled novel L. kaempferi NF-YB family genes and their characteristics, providing the basic knowledge for further in-depth studies on their roles in abiotic stress responses of L. kaempferi.

Genome-Wide Identification and Chilling Stress Analysis of the NF-Y Gene Family in Melon

The nuclear factor Y (NF-Y) transcription factor contains three subfamilies: NF-YA, NF-YB, and NF-YC. The NF-Y family have been reported to be key regulators in plant growth and stress responses. However, little attention has been given to these genes in melon (Cucumis melo L.). In this study, twenty-five NF-Ys were identified in the melon genome, including six CmNF-YAs, eleven CmNF-YBs, and eight CmNF-YCs. Their basic information (gene location, protein characteristics, and subcellular localization), conserved domains and motifs, and phylogeny and gene structure were subsequently analyzed. Results showed highly conserved motifs exist in each subfamily, which are distinct between subfamilies. Most CmNF-Ys were expressed in five tissues and exhibited distinct expression patterns. However, CmNF-YA6, CmNF-YB1/B2/B3/B8, and CmNF-YC6 were not expressed and might be pseudogenes. Twelve CmNF-Ys were induced by cold stress, indicating the NF-Y family plays a key role in melon cold tolerance. Taken together, our findings provide a comprehensive understanding of CmNF-Y genes in the development and stress response of melon and provide genetic resources for solving the practical problems of melon production.

Genome-wide identification, expression analysis, and transcriptome analysis of the IAA gene family in Zoysia japonica

BACKGROUND

AUX/IAA is an essential signaling molecule and has great physiological importance in various plants, but its function in Zoysia japonica remains unknown.

METHODS AND RESULTS

Genome-wide identification and analysis of AUX/IAA genes used bioinformatics methods to investigate the ZjIAA genes' expression of exogenous IAA hydroponics treatment for 2 h by qRT-PCR, control and exogenous IAA treated zoysia were subjected to transcriptome sequencing. ZjIAAs were distributed across the 13 subfamilies by phylogenetic analysis with Oryza sativa and Arabidopsis thaliana. Multiple sequence alignment revealed that the majority of genes were non-canonical ZjIAAs with incomplete domain. The optimal growth concentration of the IAA hormone was 0.05 mM, and the qRT-PCR analysis revealed that eight ZjIAAs were differentially expressed, with seven genes considerably upregulating and one gene significantly downregulating. The result of transcriptome sequencing revealed that 515 differentially expressed genes (DEGs) were identified, with 344 upregulated genes and 171 downregulated genes. A total of 18 genes were annotated as involved in the plant hormone signal transduction pathway. And 8 ZjIAAs exhibited distinct expressions, 7 upregulated, and only one downregulated, according to the qRT-PCR study.

CONCLUSIONS

Genome-wide identification and analysis increased the understanding of the evolution and function of the IAA family in zoysia. DEGs of control and treatment with 0.05 mM exogenous IAA hormone were investigated by transcriptome sequencing. ZjIAAs had substantial variations in the expression of associated genes, with the majority of genes upregulated and 18 genes implicated in plant hormone signal transduction.

Crucial Abiotic Stress Regulatory Network of NF-Y Transcription Factor in Plants

Nuclear Factor-Y (NF-Y), composed of three subunits NF-YA, NF-YB and NF-YC, exists in most of the eukaryotes and is relatively conservative in evolution. As compared to animals and fungi, the number of NF-Y subunits has significantly expanded in higher plants. The NF-Y complex regulates the expression of target genes by directly binding the promoter CCAAT box or by physical interaction and mediating the binding of a transcriptional activator or inhibitor. NF-Y plays an important role at various stages of plant growth and development, especially in response to stress, which attracted many researchers to explore. Herein, we have reviewed the structural characteristics and mechanism of function of NF-Y subunits, summarized the latest research on NF-Y involved in the response to abiotic stresses, including drought, salt, nutrient and temperature, and elaborated the critical role of NF-Y in these different abiotic stresses. Based on the summary above, we have prospected the potential research on NF-Y in response to plant abiotic stresses and discussed the difficulties that may be faced in order to provide a reference for the in-depth analysis of the function of NF-Y transcription factors and an in-depth study of plant responses to abiotic stress.

NF-YB family transcription factors in Arabidopsis: Structure, phylogeny, and expression analysis in biotic and abiotic stresses

Nuclear factor-Y (NF-Y) transcription factors (TFs) are conserved heterotrimeric complexes present and widespread across eukaryotes. Three main subunits make up the structural and functional aspect of the NF-Y TFs: NF-YA, NF-YB and NF-YC, which bind to the conserved CCAAT- box of the promoter region of specific genes, while also interacting with each other, thereby forming myriad combinations. The NF-YBs are expressed differentially in various tissues and plant development stages, likely impacting many of the cellular processes constitutively and under stress conditions. In this study, ten members of NF-YB family from Arabidopsis thaliana were identified and expression profiles were mined from microarray data under different biotic and abiotic conditions, revealing key insights into the involvement of this class of proteins in the cellular and biological processes in Arabidopsis. Analysis of cis-acting regulatory elements (CAREs) indicated the presence of abiotic and biotic stress-related transcription factor binding sites (TFBs), shedding light on the multifaceted roles of these TFs. Microarray data analysis inferred distinct patterns of expression in various tissues under differing treatments such as drought, cold and heat stress as well as bacterial, fungal, and viral stress, indicating their likelihood of having an expansive range of regulatory functions under native and stressed conditions; while quantitative real-time PCR (qRT-PCR) based expression analysis revealed that these TFs get real-time-modulated in a stress dependent manner. This study, overall, provides an understanding of the AtNF-YB family of TFs in their regulation and participation in various morphogenetic and defense- related pathways and can provide insights for development of transgenic plants for trait dependent studies.

The NF-Y Transcription Factor Family in Watermelon: Re-Characterization, Assembly of ClNF-Y Complexes, Hormone- and Pathogen-Inducible Expression and Putative Functions in Disease Resistance

Nuclear factor Y (NF-Y) is a heterotrimeric transcription factor that binds to the CCAAT cis-element in the promoters of target genes and plays critical roles in plant growth, development, and stress responses. In the present study, we aimed to re-characterize the ClNF-Y family in watermelon, examine the assembly of ClNF-Y complexes, and explore their possible involvement in disease resistance. A total of 25 ClNF-Y genes (7 ClNF-YAs, 10 ClNF-YBs, and 8 ClNF-YCs) were identified in the watermelon genome. The ClNF-Y family was comprehensively characterized in terms of gene and protein structures, phylogenetic relationships, and evolution events. Different types of cis-elements responsible for plant growth and development, phytohormones, and/or stress responses were identified in the promoters of the ClNF-Y genes. ClNF-YAs and ClNF-YCs were mainly localized in the nucleus, while most of the ClNF-YBs were localized in the cytoplasm of cells. ClNF-YB5, -YB6, -YB7, -YB8, -YB9, and -YB10 interacted with ClNF-YC2, -YC3, -YC4, -YC5, -YC6, -YC7, and -YC8, while ClNF-YB1 and -YB3 interacted with ClNF-YC1. A total of 37 putative ClNF-Y complexes were identified, e.g., ClNF-YA1, -YA2, -YA3, and -YA7 assembled into 13, 8, 8, and 8 ClNF-Y complexes with different ClNF-YB/-YC heterodimers. Most of the ClNF-Y genes responded with distinct expression patterns to defense hormones such as salicylic acid, methyl jasmonate, abscisic acid, and ethylene precursor 1-aminocyclopropane-1-carboxylate, and to infection by the vascular infecting fungus Fusarium oxysporum f. sp. niveum. Overexpression of ClNF-YB1, -YB8, -YB9, ClNF-YC2, and -YC7 in transgenic Arabidopsis resulted in an earlier flowering phenotype. Overexpression of ClNF-YB8 in Arabidopsis led to enhanced resistance while overexpression of ClNF-YA2 and -YC2 resulted in decreased resistance against Botrytis cinerea. Similarly, overexpression of ClNF-YA3, -YB1, and -YC4 strengthened resistance while overexpression of ClNF-YA2 and -YB8 attenuated resistance against Pseudomonas syringae pv. tomato DC3000. The re-characterization of the ClNF-Y family provides a basis from which to investigate the biological functions of ClNF-Y genes in respect of growth, development, and stress response in watermelon, and the identification of the functions of some ClNF-Y genes in disease resistance enables further exploration of the molecular mechanism of ClNF-Ys in the regulation of watermelon immunity against diverse pathogens.

Identification of potentially functional modules and diagnostic genes related to amyotrophic lateral sclerosis based on the WGCNA and LASSO algorithms

Amyotrophic lateral sclerosis (ALS) is a genetically and phenotypically heterogeneous disease results in the loss of motor neurons. Mounting information points to involvement of other systems including cognitive impairment. However, neither the valid biomarker for diagnosis nor effective therapeutic intervention is available for ALS. The present study is aimed at identifying potentially genetic biomarker that improves the diagnosis and treatment of ALS patients based on the data of the Gene Expression Omnibus. We retrieved datasets and conducted a weighted gene co-expression network analysis (WGCNA) to identify ALS-related co-expression genes. Functional enrichment analysis was performed to determine the features and pathways of the main modules. We then constructed an ALS-related model using the least absolute shrinkage and selection operator (LASSO) regression analysis and verified the model by the receiver operating characteristic (ROC) curve. Besides we screened the non-preserved gene modules in FTD and ALS-mimic disorders to distinct ALS-related genes from disorders with overlapping genes and features. Altogether, 4198 common genes between datasets with the most variation were analyzed and 16 distinct modules were identified through WGCNA. Blue module had the most correlation with ALS and functionally enriched in pathways of neurodegeneration-multiple diseases', 'amyotrophic lateral sclerosis', and 'endocytosis' KEGG terms. Further, some of other modules related to ALS were enriched in 'autophagy' and 'amyotrophic lateral sclerosis'. The 30 top of hub genes were recruited to a LASSO regression model and 5 genes (BCLAF1, GNA13, ARL6IP5, ARGLU1, and YPEL5) were identified as potentially diagnostic ALS biomarkers with validating of the ROC curve and AUC value.

Transcriptome Analysis Revealed the Molecular Response Mechanism of Non-heading Chinese Cabbage to Iron Deficiency Stress

Iron is a trace metal that is found in animals, plants, and the human body. Human iron absorption is hampered by plant iron shortage, which leads to anemia. Leafy vegetables are one of the most direct and efficient sources of iron for humans. Despite the fact that ferrotrophic disorder is common in calcareous soil, however, non-heading Chinese cabbage performs a series of reactions in response to iron deficiency stress that help to preserve iron homeostasis in vivo. In this study, we discovered that iron deficiency stress caused leaf yellowing and impeded plant development in both iron-deficient and control treatments by viewing or measuring phenotypic, chlorophyll content, and Fe²⁺ content in both iron-deficient and control treatments. We found a total of 9213 differentially expressed genes (DEGs) in non-heading Chinese cabbage by comparing root and leaf transcriptome data with iron deficiency and control treatments. For instance, 1927 DEGs co-expressed in root and leaf, including 897 up-regulated and 1030 down-regulated genes, respectively. We selected some key antioxidant genes, hormone signal transduction, iron absorption and transport, chlorophyll metabolism, and transcription factors involved in the regulation of iron deficiency stress utilizing GO enrichment, KEGG enrichment, multiple types of functional annotation, and Weighted Gene Co-expression Network Analysis (WGCNA). This study identifies prospective genes for maintaining iron homeostasis under iron-deficient stress, offering a theoretical foundation for further research into the molecular mechanisms of greater adaptation to iron-deficient stress, and perhaps guiding the development of iron-tolerant varieties.

Identification and Comprehensive Analysis of the Nuclear Factor-Y Family Genes Reveal Their Multiple Roles in Response to Nutrient Deficiencies in Brassica napus

Nuclear Factor-Y (NF-Y) transcription factors play vital roles in plant abiotic stress response. Here, the NF-Y family in Brassica napus, which is hyper-sensitive to nitrogen (N) deprivation, was comprehensively identified and systematically characterized. A total of 108 NF-Y family members were identified in B. napus and categorized into three subfamilies (38 NF-YA, 46 NF-YB and 24 NF-YC; part of the Arabidopsis NF-YC homologous genes had been lost during B. napus evolution). In addition, the expansion of the NF-Y family in B. napus was driven by whole-genome duplication and segmental duplication. Differed expression patterns of BnaNF-Ys were observed in response to multiple nutrient starvations. Thirty-four genes were regulated only in one nutrient deficient condition. Moreover, more BnaNF-YA genes were differentially expressed under nutrient limited environments compared to the BnaNF-YB and BnaNF-YC subfamilies. Sixteen hub genes responded diversely to N deprivation in five rapeseed tissues. In summary, our results laid a theoretical foundation for the follow-up functional study of the key NF-Y genes in B. napus in regulating nutrient homeostasis, especially N.

Expression of the Malus sieversii NF-YB21 Encoded Gene Confers Tolerance to Osmotic Stresses in Arabidopsis thaliana

Drought is the main environmental factor that limits the yield and quality of apples (Malus × domestica) grown in arid and semi-arid regions. Nuclear factor Ys (NF-Ys) are important transcription factors involved in the regulation of plant growth, development, and various stress responses. However, the function of NF-Y genes is poorly understood in apples. Here, we identified 43 NF-Y genes in the genome of apples and conducted an initial functional characterization of the apple NF-Y. Expression analysis of NF-Y members in M. sieversii revealed that a large number of NF-Ys were highly expressed in the roots compared with the leaves, and a large proportion of NF-Y genes responded to drought treatment. Furthermore, heterologous expression of MsNF-YB21, which was significantly upregulated by drought, led to a longer root length and, thus, conferred improved osmotic and salt tolerance in Arabidopsis. Moreover, the physiological analysis of MsNF-YB21 overexpression revealed enhanced antioxidant systems, including antioxidant enzymes and compatible solutes. In addition, genes encoding catalase (AtCAT2, AtCAT3), superoxide dismutase (AtFSD1, AtFSD3, AtCSD1), and peroxidase (AtPER12, AtPER42, AtPER47, AtPER51) showed upregulated expression in the MsNF-YB21 overexpression lines. These results for the MsNF-Y gene family provide useful information for future studies on NF-Ys in apples, and the functional analysis of MsNF-YB21 supports it as a potential target in the improvement of apple drought tolerance via biotechnological strategies.

Weighted gene co-expression network analysis identifies modules and functionally enriched pathways in the lactation process

The exponential growth in knowledge has resulted in a better understanding of the lactation process in a wide variety of animals. However, the underlying genetic mechanisms are not yet clearly known. In order to identify the mechanisms involved in the lactation process, various mehods, including meta-analysis, weighted gene co-express network analysis (WGCNA), hub genes identification, gene ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment at before peak (BP), peak (P), and after peak (AP) stages of the lactation processes have been employed. A total of 104, 85, and 26 differentially expressed genes were identified based on PB vs. P, BP vs. AP, and P vs. AP comparisons, respectively. GO and KEGG pathway enrichment analysis revealed that DEGs were significantly enriched in the "ubiquitin-dependent ERAD" and the "chaperone cofactor-dependent protein refolding" in BP vs. P and P vs. P, respectively. WGCNA identified five significant functional modules related to the lactation process. Moreover, GJA1, AP2A2, and NPAS3 were defined as hub genes in the identified modules, highlighting the importance of their regulatory impacts on the lactation process. The findings of this study provide new insights into the complex regulatory networks of the lactation process at three distinct stages, while suggesting several candidate genes that may be useful for future animal breeding programs. Furthermore, this study supports the notion that in combination with a meta-analysis, the WGCNA represents an opportunity to achieve a higher resolution analysis that can better predict the most important functional genes that might provide a more robust bio-signature for phenotypic traits, thus providing more suitable biomarker candidates for future studies.

Weighted gene co-expression network analysis of the salt-responsive transcriptomes reveals novel hub genes in green halophytic microalgae Dunaliella salina

Despite responses to salinity stress in Dunaliella salina, a unicellular halotolerant green alga, being subject to extensive study, but the underlying molecular mechanism remains unknown. Here, Empirical Bayes method was applied to identify the common differentially expressed genes (DEGs) between hypersaline and normal conditions. Then, using weighted gene co-expression network analysis (WGCNA), which takes advantage of a graph theoretical approach, highly correlated genes were clustered as a module. Subsequently, connectivity patterns of the identified modules in two conditions were surveyed to define preserved and non-preserved modules by combining the Zsummary and medianRank measures. Finally, common and specific hub genes in non-preserved modules were determined using Eigengene-based module connectivity or module membership (k_ME) measures and validation was performed by using leave-one-out cross-validation (LOOCV). In this study, the power of beta = 12 (scale-free R2 = 0.8) was selected as the soft-thresholding to ensure a scale-free network, which led to the identification of 15 co-expression modules. Results also indicate that green, blue, brown, and yellow modules are non-preserved in salinity stress conditions. Examples of enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in non-preserved modules are Sulfur metabolism, Oxidative phosphorylation, Porphyrin and chlorophyll metabolism, Vitamin B6 metabolism. Moreover, the systems biology approach was applied here, proposed some salinity specific hub genes, such as radical-induced cell death1 protein (RCD1), mitogen-activated protein kinase kinase kinase 13 (MAP3K13), long-chain acyl-CoA synthetase (ACSL), acetyl-CoA carboxylase, biotin carboxylase subunit (AccC), and fructose-bisphosphate aldolase (ALDO), for the development of metabolites accumulating strains in D. salina.

miR169 and PmRGL2 synergistically regulate the NF-Y complex to activate dormancy release in Japanese apricot (Prunus mume Sieb. et Zucc.)

This study is the first to demonstrate that GA₄-induced dormancy release is associated with the NF-Y complex, which interacts with gibberellin inhibitor RGL2 in Japanese apricot. Seasonal dormancy is not only vital for the survival in cold winter but also affects flowering of temperate fruit trees and the dormancy release depends on the accumulation of the cold temperatures (Chilling requirement-CR). To understand the mechanism of dormancy release in deciduous fruit crops, we compared miRNA sequencing data during the transition stage from paradormancy to dormancy release in the Japanese apricot and found that the miR169 family showed significant differentially up-regulated expression during dormancy induction and was down-regulated during the dormancy release periods. The 5' RACE assay and RT-qPCR validated its target gene NUCLEAR FACTOR-Y subunit A (NF-YA), which exhibited the opposite expression pattern. Further study showed that exogenous GA₄ could inhibit the expression of the gibberellic acid (GA) signal transduction suppressor PmRGL2 (RGA-LIKE 2) and promote the expression of NF-Y. Moreover, the interaction between the NF-Y family and GA inhibitor PmRGL2 was verified by the yeast-two-hybrid (Y2H) system and a bimolecular fluorescence complementarity (BiFC) experiment. These results suggest that synergistic regulation of the NF-Y and PmRGL2 complex leads to the activation of dormancy release induced by GA₄. These findings will help to elucidate the functional and regulatory roles of miR169 and NF-Y complex in seasonal bud dormancy induced by GA in Japanese apricot and provide new insights for the discovery of dormancy release mechanisms in woody plants.

Programmed cell death (PCD) control in plants: New insights from the Arabidopsis thaliana deathosome

Programmed cell death (PCD) is a genetically controlled process that leads to cell suicide in both eukaryotic and prokaryotic organisms. In plants PCD occurs during development, defence response and when exposed to adverse conditions. PCD acts controlling the number of cells by eliminating damaged, old, or unnecessary cells to maintain cellular homeostasis. Unlike in animals, the knowledge about PCD in plants is limited. The molecular network that controls plant PCD is poorly understood. Here we present a review of the current mechanisms involved with the genetic control of PCD in plants. We also present an updated version of the AtLSD1 deathosome, which was previously proposed as a network controlling HR-mediated cell death in Arabidopsis thaliana. Finally, we discuss the unclear points and open questions related to the AtLSD1 deathosome.

Transcriptome signature of two lactation stages in Ghezel sheep identifies using RNA-Sequencing

The expression of genes and their regulation during lactation in Ghezel sheep breed remains less understood. To explore the underlying molecular mechanism of the lactation process in the mammary gland, transcriptome profiles of Iranian fat-tailed Ghezel sheep breed milk at two stages, before (BF) and after peak (AF) stages of lactation were investigated. Functional impacts of differentially expressed genes (DEGs) between BF and AF stages were surveyed using Gene Ontology (GO) and Protein-Protein Interaction (PPI) network analysis. Totally, 75 DEGs were identified between BF and AF stages of lactation. The RNA-Seq results were validated by Q-RT-PCR. Gene ontology of DEGs mainly enriched in metabolic process and oxidative phosphorylation. PPI network analysis also highlighted the contribution of peroxisome proliferator-activated receptors (PPAR) signaling, oxidative phosphorylation and metabolic pathways in the lactation process. Intriguingly, the genes involved in fat metabolism dominantly down-regulated at AF stage. Our results provide new insight into transcriptional changes and add to growing body of knowledge on the lactation process in fat-tailed sheep breeds.

Systems biology approach identifies functional modules and regulatory hubs related to secondary metabolites accumulation after transition from autotrophic to heterotrophic growth condition in microalgae

Heterotrophic growth mode is among the most promising strategies put forth to overcome the low biomass and secondary metabolites productivity challenge. To shedding light on the underlying molecular mechanisms, transcriptome meta-analysis was integrated with weighted gene co-expression network analysis (WGCNA), connectivity analysis, functional enrichment, and hubs identification. Meta-analysis and Functional enrichment analysis demonstrated that most of the biological processes are up-regulated at heterotrophic growth condition, which leads to change of genetic architectures and phenotypic outcomes. WGNCA analysis of meta-genes also resulted four significant functional modules across logarithmic (LG), transition (TR), and production peak (PR) phases. The expression pattern and connectivity characteristics of the brown module as a non-preserved module vary across LG, TR, and PR phases. Functional analysis identified Carotenoid biosynthesis, Fatty acid metabolism and Methane metabolism as enriched pathways in the non-preserved module. Our integrated approach was applied here, identified some hubs, such as a serine hydroxymethyltransferase (SHMT1), which is the best candidate for development of metabolites accumulating strains in microalgae. Current study provided a new insight into underlying metabolite accumulation mechanisms and opens new avenue for the future applied studies in the microalgae field.

Integration of Cross Species RNA-seq Meta-Analysis and Machine-Learning Models Identifies the Most Important Salt Stress-Responsive Pathways in Microalga Dunaliella

Photosynthetic microalgae are potentially yielding sources of different high-value secondary metabolites. Salinity is a complex stress that influences various metabolite-related pathways in microalgae. To obtain a clear view of the underlying metabolic pathways and resolve contradictory information concerning the transcriptional regulation of Dunaliella species in salt stress conditions, RNA-seq meta-analysis along with systems levels analysis was conducted. A p-value combination technique with Fisher method was used for cross species meta-analysis on the transcriptomes of two Dunaliella salina and Dunaliellatertiolecta species. The potential functional impacts of core meta-genes were surveyed based on gene ontology and network analysis. In the current study, the integration of supervised machine-learning algorithms with RNA-seq meta-analysis was performed. The analysis shows that the lipid and nitrogen metabolism, structural proteins of photosynthesis apparatus, chaperone-mediated autophagy, and ROS-related genes are the keys and core elements of the Dunaliella salt stress response system. Cross-talk between Ca²⁺ signal transduction, lipid accumulation, and ROS signaling network in salt stress conditions are also proposed. Our novel approach opens new avenues for better understanding of microalgae stress response mechanisms and for selection of candidate gene targets for metabolite production in microalgae.