KOBAS-i: intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis

Regulation of DNA methylation during the testicular development of Shaziling pigs

DNA methylation is one of the key epigenetic regulatory mechanisms in development and spermatogenesis. However, the dynamic regulatory mechanisms of genome-wide DNA methylation during testicular development remain largely unknown. Herein, we generated a single-base resolution DNA methylome and transcriptome atlas of precocious porcine testicular tissues across three developmental stages (1, 75, and 150 days old). The results showed that the dynamic methylation patterns were directly related to the expression of the DNMT3A gene. Conjoint analysis revealed a negative regulatory pattern between promoter methylation and the positive regulation of 3'-untranslated region (3'UTR) methylation. Mechanistically, the decrease in promoter methylation affected the upregulation of meiosis-related genes, such as HORMAD1, SPO11, and SYCE1. Demethylation in the 3'UTR induced the downregulation of the INHBA gene and knockdown of INHBA inhibited cell proliferation by reducing the synthesis of activin A. These findings contribute to exploring the regulatory mechanisms of DNA methylation in testicular development.

Fatty acid metabolism-related genes are associated with flavor-presenting aldehydes in Chinese local chicken

Aldehydes are primary volatile organic compounds (VOCs) in local Chinese chicken meat and contribute green grass, fatty, citrus, and bitter almond aromas to chicken meat. To understand the genetic basis of these aldehyde VOC aromas, we used approximately 500 Chinese Jingxing Yellow (JXY) chickens to conduct genome-wide association studies (GWAS) on the flavor traits with the data of single nucleotide polymorphisms (SNPs) and insertions and deletions (INDELs). In total, 501 association variants (253 SNPs and 248 INDELs) were found to be suggestively (SNPs: p-value < 2.77e-06 and INDELs: p-value < 3.78e-05) associated with total aldehydes (the sum of nine aldehydes), hexanal, heptanal, benzaldehyde, (E,E)-2,4-nonadienal, octanal, (E)-2-decenal, nonanal, decanal, and octadecanal. Of them, six SNPs and 23 INDELs reached a genome-wide significance level (SNPs: p-value < 1.38e-07 and INDELs: p-value < 1.89e-06). Potential candidate aldehyde genes were functionally annotated for lipid metabolism, especially fatty acid-related pathways and phospholipid-related gene ontology (GO) terms. Moreover, the GWAS analysis of total aldehydes, hexanal, and nonanal generated the most significant signals, and phenotypic content differed between different genotypes at candidate gene-related loci. For total aldehydes and hexanal traits, candidate genes were annotated based on the significant and suggestive variants on chromosomes 3 and 8 with highly polymorphic linkage blocks. The following candidate genes were also identified: GALM, MAP4K3, GPCPD1, RPS6KA2, CRLS1, ASAP1, TRMT6, SDC1, PUM2, ALDH9A1, MGST3, GMEB1, MECR, LDLRAP1, GPAM and ACSL5. We also found that polyunsaturated fatty acids (PUFAs) (C18:2n6c linoleic acid and C18:3n3 linolenic acid) were significantly correlated with total aldehydes and hexanal contents. PUFAs are important aldehyde precursors, and consistently, our results suggested that candidate genes involved in fatty acid pathways and phospholipid GO terms were identified in association loci. This work provides an understanding of the genetic basis of aldehyde formation, which is a key flavor-forming compound.

Utilization of hypoxia-derived gene signatures to predict clinical outcomes and immune checkpoint blockade therapy responses in prostate cancer

Background: Increasing evidences show a clinical significance in the interaction between hypoxia and prostate cancer. However, reliable prognostic signatures based on hypoxia have not been established yet.

Methods: We screened hypoxia-related gene modules by weighted gene co-expression network analysis (WGCNA) and established a hypoxia-related prognostic risk score (HPRS) model by univariate Cox and LASSO-Cox analyses. In addition, enriched pathways, genomic mutations, and tumor-infiltrating immune cells in HPRS subgroups were analyzed and compared. HPRS was also estimated to predict immune checkpoint blockade (ICB) therapy response.

Results: A hypoxia-related 22-gene prognostic model was established. Furthermore, three independent validation cohorts showed moderate performance in predicting biochemical recurrence-free (BCR-free) survival. HPRS could be a useful tool in selecting patients who can benefit from ICB therapy. The CIBERSORT results in our study demonstrated that hypoxia might act on multiple T cells, activated NK cells, and macrophages M1 in various ways, suggesting that hypoxia might exert its anti-tumor effects by suppressing T cells and NK cells.

Conclusion: Hypoxia plays an important role in the progression of prostate cancer. The hypoxia-derived signatures are promising biomarkers to predict biochemical recurrence-free survival and ICB therapy responses in patients with prostate cancer.

GSA Central—A web platform to perform, learn, and discuss gene set analysis

Gene Set Analysis (GSA) is one of the most commonly used strategies to analyze omics data. Hundreds of GSA-related papers have been published, giving birth to a GSA field in Bioinformatics studies. However, as the field grows, it is becoming more difficult to obtain a clear view of all available methods, resources, and their quality. In this paper, we introduce a web platform called “GSA Central” which, as its name indicates, acts as a focal point to centralize GSA information and tools useful to beginners, average users, and experts in the GSA field. “GSA Central” contains five different resources: A Galaxy instance containing GSA tools (“Galaxy-GSA”), a portal to educational material (“GSA Classroom”), a comprehensive database of articles (“GSARefDB”), a set of benchmarking tools (“GSA BenchmarKING”), and a blog (“GSA Blog”). We expect that “GSA Central” will become a useful resource for users looking for introductory learning, state-of-the-art updates, method/tool selection guidelines and insights, tool usage, tool integration under a Galaxy environment, tool design, and tool validation/benchmarking. Moreover, we expect this kind of platform to become an example of a “thematic platform” containing all the resources that people in the field might need, an approach that could be extended to other bioinformatics topics or scientific fields.

Genome-Wide Association Study for Screening and Identifying Potential Shin Color Loci in Ducks

Shin color diversity is a widespread phenomenon in birds. In this study, ducks were assessed to identify candidate genes for yellow, black, and spotted tibiae. For this purpose, we performed whole-genome resequencing of an F2 population consisting of 275 ducks crossed between Runzhou crested-white ducks and Cherry Valley ducks. We obtained 12.6 Mb of single nucleotide polymorphism (SNP) data, and the three shin colors were subsequently genotyped. Genome-wide association studies (GWASs) were performed to identify candidate and potential SNPs for the three shin colors. According to the results, 2947 and 3451 significant SNPs were associated with black and yellow shins, respectively, and six potential SNPs were associated with spotted shins. Based on the SNP annotations, the MITF, EDNRB2, POU family members, and the SLC superfamily were the candidate genes regulating pigmentation. In addition, the isoforms of EDNRB2, TYR, TYRP1, and MITF-M were significantly different between the black and yellow tibiae. MITF and EDNRB2 may have synergistic roles in the regulation of melanin synthesis, and their mutations may lead to phenotypic differences in the melanin deposition between individuals. This study provides new insights into the genetic factors that may influence tibia color diversity in birds.

Discovery of Small-Molecule Degraders of the CDK9-Cyclin T1 Complex for Targeting Transcriptional Addiction in Prostate Cancer

Aberrant hyperactivation of cyclins results in carcinogenesis and therapy resistance in cancers. Direct degradation of the specific cyclin or cyclin-dependent kinase (CDK)-cyclin complex by small-molecule degraders remains a great challenge. Here, we applied the first application of hydrophobic tagging to induce degradation of CDK9-cyclin T1 heterodimer, which is required to keep productive transcription of oncogenes in cancers. LL-K9-3 was identified as a potent small-molecule degrader of CDK9-cyclin T1. Quantitative and time-resolved proteome profiling exhibited LL-K9-3 induced selective and synchronous degradation of CDK9 and cyclin T1. The expressions of androgen receptor (AR) and cMyc were reduced by LL-K9-3 in 22RV1 cells. LL-K9-3 exhibited enhanced anti-proliferative and pro-apoptotic effects compared with its parental CDK9 inhibitor SNS032 and suppressed downstream signaling of CDK9 and AR more effectively than SNS032. Moreover, LL-K9-3 inhibited AR and Myc-driven oncogenic transcriptional programs and exerted stronger inhibitory effects on several intrinsic target genes of AR than the monomeric CDK9 PROTAC (Thal-SNS032).

Long non-coding RNA Gm10561 promotes myogenesis by sponging miR-432

Skeletal myogenesis is a highly ordered process finely regulated by various factors. Long non-coding RNAs play an important regulatory role in myogenesis via multiple mechanisms. In this study, we identified the lncRNA Gm10561, which was upregulated during myogenic differentiation and is highly expressed in skeletal muscle. Knockdown of Gm10561 inhibited the proliferation and differentiation of C2C12 myoblasts in vitro and muscle growth in vivo. Overexpression of Gm10561 promoted the proliferation and differentiation of both C2C12 myoblasts and porcine muscle satellite cells. Notably, lncRNA Gm10561 is localized in the cytoplasm and competitively bound to miR-432, which directly targets MEF2C and E2F3. It was confirmed that lncRNA Gm10561 regulates the proliferation and differentiation of myoblasts by acting as a sponge of miR-432 to modulate MEF2C and E2F3 expression. Thus, the lncRNA-Gm10561-miR-432-MEF2C/E2F3 axis plays an important role in myogenesis.

Bioinformatics analysis of potential common pathogenic mechanisms for COVID-19 infection and primary Sjogren’s syndrome

Background

Accumulating evidence has revealed that the prevalence of Coronavirus 2019 (COVID-19) was significantly higher in patients with primary Sjogren’s syndrome (pSS) compared to the general population. However, the mechanism remains incompletely elucidated. This study aimed to further investigate the molecular mechanisms underlying the development of this complication.

Methods

The gene expression profiles of COVID-19 (GSE157103) and pSS (GSE40611) were downloaded from the Gene Expression Omnibus (GEO) database. After identifying the common differentially expressed genes (DEGs) for pSS and COVID-19, functional annotation, protein-protein interaction (PPI) network, module construction and hub gene identification were performed. Finally, we constructed transcription factor (TF)-gene regulatory network and TF-miRNA regulatory network for hub genes.

Results

A total of 40 common DEGs were selected for subsequent analyses. Functional analyses showed that cellular components and metabolic pathways collectively participated in the development and progression of pSS and COVID-19. Finally, 12 significant hub genes were identified using the cytoHubba plugin, including CMPK2, TYMS, RRM2, HERC5, IFI44L, IFI44, IFIT2, IFIT1, IFIT3, MX1, CDCA2 and TOP2A, which had preferable values as diagnostic markers for COVID-19 and pSS.

Conclusions

Our study reveals common pathogenesis of pSS and COVID-19. These common pathways and pivotal genes may provide new ideas for further mechanistic studies.

RNA Sequencing Reveals the Expression Profiles of circRNAs and Indicates Hsa_circ_0070562 as a Pro-osteogenic Factor in Bone Marrow-Derived Mesenchymal Stem Cells of Patients With Ankylosing Spondylitis

Recent studies have reported that circular RNAs (circRNAs) play a crucial regulatory role in a variety of human diseases. However, the roles of circRNAs in pathological osteogenesis in ankylosing spondylitis (AS) remain unclear. We conducted circRNA and miRNA expression profiling of osteogenically differentiated bone marrow-derived mesenchymal stem cells (BMSCs) of patients with AS compared with those of healthy donors (HDs) by RNA sequencing (RNA-seq). Results showed that a total of 31806 circRNAs were detected in the BMSC samples, of which 418 circRNAs were significantly differentially expressed (DE) with a fold change ≥2 and p value <0.05. Among these, 204 circRNAs were upregulated, and 214 were downregulated. GO and KEGG analyses demonstrated that the DE circRNAs were mainly involved in the regulation of biological processes of the cell matrix adhesion and the TGF-beta signaling pathway, which are closely related to AS. circRNA-miRNA interaction networks related to the TGF-beta signaling pathway were established. The results of qRT-PCR showed that has_circ_0070562 was significantly up-regulated in AS-MSCs. In vitro experiments showed that silencing of has_circ_0070562 weakened osteogenesis of AS-BMSCs. In conclusion, we identified numerous circRNAs that were dysregulated in AS-BMSCs compared with HD-BMSCs. Bioinformatic analyses suggested that these dysregulated circRNAs might play important functional roles in AS-BMSCs osteogenesis. Circ_0070562 functioned as a pro-ostegenic factor and might serve as a potential biomarker and a therapeutic target for AS.

Genome-Wide Association Analysis and Genetic Parameters for Feed Efficiency and Related Traits in Yorkshire and Duroc Pigs

Simple Summary

Genetic improvements in feed efficiency (FE) and related traits could considerably reduce pig production costs and energy consumption. Thus, we performed a genetic parameter estimation and genome-wide association study of four FE and FE-related traits, namely, average daily feed intake, average daily gain, the feed conversion ratio, and residual feed intake, of two pig breeds, Yorkshire and Duroc. The results demonstrate the genetic relationships of FE and FE-related traits with two growth traits, age and backfat thickness at 100 kg. We also identified many single-nucleotide polymorphisms (SNPs) and novel candidate genes related to these traits. In addition, we found many pathways significantly associated with FE and FE-related traits, and they are generally involved in digestive and metabolic processes. The results of this study are expected to provide a valuable reference for the genomic selection of FE and FE-related traits in pigs.

Abstract

Feed efficiency (FE) traits are key factors that can influence the economic benefits of pig production. However, little is known about the genetic architecture of FE and FE-related traits. This study aimed to identify SNPs and candidate genes associated with FE and FE-related traits, namely, average daily feed intake (ADFI), average daily gain (ADG), the feed conversion ratio (FCR), and residual feed intake (RFI). The phenotypes of 5823 boars with genotyped data (50 K BeadChip) from 1365 boars from a nucleus farm were used to perform a genome-wide association study (GWAS) of two breeds, Duroc and Yorkshire. Moreover, we performed a genetic parameter estimation for four FE and FE-related traits. The heritabilities of the FE and FE-related traits ranged from 0.13 to 0.36, and there were significant genetic correlations (−0.69 to 0.52) of the FE and FE-related traits with two growth traits (age at 100 kg and backfat thickness at 100 kg). A total of 61 significant SNPs located on eight different chromosomes associated with the four FE and FE-related traits were identified. We further identified four regions associated with FE and FE-related traits that have not been previously reported, and they may be potential novel QTLs for FE. Considering their biological functions, we finally identified 35 candidate genes relevant for FE and FE-related traits, such as the widely reported MC4R and INSR genes. A gene enrichment analysis showed that FE and FE-related traits were highly enriched in the biosynthesis, digestion, and metabolism of biomolecules. This study deepens our understanding of the genetic mechanisms of FE in pigs and provides valuable information for using marker-assisted selection in pigs to improve FE.

Establishment of a circular RNA regulatory stemness-related gene pair signature for predicting prognosis and therapeutic response in colorectal cancer

Background

Colorectal cancer (CRC) is a common malignant tumor of the digestive tract with a poor prognosis. Cancer stem cells (CSCs) affect disease outcomes and treatment responses in CRC. We developed a circular RNA (circRNA) regulatory stemness-related gene pair (CRSRGP) signature to predict CRC patient prognosis and treatment effects.

Methods

The circRNA, miRNA, and mRNA expression profiles and clinical information of CRC patients were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. CRSRGPs were established based on stemness-related genes in the competing endogenous RNA (ceRNA) network. A CRSRGP signature was generated using the least absolute shrinkage and selection operator (Lasso) and Cox regression analysis of TCGA training set. The prognosis was predicted by generating a nomogram integrating the CRSRGP signature and clinicopathologic features. The model was validated in an external validation set (GSE17536). The antitumor drug sensitivity and immunotherapy responses of CRC patients in the high-risk group (HRG) and low-risk group (LRG) were evaluated by the pRRophetic algorithm and immune checkpoint analysis.

Results

We established an 18-CRSRGP signature to predict the prognosis and treatment responses of CRC patients. In the training and external validation sets, risk scores were used to categorize CRC patients into the HRG and LRG. The Kaplan–Meier analysis showed a poor prognosis for patients in the HRG and that subgroups with different clinical characteristics had significantly different prognoses. A multivariate Cox analysis revealed that the CRSRGP signature was an independent prognostic factor. The nomogram integrating clinical features and the CRSRGP signature efficiently predicted CRC patient prognosis, outperformed the current TNM staging system, and had improved practical clinical value. Anticancer drug sensitivity predictions revealed that the tumors of patients in the HRG were more sensitive to pazopanib, sunitinib, gemcitabine, lapatinib, and cyclopamine. Analysis of immune checkpoint markers demonstrated that patients in the HRG were more likely to benefit from immunotherapy.

Conclusion

An efficient, reliable tool for evaluating CRC patient prognosis and treatment response was established based on the 18-CRSRGP signature and nomogram.

Potential of ATP5MG to Treat Metabolic Syndrome-Associated Cardiovascular Diseases

INTRODUCTION

Metabolic syndrome-associated cardiovascular disease (MetS-CVD) is a cluster of metabolism-immunity highly integrated diseases. Emerging evidence hints that mitochondrial energy metabolism may be involved in MetS-CVD development. The physiopathological role of ATP5MG, a subunit of the F0 ATPase complex, has not been fully elucidated.

METHODS

In this study, we selected ATP5MG to identify the immunity-mediated pathway and mine drugs targeting this pathway for treating MetS-CVD. Using big data from public databases, we dissected co-expressed RNA (coRNA), competing endogenous RNA (ceRNA), and interacting RNA (interRNA) genes for ATP5MG.

RESULTS

It was identified that ATP5MG may form ceRNA with COX5A through hsa-miR-142-5p and interplay with NDUFB8, SOD1, and MDH2 through RNA-RNA interaction under the immune pathway. We dug out 251 chemicals that may target this network and identified some of them as clinical drugs. We proposed five medicines for treating MetS-CVD. Interestingly, six drugs are being tested to treat COVID-19, which unexpectedly offers a new potential host-targeting antiviral strategy.

CONCLUSION

Collectively, we revealed the potential significance of the ATP5MG-centered network for developing drugs to treat MetS-CVD, which offers insights into the epigenetic regulation for metabolism-immunity highly integrated diseases.

Transcriptome Analysis Reveals Molecular Mechanisms under Salt Stress in Leaves of Foxtail Millet (Setaria italica L.)

Foxtail millet (Setaria italica L.) is an important cereal for managing future water scarcity and ensuring food security, due to its strong drought and salt stress resistance owing to its developed root system. However, the molecular responses of foxtail millet leaves to salt stress are largely unknown. In this study, seeds of 104 foxtail millet accessions were subjected to 0.17 mol·L^-1 NaCl stress during germination, and various germination-related parameters were analyzed to derive 5 salt-sensitive accessions and 13 salt-tolerant accessions. Hong Gu 2000 and Pu Huang Yu were the most salt-tolerant and salt-sensitive accessions, respectively. To determine the mechanism of the salt stress response, transcriptomic differences between the control and salt-treated groups were investigated. We obtained 2019 and 736 differentially expressed genes under salt stress in the salt-sensitive and salt-tolerant accessions, respectively. The transcription factor families bHLH, WRKY, AP2/ERF, and MYB-MYC were found to play critical roles in foxtail millet's response to salt stress. Additionally, the down-regulation of ribosomal protein-related genes causes stunted growth in the salt-sensitive accessions. The salt-tolerant accession alleviates salt stress by increasing energy production. Our findings provide novel insights into the molecular mechanism of foxtail millet's response to salt stress.

Downregulated Dual-Specificity Protein Phosphatase 1 in Ovarian Carcinoma: A Comprehensive Study With Multiple Methods

Introduction: We aimed to explore the abnormal expression of dual-specificity protein phosphatase 1 (DUSP1) and its latent molecular mechanisms in ovarian carcinoma (OVCA).

Materials and Methods: Two clinical cohorts collected from two different hospitals were used to evaluate the expression of DUSP1 protein in OVCA tissues. RNA-sequencing and microarray datasets were utilised to verify DUSP1 expression at mRNA levels in both OVCA tissues and in the peripheral blood of OVCA patients. Furthermore, an integrated calculation was performed to pool the standard mean difference (SMD) from each cohort in order to comprehensively assess the expression of DUSP1 in OVCA. Furthermore, we examined the relationship among DUSP1, tumour microenvironment (TME), and chemotherapy resistance in OVCA. Moreover, we used pathway enrichment analysis to explore the underlying mechanisms of DUSP1 in OVCA.

Results: A pooled SMD of −1.19 (95% CI [−2.00, −0.38], p = 0.004) with 1,240 samples revealed that DUSP1 was downregulated in OVCA at both mRNA and protein levels. The area under the receiver operating characteristic curve of 0.9235 indicated the downregulated DUSP1 in peripheral blood may have a non-invasive diagnostic value in OVCA. Through six algorithms, we identified that DUSP1 may related to tumour-infiltrating T cells and cancer associated fibroblasts (CAFs) in OVCA. Pathway enrichment demonstrated that DUSP1 might participate in the mitogen-activated protein kinase (MAPK) signalling pathway. Furthermore, DUSP1 may have relations with chemotherapy resistance, and a favourable combining affinity was observed in the paclitaxel-DUSP1 docking model.

Conclusion: DUSP1 was downregulated in OVCA, and this decreasing trend may affect the infiltration of CAFs. Finally, DUSP1 may have a targeting relation with paclitaxel and participate in MAPK signaling pathways.

Genomic status of yellow-breasted bunting following recent rapid population decline

Global biodiversity is facing serious threats. However, knowledge of the genomic consequences of recent rapid population declines of wild organisms is limited. Do populations experiencing recent rapid population decline have the same genomic status as wild populations that experience long-term declines? Yellow-breasted Bunting (Emberiza aureola) is a critically endangered species that has been experiencing a recent rapid population decline. To answer the question, we assembled and annotated the whole genome of Yellow-breasted Bunting. Furthermore, we found high genetic diversity, low linkage disequilibrium, and low proportion of long runs of homozygosity in Yellow-breasted Bunting, suggesting that the populations following recent rapid declines have different genomic statuses from the population that experienced long-term population decline.

Whole genome resequencing reveals the genetic contribution of Kazakh and Swiss Brown cattle to a population of Xinjiang Brown cattle

Xinjiang Brown cattle is characterized by wide environmental adaptation from its female parent (Kazakh cattle) and good production performance from its male parent (Swiss Brown cattle). However, the genetic basis underlying these characteristics hasnotbeenexplored. Here we compared 50 genomes of Xinjiang Brown cattle to the genomes of other eight breeds worldwide to analyze patterns of genetic variation in the Xinjiang Brown cattle. We found canonical genomic characteristics of cross breed with the lowest linkage disequilibrium and the highest effective population size. At the global level, Xinjiang Brown cattle had 9.88% Kazakh cattle and 90.12% Swiss Brown cattle inheritance. Our local ancestry inference revealed the segments with the excess of Kazakh cattle blood enriched in genes or pathways involved in digestion,absorption,metabolism and disease. More importantly, we also observed the completely fixed haplotypes inherited from Swiss Brown cattle harboring genes (LCORL, GHR, MEF2D, PCSK1 and MSRB3), KEGG pathways, cattle QTLs or human NHGRIGWAS catalog related to body measurement and growth traits. Our findings will not only help understand the process of cross breeding but can provide basic materials for further QTL mapping and improvement of important traits in Xinjiang Brown cattle.

Different functional connectivity optimal frequency in autism compared with healthy controls and the relationship with social communication deficits: Evidence from gene expression and behavior symptom analyses

Studies have reported that different brain regions/connections possess distinct frequency properties, which are related to brain function. Previous studies have proposed altered brain activity frequency and frequency-specific functional connectivity (FC) patterns in autism spectrum disorder (ASD), implying the varied dominant frequency of FC in ASD. However, the difference of the dominant frequency of FC between ASD and healthy controls (HCs) remains unclear. In the present study, the dominant frequency of FC was measured by FC optimal frequency, which was defined as the intermediate of the frequency bin at which the FC strength could reach the maximum. A multivariate pattern analysis was conducted to determine whether the FC optimal frequency in ASD differs from that in HCs. Partial least squares regression (PLSR) and enrichment analyses were conducted to determine the relationship between the FC optimal frequency difference of ASD/HCs and cortical gene expression. PLSR analyses were also performed to explore the relationship between FC optimal frequency and the clinical symptoms of ASD. Results showed a significant difference of FC optimal frequency between ASD and HCs. Some genes whose cortical expression patterns are related to the FC optimal frequency difference of ASD/HCs were enriched for social communication problems. Meanwhile, the FC optimal frequency in ASD was significantly related to social communication symptoms. These results may help us understand the neuro-mechanism of the social communication deficits in ASD.

PtoNF-YC9-SRMT-PtoRD26 module regulates the high saline tolerance of a triploid poplar

BACKGROUND

Sensing and responding to stresses determine the tolerance of plants to adverse environments. The triploid Chinese white poplar is widely cultivated in North China because of its adaptation to a wide range of habitats including highly saline ones. However, its triploid genome complicates any detailed investigation of the molecular mechanisms underlying its adaptations.

RESULTS

We report a haplotype-resolved genome of this triploid poplar and characterize, using reverse genetics and biochemical approaches, a MYB gene, SALT RESPONSIVE MYB TRANSCRIPTION FACTOR (SRMT), which combines NUCLEAR FACTOR Y SUBUNIT C 9 (PtoNF-YC9) and RESPONSIVE TO DESICCATION 26 (PtoRD26), to regulate an ABA-dependent salt-stress response signaling. We reveal that the salt-inducible PtoRD26 is dependent on ABA signaling. We demonstrate that ABA or salt drives PtoNF-YC9 shuttling into the nucleus where it interacts with SRMT, resulting in the rapid expression of PtoRD26 which in turn directly regulates SRMT. This positive feedback loop of SRMT-PtoRD26 can rapidly amplify salt-stress signaling. Interference with either component of this regulatory module reduces the salt tolerance of this triploid poplar.

CONCLUSION

Our findings reveal a novel ABA-dependent salt-responsive mechanism, which is mediated by the PtoNF-YC9-SRMT-PtoRD26 module that confers salt tolerance to this triploid poplar. These genes may therefore also serve as potential and important modification targets in breeding programs.

Role of PsnWRKY70 in Regulatory Network Response to Infection with Alternaria alternata (Fr.) Keissl in Populus

WRKY is an important complex family of transcription factors involved in plant immune responses. Among them, WRKY70 plays an important role in the process of the plant defense response to the invasion of pathogens. However, the defense mechanism of PsnWRKY70 is not clear in Populus nigra. In this study, we showed that PsnWRKY70-overexpression lines (OE) had fewer leaf blight symptoms than PsnWRKY70-repressing lines (RE). PsnWRKY70 activated MAP kinase cascade genes (PsnM2K4, PsnMPK3, PsnM3K18), calcium channel proteins-related genes (PsnCNG3, PsnCNGC1, PsnCNG4), and calcium-dependent protein kinases genes (PsnCDPKL, PsnCDPKW, PsnCDPKS, PsnCDPKQ). Furthermore, 129 genes of PsnWRKY70 putative genome-wide direct targets (DTGs) were identified by using transcriptome (RNA-seq) and DNA affinity purification sequencing (DAP-seq). PsnWRKY70 directly binds to the promoters of homologous genes and LRR domain proteins to promote the expression of WRKY6, WRKY18, WRKY22, and WRKY22–1, LRR domain proteins LRR8, LRR-RLK, ADR1-like 2, NB-ARC, etc. Our study suggests that PsnWRKY70 enhances the resistance of A. alternata in poplar by activating genes in both pathogen-associated molecular pattern-triggered immunity (PTI) and effector-triggered immunity (ETI).

Genomic Diversity and Selection Signatures for Weining Cattle on the Border of Yunnan-Guizhou

Weining cattle is a Chinese indigenous breed influenced by complex breeding and geographical background. The multi-ethnic breeding culture makes Weining cattle require more attention as livestock resources for its genetic diversity. Here, we used 10 Weining cattle (five newly sequenced and five downloaded) and downloaded another 48 genome data to understand the aspects of Weining cattle: genetic diversity, population structure, and cold-adapted performance. In the current study, a high level of genetic diversity was found in Weining cattle, and its breed comprised two potential ancestries, which were Bos taurus and Bos indicus. The positive selective sweep analysis in Weining cattle was analyzed using composite likelihood ratio (CLR) and nucleotide diversity (θπ), resulting in 203 overlapped genes. In addition, we studied the cold adaptation of Weining cattle by comparing with other Chinese cattle (Wannan and Wenshan cattle) by three methods (FST, θπ-ratio, and XP-EHH). Of the top 1% gene list, UBE3D and ZNF668 were analyzed, and these genes may be associated with fat metabolism and blood pressure regulation in cold adaptation. Our findings have provided invaluable information for the development and conservation of cattle genetic resources, especially in southwest China.

Identification of Ferroptotic Genes in Spinal Cord Injury at Different Time Points: Bioinformatics and Experimental Validation

Programmed cell death (PCD) is an important pathologic process after spinal cord injury (SCI). As a new type of PCD, ferroptosis is involved in the secondary SCI. However, the underlying molecular mechanism remains unclear. In this study, we validated ferroptotic phenotype in an animal model of SCI. Then, the bioinformatic analyses performed on a microarray data of SCI (GSE45006). KEGG analysis suggested that the pathways of mTOR, HIF-1, VEGF, and protein process in endoplasmic reticulum were involved in SCI-induced ferroptosis. GO analysis revealed that oxidative stress, amide metabolic process, cation transport, and cytokine production were essential biological processes in ferroptosis after SCI. We highlighted five genes including ATF-3, XBP-1, HMOX-1, DDIT-3, and CHAC-1 as ferroptotic key gene in SCI. These results contribute to exploring the ferroptotic mechanism underlying the secondary SCI and providing potential targets for clinical treatment.

Transcriptomic Profiling Revealed Signaling Pathways Associated with the Spawning of Female Zebrafish under Cold Stress

As one of the critical abiotic factors, temperature controls fish development and reproduction. However, the effects of low temperature on the transcriptional regulation of zebrafish reproduction remain largely unclear. In this study, the fecundity of zebrafish was examined after exposure to cold temperatures at 19.5 °C, 19 °C, 18.5 °C, or 18 °C. The temperature at 19 °C showed no significant influence on the fecundity of zebrafish, but temperature at 18.5 °C or 18 °C significantly blocked the spawning of females, suggesting the existence of a low temperature critical point for the spawning of zebrafish females. Based on these observations, the brains of anesthetized fish under cold stress at different cold temperatures were collected for high-throughput RNA-seq assays. Key genes, hub pathways and important biological processes responding to cold temperatures during the spawning of zebrafish were identified through bioinformatic analysis. The number of down-regulated and up-regulated genes during the temperature reduction from egg-spawning temperatures at 19.5 °C and 19 °C to non-spawning temperatures at 18.5 °C and 18 °C were 2588 and 2527 (fold change ≥ 1.5 and p-value ≤ 0.01), respectively. Venn analysis was performed to identify up- and down-regulated key genes. KEGG enrichment analysis indicated that the hub pathways overrepresented among down-regulated key genes included the GnRH signaling pathway, vascular smooth muscle contraction, C-type lectin receptor signaling pathway, phosphatidylinositol signaling system and insulin signaling pathway. GO enrichment analysis of down-regulated key genes revealed the most important biological processes inhibited under non-spawning temperatures at 18.5 °C and 18 °C were photoreceptor cell outer segment organization, circadian regulation of gene expression and photoreceptor cell maintenance. Furthermore, 99 hormone-related genes were found in the brain tissues of non-spawning and spawning groups, and GnRH signaling pathway and insulin signaling pathway were enriched from down-regulated genes related to hormones at 18.5 °C and 18 °C. Thus, these findings uncovered crucial hormone-related genes and signaling pathways controlling the spawning of female zebrafish under cold stress.

Glia Maturation Factor β as a Novel Independent Prognostic Biomarker and Potential Therapeutic Target of Kidney Renal Clear Cell Carcinoma

Kidney renal clear cell carcinoma (KIRC) has the highest mortality rate and potential for invasion among renal cancers. The diagnosis and treatment of KIRC are becoming challenging because of its diverse pathogenic mechanisms. Glia (GMFB) is a highly conserved growth and differentiation factor for glia cells and neurons, and it is closely associated with neurodegenerative diseases. However, its role in KIRC remains unknown. The present study integrated bioinformatics approaches with suitable meta-analyses to determine the position of GMFB in KIRC. There was a significant decrease in Gmfb expression in KIRC kidneys compared with normal controls. Gmfb expression was negatively associated with pathologic stage, T and M stages, and histologic grade. Univariate and multivariate analyses showed that elevated Gmfb expression was an independent factor for a favorable prognosis. Furthermore, the nomogram verified that Gmfb is a low-risk factor for KIRC. Knockdown of Gmfb in Caki-2 cells increased viability and decreased p21 and p27 levels. Overexpression of Gmfb inhibited Caki-2 cell proliferation, migration, and invasion and decreased mitochondrial membrane potential. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses considering Gmfb co-expressed differentially expressed genes (DEGs) showed that collecting duct acid secretion and mineral absorption ranked were the most important upregulated and downregulated DEGs, respectively. The upregulated hub genes for DEGs were mainly involved in nucleosome assembly, nucleosome organization, and chromatin assembly, and the downregulated hub genes were primarily associated with keratinization. The ratio of tumor-infiltrating immune cells in KIRC tissues was evaluated using CIBERSORTx. The results showed that the Gmfb expression was significantly positively correlated with macrophage M2 cells and mast resting cell infiltration levels and negatively correlated with T follicular helper, T regulatory, and B plasma cell infiltration levels. The former cell types were associated with a beneficial outcome, while the latter had a worse outcome in patients with KIRC. In summary, this study identified GMFB as a novel independent biomarker and therapeutic target for KIRC, and it provides a helpful and distinct individualized treatment strategy for KIRC with a combination of molecular targets and tumor microenvironment.

Epigenetic Alterations in Cryopreserved Human Spermatozoa: Suspected Potential Functional Defects

Background: Gene set enrichment analysis (GSEA) was conducted on raw data, and alternative splicing (AS) events were found after mRNA sequencing of human spermatozoa. In this study, we aimed to compare unknown micro-epigenetics alternations in fresh and cryopreserved spermatozoa to evaluate the effectivity of cryopreservation protocols. Methods: Spermatozoa were divided into three groups: fresh spermatozoa (group 1), cryoprotectant-free vitrified spermatozoa (group 2), and conventionally frozen spermatozoa (group 3). Nine RNA samples (three replicates in each group) were detected and were used for library preparation with an Illumina compatible kit and sequencing by the Illumina platform. Results: Three Gene Ontology (GO) terms were found to be enriched in vitrified spermatozoa compared with fresh spermatozoa: mitochondrial tRNA aminoacylation, ATP-dependent microtubule motor activity, and male meiotic nuclear division. In alternative splicing analysis, a number of unknown AS events were found, including functional gene exon skipping (SE), alternative 5′ splice sites (A5SS), alternative 3′ splice sites (A3SS), mutually exclusive exon (MXE), and retained intron (RI). Conclusions: Cryopreservation of spermatozoa from some patients can agitate epigenetic instability, including increased alternative splicing events and changes in crucial mitochondrial functional activities. For fertilization of oocytes, for such patients, it is recommended to use fresh spermatozoa whenever possible; cryopreservation of sperm is recommended to be used only in uncontested situations.

Potential mechanism of oral baicalin treating psoriasis via suppressing Wnt signaling pathway and inhibiting Th17/IL-17 axis by activating PPARγ

Psoriasis (PSO), an immune-mediated chronic inflammatory skin disease, has seriously affected the quality of patients' life. It is urgent to find effective medicines with lower costs and less side effects. Baicalin (HQG) is the main bioactive substance from Scutellaria baicalensis with effects of anti-inflammation and immunoregulation. Herein, we explored the effect of oral HQG treating PSO and its potential mechanism. Firstly, network pharmacology was used to predict that HQG may act on Estrogen, TNF-alpha (tumor necrosis factor, TNF), interleukin-17 (IL-17) signaling pathways and Th17 cell differentiation, especially the key targets including TNF, Proto-oncogene tyrosine-protein kinase Src, Peroxisome proliferator-activated receptor gamma and Matrix metalloproteinase-9. Imiquimod (IMQ)-induced mice were then used to study the effects of HQG treating PSO. HQG could significantly ameliorate the skin lesions, decrease the level of inflammatory factors and inhibit Th1/Th17 cell differentiation in IMQ-induced mice. Finally, transcriptome analysis of skin lesions integrated with the prediction of network pharmacology further demonstrated that the potential mechanism may be associated with suppressing Wnt signaling pathway and inhibiting Th17/IL-17 axis by activating PPARγ. In conclusion, this study suggested that HQG may be a promising agent for further studies in the search for therapeutic strategies to treat PSO in the future.

Regular combined training and vitamins modulated the apoptosis process in diabetic rats: Bioinformatics analysis of heart failure's differential genes expression network correlated with anti-apoptotic process

The apoptosis process could impose significantly by hyperglycemia. According to in silico language processing and high throughput raw data analysis, we recognized hub molecular mechanisms involved in the pathogenesis of diabetic hearts and suggested a new pharmaceutical approach for declining myocardial programed cell death. Fifty male Sprague-Dawley rats were classified into five groups: healthy rats as control, diabetic rats, diabetic combined resistance/endurance training, diabetic rats which consumed supplementation vitamins E and C, and the combined supplementation and training. Here, we calculated changes in gene expression based on artificial intelligence methods and evaluated gene expression in apoptotic influencing combined training and antioxidants vitamins consumption in heart injured models by streptozotocin via Real-Time PCR. Moreover, we assessed the binding affinity of the 3D structure of small molecules on macromolecule SIRT3 to a new compound pharmaceutical suggesting the decline in cell death program. The computational intelligence surveys revealed that the apoptosis process was a remarkable pathomechanism in the abnormality function of heart tissue in diabetic conditions. Furthermore, we showed that synchronizing antioxidant vitamin consumption and regular combined training could significantly decrease irreversible myocardial cell death in diabetic myocardiopathy. Hence, levels of antiapoptotic mRNA were modified in the combined training/vitamin consumption group compared with other classifications. We found that regular combined exercise and vitamin consumption could reverse the apoptosis process to enhance the survival of cardiac muscle cells in diabetes conditions. PRACTICAL APPLICATIONS: Machine learning and system biology indicated that the apoptosis process is a vital pathomechanism of hyperglycemia-induced heart failure. Sirt3/Fas/Bcl-2/Cycs and Bax, as a critical network of apoptosis, play an essential role in heart failure induced by hyperglycemia. Moreover, Type 2 diabetes and obesity increase the risk of heart failure by increasing high blood sugar levels. We calculated the binding power of the vitamins E and C on SIRT3 protein based on the drug software. In addition, this study assessed that regular combined training and vitamin consumption had an antiapoptotic effect. Also, our data might improve the hyperglycemia state.

Increased Tumor Intrinsic Growth Potential and Decreased Immune Function Orchestrate the Progression of Lung Adenocarcinoma

Background

The overall 5-year survival of lung cancer was reported to be only ~15%, with lung adenocarcinoma (LUAD) as the main pathological subtype. Before developing into invasive stages, LUAD undergoes pre-invasive stages of adenocarcinoma in situ (AIS) and minimally invasive adenocarcinoma (MIA), where surgical resection gives an excellent 5-year survival rate. Given the dramatic decline of prognosis from pre-invasive to invasive stages, a deeper understanding of key molecular changes driving the progression of LUAD is highly needed.

Methods

In this study, we performed whole-exome sequencing and RNA sequencing on surgically resected 24 AIS, 74 MIA, 99 LUAD specimens, and their adjacent paired normal tissues. Survival data were obtained by follow-up after surgery. Key molecular events were found by comparing the gene expression profiles of tumors with different stages. Finally, to measure the level of imbalance between tumor intrinsic growth potential and immune microenvironment, a tumor progressive (TP) index was developed to predict tumor progression and patients’ survival outcome and validated by external datasets.

Results

As tumors progressed to more invasive stages, they acquired higher growth potential, mutational frequency of tumor suppressor genes, somatic copy number alterations, and tumor mutation burden, along with suppressed immune function. To better predict tumor progression and patients’ outcome, TP index were built to measure the imbalance between tumor intrinsic growth potential and immune microenvironment. Patients with a higher TP index had significantly worse recurrence-free survival [Hazard ratio (HR), 10.47; 95% CI, 3.21–34.14; p < 0.0001] and overall survival (OS) [Hazard ratio (HR), 4.83e8; 95% CI, 0–Inf; p = 0.0013]. We used The Cancer Genome Atlas (TCGA)-LUAD dataset for validation and found that patients with a higher TP index had significantly worse OS (HR, 1.10; 95% CI, 0.83–1.45; p = 0.048), demonstrating the prognostic value of the TP index for patients with LUAD.

Conclusions

The imbalance of tumor intrinsic growth potential and immune function orchestrate the progression of LUAD, which can be measured by TP index. Our study provided new insights into predicting survival of patients with LUAD and new target discovery for LUAD through assessing the imbalance between tumor intrinsic growth potential and immune function.

Drought Stress Exacerbates Fungal Colonization and Endodermal Invasion and Dampens Defense Responses to Increase Dry Root Rot in Chickpea

Drought plays a central role in increasing the incidence and severity of dry root rot (DRR) disease in chickpea. This is an economically devastating disease, compromising chickpea yields particularly severely in recent years due to erratic rainfall patterns. Macrophomina phaseolina (formerly Rhizoctonia bataticola) is the causal agent of DRR disease in the chickpea plant. The infection pattern in chickpea roots under well-watered conditions and drought stress are poorly understood at present. This study provides detailed disease symptomatology and the characteristics of DRR fungus at morphological and molecular levels. Using microscopy techniques, the infection pattern of DRR fungus in susceptible chickpea roots was investigated under well-watered and drought-stress conditions. Our observations suggested that drought stress intensifies the progression of already ongoing infection by weakening the endodermal barrier and overall defense. Transcriptomic analysis suggested that the plant's innate immune defense program is downregulated in infected roots when subjected to drought stress. Furthermore, genes involved in hormonal regulation are differentially expressed under drought stress. These findings provide hints in terms of potential chickpea genes to target in crop improvement programs to develop climate-change-resilient cultivars.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

SPI1-related protein inhibits cervical cancer cell progression and prevents macrophage cell migration

AIM

The functions and molecular mechanisms of SPI1-related protein (SPIB) were examined in cervical cancer (CC) cells.

METHODS

Genes related to miscarriage and prognosis in CC were identified by Kaplan-Meier and differential expression analysis, respectively. Cell proliferation, apoptosis, migration, and invasion were examined by cell counting kit-8, flow cytometry, transwell migration, and transwell invasion assays, respectively. The potential functions and molecular mechanisms of SPIB in CC were speculated by gene set enrichment analysis (GSEA) analysis. The mRNA and protein levels of genes were examined by RT-qPCR and western blot assays, respectively. The effect of SPIB on macrophage cells was tested by macrophage recruitment assay and bioinformatics analysis.

RESULTS

A total of 753 dysregulated genes were identified in 88 TCGA CC samples with a history of one or more miscarriages versus 208 CC samples with no miscarriage history. Also, 91 genes related to CC prognosis were identified. SPIB, a gene related to both miscarriage and CC prognosis, inhibited Hela cell proliferation, migration, and invasion, and facilitated Hela cell apoptosis. GSEA analysis disclosed that SPIB might play vital roles in immunity, chemokine signaling pathway, and macrophage chemotaxis/activation in CC. Moreover, SPIB inhibited C-X-C motif chemokine ligand 8 (CXCL8), C-C motif chemokine ligand 17 (CCL17), and C-C motif chemokine ligand 25 (CCL25) expression in Hela cells, and SPIB overexpression in Hela cells hampered THP-1 cell migration. Higher SPIB expression was associated with less M2 macrophage infiltration in CC.

CONCLUSIONS

SPIB inhibited CC-cell progression and hindered macrophage cell migration in CC.

Time Course Analysis of Transcriptome in Human Myometrium Depending on Labor Duration and Correlating With Postpartum Blood Loss

The maintenance of coordinated powerful episodic contractions of the uterus is the crucial factor for normal labor. The uterine contractility is gradually enhanced with the progression of labor, which is related to the gene expression of the myometrium. Competing endogenous RNA (ceRNA) can also regulate the gene expression. To better understand the role of ceRNA network in labor, transcriptome sequencing was performed on the myometrium of 17 parturients at different labor durations (0–24 h). From this, expression levels of mRNA, long non-coding RNA (lncRNA), circular RNA (circRNA), and microRNA (miRNA) were correlated with labor duration. Then, targeting relationships between mRNAs, lncRNAs, circRNAs, and miRNAs were predicted, and the ceRNA regulatory network was established. The mRNA expression patterns associated with cervical dilation and postpartum bleeding were further investigated. This analysis identified 932 RNAs positively correlated with labor duration (859 mRNAs, 28 lncRNAs, and 45 circRNAs) and 153 RNAs negatively correlated with labor duration (122 mRNAs, 28 lncRNAs, and 3 miRNAs). These mRNAs were involved in protein metabolism, transport, and cytoskeleton functions. According to the targeting relationship among these ceRNAs and mRNAs, a ceRNA network consisting of 3 miRNAs, 72 mRNAs, 2 circRNAs, and 1 lncRNA was established. In addition, two mRNA expression patterns were established using time-series analysis of mRNA expression in different phases of cervical dilation. A ceRNA network analysis for blood loss was performed; postpartum bleeding was closely related to inflammatory response, angiogenesis, and hemostasis. This study identified human myometrial transcriptome and established the ceRNA regulatory network depending on labor duration and highlighted the dynamic changes that occur at ceRNAs during parturition, which need to be considered more in the future to better understand how changes in gene expression are relevant to functional changes in human myometrium at labor.

Dissecting the Genetic Structure of Maize Leaf Sheaths at Seedling Stage by Image-Based High-Throughput Phenotypic Acquisition and Characterization

The rapid development of high-throughput phenotypic detection techniques makes it possible to obtain a large number of crop phenotypic information quickly, efficiently, and accurately. Among them, image-based phenotypic acquisition method has been widely used in crop phenotypic identification and characteristic research due to its characteristics of automation, non-invasive, non-destructive and high throughput. In this study, we proposed a method to define and analyze the traits related to leaf sheaths including morphology-related, color-related and biomass-related traits at V6 stage. Next, we analyzed the phenotypic variation of leaf sheaths of 418 maize inbred lines based on 87 leaf sheath-related phenotypic traits. In order to further analyze the mechanism of leaf sheath phenotype formation, 25 key traits (2 biomass-related, 19 morphology-related and 4 color-related traits) with heritability greater than 0.3 were analyzed by genome-wide association studies (GWAS). And 1816 candidate genes of 17 whole plant leaf sheath traits and 1,297 candidate genes of 8 sixth leaf sheath traits were obtained, respectively. Among them, 46 genes with clear functional descriptions were annotated by single nucleotide polymorphism (SNPs) that both Top1 and multi-method validated. Functional enrichment analysis results showed that candidate genes of leaf sheath traits were enriched into multiple pathways related to cellular component assembly and organization, cell proliferation and epidermal cell differentiation, and response to hunger, nutrition and extracellular stimulation. The results presented here are helpful to further understand phenotypic traits of maize leaf sheath and provide a reference for revealing the genetic mechanism of maize leaf sheath phenotype formation.

ARG1 as a promising biomarker for sepsis diagnosis and prognosis: evidence from WGCNA and PPI network

Background

Sepsis is a life-threatening multi-organ dysfunction caused by the dysregulated host response to infection. Sepsis remains a major global concern with high mortality and morbidity, while management of sepsis patients relies heavily on early recognition and rapid stratification. This study aims to identify the crucial genes and biomarkers for sepsis which could guide clinicians to make rapid diagnosis and prognostication.

Methods

Preliminary analysis of multiple global datasets, including 170 samples from patients with sepsis and 110 healthy control samples, revealed common differentially expressed genes (DEGs) in peripheral blood of patients with sepsis. After Gene Oncology (GO) and pathway analysis, the Weighted Gene Correlation Network Analysis (WGCNA) was used to screen for genes most related with clinical diagnosis. Also, the Protein-Protein Interaction Network (PPI Network) was constructed based on the DEGs and the hub genes were found. The results of WGCNA and PPI network were compared and one shared gene was discovered. Then more datasets of 728 experimental samples and 355 control samples were used to prove the diagnostic and prognostic value of this gene. Last, we used real-time PCR to confirm the bioinformatic results.

Results

Four hundred forty-four common differentially expressed genes in the blood of sepsis patients from different ethnicities were identified. Fifteen genes most related with clinical diagnosis were found by WGCNA, and 24 hub genes with most node degrees were identified by PPI network. ARG1 turned out to be the unique overlapped gene. Further analysis using more datasets showed that ARG1 was not only sharply up-regulated in sepsis than in healthy controls, but also significantly high-expressed in septic shock than in non-septic shock, significantly high-expressed in severe or lethal sepsis than in uncomplicated sepsis, and significantly high-expressed in non-responders than in responders upon early treatment. These all demonstrate the performance of ARG1 as a key biomarker. Last, the up-regulation of ARG1 in the blood was confirmed experimentally.

Conclusions

We identified crucial genes that may play significant roles in sepsis by WGCNA and PPI network. ARG1 was the only overlapped gene in both results and could be used to make an accurate diagnosis, discriminate the severity and predict the treatment response of sepsis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41065-022-00240-1.

Next-Generation Sequencing analysis discloses genes implicated in equine endometrosis that may lead to tumorigenesis

Endometrosis is a periglandular fibrosis associated with dysfunction of affected glandular epithelial cells that is the most common cause of reduced fertility in mares, although it is not fully understood. The etiology of the disease is still partially unknown. This study focuses on understanding the genetic mechanisms potentially underlying endometrosis in mares using the Next Generation Sequencing (NGS) technique. Endometrial samples, used in the study, were obtained in the anestrus phase both from healthy mares and those diagnosed with endometrosis. The NGS data were analyzed for gene involvement in biological processes and pathways (e.g. STAR, KOBAS-I, STRING, and ClustVis software). Bioinformatic analysis revealed differential expression of 55 transcripts. In tissues with endometrosis, most genes displayed upregulated expression. The protein-protein interaction analysis disclosed a substantial transcript network including transcripts related to metabolism e.g. sulfur metabolism (SELENBP1), ovarian steroidogenesis, steroid hormone biosynthesis, and chemical carcinogenesis (CYP1B1), COXs (COX4I1, COX3, UQCRFS1) as well as transcripts related to immune response e.g. MMP7, JCHAIN, PIGR, CALR, B2M, FCGRT. Interestingly, the latter has been previously linked with various pathologies including cancers in the female reproductive system. In conclusion, this study evaluated genes that are not directly impacted by sex hormone feedback, but that create a metabolic and immune environment in tissues, thus influencing fertility and pregnancy in mares with endometrosis. Moreover, some of the identified genes may be implicated in tumorigenesis of endometrial lesions. These data may be useful as a starting point in further research, such as the development of targeted strategies for rapid diagnosis and/or prevention of this pathology based on gene and protein-protein interactions.

Gene expression variation explains maize seed germination heterosis

BACKGROUND

Heterosis has been extensively utilized in plant breeding, however, the underlying molecular mechanism remains largely elusive. Maize (Zea mays), which exhibits strong heterosis, is an ideal material for studying heterosis.

RESULTS

In this study, there is faster imbibition and development in reciprocal crossing Zhengdan958 hybrids than in their parent lines during seed germination. To investigate the mechanism of heterosis of maize germination, comparative transcriptomic analyses were conducted. The gene expression patterns showed that 1324 (47.27%) and 1592 (66.44%) of the differential expression genes between hybrids and either parental line display parental dominance up or higher levels in the reciprocal cross of Zhengdan958, respectively. Notably, these genes were mainly enriched in metabolic pathways, including carbon metabolism, glycolysis/gluconeogenesis, protein processing in endoplasmic reticulum, etc. CONCLUSION: Our results provide evidence for the higher expression level genes in hybrid involved in metabolic pathways acting as main contributors to maize seed germinating heterosis. These findings provide new insights into the gene expression variation of maize embryos and improve the understanding of maize seed germination heterosis.

Cross Analysis of Genomic-Pathologic Features on Multiple Primary Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a prevalent malignancy cancer worldwide with a poor prognosis. Hepatic resection is indicated as a potentially curative option for HCC patients in the early stage. However, due to multiple nodules, it leads to clinical challenges for surgical management. Approximately 41%–75% of HCC cases are multifocal at initial diagnosis, which may arise from multicentric occurrence (MO-HCC) or intrahepatic metastasis (IM-HCC) pattern with significantly different clinical outcomes. Effectively differentiating the two mechanisms is crucial to prioritize the allocation of surgery for multifocal HCC. In this study, we collected a multifocal hepatocellular carcinoma cohort of 17 patients with a total of 34 samples. We performed whole-exome sequencing and staining of pathological HE sections for each lesion. Reconstruction of the clonal evolutionary pattern using genome mutations showed that the intrahepatic metastogenesis pattern had a poorer survival performance than independent origins, with variants in the TP53, ARID1A, and higher CNV variants occurring more significantly in the metastatic pattern. Cross-modality analysis with pathology showed that molecular classification results were consistent with pathology results in 70.6% of patients, and we found that pathology results could further complement the classification for undefined patterns of occurrence. Based on these results, we propose a model to differentiate the pattern of multifocal hepatocellular carcinoma based on the pathological results and genome mutations information, which can provide guidelines for diagnosing and treating multifocal hepatocellular carcinoma.

Changes in Higher-Order Chromosomal Structure of Klebsiella pneumoniae Under Simulated Microgravity

Our previous work have shown that certain subpopulations of Klebsiella pneumoniae exhibit significant phenotypic changes under simulated microgravity (SMG), including enhanced biofilm formation and cellulose synthesis, which may be evoked by changes in gene expression patterns. It is well known that prokaryotic cells genomic DNA can be hierarchically organized into different higher-order three-dimensional structures, which can highly influence gene expression. It is remain elusive whether phenotypic changes induced by SMG in the subpopulations of K. pneumoniae are driven by genome higher-order structural changes. Here, we investigated the above-mentioned issue using the wild-type (WT) K. pneumoniae (WT was used as a control strain and continuously cultivated for 2 weeks under standard culture conditions of normal gravity) and two previous identified subpopulations (M1 and M2) obtained after 2 weeks of continuous incubation in a SMG device. By the combination of genome-wide chromosome conformation capture (Hi-C), RNA-seq and whole-genome methylation (WGS) analyses, we found that the along with the global chromosome interactions change, the compacting extent of M1, M2 subpopulations were much looser under SMG and even with an increase in active, open chromosome regions. In addition, transcriptome data showed that most differentially expressed genes (DEGs) were upregulated, whereas a few DEGs were downregulated in M1 and M2. The functions of both types DEGs were mainly associated with membrane fractions. Additionally, WGS analysis revealed that methylation levels were lower in M1 and M2. Using combined analysis of multi-omics data, we discovered that most upregulated DEGs were significantly enriched in the boundary regions of the variable chromosomal interaction domains (CIDs), in which genes regulating biofilm formation were mainly located. These results suggest that K. pneumoniae may regulate gene expression patterns through DNA methylation and changes in genome structure, thus resulting in new phenotypes in response to altered gravity.

Caterpillar-Induced Rice Volatile (E)-β-Farnesene Impairs the Development and Survival of Chilo suppressalis Larvae by Disrupting Insect Hormone Balance

Significant research progress has recently been made on establishing the roles of tps46 in rice defense. (E)-β-farnesene (Eβf) is a major product of tps46 activity but its physiological functions and potential mechanisms against Chilo suppressalis have not yet been clarified. In the present study, C. suppressalis larvae were artificially fed a diet containing 0.8 g/kg Eβf and the physiological performance of the larvae was evaluated. In response to Eβf treatment, the average 2nd instar duration significantly increased from 4.78 d to 6.31 d while that of the 3rd instar significantly increased from 5.70 d to 8.00 d compared with the control. There were no significant differences between the control and Eβf-fed 4th and 5th instars in terms of their durations. The mortalities of the 2nd and 3rd Eβf-fed instars were 21.00-fold and 6.39-fold higher, respectively, than that of the control. A comparative transcriptome analysis revealed that multiple differentially expressed genes are involved in insect hormone biosynthesis. An insect hormone assay on the 3rd instars disclosed that Eβf disrupted the balance between the juvenile hormone and ecdysteroid levels. Eβf treatment increased the juvenile hormones titers but not those of the ecdysteroids. The qPCR results were consistent with those of the RNA-Seq. The foregoing findings suggested that Eβf impairs development and survival in C. suppressalis larvae by disrupting their hormone balance. Moreover, Eβf altered the pathways associated with carbohydrate and xenobiotic metabolism as well as those related to cofactors and vitamins in C. suppressalis larvae. The discoveries of this study may contribute to the development and implementation of an integrated control system for C. suppressalis infestations in rice.

10-Hydroxydec-2-Enoic Acid Reduces Hydroxyl Free Radical-Induced Damage to Vascular Smooth Muscle Cells by Rescuing Protein and Energy Metabolism

10-Hydroxydec-2-enoic acid (10-HDA), an unsaturated hydroxyl fatty acid from the natural food royal jelly, can protect against cell and tissue damage, yet the underlying mechanisms are still unexplored. We hypothesized that the neutralization of the hydroxyl free radical (•OH), the most reactive oxygen species, is an important factor underlying the cytoprotective effect of 10-HDA. In this study, we found that the •OH scavenging rate by 10-HDA (2%, g/ml) was more than 20%, which was achieved through multiple-step oxidization of the -OH group and C=C bond of 10-HDA. Moreover, 10-HDA significantly enhanced the viability of vascular smooth muscle cells (VSMCs) damaged by •OH (P < 0.01), significantly attenuated •OH-derived malondialdehyde production that represents cellular lipid peroxidation (P < 0.05), and significantly increased the glutathione levels in •OH-stressed VSMCs (P < 0.05), indicating the role of 10-HDA in reducing •OH-induced cytotoxicity. Further proteomic analyses of VSMCs identified 195 proteins with decreased expression by •OH challenge that were upregulated by 10-HDA rescue and were primarily involved in protein synthesis (such as translation, protein transport, ribosome, and RNA binding) and energy metabolism (such as fatty acid degradation and glycolysis/gluconeogenesis). Taken together, these findings indicate that 10-HDA can effectively promote cell survival by antagonizing •OH-induced injury in VSMCs. To the best of our knowledge, our results provide the first concrete evidence that 10-HDA-scavenged •OH could be a potential pharmacological application for maintaining vascular health.

Network-Based Analysis Reveals Novel Biomarkers in Peripheral Blood of Patients With Preeclampsia

WGCNA is a potent systems biology approach that explains the connection of gene expression based on a microarray database, which facilitates the discovery of disease therapy targets or potential biomarkers. Preeclampsia is a kind of pregnancy-induced hypertension caused by complex factors. The disease’s pathophysiology, however, remains unknown. The focus of this research is to utilize WGCNA to identify susceptible modules and genes in the peripheral blood of preeclampsia patients. Obtain the whole gene expression data of GSE48424 preeclampsia patients and normal pregnant women from NCBI’s GEO database. WGCNA is used to construct a gene co-expression network by calculating correlation coefficients between modules and phenotypic traits, screening important modules, and filtering central genes. To identify hub genes, we performed functional enrichment analysis, pathway analysis, and protein-protein interaction (PPI) network construction on key genes in critical modules. Then, the genetic data file GSE149437 and clinical peripheral blood samples were used as a validation cohort to determine the diagnostic value of these key genes. Nine gene co-expression modules were constructed through WGCNA analysis. Among them, the blue module is significantly related to preeclampsia and is related to its clinical severity. Thirty genes have been discovered by using the intersection of the genes in the blue module and the DEGs genes as the hub genes. It was found that HDC, MS4A2, and SLC18A2 scored higher in the PPI network and were identified as hub genes. These three genes were also differentially expressed in peripheral blood validation samples. Based on the above three genes, we established the prediction model of peripheral blood markers of preeclampsia and drew the nomogram and calibration curve. The ROC curves were used in the training cohort GSE48424 and the validation cohort GSE149437 to verify the predictive value of the above model. Finally, it was confirmed in the collected clinical peripheral blood samples that MS4A2 was differentially expressed in the peripheral blood of early-onset and late-onset preeclampsia, which is of great significance. This study provides a new biomarker and prediction model for preeclampsia.

An integrative pan-cancer analysis of molecular characteristics and oncogenic role of mitochondrial creatine kinase 1A (CKMT1A) in human tumors

In recent years, several studies have suggested that mitochondrial creatine kinase 1A (CKMT1A) plays a key role in various cancer types. However, there is still a lack of systematic understanding of the contribution of CKMT1A in different types of cancer. Therefore, this study aims to explore the potential role of CKMT1A in human tumors. Firstly, we evaluated the expression level of CKMT1A in 33 types of tumors. Secondly, we used the GEPIA2 and Kaplan–Meier plotter to explore the relationship between CKMT1A expression and survival prognosis. Furthermore, the genetic alterations of CKMT1A were analyzed by the cBioPortal web. In addition, we performed immune infiltration analysis and gene enrichment pathway analysis. CKMT1A was highly expressed in most types of cancers and there was a significant correlation between CKMT1A expression and the prognosis of patients for certain tumors. Non-Small Cell Lung Cancer cases with altered CKMT1A showed a poorer overall survival. CKMT1A expression was negatively correlated with the infiltration of cancer-associated fibroblasts in most tumors. We also found that its expression was negatively associated with CD8⁺ T-cell infiltration in several tumors. Furthermore, enrichment analysis revealed that “Glycolysis/ Gluconeogenesis” and “metabolic pathways” functions were involved in the functional mechanism of CKMT1A. Taken together, our studies will provide a relatively clear and integrative understanding of the role of CKMT1A across different tumors. All these findings will lay a solid foundation for further molecular assays of CKMT1A in tumorigenesis and provide the rationale for developing novel therapeutic strategies.

TRPV1 SUMOylation suppresses itch by inhibiting TRPV1 interaction with H1 receptors

The molecular mechanism underlying the functional interaction between H1R and TRPV1 remains unclear. We show here that H1R directly binds to the carboxy-terminal region of TRPV1 at residues 715-725 and 736-749. Cell-penetrating peptides containing these sequences suppress histamine-induced scratching behavior in a cheek injection model. The H1R-TRPV1 binding is kept at a minimum at rest in mouse trigeminal neurons due to TRPV1 SUMOylation and it is enhanced upon histamine treatment through a transient TRPV1 deSUMOylation. The knockin of the SUMOylation-deficient TRPV1^K823R mutant in mice leads to constitutive enhancement of H1R-TRPV1 binding, which exacerbates scratching behaviors induced by histamine. Conversely, SENP1 conditional knockout in sensory neurons enhances TRPV1 SUMOylation and suppresses the histamine-induced scratching response. In addition to interfering with binding, TRPV1 SUMOylation promotes H1R degradation through ubiquitination. Our work unveils the molecular mechanism of histaminergic itch by which H1R directly binds to deSUMOylated TRPV1 to facilitate the transduction of the pruritogen signal to the scratching response.

Genome-Wide Association Study of Feed Efficiency Related Traits in Ducks

Simple Summary

Good feed efficiency (FE) is an important trait to ensure the economic output of the livestock and poultry industries. Herein, a genome-wide association study was conducted to identify potential variants and genes associated with seven FE measures in ducks. Genomic DNA samples of 308 ducks were collected and sequenced. All animals were evaluated concerning body weight gain (BWG), feed intake (FI), residual feed intake (RFI), feed conversion ratio (FCR), and weight at 21 (BW21) and 42 days of age (BW42). Overall, 4 (FCR), 3 (FI), 36 (RFI), 6 (BWG), 8 (BW21), and 10 (BW42) single nucleotide polymorphisms (SNPs) were significantly associated with these FE traits, respectively. Moreover, candidate genes close to the identified variants were found to be mainly involved in key pathways and terms related to metabolism. In summary, these findings improve our understanding of poultry genetics and provide new foundations for breeding programs aimed at maximizing the economic potential of duck breeding and farming.

Abstract

Feed efficiency (FE) is the most important economic trait in the poultry and livestock industry. Thus, genetic improvement of FE may result in a considerable reduction of the cost and energy burdens. As genome-wide association studies (GWASs) can help identify candidate variants influencing FE, the present study aimed to analyze the phenotypic correlation and identify candidate variants of the seven FE traits in ducks. All traits were found to have significant positive correlations with varying degrees. In particular, residual feed intake presented correlation coefficients of 0.61, 0.54, and 0.13 with feed conversion ratio, and feed intake, respectively. Furthermore, data from seven FE-related GWAS revealed 4 (FCR), 3 (FI), 36 (RFI), 6 (BWG), 8 (BW21), and 10 (BW42) SNPs were significantly associated with body weight gain, feed intake, residual feed intake, feed conversion ratio, and weight at 21 and 42 days, respectively. Candidate SNPs of seven FE trait-related genes were involved in galactose metabolism, starch, propanoate metabolism, sucrose metabolism and etc. Taken together, these findings provide insight into the genetic mechanisms and genes involved in FE-related traits in ducks. However, further investigations are warranted to further validate these findings.

Identification of Potential Biomarkers for Pan-Cancer Diagnosis and Prognosis Through the Integration of Large-Scale Transcriptomic Data

Cancer is one of the leading causes of death worldwide, bringing a significant burden to human health and society. Accurate cancer diagnosis and biomarkers that can be used as robust therapeutic targets are of great importance as they facilitate early and effective therapies. Shared etiology among cancers suggests the existence of pan-cancer biomarkers, performance of which could benefit from the large sample size and the heterogeneity of the studied patients. In this study, we conducted a systematic RNA-seq study of 9,213 tumors and 723 para-cancerous tissue samples of 28 solid tumors from the Cancer Genome Atlas (TCGA) database, and 7,008 normal tissue samples from the Genotype-Tissue Expression (GTEx) database. By differential gene expression analysis, we identified 214 up-regulated and 186 downregulated differentially expressed genes (DEGs) in more than 80% of the studied tumors, respectively, and obtained 20 highly linked up- and downregulated hub genes from them. These markers have rarely been reported in multiple tumors simultaneously. We further constructed pan-cancer diagnostic models to classify tumors and para-cancerous tissues using 10 up-regulated hub genes with an AUC of 0.894. Survival analysis revealed that these hub genes were significantly associated with the overall survival of cancer patients. In addition, drug sensitivity predictions for these hub genes in a variety of tumors obtained several broad-spectrum anti-cancer drugs targeting pan-cancer. Furthermore, we predicted immunotherapy sensitivity for cancers based on tumor mutational burden (TMB) and the expression of immune checkpoint genes (ICGs), providing a theoretical basis for the treatment of tumors. In summary, we identified a set of biomarkers that were differentially expressed in multiple types of cancers, and these biomarkers can be potentially used for diagnosis and used as therapeutic targets.

Identification of novel regulatory pathways across normal human bronchial epithelial cell lines (NHBEs) and peripheral blood mononuclear cell lines (PBMCs) in COVID-19 patients using transcriptome analysis

The SARS-CoV-2 is one of the most infectious and deadly coronaviruses, which has gripped the world, causing the COVID-19 pandemic. Despite the numerous studies being conducted on this virus, many uncertainties are with the disease. This is exacerbated by the speedy mutations acquired by the viral strain, which enables the disease to present itself differently in different people, introducing new factors of uncertainty. This study aims at the identification of regulatory pathways across two cell lines, namely, the peripheral blood mononuclear cell line (PBMC) and the normal human bronchial epithelial (NHBE) cell line. Both the above-mentioned cell lines were considered because they support viral replication. Furthermore, the NHBE cell line captures vital changes in the lungs, which are the main organs affected by the COVID-19 patients, and the PBMC cell line is closely linked to the body's immune system. RNA-Seq analysis, differential gene expression and gene set enrichment analysis for pathway identification were followed. Pathway analysis throws light upon the various systems affected in the body due to the COVID-19. Gene regulatory networks associated with the significant pathways were also designed. These networks aid in identifying various gene targets, along with their interactions. Studying the functionality of the pathways and the gene interactions associated with them, aided by long COVID studies, will provide immense clarity about the current COVID-19 scenario. In the long term, this will help in the design of therapeutic approaches against the SARS-CoV-2 and can also contribute to drug repurposing studies. Ultimately, this study identifies and analyses the relationship of various undiscovered or lesser explored pathways in the human body to the SARS-CoV-2 and establish a clearer picture of the association to help streamline further studies and approaches.

Genome-Wide Identification and Characterization of Long Non-Coding RNAs in Longissimus dorsi Skeletal Muscle of Shandong Black Cattle and Luxi Cattle

There is an increasing understanding of the possible regulatory role of long non-coding RNAs (LncRNA). Studies on livestock have mainly focused on the regulation of cell differentiation, fat synthesis, and embryonic development. However, there has been little study of skeletal muscle of domestic animals and the potential role of lncRNA. In this study, the transcriptome numbers of longissimus muscle of different beef cattle (Shandong black catle and Luxi catlle) were used to construct muscle related lncRNAs-miRNA-mRNA interaction network through bioinformatics analysis. This is helpful to clarify the molecular mechanism of bovine muscle development, and can be used to promote animal husbandry and improve animal husbandry production. According to the screening criteria of |FC|≧2 and q < 0.05, a total of 1,415 transcripts (of which 480 were LncRNAs) were differentially expressed (q < 0.05) in the different breeds. Further, we found that the most differentially expressed LncRNAs were found on chromosome 9, in which the differentially expressed LncRNAs targeted 1,164 protein coding genes (MYORG, Wnt4, PAK1, ADCY7,etc) (upstream and downstream<50 Kb). In addition, Pearson’s correlation coefficients of co-expression levels indicated a potential trans regulatory relationship between the differentially expressed LncRNAs and 43844 mRNAs (r > 0.9). The identified co-expressed mRNAs (MYORG, Dll1, EFNB2, SOX6, MYOCD, and MYLK3) are related to the formation of muscle structure, and enriched in muscle system process, strained muscle cell differentiation, muscle cell development, striated muscle tissue development, calcium signaling, and AMPK signaling. Additionally, we also found that some LncRNAs (LOC112444238, LOC101903367, LOC104975788, LOC112441863, LOC112449549, and LOC101907194) may interact with miRNAs related to cattle muscle growth and development. Based on this, we constructed a LncRNAs-miRNA-mRNA interaction network as the putative basis for biological regulation in cattle skeletal muscle. Interestingly, a candidate differential LncRNA (LOC104975788) and a protein-coding gene (Pax7) contain miR-133a binding sites and binding was confirmed by luciferase reporter assay. LOC104975788 may combined miR-133a competitively with Pax7, thus relieving the inhibitory effect of miR-133a on Pax7 to regulate skeletal muscle development. These results will provide the theoretical basis for further study of LncRNA regulation and activity in different cattle breeds.

L-5-hydroxytryptophan promotes antitumor immunity by inhibiting PD-L1 inducible expression

BACKGROUND

The repression or downregulation of programmed death-ligand 1 (PD-L1) can release its inhibition of T cells and activate antitumor immune responses. Although PD-1 and PD-L1 antibodies are promising treatments for diverse tumor types, their inherent disadvantages and immune-related adverse events remain significant issues. The development of small molecule inhibitors targeting the interaction surface of PD-1 and PD-L1 has been reviving, yet many challenges remain. To address these issues, we aimed to find small molecules with durable efficacy and favorable biosafety that alter PD-L1 surface expression and can be developed into a promising alternative and complementary therapy for existing anti-PD-1/PD-L1 therapies.

METHODS

Cell-based screen of 200 metabolic molecules using a high-throughput flow cytometry assay of PD-L1 surface expression was conducted, and L-5-hydroxytryptophan (L-5-HTP) was found to suppress PD-L1 expression induced by interferon gamma (IFN-γ). Inhibition of PD-L1 induction and antitumor effect of L-5-HTP were evaluated in two syngeneic mouse tumor models. Flow cytometry was performed to investigate the change in the tumor microenvironment caused by L-5-HTP treatment.

RESULTS

We discovered that L-5-HTP suppressed IFN-γ-induced PD-L1 expression in tumor cells transcriptionally, and this effect was directly due to itself. Mechanistically, L-5-HTP inhibited IFN-γ-induced expression of RTK ligands and thus suppressed phosphorylation-mediated activation of RTK receptors and the downstream MEK/ERK/c-JUN signaling cascade, leading to decreased PD-L1 induction. In syngeneic mouse tumor models, treatment with 100 mg/kg L-5-HTP (intraperitoneal) inhibited PD-L1 expression and exhibited antitumor effect. L-5-HTP upregulated the ratio of granzyme B+ CD8+ activated cytotoxic T cells. An intact immune system and PD-L1 expression was critical for L-5-HTP to exert its antitumor effects. Furthermore, L-5-HTP acted synergistically with PD-1 antibody to improve anticancer effect.

CONCLUSION

Our study illustrated L-5-HTP's inhibitory effect on PD-L1 induction stimulated by IFN-γ in tumor cells and also provided insight into repurposing L-5-HTP for use in tumor immunotherapy.

Identification of ST1 reveals a selection involving hitchhiking of seed morphology and oil content during soybean domestication

Seed morphology and quality of cultivated soybean (Glycine max) have changed dramatically during domestication from their wild relatives, but their relationship to selection is poorly understood. Here, we describe a semi-dominant locus, ST1 (Seed Thickness 1), affecting seed thickness and encoding a UDP-D-glucuronate 4-epimerase, which catalyses UDP-galacturonic acid production and promotes pectin biosynthesis. Interestingly, this morphological change concurrently boosted seed oil content, which, along with up-regulation of glycolysis biosynthesis modulated by ST1, enabled soybean to become a staple oil crop. Strikingly, ST1 and an inversion controlling seed coat colour formed part of a single selective sweep. Structural variation analysis of the region surrounding ST1 shows that the critical mutation in ST1 existed in earlier wild relatives of soybean and the region containing ST1 subsequently underwent an inversion, which was followed by successive selection for both traits through hitchhiking during selection for seed coat colour. Together, these results provide direct evidence that simultaneously variation for seed morphology and quality occurred earlier than variation for seed coat colour during soybean domestication. The identification of ST1 thus sheds light on a crucial phase of human empirical selection in soybeans and provides evidence that our ancestors improved soybean based on taste.

ABCA8 inhibits breast cancer cell proliferation by regulating the AMP activated protein kinase/mammalian target of rapamycin signaling pathway

ATP-binding cassette (ABC) subfamily A member 8 (ABCA8) has been reported to play a vital role in cancer development. Our study aimed to explore the role and the molecular mechanism of ABCA8 in breast cancer (BC) progression. GSE65194, GSE15852, and GSE45827 datasets were used to identify differentially expressed genes (DEGs) in BC. The diagnosis and prognosis value were determined using ROC curve analysis and Kaplan-Meier plotter, respectively. The relationship between ABCA8 expression and clinicopathological features in BC was analyzed by TCGA. Co-expressed genes of ABCA8 in BC were screened out through GEPIA and subjected to KEGG pathway enrichment analysis. Cell proliferation was evaluated by CCK-8 and EdU incorporation assays. Proliferating cell nuclear antigen (PCNA) expression and the changes of the AMP activated protein kinase (AMPK)/mammalian target of rapamycin (mTOR) pathway were measured by western blot analysis. Totally 4 overlapping DEGs were identified and all reduced in BC samples. ABCA8 with high diagnostic and prognostic values was selected for further exploration. Low ABCA8 expression was correlated with clinicopathological features in BC patients. ABCA8 overexpression inhibited BC cell proliferation. The top 20 co-expressed genes of ABCA8 were identified by GEPIA and significantly enriched in AMPK signaling pathway. Inhibition of AMPK/mTOR pathway reversed the suppressive effect of ABCA8 on BC cell growth. These results suggested that ABCA8 overexpression repressed BC cell proliferation through regulating the AMPK/mTOR signaling pathway.