Identification of overlay differentially expressed genes in both rats and goats with blast lung injury through comparative transcriptomics / 中华创伤杂志(英文版)

PURPOSE@#To identify the potential target genes of blast lung injury (BLI) for the diagnosis and treatment.@*METHODS@#This is an experimental study. The BLI models in rats and goats were established by conducting a fuel-air explosive power test in an unobstructed environment, which was subsequently validated through hematoxylin-eosin staining. Transcriptome sequencing was performed on lung tissues from both goats and rats. Differentially expressed genes were identified using the criteria of q ≤ 0.05 and |log2 fold change| ≥ 1. Following that, enrichment analyses were conducted for gene ontology and the Kyoto Encyclopedia of Genes and Genomes pathways. The potential target genes were further confirmed through quantitative real-time polymerase chain reaction and enzyme linked immunosorbent assay.@*RESULTS@#Observations through microscopy unveiled the presence of reddish edema fluid, erythrocytes, and instances of focal or patchy bleeding within the alveolar cavity. Transcriptome sequencing analysis identified a total of 83 differentially expressed genes in both rats and goats. Notably, 49 genes exhibited a consistent expression pattern, with 38 genes displaying up-regulation and 11 genes demonstrating down-regulation. Enrichment analysis highlighted the potential involvement of the interleukin-17 signaling pathway and vascular smooth muscle contraction pathway in the underlying mechanism of BLI. Furthermore, the experimental findings in both goats and rats demonstrated a strong association between BLI and several key genes, including anterior gradient 2, ankyrin repeat domain 65, bactericidal/permeability-increasing fold containing family A member 1, bactericidal/permeability-increasing fold containing family B member 1, and keratin 4, which exhibited up-regulation.@*CONCLUSIONS@#Anterior gradient 2, ankyrin repeat domain 65, bactericidal/permeability-increasing fold containing family A member 1, bactericidal/permeability-increasing fold containing family B member 1, and keratin 4 hold potential as target genes for the prognosis, diagnosis, and treatment of BLI.

Spatial transcriptomics reveals that metabolic characteristics define the tumor immunosuppression microenvironment via iCAF transformation in oral squamous cell carcinoma / 国际口腔科学杂志·英文版

Tumor progression is closely related to tumor tissue metabolism and reshaping of the microenvironment. Oral squamous cell carcinoma (OSCC), a representative hypoxic tumor, has a heterogeneous internal metabolic environment. To clarify the relationship between different metabolic regions and the tumor immune microenvironment (TME) in OSCC, Single cell (SC) and spatial transcriptomics (ST) sequencing of OSCC tissues were performed. The proportion of TME in the ST data was obtained through SPOTlight deconvolution using SC and GSE103322 data. The metabolic activity of each spot was calculated using scMetabolism, and k-means clustering was used to classify all spots into hyper-, normal-, or hypometabolic regions. CD4T cell infiltration and TGF-β expression is higher in the hypermetabolic regions than in the others. Through CellPhoneDB and NicheNet cell-cell communication analysis, it was found that in the hypermetabolic region, fibroblasts can utilize the lactate produced by glycolysis of epithelial cells to transform into inflammatory cancer-associated fibroblasts (iCAFs), and the increased expression of HIF1A in iCAFs promotes the transcriptional expression of CXCL12. The secretion of CXCL12 recruits regulatory T cells (Tregs), leading to Treg infiltration and increased TGF-β secretion in the microenvironment and promotes the formation of a tumor immunosuppressive microenvironment. This study delineates the coordinate work axis of epithelial cells-iCAFs-Tregs in OSCC using SC, ST and TCGA bulk data, and highlights potential targets for therapy.

Single-Cell Mapping of Brain Myeloid Cell Subsets Reveals Key Transcriptomic Changes Favoring Neuroplasticity after Ischemic Stroke / 神经科学通报·英文版

Interactions between brain-resident and peripheral infiltrated immune cells are thought to contribute to neuroplasticity after cerebral ischemia. However, conventional bulk sequencing makes it challenging to depict this complex immune network. Using single-cell RNA sequencing, we mapped compositional and transcriptional features of peri-infarct immune cells. Microglia were the predominant cell type in the peri-infarct region, displaying a more diverse activation pattern than the typical pro- and anti-inflammatory state, with axon tract-associated microglia (ATMs) being associated with neuronal regeneration. Trajectory inference suggested that infiltrated monocyte-derived macrophages (MDMs) exhibited a gradual fate trajectory transition to activated MDMs. Inter-cellular crosstalk between MDMs and microglia orchestrated anti-inflammatory and repair-promoting microglia phenotypes and promoted post-stroke neurogenesis, with SOX2 and related Akt/CREB signaling as the underlying mechanisms. This description of the brain's immune landscape and its relationship with neurogenesis provides new insight into promoting neural repair by regulating neuroinflammatory responses.

Distinct mononuclear diploid cardiac subpopulation with minimal cell-cell communications persists in embryonic and adult mammalian heart / 医学前沿

A small proportion of mononuclear diploid cardiomyocytes (MNDCMs), with regeneration potential, could persist in adult mammalian heart. However, the heterogeneity of MNDCMs and changes during development remains to be illuminated. To this end, 12 645 cardiac cells were generated from embryonic day 17.5 and postnatal days 2 and 8 mice by single-cell RNA sequencing. Three cardiac developmental paths were identified: two switching to cardiomyocytes (CM) maturation with close CM-fibroblast (FB) communications and one maintaining MNDCM status with least CM-FB communications. Proliferative MNDCMs having interactions with macrophages and non-proliferative MNDCMs (non-pMNDCMs) with minimal cell-cell communications were identified in the third path. The non-pMNDCMs possessed distinct properties: the lowest mitochondrial metabolisms, the highest glycolysis, and high expression of Myl4 and Tnni1. Single-nucleus RNA sequencing and immunohistochemical staining further proved that the Myl4+Tnni1+ MNDCMs persisted in embryonic and adult hearts. These MNDCMs were mapped to the heart by integrating the spatial and single-cell transcriptomic data. In conclusion, a novel non-pMNDCM subpopulation with minimal cell-cell communications was unveiled, highlighting the importance of microenvironment contribution to CM fate during maturation. These findings could improve the understanding of MNDCM heterogeneity and cardiac development, thus providing new clues for approaches to effective cardiac regeneration.

New strategies for the treatment of carcinoma of unknown primary / 中华肿瘤杂志

Carcinoma of unknown primary (CUP) is a kind of metastatic tumor whose primary origin cannot be identified after adequate examination and evaluation. The main treatment modality of CUP is empiric chemotherapy, and the median overall survival time is less than 1 year. Compared with immunohistochemistry, novel method based on gene expression profiling have improved the sensitivity and specificity of CUP detection, but its guiding value for treatment is still controversial. The approval of immune checkpoint inhibitors and pan-cancer antitumor agents has improved the prognosis of patients with CUP, and targeted therapy and immunotherapy based on specific molecular characteristics are the main directions of future research. Given the high heterogeneity and unique clinicopathological characteristics of CUP, "basket trial" is more suitable for clinical trial design in CUP.

STAT3 as a candidate transcriptomic prognosticator of sepsis severity levels

Background@#Sepsis is a life-threatening multiple-organ dysfunction caused by a dysregulated host response to infection and is the leading cause of death in non-cardiac intensive care facilities. Early reliable prediction of sepsis outcomes leads to cost-efficient resource allocation and therapeutic strategies. However, there are still no reliable markers to predict the outcome of patients at the initial stage of sepsis. Analyzing transcription profiles enables researchers to predict early outcomes using transcripts and their expression patterns. Transcriptomic profiling of septic patients has been done recently; however, analysis of prognostic outcomes is still scarce.@*Objective@#This study aimed to determine transcriptional indicators that may be useful in the prognosis of the severity of sepsis.@*Methods@#This is a prospective cohort study of Filipino patients admitted for sepsis at the national tertiary referral hospital in Manila, Philippines. We conducted differentially expressed gene analysis, network analyses, and area under the curve study of publicly available datasets of surviving vs. non-surviving sepsis patients to identify candidate prognosticator markers. Quantitative PCR was used to characterize the expression of each marker. A model using ordinal logistic regression analysis was done to determine which among the markers can best predict the outcome of sepsis severity.@*Results@#We identified ACTB, RAC1, STAT3, and UBQLN1 as candidate mRNA prognosticators. The expression of STAT3, a gene involved in immunosuppression, is inversely correlated with the severity of sepsis.@*Conclusion@#Transcriptomic markers such as STAT3 can predict the severity of patients with sepsis. Early detection of its inverse expression may prompt early and more aggressive management of patients.

Identification and validation of novel biomarkers for cold-dampness syndrome of rheumatoid arthritis based on integration of multiple bioinformatics methods / 中国中药杂志

This study aims to identify the novel biomarkers of cold-dampness syndrome(RA-Cold) of rheumatoid arthritis(RA) by gene set enrichment analysis(GSEA), weighted gene correlation network analysis(WGCNA), and clinical validation. Firstly, transcriptome sequencing was carried out for the whole blood samples from RA-Cold patients, RA patients with other traditional Chinese medicine(TCM) syndromes, and healthy volunteers. The differentially expressed gene(DEG) sets of RA-Cold were screened by comparison with the RA patients with other TCM syndromes and healthy volunteers. Then, GSEA and WGCNA were carried out to screen the key DEGs as candidate biomarkers for RA-Cold. Experimentally, the expression levels of the candidate biomarkers were determined by RT-qPCR for an independent clinical cohort(not less than 10 cases/group), and the clinical efficacy of the candidates was assessed using the receiver operating characteristic(ROC) curve. The results showed that 3 601 DEGs associated with RA-Cold were obtained, including 106 up-regulated genes and 3 495 down-regulated genes. The DEGs of RA-Cold were mainly enriched in the pathways associated with inflammation-immunity regulation, hormone regulation, substance and energy metabolism, cell function regulation, and synovial pannus formation. GSEA and WGCNA showed that recombinant proteasome 26S subunit, ATPase 2(PSMC2), which ranked in the top 50% in terms of coefficient of variation, representativeness of pathway, and biological modules, was a candidate biomarker of RA-Cold. Furthermore, the validation results based on the clinical independent sample set showed that the F1 value, specificity, accuracy, and precision of PSMC2 for RA-Cold were 70.3%, 61.9%, 64.5%, and 81.3%, respectively, and the area under the curve(AUC) value was 0.96. In summary, this study employed the "GSEA-WGCNA-validation" integrated strategy to identify novel biomarkers of RA-Cold, which helped to improve the TCM clinical diagnosis and treatment of core syndromes in RA and provided an experimental basis for TCM syndrome differentiation.

Screening and expression analysis of transcription factors involved in genuineness of Codonopsis pilosula in Shanxi / 中国中药杂志

This study aims to mine the transcription factors that affect the genuineness of Codonopsis pilosula in Shanxi based on the transcriptome data of C. pilosula samples collected from Shanxi and Gansu, and then analyze the gene expression patterns, which will provide a theoretical basis for the molecular assisted breeding of C. pilosula. Gene ontology(GO) functional annotation, conserved motif prediction, and gene expression pattern analysis were performed for the differential transcription factors predicted based on the transcriptome data of C. pilosula from different habitats. A total of 61 differentially expressed genes(DEGs) were screened out from the transcriptome data. Most of the DEGs belonged to AP2/ERF-ERF family, with the conserved motif of [2X]-[LG]-[3X]-T-[3X]-[AARAYDRAA]-[3X]-[RG]-[2X]-A-[2X]-[NFP]. Forty-three of the DEGs showed significantly higher gene expression in C. pilosula samples from Shanxi than in the samples from Gansu, including 11 genes in the AP2/ERF-ERF family, 5 genes in the NAC fa-mily, 1 gene in the bHLH family, and 2 genes in the RWP-RK family, while 18 transcription factors showed higher expression levels in the samples from Gansu. GO annotation predicted that most of the DEGs were enriched in GO terms related to transcriptional binding activity(103), metabolic process(26), and stress response(23). The expression of transcription factor genes, CpNAC92, CpNAC100, CpbHLH128, and CpRAP2-7 was higher in the samples from Shanxi and in the roots of C. pilosula. CpNAC92, CpbHLH128, and CpRAP2-7 responded to the low temperature, temperature difference, and iron stresses, while CpNAC100 only responded to low temperature and iron stresses. The screening and expression analysis of the specific transcription factors CpNAC92, CpNAC100, CpbHLH128, and CpRAP2-7 in C. pilosula in Shanxi laid a theoretical foundation for further research on the mechanism of genuineness formation of C. pilosula.

Transcriptional regulation mechanism of differential accumulation of flavonoids in leaves and roots of Sarcandra glabra based on metabonomics and transcriptomics / 中国中药杂志

This study aims to explore the molecular regulation mechanism of the differential accumulation of flavonoids in the leaves and roots of Sarcandra glabra. Liquid chromatography-mass spectrometry(LC-MS) and high-throughput transcriptome sequencing(RNA-seq) were employed to screen out the flavonoid-related differential metabolites and differentially expressed genes(DEGs) encoding key metabolic enzymes. Eight DEGs were randomly selected for qRT-PCR verification. The results showed that a total of 37 flavonoid-related differential metabolites between the leaves and roots of S. glabra were obtained, including pinocembrin, phlorizin, na-ringenin, kaempferol, leucocyanidin, and 5-O-caffeoylshikimic acid. The transcriptome analysis predicted 36 DEGs associated with flavonoids in the leaves and roots of S. glabra, including 2 genes in the PAL pathway, 3 genes in the 4CL pathway, 2 genes in the CHS pathway, 4 genes in the CHI pathway, 2 genes in the FLS pathway, 1 gene in the DFR pathway, 1 gene in the CYP73A pathway, 1 gene in the CYP75B1 pathway, 3 genes in the PGT1 pathway, 6 genes in the HCT pathway, 2 genes in the C3'H pathway, 1 gene in the CCOAOMT pathway, 1 gene in the ANR pathway, 1 gene in the LAR pathway, 2 genes in the 3AT pathway, 1 gene in the BZ1 pathway, 2 genes in the IFTM7 pathway, and 1 gene in the CYP81E9 pathway. Six transcription factors, including C2H2, bHLH, and bZIP, were involved in regulating the differential accumulation of flavonoids in the leaves and roots of S. glabra. The qRT-PCR results showed that the up-or down-regulated expression of the 8 randomly selected enzyme genes involved in flavonoid synthesis in the leaves and roots of S. glabra was consistent with the transcriptome sequencing results. This study preliminarily analyzed the transcriptional regulation mechanism of differential accumulation of flavonoids in the leaves and roots of S. glabra, laying a foundation for further elucidating the regulatory effects of key enzyme genes and corresponding transcription factors on the accumulation of flavonoids in S. glabra.

Selection and validation of reference genes for quantitative real-time PCR analysis in Paeonia veitchii / 中国中药杂志

Paeonia veitchii and P. lactiflora are both original plants of the famous Chinese medicinal drug Paeoniae Radix Rubra in the Chinese Pharmacopoeia. They have important medicinal value and great potential in the flower market. The selection of stable and reliable reference genes is a necessary prerequisite for molecular research on P. veitchii. In this study, two reference genes, Actin and GAPDH, were selected as candidate genes from the transcriptome data of P. veitchii. The expression levels of the two candidate genes in different tissues(phloem, xylem, stem, leaf, petiole, and ovary) and different growth stages(bud stage, flowering stage, and dormant stage) of P. veitchii were detected using real-time fluorescence quantitative technology(qRT-PCR). Then, the stability of the expression of the two reference genes was comprehensively analyzed using geNorm, NormFinder, BestKeeper, ΔCT, and RefFinder. The results showed that the expression patterns of Actin and GAPDH were stable in different tissues and growth stages of P. veitchii. Furthermore, the expression levels of eight genes(Pv-TPS01, Pv-TPS02, Pv-CYP01, Pv-CYP02, Pv-CYP03, Pv-BAHD01, Pv-UGT01, and Pv-UGT02) in different tissues were further detected based on the transcriptome data of P. veitchii. The results showed that when Actin and GAPDH were used as reference genes, the expression trends of the eight genes in different tissues of P. veitchii were consistent, validating the reliability of Actin and GAPDH as reference genes for P. veitchii. In conclusion, this study finds that Actin and GAPDH can be used as reference genes for studying gene expression levels in different tissues and growth stages of P. veitchii.

Spatial transcriptome analysis of long non-coding RNAs reveals tissue specificity and functional roles in cancer / 浙江大学学报（英文版）（B辑：生物医学和生物技术）

Long non-coding RNAs (lncRNAs) play a significant role in maintaining tissue morphology and functions, and their precise regulatory effectiveness is closely related to expression patterns. However, the spatial expression patterns of lncRNAs in humans are poorly characterized. Here, we constructed five comprehensive transcriptomic atlases of human lncRNAs covering thousands of major tissue samples in normal and disease states. The lncRNA transcriptomes exhibited high consistency within the same tissues across resources, and even higher complexity in specialized tissues. Tissue-elevated (TE) lncRNAs were identified in each resource and robust TE lncRNAs were refined by integrative analysis. We detected 1 to 4684 robust TE lncRNAs across tissues; the highest number was in testis tissue, followed by brain tissue. Functional analyses of TE lncRNAs indicated important roles in corresponding tissue-related pathways. Moreover, we found that the expression features of robust TE lncRNAs made them be effective biomarkers to distinguish tissues; TE lncRNAs also tended to be associated with cancer, and exhibited differential expression or were correlated with patient survival. In summary, spatial classification of lncRNAs is the starting point for elucidating the function of lncRNAs in both maintenance of tissue morphology and progress of tissue-constricted diseases.

Expression analysis of transcription factors in sugarcane during cold stress / Análise de expressão de fatores de transcrição em cana-de-açúcar sob estresse por frio

Abstract Transcription factors (TF) are a wide class of genes in plants, and these can regulate the expression of other genes in response to various environmental stresses (biotic and abiotic). In the current study, transcription factor activity in sugarcane was examined during cold stress. Initially, RNA transcript reads of two sugarcane cultivars (ROC22 and GT08-1108) under cold stress were downloaded from SRA NCBI database. The reads were aligned into a reference genome and the differential expression analyses were performed with the R/Bioconductor edgeR package. Based on our analyses in the ROC22 cultivar, 963 TF genes were significantly upregulated under cold stress among a total of 5649 upregulated genes, while 293 TF genes were downregulated among a total of 3,289 downregulated genes. In the GT08-1108 cultivar, 974 TF genes were identified among 5,649 upregulated genes and 283 TF genes were found among 3,289 downregulated genes. Most transcription factors were annotated with GO categories related to protein binding, transcription factor binding, DNA-sequence-specific binding, transcription factor complex, transcription factor activity in RNA polymerase II, the activity of nucleic acid binding transcription factor, transcription corepressor activity, sequence-specific regulatory region, the activity of transcription factor of RNA polymerase II, transcription factor cofactor activity, transcription factor activity from plastid promoter, transcription factor activity from RNA polymerase I promoter, polymerase II and RNA polymerase III. The findings of above results will help to identify differentially expressed transcription factors during cold stress. It also provides a comprehensive analysis of the regulation of the transcription activity of many genes. Therefore, this study provides the molecular basis for improving cold tolerance in sugarcane and other economically important grasses.

Resumo Fatores de transcrição (FT) são uma ampla classe de genes em plantas e podem regular a expressão de outros genes em resposta a vários estresses ambientais (estresses bióticos e abióticos). No presente estudo, a atividade do fator de transcrição na cana-de-açúcar foi examinada durante o estresse pelo frio. Inicialmente, as leituras de transcrição de RNA de duas cultivares de cana-de-açúcar (ROC22 e GT08-1108) sob estresse frio foram baixadas do banco de dados SRA NCBI. As leituras foram alinhadas em um genoma de referência e as análises de expressão diferencial foram realizadas com o pacote R / Bioconductor edgeR. Com base em nossas análises no cultivar ROC22, 963 genes TF foram significativamente regulados positivamente sob estresse pelo frio entre um total de 5.649 genes regulados positivamente, enquanto 293 genes TF foram regulados negativamente entre um total de 3.289 genes regulados negativamente. No cultivar GT08-1108, 974 genes TF foram identificados entre 5.649 genes regulados positivamente e 283 genes TF foram encontrados entre 3.289 genes regulados negativamente. Os fatores de transcrição, em sua maioria, foram anotados com categorias GO relacionadas à ligação de proteína, ligação de fator de transcrição, ligação específica de sequência de DNA, complexo de fator de transcrição, atividade de fator de transcrição em RNA polimerase II, atividade de fator de transcrição de ligação de ácido nucleico, atividade de corepressor de transcrição, sequência específica da região reguladora, atividade do fator de transcrição da RNA polimerase II, atividade do cofator do fator de transcrição, atividade do fator de transcrição do promotor do plastídio, atividade do fator de transcrição do promotor da RNA polimerase I, polimerase II e RNA polimerase III. As descobertas dos resultados acima ajudarão a identificar fatores de transcrição expressos diferencialmente durante o estresse pelo frio. Ele também fornece uma análise abrangente da regulação da atividade de transcrição de muitos genes. Portanto, este estudo fornece base molecular para melhorar a tolerância ao frio em cana-de-açúcar e outras gramíneas economicamente importantes.

Genome-wide identification of bZIP family genes and screening of candidate AarbZIPs involved in terpenoid biosynthesis in Artemisia argyi / 中国中药杂志

Artemisia argyi is an important medicinal and economic plant in China, with the effects of warming channels, dispersing cold, and relieving pain, inflammation, and allergy. The essential oil of this plant is rich in volatile terpenoids and widely used in moxi-bustion and healthcare products, with huge market potential. The bZIP transcription factors compose a large family in plants and are involved in the regulation of plant growth and development, stress response, and biosynthesis of secondary metabolites such as terpenoids. However, little is known about the bZIPs and their roles in A. argyi. In this study, the bZIP transcription factors in the genome of A. argyi were systematically identified, and their physicochemical properties, phylogenetic relationship, conserved motifs, and promoter-binding elements were analyzed. Candidate AarbZIP genes involved in terpenoid biosynthesis were screened out. The results showed that a total of 156 AarbZIP transcription factors were identified at the genomic level, with the lengths of 99-618 aa, the molecular weights of 11.7-67.8 kDa, and the theoretical isoelectric points of 4.56-10.16. According to the classification of bZIPs in Arabidopsis thaliana, the 156 AarbZIPs were classified into 12 subfamilies, and the members in the same subfamily had similar conserved motifs. The cis-acting elements of promoters showed that AarbZIP genes were possibly involved in light and hormonal pathways. Five AarbZIP genes that may be involved in the regulation of terpenoid biosynthesis were screened out by homologous alignment and phylogenetic analysis. The qRT-PCR results showed that the expression levels of the five AarbZIP genes varied significantly in different tissues of A. argyi. Specifically, AarbZIP29 and AarbZIP55 were highly expressed in the leaves and AarbZIP81, AarbZIP130, and AarbZIP150 in the flower buds. This study lays a foundation for the functional study of bZIP genes and their regulatory roles in the terpenoid biosynthesis in A. argyi.

Genome-wide identification of the BmAKR gene family in the silkworm (Bombyx mori) and their expression analysis in diapause eggs and nondiapause eggs / 生物工程学报

The aldo-keto reductase super family (AKRs) has a wide range of substrate specificity. However, the systematic identification of insect AKR gene family members has not been reported. In this study, bioinformatics methods were used to predict the phylogenetic evolution, physical and chemical properties, chromosome location, conserved motifs, and gene structure of AKR family members in Bombyx mori (BmAKR). Transcriptome data or quantitative real time polymerase chain reaction (qRT-PCR) were used to analyze the expression level of BmAKR genes during different organizational periods and silkworm eggs in different developmental states. Moreover, Western blotting was used to detect the expression level of the BmAKR in silkworm eggs. The results showed that 11 BmAKR genes were identified. These genes were distributed on 4 chromosomes of the silkworm genome, all of which had the (α/β) 8-barrel motif, and their physical and chemical characteristics were relatively similar. Phylogenetic analysis showed that the BmAKR genes could be divided into 2 subgroups (AKR1 and AKR2). Transcriptome data analysis showed that the expression of BmAKR genes were quite different in different tissues and periods. Moreover, the expression analysis of BmAKR genes in silkworm eggs showed that some genes were expressed significantly higher in nondiapause eggs than in diapause eggs; but another gene, BmAKR1-1, was expressed significantly higher in diapause eggs than in nondiapause eggs. The detection of protein level found that the difference trend of BmAKR1-1 in diapause eggs and non-diapause eggs was consistent with the results of qRT-PCR. In conclusion, BmAKR1-1 was screened out as candidates through the identification and analysis of the BmAKR genes in silkworm, which may regulate silkworm egg development is worthy of further investigation.

Exploratory research on the probable shared molecular mechanism and transcription factors between chronic periodontitis and chronic obstructive pulmonary disease / 华西口腔医学杂志

OBJECTIVES@#To investigate possible cross-talk genes, associated pathways, and transcription factors between chronic periodontitis (CP) and chronic obstructive pulmonary disease (COPD).@*METHODS@#The gene expression profiles of CP (GSE10334 and GSE16134) and COPD (GSE76925) were downloaded from the GEO database. Differential expression and functional clustering analyses were performed. The protein‑protein interaction (PPI) network was constructed. The core cross-talk genes were filtered using four topological analysis algorithms and modular segmentation. Then, functional clustering analysis was performed again.@*RESULTS@#GSE10334 detected 164 differentially expressed genes (DEGs) (119 upregulated and 45 downregulated). GSE16134 identified 208 DEGs (154 upregulated and 54 downregulated). GSE76925 identified 1 408 DEGs (557 upregulated and 851 downregulated). The PPI network included 21 nodes and 20 edges. The final screening included seven cross-talk genes: CD79A, FCRLA, CD19, IRF4, CD27, SELL, and CXCL13. Relevant pathways included primary immunodeficiency, the B-cell receptor signaling pathway, and cytokine-cytokine receptor interaction.@*CONCLUSIONS@#This study indicates the probability of shared pathophysiology between CP and COPD, and their cross-talk genes, associated pathways, and transcription factors may offer novel concepts for future mechanistic investigations.

Characterization of candidate factors associated with the metastasis and progression of high-grade serous ovarian cancer / 中华医学杂志(英文版)

BACKGROUND@#High-grade serous ovarian cancer (HGSOC) is the biggest cause of gynecological cancer-related mortality because of its extremely metastatic nature. This study aimed to explore and evaluate the characteristics of candidate factors associated with the metastasis and progression of HGSOC.@*METHODS@#Transcriptomic data of HGSOC patients' samples collected from primary tumors and matched omental metastatic tumors were obtained from three independent studies in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were selected to evaluate the effects on the prognosis and progression of ovarian cancer using data from The Cancer Genome Atlas (TCGA) database. Hub genes' immune landscapes were estimated by the Tumor Immune Estimation Resource (TIMER) database. Finally, using 25 HGSOC patients' cancer tissues and 10 normal fallopian tube tissues, immunohistochemistry (IHC) was performed to quantify the expression levels of hub genes associated with International Federation of Gynecology and Obstetrics (FIGO) stages.@*RESULTS@#Fourteen DEGs, ADIPOQ , ALPK2 , BARX1 , CD37 , CNR2 , COL5A3 , FABP4 , FAP , GPR68 , ITGBL1 , MOXD1 , PODNL1 , SFRP2 , and TRAF3IP3 , were upregulated in metastatic tumors in every database while CADPS , GATA4 , STAR , and TSPAN8 were downregulated. ALPK2 , FAP , SFRP2 , GATA4 , STAR , and TSPAN8 were selected as hub genes significantly associated with survival and recurrence. All hub genes were correlated with tumor microenvironment infiltration, especially cancer-associated fibroblasts and natural killer (NK) cells. Furthermore, the expression of FAP and SFRP2 was positively correlated with the International Federation of Gynecology and Obstetrics (FIGO) stage, and their increased protein expression levels in metastatic samples compared with primary tumor samples and normal tissues were confirmed by IHC ( P = 0.0002 and P = 0.0001, respectively).@*CONCLUSIONS@#This study describes screening for DEGs in HGSOC primary tumors and matched metastasis tumors using integrated bioinformatics analyses. We identified six hub genes that were correlated with the progression of HGSOC, particularly FAP and SFRP2 , which might provide effective targets to predict prognosis and provide novel insights into individual therapeutic strategies for HGSOC.

Altered microRNA expression profiles of human spermatozoa in normal fertile men of different ages / 亚洲男科学杂志(英文版)

MicroRNAs (miRNAs) are mediators of the aging process. The purpose of this work was to analyze the miRNA expression profiles of spermatozoa from men of different ages with normal fertility. Twenty-seven donors were divided into three groups by age (Group A, n = 8, age: 20-30 years; Group B, n = 10, age: 31-40 years; and Group C, n = 9, age: 41-55 years) for high-throughput sequencing analysis. Samples from 65 individuals (22, 22, and 21 in Groups A, B, and C, respectively) were used for validation by quantitative real-time polymerase chain reaction (qRT-PCR). A total of 2160 miRNAs were detected: 1223 were known, 937 were newly discovered and unnamed, of which 191 were expressed in all donors. A total of 7, 5, and 17 differentially expressed microRNAs (DEMs) were found in Group A vs B, Group B vs C, and Group A vs C comparisons, respectively. Twenty-two miRNAs were statistically correlated with age. Twelve miRNAs were identified as age-associated miRNAs, including hsa-miR-127-3p, mmu-miR-5100_L+2R-1, efu-miR-9226_L-2_1ss22GA, cgr-miR-1260_L+1, hsa-miR-652-3p_R+1, pal-miR-9993a-3p_L+2R-1, hsa-miR-7977_1ss6AG, hsa-miR-106b-3p_R-1, hsa-miR-186-5p, PC-3p-59611_111, hsa-miR-93-3p_R+1, and aeca-mir-8986a-p5_1ss1GA. There were 9165 target genes of age-associated miRNAs. Gene Ontology (GO) analysis of the target genes identified revealed enrichment of protein binding, membrane, cell cycle, and so on. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of age-related miRNAs for target genes revealed 139 enriched pathways, such as signaling pathways regulating stem cell pluripotency, metabolic pathways, and the Hippo signaling pathway. This suggests that miRNAs play a key role in male fertility changes with increasing age and provides new evidence for the study of the mechanism of age-related male fertility decline.

Bioinformatics analysis and identification to immune-related markers of osteoporosis / 细胞与分子免疫学杂志

Objective To identify immune-related dysregulation mechanisms and potential diagnostic predictive biomarkers in osteoporosis. Methods Gene expression data for both osteoporosis and control populations were retrieved from the GSE35958 and GSE56815 datasets. Immune-related differentially expressed genes (DEGs) were obtained by screening DEGs and were compared with the immunology database and analysis portal (ImmPort) database. Enrichment analysis of these immune-related DEGs was conducted using the Clusterprofiler software package. A protein-protein interaction network was built with the STRING database, which is a search tool for finding interacting genes/proteins, and the top 10 genes with the highest network connectivity were identified as candidate genes. Subsequently, the diagnostic predictive effect of candidate genes was evaluated using receiver operating characteristic (ROC) curves, logistic regression, and column plots. Finally, PCR and Western blot analysis were applied to detect the differential expression of these genes in bone marrow tissue of patients with osteoporosis. Results A total of 138 immune-related DEGs were obtained through intersection analysis. The results of the enrichment analysis indicated that these genes were involved in biological functions such as immune inflammation and signaling pathways including T cell receptors, mitogen activated protein kinase (MAPK), rat sarcoma virus oncogene homologs (Ras), osteoclast differentiation, and B cell receptors. In addition, among the candidate genes, upregulated vascular endothelial growth factor A (VEGFA) and epidermal growth factor receptor (EGFR) and downregulated AKT1, SRC, and JUN in osteoporosis showed the highest connectivity. Among them, VEGFA, EGFR, JUN, and AKT1 demonstrated the best diagnostic predictive value. Conclusion The screening of immune-related DEGs will enhance the understanding of osteoporosis and facilitate the development of immunotherapy targets.

Screening of Genes Co-Associated with Sudden Infant Death Syndrome and Infectious Sudden Death in Infancy and Bioinformatics Analysis of Their Regulatory Networks / 法医学杂志

OBJECTIVES@#The common differentially expressed mRNAs in brain, heart and liver tissues of deceased sudden infant death syndrome (SIDS) and infectious sudden death in infancy (ISDI) confirmed by autopsy was screened by bioinformatics to explore the common molecular markers and pathogenesis of SIDS and ISDI.@*METHODS@#The datasets of GSE70422 and GSE136992 were downloaded, the limma of R software was used to screen differentially expressed mRNA in different tissue samples of SIDS and ISDI decedents for overlapping analysis. The clusterProfiler of R software was used to conduct gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. The protein-protein interaction (PPI) network was constructed by STRING database, while the hub gene was screened by cytoHubba plug-in.@*RESULTS@#Compared with the control group, there were 19 significant differentially expressed genes in the tissue samples of SIDS and ISDI decedents, among which 16 in the heart tissue and 3 in the liver tissue, and the astrotactin 1 (ASTN1) gene expression difference in the heart tissue was most significant. The PPI network identified Ras homolog family member A (RHOA), integrin subunit alpha 1 (ITGA1), and H2B clustered histone 5 (H2BC5) were hub genes. The analysis of GO and KEGG showed that differentially expressed genes were enriched in the molecular pathways of actin cytoskeleton regulation, focal adhesion and response to mycophenolic acid.@*CONCLUSIONS@#ASTN1, RHOA and ITGA1 may participate in the development of SIDS and ISDI. The enrichment of differentially expressed genes in immune and inflammatory pathways suggests a common molecular regulatory mechanism between SIDS and ISDI. These findings are expected to provide new biomarkers for molecular anatomy and forensic identification of SIDS and ISDI.

Transcriptional profiling of dental sensory and proprioceptive trigeminal neurons using single-cell RNA sequencing / 国际口腔科学杂志·英文版

Dental primary afferent (DPA) neurons and proprioceptive mesencephalic trigeminal nucleus (MTN) neurons, located in the trigeminal ganglion and the brainstem, respectively, are essential for controlling masticatory functions. Despite extensive transcriptomic studies on various somatosensory neurons, there is still a lack of knowledge about the molecular identities of these populations due to technical challenges in their circuit-validated isolation. Here, we employed high-depth single-cell RNA sequencing (scRNA-seq) in combination with retrograde tracing in mice to identify intrinsic transcriptional features of DPA and MTN neurons. Our transcriptome analysis revealed five major types of DPA neurons with cell type-specific gene enrichment, some of which exhibit unique mechano-nociceptive properties capable of transmitting nociception in response to innocuous mechanical stimuli in the teeth. Furthermore, we discovered cellular heterogeneity within MTN neurons that potentially contribute to their responsiveness to mechanical stretch in the masseter muscle spindles. Additionally, DPA and MTN neurons represented sensory compartments with distinct molecular profiles characterized by various ion channels, receptors, neuropeptides, and mechanoreceptors. Together, our study provides new biological insights regarding the highly specialized mechanosensory functions of DPA and MTN neurons in pain and proprioception.