Inferring upstream regulatory genes of FOXP3 in human regulatory T cells from time-series transcriptomic data.

The discovery of upstream regulatory genes of a gene of interest still remains challenging. Here we applied a scalable computational method to unbiasedly predict candidate regulatory genes of critical transcription factors by searching the whole genome. We illustrated our approach with a case study on the master regulator FOXP3 of human primary regulatory T cells (Tregs). While target genes of FOXP3 have been identified, its upstream regulatory machinery still remains elusive. Our methodology selected five top-ranked candidates that were tested via proof-of-concept experiments. Following knockdown, three out of five candidates showed significant effects on the mRNA expression of FOXP3 across multiple donors. This provides insights into the regulatory mechanisms modulating FOXP3 transcriptional expression in Tregs. Overall, at the genome level this represents a high level of accuracy in predicting upstream regulatory genes of key genes of interest.

Nano zinc oxide mediated resuscitation of aged Cajanus cajan via modulating aquaporin, cell cycle regulatory genes and hormonal responses.

KEY MESSAGE: Nanoparticle pretreatment improved the health of aged Cajanus cajan seeds viz., regulation of redox status, gene expression, and restoration of hormonal homeostasis. Ageing deteriorates the quality of seeds by lowering their vigor and viability, and terminating with loss of germination. These days, nanotechnology has been seen to revolutionize the agricultural sectors, and particularly nano zinc oxide (nZnO) has gained considerable interests due to its distinctive properties. The aim of the present work was to decipher the possibilities of using nZnO to rejuvenate accelerated aged (AA) seeds of Cajanus cajan. Both chemically (CnZnO) and green (GnZnO; synthesized using Moringa oleifera) fabricated nZnOs were characterized via standard techniques to interpret their purity, size, and shape. Experimental results revealed erratic germination with a decline in viability and membrane stability as outcomes of reactive oxygen intermediate (ROI) buildup in AA seeds. Application of nZnO substantially rebated the accrual of ROI, along with enhanced production of antioxidants, α-amylase activity, total sugar, protein and DNA content. Higher level of zinc was assessed qualitatively/ histologically and quantitatively in nZnO pulsed AA seeds, supporting germination without inducing toxicity. Meantime, augmentation in the gibberellic acid with a simultaneous reduction in the abscisic acid level were noted in nZnO invigorated seeds than that determined in the AA seeds, suggesting possible involvement of ROI in hormonal signalling. Furthermore, nZnO-subjected AA seeds unveiled differential expression of aquaporins and cell cycle regulatory genes. Summarizing, among CnZnO and GnZnO, later one holds better potential for a revival of AA seeds of Cajanus cajan by providing considerable tolerance against ageing-associated deterioration via recouping the cellular redox homeostasis, hormonal signaling, and alteration in expression patterns of aquaporin and cell cycle regulatory genes.

Construction and verification of a histone deacetylases-related prognostic signature model for colon cancer.

Histone deacetylases (HDACs) contribute significantly to the initiation, progression, and prognosis of colorectal adenocarcinoma (COAD). Additionally, HDACs regulate the tumor microenvironment, immune escape, and tumor stem cells, and are closely linked to COAD prognosis. We developed a prognostic model for COAD that incorporates HDACs to evaluate their specific roles. The COAD dataset containing clinical and mutation data was collected using the TCGA and GEO databases to obtain genes associated with HDAC. LASSO analysis and univariate and multivariate Cox regression analysis were used to determine the presence of prognostic genes. Multivariate Cox analysis was also used to determine risk scores for HDAC-related features. Furthermore, genomic alterations, immune infiltration, and drug response were compared between high- and low-risk groups. Cellular experiments validated the potential regulatory role of BRD3 on COAD proliferation, migration, and apoptosis. The median risk scores, calculated based on the characteristics, demonstrated a more significant prognostic improvement in patients in the low-risk group. Furthermore, HDAC-related features were identified as important independent prognostic factors for patients with COAD. Additionally, genomic mutation status, immune infiltration, and function, as well as response to immunotherapy and chemotherapy, were found to be associated with risk scores. Subgroup analyses indicate that anti-PD-1 therapy may be beneficial for patients in the low-risk group. Additionally, a decrease in risk score was associated with a decrease in immune infiltration. Finally, HCT116 and HT29 cells exhibited inhibition of BRD3 gene proliferation and migration, as well as promotion of apoptosis. In patients with COAD, HDAC-related characteristics may be useful in predicting survival and selecting treatment.

Small Non-Coding RNAs and Their Role in Locoregional Metastasis and Outcomes in Early-Stage Breast Cancer Patients.

Deregulation of small non-coding RNAs (sncRNAs) has been associated with the onset of metastasis. We evaluated the expression of sncRNAs in patients with early-stage breast cancer, performing RNA sequencing in 60 patients for whom tumor and sentinel lymph node (SLN) samples were available, and conducting differential expression, gene ontology, enrichment and survival analyses. Sequencing annotation classified most of the sncRNAs into small nucleolar RNA (snoRNAs, 70%) and small nuclear RNA (snRNA, 13%). Our results showed no significant differences in sncRNA expression between tumor or SLNs obtained from the same patient. Differential expression analysis showed down-regulation (n = 21) sncRNAs and up-regulation (n = 2) sncRNAs in patients with locoregional metastasis. The expression of SNHG5, SNORD90, SCARNA2 and SNORD78 differentiated luminal A from luminal B tumors, whereas SNORD124 up-regulation was associated with luminal B HER2+ tumors. Discriminating analysis and receiver-operating curve analysis revealed a signature of six snoRNAs (SNORD93, SNORA16A, SNORD113-6, SNORA7A, SNORA57 and SNORA18A) that distinguished patients with locoregional metastasis and predicted patient outcome. Gene ontology and Reactome pathway analysis showed an enrichment of biological processes associated with translation initiation, protein targeting to specific cell locations, and positive regulation of Wnt and NOTCH signaling pathways, commonly involved in the promotion of metastases. Our results point to the potential of several sncRNAs as surrogate markers of lymph node metastases and patient outcome in early-stage breast cancer patients. Further preclinical and clinical studies are required to understand the biological significance of the most significant sncRNAs and to validate our results in a larger cohort of patients.

Iron metabolism disorder and multiple sclerosis: a comprehensive analysis.

Background: Multiple sclerosis (MS) is the most common chronic inflammatory disease of the central nervous system. Currently, the pathological mechanisms of MS are not fully understood, but research has suggested that iron metabolism disorder may be associated with the onset and clinical manifestations of MS. Methods and materials: The study utilized publicly available databases and bioinformatics techniques for gene expression data analysis, including differential expression analysis, weighted correlation network analysis, gene enrichment analysis, and construction of logistic regression models. Subsequently, Mendelian randomization was used to assess the causal relationship between different iron metabolism markers and MS. Results: This study identified IREB2, LAMP2, ISCU, ATP6V1G1, ATP13A2, and SKP1 as genes associated with multiple sclerosis (MS) and iron metabolism, establishing their multi-gene diagnostic value for MS with an AUC of 0.83. Additionally, Mendelian randomization analysis revealed a potential causal relationship between transferrin saturation and MS (p=2.22E-02; OR 95%CI=0.86 (0.75, 0.98)), as well as serum transferrin and MS (p=2.18E-04; OR 95%CI=1.22 (1.10, 1.36)). Conclusion: This study comprehensively explored the relationship between iron metabolism and MS through integrated bioinformatics analysis and Mendelian randomization methods. The findings provide important insights for further research into the role of iron metabolism disorder in the pathogenesis of MS and offer crucial theoretical support for the treatment of MS.

Gene Expression Patterns Associated with Survival in Glioblastoma.

The aim of this study was to investigate gene expression alterations associated with overall survival (OS) in glioblastoma (GBM). Using the Nanostring nCounter platform, we identified four genes (COL1A2, IGFBP3, NGFR, and WIF1) that achieved statistical significance when comparing GBM with non-neoplastic brain tissue. The four genes were included in a multivariate Cox Proportional Hazard model, along with age, extent of resection, and O6-methylguanine-DNA methyltransferase (MGMT) promotor methylation, to create a unique glioblastoma prognostic index (GPI). The GPI score inversely correlated with survival: patient with a high GPI had a median OS of 7.5 months (18-month OS = 9.7%) whereas patients with a low GPI had a median OS of 20.1 months (18-month OS = 54.5%; log rank p-value = 0.004). The GPI score was then validated in 188 GBM patients from The Cancer Genome Atlas (TCGA) from a national data base; similarly, patients with a high GPI had a median OS of 10.5 months (18-month OS = 12.4%) versus 16.9 months (18-month OS = 41.5%) for low GPI (log rank p-value = 0.0003). We conclude that this novel mRNA-based prognostic index could be useful in classifying GBM patients into risk groups and refine prognosis estimates to better inform treatment decisions or stratification into clinical trials.

Tumor lymphangiogenesis index reveals the immune landscape and immunotherapy response in lung adenocarcinoma.

Background: Lymphangiogenesis (LYM) has an important role in tumor progression and is strongly associated with tumor metastasis. However, the clinical application of LYM has not progressed as expected. The potential value of LYM needs to be further developed in lung adenocarcinoma (LUAD) patients. Methods: The Sequencing data and clinical characteristics of LUAD patients were downloaded from The Cancer Genome Atlas and GEO databases. Multiple machine learning algorithms were used to screen feature genes and develop the LYM index. Immune cell infiltration, immune checkpoint expression, Tumor Immune Dysfunction and Exclusion (TIDE) algorithm and drug sensitivity analysis were used to explore the correlation of LYM index with immune profile and anti-tumor therapy. Results: We screened four lymphangiogenic feature genes (PECAM1, TIMP1, CXCL5 and PDGFB) to construct LYM index based on multiple machine learning algorithms. We divided LUAD patients into the high LYM index group and the low LYM index group based on the median LYM index. LYM index is a risk factor for the prognosis of LUAD patients. In addition, there was a significant difference in immune profile between high LYM index and low LYM index groups. LUAD patients in the low LYM index group seemed to benefit more from immunotherapy based on the results of TIDE algorithm. Conclusion: Overall, we confirmed that the LYM index is a prognostic risk factor and a valuable predictor of immunotherapy response in LUAD patients, which provides new evidence for the potential application of LYM.

Identification of PTPN12 Phosphatase as a Novel Negative Regulator of Hippo Pathway Effectors YAP/TAZ in Breast Cancer.

The Hippo pathway plays crucial roles in governing various biological processes during tumorigenesis and metastasis. Within this pathway, upstream signaling stimuli activate a core kinase cascade, involving MST1/2 and LATS1/2, that subsequently phosphorylates and inhibits the transcriptional co-activators YAP and its paralog TAZ. This inhibition modulates the transcriptional regulation of downstream target genes, impacting cell proliferation, migration, and death. Despite the acknowledged significance of protein kinases in the Hippo pathway, the regulatory influence of protein phosphatases remains largely unexplored. In this study, we conducted the first gain-of-functional screen for protein tyrosine phosphatases (PTPs) regulating the Hippo pathway. Utilizing a LATS kinase biosensor (LATS-BS), a YAP/TAZ activity reporter (STBS-Luc), and a comprehensive PTP library, we identified numerous novel PTPs that play regulatory roles in the Hippo pathway. Subsequent experiments validated PTPN12, a master regulator of oncogenic receptor tyrosine kinases (RTKs), as a previously unrecognized negative regulator of the Hippo pathway effectors, oncogenic YAP/TAZ, influencing breast cancer cell proliferation and migration. In summary, our findings offer valuable insights into the roles of PTPs in the Hippo signaling pathway, significantly contributing to our understanding of breast cancer biology and potential therapeutic strategies.

Bioinformatics-based analysis of the relationship between plasminogen regulatory genes and photoaging.

BACKGROUND: Ultraviolet radiation causes skin photoaging by producing a variety of enzymes, which impact both skin health and hinder beauty. Currently, the early diagnosis and treatment of photoaging remain a challenge. Bioinformatics analysis has strong advantages in exploring core genes and the biological pathways of photoaging. AIMS: To screen and validate key risk genes associated with plasminogen in photoaging and to identify potential target genes for photoaging. METHODS: Two human transcriptome datasets were obtained by searching the Gene Expression Omnibus (GEO) database, and the mRNAs in the GSE131789 dataset were differentially analyzed, and then the weighted gene co-expression network analysis (WGCNA) was performed to find out the strongest correlations. Template genes, interaction analysis of differentially expressed genes (DEGs), modular genes with the most WGCNA correlations, and genecard database genes related to plasminogen were performed, and further Kyoto genes and Genome Encyclopedia (KEGG) pathway analysis. Two different algorithms, least absolute shrinkage and selection operator (LASSO) and support vector machines-recursive feature elimination (SVM-RFE), were used to find key genes. Then the data set (GSE206495) was validated and analyzed. Real-time PCR was performed to validate the expression of key genes through in vitro cellular experiments. RESULTS: IFI6, IFI44L, HRSP12, and BMP4 were screened from datasets as key genes for photoaging and further analysis showed that these genes have significant diagnostic value for photoaging. CONCLUSION: IFI6, IFI44L, HRSP12, and BMP4 play a key role in the pathogenesis of photoaging, and serve as promising potential predictive biomarkers for photoaging.

[Preclinical diagnostics of von Hippel-Lindau syndrome in a child].

The description of the child aged 5 months with the von Hippel-Lindau syndrome without any manifestations of this syndrome is presented. The reason for the molecular genetic examination was the presence of cases of this syndrome in the family (mother and sister). The heterozygous variant c.355T>C p.F119L was found in the VHL gene. An objective examination revealed no pathology. A comprehensive laboratory and instrumental examination aimed at searching for components of the von Hippel-Lindau syndrome, including a blood test for metanephrines and normetanephrines, ultrasound of the abdominal organs, examination of the fundus, also did not reveal any abnormalities. Given the results of molecular genetic diagnosis, the child remains under observation and will undergo regular examinations to identify components of the von Hippel-Lindau syndrome, including blood/urine tests for normetanephrines.

Clinical and Genetic Analysis of Multiple Idiopathic Cervical Root Resorption.

OBJECTIVE: To explore the genetic background and clinical phenotypes of multiple idiopathic cervical root resorption (MICRR) in a Chinese family. METHODS: The proband and his three family members were clinically examined and had radiographs taken with a radiovisiography (RVG) system and CBCT to define the diagnosis of MICRR. Genomic DNA (gDNA) was extracted from peripheral blood samples of the patient, his father, mother and younger sister for whole exome sequencing (WES). The pathogenicity of rare variants with minor allele frequency (MAF) less than 0.005 were analysed following possible inheritance patterns, predicted results from 12 software programs, the American College of Medical Genetics (ACMG) 2015 criteria, and information from ClinVar, OMIM and HGMD databases as well as gene function. RESULTS: The proband presented the typical MICRR phenotypes such as thin cervical pulp wall and apple core-like lesions in radiographs. Following the recessive inheritance pattern, WES analysis identified SHROOM2, SYTL5, MAGED1 and FLNA with a higher chance of causing MICRR. Four genes with compound heterozygous variants and another 27 genes with de novo variants either in autosomal-dominant or autosomal-recessive pattern were also found to have the potential pathogenicity. CONCLUSION: A total of 35 novel potential pathogenic genes were found to be associated with MICRR from a Chinese family through WES. The new genetic background of MICRR may be helpful for clinical and molecular diagnosis.

The integrated bioinformatic analysis identifies immune microenvironment-related potential biomarkers for patients with gestational diabetes mellitus.

Background: Gestational diabetes mellitus (GDM), a transient disease, may lead to short- or long-term adverse influences on maternal and fetal health. Therefore, its potential functions, mechanisms and related molecular biomarkers must be comprehended for the control, diagnosis and treatment of GDM. Methods: The differentially expressed genes (DEGs) were identified using GSE49524 and GSE87295 associated with GDM from the Gene Expression Omnibus database, followed by function enrichment analysis, protein-protein interactions network construction, hub DEGs mining, diagnostic value evaluation and immune infiltration analysis. Finally, hub DEGs, the strongest related to immune infiltration, were screened as immune-related biomarkers. Results: A hundred and seven DEGs were identified between patients with GDM and healthy individuals. Six hub genes with high diagnostic values, including ALDH1A1, BMP4, EFNB2, MME, PLAUR and SLIT2, were identified. Among these, two immune-related genes (PLAUR and SLIT2) with the highest absolute correlation coefficient were considered immune-related biomarkers in GDM. Conclusion: Our study provides a comprehensive analysis of GDM, which would provide a foundation for the development of diagnosis and treatment of GDM.

Whole genome discovery of regulatory genes responsible for the response of chicken to heat stress.

Long noncoding RNAs (lncRNAs) are functional bridges connecting the genome with phenotypes by interacting with DNA, mRNA, and proteins. Using publically available acute heat stress (AHS)-related RNA-seq data, we discovered novel lncRNAs and tested their association with AHS along with ~ 8800 known lncRNAs and ~ 28,000 mRNA transcripts. Our pipeline discovered a total of 145 potentially novel-lncRNAs. One of them (Fishcomb_p-value = 0.06) along with another novel transcript (annotated as protein-coding; Fishcomb_p-value = 0.03) were identified as significantly associated with AHS. We found five known-lncRNAs and 134 mRNAs transcripts that were significantly associated with AHS. Four novel lncRNAs interact cis-regulated with 12 mRNA transcripts and are targeted by 11 miRNAs. Also six meta-lncRNAs associate with 134 meta-mRNAs through trans-acting co-expression, each targeted by 15 and 216 miRNAs, respectively. Three of the known-lncRNAs significantly co-expressed with almost 97 of the significant mRNAs (Pearson correlation p-value < 0.05). We report the mentioned three known-lncRNAs (ENSGALT00000099876, ENSGALT00000107573, and ENSGALT00000106323) as the most, significantly regulatory elements of AHS in chicken. It can be concluded that in order to alleviate the adverse effects of AHS on chicken, the manipulation of the three regulatory lncRNAs could lead to a more desirable result than the manipulation of the most significant mRNAs.

The immune cell infiltration-associated molecular subtypes and gene signature predict prognosis for osteosarcoma patients.

Host immune dysregulation involves in the initiation and development of osteosarcoma (OS). However, the exact role of immune cells in OS remains unknown. We aimed to distinguish the molecular subtypes and establish a prognostic model in OS patients based on immunocyte infiltration. The gene expression profile and corresponding clinical feature of OS patients were obtained from TARGET and GSE21257 datasets. MCP-counter and univariate Cox regression analyses were applied to identify immune cell infiltration-related molecular subgroups. Functional enrichment analysis and immunocyte infiltration analysis were performed between two subgroups. Furthermore, Cox regression and LASSO analyses were performed to establish the prognostic model for the prediction of prognosis and metastasis in OS patients. The subgroup with low infiltration of monocytic lineage (ML) was related to bad prognosis in OS patients. 435 DEGs were screened between the two subgroups. Functional enrichment analysis revealed these DEGs were involved in immune- and inflammation-related pathways. Three important genes (including TERT, CCDC26, and IL2RA) were identified to establish the prognostic model. The risk model had good prognostic performance for the prediction of metastasis and overall survival in OS patients. A novel stratification system was established based on ML-related signature. The risk model could predict the metastasis and prognosis in OS patients. Our findings offered a novel sight for the prognosis and development of OS.

Genetic architecture of the structural connectome.

Myelinated axons form long-range connections that enable rapid communication between distant brain regions, but how genetics governs the strength and organization of these connections remains unclear. We perform genome-wide association studies of 206 structural connectivity measures derived from diffusion magnetic resonance imaging tractography of 26,333 UK Biobank participants, each representing the density of myelinated connections within or between a pair of cortical networks, subcortical structures or cortical hemispheres. We identify 30 independent genome-wide significant variants after Bonferroni correction for the number of measures studied (126 variants at nominal genome-wide significance) implicating genes involved in myelination (SEMA3A), neurite elongation and guidance (NUAK1, STRN, DPYSL2, EPHA3, SEMA3A, HGF, SHTN1), neural cell proliferation and differentiation (GMNC, CELF4, HGF), neuronal migration (CCDC88C), cytoskeletal organization (CTTNBP2, MAPT, DAAM1, MYO16, PLEC), and brain metal transport (SLC39A8). These variants have four broad patterns of spatial association with structural connectivity: some have disproportionately strong associations with corticothalamic connectivity, interhemispheric connectivity, or both, while others are more spatially diffuse. Structural connectivity measures are highly polygenic, with a median of 9.1 percent of common variants estimated to have non-zero effects on each measure, and exhibited signatures of negative selection. Structural connectivity measures have significant genetic correlations with a variety of neuropsychiatric and cognitive traits, indicating that connectivity-altering variants tend to influence brain health and cognitive function. Heritability is enriched in regions with increased chromatin accessibility in adult oligodendrocytes (as well as microglia, inhibitory neurons and astrocytes) and multiple fetal cell types, suggesting that genetic control of structural connectivity is partially mediated by effects on myelination and early brain development. Our results indicate pervasive, pleiotropic, and spatially structured genetic control of white-matter structural connectivity via diverse neurodevelopmental pathways, and support the relevance of this genetic control to healthy brain function.

Comprehensive analysis indicated that NDE1 is a potential biomarker for pan-cancer and promotes bladder cancer progression.

BACKGROUND: The nuclear distribution E homologue 1 (NDE1) is a crucial dynein binding partner. The NDE1 protein has the potential to disrupt the normal functioning of centrosomes, leading to a compromised ability to generate spindles and ensure precise separation of chromosomes during cell division. The potential consequences of this phenomenon include genomic instability, malignant transformation and the proliferation of neoplastic growths. However, studies examining the connection between NDE1 and cancer is still very rare. METHODS: The expression level, prognostic impact, gene change, DNA methylation, protein interaction, mRNA m6A modification, ceRNA network, associated gene and function enrichment, and immune-related effects of NDE1 in pan-cancer were examined using a range of online analytic tools and the R software package. The CCK-8 test, transwell assay, scratch assay and colony formation assay were used to confirm the effects of NDE1 on the proliferation, invasion and metastasis of bladder cancer cells. RESULTS: Numerous tumour types have elevated NDE1, which is linked to a bad prognosis. NDE1 is an excellent diagnostic tool for many different types of cancer. Numerous malignancies have been linked to genetic changes in NDE1. NDE1 was connected to TMB, MSI, several immunological checkpoint genes and immune cell infiltration. NDE1 is linked to a number of immunological subtypes. NDE1 could affect how well immunotherapy works to treat different types of cancer. NDE1 was mostly associated with cell cycle, chromosomal segregation, DNA replication and mitotic segregation, according to GO and KEGG analyses. NDE1 physically binds to PAFAH1B1 and DCTN1, respectively. The proliferation, invasion and metastasis of bladder cancer cells may be prevented by NDE1 knockdown. Furthermore, knockdown of NDE1 promoted the apoptosis of bladder cancer cells. CONCLUSION: High expression of NDE1 is present in a variety of tumours, which is linked to a bad prognosis for cancer. Knockdown of NDE1 inhibited the proliferation, invasion and metastasis of bladder cancer cells, and promoted the apoptosis. For a number of malignancies, NDE1 may be a biomarker for immunotherapy and prognosis.

A Novel Signature Based on m6A Regulator-Mediated Genes Along Glycolytic Pathway Predicts Prognosis and Immunotherapy Responses of Gastric Cancer Patients.

N6-methyladenosine (m6A) modification is involved in many aspects of gastric cancer (GC). Moreover, m6A and glycolysis-related genes (GRGs) play important roles in immunotherapeutic and prognostic implication of GC. However, GRGs involved in m6A regulation have never been analyzed comprehensively in GC. Herein, the study aims to identify and validate a novel signature based on m6A-related GRGs in GC patients. Therefore, a m6A-related GRGs signature is established, which can predict the survival of patients with GC and remain an independent prognostic factor in multivariate analyses. Clinical significance of the model is well validated in internal cohort and independent validation cohort. In addition, the expression levels of risk model-related GRGs in clinical samples are validated. Consistent with the database results, all model genes are up-regulated in expression except DCN. After regrouping the patients based on this risk model, the study can effectively distinguish between them in respect to immune-cell infiltration microenvironment and immunotherapeutic response. Additionally, candidate drugs targeting risk model-related GRGs are confirmed. Finally, a nomogram combining risk scores and clinical parameters is created, and calibration plots show that the nomogram can accurately predict survival. This risk model can serve as a reliable assessment tool for predicting prognosis and immunotherapeutic responses in GC patients.

Identification of candidate regulatory genes for intramuscular fatty acid composition in pigs by transcriptome analysis.

BACKGROUND: Intramuscular fat (IMF) content and its fatty acid (FA) composition are typically controlled by several genes, each with a small effect. In the current study, to pinpoint candidate genes and putative regulators involved in FA composition, we performed a multivariate integrative analysis between intramuscular FA and transcriptome profiles of porcine longissimus dorsi (LD) muscle. We also carried out a combination of network, regulatory impact factor (RIF), in silico prediction of putative target genes, and functional analyses to better support the biological relevance of our findings. RESULTS: For this purpose, we used LD RNA-Seq and intramuscular FA composition profiles of 129 Iberian × Duroc backcrossed pigs. We identified 378 correlated variables (13 FA and 365 genes), including six FA (C20:4n-6, C18:2n-6, C20:3n-6, C18:1n-9, C18:0, and C16:1n-7) that were among the most interconnected variables in the predicted network. The detected FA-correlated genes include genes involved in lipid and/or carbohydrate metabolism or in regulation of IMF deposition (e.g., ADIPOQ, CHUK, CYCS, CYP4B1, DLD, ELOVL6, FBP1, G0S2, GCLC, HMGCR, IDH3A, LEP, LGALS12, LPIN1, PLIN1, PNPLA8, PPP1R1B, SDR16C5, SFRP5, SOD3, SNW1, and TFRC), meat quality (GALNT15, GOT1, MDH1, NEU3, PDHA1, SDHD, and UNC93A), and transport (e.g., EXOC7 and SLC44A2). Functional analysis highlighted 54 over-represented gene ontology terms, including well-known biological processes and pathways that regulate lipid and carbohydrate metabolism. RIF analysis suggested a pivotal role for six transcription factors (CARHSP1, LBX1, MAFA, PAX7, SIX5, and TADA2A) as putative regulators of gene expression and intramuscular FA composition. Based on in silico prediction, we identified putative target genes for these six regulators. Among these, TADA2A and CARHSP1 had extreme RIF scores and present novel regulators in pigs. In addition, the expression of TADA2A correlated (either positively or negatively) with C20:4n-6, C18:2n-6, C20:3n-6, C18:1n-9, and that of CARHSP1 correlated (positively) with the C16:1n-7 lipokine. We also found that these two transcription factors share target genes that are involved in lipid metabolism (e.g., GOT1, PLIN1, and TFRC). CONCLUSIONS: This integrative analysis of muscle transcriptome and intramuscular FA profile revealed valuable information about key candidate genes and potential regulators for FA and lipid metabolism in pigs, among which some transcription factors are proposed to control gene expression and modulate FA composition differences.

Integrated Transcriptomic and Proteomic Analysis Identifies Novel Regulatory Genes Associated with Plant Growth Regulator-Induced Astringency in Grape Berries.

Astringency influences the sensory characteristics and flavor quality of table grapes. We tested the astringency sensory attributes of berries and investigated the concentration of flavan-3-ols/proanthocyanidins (PAs) in skins after the application of the plant growth regulators CPPU and GA3 to the flowers and young berries of the "Summer Black" grape. Our results showed that CPPU and GA3 applications increase sensory astringency perception scores and flavan-3-ol/proanthocyanidin concentrations. Using integrated transcriptomic and proteomic analysis, differentially expressed transcripts and proteins associated with growth regulator treatment were identified, including those for flavonoid biosynthesis that contribute to the changes in sensory astringency levels. Transient overexpression of candidate astringency-related regulatory genes in grape leaves revealed that VvWRKY71, in combination with VvMYBPA1 and VvMYC1, could promote the biosynthesis of proanthocyanidins, while overexpression of VvNAC83 reduced the accumulation of proanthocyanidins. However, in transient promoter studies in Nicotiana benthamiana, VvWRKY71 repressed the promoter of VvMYBPA2, while VvNAC83 had no significant effect on the promoter activity of four PA-related genes, and VvMYBPA1 was shown to activate its own promoter. This study provides new insights into the molecular mechanisms of sensory astringency formation induced by plant growth regulators in grape berries.

Genome-wide identification of the N6-methyladenosine regulatory genes reveals NtFIP37B increases drought resistance of tobacco (Nicotiana tabacum L.).

BACKGROUND: N6-methyladenosine (m6A) is one of the common internal RNA modifications found in eukaryotes. The m6A modification can regulate various biological processes in organisms through the modulation of alternative splicing, alternative polyadenylation, folding, translation, localization, transport, and decay of multiple types of RNA, without altering the nucleotide sequence. The three components involved in m6A modification, namely writer, eraser, and reader, mediate the abundance of RNA m6A modification through complex collaborative actions. Currently, research on m6A regulatory genes in plants is still in its infancy. RESULTS: In this study, we identified 52 candidate m6A regulatory genes in common tobacco (Nicotiana tabacum L.). Gene structure, conserved domains, and motif analysis showed structural and functional diversity among different subgroups of tobacco m6A regulatory genes. The amplification of m6A regulatory genes were mainly driven by polyploidization and dispersed duplication, and duplicated genes evolved through purified selection. Based on the potential regulatory network and expression pattern analysis of m6A regulatory genes, a significant number of m6A regulatory genes might play important roles in growth, development, and stress response processes. Furthermore, we have confirmed the critical role of NtFIP37B, an m6A writer gene in tobacco, in enhancing drought resistance. CONCLUSIONS: This study provides useful information for better understanding the evolution of m6A regulatory genes and the role of m6A modification in tobacco stress response, and lays the foundation for further elucidating the function of m6A regulatory genes in tobacco.