RTCpredictor: identification of read-through chimeric RNAs from RNA sequencing data.

Read-through chimeric RNAs are being recognized as a means to expand the functional transcriptome and contribute to cancer tumorigenesis when mis-regulated. However, current software tools often fail to predict them. We have developed RTCpredictor, utilizing a fast ripgrep tool to search for all possible exon-exon combinations of parental gene pairs. We also added exonic variants allowing searches containing common SNPs. To our knowledge, it is the first read-through chimeric RNA specific prediction method that also provides breakpoint coordinates. Compared with 10 other popular tools, RTCpredictor achieved high sensitivity on a simulated and three real datasets. In addition, RTCpredictor has less memory requirements and faster execution time, making it ideal for applying on large datasets.

Rhabdomyosarcomas are oncogene addicted to the activation of AVIL.

Rhabdomyosarcoma (RMS) is one of the most common pediatric soft-tissue cancer. Previously, we discovered a gene fusion, MARS-AVIL formed by chromosomal inversion in RMS. Suspecting that forming a fusion with a housekeeping gene may be one of the mechanisms to dysregulate an oncogene, we investigated AVIL expression and its role in RMS. We first showed that MARS-AVIL translates into an in-frame fusion protein, which is critical for RMS cell tumorigenesis. Besides forming a gene fusion with the housekeeping gene, MARS, the AVIL locus is often amplified, and its RNA and protein expression are overexpressed in the majority of RMSs. Tumors with AVIL dysregulation exhibit evidence of oncogene addiction: Silencing MARS-AVIL in cells harboring the fusion, or silencing AVIL in cells with AVIL overexpression, nearly eradicated the cells in culture, as well as inhibited in vivo xenograft growth in mice. Conversely, gain-of-function manipulations of AVIL led to increased cell growth and migration, enhanced foci formation in mouse fibroblasts, and most importantly transformed mesenchymal stem cells in vitro and in vivo. Mechanistically, AVIL seems to serve as a converging node functioning upstream of two oncogenic pathways, PAX3-FOXO1 and RAS, thus connecting two types of RMS associated with these pathways. Interestingly, AVIL is overexpressed in other sarcoma cells as well, and its expression correlates with clinical outcomes, with higher levels of AVIL expression being associated with worse prognosis. AVIL is a bona fide oncogene in RMS, and RMS cells are addicted to its activity.

Exosomes from Intrahepatic Cholestasis of Pregnancy Induce Cell Apoptosis Through the miRNA-6891-5p/YWHAE Pathway.

BACKGROUND: Intrahepatic cholestasis of pregnancy (ICP) is associated with adverse pregnancy outcomes; however, the underlying mechanisms are not fully understood. AIMS: This study aimed to determine the role of exosomal miR-6891-5p in placental trophoblast dysfunction in ICP and identify new biomarkers for ICP diagnosis. METHODS: Serum samples were collected from ICP patients and healthy pregnant women, and serum exosomes were extracted and identified. Fluorescent dye labeling of exosomes and cell-verified cell phagocytosis were performed. In vitro experiments were conducted by adding taurocholic acid to simulate the ICP environment. Cell proliferation and apoptosis levels were detected using flow cytometry and the cell counting kit-8 assay. Mimics were constructed to overexpress miR-6891-5p in cells, and the binding site between miR-6891-5p and YWHAE was verified using luciferase reporter genes. RESULTS: miR-6891-5p expression was significantly decreased in serum exosomes of ICP patients. Co-culturing with exosomes derived from ICP patients' serum (ICP-Exos) decreased HTR-8/SVeno cell proliferation and increased apoptosis levels. miR-6891-5p upregulation in HTR-8/SVeno cells significantly increased cell viability and reduced cell apoptosis levels, as determined by the cell counting kit-8 assay and flow cytometry. A double luciferase assay confirmed that miR-6891-5p affected the expression of the downstream YWHAE protein. CONCLUSIONS: This study indicates that serum exosomes from ICP patients can impact the apoptosis of placental trophoblast HTR-8/SVeno cells through the miR-6891-5P/YWHAE pathway and can serve as specific molecular markers for ICP diagnosis.

Gold nanoparticle dimer-based immunochromatography for in situ ultrasensitive detection of porcine epidemic diarrhea virus.

The low detection sensitivity of lateral-flow immunochromatography assay (LFIA) based on spherical gold nanoparticle (AuNP) limits its wide applications. In the present study, AuNP dimers with strong plasma scattering and robust signal output were synthesized via the Ag ion soldering (AIS) strategy and used as labeled probes in LFIA to boost the sensitivity without any extra operation process and equipment. The established LFIA exhibited high sensitivity with a limit of detection (LOD) of 2.0 × 102 TCID50/mL for PEDV, which provides 50 times higher sensitivity than commercial LFIA based on spherical colloidal gold. In addition, the AuNP dimer-based LFIA showed strong specificity, good reproducibility, high stability, and good accordance to reverse transcription polymer chain reaction (RT-PCR) when detecting 109 clinical samples. Thus, the AuNP dimers is a promising probe for LFIA and the developed AuNP dimer-based LFIA is suitable for the rapid detection of PEDV in the field.

Investigation of interfacial adsorption between microplastics and methylparaben in aqueous solution.

Microplastics and parabens are considered to be a global contaminants, especially in the aquatic ecosystem. The interfacial interaction between four types of microplastics including polystyrene, polyethylene, polyethylene terephthalate, and polyvinyl chloride, and methylparaben were investigated in this study. The results showed that molecular layer dominates the adsorption, with the rate significantly affected by both internal diffusion and external diffusion. Among the four types, polystyrene and polyvinyl chloride showed the smallest and biggest adsorption capability, with the values were 0.656 and 1.269 mg g-1, respectively. For the adsorption capability, smaller particle size and higher pH value possessed positive effects. However, the existence of metal ions could inhibit the adsorption process, except for a weak promotion at low salinity. Physical adsorption effects, such as electrostatic interaction, hydrogen bond formation, and covalent bond formation, had been identified that dominated the adsorption. This finding could be served as a speculative foundation for the further study of the toxicity, migration, and ecological risk assessment of microplastics in aquatic ecosystem.

Dynamic landscape and evolution of m6A methylation in human.

m6A is a prevalent internal modification in mRNAs and has been linked to the diverse effects on mRNA fate. To explore the landscape and evolution of human m6A, we generated 27 m6A methylomes across major adult tissues. These data reveal dynamic m6A methylation across tissue types, uncover both broadly or tissue-specifically methylated sites, and identify an unexpected enrichment of m6A methylation at non-canonical cleavage sites. A comparison of fetal and adult m6A methylomes reveals that m6A preferentially occupies CDS regions in fetal tissues. Moreover, the m6A sub-motifs vary between fetal and adult tissues or across tissue types. From the evolutionary perspective, we uncover that the selection pressure on m6A sites varies and depends on their genic locations. Unexpectedly, we found that ∼40% of the 3'UTR m6A sites are under negative selection, which is higher than the evolutionary constraint on miRNA binding sites, and much higher than that on A-to-I RNA modification. Moreover, the recently gained m6A sites in human populations are clearly under positive selection and associated with traits or diseases. Our work provides a resource of human m6A profile for future studies of m6A functions, and suggests a role of m6A modification in human evolutionary adaptation and disease susceptibility.

Binaural processing deficit and cognitive impairment in Alzheimer's disease.

Speech comprehension in noisy environments depends on central auditory functions, which are vulnerable in Alzheimer's disease (AD). Binaural processing exploits two ear sounds to optimally process degraded sound information; its characteristics are poorly understood in AD. We studied behavioral and electrophysiological alterations in binaural processing among 121 participants (AD = 27; amnestic mild cognitive impairment [aMCI] = 33; subjective cognitive decline [SCD] = 30; cognitively normal [CN] = 31). We observed impairment of binaural processing in AD and aMCI, and detected a U-shaped curve change in phase synchrony (declining from CN to SCD and to aMCI, but increasing from aMCI to AD). This improvement in phase synchrony accompanying more severe cognitive stages could reflect neural adaptation for binaural processing. Moreover, increased phase synchrony is associated with worse memory during the stages when neural adaptation apparently occurs. These findings support a hypothesis that neural adaptation for binaural processing deficit may exacerbate cognitive impairment, which could help identify biomarkers and therapeutic targets in AD.

The landscape of miRNA editing in animals and its impact on miRNA biogenesis and targeting.

Adenosine-to-inosine (A-to-I) RNA editing regulates miRNA biogenesis and function. To date, fewer than 160 miRNA editing sites have been identified. Here, we present a quantitative atlas of miRNA A-to-I editing through the profiling of 201 pri-miRNA samples and 4694 mature miRNA samples in human, mouse, and Drosophila. We identified 4162 sites present in ∼80% of the pri-miRNAs and 574 sites in mature miRNAs. miRNA editing is prevalent in many tissue types in human. However, high-level editing is mostly found in neuronal tissues in mouse and Drosophila Interestingly, the edited miRNAs in neuronal and non-neuronal tissues in human gain two distinct sets of new targets, which are significantly associated with cognitive and organ developmental functions, respectively. Furthermore, we reveal that miRNA editing profoundly affects asymmetric strand selection. Altogether, these data provide insight into the impact of RNA editing on miRNA biology and suggest that miRNA editing has recently gained non-neuronal functions in human.

One-Pot Methylenation-Cyclization Employing Two Molecules of CO2 with Arylamines and Enaminones.

One-pot methylenation-cyclization employing two molecules of CO2 with enaminones and primary aromatic amines was discussed for the first time to access cyclized products. This 1,5,7-triazabicyclo[4.4.0]dec-5-ene and ZnCl2-catalyzed procedure was characterized by the selective conversion of two molecules of CO2 into methylene groups in a multicomponent cyclization reaction. According to the computational study and control experiments, the reaction might proceed through the generation of bis(silyl)acetal and condensation of arylamine and aza-Diels-Alder processes. Moreover, the resulting products will probably be potential organic building blocks with adjustable photophysical properties.

Synthesis of fused-tetrahydropyrimidines: one-pot methylenation-cyclization utilizing two molecules of CO2.

A methylenation-cyclization reaction, employing cyclic enaminones with primary aromatic amines and two molecules of CO2, furnishing fused-tetrahydropyrimidines, is discussed. In this Cs2CO3 and ZnI2 catalyzed one-pot two-step procedure, two molecules of CO2 were selectively converted to methylene groups. The multi-component reaction might proceed through the formation of bis(silyl)acetal which was followed by condensation and further aza-Diels-Alder reaction. Hydroquinazoline, hydrocyclopenta[d]pyrimidine and hydroindeno[1,2-d]pyrimidine derivatives could be prepared with CO2 as the C1 source, effectively.

Functional heritage: the evolution of chimeric RNA into a gene.

Once believed to be unique features of neoplasia, chimeric RNAs are now being discovered in normal physiology. We speculated that some chimeric RNAs may be functional precursors of genes, and that forming chimeric RNA at the transcriptional level may be a 'trial' mechanism before the functional element is fixed into the genome. Supporting this idea, we identified a chimeric RNA, HNRNPA1L2-SUGT1 (H-S), whose sequence is highly similar to that of a 'pseudogene' MRPS31P5. Sequence analysis revealed that MRPS31P5 transcript is more similar to H-S chimeric RNA than its 'parent' gene, MRPS31. Evolutionarily, H-S precedes MRPS31P5, as it can be detected bioinformatically and experimentally in marmosets, which do not yet possess MRPS31P5 in their genome. Conversely, H-S is minimally expressed in humans, while instead, MRPS31P5 is abundantly expressed. Silencing H-S in marmoset cells resulted in similar phenotype as silencing MRPS31P5 in human cells. In addition, whole transcriptome analysis and candidate downstream target validation revealed common signalling pathways shared by the two transcripts. Interestingly, H-S failed to rescue the phenotype caused by silencing MPRS31P5 in human and rhesus cells, whereas MRPS31P5 can at least partially rescue the phenotype caused by silencing H-S in marmoset cells, suggesting that MRPS31P5 may have further evolved into a distinct entity. Thus, multiple lines of evidence support that MRPS31P5 is not truly a pseudogene of MRPS31, but a likely functional descendent of H-S chimera. Instead being a gene fusion product, H-S is a product of cis-splicing between adjacent genes, while MRPS31P5 is likely produced by genome rearrangement.

Elevated CO2 concentration enhances drought resistance of soybean by regulating cell structure, cuticular wax synthesis, photosynthesis, and oxidative stress response.

The atmospheric [CO2] and the frequency and intensity of extreme weather events such as drought are increased, leading to uncertainty to soybean production. Elevated [CO2] (eCO2) partially mitigates the adverse effects of drought stress on crop growth and photosynthetic performance, but the mitigative mechanism is not well understood. In this study, soybean seedlings under drought stress simulated by PEG-6000 were grown in climate chambers with different [CO2] (400 µmol mol-1 and 700 µmol mol-1). The changes in anatomical structure, wax content, photosynthesis, and antioxidant enzyme were investigated by the analysis of physiology and transcriptome sequencing (RNA-seq). The results showed that eCO2 increased the thickness of mesophyll cells and decreased the thickness of epidermal cells accompanied by reduced stomatal conductance, thus reducing water loss in soybean grown under drought stress. Meanwhile, eCO2 up-regulated genes related to wax anabolism, thus producing more epidermal wax. Under drought stress, eCO2 increased net photosynthetic rate (PN), ribulose-1,5-bisphosphate carboxylase/oxygenase activity, and alerted the gene expressions in photosynthesis. The increased sucrose synthesis and decreased sucrose decomposition contributed to the progressive increase in the soluble saccharide contents under drought stress with or without eCO2. In addition, eCO2 increased the expressions of genes associated with peroxidase (POD) and proline (Pro), thus enhancing POD activity and Pro content and improving the drought resistance in soybean. Taken together, these findings deepen our understanding of the effects of eCO2 on alleviating drought stress in soybean and provide potential target genes for the genetic improvement of drought tolerance in soybean.

The role of ERp29/FOS/EMT pathway in excessive apoptosis of placental trophoblast cells in intrahepatic cholestasis of pregnancy.

BACKGROUND: Abnormal bile acid metabolism leading to changes in placental function during pregnancy. To determine whether endoplasmic reticulum protein 29 (ERp29) can mediate the pregnancy effects of cholestasis by altering the level of trophoblast cell apoptosis. METHODS: ERp29 in serum of 66 intrahepatic cholestasis of pregnancy (ICP) pregnant women and 74 healthy were detected by ELISA. Subcutaneous injection of ethinyl estradiol (E2) was used to induce ICP in pregnant rats. Taurocholic acid (TCA) was used to simulate the ICP environment, and TGF-ß1 was added to induce the epithelial mesenchymal transformation (EMT) process. The scratch, migration, and invasion test were used to detect the EMT process. ERp29 overexpression/knockdown vector were constructed and transfected to verify the role of ERp29 in the EMT process. Downstream gene was obtained through RNA-seq. RESULTS: Compared with the healthy pregnant women, the expression levels of ERp29 in serum of ICP pregnancy women were significantly increased (P < 0.001). ERp29 in the placenta tissue of the ICP pregnant rats increased significantly, and the level of apoptosis increased. The placental tissues of the ICP had high expression of E-cadherin and low expression of N-cadherin, snail1, vimentin. After HTR-8/SVneo cells were induced by TCA, EMT was inhibited, while the ERp29 increased. Cell and animal experiments showed that, knockdown of ERp29 reduced the inhibition of EMT, the ICP progress was alleviated. Overexpression of FOS salvaged the inhibitory effects of ERp29 on cell EMT. DISCUSSION: The high level of ERp29 in placental trophoblast cells reduced FOS mRNA levels, inhibited the EMT process and aggravated the occurrence and development of ICP.

Effect of Glycolipids Application Combined with Nitrogen Fertilizer Reduction on Maize Nitrogen Use Efficiency and Yield.

Microbial-driven N turnover is important in regulating N fertilizer use efficiency through the secretion of metabolites like glycolipids. Currently, our understanding of the potential of glycolipids to partially reduce N fertilizer use and the effects of glycolipids on crop yield and N use efficiency is still limited. Here, a three-year in situ field experiment was conducted with seven treatments: no fertilization (CK); chemical N, phosphorus and potassium (NPK); NPK plus glycolipids (N+PKT); and PK plus glycolipids with 10% (0.9 N+PKT), 20% (0.8 N+PKT), 30% (0.7 N+PKT), and 100% (PKT) N reduction. Compared with NPK, glycolipids with 0-20% N reduction did not significantly reduce maize yields, and also increased N uptake by 6.26-11.07%, but no significant changes in grain or straw N uptake. The N resorption efficiency under 0.9 N+PKT was significantly greater than that under NPK, while the apparent utilization rates of N fertilizer and partial factor productivity of N under 0.9 N+PKT were significantly greater than those under NPK. Although 0.9 N+PKT led to additional labor and input costs, compared with NPK, it had a greater net economic benefit. Our study demonstrates the potential for using glycolipids in agroecosystem management and provides theoretical support for optimizing fertilization strategies.

A protein-encoding CCDC7 circular RNA inhibits the progression of prostate cancer by up-regulating FLRT3.

Circular RNAs (circRNAs) are a family of endogenous RNAs that have become a focus of biological research in recent years. Emerging evidence has revealed that circRNAs exert biological functions by acting as transcriptional regulators, microRNA sponges, and binding partners with RNA-binding proteins. However, few studies have identified coding circRNAs, which may lead to a hidden repertoire of proteins. In this study, we unexpectedly discovered a protein-encoding circular RNA circCCDC7(15,16,17,18,19) while we were searching for prostate cancer related chimeric RNAs. circCCDC7(15,16,17,18,19) is derived from exon 19 back spliced to exon 15 of the CCDC7 gene. It is significantly downregulated in patients with high Gleason score. Prostate cancer patients with decreased circCCDC7(15,16,17,18,19) expression have a worse prognosis, while linear CCDC7 had no such association. Overexpressed circCCDC7(15,16,17,18,19) inhibited prostate cancer cell migration, invasion, and viability, supporting classification of circCCDC7(15,16,17,18,19) as a bona fide tumor suppressor gene. We provide evidence that its tumor suppressive activity is driven by the protein it encodes, and that circCCDC7(15,16,17,18,19) encodes a secretory protein. Consistently, conditioned media from circCCDC7(15,16,17,18,19) overexpressing cells has the same tumor suppressive activity. We further demonstrate that the tumor suppressive activity of circCCDC7(15,16,17,18,19) is at least partially mediated by FLRT3, whose expression also negatively correlates with Gleason score and clinical prognosis. In conclusion, circCCDC7(15,16,17,18,19) functions as a tumor suppressor in prostate cancer cells through the circCCDC7-180aa secretory protein it encodes, and is a promising therapeutic peptide for prostate cancer.

UBA1-CDK16 : A Sex-Specific Chimeric RNA and Its Role in Immune Sexual Dimorphism.

RNA processing mechanisms, such as alternative splicing and RNA editing, have been recognized as critical means to expand the transcriptome. Chimeric RNAs formed by intergenic splicing provide another potential layer of RNA diversification. By analyzing a large set of RNA-Seq data and validating results in over 1,200 blood samples, we identified UBA1-CDK16 , a female-specific chimeric transcript. Intriguingly, both parental genes, are expressed in males and females. Mechanistically, UBA1-CDK16 is produced by cis-splicing between the two adjacent X-linked genes, originating from the inactive X chromosome. A female-specific chromatin loop, formed between the junction sites, facilitates the alternative splicing of its readthrough precursor. This unique chimeric transcript exhibits evolutionary conservation, evolving to be female-specific from non-human primates to humans. Furthermore, our investigation reveals that UBA1-CDK16 is enriched in the myeloid lineage and plays a regulatory role in myeloid differentiation. Notably, female COVID-19 patients who tested negative for this chimeric transcript displayed higher counts of neutrophils, highlighting its potential role in disease pathogenesis. These findings support the notion that chimeric RNAs represent a new repertoire of transcripts that can be regulated independently from the parental genes, and a new class of RNA variance with potential implications in sexual dimorphism and immune responses.

Association of maternal mineral status with the risk of preterm birth: a retrospective cohort study.

Background: There has been a gradual increase in the proportion of preterm birth in China during the past several decades. Maternal malnutrition is a significant determinant for preterm birth. Nevertheless, comprehensive studies investigating serum mineral levels during pregnancy associated with preterm birth remain scarce. This study aims to assess the associations between maternal serum mineral levels and the risk of preterm birth. Methods: This retrospective cohort study of 18,048 pregnant women used data from a tertiary hospital in China from January 2016 to December 2022. Demographic data and serum mineral concentrations in the second and third trimesters of mothers were collected from the hospital information system. Analysis was performed using restricted cubic splines and logistic regression models. Results: The proportion of preterm birth in this study was 6.01%. Phosphorus [P for overall = 0.005; P for nonlinear = 0.490; OR (95%CI) = 1.11 (1.04, 1.18)] and chlorine [P for overall = 0.002; P for nonlinear = 0.058; OR (95%CI) = 1.11 (1.03, 1.19)] showed a significant positive correlation with preterm birth in a linear fashion. Furthermore, serum levels of potassium (P for nonlinear <0.001), sodium (P for nonlinear = 0.004), and magnesium (P for nonlinear <0.001) exhibited non-linear relationships with the risk of preterm birth. Conclusion: Serum levels of some minerals during pregnancy were associated with the risk of preterm birth among pregnant women. In addition to commonly recognized micronutrients such as folic acid, iron, and vitamin D, healthcare providers should also pay attention to the levels of these minerals during pregnancy.

Roles of non-coding RNAs in eye development and diseases.

The prevalence of ocular disorders is dramatically increasing worldwide, especially those that cause visual impairment and permanent loss of vision, including cataract, glaucoma, age-related macular degeneration, and diabetic retinopathy. Extensive evidence has shown that ncRNAs are key regulators in various biogenesis and biological functions, controlling gene expression related to histogenesis and cell differentiation in ocular tissues. Aberrant expression and function of ncRNA can lead to dysfunction of visual system and mediate progression of eye disorders. Here, we mainly offer an overview of the role of precise modulation of ncRNAs in eye development and function in patients with eye diseases. We also highlight the challenges and future perspectives in conducting ncRNA studies, focusing specifically on the role of ncRNAs that may hold expanded promise for their diagnostic and therapeutic applications in various eye diseases. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA in Disease and Development > RNA in Disease RNA in Disease and Development > RNA in Development.

Physiological and transcriptome analysis of response of soybean (Glycine max) to cadmium stress under elevated CO2 concentration.

The continuous accumulation of Cd has long-lasting detrimental effects on plant growth and food safety. Although elevated CO2 concentration (EC) has been reported to reduce Cd accumulation and toxicity in plants, evidence on the functions of elevated CO2 concentration and its mechanisms in the possible alleviation of Cd toxicity in soybean are limited. Here, we used physiological and biochemical methods together with transcriptomic comparison to explore the effects of EC on Cd-stressed soybean. Under Cd stress, EC significantly increased the weight of roots and leaves, promoted the accumulations of proline, soluble sugars, and flavonoid. In addition, the enhancement of GSH activity and GST gene expressions promoted Cd detoxification. These defensive mechanisms reduced the contents of Cd2+, MDA, and H2O2 in soybean leaves. The up-regulation of genes encoding phytochelatin synthase, MTPs, NRAMP, and vacuoles protein storage might play vital roles in the transportation and compartmentalization process of Cd. The MAPK and some transcription factors such as bHLH, AP2/ERF, and WRKY showed changed expressions and might be engaged in mediation of stress response. These findings provide a boarder view on the regulatory mechanism of EC on Cd stress and provide numerous potential target genes for future engineering of Cd-tolerant cultivars in soybean breeding programs under climate changes scenarios.

Identification and validation of potential hypoxia-related genes associated with coronary artery disease.

Introduction: Coronary artery disease (CAD) is one of the most life-threatening cardiovascular emergencies with high mortality and morbidity. Increasing evidence has demonstrated that the degree of hypoxia is closely associated with the development and survival outcomes of CAD patients. However, the role of hypoxia in CAD has not been elucidated. Methods: Based on the GSE113079 microarray dataset and the hypoxia-associated gene collection, differential analysis, machine learning, and validation of the screened hub genes were carried out. Results: In this study, 54 differentially expressed hypoxia-related genes (DE-HRGs), and then 4 hub signature genes (ADM, PPFIA4, FAM162A, and TPBG) were identified based on microarray datasets GSE113079 which including of 93 CAD patients and 48 healthy controls and hypoxia-related gene set. Then, 4 hub genes were also validated in other three CAD related microarray datasets. Through GO and KEGG pathway enrichment analyses, we found three upregulated hub genes (ADM, PPFIA4, TPBG) were strongly correlated with differentially expressed metabolic genes and all the 4 hub genes were mainly enriched in many immune-related biological processes and pathways in CAD. Additionally, 10 immune cell types were found significantly different between the CAD and control groups, especially CD8 T cells, which were apparently essential in cardiovascular disease by immune cell infiltration analysis. Furthermore, we compared the expression of 4 hub genes in 15 cell subtypes in CAD coronary lesions and found that ADM, FAM162A and TPBG were all expressed at higher levels in endothelial cells by single-cell sequencing analysis. Discussion: The study identified four hypoxia genes associated with coronary heart disease. The findings provide more insights into the hypoxia landscape and, potentially, the therapeutic targets of CAD.