Genome-wide RNA Sequencing Research Articles

Genome-Wide Analysis in the Study of the Fetal Growth Restriction Pathogenetics

Fetal growth restriction is a complication of pregnancy that defined as the inability of the fetus to realize its genetically determined growth potential. Despite the high social and medical significance of this problem the exact pathogenesis of fetal growth restriction is not known by now. Therefore, the analysis of the molecular genetics mechanisms of this pathology within the framework of approaches using modern high-performance technologies of next generation sequencing is of undoubted interest. In this review we focused on the analysis of data obtained in studies of the fetal growth restriction genetics component. The authors of these researches used next generation sequencing technologies and carried out whole transcriptome profiling. The results of the genes expression genome-wide analysis in placental tissue allow us to identify 1430 differentially expressed genes between fetal growth restriction and normal pregnancy, of which only 1% were found in at least two studies. These differentially expressed genes are involved in the Wnt/β-catenin signaling pathway that plays an important role in cell migration, neural pattern formation and organogenesis during embryonic development. Common genes are associated with both obstetric and gynecological diseases, as well as with various somatic conditions from the groups of neurodegenerative, cardiovascular diseases and mental disorders, which probably reflects their involvement in the development of postnatal consequences of fetal growth restriction. The results of our work don‘t point only to potential molecular mechanisms and key genes underlying fetal growth restriction, but also indicate the important role of gene-gene communications in this pathology implementation: about 30% of all identified differentially expressed genes products interact with each other within the same gene network. In general, genome-wide RNA sequencing combined with the analysis of protein-protein interactions represents a promising direction in research on the development and functioning of the placenta, as well as the identification of placental insufficiency diseases genetics mechanisms, including fetal growth restriction.

Transcriptional regulators ensuring specific gene expression and decision-making at high TGFβ doses.

TGFβ-signaling regulates cancer progression by controlling cell division, migration, and death. These outcomes are mediated by gene expression changes, but the mechanisms of decision-making toward specific fates remain unclear. Here, we combine SMAD transcription factor imaging, genome-wide RNA sequencing, and morphological assays to quantitatively link signaling, gene expression, and fate decisions in mammary epithelial cells. Fitting genome-wide kinetic models to our time-resolved data, we find that most of the TGFβ target genes can be explained as direct targets of SMAD transcription factors, whereas the remainder show signs of complex regulation, involving delayed regulation and strong amplification at high TGFβ doses. Knockdown experiments followed by global RNA sequencing revealed transcription factors interacting with SMADs in feedforward loops to control delayed and dose-discriminating target genes, thereby reinforcing the specific epithelial-to-mesenchymal transition at high TGFβ doses. We identified early repressors, preventing premature activation, and a late activator, boosting gene expression responses for a sufficiently strong TGFβ stimulus. Taken together, we present a global view of TGFβ-dependent gene regulation and describe specificity mechanisms reinforcing cellular decision-making.

EF-hand calcium sensor, EfhP, controls transcriptional regulation of iron uptake by calcium in Pseudomonas aeruginosa.

The human pathogen Pseudomonas aeruginosa (Pa) poses a major risk for a range of severe infections, particularly lung infections in patients suffering from cystic fibrosis (CF). As previously reported, the virulent behavior of this pathogen is enhanced by elevated levels of Ca2+ that are commonly present in CF nasal and lung fluids. In addition, a Ca2+-binding EF-hand protein, EfhP (PA4107), was partially characterized and shown to be critical for the Ca2+-regulated virulence in P. aeruginosa. Here, we describe the rapid (10 min, 60 min), and adaptive (12 h) transcriptional responses of PAO1 to elevated Ca2+ detected by genome-wide RNA sequencing and show that efhP deletion significantly hindered both rapid and adaptive Ca2+ regulation. The most differentially regulated genes included multiple Fe sequestering mechanisms, a large number of extracytoplasmic function sigma factors (ECFσ), and several virulence factors, such as the production of pyocins. The Ca2+ regulation of Fe uptake was also observed in CF clinical isolates and appeared to involve the global regulator Fur. In addition, we showed that the efhP transcription is controlled by Ca2+ and Fe, and this regulation required a Ca2+-dependent two-component regulatory system CarSR. Furthermore, the efhP expression is significantly increased in CF clinical isolates and upon pathogen internalization into epithelial cells. Overall, the results established for the first time that Ca2+ controls Fe sequestering mechanisms in P. aeruginosa and that EfhP plays a key role in the regulatory interconnectedness between Ca2+ and Fe signaling pathways, the two distinct and important signaling pathways that guide the pathogen's adaptation to the host.IMPORTANCEPseudomonas aeruginosa (Pa) poses a major risk for severe infections, particularly in patients suffering from cystic fibrosis (CF). For the first time, kinetic RNA sequencing analysis identified Pa rapid and adaptive transcriptional responses to Ca2+ levels consistent with those present in CF respiratory fluids. The most highly upregulated processes include iron sequestering, iron starvation sigma factors, and self-lysis factors pyocins. An EF-hand Ca2+ sensor, EfhP, is required for at least 1/3 of the Ca2+ response, including the majority of the iron uptake mechanisms and the production of pyocins. Transcription of efhP itself is regulated by Ca2+ and Fe, and increases during interactions with host epithelial cells, suggesting the protein's important role in Pa infections. The findings establish the regulatory interconnectedness between Ca2+ and iron signaling pathways that shape Pa transcriptional responses. Therefore, understanding Pa's transcriptional response to Ca2+ and associated regulatory mechanisms will serve in the development of future therapeutics targeting Pa's dangerous infections.

RUNX2 is stabilised by TAZ and drives pulmonary artery calcification and lung vascular remodelling in pulmonary hypertension due to left heart disease.

Calcification is common in chronic vascular disease, yet its role in pulmonary hypertension due to left heart disease is unknown. Here, we probed for the role of runt-related transcription factor-2 (RUNX2), a master transcription factor in osteogenesis, and its regulation by the HIPPO pathway transcriptional coactivator with PDZ-binding motif (TAZ) in the osteogenic reprogramming of pulmonary artery smooth muscle cells and vascular calcification in patients with pulmonary hypertension due to left heart disease. We similarly examined its role using an established rat model of pulmonary hypertension due to left heart disease induced by supracoronary aortic banding. Pulmonary artery samples were collected from patients and rats with pulmonary hypertension due to left heart disease. Genome-wide RNA sequencing was performed, and pulmonary artery calcification assessed. Osteogenic signalling via TAZ and RUNX2 was delineated by protein biochemistry. In vivo, the therapeutic potential of RUNX2 or TAZ inhibition by CADD522 or verteporfin was tested in the rat model. Gene ontology term analysis identified significant enrichment in ossification and osteoblast differentiation genes, including RUNX2, in pulmonary arteries of patients and lungs of rats with pulmonary hypertension due to left heart disease. Pulmonary artery calcification was evident in both patients and rats. Both TAZ and RUNX2 were upregulated and activated in pulmonary artery smooth muscle cells of patients and rats. Co-immunoprecipitation revealed a direct interaction of RUNX2 with TAZ in pulmonary artery smooth muscle cells. TAZ inhibition or knockdown decreased RUNX2 abundance due to accelerated RUNX2 protein degradation rather than reduced de novo synthesis. Inhibition of either TAZ or RUNX2 attenuated pulmonary artery calcification, distal lung vascular remodelling and pulmonary hypertension development in the rat model. Pulmonary hypertension due to left heart disease is associated with pulmonary artery calcification that is driven by TAZ-dependent stabilisation of RUNX2, causing osteogenic reprogramming of pulmonary artery smooth muscle cells. The TAZ-RUNX2 axis may present a therapeutic target in pulmonary hypertension due to left heart disease.

Multi-scalar data integration decoding risk genes for chronic kidney disease

BackgroundChronic Kidney Disease (CKD) impacts over 10% of the global population, and recent advancements in high-throughput analytical technologies are uncovering the complex physiology underlying this condition. By integrating Genome-Wide Association Studies (GWAS), RNA sequencing (RNA-seq/RNA array), and single-cell RNA sequencing (scRNA-seq) data, our study aimed to explore the genes and cell types relevant to CKD traits.MethodsGWAS summary data for end-stage renal failure (ESRD) and decreased eGFR (CKD) with or without diabetes and (micro)proteinuria were obtained from the GWAS Catalog and the UK Biobank (UKB) database. Two gene Expression Omnibus (GEO) transcriptome datasets were used to establish glomerular and tubular gene expression differences between CKD patients and healthy individuals. Two scRNA-seq datasets were utilized to obtain the expression of key genes at the single-cell level. The expression profile, differentially expressed genes (DEGs), gene-gene interaction, and pathway enrichment were analysed for these CKD risk genes.ResultsA total of 779 distinct SNPs were identified from GWAS across different CKD traits, involving 681 genes. While many of these risk genes are specific to the CKD traits of renal failure, decreased eGFR, and (micro)proteinuria, they share common pathways, including extracellular matrix (ECM). ECM modeling was enriched in upregulated glomerular and tubular DEGs from CKD kidneys compared to healthy controls, with the expression of relevant collagen genes, such as COL1A2, prevalent in fibroblasts/myofibroblasts. Additionally, immune responses, including T cell differentiation, were dysregulated in CKD kidneys. The late podocyte signature gene THSD7A was enriched in podocytes but downregulated in CKD. We also highlighted that the regulated risk genes of CKD are mainly expressed in tubular cells and immune cells in the kidney.ConclusionsOur integrated analysis highlight the genes, pathways, and relevant cell types associational with the pathogenesis of kidney traits, as a basis for further mechanistic studies to understand the pathogenesis of CKD.

The role of DNA methylation in chondrogenesis of human iPSCs as a stable marker of cartilage quality

BackgroundLack of insight into factors that determine purity and quality of human iPSC (hiPSC)-derived neo-cartilage precludes applications of this powerful technology toward regenerative solutions in the clinical setting. Here, we set out to generate methylome-wide landscapes of hiPSC-derived neo-cartilages from different tissues-of-origin and integrated transcriptome-wide data to identify dissimilarities in set points of methylation with associated transcription and the respective pathways in which these genes act.MethodsWe applied in vitro chondrogenesis using hiPSCs generated from two different tissue sources: skin fibroblasts and articular cartilage. Upon differentiation toward chondrocytes, these are referred to as hFiCs and hCiC, respectively. Genome-wide DNA methylation and RNA sequencing datasets were generated of the hiPSC-derived neo-cartilages, and the epigenetically regulated transcriptome was compared to that of neo-cartilage deposited by human primary articular cartilage (hPAC).ResultsMethylome-wide landscapes of neo-cartilages of hiPSCs reprogrammed from two different somatic tissues were 85% similar to that of hPACs. By integration of transcriptome-wide data, differences in transcriptionally active CpGs between hCiC relative to hPAC were prioritized. Among the CpG-gene pairs lower expressed in hCiCs relative to hPACs, we identified genes such as MGP, GDF5, and CHAD enriched in closely related pathways and involved in cartilage development that likely mark phenotypic differences in chondrocyte states. Vice versa, among the CpG-gene pairs higher expressed, we identified genes such as KIF1A or NKX2-2 enriched in neurogenic pathways and likely reflecting off target differentiation.ConclusionsWe did not find significant variation between the neo-cartilages derived from hiPSCs of different tissue sources, suggesting that application of a robust differentiation protocol such as we applied here is more important as compared to the epigenetic memory of the cells of origin. Results of our study could be further exploited to improve quality, purity, and maturity of hiPSC-derived neo-cartilage matrix, ultimately to realize introduction of sustainable, hiPSC-derived neo-cartilage implantation into clinical practice.

Circulating microRNAs and isomiRs as biomarkers for the initial insult and epileptogenesis in four experimental epilepsy models: The EPITARGET study.

Structural epilepsies can manifest months or years after the occurrence of an initial epileptogenic insult, making them amenable for secondary prevention. However, development of preventive treatments has been challenged by a lack of biomarkers for identifying the subset of individuals with the highest risk of epilepsy after the epileptogenic insult. Four different rat models of epileptogenesis were investigated to identify differentially expressed circulating microRNA (miRNA) and isomiR profiles as biomarkers for epileptogenesis. Plasma samples were collected on day 2 and day 9 during the latency period from animals that did or did not develop epilepsy during long-term video-electroencephalographic monitoring. miRNAs and isomiRs were identified and measured in an unsupervised manner, using a genome-wide small RNA sequencing platform. Receiver operating characteristic analysis was performed to determine the performance of putative biomarkers. Two days after an epileptogenic insult, alterations in the levels of several plasma miRNAs and isomiRs predicted epileptogenesis in a model-specific manner. One miRNA, miR-3085, showed good sensitivity (but low specificity) as a prognostic biomarker for epileptogenesis in all four models (area under the curve = .729, sensitivity = 83%, specificity = 64%, p < .05). Identified plasma miRNAs and isomiRs are mostly etiology-specific rather than common prognostic biomarkers of epileptogenesis. These data imply that in preclinical and clinical studies, it may be necessary to identify specific biomarkers for different epilepsy etiologies. Importantly, circulating miRNAs like miR-3085 with high negative predictive value for epileptogenesis in different etiologies could be useful candidates for initial screening purposes of epileptogenesis risk.

Genetic regulation of microRNAs in the older adult brain and their contribution to neuropsychiatric conditions.

MicroRNAs are essential post-transcriptional regulators of gene expression and involved in many biological processes; however, our understanding of their genetic regulation and role in brain illnesses is limited. Here, we mapped brain microRNA expression quantitative trait loci (miR-QTLs) using genome-wide small RNA sequencing profiles from dorsolateral prefrontal cortex (dlPFC) samples of 604 older adult donors of European ancestry. miR-QTLs were identified for 224 miRNAs (48% of 470 tested miRNAs) at false discovery rate < 1%. We found that miR-QTLs were enriched in brain promoters and enhancers, and that intragenic miRNAs often did not share QTLs with their host gene. Additionally, we integrated the brain miR-QTLs with results from 16 GWAS of psychiatric and neurodegenerative diseases using multiple independent integration approaches and identified four miRNAs that contribute to the pathogenesis of bipolar disorder, major depression, post-traumatic stress disorder, schizophrenia, and Parkinson's disease. This study provides novel insights into the contribution of miRNAs to the complex biological networks that link genetic variation to disease.

The methyltransferase SETD3 regulates mRNA alternative splicing through interacting with hnRNPK

The methyltransferase SETD3 is an enzyme essential for catalyzing histidine-73 methylation on β-Actin, thereby promoting its polymerization and regulating muscle contraction. Although increasing evidence suggests that SETD3 is involved in multiple physiological or pathological events, its biological functions remain incompletely understood. In this study, we utilize in situ proximity labeling combined with mass spectrometry analysis to detect potential interacting partners of SETD3. Unexpectedly, we find that many splicing factors are associated with SETD3. Genome-wide RNA sequencing reveals that SETD3 regulates pre-mRNA splicing events, predominantly influencing exon skipping. Biochemical and bioinformatic analyses suggest that SETD3 interacts with hnRNPK, and they collaboratively regulate exon skipping in a common subset of genes. Functionally, we demonstrate that SETD3 and hnRNPK are required for retention of exon 7 skipping in the FNIP1 gene. This promotes FNIP1-mediated nuclear translocation of the transcription factor TFEB and the subsequent induction of lysosomal and mitochondrial biogenesis. Overall, this study uncovers a novel function of SETD3 in modulating mRNA exon splicing.

Paecilomyces variotii extracts promote growth and alleviate nutrient deficiency symptoms in apple

It is essential to find environment-friendly agrochemicals to cope with the problems of nutrient imbalance, fruit quality decline, and physiological disorders during apple fruit production, which is beneficial for improving the quality and yield of apple. A natural extract from Paecilomyces variotii (ZNC), an endophytic fungus, has been used widely to enhance crop performance. However, an understanding of the mechanism underlying ZNC-triggered growth and alleviation of nutrient deficiency-associated symptoms in apples remains elusive. Here, the photosynthesis, leaf growth, and fruit quality were enhanced by adding ZNC. In addition, ZNC relieved nutrient deficiency-related symptoms promoted the differentiation of root morphology and vitality, and reduced the accumulation of osmoprotectants and reactive oxygen species, thereby promoting growth under normal and nutrient-deficient conditions. Finally, genome-wide RNA sequencing reveals the ZNC-regulated mechanisms involved in hormone and metal ion pathways. Our study reveals the role of ZNC in promoting growth and improving the quality of apple, providing a new direction for reducing the use of chemical fertilizers and pesticides.

Integrated Genome-Wide Analysis of DNA Methylation and Gene Expression in the Hippocampi of 5xFAD Alzheimer's Disease Mouse Model.

DNA methylation forms 5-methylcytosine and its regulation in the hippocampus is critical for learning and memory. Indeed, dysregulation of DNA methylation is associated with neurological diseases. Alzheimer's disease (AD) is the predominant of dementia and a neurodegenerative disorder. We examined the learning and memory function in 3- and 9-month-old wild-type and 5xfamiliar Alzheimer's disease (5xFAD) transgenic mice by performing the object recognition memory and Y-maze tests, and identified the hippocampal amyloid beta burden. To investigate the epigenetically regulated genes involved in the development or neuropathology of AD, we performed genome-wide DNA methylation sequencing and RNA sequencing analyses in the hippocampus of 9-month-old wild-type and 5xFAD tg mice. To validate the genes inversely regulated by epigenetics, we confirmed their methylation status and mRNA levels. At 9 months of age, 5xFAD tg mice showed significant cognitive impairment and amyloid-beta plaques in the hippocampus. DNA methylation sequencing identified a total of 13,777 differentially methylated regions, including 4484 of hyper- and 9293 of hypomethylated regions, that are associated with several gene ontology (GO) terms including 'nervous system development' and 'axon guidance'. In RNA sequencing analysis, we confirmed a total of 101 differentially expressed genes, including 52 up- and 49 downregulated genes, associated with GO functions such as 'positive regulation of synaptic transmission, glutamatergic' and 'actin filament organization'. Through further integrated analysis of DNA methylation and RNA sequencing, three epigenetically regulated genes were selected: thymus cell antigen 1, theta (Thy1), myosin VI (Myo6), and filamin A-interacting protein 1-like (Filip1l). The methylation level of Thy1 decreased and its mRNA levels increased, whereas that of Myo6 and Filip1l increased and their mRNA levels decreased. The common functions of these three genes may be associated with the neural cytoskeleton and synaptic plasticity. We suggest that the candidate genes epigenetically play a role in AD-associated neuropathology (i.e., amyloid-beta plaques) and memory deficit by influencing neural structure and synaptic plasticity. Furthermore, counteracting dysregulated epigenetic changes may delay or ameliorate AD onset or symptoms.

Smoking-informed methylation and expression QTLs in human brain and colocalization with smoking-associated genetic loci

Smoking is a leading cause of preventable morbidity and mortality. Smoking is heritable, and genome-wide association studies (GWASs) of smoking behaviors have identified hundreds of significant loci. Most GWAS-identified variants are noncoding with unknown neurobiological effects. We used genome-wide genotype, DNA methylation, and RNA sequencing data in postmortem human nucleus accumbens (NAc) to identify cis-methylation/expression quantitative trait loci (meQTLs/eQTLs), investigate variant-by-cigarette smoking interactions across the genome, and overlay QTL evidence at smoking GWAS-identified loci to evaluate their regulatory potential. Active smokers (N = 52) and nonsmokers (N = 171) were defined based on cotinine biomarker levels and next-of-kin reporting. We simultaneously tested variant and variant-by-smoking interaction effects on methylation and expression, separately, adjusting for biological and technical covariates and correcting for multiple testing using a two-stage procedure. We found >2 million significant meQTL variants (padj < 0.05) corresponding to 41,695 unique CpGs. Results were largely driven by main effects, and five meQTLs, mapping to NUDT12, FAM53B, RNF39, and ADRA1B, showed a significant interaction with smoking. We found 57,683 significant eQTL variants for 958 unique eGenes (padj < 0.05) and no smoking interactions. Colocalization analyses identified loci with smoking-associated GWAS variants that overlapped meQTLs/eQTLs, suggesting that these heritable factors may influence smoking behaviors through functional effects on methylation/expression. One locus containing MUSTN1 and ITIH4 colocalized across all data types (GWAS, meQTL, and eQTL). In this first genome-wide meQTL map in the human NAc, the enriched overlap with smoking GWAS-identified genetic loci provides evidence that gene regulation in the brain helps explain the neurobiology of smoking behaviors.

MEK/ERK signaling drives the transdifferentiation of supporting cells into functional hair cells by modulating the Notch pathway.

Loss of cochlear hair cells (HCs) leads to permanent hearing loss in mammals, and regenerative medicine is regarded as an ideal strategy for hearing recovery. Limited genetic and pharmaceutical approaches for HC regeneration have been established, and the existing strategies cannot achieve recovery of auditory function. A promising target to promote HC regeneration is MEK/ERK signaling because dynamic shifts in its activity during the critical stages of inner ear development have been observed. Here, we first showed that MEK/ERK signaling is activated specifically in supporting cells (SCs) after aminoglycoside-induced HC injury. We then selected 4 MEK/ERK signaling inhibitors, and PD0325901 (PD03) was found to induce the transdifferentiation of functional supernumerary HCs from SCs in the neonatal mammalian cochlear epithelium. We next found that PD03 facilitated the generation of HCs in inner ear organoids. Through genome-wide high-throughput RNA sequencing and verification, we found that the Notch pathway is the downstream target of MEK/ERK signaling. Importantly, delivery of PD03 into the inner ear induced mild HC regeneration in vivo. Our study thus reveals the importance of MEK/ERK signaling in cell fate determination and suggests that PD03 might serve as a new approach for HC regeneration.

Assessing the ability of dihydrosterculic acid to activate PPARα

Peroxisome proliferator-activated receptor α (PPARα) is a ligand-activated transcription factor associated with the expression of specific target genes involved in fatty acid β-oxidation. Activation of PPARα could potentially be utilized to treat diseases associated with the accumulation of hepatic lipids such as non-alcoholic fatty liver disease (NAFLD), which currently impacts approximately 25% of the world’s population. At present, a class of synthetic ligands called fibrates can be prescribed to help manage NAFLD but are weak agonists for PPARα and therefore have limited effcacy. Consequently, alternative ligands for PPARα must be established and no natural endogenous ligands have been conclusively identified. This study expands on our previous work and aims to identify the direct effects of dihydrosterculic acid (DHSA), a cyclopropyl fatty acid found in cottonseed oil (CSO), on PPARα activity. To identify downstream effects of DHSA, male mice were placed on a CSO- (DHSA-containing) or isocaloric oil-enriched diet for 6 weeks and tissues were sent for genome-wide RNA sequencing analysis. CSO-fed mice demonstrated increased PPARα gene expression in the liver and concomitant decreased triglyceride in plasma (-75.1±11.5 milligrams/deciliter (mg/dL); p&lt;0.01) and liver (-86.7±71.7 mg/gram (g) protein; p=0.1) and decreased free fatty acid in liver (-16.5±3.4 micromolar (μM)/g protein; p&lt;0.01) relative to the oil-enriched control diet. Confirmational tissue analysis of liver samples show an increase in PPARα target gene expression in the DHSA containing groups compared to control diet groups, indicating that DHSA could potentially serve as an endogenous ligand for PPARα to increase its transcriptional activity to increase fat oxidation. To further test this hypothesis, female PPARα knockout mice were fed a DHSA containing diet or control diet for 4 weeks (n=8). Following the 4-week period, the chow-fed wild type females gained more weight (1.3±0.6 g; p&lt;0.05) than the PPARα knockout mice. Otherwise, there were no significant differences in body weight or food intake between diets or genotypes. Additional in vitro studies show increased PPARα target gene expression following treatments with DHSA, specifically for CPT1a (2.7±0.2 fold increase; p&lt;0.01), ACOX1 (1.3±0.1 fold increase; p&lt;0.05), and ACADVL (1.6±0.1 fold increase, p&lt;0.01), which encode key enzymes in lipid oxidation pathways. SCD1, the rate limiting enzyme in lipogenesis, was also downregulated (0.5±0.03 fold decrease, p&lt;0.01) in response to a 24-hour treatment with DHSA indicating a reduction in endogenous lipogenesis. From these results, we believe that DHSA could function through a PPARα-dependent mechanism. Further analysis of the tissues collected from the PPARα knockout mouse model will aid in determining the mechanism of DHSA as well as identifying its potential for the treatment of excess lipid accumulation. This project is funded by Cotton Incorporated. This is the full abstract presented at the American Physiology Summit 2024 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

Spatiotemporal ATF3 Expression Determines VSMC Fate in Abdominal Aortic Aneurysm.

Abdominal aortic aneurysm (AAA) is a catastrophic disease with little effective therapy, likely due to the limited understanding of the mechanisms underlying AAA development and progression. ATF3 (activating transcription factor 3) has been increasingly recognized as a key regulator of cardiovascular diseases. However, the role of ATF3 in AAA development and progression remains elusive. Genome-wide RNA sequencing analysis was performed on the aorta isolated from saline or Ang II (angiotensin II)-induced AAA mice, and ATF3 was identified as the potential key gene for AAA development. To examine the role of ATF3 in AAA development, vascular smooth muscle cell-specific ATF3 knockdown or overexpressed mice by recombinant adeno-associated virus serotype 9 vectors carrying ATF3, or shRNA-ATF3 with SM22α (smooth muscle protein 22-α) promoter were used in Ang II-induced AAA mice. In human and murine vascular smooth muscle cells, gain or loss of function experiments were performed to investigate the role of ATF3 in vascular smooth muscle cell proliferation and apoptosis. In both Ang II-induced AAA mice and patients with AAA, the expression of ATF3 was reduced in aneurysm tissues but increased in aortic lesion tissues. The deficiency of ATF3 in vascular smooth muscle cell promoted AAA formation in Ang II-induced AAA mice. PDGFRB (platelet-derived growth factor receptor β) was identified as the target of ATF3, which mediated vascular smooth muscle cell proliferation in response to TNF-alpha (tumor necrosis factor-α) at the early stage of AAA. ATF3 suppressed the mitochondria-dependent apoptosis at the advanced stage by upregulating its direct target BCL2. Our chromatin immunoprecipitation results also demonstrated that the recruitment of NFκB1 and P300/BAF/H3K27ac complex to the ATF3 promoter induces ATF3 transcription via enhancer activation. NFKB1 inhibitor (andrographolide) inhibits the expression of ATF3 by blocking the recruiters NFKB1 and ATF3-enhancer to the ATF3-promoter region, ultimately leading to AAA development. Our results demonstrate a previously unrecognized role of ATF3 in AAA development and progression, and ATF3 may serve as a novel therapeutic and prognostic marker for AAA.

Identification of potential candidate genes associated with milk protein differences in Holstein cows: a meta-analysis integrating GWAS and RNA-Seq transcriptome

Despite the identification of candidate genes influencing milk protein, the connections between genes and regulatory pathways remains elusive. This study aimed integrate findings from genome-wide association studies (GWAS) and RNA sequencing (RNA-Seq) through meta-analysis to pinpoint single nucleotide polymorphisms (SNPs) and genes responsible for high and low protein yield in cows. Previous GWAS and RNA-Seq analyses had identified 663 SNPs and 1106 genes ( P &lt; 0.05). Twenty SNPs from GWAS, 10 genes from RNA-Seq, and 49 SNP/gene associations from both datasets, were identified using meta-analysis. Meta-analysis validated several SNPs previously identified through GWAS, such as rs135549651 ( P = 2.6 × 10−256), rs109146371 ( P = 3.1 × 10−208), rs109350371 ( P = 4.0 × 10−207), and rs109774038 ( P = 8.6 × 10−587). Genes identified in RNA-Seq experiments, including NR4A1 ( P = 3.2 × 10−7), ATF3 ( P = 9.6 × 10−7), CDH16 ( P = 9.9 × 10−7), VEGFA ( P = 1.0 × 10−6), and SAA3 ( P = 7.3 × 10−11), were confirmed. The combined GWAS and RNA-Seq datasets highlighted CCND2 ( P = 8.9 × 10−111), MAPK15 ( P = 1.3 × 10−151), and CPSF1 ( P = 1.2 × 10−306) as the most significant genes. Additionally, significant gene ontology (GO) terms, including ionizing radiation ( P = 1.5 × 10−4), nuclear pore cytoplasmic filaments ( P = 9.4 × 10−5), and phenylalanine 4-monooxygenase activity ( P = 1.4 × 10−5), were identified. In conclusion, the integration of GWAS and RNA-Seq, coupled with GO enrichment, allowed identification of candidate SNPs and genes with higher accuracy. These findings improve our knowledge about genomic architecture of milk protein and enhance evaluation of Holstein cows.

Type I interferon blockade with anifrolumab in patients with systemic lupus erythematosus modulates key immunopathological pathways in a gene expression and proteomic analysis of two phase 3 trials

IntroductionAnifrolumab is a type I interferon (IFN) receptor 1 (IFNAR1) blocking antibody approved for treating patients with systemic lupus erythematosus (SLE). Here, we investigated the immunomodulatory mechanisms of anifrolumab using...

Abstract 6389: An exosome-based liquid biopsy powered by machine learning predicts progression-free and overall survival before commencing first-line EGFR inhibitors for metastatic colorectal cancer: The EXONERATE Study

Abstract Introduction: We developed a liquid biopsy assay based on exosomal and cell-free microRNA (exo- and cf-miRNA, respectively) for patients with chemotherapy-naïve, metastatic, unresectable, RAS and BRAF wild-type colorectal cancer (mCRC), to predict progression-free (PFS) and overall survival (OS) before the commencement of first-line treatment with EGFR inhibitors (EXONERATE). Methods: Patients with mCRC were consecutively enrolled in two open-label, nationwide clinical trials and received first-line cetuximab or panitumumab. Cf- and exo-miRNA expression was quantified using RT-qPCR in plasma samples obtained before exposure to EGFR inhibitors. The prespecified primary endpoint was 12-month PFS, hierarchically tested in left-sided mCRC, right-sided mCRC, and the entire cohort. As secondary endpoints, the EXONERATE status was correlated with PFS and OS hazard ratios (HR). Results: After excluding those with RAS/BRAF mutations and unavailable samples, the complete analysis dataset comprised 120 patients (94 with left-sided and 26 with right-sided mCRC; 87 received cetuximab and 33 panitumumab). As expected, right-sided mCRC patients showed inferior PFS and OS outcomes (HR 1.80 [CI95% = 1.13 - 2.80] and 1.80 [CI95% = 1.00 - 3.27], respectively). Genome-wide small RNA sequencing on 20 patient samples identified 12 cf- and 14 exo-miRNA candidates, then narrowed to 8 and 9 based on best performance. We employed machine learning algorithms to develop two liquid biopsy assays, one based on cf- and the other on exo-miRNAs, predicting 12mPFS with AUC values of 0.84 and 0.81, respectively. These were then combined to establish the EXONERATE assay for left-sided mCRC (AUC = 0.89), successfully validated in right-sided mCRC patients (AUC= 0.84). Among left-sided mCRC patients, 35% were EXONERATE-high, exhibiting shorter median PFS (mPFS = 9.5 vs. 18.5 months, p&amp;lt;.001) and a higher mortality rate (HR 0.21 [CI95% = 0.10 - 0.42]). In the right-sided mCRC cohort, the EXONERATE-high panel robustly identified patients with shorter mPFS (8.6 vs. 23.7 months, p=.0042) and predicted 12mPFS with high accuracy (sensitivity 100%, NPV 81%, PPV 100%). Although in univariate analysis both the EXONERATE panel and sidedness showed associations with PFS and OS, in multivariate analysis only the EXONERATE panel could fully account for the differences in PFS and OS (HR 0.22 and HR 0.19, respectively). Hence, we employed the EXONERATE panel agnostically of tumor sidedness and observed that EXONERATE-high patients demonstrated poorer survival outcomes in mPFS (8.9 vs. 19.6 months, HR 0.27 [CI95% = 0.17 - 0.41]) and mOS (22.6 vs. &amp;gt;60 months, HR 0.18 [CI95% = 0.09 - 0.37]). Conclusion: The EXONERATE panel robustly predicted PFS and OS outcomes in both right- and left-sided mCRC patients before receiving either EGFR inhibitor. Citation Format: Xu Caiming, Alessandro Mannucci, Roy Souvick, Esposito Francis, Oliveres Helena, Alonso-Orduña Vicente, Escudero Pilar, Fernández-Martos Carlos, Salud Antonieta, Gallego Javier, Rodriguez Jose Ramon, Martín-Richard Marta M. Marta, Fernández-Plana Julen, Manzano Hermini, Aparicio Jorge, Feliu Jaime, Maurel Joan, Ajay Goel. An exosome-based liquid biopsy powered by machine learning predicts progression-free and overall survival before commencing first-line EGFR inhibitors for metastatic colorectal cancer: The EXONERATE Study [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6389.

Genetic background of hematological parameters in Holstein cattle based on genome-wide association and RNA sequencing analyses

Hematological parameters refer to the assessment of changes in the number and distribution of blood cells, including leukocytes (LES), erythrocytes (ERS), and platelets (PLS), which are essential for the early diagnosis of hematological system disorders and other systemic diseases in livestock. In this context, the primary objectives of this study were to investigate the genomic background of 19 hematological parameters in Holstein cattle, focusing on LES, ERS, and PLS blood components. Genetic and phenotypic (co)variances of hematological parameters were calculated based on the Average Information Restricted Maximum Likelihood (AIREML) method and 1,610 genotyped individuals and 5,499 hematological parameter records from 4,543 cows. Furthermore, we assessed the genetic relationship between these hematological parameters and other economically important traits in dairy cattle breeding programs. We also carried out genome-wide association studies and candidate gene analyses. Blood samples from 21 primiparous cows were used to identify candidate genes further through RNA sequencing (RNA-seq) analyses. Hematological parameters generally exhibited low-to-moderate heritabilities ranging from 0.01 to 0.29, with genetic correlations between them ranging from -0.88 ± 0.09 (between mononuclear cell ratio and lymphocyte cell ratio) to 0.99 ± 0.01 (between white blood cell count and granulocyte cell count). Furthermore, low to moderate approximate genetic correlations between hematological parameters with one longevity, 4 fertility, and 5 health traits were observed. One-hundred-and-99 significant single nucleotide polymorphisms (SNP) located primarily on the Bos taurus autosomes (BTA) BTA4, BTA6, and BTA8 were associated with 16 hematological parameters. Based on the RNA-seq analyses, 6,687 genes were significantly downregulated and 4,119 genes were upregulated when comparing 2 groups of cows with high and low phenotypic values. By integrating genome-wide association studies (GWAS), RNA-seq, and previously published results, the main candidate genes associated with hematological parameters in Holstein cattle were ACRBP, ADAMTS3, CANT1, CCM2L, CNN3, CPLANE1, GPAT3, GRIP2, PLAGL2, RTL6, SOX4, WDFY3, and ZNF614. Hematological parameters are heritable and moderately to highly genetically correlated among themselves. The large number of candidate genes identified based on GWAS and RNA-seq indicate the polygenic nature and complex genetic determinism of hematological parameters in Holstein cattle.

CHD7 and SOX2 act in a common gene regulatory network during mammalian semicircular canal and cochlear development

Inner ear morphogenesis requires tightly regulated epigenetic and transcriptional control of gene expression. CHD7, an ATP-dependent chromodomain helicase DNA-binding protein, and SOX2, an SRY-related HMG box pioneer transcription factor, are known to contribute to vestibular and auditory system development, but their genetic interactions in the ear have not been explored. Here, we analyzed inner ear development and the transcriptional regulatory landscapes in mice with variable dosages of Chd7 and/or Sox2. We show that combined haploinsufficiency for Chd7 and Sox2 results in reduced otic cell proliferation, severe malformations of semicircular canals, and shortened cochleae with ectopic hair cells. Examination of mice with conditional, inducible Chd7 loss by Sox2CreER reveals a critical period (~E9.5) of susceptibility in the inner ear to combined Chd7 and Sox2 loss. Data from genome-wide RNA-sequencing and CUT&Tag studies in the otocyst show that CHD7 regulates Sox2 expression and acts early in a gene regulatory network to control expression of key otic patterning genes, including Pax2 and Otx2. CHD7 and SOX2 directly bind independently and cooperatively at transcription start sites and enhancers to regulate otic progenitor cell gene expression. Together, our findings reveal essential roles for Chd7 and Sox2 in early inner ear development and may be applicable for syndromic and other forms of hearing or balance disorders.