Search

Abstract Mutations in mitochondrial energy-producing genes lead to a heterogeneous group of untreatable disorders known as primary mitochondrial diseases (MD). Leigh syndrome (LS) is the most common pediatric MD and is characterized by progressive neuromuscular affectation and premature death. Here, we show that daily cannabidiol (CBD) administration significantly extends lifespan and ameliorates pathology in two LS mouse models, and improves cellular function in fibroblasts from LS patients. CBD delays motor decline and neurodegenerative signs, improves social deficits and breathing abnormalities, decreases thermally induced seizures, and improves neuropathology in affected brain regions. Mechanistically, we identify peroxisome proliferator-activated receptor gamma (PPARγ) as a key nuclear receptor mediating CBD’s beneficial effects, while also providing proof of dysregulated PPARγ expression and activity as a common feature in both mouse neurons and fibroblasts from LS patients. Taken together, our results provide the first evidence for CBD as a potential treatment for LS.

Background & Aims: Growing evidence suggests an increased prevalence of metabolic dysfunction-associated steatotic liver disease (MASLD) in the context of immune-mediated inflammatory diseases (IMIDs). We aimed to clinically and mechanistically characterize steatotic liver disease (SLD) in a prospective cohort of patients with IMID compared to controls. Methods: Cross-sectional, case-control study including a subset of patients with IMID. Controls from the general population were age-, sex-, type 2 diabetes-, and BMI-matched at a 1:2 ratio. SLD was established using controlled attenuation parameter. Liver biopsies were obtained when significant liver fibrosis was suspected. Total RNA was extracted from freshly frozen cases and analyzed by RNA-seq. Differential gene expression was performed with ‘limma-voom’. Gene set-enrichment analysis was performed using the fgsea R package with a preranked “limma t-statistic” gene list. Results: A total of 1,456 patients with IMID and 2,945 controls were included. Advanced SLD (liver stiffness measurement ≥9.7 kPa) (13.46% vs. 3.79%; p

IntroductionUnderstanding the adaptive capacity to current climate change of drought-sensitive tree species is mandatory, given their limited prospect of migration and adaptation as long-lived, sessile organisms. Knowledge about the molecular and eco-physiological mechanisms that control drought resilience is thus key, since water shortage appears as one of the main abiotic factors threatening forests ecosystems. However, our current background is scarce, especially in conifers, due to their huge and complex genomes.MethodsHere we investigated the eco-physiological and transcriptomic basis of drought response of the climate change-threatened conifer Cedrus atlantica. We studied C. atlantica seedlings from two locations with contrasting drought conditions to investigate a local adaptation. Seedlings were subjected to experimental drought conditions, and were monitored at immediate (24 hours) and extended (20 days) times. In addition, post-drought recovery was investigated, depicting two contrasting responses in both locations (drought resilient and non-resilient). Single nucleotide polymorphisms (SNPs) were also studied to characterize the genomic basis of drought resilience and investigate a rapid local adaptation of C. atlantica.ResultsDe novo transcriptome assembly was performed for the first time in this species, providing differences in gene expression between the immediate and extended treatments, as well as among the post-drought recovery phenotypes. Weighted gene co-expression network analysis showed a regulation of stomatal closing and photosynthetic activity during the immediate drought, consistent with an isohydric dynamic. During the extended drought, growth and flavonoid biosynthesis inhibition mechanisms prevailed, probably to increase root-to-shoot ratio and to limit the energy-intensive biosynthesis of secondary metabolites. Drought sensitive individuals failed in metabolism and photosynthesis regulation under drought stress, and in limiting secondary metabolite production. Moreover, genomic differences (SNPs) were found between drought resilient and sensitive seedlings, and between the two studied locations, which were mostly related to transposable elements.DiscussionThis work provides novel insights into the transcriptomic basis of drought response of C. atlantica, a set of candidate genes mechanistically involved in its drought sensitivity and evidence of a rapid local adaptation. Our results may help guide conservation programs for this threatened conifer, contribute to advance drought-resilience research and shed light on trees' adaptive potential to current climate change.

Abstract Acute infection and chronic infection are the two most common fates of pathogenic virus infections. While several factors that contribute to these fates are described, the critical control points and the mechanisms that underlie infection fate regulation are incompletely understood. Using the acute and chronic lymphocytic choriomeningitis virus (LCMV) infection model of mice, we find that the early dynamic pattern of the IFN-I response is a differentiating trait between both infection fates. Acute-infected mice generate a 2-wave IFN-I response while chronic-infected mice generate only a 1-wave response. The underlying cause is a temporal difference in CD8 T cell-mediated killing of splenic marginal zone CD169+ macrophages. It occurs later in acute infection and thus enables CD169+ marginal zone macrophages to produce the 2nd IFN-I wave. This is required for subsequent immune events including induction of inflammatory macrophages, generation of effector CD8+ T cells and virus clearance. Importantly, these benefits come at a cost for the host in the form of spleen fibrosis. Due to an earlier marginal zone destruction, these ordered immune events are deregulated in chronic infection. Our findings demonstrate the critical importance of kinetically well-coordinated sequential immune events for acute infection control and highlights that it may come at a cost for the host organism.

Abstract Accumulation of lipid-laden macrophages within the arterial neointima is a critical step in atherosclerotic plaque formation. Here, we show that reduced levels of the cellular plasticity factor ZEB1 in macrophages increase atherosclerotic plaque formation and the chance of cardiovascular events. Compared to control counterparts (Zeb1 WT/Apoe KO), male mice with Zeb1 ablation in their myeloid cells (Zeb1 ∆M/Apoe KO) have larger atherosclerotic plaques and higher lipid accumulation in their macrophages due to delayed lipid traffic and deficient cholesterol efflux. Zeb1 ∆M/Apoe KO mice display more pronounced systemic metabolic alterations than Zeb1 WT/Apoe KO mice, with higher serum levels of low-density lipoproteins and inflammatory cytokines and larger ectopic fat deposits. Higher lipid accumulation in Zeb1 ∆M macrophages is reverted by the exogenous expression of Zeb1 through macrophage-targeted nanoparticles. In vivo administration of these nanoparticles reduces atherosclerotic plaque formation in Zeb1 ∆M/Apoe KO mice. Finally, low ZEB1 expression in human endarterectomies is associated with plaque rupture and cardiovascular events. These results set ZEB1 in macrophages as a potential target in the treatment of atherosclerosis.

BackgroundThe BCN02-trial combined therapeutic vaccination with a viral latency reversing agent (romidepsin, RMD) in HIV-1-infected individuals and included a monitored antiretroviral pause (MAP) as an efficacy read-out identifying individuals with an early or late (< or > 4weeks) viral-rebound. Integrated -omics analyses were applied prior treatment interruption to identify markers of virus control during MAP.MethodsPBMC, whole-genome DNA methylation and transcriptomics were assessed in 14 BCN02 participants, including 8 Early and 4 Late viral-rebound individuals. Chromatin state, histone marks and integration analysis (histone-3 acetylation (H3Ac), viral load, proviral levels and HIV-specific T cells responses) were included. REDUC-trial samples (n = 5) were included as a control group for RMD administration alone.FindingsDNA methylation imprints after receiving the complete intervention discriminated Early versus Late viral-rebound individuals before MAP. Also, differential chromatin accessibility and histone marks at DNA methylation level were detected. Importantly, the differential DNA methylation positions (DMPs) between Early and Late rebounders before MAP were strongly associated with viral load, proviral levels as well as the HIV-specific T-cell responses. Most of these DMPs were already present prior to the intervention and accentuated after RMD infusion.InterpretationThis study identifies host DNA methylation profiles and epigenetic cascades that are predictive of subsequent virus control in a kick-and-kill HIV cure strategy.FundingEuropean Union Horizon 2020 Framework Programme for Research and Innovation under Grant Agreement N°681137-EAVI2020 and N°847943-MISTRAL, the Ministerio de Ciencia e Innovación (SAF2017_89726_R), and the National Institutes of Health-National Institute of Allergy and Infectious Diseases Program Grant P01-AI131568.

Abstract Acute inflammation can either resolve through immunosuppression or persist, leading to chronic inflammation. These transitions are driven by distinct molecular and metabolic reprogramming of immune cells. The anti-diabetic drug Metformin inhibits acute and chronic inflammation through mechanisms still not fully understood. Here, we report that the anti-inflammatory and reactive-oxygen-species-inhibiting effects of Metformin depend on the expression of the plasticity factor ZEB1 in macrophages. Using mice lacking Zeb1 in their myeloid cells and human patient samples, we show that ZEB1 plays a dual role, being essential in both initiating and resolving inflammation by inducing macrophages to transition into an immunosuppressed state. ZEB1 mediates these diverging effects in inflammation and immunosuppression by modulating mitochondrial content through activation of autophagy and inhibition of mitochondrial protein translation. During the transition from inflammation to immunosuppression, Metformin mimics the metabolic reprogramming of myeloid cells induced by ZEB1. Mechanistically, in immunosuppression, ZEB1 inhibits amino acid uptake, leading to downregulation of mTORC1 signalling and a decrease in mitochondrial translation in macrophages. These results identify ZEB1 as a driver of myeloid cell metabolic plasticity, suggesting that targeting its expression and function could serve as a strategy to modulate dysregulated inflammation and immunosuppression.

Duchenne and Becker muscular dystrophies, caused by pathogenic variants in DMD, are the most common inherited neuromuscular conditions in childhood. These diseases follow an X-linked recessive inheritance pattern, and mainly males are affected. The most prevalent pathogenic variants in the DMD gene are copy number variants (CNVs), and most patients achieve their genetic diagnosis through Multiplex Ligation-dependent Probe Amplification (MLPA) or exome sequencing. Here, we investigated a female patient presenting with muscular dystrophy who remained genetically undiagnosed after MLPA and exome sequencing. RNA sequencing (RNAseq) from the patient’s muscle biopsy identified an 85% reduction in DMD expression compared to 116 muscle samples included in the cohort. A de novo balanced translocation between chromosome 17 and the X chromosome (t(X;17)(p21.1;q23.2)) disrupting the DMD and BCAS3 genes was identified through trio whole genome sequencing (WGS). The combined analysis of RNAseq and WGS played a crucial role in the detection and characterisation of the disease-causing variant in this patient, who had been undiagnosed for over two decades. This case illustrates the diagnostic odyssey of female DMD patients with complex structural variants that are not detected by current panel or exome sequencing analysis.

Abstract Background New advances in high-throughput technologies have allowed for the massive analysis of genomic data, providing new opportunities for the characterization of the transcriptome architectures. Recent studies in pigs have employed RNA-Seq to explore the transcriptome of different tissues in a reduced number of animals. The main goal of this study was the identification of differentially-expressed genes in the liver of Iberian x Landrace crossbred pigs showing extreme phenotypes for intramuscular fatty acid composition using RNA-Seq. Results The liver transcriptomes of two female groups (H and L) with phenotypically extreme intramuscular fatty acid composition were sequenced using RNA-Seq. A total of 146 and 180 unannotated protein-coding genes were identified in intergenic regions for the L and H groups, respectively. In addition, a range of 5.8 to 7.3% of repetitive elements was found, with SINEs being the most abundant elements. The expression in liver of 186 (L) and 270 (H) lncRNAs was also detected. The higher reproducibility of the RNA-Seq data was validated by RT-qPCR and porcine expression microarrays, therefore showing a strong correlation between RT-qPCR and RNA-Seq data (ranking from 0.79 to 0.96), as well as between microarrays and RNA-Seq (r=0.72). A differential expression analysis between H and L animals identified 55 genes differentially-expressed between groups. Pathways analysis revealed that these genes belong to biological functions, canonical pathways and three gene networks related to lipid and fatty acid metabolism. In concordance with the phenotypic classification, the pathways analysis inferred that linolenic and arachidonic acids metabolism was altered between extreme individuals. In addition, a connection was observed among the top three networks, hence suggesting that these genes are interconnected and play an important role in lipid and fatty acid metabolism. Conclusions In the present study RNA-Seq was used as a tool to explore the liver transcriptome of pigs with extreme phenotypes for intramuscular fatty acid composition. The differential gene expression analysis showed potential gene networks which affect lipid and fatty acid metabolism. These results may help in the design of selection strategies to improve the sensorial and nutritional quality of pork meat.

Abstract Background Performing high throughput sequencing on samples pooled from different individuals is a strategy to characterize genetic variability at a small fraction of the cost required for individual sequencing. In certain circumstances some variability estimators have even lower variance than those obtained with individual sequencing. SNP calling and estimating the frequency of the minor allele from pooled samples, though, is a subtle exercise for at least three reasons. First, sequencing errors may have a much larger relevance than in individual SNP calling: while their impact in individual sequencing can be reduced by setting a restriction on a minimum number of reads per allele, this would have a strong and undesired effect in pools because it is unlikely that alleles at low frequency in the pool will be read many times. Second, the prior allele frequency for heterozygous sites in individuals is usually 0.5 (assuming one is not analyzing sequences coming from, e.g. cancer tissues), but this is not true in pools: in fact, under the standard neutral model, singletons (i.e. alleles of minimum frequency) are the most common class of variants because P(f) ∝ 1/f and they occur more often as the sample size increases. Third, an allele appearing only once in the reads from a pool does not necessarily correspond to a singleton in the set of individuals making up the pool, and vice versa, there can be more than one read – or, more likely, none – from a true singleton. Results To improve upon existing theory and software packages, we have developed a Bayesian approach for minor allele frequency (MAF) computation and SNP calling in pools (and implemented it in a program called snape): the approach takes into account sequencing errors and allows users to choose different priors. We also set up a pipeline which can simulate the coalescence process giving rise to the SNPs, the pooling procedure and the sequencing. We used it to compare the performance of snape to that of other packages. Conclusions We present a software which helps in calling SNPs in pooled samples: it has good power while retaining a low false discovery rate (FDR). The method also provides the posterior probability that a SNP is segregating and the full posterior distribution of f for every SNP. In order to test the behaviour of our software, we generated (through simulated coalescence) artificial genomes and computed the effect of a pooled sequencing protocol, followed by SNP calling. In this setting, snape has better power and False Discovery Rate (FDR) than the comparable packages samtools, PoPoolation, Varscan : for N = 50 chromosomes, snape has power ≈ 35%and FDR ≈ 2.5%. snape is available at http://code.google.com/p/snape-pooled/ (source code and precompiled binaries).

Abstract Background Although RNA-seq greatly advances our understanding of complex transcriptome landscapes, such as those found in mammals, complete RNA-seq studies in livestock and in particular in the pig are still lacking. Here, we used high-throughput RNA sequencing to gain insight into the characterization of the poly-A RNA fraction expressed in pig male gonads. An expression analysis comparing different mapping approaches and detection of allele specific expression is also discussed in this study. Results By sequencing testicle mRNA of two phenotypically extreme pigs, one Iberian and one Large White, we identified hundreds of unannotated protein-coding genes (PcGs) in intergenic regions, some of them presenting orthology with closely related species. Interestingly, we also detected 2047 putative long non-coding RNA (lncRNA), including 469 with human homologues. Two methods, DEGseq and Cufflinks, were used for analyzing expression. DEGseq identified 15% less expressed genes than Cufflinks, because DEGseq utilizes only unambiguously mapped reads. Moreover, a large fraction of the transcriptome is made up of transposable elements (14500 elements encountered), as has been reported in previous studies. Gene expression results between microarray and RNA-seq technologies were relatively well correlated (r = 0.71 across individuals). Differentially expressed genes between Large White and Iberian showed a significant overrepresentation of gamete production and lipid metabolism gene ontology categories. Finally, allelic imbalance was detected in ~ 4% of heterozygous sites. Conclusions RNA-seq is a powerful tool to gain insight into complex transcriptomes. In addition to uncovering many unnanotated genes, our study allowed us to determine that a considerable fraction is made up of long non-coding transcripts and transposable elements. Their biological roles remain to be determined in future studies. In terms of differences in expression between Large White and Iberian pigs, these were largest for genes involved in spermatogenesis and lipid metabolism, which is consistent with phenotypic extreme differences in prolificacy and fat deposition between these two breeds.

Liver is a unique organ in displaying a reparative and regenerative response after acute/chronic damage or partial hepatectomy, when all the cell types must proliferate to re-establish the liver mass. The NADPH oxidase NOX4 mediates Transforming Growth Factor-beta (TGF-β) actions, including apoptosis in hepatocytes and activation of stellate cells to myofibroblasts. Aim of this work was to analyze the impact of NOX4 in liver regeneration by using two mouse models where Nox4 was deleted: 1) general deletion of Nox4 (NOX4-/-) and 2) hepatocyte-specific deletion of Nox4 (NOX4hepKO). Liver regeneration was analyzed after 2/3 partial hepatectomy (PH). Results indicated an earlier recovery of the liver-to-body weight ratio in both NOX4-/- and NOX4hepKO mice and an increased survival, when compared to corresponding WT mice. The regenerative hepatocellular fat accumulation and the parenchyma organization recovered faster in NOX4 deleted livers. Hepatocyte proliferation, analyzed by Ki67 and phospho-Histone3 immunohistochemistry, was accelerated and increased in NOX4 deleted mice, coincident with an earlier and increased Myc expression. Primary hepatocytes isolated from NOX4 deleted mice showed higher proliferative capacity and increased expression of Myc and different cyclins in response to serum. Transcriptomic analysis through RNA-seq revealed significant changes after PH in NOX4-/- mice and support a relevant role for Myc in a node of regulation of proliferation-related genes. Interestingly, RNA-seq also revealed changes in the expression of genes related to activation of the TGF-β pathway. In fact, levels of active TGF-β1, phosphorylation of Smads and levels of its target p21 were lower at 24 h in NOX4 deleted mice. Nox4 did not appear to be essential for the termination of liver regeneration in vivo, neither for the in vitro hepatocyte response to TGF-β1 in terms of growth inhibition, which suggest its potential as therapeutic target to improve liver regeneration, without adverse effects.

Abstract Latency is a major barrier towards virus elimination in HIV-1-infected individuals. Yet, the mechanisms that contribute to the maintenance of HIV-1 latency are incompletely understood. Here we describe the Schlafen 12 protein (SLFN12) as an HIV-1 restriction factor that establishes a post-transcriptional block in HIV-1-infected cells and thereby inhibits HIV-1 replication and virus reactivation from latently infected cells. The inhibitory activity is dependent on the HIV-1 codon usage and on the SLFN12 RNase active sites. Within HIV-1-infected individuals, SLFN12 expression in PBMCs correlated with HIV-1 plasma viral loads and proviral loads suggesting a link with the general activation of the immune system. Using an RNA FISH-Flow HIV-1 reactivation assay, we demonstrate that SLFN12 expression is enriched in infected cells positive for HIV-1 transcripts but negative for HIV-1 proteins. Thus, codon-usage dependent translation inhibition of HIV-1 proteins participates in HIV-1 latency and can restrict the amount of virus release after latency reversal.

Abstract Background Chronic stress can produce a severe negative impact on health not only in the exposed individuals but also in their offspring. Indeed, chronic stress may be contributing to the current worldwide scenario of increasing infertility and decreasing gamete quality in human populations. Here, we evaluate the effect of chronic stress on behavior and male reproductive parameters in zebrafish. Our goal is to provide information on the impact that chronic stress has at molecular, histological, and physiological level in a vertebrate model species. Results We evaluated the effects of a 21-day chronic stress protocol covering around three full waves of spermatogenesis in Danio rerio adult males. The induction of chronic stress produced anxiety-like behavior in stressed males as assessed by a novel tank test. At a molecular level, the induction of chronic stress consistently resulted in the overexpression of two genes related to endoplasmic reticulum (ER) stress in the brain. Gene set enrichment analysis (GSEA) of testes suggested a dysregulation of the nonsense-mediated decay (NMD) pathway, which was also confirmed on qPCR analysis. Histological analysis of the testicle did not show significant differences in terms of the relative proportions of each germ-cell type; however, the quality of sperm from stressed males was compromised in terms of motility. RNA-seq analysis in stress-derived larval progenies revealed molecular alterations, including those predicted to affect translation initiation, DNA repair, cell cycle control, and response to stress. Conclusions Induction of chronic stress during a few cycles of spermatogenesis in the vertebrate zebrafish model affects behavior, gonadal gene expression, final gamete quality, and progeny. The NMD surveillance pathway (a key cellular mechanism that regulates the stability of both normal and mutant transcripts) is severely affected in the testes by chronic stress and therefore the control and regulation of RNAs during spermatogenesis may be affected altering the molecular status in the progeny. Graphical Abstract

Abstract Background Inclusion body myositis (IBM) is an inflammatory myopathy clinically characterized by proximal and distal muscle weakness, with inflammatory infiltrates, rimmed vacuoles and mitochondrial changes in muscle histopathology. There is scarce knowledge on IBM aetiology, and non‐established biomarkers or effective treatments are available, partly due to the lack of validated disease models. Methods We have performed transcriptomics and functional validation of IBM muscle pathological hallmarks in fibroblasts from IBM patients (n = 14) and healthy controls (n = 12), paired by age and sex. The results comprise an mRNA‐seq, together with functional inflammatory, autophagy, mitochondrial and metabolic changes between patients and controls. Results Gene expression profile of IBM vs control fibroblasts revealed 778 differentially expressed genes (P‐value adj

Highlights Fatty acid metabolism and CPT1A in AgRP neurons show marked sex-dependent differences in the control of feeding. Cpt1a gene deletion in AgRP neurons increases energy expenditure in females but not in males. CPT1A in AgRP neurons is involved in the control of thirst and fluid homeostasis. Cpt1a gene deletion in AgRP neurons induces morphological, mitochondrial, and gene expression alterations in a sex-dependent manner.

African swine fever (ASF) is a deadly disease of swine currently causing a worldwide pandemic, leading to severe economic consequences for the porcine industry. The control of disease spread is hampered by the limitation of available effective vaccines. Live attenuated vaccines (LAVs) are currently the most advanced vaccine prototypes, providing strong protection against ASF. However, the significant advances achieved using LAVs must be complemented with further studies to analyze vaccine-induced immunity. Here, we characterized the onset of cross-protective immunity triggered by the LAV candidate BA71ΔCD2. Intranasally vaccinated pigs were challenged with the virulent Georgia 2007/1 strain at days 3, 7 and 12 postvaccination. Only the animals vaccinated 12 days before the challenge had effectively controlled infection progression, showing low virus loads, minor clinical signs and a lack of the unbalanced inflammatory response characteristic of severe disease. Contrarily, the animals vaccinated 3 or 7 days before the challenge just showed a minor delay in disease progression. An analysis of the humoral response and whole blood transcriptome signatures demonstrated that the control of infection was associated with the presence of virus-specific IgG and a cytotoxic response before the challenge. These results contribute to our understanding of protective immunity induced by LAV-based vaccines, encouraging their use in emergency responses in ASF-affected areas.

Despite being in the same pathway, mutations of KRAS and BRAF in colorectal carcinomas (CRCs) determine distinct progression courses. ZEB1 induces an epithelial-to-mesenchymal transition (EMT) and is associated with worse progression in most carcinomas. Using samples from patients with CRC, mouse models of KrasG12D and BrafV600E CRC, and a Zeb1-deficient mouse, we show that ZEB1 had opposite functions in KRAS- and BRAF-mutant CRCs. In KrasG12D CRCs, ZEB1 was correlated with a worse prognosis and a higher number of larger and undifferentiated (mesenchymal or EMT-like) tumors. Surprisingly, in BrafV600E CRC, ZEB1 was associated with better prognosis; fewer, smaller, and more differentiated (reduced EMT) primary tumors; and fewer metastases. ZEB1 was positively correlated in KRAS-mutant CRC cells and negatively in BRAF-mutant CRC cells with gene signatures for EMT, cell proliferation and survival, and ERK signaling. On a mechanistic level, ZEB1 knockdown in KRAS-mutant CRC cells increased apoptosis and reduced clonogenicity and anchorage-independent growth; the reverse occurred in BRAFV600E CRC cells. ZEB1 is associated with better prognosis and reduced EMT signature in patients harboring BRAF CRCs. These data suggest that ZEB1 can function as a tumor suppressor in BRAF-mutant CRCs, highlighting the importance of considering the KRAS/BRAF mutational background of CRCs in therapeutic strategies targeting ZEB1/EMT.

Summary: Human embryonic stem cells (hESCs) can differentiate into any cell lineage. Here, we report that ZEB1 and ZEB2 promote and inhibit mesodermal-to-myogenic specification of hESCs, respectively. Knockdown and/or overexpression experiments of ZEB1, ZEB2, or PAX7 in hESCs indicate that ZEB1 is required for hESC Nodal/Activin-mediated mesodermal specification and PAX7+ human myogenic progenitor (hMuP) generation, while ZEB2 inhibits these processes. ZEB1 downregulation induces neural markers, while ZEB2 downregulation induces mesodermal/myogenic markers. Mechanistically, ZEB1 binds to and transcriptionally activates the PAX7 promoter, while ZEB2 binds to and activates the promoter of the neural OTX2 marker. Transplanting ZEB1 or ZEB2 knocked down hMuPs into the muscles of a muscular dystrophy mouse model, showing that hMuP engraftment and generation of dystrophin-positive myofibers depend on ZEB1 and are inhibited by ZEB2. The mouse model results suggest that ZEB1 expression and/or downregulating ZEB2 in hESCs may also enhance hESC regenerative capacity for human muscular dystrophy therapy.

Embryonic stem cell (ESC) derivation from single blastomeres of 8-cell mouse embryos results in lower derivation rates than that from whole blastocysts, raising a biological question about the developmental potential of sister blastomeres. We aimed to assess the ability of 8-cell blastomeres to produce epiblast cells and ESC lines after isolation, and the properties of the resulting lines. Our results revealed unequal competence among sister blastomeres to produce ESC lines. At least half of the blastomeres possess a lower potential to generate ESCs, although culture conditions and blastomeres plasticity can redirect their non-pluripotent fate towards the epiblast lineage, allowing us to generate up to seven lines from the same embryo. Lines originated from the same embryo segregated into two groups according to their transcriptional signatures. While the expression of genes related to pluripotency and development was higher in one group, no differences were found in their trilineage differentiation ability. These results may help to improve our understanding of the ESC derivation process from single blastomeres and cell fate determination in the preimplantation mouse embryos.

C-reactive protein (CRP) is an evolutionary highly conserved protein. Like humans, CRP acts as a major acute phase protein in pigs. While CRP regulatory mechanisms have been extensively studied in humans, little is known about the molecular mechanisms that control pig CRP gene expression. The main goal of the present work was to study the regulatory mechanisms and identify functional genetic variants regulating CRP gene expression and CRP blood levels in pigs. The characterization of the porcine CRP proximal promoter region revealed a high level of conservation with both cow and human promoters, sharing binding sites for transcription factors required for CRP expression. Through genome-wide association studies and fine mapping, the most associated variants with both mRNA and protein CRP levels were localized in a genomic region 39.3 kb upstream of CRP. Further study of the region revealed a highly conserved putative enhancer that contains binding sites for several transcriptional regulators such as STAT3, NF-kB or C/EBP-β. Luciferase reporter assays showed the necessity of this enhancer-promoter interaction for the acute phase induction of CRP expression in liver, where differences in the enhancer sequences significantly modified CRP activity. The associated polymorphisms disrupted the putative binding sites for HNF4α and FOXA2 transcription factors. The high correlation between HNF4α and CRP expression levels suggest the participation of HNF4α in the regulatory mechanism of porcine CRP expression through the modification of its binding site in liver. Our findings determine, for the first time, the relevance of a distal regulatory element essential for the acute phase induction of porcine CRP in liver and identify functional polymorphisms that can be included in pig breeding programs to improve immunocompetence.

Abstract Fatty acids (FAs) play an essential role as mediators of cell signaling and signal transduction, affecting metabolic homeostasis and determining meat quality in pigs. However, FAs are transformed by the action of several genes, such as those encoding desaturases and elongases of FAs in lipogenic tissues. The aim of the current work was to identify candidate genes, biological processes, and pathways involved in the modulation of intramuscular FA profile from longissimus dorsi muscle. FA profile by gas chromatography of methyl esters and gene expression by RNA-Seq were determined in 129 Iberian × Duroc backcrossed pigs. An association analysis between the muscle transcriptome and its FA profile was performed, followed by a concordance and functional analysis. Overall, a list of well-known (e.g., PLIN1, LEP, ELOVL6, SC5D, NCOA2, ACSL1, MDH1, LPL, LGALS12, TFRC, GOT1, and FBP1) and novel (e.g., TRARG1, TANK, ENSSSCG00000011196, and ENSSSCG00000038429) candidate genes was identified, either in association with specific or several FA traits. Likewise, several of these genes belong to biological processes and pathways linked to energy, lipid, and carbohydrate metabolism, which seem determinants in the modulation of FA compositions. This study can contribute to elucidate the complex relationship between gene expression and FA profile in pig muscle.

The vestibular ganglion contains primary sensory neurons that are postsynaptic to the transducing hair cells (HC) and project to the central nervous system. Understanding the response of these neurons to HC stress or loss is of great interest as their survival and functional competence will determine the functional outcome of any intervention aiming at repair or regeneration of the HCs. We have shown that subchronic exposure to the ototoxicant 3,3′-iminodipropionitrile (IDPN) in rats and mice causes a reversible detachment and synaptic uncoupling between the HCs and the ganglion neurons. Here, we used this paradigm to study the global changes in gene expression in vestibular ganglia using RNA-seq. Comparative gene ontology and pathway analyses of the data from both model species indicated a robust downregulation of terms related to synapses, including presynaptic and postsynaptic functions. Manual analyses of the most significantly downregulated transcripts identified genes with expressions related to neuronal activity, modulators of neuronal excitability, and transcription factors and receptors that promote neurite growth and differentiation. For choice selected genes, the mRNA expression results were replicated by qRT-PCR, validated spatially by RNA-scope, or were demonstrated to be associated with decreased expression of the corresponding protein. We conjectured that decreased synaptic input or trophic support on the ganglion neurons from the HC was triggering these expression changes. To support this hypothesis, we demonstrated decreased expression of BDNF mRNA in the vestibular epithelium after subchronic ototoxicity and also downregulated expression of similarly identified genes (e.g Etv5, Camk1g, Slc17a6, Nptx2, Spp1) after HC ablation with another ototoxic compound, allylnitrile. We conclude that vestibular ganglion neurons respond to decreased input from HCs by decreasing the strength of all their synaptic contacts, both as postsynaptic and presynaptic players.

Abstract In non-mammalian vertebrates, the molecular mechanisms involved in the transformation of haploid germ cells (HGCs) into spermatozoa (spermiogenesis) are largely unknown. Here, we investigated this process in the marine teleost gilthead seabream (Sparus aurata) through the examination of the changes in the transcriptome between cell-sorted HGCs and ejaculated sperm (SPZEJ). Samples were collected under strict quality controls employing immunofluorescence microscopy as well as by determining the sperm motion kinematic parameters by computer-assisted sperm analysis. Deep sequencing by RNA-seq identified a total of 7286 differentially expressed genes (DEGs) (p-value

Abstract Low‐flow vascular malformations are congenital overgrowths composed of abnormal blood vessels potentially causing pain, bleeding and obstruction of different organs. These diseases are caused by oncogenic mutations in the endothelium, which result in overactivation of the PI3K/AKT pathway. Lack of robust in vivo preclinical data has prevented the development and translation into clinical trials of specific molecular therapies for these diseases. Here, we demonstrate that the Pik3caH1047R activating mutation in endothelial cells triggers a transcriptome rewiring that leads to enhanced cell proliferation. We describe a new reproducible preclinical in vivo model of PI3K‐driven vascular malformations using the postnatal mouse retina. We show that active angiogenesis is required for the pathogenesis of vascular malformations caused by activating Pik3ca mutations. Using this model, we demonstrate that the AKT inhibitor miransertib both prevents and induces the regression of PI3K‐driven vascular malformations. We confirmed the efficacy of miransertib in isolated human endothelial cells with genotypes spanning most of human low‐flow vascular malformations.

IntroductionTumor infiltrating lymphocytes (TILs) are known to be a prognostic and predictive biomarker in breast cancer, particularly in triple negative breast cancer (TNBC) patients. International guidelines have been proposed to evaluate them in the clinical setting as a continuous variable, without a clear defined cut-off. However, there are scenarios where the immune infiltration is heterogeneous that some areas of the patient’s tumour have high numbers of TILs while other areas completely lack them. This spontaneous presentation of a heterogeneous immune infiltration could be a great opportunity to study why some tumours present TILs at diagnosis but others do not, while eliminating inter patient’s differences.MethodsIn this study, we have identified five TNBC patients that showed great TIL heterogeneity, with areas of low (≤5%) and high (≥50%) numbers of TILs in their surgical specimens. To evaluate immune infiltration heterogeneity, we performed and analyzed bulk RNA-sequencing in three independent triplicates from the high and low TIL areas of each patient.ResultsGene expression was homogeneous within the triplicates in each area but was remarkable different between TILs regions. These differences were not only due to the presence of TILs as there were other non-inflammatory genes and pathways differentially expressed between the two areas.DiscussionThis highlights the importance of intratumour heterogeneity driving the immune infiltration, and not patient’s characteristics like the HLA phenotype, germline DNA or immune repertoire.

IntroductionUnderstanding the adaptive capacity to current climate change of drought-sensitive tree species is mandatory, given their limited prospect of migration and adaptation as long-lived, sessile organisms. Knowledge about the molecular and eco-physiological mechanisms that control drought resilience is thus key, since water shortage appears as one of the main abiotic factors threatening forests ecosystems. However, our current background is scarce, especially in conifers, due to their huge and complex genomes.MethodsHere we investigated the eco-physiological and transcriptomic basis of drought response of the climate change-threatened conifer Cedrus atlantica. We studied C. atlantica seedlings from two locations with contrasting drought conditions to investigate a local adaptation. Seedlings were subjected to experimental drought conditions, and were monitored at immediate (24 hours) and extended (20 days) times. In addition, post-drought recovery was investigated, depicting two contrasting responses in both locations (drought resilient and non-resilient). Single nucleotide polymorphisms (SNPs) were also studied to characterize the genomic basis of drought resilience and investigate a rapid local adaptation of C. atlantica.ResultsDe novo transcriptome assembly was performed for the first time in this species, providing differences in gene expression between the immediate and extended treatments, as well as among the post-drought recovery phenotypes. Weighted gene co-expression network analysis showed a regulation of stomatal closing and photosynthetic activity during the immediate drought, consistent with an isohydric dynamic. During the extended drought, growth and flavonoid biosynthesis inhibition mechanisms prevailed, probably to increase root-to-shoot ratio and to limit the energy-intensive biosynthesis of secondary metabolites. Drought sensitive individuals failed in metabolism and photosynthesis regulation under drought stress, and in limiting secondary metabolite production. Moreover, genomic differences (SNPs) were found between drought resilient and sensitive seedlings, and between the two studied locations, which were mostly related to transposable elements.DiscussionThis work provides novel insights into the transcriptomic basis of drought response of C. atlantica, a set of candidate genes mechanistically involved in its drought sensitivity and evidence of a rapid local adaptation. Our results may help guide conservation programs for this threatened conifer, contribute to advance drought-resilience research and shed light on trees’ adaptive potential to current climate change.

The emerging whitefly-transmitted crinivirus tomato chlorosis virus (ToCV) causes substantial economic losses by inducing yellow leaf disorder in tomato crops. This study explores potential resistance mechanisms by examining early-stage molecular responses to ToCV. A time-course transcriptome analysis compared naïve, mock, and ToCV-infected plants at 2, 7, and 14 days post-infection (dpi). Gene expression changes were most notable at 2 and 14 dpi, likely corresponding to whitefly feeding and viral infection. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed key genes and pathways associated with ToCV infection, including those related to plant immunity, flavonoid and steroid biosynthesis, photosynthesis, and hormone signaling. Additionally, virus-derived small interfering RNAs (vsRNAs) originating from ToCV predominantly came from RNA2 and were 22 nucleotides in length. Furthermore, two genes involved in plant immunity, Hsp90 (heat shock protein 90) and its co-chaperone Sgt1 (suppressor of the G2 allele of Skp1) were targeted through viral-induced gene silencing (VIGS), showing a potential contribution to basal resistance against viral infections since their reduction correlated with increased ToCV accumulation. This study provides insights into tomato plant responses to ToCV, with potential implications for developing effective disease control strategies.

Recent research has provided compelling evidence demonstrating that paternal exposure to different stressors can influence their offspring’s phenotypes. We hypothesized that paternal stress can negatively impact the progeny, altering different miRs and triggering different physiological alterations that could compromise offspring development. To investigate this, we exposed zebrafish male siblings to a chronic stress protocol for 21 days. We performed RNA-sequencing (RNA-seq) analyses to identify differentially expressed small noncoding RNAs in 7-day postfertilization (dpf) larvae derived from paternally stressed males crossed with control females compared with the control progeny. We found a single miRNA differentially expressed—miR-29a—which was validated in larva and was also tested in the sperm, testicles, and brain of the stressed progenitors. We observed a vertical transmission of chronic stress to the unexposed larvae, reporting novel consequences of paternally inherited chronic stress at a molecular level. The deregulation of mi-R29a in those larvae could affect relevant biological processes affecting development, morphogenesis, or neurogenesis, among others. Additionally, these disruptions were associated with reduced rates of survival and hatching in the affected offspring.

Omega-3 (n-3) and omega-6 (n-6) polyunsaturated fatty acids (PUFAs) are essential fatty acids with antagonistic inflammatory functions that play vital roles in metabolic health and immune response. Current commercial swine diets tend to over-supplement with n-6 PUFAs, which may increase the likelihood of developing inflammatory diseases and affect the overall well-being of the animals. However, it is still poorly understood how n-6/n-3 PUFA ratios affect the porcine transcriptome expression and how messenger RNAs (mRNAs) and microRNAs (miRNAs) might regulate biological processes related to PUFA metabolism. On account of this, we selected a total of 20 Iberian × Duroc crossbred pigs with extreme values for n-6/n-3 FA ratio (10 high vs 10 low), and longissimus dorsi muscle samples were used to identify differentially expressed mRNAs and miRNAs. The observed differentially expressed mRNAs were associated to biological pathways related to muscle growth and immunomodulation, while the differentially expressed microRNAs (ssc-miR-30a-3p, ssc-miR-30e-3p, ssc-miR-15b and ssc-miR-7142-3p) were correlated to adipogenesis and immunity. Relevant miRNA-to-mRNA regulatory networks were also predicted (i.e., mir15b to ARRDC3; mir-7142-3p to METTL21C), and linked to lipolysis, obesity, myogenesis, and protein degradation. The n-6/n-3 PUFA ratio differences in pig skeletal muscle revealed genes, miRNAs and enriched pathways involved in lipid metabolism, cell proliferation and inflammation.