3t-seq: automatic gene expression analysis of single-copy genes, transposable elements, and tRNAs from RNA-seq data.

RNA sequencing is the gold-standard method to quantify transcriptomic changes between two conditions. The overwhelming majority of data analysis methods available are focused on polyadenylated RNA transcribed from single-copy genes and overlook transcripts from repeated sequences such as transposable elements (TEs). These self-autonomous genetic elements are increasingly studied, and specialized tools designed to handle multimapping sequencing reads are available. Transfer RNAs are transcribed by RNA polymerase III and are essential for protein translation. There is a need for integrated software that is able to analyze multiple types of RNA. Here, we present 3t-seq, a Snakemake pipeline for integrated differential expression analysis of transcripts from single-copy genes, TEs, and tRNA. 3t-seq produces an accessible report and easy-to-use results for downstream analysis starting from raw sequencing data and performing quality control, genome mapping, gene expression quantification, and statistical testing. It implements three methods to quantify TEs expression and one for tRNA genes. It provides an easy-to-configure method to manage software dependencies that lets the user focus on results. 3t-seq is released under MIT license and is available at https://github.com/boulardlab/3t-seq.

Plasticity mechanisms of genetically distinct Purkinje cells.

Despite its uniform appearance, the cerebellar cortex is highly heterogeneous in terms of structure, genetics and physiology. Purkinje cells (PCs), the principal and sole output neurons of the cerebellar cortex, can be categorized into multiple populations that differentially express molecular markers and display distinctive physiological features. Such features include action potential rate, but also their propensity for synaptic and intrinsic plasticity. However, the precise molecular and genetic factors that correlate with the differential physiological properties of PCs remain elusive. In this article, we provide a detailed overview of the cellular mechanisms that regulate PC activity and plasticity. We further perform a pathway analysis to highlight how molecular characteristics of specific PC populations may influence their physiology and plasticity mechanisms.

Nutri-cereal tissue-specific transcriptome atlas during development: Functional integration of gene expression to identify mineral uptake pathways in little millet (Panicum sumatrense).

Little millet (Panicum sumatrense Roth ex Roem. & Schult.) is an essential minor millet of southeast Asia and Africa's temperate and subtropical regions. The plant is stress-tolerant, has a short life cycle, and has a mineral-rich nutritional profile associated with unique health benefits. We report the developmental gene expression atlas of little millet (genotype JK-8) from ten tissues representing different stages of its life cycle, starting from seed germination and vegetative growth to panicle maturation. The developmental transcriptome atlas led to the identification of 342 827 transcripts. The BUSCO analysis and comparison with the transcriptomes of related species confirm that this study presents high-quality, in-depth coverage of the little millet transcriptome. In addition, the eFP browser generated here has a user-friendly interface, allowing interactive visualizations of tissue-specific gene expression. Using these data, we identified transcripts, the orthologs of which in Arabidopsis and rice are involved in nutrient acquisition, transport, and response pathways. The comparative analysis of the expression levels of these transcripts holds great potential for enhancing the mineral content in crops, particularly zinc and iron, to address the issue of "hidden hunger" and to attain nutritional security, making it a valuable asset for translational research.

A chromosome-level genome for the nudibranch gastropod Berghia stephanieae helps parse clade-specific gene expression in novel and conserved phenotypes.

BACKGROUND: How novel phenotypes originate from conserved genes, processes, and tissues remains a major question in biology. Research that sets out to answer this question often focuses on the conserved genes and processes involved, an approach that explicitly excludes the impact of genetic elements that may be classified as clade-specific, even though many of these genes are known to be important for many novel, or clade-restricted, phenotypes. This is especially true for understudied phyla such as mollusks, where limited genomic and functional biology resources for members of this phylum have long hindered assessments of genetic homology and function. To address this gap, we constructed a chromosome-level genome for the gastropod Berghia stephanieae (Valdés, 2005) to investigate the expression of clade-specific genes across both novel and conserved tissue types in this species. RESULTS: The final assembled and filtered Berghia genome is comparable to other high-quality mollusk genomes in terms of size (1.05 Gb) and number of predicted genes (24,960 genes) and is highly contiguous. The proportion of upregulated, clade-specific genes varied across tissues, but with no clear trend between the proportion of clade-specific genes and the novelty of the tissue. However, more complex tissue like the brain had the highest total number of upregulated, clade-specific genes, though the ratio of upregulated clade-specific genes to the total number of upregulated genes was low. CONCLUSIONS: Our results, when combined with previous research on the impact of novel genes on phenotypic evolution, highlight the fact that the complexity of the novel tissue or behavior, the type of novelty, and the developmental timing of evolutionary modifications will all influence how novel and conserved genes interact to generate diversity.

The expression spectrum of yak epididymal epithelial cells reveals the functional diversity of caput, corpus and cauda regions.

Sperm undergo a series of changes in the epididymis region before acquiring the ability to move and fertilize, and the identification of genes expressed in a region-specific manner in the epididymis provides a valuable insight into functional differences between regions. We collected epididymal tissue from three yaks and cultured epithelial cells from the caput, corpus and cauda regions of the yak epididymis using the tissue block method. RNA sequencing analysis (RNA-seq) technology was used to detect gene expression in yak epididymal caput, corpus and cauda epithelial cells. The results showed that the DEGs were highest in the caput vs. corpus comparison, and lowest in the corpus vs. cauda comparison. Six DEGs were verified by real-time fluorescence quantitative PCR (qRT-PCR), consistent with transcriptome sequencing results. The significantly enriched DNA replication pathway in the caput vs. corpus was coordinated with cell proliferation, while upregulated DEGs such as POLD1 and MCM4 were found in the DNA replication pathway. The AMPK signaling pathway was found significantly enriched in the caput vs cauda, suggesting its involvement in sperm maturation and capacitation. The TGF beta signaling pathway was screened in the corpus vs cauda and is crucial for mammalian reproductive regulation. Upregulated DEGs (TGFB3, INHBA, INHBB) are involved in the TGF beta signaling pathway. This study provides a reference for culturing yak epididymal epithelial cells in vitro, and elucidates the transcriptional profiles of epithelial cells in different segments of the epididymis, revealing the regulatory and functional differences between different segments, providing basic data for exploring the molecular mechanism of yak sperm maturation and improving the reproductive capacity of high-altitude mammals.

Curare and GenExVis: a versatile toolkit for analyzing and visualizing RNA-Seq data.

Even though high-throughput transcriptome sequencing is routinely performed in many laboratories, computational analysis of such data remains a cumbersome process often executed manually, hence error-prone and lacking reproducibility. For corresponding data processing, we introduce Curare, an easy-to-use yet versatile workflow builder for analyzing high-throughput RNA-Seq data focusing on differential gene expression experiments. Data analysis with Curare is customizable and subdivided into preprocessing, quality control, mapping, and downstream analysis stages, providing multiple options for each step while ensuring the reproducibility of the workflow. For a fast and straightforward exploration and visualization of differential gene expression results, we provide the gene expression visualizer software GenExVis. GenExVis can create various charts and tables from simple gene expression tables and DESeq2 results without the requirement to upload data or install software packages. In combination, Curare and GenExVis provide a comprehensive software environment that supports the entire data analysis process, from the initial handling of raw RNA-Seq data to the final DGE analyses and result visualizations, thereby significantly easing data processing and subsequent interpretation.

CrustyBase v.2.0: new features and enhanced utilities to support open science.

BACKGROUND: Transcriptomes present a rich, multi-dimensional subset of genomics data. They provide broad insights into genetic sequence, and more significantly gene expression, across biological samples. This technology is frequently employed for describing the genetic response to experimental conditions and has created vast libraries of datasets which shed light on gene function across different tissues, diseases, diets and developmental stages in many species. However, public accessibility of these data is impeded by a lack of suitable software interfaces and databases with which to locate and analyse them. BODY: Here we present an update on the status of CrustyBase.org, an online resource for analysing and sharing crustacean transcriptome datasets. Since its release in October 2020, the resource has provided many thousands of transcriptome sequences and expression profiles to its users and received 19 new dataset imports from researchers across the globe. In this article we discuss user analytics which point towards the utilization of this resource. The architecture of the application has proven robust with over 99.5% uptime and effective reporting of bugs through both user engagement and the error logging mechanism. We also introduce several new features that have been developed as part of a new release of CrustyBase.org. Two significant features are described in detail, which allow users to navigate through transcripts directly by submission of transcript identifiers, and then more broadly by searching for encoded protein domains by keyword. The latter is a novel and experimental feature, and grants users the ability to curate gene families from any dataset hosted on CrustyBase in a matter of minutes. We present case studies to demonstrate the utility of these features. CONCLUSION: Community engagement with this resource has been very positive, and we hope that improvements to the service will further enable the research of users of the platform. Web-based platforms such as CrustyBase have many potential applications across life science domains, including the health sector, which are yet to be realised. This leads to a wider discussion around the role of web-based resources in facilitating an open and collaborative research community.

Transcriptomic analysis of Anopheles gambiae from Benin reveals overexpression of salivary and cuticular proteins associated with cross-resistance to pyrethroids and organophosphates.

BACKGROUND: Insecticide resistance (IR) is one of the major threats to malaria vector control programs in endemic countries. However, the mechanisms underlying IR are poorly understood. Thus, investigating gene expression patterns related to IR can offer important insights into the molecular basis of IR in mosquitoes. In this study, RNA-Seq was used to characterize gene expression in Anopheles gambiae surviving exposure to pyrethroids (deltamethrin, alphacypermethrin) and an organophosphate (pirimiphos-methyl). RESULTS: Larvae of An. gambiae s.s. collected from Bassila and Djougou in Benin were reared to adulthood and phenotyped for IR using a modified CDC intensity bottle bioassay. The results showed that mosquitoes from Djougou were more resistant to pyrethroids (5X deltamethrin: 51.7% mortality; 2X alphacypermethrin: 47.4%) than Bassila (1X deltamethrin: 70.7%; 1X alphacypermethrin: 77.7%), while the latter were more resistant to pirimiphos-methyl (1.5X: 48.3% in Bassila and 1X: 21.5% in Djougou). RNA-seq was then conducted on resistant mosquitoes, non-exposed mosquitoes from the same locations and the laboratory-susceptible An. gambiae s.s. Kisumu strain. The results showed overexpression of detoxification genes, including cytochrome P450s (CYP12F2, CYP12F3, CYP4H15, CYP4H17, CYP6Z3, CYP9K1, CYP4G16, and CYP4D17), carboxylesterase genes (COEJHE5E, COE22933) and glutathione S-transferases (GSTE2 and GSTMS3) in all three resistant mosquito groups analyzed. Genes encoding cuticular proteins (CPR130, CPR10, CPR15, CPR16, CPR127, CPAP3-C, CPAP3-B, and CPR76) were also overexpressed in all the resistant groups, indicating their potential role in cross resistance in An. gambiae. Salivary gland protein genes related to 'salivary cysteine-rich peptide' and 'salivary secreted mucin 3' were also over-expressed and shared across all resistant groups. CONCLUSION: Our results suggest that in addition to metabolic enzymes, cuticular and salivary gland proteins could play an important role in cross-resistance to multiple classes of insecticides in Benin. These genes warrant further investigation to validate their functional role in An. gambiae resistance to insecticides.

Transcriptomic responses of Mediterranean sponges upon encounter with symbiont microbial consortia.

BACKGROUND: Sponges (phylum Porifera) constantly interact with microbes. They graze on microbes from the water column by filter-feeding and they harbor symbiotic partners within their bodies. In experimental setups, sponges take up symbionts at lower rates compared with seawater microbes. This suggests that sponges have the capacity to differentiate between microbes and preferentially graze in non-symbiotic microbes, although the underlying mechanisms of discrimination are still poorly understood. Genomic studies showed that, compared to other animal groups, sponges present an extended repertoire of immune receptors, in particular NLRs, SRCRs, and GPCRs, and a handful of experiments showed that sponges regulate the expression of these receptors upon encounter with microbial elicitors. We hypothesize that sponges may rely on differential expression of their diverse repertoire of poriferan immune receptors to sense different microbial consortia while filter-feeding. To test this, we characterized the transcriptomic response of two sponge species, Aplysina aerophoba and Dysidea avara, upon incubation with microbial consortia extracted from A. aerophoba in comparison with incubation with seawater microbes. The sponges were sampled after 1 h, 3 h, and 5 h for RNA-Seq differential gene expression analysis. RESULTS: D. avara incubated with A. aerophoba-symbionts regulated the expression of genes related to immunity, ubiquitination, and signaling. Within the set of differentially-expressed immune genes we identified different families of Nucleotide Oligomerization Domain (NOD)-Like Receptors (NLRs). These results represent the first experimental evidence that different types of NLRs are involved in microbial discrimination in a sponge. In contrast, the transcriptomic response of A. aerophoba to its own symbionts involved comparatively fewer genes and lacked genes encoding for immune receptors. CONCLUSION: Our work suggests that: (i) the transcriptomic response of sponges upon microbial exposure may imply "fine-tuning" of baseline gene expression as a result of their interaction with microbes, (ii) the differential response of sponges to microbial encounters varied between the species, probably due to species-specific characteristics or related to host's traits, and (iii) immune receptors belonging to different families of NLR-like genes played a role in the differential response to microbes, whether symbionts or food bacteria. The regulation of these receptors in sponges provides further evidence of the potential role of NLRs in invertebrate host-microbe interactions. The study of sponge responses to microbes exemplifies how investigating different animal groups broadens our knowledge of the evolution of immune specificity and symbiosis.

Identification of candidate genes associated with host-seeking behavior in the parasitoid wasp Diachasmimorpha longicaudata.

BACKGROUND: Diachasmimorpha longicaudata is a hymenopteran fruit fly endoparasitoid. Females of this species find their hosts for oviposition by using complex sensorial mechanisms in response to physical and chemical stimuli associated with the host and host habitat. Ecological and behavioral aspects related to host-seeking behavior for oviposition have been extensively studied in D. longicaudata, including the identification of volatile organic compounds acting as attractants to females. In this sense, molecular mechanisms of chemoreception have been explored in this species, including a preliminary characterization of odorant-binding proteins (OBPs), chemosensory proteins (CSPs) and odorant receptors (ORs), among other proteins. Functional assays on OBP and CSP have been conducted as a first approach to identify molecular mechanisms associated with the female host-seeking behavior for oviposition. The aims of the present study were to identify the D. longicaudata sensory gene repertoire expressed in the antenna of sexually mature and mated individuals of both sexes, and subsequently, characterize transcripts differentially expressed in the antennae of females to identify candidate genes associated with the female host-seeking behavior for oviposition. RESULTS: A total of 33,745 predicted protein-coding sequences were obtained from a de novo antennal transcriptome assembly. Ten sensory-related gene families were annotated as follows: 222 ORs, 44 ionotropic receptors (IRs), 25 gustatory receptors (GRs), 9 CSPs, 13 OBPs, 2 ammonium transporters (AMTs), 8 pickpocket (PPKs) receptors, 16 transient receptor potential (TRP) channels, 12 CD36/SNMPs and 3 Niemann-Pick type C2 like proteins (NPC2-like). The differential expression analysis revealed 237 and 151 transcripts up- and downregulated, respectively, between the female and male antennae. Ninety-seven differentially expressed transcripts corresponded to sensory-related genes including 88 transcripts being upregulated (87 ORs and one TRP) and nine downregulated (six ORs, two CSPs and one OBP) in females compared to males. CONCLUSIONS: The sensory gene repertoire of D. longicaudata was similar to that of other taxonomically related parasitoid wasps. We identified a high number of ORs upregulated in the female antenna. These results may indicate that this gene family has a central role in the chemoreception of sexually mature females during the search for hosts and host habitats for reproductive purposes.