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.

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.

The pattern of alternative splicing and DNA methylation alteration and their interaction in linseed (Linum usitatissimum L) response to repeated drought stresses

BACKGROUND: Drought stress has significantly hampered agricultural productivity worldwide and can also result in modifications to DNA methylation levels. However, the dynamics of DNA methylation and its association with the changes in gene transcription and alternative splicing (AS) under drought stress are unknown in linseed, which is frequently cultivated in arid and semiarid regions. RESULTS: We analysed AS events and DNA methylation patterns in drought-tolerant (Z141) and drought-sensitive (NY-17) linseed under drought stress (DS) and repeated drought stress (RD) treatments. We found that the number of intron-retention (IR) and alternative 3' splice site (Alt3'SS) events were significantly higher in Z141 and NY-17 under drought stress. We found that the linseed response to the DS treatment was mainly regulated by transcription, while the response to the RD treatment was coregulated by transcription and AS. Whole genome-wide DNA methylation analysis revealed that drought stress caused an increase in the overall methylation level of linseed. Although we did not observe any correlation between differentially methylated genes (DMGs) and differentially spliced genes (DSGs) in this study, we found that the DSGs whose gene body region was hypermethylated in Z141 and hypomethylated in NY-17 were enriched in abiotic stress response Gene Ontology (GO) terms. This finding implies that gene body methylation plays an important role in AS regulation in some specific genes. CONCLUSION: Our study is the first comprehensive genome-wide analysis of the relationship between linseed methylation changes and AS under drought and repeated drought stress. Our study revealed different interaction patterns between differentially expressed genes (DEGs) and DSGs under DS and RD treatments and differences between methylation and AS regulation in drought-tolerant and drought-sensitive linseed varieties. The findings will probably be of interest in the future. Our results provide interesting insights into the association between gene expression, AS, and DNA methylation in linseed under drought stress. Differences in these associations may account for the differences in linseed drought tolerance.

The effects of repeated doses of xylazine-ketamine and medetomidineketamine anesthesia on DNA damage in the liver and kidney

Purpose: This study evaluated the DNA damage caused by repeated doses of xylazine-ketamine and medetomidine-ketamine anesthesia in the liver and kidneys. Methods: In this study, 60 rats were used. The rats were divided into group 1 (xylazine-ketamine), and group 2 (medetomidine-ketamine), and these anesthetic combinations were administered to the rats at repeated doses with 30-min intervals. The effects of these anesthetic agents on the tumor necrosis factor-alpha gene for DNA damage were investigated. Results: According to the gene expression results, it was observed that a single dose of xylazine-ketamine was 2.9-fold expressed, while first and second repeat doses did not show significant changes in expression levels. However, in the case of the third repetition, it was observed to be 3.8-fold overexpressed. In the case of medetomidine-ketamine administration, it was observed that a single-dose application resulted in a 1.04-fold expression, while the first and the third repeat doses showed a significant down expression. The samples from the second repeat dose administration group were found to have insignificant levels of expression. Conclusions: This study can contribute to understanding the safe anesthetic combination in research and operations in which xylazine-ketamine and medetomidine-ketamine combinations are used.

O impacto da exposição química sobre o transcriptoma de trabalhadores formais e informais / The impact of chemical exposure on the transcriptome of formal and informal workers

Introdução: Exposições químicas podem variar em populações geograficamente distintas, com diferentes hábitos, estilos de vida e características individuais. Alguns elementos químicos encontrados no ambiente são capazes de alterar a expressão gênica humana. Objetivos: a) quantificar as concentrações de elementos potencialmente tóxicos (EPTs: As, Cd, Cr, Cu, Hg, Mn, Ni, Pb, Sb, Sn, e Zn) na urina da população de Limeira, e nas peças de joias e pó de solda; b) quantificar as concentrações de EPTs no sangue dos participantes de Limeira e Volta Redonda; c) avaliar os riscos de doenças associadas à exposição ocupacional; d) avaliar o impacto da exposição ocupacional na expressão gênica de trabalhadores formais e informais. Métodos: O grupo Exposto foi composto por trabalhadores informais que realizam soldagem de joias e bijuterias em ambiente domiciliar na cidade de Limeira, SP; e por trabalhadores formais da Companhia Siderúrgica Nacional, em Volta Redonda, RJ. O grupo Controle incluiu moradores dos mesmos bairros dos trabalhadores, mas que não desenvolviam nenhuma atividade diretamente relacionada à exposição química. Amostras de sangue foram coletadas para quantificação de glicose, insulina, perfil lipídico, EPTs e para análise transcriptômica. Em Limeira, também foi quantificada a concentração de EPTs na urina. Para transcriptômica, o RNA foi extraído e hibridizado com Agilent SurePrint G3 Human Gene Expression 8x60K v2 Microarray. O pré-processamento, análise estatística e de vias de interesse foram realizados no software GeneSpring GX. Todos os participantes preencheram questionários sobre hábitos, percepção de risco, morbidade referida e exposição ocupacional. A associação entre exposição a EPTs e desfechos de saúde foi testada por modelo de regressão de Poisson multivariado. Resultados: Nos trabalhadores informais, foram detectados 16 genes superexpressos e 33 subexpressos em comparação com os controles (fold-change > 2). A análise de vias indicou genes enriquecidos em vias do processo inflamatório (quimiocinas MAPK, receptor Toll-like e NF-kappa B). Nos trabalhadores formais, foram encontrados 20 genes superexpressos e 50 subexpressos, com vias relacionadas à resposta imune e ao processo de aterosclerose. O único gene diferencialmente expresso (DEG) em comum nas duas populações foi o IFI27 relacionado na literatura a diferentes tipos de câncer. A produção informal de joias de Limeira foi associada a exposição dos trabalhadores ao Cd, com concentrações significativamente maiores na urina e no sangue dos trabalhadores comparado aos controles. Além disso, foi observada uma associação positiva entre as concentrações de Cd no sangue e a glicemia. As concentrações de As e Pb também foram maiores no sangue dos trabalhadores informais comparado aos controles, sendo que participantes com concentrações de Pb superiores a 2,6 µg dL-1 apresentaram prevalência de manifestações neurológicas 2,3 vezes maior (IC 95%: 1,17 - 4,58; p = 0,02). Não foram observadas diferenças significativas nos EPTs entre os grupos de Volta Redonda, provavelmente devido ao uso de equipamentos de proteção individual e à poluição ambiental na região. Conclusão: As diferenças na expressão gênica relacionadas à exposição ocupacional estão associadas, principalmente, à inflamação e à resposta imune. Os resultados sugerem que a exposição ocupacional prolongada a elementos tóxicos pode levar a consequências negativas para a saúde, como por exemplo, um aumento da prevalência de manifestações neurológicas. Os resultados exploratórios desta tese são um ponto de partida para estudos em populações sensíveis e pouco estudadas, especialmente, de países em desenvolvimento. Análises adicionais devem ser realizadas para investigar efeitos diretos e validar associações causais.

Introduction: Chemical exposures may vary in geographically distinct populations, with different habits, lifestyles, and individual characteristics. Objectives: a) to determine potentially toxic elements' (EPTs: As, Cd, Cr, Cu, Hg, Mn, Ni, Pb, Sb, Sn, and Zn) in the urine of the population of Limeira, and in jewelry pieces and soldering powder; b) determine the PTEs concentrations in the participants' blood from Limeira and Volta Redonda municipalities; c) investigate disease risks associated with occupational exposure; d) to evaluate the impact of occupational exposure on gene expression profile. Methods: The exposed group was composed of informal workers who perform soldering of jewelry inside their houses in the city of Limeira, SP; and formal workers from a steel company in the city of Volta Redonda, RJ. Control participants were recruited from the same neighborhoods without occupational chemical exposure. Blood samples were collected for blood glucose, insulin, lipid profile, and PTE determinations, and for transcriptomic analysis. In Limeira, PTE concentration in urine was also determined. RNA was extracted and hybridized to Agilent SurePrint G3 Human Gene Expression 8x60K v2 Microarray for transcriptomics analysis. Pre-processing, statistical, and pathway analysis were performed in GeneSpring GX software. All participants completed questionnaires about household risk, reported morbidity, and occupational exposure. The association between PTEs exposure and health outcomes was tested by a multivariable robust Poisson regression model. Results: 16 up- and 33 down-regulated genes (fold-change > 2) were observed in the informal workers. Pathway analysis revealed genes enriched in inflammatory process (MAPK, Toll-like receptor, and NF-kappa B chemokine signaling pathways). In formal workers, 20 up- and 50 down-regulated genes were found with pathways related to immune response and atherosclerosis development. The gene IFI27 which has been associated with various types of cancer was the only one commonly differentially expressed between informal and formal workers. Informal jewelry production in Limeira increased exposure to Cd, with significantly higher concentrations in the urine and blood of informal exposed workers compared to controls. Furthermore, a positive association was observed between blood Cd concentrations and glycemia. The blood concentration of As and Pb were also Participants with Pb concentrations higher than 2.6 µg dL-1 showed a prevalence of neurological manifestations 2.3 times higher (95% CI: 1.17 - 4.58; p = 0.02) than those with lower lead concentrations. No significant differences were observed between formal workers from Volta Redonda and their control group, probably, because of the use of individual protection equipment and the environmental pollution in the region. Conclusion: Differences in gene expression related to occupational exposure are mainly associated with inflammation and immune response. The results suggest that prolonged occupational exposure to toxic elements could lead to negative health outcomes, such as higher prevalence of neurological manifestations. These exploratory results are a starting point for analysis in sensitive and understudied populations, especially in developing countries. Further analysis should be carried out to investigate its direct effects and to validate causal associations.

Avaliação dos aspectos moleculares e sua relevância na resposta terapêutica e prognóstico de pacientes com câncer de bexiga submetidos a quimioterapia com cisplatina / Evaluation of molecular aspects and their relevance in therapeutic response and prognosis of patients with bladder cancer undergoing chemotherapy with cisplatin

INTRODUÇÃO: O câncer de bexiga é a 10a neoplasia maligna mais comum no mundo. O tipo músculo-invasivo tem como tratamento padrão-ouro combinações de quimioterapia e cirurgia. Muitos pacientes, no entanto, não respondem a quimioterapia e são submetidos a uma terapia ineficaz. Identificar quais pacientes não irão se beneficiar da terapia pode reduzir toxicidades e custos. A base de cisplatina, a quimioterapia no tumor de bexiga age através de dano ao DNA. Levantamos a hipótese que pacientes que respondem apresentam mutação ou baixa expressão de genes de reparo ao DNA. OBJETIVO: Comparar os aspectos moleculares entre 2 grupos de pacientes com tumor de bexiga submetidos a quimioterapia: um que respondeu adequadamente e outro que não respondeu. Identificar se os genes diferencialmente expressos ou mutados estão associados ao reparo do DNA e, consequentemente, com a resposta terapêutica. MÉTODOS: O estudo foi realizado através de análise de dados moleculares e clínicos de pacientes com tumor de bexiga músculo-invasivo submetidos a quimioterapia com cisplatina. Os pacientes foram divididos em dois grupos de acordo com a resposta a quimioterapia: Respondedores e Não respondedores. Realizamos comparação do perfil de expressão gênica, regulação epigenética e mutações entre grupos. Os genes diferencialmente expressos e/ou mutados foram avaliados se estão associados com as vias de reparo do DNA e, consequentemente, com a reposta as terapias sistêmicas. Por fim, realizamos análise de sobrevida a fim de avaliar correlação entre expressão dos genes identificados com a sobrevida global dos pacientes. RESULTADOS: A casuística foi composta por 63 pacientes. 41 apresentaram resposta a quimioterapia e 22 não responderam. Após análise de expressão gênica diferencial, identificamos expressão significativamente menor dos genes DDB1 e PRPF9, associados as vias de reparo do DNA, no grupo de pacientes respondedores. Na análise de sobrevida, identificamos que a alta expressão do gene DDB1 esteve significativamente associada a uma sobrevida global pior. Na comparação de padrão mutacional, identificamos uma taxa de mutação, significativamente maior, do gene CDK12, associado ao reparo do DNA, no grupo de pacientes respondedores. Esses padrões de expressão ou mutação identificados nestes 3 genes de reparo de dano ao DNA estiveram fortemente associados a uma resposta satisfatória a quimioterapia. Dessa forma, são potenciais biomarcadores capazes de predizer resposta a quimioterapia em câncer de bexiga. CONCLUSÕES: Identificamos potenciais biomarcadores associados com a resposta terapêutica e sobrevida após uso de cisplatina em pacientes com tumor de bexiga. Estes, após validações, potencialmente poderão ser agrupado em escores capazes de predizer resposta a quimioterapia no câncer de bexiga e auxiliar na tomada de decisão, implicando em redução de toxicidades e custos.

INTRODUCTION: Bladder cancer is the 10th most common malignancy in world. The gold standard treatment for muscle-invasive type is combinations of chemotherapy and surgery. Many patients, however, do not respond to chemotherapy and are subjected to ineffective therapy. Identifying which patients will not benefit from therapy can reduce toxicities and costs. Cisplatin-based chemotherapy in bladder tumors acts through DNA damage. We hypothesize that patients who respond have mutation or low expression of DNA repair genes. OBJECTIVE: To compare the molecular aspects between 2 groups of patients with bladder tumors undergoing chemotherapy: one that responded adequately and other that did not respond. Identify whether differentially expressed or mutated genes are associated with DNA repair and, consequently, with therapeutic response. METHODS: The study was carried out through the analysis of molecular and clinical data from patients with muscle-invasive bladder tumors undergoing chemotherapy with cisplatin. Patients were divided into two groups according to their response to chemotherapy: Responders and Non-responders. We compared gene expression profile, epigenetic regulation and mutations between groups. Differentially expressed and/or mutated genes were assessed whether they are associated with DNA repair pathways and, consequently, with response to systemic therapies. Finally, we performed a survival analysis in order to evaluate the correlation between the expression of the identified genes and the overall survival of the patients. RESULTS: The sample consisted of 63 patients. 41 responded to chemotherapy and 22 did not respond. After differential gene expression analysis, we identified significantly lower expression of the DDB1 and PRPF9 genes, associated with DNA repair pathways, in Responders group. In survival analysis, we identified that high expression of the DDB1 gene was significantly associated with worse overall survival. When comparing the mutational pattern, we identified a significantly higher mutation rate in the CDK12 gene, associated with DNA repair, in Responders group. These expression or mutation patterns identified in these 3 DNA damage repair genes were strongly associated with a satisfactory response to chemotherapy. Therefore, they are potential biomarkers capable of predicting response to chemotherapy in bladder cancer. CONCLUSIONS: We identified potential biomarkers associated with therapeutic response and survival after the use of cisplatin in patients with bladder tumors. These, after validation, could potentially be grouped into scores capable of predicting response to chemotherapy in bladder cancer and assisting in decisionmaking, resulting in a reduction in toxicities and costs.

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.

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.

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.

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.

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.