Genome-wide Transcriptome Analysis Research Articles

Integration of CTCF loops, methylome, and transcriptome in differentiating LUHMES as a model for imprinting dynamics of the 15q11-q13 locus in human neurons.

Human cell line models, including the neuronal precursor line LUHMES, are important for investigating developmental transcriptional dynamics within imprinted regions, particularly the 15q11-q13 Angelman (AS) and Prader-Willi (PWS) syndrome locus. AS results from loss of maternal UBE3A in neurons, where the paternal allele is silenced by a convergent antisense transcript UBE3A-ATS, a lncRNA that terminates at PWAR1 in non-neurons. qRT-PCR analysis confirmed the exclusive and progressive increase in UBE3A-ATS in differentiating LUHMES neurons, validating their use for studying UBE3A silencing. Genome-wide transcriptome analyses revealed changes to 11 834 genes during neuronal differentiation, including the upregulation of most genes within the 15q11-q13 locus. To identify dynamic changes in chromatin loops linked to transcriptional activity, we performed a HiChIP validated by 4C, which identified two neuron-specific CTCF loops between MAGEL2-SNRPN and PWAR1-UBE3A. To determine if allele-specific differentially methylated regions (DMR) may be associated with CTCF loop anchors, whole genome long-read nanopore sequencing was performed. We identified a paternally hypomethylated DMR near the SNRPN upstream loop anchor exclusive to neurons and a paternally hypermethylated DMR near the PWAR1 CTCF anchor exclusive to undifferentiated cells, consistent with increases in neuronal transcription. Additionally, DMRs near CTCF loop anchors were observed in both cell types, indicative of allele-specific differences in chromatin loops regulating imprinted transcription. These results provide an integrated view of the 15q11-q13 epigenetic landscape during LUHMES neuronal differentiation, underscoring the complex interplay of transcription, chromatin looping, and DNA methylation. They also provide insights for future therapeutic approaches for AS and PWS.

The missing link: ARID1B non-truncating variants causing Coffin-Siris syndrome due to protein aggregation

ARID1B is the most frequently mutated gene in Coffin-Siris syndrome (CSS). To date, the vast majority of causative variants reported in ARID1B are truncating, leading to nonsense-mediated mRNA decay. In the absence of experimental data, only few ARID1B amino acid substitutions have been classified as pathogenic, mainly based on clinical data and their de novo occurrence, while most others are currently interpreted as variants of unknown significance. The present study substantiates the pathogenesis of ARID1B non-truncating/NMD-escaping variants located in the SMARCA4-interacting EHD2 and DNA-binding ARID domains. Overexpression assays in cell lines revealed that the majority of EHD2 variants lead to protein misfolding and formation of cytoplasmic aggresomes surrounded by vimentin cage-like structures and co-localizing with the microtubule organisation center. ARID domain variants exhibited not only aggresomes, but also nuclear aggregates, demonstrating robust pathological effects. Protein levels were not compromised, as shown by quantitative western blot analysis. In silico structural analysis predicted the exposure of amylogenic segments in both domains due to the nearby variants, likely causing this aggregation. Genome-wide transcriptome and methylation analysis in affected individuals revealed expression and methylome patterns consistent with those of the pathogenic haploinsufficiency ARID1B alterations in CSS cases. These results further support pathogenicity and indicate two approaches for disambiguation of such variants in everyday practice. The few affected individuals harbouring EHD2 non-truncating variants described to date exhibit mild CSS clinical traits. In summary, this study paves the way for the re-evaluation of previously unclear ARID1B non-truncating variants and opens a new era in CSS genetic diagnosis.

Single-cell RNA sequencing reveals that NRF2 regulates vascular smooth muscle cell phenotypic switching in abdominal aortic aneurysm.

Abdominal aortic aneurysm (AAA) is a life-threatening disease characterized by extensive membrane destruction in the vascular wall that is closely associated with vascular smooth muscle cell (VSMC) phenotypic switching. A thorough understanding of the changes in regulatory factors during VSMC phenotypic switching is essential for managing AAA therapy. In this study, we revealed the impact of NRF2 on the modulation of VSMC phenotype and the development of AAA based on single-cell RNA sequencing analysis. By utilizing a murine model of VSMC-specific knockout of nuclear factor E2-related factor 2 (NRF2), we observed that the absence of NRF2 in VSMCs exacerbated AAA formation in an angiotensin II-induced AAA model. The downregulation of NRF2 promoted VSMC phenotypic switching, leading to an enhanced inflammatory response. Through genome-wide transcriptome analysis and loss- or gain-of-function experiments, we discovered that NRF2 upregulated the expression of VSMC contractile phenotype-specific genes by facilitating microRNA-145 (miR-145) expression. Our data identified NRF2 as a novel regulator involved in maintaining the VSMC contractile phenotype while also influencing AAA formation through an miR-145-dependent regulatory mechanism.

Gene expression patterns associated with PFOA exposure in Czech Young men and women

Perfluorooctanoic acid (PFOA), a member of per- and polyfluoroalkyl substances (PFASs), has been widely used in manufacturing for decades. Currently, PFOA is strictly regulated, but due to its high stability and persistence, it is detected in both environmental as well as in human matrices. To elucidate mechanisms of PFOA toxicity in humans, we determined the genome-wide transcriptomic changes of peripheral blood mononuclear cells (PBMC) responding to PFOA exposure in a sex-stratified analysis. This work employed samples from 145 female and 143 male participants of the CELSPAC: YA study to characterize PFOA-associated transcripts in a broader context using computational analysis. PFOA-associated gene expression differed significantly between men and women, as only 2% of mapped genes were expressed in both sexes. Disease-specific enrichment analysis revealed cancer and immune-related disease terms as those most enriched in male and female populations. Patterns of enriched terms within the gene set enrichment analysis indicated three main targets of PFOA toxicity: i) lipid metabolism for women; ii) cell cycle regulation for men; and iii) immune system response for both sexes. In summary, our genome-wide transcriptomics analysis described sex-specific differences in PFOA-associated gene expression and provided evidence about biological pathways underlying PFOA toxicity in humans.

Genetics, transcriptomics, metagenomics, and metabolomics in the pathogenesis and prediction of atrial fibrillation.

The primary cellular substrates of atrial fibrillation (AF) and the mechanisms underlying AF onset remain poorly characterized and therefore, its risk assessment lacks precision. While the use of omics may enable discovery of novel AF risk factors and narrow down the cellular pathways involved in AF pathogenesis, the work is far from complete. Large-scale genome-wide association studies and transcriptomic analyses that allow an unbiased, non-candidate-gene-based delineation of molecular changes associated with AF in humans have identified at least 150 genetic loci associated with AF. However, only few of these loci have been thoroughly mechanistically dissected, indicating that much remains to be discovered for targeted diagnostics and therapeutics. Metabolomics and metagenomics, on the other hand, add to the understanding of AF downstream of the primary substrate and integrate the signalling of environmental and host factors, respectively. These two rapidly developing fields have already provided several correlates of prevalent and incident AF that require additional validation in external cohorts and experimental studies. In this review, we take a look at the recent developments in genetics, transcriptomics, metagenomics, and metabolomics and how they may aid in improving the discovery of AF risk factors and shed light into the molecular mechanisms leading to AF onset.

Three-dimensional environment sensitizes pancreatic cancer cells to the anti-proliferative effect of budesonide by reprogramming energy metabolism

BackgroundPancreatic ductal adenocarcinoma (PDAC) is the most lethal cancer with an aggressive metastatic phenotype and very poor clinical prognosis. Interestingly, a lower occurrence of PDAC has been described in individuals with severe and long-standing asthma. Here we explored the potential link between PDAC and the glucocorticoid (GC) budesonide, a first-line therapy to treat asthma.MethodsWe tested the effect of budesonide and the classical GCs on the morphology, proliferation, migration and invasiveness of patient-derived PDAC cells and pancreatic cancer cell lines, using 2D and 3D cultures in vitro. Furthermore, a xenograft model was used to investigate the effect of budesonide on PDAC tumor growth in vivo. Finally, we combined genome-wide transcriptome analysis with genetic and pharmacological approaches to explore the mechanisms underlying budesonide activities in the different environmental conditions.ResultsWe found that in 2D culture settings, high micromolar concentrations of budesonide reduced the mesenchymal invasive/migrating features of PDAC cells, without affecting proliferation or survival. This activity was specific and independent of the Glucocorticoid Receptor (GR). Conversely, in a more physiological 3D environment, low nanomolar concentrations of budesonide strongly reduced PDAC cell proliferation in a GR-dependent manner. Accordingly, we found that budesonide reduced PDAC tumor growth in vivo. Mechanistically, we demonstrated that the 3D environment drives the cells towards a general metabolic reprogramming involving protein, lipid, and energy metabolism (e.g., increased glycolysis dependency). This metabolic change sensitizes PDAC cells to the anti-proliferative effect of budesonide, which instead induces opposite changes (e.g., increased mitochondrial oxidative phosphorylation). Finally, we provide evidence that budesonide inhibits PDAC growth, at least in part, through the tumor suppressor CDKN1C/p57Kip2.ConclusionsCollectively, our study reveals that the microenvironment influences the susceptibility of PDAC cells to GCs and provides unprecedented evidence for the anti-proliferative activity of budesonide on PDAC cells in 3D conditions, in vitro and in vivo. Our findings may explain, at least in part, the reason for the lower occurrence of pancreatic cancer in asthmatic patients and suggest a potential suitability of budesonide for clinical trials as a therapeutic approach to fight pancreatic cancer.

Salicylic Acid Is Involved in the Growth Inhibition Caused by Excessive Ammonium in Oilseed Rape (Brassica napus L.).

Rapeseed (Brassica napus L.) is extremely sensitive to excessive NH4+ toxicity. There remains incomplete knowledge of the causal factors behind the growth suppression in NH4+-nourished plants, with limited studies conducted specifically on field crop plants. In this study, we found that NH4+ toxicity significantly increased salicylic acid (SA) accumulation by accelerating the conversion of SA precursors. Moreover, exogenous SA application significantly aggravated NH4+ toxicity symptoms in the rapeseed shoots. Genome-wide differential transcriptomic analysis showed that NH4+ toxicity increased the expression of genes involved in the biosynthesis, transport, signaling transduction, and conversion of SA. SA treatment significantly increased shoot NH4+ concentrations by reducing the activities of glutamine synthase and glutamate synthase in NH4+-treated rapeseed plants. The application of an SA biosynthesis inhibitor, ABT, alleviated NH4+ toxicity symptoms. Furthermore, SA induced putrescine (Put) accumulation, resulting in an elevated ratio of Put to [spermidine (Spd) + spermine (Spm)] in the NH4+-treated plants, while the opposite was true for ABT. The application of exogenous Put and its biosynthesis inhibitor DFMA induced opposite effects on NH4+ toxicity in rapeseed shoots. These results indicated that the increased endogenous SA contributed noticeably to the toxicity caused by the sole NH4+-N supply in rapeseed shoots. This study provided fresh perspectives on the mechanism underlying excessive NH4+-induced toxicity and the corresponding alleviating strategies in plants.

Whole Genome Scan Uncovers Candidate Genes Related to Milk Production Traits in Barka Cattle.

In this study, our primary aim was to explore the genomic landscape of Barka cattle, a breed recognized for high milk production in a semi-arid environment, by focusing on genes with known roles in milk production traits. We employed genome-wide analysis and three selective sweep detection methods (ZFST, θπ ratio, and ZHp) to identify candidate genes associated with milk production and composition traits. Notably, ACAA1, P4HTM, and SLC4A4 were consistently identified by all methods. Functional annotation highlighted their roles in crucial biological processes such as fatty acid metabolism, mammary gland development, and milk protein synthesis. These findings contribute to understanding the genetic basis of milk production in Barka cattle, presenting opportunities for enhancing dairy cattle production in tropical climates. Further validation through genome-wide association studies and transcriptomic analyses is essential to fully exploit these candidate genes for selective breeding and genetic improvement in tropical dairy cattle.

Genome-Wide Transcriptomic and Metabolomic Analyses Unveiling the Defence Mechanisms of Populus tremula against Sucking and Chewing Insect Herbivores.

Plants and insects coevolved as an evolutionarily successful and enduring association. The molecular arms race led to evolutionary novelties regarding unique mechanisms of defence and detoxification in plants and insects. While insects adopt mechanisms to conquer host defence, trees develop well-orchestrated and species-specific defence strategies against insect herbivory. However, current knowledge on the molecular underpinnings of fine-tuned tree defence responses against different herbivore insects is still restricted. In the current study, using a multi-omics approach, we unveiled the defence response of Populus tremula against aphids (Chaitophorus populialbae) and spongy moths (Lymantria dispar) herbivory. Comparative differential gene expression (DGE) analyses revealed that around 272 and 1203 transcripts were differentially regulated in P. tremula after moth and aphid herbivory compared to uninfested controls. Interestingly, 5716 transcripts were differentially regulated in P. tremula between aphids and moth infestation. Further investigation showed that defence-related stress hormones and their lipid precursors, transcription factors, and signalling molecules were over-expressed, whereas the growth-related counterparts were suppressed in P. tremula after aphid and moth herbivory. Metabolomics analysis documented that around 37% of all significantly abundant metabolites were associated with biochemical pathways related to tree growth and defence. However, the metabolic profiles of aphid and moth-fed trees were quite distinct, indicating species-specific response optimization. After identifying the suitable reference genes in P. tremula, the omics data were further validated using RT-qPCR. Nevertheless, our findings documented species-specific fine-tuning of the defence response of P. tremula, showing conservation on resource allocation for defence overgrowth under aphid and moth herbivory. Such findings can be exploited to enhance our current understanding of molecular orchestration of tree responses against herbivory and aid in developing insect pest resistance P. tremula varieties.

Differential preoperative transcriptomic profiles are related to the quality of left ventricle reverse remodelling one year after valve replacement in patients with severe aortic valve stenosis

Abstract Funding Acknowledgements Type of funding sources: Public Institution(s). Main funding source(s): Spain´s Ministerio de Economía y Competitividad [(PI21/00084, CIBERCV (CB16/11/00264), co-funded by Fondo Europeo de Desarrollo Regional (FEDER)] Instituto de Investigación Sanitaria Marqués de Valdecilla, IDIVAL (INNVAL21/24; PREVAL22/09) Background Aortic valve stenosis (AS) promotes a progressive remodelling of the left ventricle (LV) characterised by cardiomyocyte hypertrophy and increased deposition of extracellular matrix elements. These adaptive plastic changes may evolve towards heart failure over time. The extent of LV remodelling and its persistence after valve replacement are significant prognostic factors for the long-term outcome of AS patients. The biological processes that determine the reverse remodelling rate after surgery of AS patients are unknown. Purpose To ascertain preoperative LV transcriptomic determinants of favourable reverse remodelling one year after surgery in AS patients. To this end, a genome-wide transcriptomic analysis was performed in LV biopsies from AS patients and surgical controls. Methods and Results The study was performed in AS (n=24) and control patients (n=12) (age and sex matched) followed up for one year after surgery. RNA-seq was performed on LV biopsies obtained during surgery. AS patients were stratified according to the quality of the reverse remodelling at one year in good (LV mass index ≤51 g/m2.7, n=12) or poor (LV mass index >51 g/m2.7, n=12) performers. RNA-seq libraries were prepared and sequenced using an Illumina Next Seq Unit. Resulting reads were mapped against GRCh38.99 version of the human transcriptome. The results from the principal component analysis revealed a distinctive expression gene pattern for the AS condition, accounting for 60.1 % of the variance, with no difference for the sex condition (8.7 %). Differential gene expression was considered significant when yielding a Benjamini-Hochberg adjusted p-value < 0.05. Compared to controls, 339 and 443 genes were differentially expressed in AS with good and poor reverse remodelling, respectively, with 509 overlapping. Gene Set Enrichment Analysis (GSEA) of differentially expressed genes in patients with poor reverse remodelling compared to control patients, showed a significant enrichment (NES > 1.5) in genes associated with angiogenesis, epithelial-mesenchymal transition (EMT), and heart growth. The down-regulated genes (NES < 0) were associated with transcription processes, such as histone modification and RNA processing-related genes, highlighting the significance of epigenetic mechanisms in AS patients exhibiting poor reverse remodelling. In contrast, such transcriptomic profile patterns were not observed in AS with good-reverse remodelling. Conclusion Aortic stenosis patients with good or poor reverse remodelling one year after aortic valve replacement exhibited distinct preoperative LV transcriptome profile. Gene ontology enrichment analysis revealed angiogenesis, EMT, and epigenetic processes as critical biological functions associated with poor reverse remodelling in AS patients.

Genome-wide investigation of HD-ZIP gene family and functional characterization of BnaHDZ149 and BnaHDZ22 in salt and drought response in Brassica napus L.

HD-ZIP proteins comprise a plant-specific transcription factor family, which play pivotal roles in plant development and adaptation to ever-changing environment. Although HD-ZIP family members have been identified in some plant species, so far our knowledge about HD-ZIP genes in rapeseed is still limited. In this study, 178 Brassica napus HD-ZIP (BnaHDZ) family members were identified in the rapeseed genome. The phylogenetic relationship, chromosomal locations, intron-exon structures, motif composition, and expression patterns of the BnaHDZ members were analyzed. The BnaHDZ family can be phylogenetically divided into four categories (Ⅰ, Ⅱ, Ⅲ and Ⅳ). Genome-wide transcriptome analysis revealed that most of the HD-ZIP I members respond to at least one abiotic stress. Two closely homologous stress-responsive HD-ZIP Ⅰ genes, BnaHDZ22 and BnaHDZ149, were identified to be involved in drought and salt responses, and selected for further functional characterization. Overexpressing BnaHDZ149 in rapeseed increased salt sensitivity of the transgenic plants, whereas overexpressing BnaHDZ22 increased sensitivity of the transgenic plants to polyethylene glycol (PEG)-simulated drought stress. This research provides not only a comprehensive landscape of BnaHDZ genes, but also a theoretical basis for elucidating the molecular mechanism of the abiotic stress responses of the HD-ZIP family in rapeseed.

Genome-wide transcriptome analyses reveal changes in glutathione-overproducing yeast obtained by ARTP mutagenesis for rice wine brewing

The atmospheric and room temperature plasma (ARTP) technique efficiently enhances yeast glutathione (GSH) yield. This study focused on obtaining a GSH over-producing strain, Saccharomyces cerevisiae SCZ40, through mutagenesis and stability testing. SCZ40 exhibited a remarkable 71.1% increase in GSH content during rice wine fermentation. Scanning electron microscopy (SEM) and whole-genome sequencing identified ARTP-induced changes in cell morphology and alterations in the protein-serine/threonine kinase and isocitrate dehydrogenase pathways linked to GSH synthesis. Transcriptomic analysis revealed 507 differentially expressed genes, highlighting the impact of ARTP on S. cerevisiae. Up/down-regulated genes associated with glycolysis/gluconeogenesis suggested a shift in glucose distribution favoring the pentose phosphate pathway and positively influencing GSH production. Pathway analysis further elucidated changes in glycolysis/gluconeogenesis, tricarboxylic acid cycle (TCA), and GSH biosynthesis pathways. Validation through quantitative real-time polymerase chain reaction (qRT-PCR) confirmed RNA-seq data reliability.

Genome-wide association study and transcriptome analysis reveal natural variation of key genes regulation flowering time in rapeseed

Genome-wide association study and transcriptome analysis reveal natural variation of key genes regulation flowering time in rapeseed

Comparison of the stage-dependent mitochondrial changes in response to pressure overload between the diseased right and left ventricle in the rat.

The right ventricle (RV) differs developmentally, anatomically and functionally from the left ventricle (LV). Therefore, characteristics of LV adaptation to chronic pressure overload cannot easily be extrapolated to the RV. Mitochondrial abnormalities are considered a crucial contributor in heart failure (HF), but have never been compared directly between RV and LV tissues and cardiomyocytes. To identify ventricle-specific mitochondrial molecular and functional signatures, we established rat models with two slowly developing disease stages (compensated and decompensated) in response to pulmonary artery banding (PAB) or ascending aortic banding (AOB). Genome-wide transcriptomic and proteomic analyses were used to identify differentially expressed mitochondrial genes and proteins and were accompanied by a detailed characterization of mitochondrial function and morphology. Two clearly distinguishable disease stages, which culminated in a comparable systolic impairment of the respective ventricle, were observed. Mitochondrial respiration was similarly impaired at the decompensated stage, while respiratory chain activity or mitochondrial biogenesis were more severely deteriorated in the failing LV. Bioinformatics analyses of the RNA-seq. and proteomic data sets identified specifically deregulated mitochondrial components and pathways. Although the top regulated mitochondrial genes and proteins differed between the RV and LV, the overall changes in tissue and cardiomyocyte gene expression were highly similar. In conclusion, mitochondrial dysfuntion contributes to disease progression in right and left heart failure. Ventricle-specific differences in mitochondrial gene and protein expression are mostly related to the extent of observed changes, suggesting that despite developmental, anatomical and functional differences mitochondrial adaptations to chronic pressure overload are comparable in both ventricles.

Genome-wide transcriptome analysis reveals that Bisphenol A activates immune responses in skeletal muscle

Cumulative human exposure to the environmental toxin, bisphenol A (BPA), has raised important health concerns in recent decades. However, the direct genomic regulation of BPA in skeletal muscles and its clinical significance are poorly understood. Therefore, we conducted a genome-wide transcriptome analysis after daily oral administration of BPA at the lowest observed adverse-effect level (LOAEL, 50 mg/kg) in male mice for six weeks to explore the gene-expression regulations in skeletal muscle induced by BPA. The primary Gene Ontology terms linked to BPA-dependent, differentially expressed genes at LOAEL comprised adaptive-immune response, positive regulation of T cell activation, and immune system process. The gene-set enrichment analysis disclosed increased complement-associated genes [complement components 3 (C3) and 4B, complement factor D, complement receptor 2, and immunoglobulin lambda constant 2] in the group administered with BPA, with a false-discovery rate of <0.05. Subsequent validation analysis conducted in BPA-fed animal skeletal muscle tissue and in vitro experiments confirmed that BPA induced immune activation, as evidenced by increased levels of C3 and C4α proteins in mice, C2C12 myoblasts, and mouse skeletal muscle cells. In addition, BPA markedly upregulated the transcription of tumor necrosis factor-α (Tnfα) in C2C12 myoblasts and mouse skeletal muscle cells, which was substantially inhibited by 5z-7-oxozeanol and parthenolide, providing further evidence of BPA-induced inflammation in muscle cells. Our bioinformatics and subsequent animal and in vitro validations demonstrate that BPA can activate inflammation in skeletal muscle, which could be a risk factor underlying chronic muscle weakness and wastage.

Genome-wide DNA methylation and transcriptomic analysis of liver tissues subjected to early ischemia/reperfusion injury upon human liver transplantation

Introduction and ObjectivesEpigenetic changes represent a mechanism connecting external stresses with long-term modifications of gene expression programs. In solid organ transplantation, ischemia-reperfusion injury (IRI) appears to induce epigenomic changes in the graft, although the currently available data are extremely limited. The present study aimed to characterize variations in DNA methylation and their effects on the transcriptome in liver transplantation from brain-dead donors. Patients and Methods12 liver grafts were evaluated through serial biopsies at different timings in the procurement-transplantation process: T0 (warm procurement, in donor), T1 (bench surgery), and T2 (after reperfusion, in recipient). DNA methylation (DNAm) and transcriptome profiles of biopsies were analyzed using microarrays and RNAseq. ResultsSignificant variations in DNAm were identified, particularly between T2 and T0. Functional enrichment of the best 1000 ranked differentially methylated promoters demonstrated that 387 hypermethylated and 613 hypomethylated promoters were involved in spliceosomal assembly and response to biotic stimuli, and inflammatory immune responses, respectively. At the transcriptome level, T2 vs. T0 showed an upregulation of 337 and downregulation of 61 genes, collectively involved in TNF-α, NFKB, and interleukin signaling. Cell enrichment analysis individuates macrophages, monocytes, and neutrophils as the most significant tissue-cell type in the response. ConclusionsIn the process of liver graft procurement-transplantation, IRI induces significant epigenetic changes that primarily act on the signaling pathways of inflammatory responses dependent on TNF-α, NFKB, and interleukins. Our DNAm datasets are the early IRI methylome literature and will serve as a launch point for studying the impact of epigenetic modification in IRI.

Integrated genome-wide association and transcriptomic analysis to identify receptor kinase genes to stripe rust resistance in wheat germplasm from southwestern China

Stripe rust of wheat, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most important diseases of wheat worldwide. Identification of new and elite Pst-resistance loci or genes has the potential to enhance overall resistance to this pathogen. Here, we conducted an integrated genome-wide association study (GWAS) and transcriptomic analysis to screen for loci associated with resistance to stripe rust in 335 accessions from Yunnan, including 311 landraces and 24 cultivars. Based on the environmental phenotype, we identified 113 protein kinases significantly associated with Pst resistance using mixed linear model (MLM) and generalized linear model (GLM) models. Transcriptomic analysis revealed that 52 of 113 protein kinases identified by GWAS were up and down regulated in response to Pst infection. Among these genes, a total of 15 receptor kinase genes were identified associated with Pst resistance. 11 candidate genes were newly discovered in Yunnan wheat germplasm. Our results revealed that resistance alleles to stripe rust were accumulated in Yunnan wheat germplasm, implying direct or indirect selection for improving stripe rust resistance in elite wheat breeding programs.

Crosstalk interactions between transcription factors ERRα and PPARα assist PPARα-mediated gene expression

ObjectiveThe peroxisome proliferator-activated receptor α (PPARα) is a transcription factor driving target genes involved in fatty acid β-oxidation. To what extent various PPARα interacting proteins may assist its function as a transcription factor is incompletely understood. An ORFeome-wide unbiased mammalian protein–protein interaction trap (MAPPIT) using PPARα as bait revealed a PPARα-ligand-dependent interaction with the orphan nuclear receptor estrogen-related receptor α (ERRα). The goal of this study was to characterize the nature of the interaction in depth and to explore whether it was of physiological relevance. MethodsWe used orthogonal protein–protein interaction assays and pharmacological inhibitors of ERRα in various systems to confirm a functional interaction and study the impact of crosstalk mechanisms. To characterize the interaction surfaces and contact points we applied a random mutagenesis screen and structural overlays. We pinpointed the extent of reciprocal ligand effects of both nuclear receptors via coregulator peptide recruitment assays. On PPARα targets revealed from a genome-wide transcriptome analysis, we performed an ERRα chromatin immunoprecipitation analysis on both fast and fed mouse livers. ResultsRandom mutagenesis scanning of PPARα's ligand-binding domain and coregulator profiling experiments supported the involvement of (a) bridging coregulator(s), while recapitulation of the interaction in vitro indicated the possibility of a trimeric interaction with RXRα. The PPARα·ERRα interaction depends on 3 C-terminal residues within helix 12 of ERRα and is strengthened by both PGC1α and serum deprivation. Pharmacological inhibition of ERRα decreased the interaction of ERRα to ligand-activated PPARα and revealed a transcriptome in line with enhanced mRNA expression of prototypical PPARα target genes, suggesting a role for ERRα as a transcriptional repressor. Strikingly, on other PPARα targets, including the isolated PDK4 enhancer, ERRα behaved oppositely. Chromatin immunoprecipitation analyses demonstrate a PPARα ligand-dependent ERRα recruitment onto chromatin at PPARα-binding regions, which is lost following ERRα inhibition in fed mouse livers. ConclusionsOur data support the coexistence of multiple layers of transcriptional crosstalk mechanisms between PPARα and ERRα, which may serve to finetune the activity of PPARα as a nutrient-sensing transcription factor.

Genetic characterization of a locus responsible for low pungency using EMS-induced mutants in Capsicum annuum L.

Key messageThe pepper mutants (‘221-2-1a’ and ‘1559-1-2h’) with very low pungency were genetically characterized. The Pun4 locus, responsible for the reduced pungency of the mutant fruits, was localized to a 208 Mb region on chromosome 6. DEMF06G16460, encoding 3-ketoacyl-CoA synthase, was proposed as a strong candidate gene based on the genetic analyses of bulked segregants, DEG, and expression analyses.Capsaicinoids are unique alkaloids present in pepper (Capsicum spp.), synthesized through the condensation of by-products from the phenylpropanoid and branched-chain fatty acid pathways, and accumulating in the placenta. In this study, we characterized two allelic ethyl methanesulfonate-induced mutant lines with extremely low pungency (‘221-2-1a’ and ‘1559-1-2h′). These mutants, derived from the pungent Korean landrace 'Yuwolcho,' exhibited lower capsaicinoid content than Yuwolcho but still contained a small amount of capsaicinoid with functional capsaicinoid biosynthetic genes. Genetic crosses between the mutants and Yuwolcho or pungent lines indicated that a single recessive mutation was responsible for the low-pungency phenotype of mutant 221-2-1a; we named the causal locus Pungency 4 (Pun4). To identify Pun4, we combined genome-wide polymorphism analysis and transcriptome analysis with bulked-segregant analysis. We narrowed down the location of Pun4 to a 208-Mb region on chromosome 6 containing five candidate genes, of which DEMF06G16460, encoding a 3-ketoacyl-CoA synthase associated with branched-chain fatty acid biosynthesis, is the most likely candidate for Pun4. The expression of capsaicinoid biosynthetic genes in placental tissues in Yuwolcho and the mutant was consistent with the branched-chain fatty acid pathway playing a pivotal role in the lower pungency observed in the mutant. We also obtained a list of differentially expressed genes in placental tissues between the mutant and Yuwolcho, from which we selected candidate genes using gene co-expression analysis. In summary, we characterized the capsaicinoid biosynthesis-related locus Pun4 through integrated of genetic, genomic, and transcriptome analyses. These findings will contribute to our understanding of capsaicinoid biosynthesis in pepper.

GDF5 as a rejuvenating treatment for age-related neuromuscular failure.

Sarcopenia involves a progressive loss of skeletal muscle force, quality and mass during ageing, which results in increased inability and death; however, no cure has been established thus far. Growth differentiation factor 5 (GDF5) has been described to modulate muscle mass maintenance in various contexts. For our proof of concept, we overexpressed GDF5 by AAV vector injection in Tibialis Anterior (TA) muscle of adult aged (20 months) mice and performed molecular and functional analysis of skeletal muscle. We analysed human Vastus Lateralis muscle biopsies from adult young (21-42 years) and aged (77-80 years) donors, quantifying the molecular markers modified by GDF5 overexpression (OE) in mouse muscle. We validated the major effects of GDF5 overexpression using human immortalized myotubes and Schwann Cells (SCs). We established a pre-clinical study by treating chronically (for 4 months) aged mice using recombinant GDF5 protein (rGDF5) in systemic administration and evaluated the long-term effect of this treatment on muscle mass and function. Here, we demonstrated that GDF5 OE in the old TAs promoted an increase of 16.5% of muscle weight (P = 0.0471) associated with a higher percentage of 5000-6000 µm2 large fibres (P = 0.0211), without the induction of muscle regeneration. Muscle mass gain was associated with an amelioration of 26.8% of rate of force generation (P = 0.0330) and a better neuromuscular connectivity (P = 0.0098). Moreover, GDF5 OE preserved neuromuscular junction (NMJ) morphology (38.5% of nerve terminal area increase, P < 0.0001) and stimulated the expression of re-innervation-related genes, in particular markers of SCs (fold change 3.19 for S100b gene expression, P = 0.0101). To further characterize the molecular events induced by GDF5 OE during ageing, we performed a genome-wide transcriptomic analysis of treated muscles and showed that this factor leads to a "rejuvenating" transcriptomic signature in aged mice, as 42% of the transcripts dysregulated by ageing reverted to youthful expression levels upon GDF5 OE (P < 0.05). Towards a pre-clinical approach, we performed a long-term systemic treatment using rGDF5 and showed its effectiveness in counteracting age-related muscle wasting, improving muscle function (17,8% of absolute maximal force increase, P = 0.0079), ensuring neuromuscular connectivity and preventing NMJ degeneration (7,96% of AchR area increase, P = 0.0125). In addition, in human muscle biopsies, we found the same age-related alterations than those observed in mice and improved by GDF5 and reproduced its major effects on human cells, suggesting this treatment as efficient in humans. Overall, these data provide a foundation to examine the curative potential of GDF5 drug in clinical trials for sarcopenia and, eventually, other neuromuscular diseases.