Molecular Insight of Plants Response to Drought Stress: Perspectives and New Insights towards Food Security.

In the face of climate-induced challenges, understanding the intricate molecular mechanisms underlying drought tolerance in plants has become imperative [...].

Variation in WIDTH OF LEAF AND GRAIN contributes to grain and leaf size by controlling LARGE2 stability in rice.

Grain and flag leaf size are two important agronomic traits that influence grain yield in rice (Oryza sativa). Many QTLs and genes that regulate these traits individually have been identified, however, few QTLs and genes that simultaneously control these two traits have been identified. In this study, we conducted a genome-wide association analysis in rice and detected a major locus, WIDTH OF LEAF AND GRAIN (WLG), that associated with both grain and flag leaf width. WLG encodes a RING-domain E3 ubiquitin ligase. WLGhap.B, which possesses five SNP variations compared to WLGhap.A, encodes a protein with enhanced ubiquitination activity that confers increased rice leaf width and grain size, whereas mutation of WLG leads to narrower leaves and smaller grains. Both WLGhap.A and WLGhap.B interact with LARGE2, a HETC-type E3 ligase, however, WLGhap.B exhibits stronger interaction with LARGE2, thus higher ubiquitination activity towards LARGE2 compared with WLGhap.A. Lysine1021 is crucial for the ubiquitination of LARGE2 by WLG. Loss-of-function of LARGE2 in wlg-1 phenocopies large2-c in grain and leaf width, suggesting that WLG acts upstream of LARGE2. These findings reveal the genetic and molecular mechanism by which the WLG-LARGE2 module mediates grain and leaf size in rice, and suggest the potential of WLGhap.B in improving rice yield.

Analysis of blood biochemistry and non-targeted metabolomics of endometritis in dairy cows.

The incidence of bovine endometritis, which has a negative impact on the reproduction of dairy cows, has been recently increasing. In this study, the differential markers and metabolites of healthy cows and cows with endometritis were analyzed by measuring blood biochemical indicators and immune factors using biochemical and enzyme-linked immunosorbent assay kits combined with nontargeted metabolomics. The LC-QTOF platform was used to evaluate the serum metabolomics of healthy cows and cows with endometritis after 21-27 days of calving. The results showed that glucose, free fatty acid, calcium, sodium, albumin, and alanine aminotransferase levels were significantly lower in the serum of cows with endometritis than in healthy cows (P < 0.05). However, the serum potassium, interleukin-1, interleukin-6, and tumor necrosis factor levels were significantly higher in cows with endometritis (P < 0.05). In addition, the serum metabolome data analysis of the two groups showed that the expression of 468 metabolites was significantly different (P < 0.05), of which 291 were upregulated and 177 were downregulated. These metabolites were involved in 78 metabolic pathways, including amino acid, nucleotide, carbohydrate, lipid, and vitamin metabolism pathways; signal transduction pathways, and other biological pathways. Taken together, negative energy balance and immune activation, which are related to local abnormalities in amino acid, lipid, and carbohydrate metabolism, were the important causes of endometritis in dairy cows. Metabolites such as glucose, carnosine, dehydroascorbic acid, L-malic acid, tetrahydrofolic acid, and UDP-glucose may be used as key indicators in the hematological diagnosis and treatment of endometritis in dairy cows.

Estimation of genetic parameters of Holstein calf survival.

Calf survival is not only an animal welfare issue but also helps to avoid huge losses in economic and genetic material due to calf mortality. Therefore, improving calf survival is essential in dairy breeding. The objective of this study was to explore the factors affecting the survival of Holstein calves in the Ningxia Region and to estimate the genetic parameters of calves using linear models and threshold models. Descriptive statistics were made for 43,847 Holstein calves born from 2018 to 2022 in Ningxia. The number of calves that died at 2-30 d was the highest, the survival rate was the lowest at 451-750 d, followed by 61-180 d and 2-30 d. Studies on the survival rates of calves born in different months have found that calves born in April have the lowest survival rates and calves born in October and December have higher survival rates. Calves born in autumn, third parity, and singleton calves are more likely to survive. The heritability of calf survival traits ranged from 0.002 ~ 0.136. Thus survival is a low heritability trait. Genetic correlation between different survival stages ranged from 0.3991 (2-30 d to 451-750 d) to 0.9985 (361-450 d to 451-750 d), the phenotypic correlation ranged from 0.1476 (2-30 d to 451-750 d) to 0.9582 (361-450 d to 451-750 d). The low genetic correlation between early and late survival suggests that survival in early and late stages may be influenced by different genetic factors. This study is helpful to understand the survival status of Holstein calves and provide a theoretical basis for improving the survival rate of calves.

Polymorphism of ADAM12, DPP6 and PRKN genes and their associations with milk production traits in Holstein.

Milk production traits as the most important economic traits of dairy cows, they directly reflect the benefits of breeding and the economic benefits of pasture. In this study, A disintegrin and metalloproteinase-12 (ADAM12), Parkinson's disease gene 2 (PRKN) and dipeptidyl peptidase-like protein subtype 6 (DPP6) polymorphism in 384 Chinese Holstein cows were detected by time-of-flight mass spectrometry and through statistical analysis using software such as Popgene 32, SAS 9.4 and Origin 2022, the relationship between single nucleotide polymorphisms (SNPs) of three genes with four milk production traits such as daily milk yield (DMY), milk fat percentage (MFP), milk protein percentage (MPP) and somatic cell score (SCS) was verified at molecular level. The results showed that four polymorphic loci (116,467,133, 116,604,487, 116,618,268 and 116,835,111) of DPP6 gene, two polymorphic loci (97,665,052 and 97,159,837) of PRKN gene and two polymorphic loci (45,542,714 and 45,553,888) of ADAM12 gene were detected. PRKN-97665052, DPP6-116467133, ADAM12-45553888, DPP6-116604487 and DPP6-116835111 were all in Hardy-Weinberg equilibrium state (p > .05). ADAM12-45542714, PRKN-97159837 and PRKN-97665052 were moderately polymorphic (0.25 ≤ PIC <0.50) in Holstein. It is evident that the selection potential and genetic variation of these five loci are relatively large, and the genetic richness is relatively high. The correlation analysis of different genotypes between these eight loci and milk production traits of Holstein showed that ADAM12-45542714 and DPP6-116835111 (p < .01) had an extremely significant effects on the DMY of Chinese Holstein in Ningxia, while PRKN-97665052 had an extremely significant effect on MFP (p < .01). The effect of PRKN-97665052 and DPP6-116467133 on MPP of Holstein were extremely significant (p < .01). DPP6-116618268 had an extremely significant effect on the SCS of Holstein in Ningxia (p < .01), and AA genotype individuals showed a higher SCS than GG genotype individuals; the other two loci (ADAM12-45553888 and DPP6-116604487) had no significant effects on milk production traits of Holstein (p > .05). In addition, through the joint analysis of DPP6, PRKN and ADAM12 gene loci, it was found that the interaction effect between the three gene loci could significantly affect the DMY, SCS (p < .01) and MPP (p < .05). In conclusion, several different loci of DPP6, PRKN and ADAM12 genes can affect the milk production traits of Holstein to different degrees. PRKN, DPP6 and ADAM12 genes can be used as potential candidate genes for milk production traits of Holstein for marker-assisted selection, providing theoretical basis for breeding of Holstein.

Control of rice ratooning ability by a nucleoredoxin that inhibits histidine kinase dimerization to attenuate cytokinin signaling in axillary buds.

Rice ratooning, the fast outgrowth of dormant buds on stubble, is an important cropping practice in rice production. However, the low ratooning ability (RA) of most rice varieties restricts the application of this cost-efficient system, and the genetic basis of RA remains unknown. In this study, we dissected the genetic architecture of RA by a genome-wide association study in a natural rice population. Rice ratooning ability 3 (RRA3), encoding a hitherto not characterized nucleoredoxin involved in reduction of disulfide bonds, was identified as the causal gene of a major locus controlling RA. Overexpression of RRA3 in rice significantly accelerated leaf senescence and reduced RA, whereas knockout of RRA3 significantly delayed leaf senescence and increased RA and ratoon yield. We demonstrated that RRA3 interacts with Oryza sativa histidine kinase 4 (OHK4), a cytokinin receptor, and inhibits the dimerization of OHK4 through disulfide bond reduction. This inhibition ultimately led to decreased cytokinin signaling and reduced RA. In addition, variations in the RRA3 promoter were identified to be associated with RA. Introgression of a superior haplotype with weak expression of RRA3 into the elite rice variety Guichao 2 significantly increased RA and ratoon yield by 23.8%. Collectively, this study not only uncovers an undocumented regulatory mechanism of cytokinin signaling through de-dimerization of a histidine kinase receptor-but also provides an eximious gene with promising value for ratoon rice breeding.

ABF1 Positively Regulates Rice Chilling Tolerance via Inducing Trehalose Biosynthesis.

Chilling stress seriously limits grain yield and quality worldwide. However, the genes and the underlying mechanisms that respond to chilling stress remain elusive. This study identified ABF1, a cold-induced transcription factor of the bZIP family. Disruption of ABF1 impaired chilling tolerance with increased ion leakage and reduced proline contents, while ABF1 over-expression lines exhibited the opposite tendency, suggesting that ABF1 positively regulated chilling tolerance in rice. Moreover, SnRK2 protein kinase SAPK10 could phosphorylate ABF1, and strengthen the DNA-binding ability of ABF1 to the G-box cis-element of the promoter of TPS2, a positive regulator of trehalose biosynthesis, consequently elevating the TPS2 transcription and the endogenous trehalose contents. Meanwhile, applying exogenous trehalose enhanced the chilling tolerance of abf1 mutant lines. In summary, this study provides a novel pathway 'SAPK10-ABF1-TPS2' involved in rice chilling tolerance through regulating trehalose homeostasis.

Regulation of FLC nuclear import by coordinated action of the NUP62-subcomplex and importin ß SAD2.

Flowering locus C (FLC) is a central transcriptional repressor that controls flowering time. However, how FLC is imported into the nucleus is unknown. Here, we report that Arabidopsis nucleoporins 62 (NUP62), NUP58, and NUP54 composed NUP62-subcomplex modulates FLC nuclear import during floral transition in an importin α-independent manner, via direct interaction. NUP62 recruits FLC to the cytoplasmic filaments and imports it into the nucleus through the NUP62-subcomplex composed central channel. Importin ß supersensitive to ABA and drought 2 (SAD2), a carrier protein, is critical for FLC nuclear import and flower transition, which facilitates FLC import into the nucleus mainly through the NUP62-subcomplex. Proteomics, RNA-seq, and cell biological analyses indicate that the NUP62-subcomplex mainly mediates the nuclear import of cargos with unconventional nuclear localization sequences (NLSs), such as FLC. Our findings illustrate the mechanisms of the NUP62-subcomplex and SAD2 on FLC nuclear import process and floral transition, and provide insights into the role of NUP62-subcomplex and SAD2 in protein nucleocytoplasmic transport in plants.

Cell type-specific proteomics uncovers a RAF15-SnRK2.6/OST1 kinase cascade in guard cells.

Multicellular organisms such as plants contain various cell types with specialized functions. Analyzing the characteristics of each cell type reveals specific cell functions and enhances our understanding of organization and function at the organismal level. Guard cells (GCs) are specialized epidermal cells that regulate the movement of the stomata and gaseous exchange, and provide a model genetic system for analyzing cell fate, signaling, and function. Several proteomics analyses of GC are available, but these are limited in depth. Here we used enzymatic isolation and flow cytometry to enrich GC and mesophyll cell protoplasts and perform in-depth proteomics in these two major cell types in Arabidopsis leaves. We identified approximately 3,000 proteins not previously found in the GC proteome and more than 600 proteins that may be specific to GC. The depth of our proteomics enabled us to uncover a guard cell-specific kinase cascade whereby Raf15 and Snf1-related kinase2.6 (SnRK2.6)/OST1(open stomata 1) mediate abscisic acid (ABA)-induced stomatal closure. RAF15 directly phosphorylated SnRK2.6/OST1 at the conserved Ser175 residue in its activation loop and was sufficient to reactivate the inactive form of SnRK2.6/OST1. ABA-triggered SnRK2.6/OST1 activation and stomatal closure was impaired in raf15 mutants. We also showed enrichment of enzymes and flavone metabolism in GC, and consistent, dramatic accumulation of flavone metabolites. Our study answers the long-standing question of how ABA activates SnRK2.6/OST1 in GCs and represents a resource potentially providing further insights into the molecular basis of GC and mesophyll cell development, metabolism, structure, and function.

Bna.EPF2 Enhances Drought Tolerance by Regulating Stomatal Development and Stomatal Size in Brassica napus.

Drought stress severely affects global plant growth and production. The enhancement of plant water-use efficiency (WUE) and drought tolerance by the manipulation of the stomata is an effective strategy to deal with water shortage. However, increasing the WUE and drought tolerance by manipulation on the stomata has rarely been tested in Brassica napus. Here, we isolated Bna.EPF2, an epidermal patterning factor (EPF) in Brassica napus (ecotype Westar), and identified its role in drought performance. Bna.EPF2 overexpression lines had a reduction average of 19.02% in abaxial stomatal density and smaller stomatal pore size, leading to approximately 25% lower transpiration, which finally resulted in greater instantaneous WUE and enhanced drought tolerance. Interestingly, the reduction in stomatal density did not affect the CO2 assimilation or yield-related agronomic traits in Bna.EPF2 overexpression plants. Together with the complementation of Bna.EPF2 significantly decreasing the stomatal density of Arabidopsis epf2, and Bna.EPF2 being expressed in mature guard cells, these results suggest that Bna.EPF2 not only functions in stomatal density development, but also in stomatal dimension in Brassicas. Taken together, our results suggest that Bna.EPF2 improves WUE and drought tolerance by the regulation of stomatal density and stomatal size in Brassica without growth and yield penalty, and provide insight into the manipulation of this gene in the breeding of drought tolerant plants with increased production under water deficit conditions.

PECTIN METHYLESTERASE INHIBITOR18 functions in stomatal dynamics and stomatal dimension.

Pectin methylesterification in guard cell (GC) walls plays an important role in stomatal development and stomatal response to external stimuli, and pectin methylesterase inhibitors (PMEIs) modulate pectin methylesterification by inhibition of pectin methylesterase (PME). However, the function of PMEIs has not been reported in stomata. Here, we report the role of Arabidopsis (Arabidopsis thaliana) PECTIN METHYLESTERASE INHIBITOR18 in stomatal dynamic responses to environmental changes. PMEI18 mutation increased pectin demethylesterification and reduced pectin degradation, resulting in increased stomatal pore size, impaired stomatal dynamics, and hypersensitivity to drought stresses. In contrast, overexpression of PMEI18 reduced pectin demethylesterification and increased pectin degradation, causing more rapid stomatal dynamics. PMEI18 interacted with PME31 in plants, and in vitro enzymatic assays demonstrated that PMEI18 directly inhibits the PME activity of PME31 on pectins. Genetic interaction analyses suggested that PMEI18 modulates stomatal dynamics mainly through inhibition of PME31 on pectin methylesterification in cell walls. Our results provide insight into the molecular mechanism of the PMEI18-PME31 module in stomatal dynamics and highlight the role of PMEI18 and PME31 in stomatal dynamics through modulation of pectin methylesterification and distribution in GC walls.

IQ67 DOMAIN protein 21 is critical for indentation formation in pavement cell morphogenesis.

In plants, cortical microtubules anchor to the plasma membrane in arrays and play important roles in cell shape. However, the molecular mechanism of microtubule binding proteins, which connect the plasma membrane and cortical microtubules in cell morphology remains largely unknown. Here, we report that a plasma membrane and microtubule dual-localized IQ67 domain protein, IQD21, is critical for cotyledon pavement cell (PC) morphogenesis in Arabidopsis. iqd21 mutation caused increased indentation width, decreased lobe length, and similar lobe number of PCs, whereas IQD21 overexpression had a different effect on cotyledon PC shape. Weak overexpression led to increased lobe number, decreased indentation width, and similar lobe length, while moderate or great overexpression resulted in decreased lobe number, indentation width, and lobe length of PCs. Live-cell observations revealed that IQD21 accumulation at indentation regions correlates with lobe initiation and outgrowth during PC development. Cell biological and genetic approaches revealed that IQD21 promotes transfacial microtubules anchoring to the plasma membrane via its polybasic sites and bundling at the indentation regions in both periclinal and anticlinal walls. IQD21 controls cortical microtubule organization mainly through promoting Katanin 1-mediated microtubule severing during PC interdigitation. These findings provide the genetic evidence that transfacial microtubule arrays play a determinant role in lobe formation, and the insight into the molecular mechanism of IQD21 in transfacial microtubule organization at indentations and puzzle-shaped PC development.

Multiple cyclic nucleotide-gated channels function as ABA-activated Ca2+ channels required for ABA-induced stomatal closure in Arabidopsis.

Abscisic acid (ABA)-activated inward Ca2+-permeable channels in the plasma membrane (PM) of guard cells are required for the initiation and regulation of ABA-specific cytosolic Ca2+ signaling and stomatal closure in plants. But the identities of the PM Ca2+ channels are still unknown. We hypothesized that the ABA-activated Ca2+ channels consist of multiple CYCLIC NUCLEOTIDE-GATED CHANNEL (CNGC) proteins from the CNGC family, which is known as a Ca2+-permeable channel family in Arabidopsis (Arabidopsis thaliana). In this research, we observed high expression of multiple CNGC genes in Arabidopsis guard cells, namely CNGC5, CNGC6, CNGC9, and CNGC12. The T-DNA insertional loss-of-function quadruple mutant cngc5-1 cngc6-2 cngc9-1 cngc12-1 (hereafter c5/6/9/12) showed a strong ABA-insensitive phenotype of stomatal closure. Further analysis revealed that ABA-activated Ca2+ channel currents were impaired, and ABA-specific cytosolic Ca2+ oscillation patterns were disrupted in c5/6/9/12 guard cells compared with in wild-type guard cells. All ABA-related phenotypes of the c5/6/9/12 mutant were successfully rescued by the expression of a single gene out of the four CNGCs under the respective native promoter. Thus, our findings reveal a type of ABA-activated PM Ca2+ channel comprising multiple CNGCs, which is essential for ABA-specific Ca2+ signaling of guard cells and ABA-induced stomatal closure in Arabidopsis.

Identification of core genes affecting IMF deposition in bovine.

Intramuscular fat (IMF) content is an important economic factor in beef production. However, knowledge on the key factors controlling bovine IMF is limited. In this study, using weighted gene co-expression network analysis (WGCNA), nine modules were identified and the number of transcripts in these modules ranged from 36 to 3191. Two modules were found to be significantly associated with fat deposition and three genes (TCAP, MYH7, and TNNC1) were further identified by Protein-protein interaction (PPI), which may be the hub genes regulating bovine IMF deposition. In addition, considering the genetic variation, the PCK1 gene was found by functional enrichment analysis of overlapping genes, which was previously reported to be involved in IMF deposition. We noted that the core promoter region of buffalo PCK1 binds to transcription factors involved in lipid metabolism while cattle PCK1 binds transcription factors involved in muscle development. The results suggest that PCK1 participated in IMF deposition of buffalo and cattle in different ways. In summary, gene expression networks and new candidate genes associated with IMF deposition identified in this study. This would lay the foundation for further research into the molecular regulatory mechanisms underlying bovine IMF deposition.

Identification of critical lncRNAs for milk fat metabolism in dairy cows using WGCNA and the construction of a ceRNAs network.

As key regulators, long non-coding RNAs (lncRNAs) play a crucial role in the ruminant mammary gland. However, the function of lncRNAs in milk fat synthesis from dairy cows is largely unknown. In this study, we used the weighted gene co-expression network analysis (WGCNA) to comprehensive analyze the expression profile data of lncRNAs from the group's previous Illumina PE150 sequencing results based on bovine mammary epithelial cells from high- and low-milk-fat-percentage (MFP) cows, and identify core_lncRNAs significantly associated with MFP by module membership (MM) and gene significance (GS). Functional enrichment analysis (Gene Ontology, Kyoto Encyclopedia of Genes and Genomes) of core_lncRNA target genes (co-localization and co-expression) was performed to screen potential lncRNAs regulating milk fat metabolism and further construct an interactive regulatory network of lipid metabolism-related competing endogenous RNAs (ceRNAs). A total of 4876 lncRNAs were used to construct the WGCNA. The MEdarkturquoise module among the 19 modules obtained was significantly associated with MFP (r = 0.78, p-value <0.05) and contained 64 core_lncRNAs (MM > 0.8, GS > 0.4). Twenty-four lipid metabolism-related lncRNAs were identified by core_lncRNA target gene enrichment analysis. TCONS_00054233, TCONS_00152292, TCONS_00048619, TCONS_00033839, TCONS_00153791 and TCONS_00074642 were key candidate lncRNAs for regulating milk fat synthesis. The 22 ceRNAs most likely to be involved in milk fat metabolism were constructed by interaction network analysis, and TCONS_00133813 and bta-miR-2454-5p were located at the network's core. TCONS_00133813_bta-miR-2454-5p_TNFAIP3, TCONS_00133813_bta-miR-2454-5p_ARRB1 and TCONS_00133813_bta-miR-2454-5p_PIK3R1 are key candidate ceRNAs associated with milk fat metabolism. This study provides a framework for the co-expression module of MFP-related lncRNAs in ruminants, identifies several major lncRNAs and ceRNAs that influence milk fat synthesis, and provides a new understanding of the complex biology of milk fat synthesis in dairy cows.

Two Arabidopsis MYB-SHAQKYF transcription repressors regulate leaf wax biosynthesis via transcriptional suppression on DEWAX.

Plant cuticular wax accumulation limits nonstomatal transpiration and is regulated by external environmental stresses. DEWAX (DECREASE WAX BIOSYNTHESIS) plays a vital role in diurnal wax biosynthesis. However, how DEWAX expression is controlled and the molecular mechanism of wax biosynthesis regulated by the diurnal cycle remains largely unknown. Here, we identified two Arabidopsis MYB-SHAQKYF transcription factors, MYS1 and MYS2, as new regulators in wax biosynthesis and drought tolerance. Mutations of both MYS1 and MYS2 caused significantly reduced leaf wax, whereas overexpression of MYS1 or MYS2 increased leaf wax biosynthesis and enhanced drought tolerance. Our results demonstrated that MYS1 and MYS2 act as transcription repressors and directly suppress DEWAX expression via ethylene response factor-associated amphiphilic repression motifs. Genetic interaction analysis with DEWAX, SPL9 (SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 9), and CER1 (ECERIFERUM 1) in wax biosynthesis and under drought stresses demonstrated that MYS1 and MYS2 act upstream of the DEWAX-SPL9 module, thus regulating CER1 expression. Expression analysis suggested that the diurnal expression pattern of DEWAX is partly regulated by MYS1 and MYS2. Our findings demonstrate the roles of two unidentified transcription repressors, MYS1 and MYS2, in wax biosynthesis and provide insights into the mechanism of diurnal cycle-regulated wax biosynthesis.

Weighted Gene Co-Expression Network Analysis Identifies Key Modules and Central Genes Associated With Bovine Subcutaneous Adipose Tissue.

Background: Fat deposition is an important economic trait in livestock and poultry production. However, the relationship between various genes and signal pathways of fat deposition is still unclear to a large extent. The purpose of this study is to analyze the potential molecular targets and related molecular pathways in bovine subcutaneous adipose tissue. Results: We downloaded the GSE116775 microarray dataset from Gene Expression Omnibus (GEO). The weighted gene co-expression network (WGCNA) was used to analyze the gene expression profile, and the key gene modules with the highest correlation with subcutaneous adipose tissue were identified, and the functional enrichment of the key modules was analyzed. Then, the "real" Hub gene was screened by in-module analysis and protein-protein interaction network (PPI), and its expression level in tissue samples and adipocytes was verified. The study showed that a total of nine co-expression modules were identified, and the number of genes in these modules ranged from 101 to 1,509. Among them, the blue module is most closely related to subcutaneous adipose tissue, containing 1,387 genes. These genes were significantly enriched in 10 gene ontologies including extracellular matrix organization, biological adhesion, and collagen metabolic process, and were mainly involved in pathways including ECM-receptor interaction, focal adhesion, cAMP signaling pathway, PI3K-AKT signaling pathway, and regulation of lipolysis in adipocytes. In the PPI network and coexpression network, five genes (CAV1, ITGA5, COL5A1, ABL1, and HSPG2) were identified as "real" Hub genes. Analysis of Hub gene expression by dataset revealed that the expression of these Hub genes was significantly higher in subcutaneous adipose tissue than in other tissues. In addition, real-time fluorescence quantitative PCR (qRT-PCR) analysis based on tissue samples and adipocytes also confirmed the above results. Conclusion: In this study, five key genes related to subcutaneous adipose tissue were discovered, which laid a foundation for further study of the molecular regulation mechanism of subcutaneous adipose tissue development and adipose deposition.

Identification of Key Genes Associated With Early Calf-Hood Nutrition in Subcutaneous and Visceral Adipose Tissues by Co-Expression Analysis.

Background: Substantive evidence has confirmed that nutrition state is associated with health risk and the onset of pubertal and metabolic profile. Due to heterogeneity, adipose tissues in different anatomical positions tend to show various metabolic mechanisms for nutrition. To date, the complicated molecular mechanisms of early calf-hood nutrition on bovine adipose tissue are still largely unknown. This study aimed to identify key genes and functionally enriched pathways associated with early calf-hood nutrition in visceral and subcutaneous adipose tissue. Results: The RNA-seq data of visceral and subcutaneous adipose tissues of calves feeding on low and high dietary nutrition for more than 100 days were downloaded and analyzed by weighted gene co-expression network analysis (WGCNA). Two modules that positively associated with a low plane of nutrition diet and two modules with a high plane of nutrition diet were identified in the subcutaneous adipose tissue. The blue and yellow modules, most closely associated with low and high nutrition, were selected for the functional enrichment analysis and exploration of hub genes. The results showed that genes in the blue module were significantly enriched in pathways that related to fat metabolism, reproduction, and cell communication. Genes in the yellow module were enriched in pathways related to fat metabolism, reproduction, cell proliferation, and senescence. Meanwhile, the blue and brown modules in visceral adipose tissue were most closely associated with low and high nutrition, respectively. Notably, genes of the blue module were significantly enriched in pathways related to substance metabolism, and genes in the brown module were significantly enriched in energy metabolism and disease pathways. Finally, key genes in subcutaneous adipose tissue for low nutrition (PLCG1, GNA11, and ANXA5) and high nutrition (BUB1B, ASPM, RRM2, PBK, NCAPG, and MKI67), and visceral adipose tissue for low nutrition (RPS5, RPL4, RPL14, and RPLP0) and high nutrition (SDHA and AKT1) were obtained and verified. Conclusion: The study applied WGCNA to identify hub genes and functionally enriched pathways in subcutaneous and visceral adipose tissue and provided a basis for studying the effect of early calf-hood nutrition on the two adipose tissue types.