Developmental RNA-Seq transcriptomics of haploid germ cells and spermatozoa uncovers novel pathways associated with teleost spermiogenesis

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 < 0.01) between both cell types, of which nearly half were upregulated in SPZEJ compared to HCGs. In addition, approximately 9000 long non-coding RNAs (lncRNAs) were found, of which 56% were accumulated or emerged de novo in SPZEJ. The upregulated transcripts are involved in transcriptional and translational regulation, chromatin and cytoskeleton organization, metabolic processes such as glycolysis and oxidative phosphorylation, and also include a number of ion and water channels, exchangers, transporters and receptors. Pathway analysis conducted on DEGs identified 37 different signaling pathways enriched in SPZEJ, including 13 receptor pathways, from which the most predominant correspond to the chemokine and cytokine, gonadotropin-releasing hormone receptor and platelet derived growth factor signaling pathways. Our data provide new insight into the mRNA and lncRNA cargos of teleost spermatozoa and uncover the possible involvement of novel endocrine mechanisms during the differentiation and maturation of spermatozoa.info:eu-repo/semantics/publishedVersio

Alteration in the Culex pipiens transcriptome reveals diverse mechanisms of the mosquito immune system implicated upon Rift Valley fever phlebovirus exposure

Rift Valley fever phlebovirus (RVFV) causes an emerging zoonotic disease and is mainly transmitted by Culex and Aedes mosquitoes. While Aedes aegypti-dengue virus (DENV) is the most studied model, less is known about the genes involved in infection-responses in other mosquito-arboviruses pairing. The main objective was to investigate the molecular responses of Cx. pipiens to RVFV exposure focusing mainly on genes implicated in innate immune responses. Mosquitoes were fed with blood spiked with RVFV. The fully-engorged females were pooled at 3 different time points: 2 hours post-exposure (hpe), 3- and 14-days post-exposure (dpe). Pools of mosquitoes fed with non-infected blood were also collected for comparisons. Total RNA from each mosquito pool was subjected to RNA-seq analysis and a de novo transcriptome was constructed. A total of 451 differentially expressed genes (DEG) were identified. Most of the transcriptomic alterations were found at an early infection stage after RVFV exposure. Forty-eight DEG related to immune infection-response were characterized. Most of them were related with the RNAi system, Toll and IMD pathways, ubiquitination pathway and apoptosis. Our findings provide for the first time a comprehensive view on Cx. pipiens-RVFV interactions at the molecular level. The early depletion of RNAi pathway genes at the onset of the RVFV infection would allow viral replication in mosquitoes. While genes from the Toll and IMD immune pathways were altered in response to RVFV none of the DEG were related to the JAK/STAT pathway. The fact that most of the DEG involved in the Ubiquitin-proteasome pathway (UPP) or apoptosis were found at an early stage of infection would suggest that apoptosis plays a regulatory role in infected Cx. pipiens midguts. This study provides a number of target genes that could be used to identify new molecular targets for vector control.info:eu-repo/semantics/publishedVersio

A method for multiplexed full-length single-molecule sequencing of the human mitochondrial genome

Data processing; DNA sequencing; Genomic analysisTratamiento de datos; Secuenciación de ADN; Análisis genómicoTractament de dades; Seqüenciació d'ADN; Anàlisi genòmicaMethods to reconstruct the mitochondrial DNA (mtDNA) sequence using short-read sequencing come with an inherent bias due to amplification and mapping. They can fail to determine the phase of variants, to capture multiple deletions and to cover the mitochondrial genome evenly. Here we describe a method to target, multiplex and sequence at high coverage full-length human mitochondrial genomes as native single-molecules, utilizing the RNA-guided DNA endonuclease Cas9. Combining Cas9 induced breaks, that define the mtDNA beginning and end of the sequencing reads, as barcodes, we achieve high demultiplexing specificity and delineation of the full-length of the mtDNA, regardless of the structural variant pattern. The long-read sequencing data is analysed with a pipeline where our custom-developed software, baldur, efficiently detects single nucleotide heteroplasmy to below 1%, physically determines phase and can accurately disentangle complex deletions. Our workflow is a tool for studying mtDNA variation and will accelerate mitochondrial research.This research has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 824110 – EASI-Genomics (I.G.G.) and the ERC Synergy project BCLL@las under grant agreement No 810287 (I.G.G.). Institutional support was from the Spanish Instituto de Salud Carlos III, Fondo de Investigaciones Sanitarias and cofunded with ERDF funds (PI19/01772). We acknowledge the institutional support of the Spanish Ministry of Science and Innovation through the Instituto de Salud Carlos III and the 2014–2020 Smart Growth Operating Program, to the EMBL partnership and institutional co-financing with the European Regional Development Fund (MINECO/FEDER, BIO2015-71792-P). We also acknowledge the support of the Centro de Excelencia Severo Ochoa, and the Generalitat de Catalunya through the Departament de Salut, Departament d’Empresa i Coneixement and the CERCA Programme to the institute

Developmental RNA-Seq transcriptomics of haploid germ cells and spermatozoa uncovers novel pathways associated with teleost spermiogenesis

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 (SPZ). 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 < 0.01) between both cell types, of which nearly half were upregulated in SPZ compared to HCGs. In addition, approximately 9000 long non-coding RNAs (lncRNAs) were found, of which 56% were accumulated or emerged de novo in SPZ. The upregulated transcripts are involved in transcriptional and translational regulation, chromatin and cytoskeleton organization, metabolic processes such as glycolysis and oxidative phosphorylation, and also include a number of ion and water channels, exchangers, transporters and receptors. Pathway analysis conducted on DEGs identified 37 different signaling pathways enriched in SPZ, including 13 receptor pathways, from which the most predominant correspond to the chemokine and cytokine, gonadotropin-releasing hormone receptor and platelet derived growth factor signaling pathways. Our data provide new insight into the mRNA and lncRNA cargos of teleost spermatozoa and uncover the possible involvement of novel endocrine mechanisms during the differentiation and maturation of spermatozoa

Epigenetic regulation of gene expression in response to a changing environment in CHO cell batch culture

Chinese Hamster Ovary (CHO) cells have been the workhorse for industrial production of recombinant therapeutic proteins since 1987. Variations in cellular environment and phenotypes that occur throughout the bioprocess can bring about significant changes in productivity and quality of recombinant proteins. This can potentially lead to rejection of the production lot. Hence, there is interest in an in-depth understanding of cell-line behavior and control to achieve more predictable and reliable process performance. Biological systems undergo dynamic changes over time, where individual genes are turned “on”, “off” or “paused” as and when required. So far, there is very little information available for CHO cell lines, that elucidates the effect of dynamic epigenetic regulation on temporal expression of genes in response to altered substrate availability and culture conditions. While DNA methylation levels around TSS induce either expression or silencing of genes, transcriptional regulation is primarily considered to be an interplay of transcription factors and chromatin modifiers. On top of these, there is a rapid increase in indications that connects phase-specific long non-coding RNAs (lncRNAs) in transcriptional and post-transcriptional gene regulation. Unfortunately, the mechanism of interaction of these lncRNAs with coding genes have not been studied extensively. In this study, the gene transcription dynamics throughout a batch culture of CHO cells was examined by analyzing expression profiles and histone modifications in regular 12-24 hour intervals. Chromatin states and differential methylation profiles were used to understand the role of epigenetic modifiers in the regulation of gene expression. A good correlation between expression level and absence of DNA-methylation in the promoter regions was observed. Genes having all essential active chromatin marks - specific for promoter activity, genic enhancer and active transcription, also showed significantly high positive correlation between the changes in expression levels and histone marks. Both transcription and chromatin modifications during different growth phases were found to be highly dynamic. Clusters of genes showing similar trends of expression depict gradual and continuous adaptation to the changing substrate concentrations. Less narrowly spaced temporal analyses would have prevented detection of critical regulators involved in transient changes during the batch culture. Here, we also report a plausible mode of interaction of lncRNAs with the coding genes mediated by RNA-DNA-DNA triplex formations. Based on the identified interactions, we could predict possible gene targets and the target sites for the expressed lncRNAs and show high level of correlation of expression levels between interacting pairs. To the best of our knowledge this is the first and most comprehensive report of genome wide transcriptional regulation by epigenetic modifiers for CHO. Please click Additional Files below to see the full abstract

European sea bass brain DLB-1 cell line is susceptible to nodavirus: A transcriptomic study

Viral diseases are responsible for high rates of mortality and subsequent economic losses in modern aquaculture. The nervous necrosis virus (NNV) produces viral encephalopathy and retinopathy (VER), which affects the fish central nervous system. It is considered one of the most serious viral diseases in marine aquaculture, the European sea bass (Dicentrarchus labrax) being amongst the most susceptible. We have evaluated the European sea bass brain derived cell line (DLB-1) susceptibility to NNV genotypes and evaluated its transcriptomic profile. DLB-1 cells supported NNV gene transcription and replication since strains belonging to the four NNV genotypes produce cytopathic effects. Afterwards, DLB-1 cells were infected with an RGNNV strain, the one which showed the highest replication, for 12 and 72 h and an RNA-seq analysis was performed to identify potential genes involved in the host-NNV interactions. Differential expression analysis showed the up-regulation of many genes related to immunity, heat-shock proteins or apoptosis but not to proteasome or autophagy processes. These data suggest that the immune response, mainly the interferon (IFN) pathway, is not powerful enough to abrogate the infection, and cells finally suffer stress and die by apoptosis liberating infective particles. GO enrichment also revealed, for the first time, the down-regulation of terms related to brain/neuron biology indicating molecular mechanisms causing the pathogenic effect of NNV. This study opens the way to understand key elements in sea bass brain and NNV interactions.Versión del edito

Cross-protection against African swine fever virus upon intranasal vaccination is associated with an adaptive-innate immune crosstalk

African swine fever virus (ASFV) is causing a worldwide pandemic affecting the porcine industry and leading to important global economic consequences. The virus causes a highly lethal hemorrhagic disease in wild boars and domestic pigs. Lack of effective vaccines hampers the control of virus spread, thus increasing the pressure on the scientific community for urgent solutions. However, knowledge on the immune components associated with protection is very limited. Here we characterized the in vitro recall response induced by immune cells from pigs intranasally vaccinated with the BA71ΔCD2 deletion mutant virus. Vaccination conferred dose-dependent cross-protection associated with both ASFV-specific antibodies and IFNγ-secreting cells. Importantly, bulk and single-cell transcriptomics of blood and lymph node cells from vaccinated pigs revealed a positive feedback from adaptive to innate immunity. Indeed, activation of Th1 and cytotoxic T cells was concomitant with a rapid IFNγ-dependent triggering of an inflammatory response characterized by TNF-producing macrophages, as well as CXCL10-expressing lymphocytes and cross-presenting dendritic cells. Altogether, this study provides a detailed phenotypic characterization of the immune cell subsets involved in cross-protection against ASFV, and highlights key functional immune mechanisms to be considered for the development of an effective ASF vaccine.info:eu-repo/semantics/publishedVersio

Cross-protection against African swine fever virus upon intranasal vaccination is associated with an adaptive-innate immune crosstalk

African swine fever virus (ASFV) is causing a worldwide pandemic affecting the porcine industry and leading to important global economic consequences. The virus causes a highly lethal hemorrhagic disease in wild boars and domestic pigs. Lack of effective vaccines hampers the control of virus spread, thus increasing the pressure on the scientific community for urgent solutions. However, knowledge on the immune components associated with protection is very limited. Here we characterized the in vitro recall response induced by immune cells from pigs intranasally vaccinated with the BA71ΔCD2 deletion mutant virus. Vaccination conferred dose-dependent cross-protection associated with both ASFV-specific antibodies and IFNγ-secreting cells. Importantly, bulk and single-cell transcriptomics of blood and lymph node cells from vaccinated pigs revealed a positive feedback from adaptive to innate immunity. Indeed, activation of Th1 and cytotoxic T cells was concomitant with a rapid IFNγ-dependent triggering of an inflammatory response characterized by TNF-producing macrophages, as well as CXCL10-expressing lymphocytes and cross-presenting dendritic cells. Altogether, this study provides a detailed phenotypic characterization of the immune cell subsets involved in cross-protection against ASFV, and highlights key functional immune mechanisms to be considered for the development of an effective ASF vaccine. African swine fever (ASF) pandemic is currently the number one threat for the porcine industry worldwide. Lack of treatments hampers its control, and the insufficient knowledge regarding the immune effector mechanisms required for protection hinders rational vaccine design. Here we present the first comprehensive study characterizing the complex cellular immune response involved in cross-protection against ASF. We show that, upon in vitro reactivation, cells from immune pigs induce a Th1-biased recall response that in turn enhances the antiviral innate response. Our results suggest that this positive feedback regulation of innate immunity plays a key role in the early control of ASF virus infection. Altogether, this work represents a step forward in the understanding of ASF immunology and provide critical immune components that should be considered to more rationally design future ASF vaccines

Differential expression of long non-coding RNAs are related to proliferation and histological diversity in follicular lymphomas

Long non‐coding RNAs (lncRNAs) comprise a family of non‐coding transcripts that are emerging as relevant gene expression regulators of different processes, including tumour development. To determine the possible contribution of lncRNA to the pathogenesis of follicular lymphoma (FL) we performed RNA‐sequencing at high depth sequencing in primary FL samples ranging from grade 1‐3A to aggressive grade 3B variants using unpurified (n = 16) and purified (n = 12) tumour cell suspensions from nodal samples. FL grade 3B had a significantly higher number of differentially expressed lncRNAs (dif‐lncRNAs) with potential target coding genes related to cell cycle regulation. Nine out of the 18 selected dif‐lncRNAs were validated by quantitative real time polymerase chain reaction in an independent series (n = 43) of FL. RP4‐694A7.2 was identified as the top deregulated lncRNA potentially involved in cell proliferation. RP4‐694A7.2 silencing in the WSU‐FSCCL FL cell line reduced cell proliferation due to a block in the G1/S phase. The relationship between RP4‐694A7.2 and proliferation was confirmed in primary samples as its expression levels positively related to the Ki‐67 proliferation index. In summary, lncRNAs are differentially expressed across the clinico‐biological spectrum of FL and a subset of them, related to cell cycle, may participate in cell proliferation regulation in these tumours

De novo transcriptome sequencing and gene co-expression reveal a genomic basis for drought sensitivity and evidence of a rapid local adaptation on Atlas cedar (Cedrus atlantica)

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