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.

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.

A Phylogenetic Framework to Simulate Synthetic Interspecies RNA-Seq Data.

Interspecies RNA-Seq datasets are increasingly common, and have the potential to answer new questions about the evolution of gene expression. Single-species differential expression analysis is now a well-studied problem that benefits from sound statistical methods. Extensive reviews on biological or synthetic datasets have provided the community with a clear picture on the relative performances of the available methods in various settings. However, synthetic dataset simulation tools are still missing in the interspecies gene expression context. In this work, we develop and implement a new simulation framework. This tool builds on both the RNA-Seq and the phylogenetic comparative methods literatures to generate realistic count datasets, while taking into account the phylogenetic relationships between the samples. We illustrate the usefulness of this new framework through a targeted simulation study, that reproduces the features of a recently published dataset, containing gene expression data in adult eye tissue across blind and sighted freshwater crayfish species. Using our simulated datasets, we perform a fair comparison of several approaches used for differential expression analysis. This benchmark reveals some of the strengths and weaknesses of both the classical and phylogenetic approaches for interspecies differential expression analysis, and allows for a reanalysis of the crayfish dataset. The tool has been integrated in the R package compcodeR, freely available on Bioconductor.

Effects of rpl1001 Gene Deletion on Cell Division of Fission Yeast and Its Molecular Mechanism.

The rpl1001 gene encodes 60S ribosomal protein L10, which is involved in intracellular protein synthesis and cell growth. However, it is not yet known whether it is involved in the regulation of cell mitosis dynamics. This study focuses on the growth, spore production, cell morphology, the dynamics of microtubules, chromosomes, actin, myosin, and mitochondria of fission yeast (Schizosaccharomyces pombe) to investigate the impact of rpl1001 deletion on cell mitosis. RNA-Seq and bioinformatics analyses were also used to reveal key genes, such as hsp16, mfm1 and isp3, and proteasome pathways. The results showed that rpl1001 deletion resulted in slow cell growth, abnormal spore production, altered cell morphology, and abnormal microtubule number and length during interphase. The cell dynamics of the rpl1001Δ strain showed that the formation of a monopolar spindle leads to abnormal chromosome segregation with increased rate of spindle elongation in anaphase of mitosis, decreased total time of division, prolonged formation time of actin and myosin loops, and increased expression of mitochondrial proteins. Analysis of the RNA-Seq sequencing results showed that the proteasome pathway, up-regulation of isp3, and down-regulation of mfm1 and mfm2 in the rpl1001Δ strain were the main factors underpinning the increased number of spore production. Also, in the rpl1001Δ strain, down-regulation of dis1 caused the abnormal microtubule and chromosome dynamics, and down-regulation of hsp16 and pgk1 were the key genes affecting the delay of actin ring and myosin ring formation. This study reveals the effect and molecular mechanism of rpl1001 gene deletion on cell division, which provides the scientific basis for further clarifying the function of the Rpl1001 protein in cell division.

De novo transcriptome assembly and global analysis of differential gene expression of aphid tolerant wild mustard Rorippa indica (L.) Hiern infested by mustard aphid Lipaphis Erysimi (L.) Kaltenbach.

Rapeseed-mustard, the oleiferous Brassica species are important oilseed crops cultivated all over the globe. Mustard aphid Lipaphis erysimi (L.) Kaltenbach is a major threat to the cultivation of rapeseed-mustard. Wild mustard Rorippa indica (L.) Hiern shows tolerance to mustard aphids as a nonhost and hence is an important source for the bioprospecting of potential resistance genes and defense measures to manage mustard aphids sustainably. We performed mRNA sequencing of the R. indica plant uninfested and infested by the mustard aphids, harvested at 24 hours post-infestation. Following quality control, the high-quality reads were subjected to de novo assembly of the transcriptome. As there is no genomic information available for this potential wild plant, the raw reads will be useful for further bioinformatics analysis and the sequence information of the assembled transcripts will be helpful to design primers for the characterization of specific gene sequences. In this study, we also used the generated resource to comprehensively analyse the global profile of differential gene expression in R. indica in response to infestation by mustard aphids. The functional enrichment analysis of the differentially expressed genes reveals a significant immune response and suggests the possibility of chitin-induced defense signaling.

Ameliorative impacts of gamma-aminobutyric acid (GABA) on seedling growth, physiological biomarkers, and gene expression in eight wheat (Triticum aestivum L.) cultivars under salt stress.

Plants spontaneously accumulate γ-aminobutyric acid (GABA), a nonprotein amino acid, in response to various stressors. Nevertheless, there is limited knowledge regarding the precise molecular mechanisms that plants employ to cope with salt stress. The objective of this study was to investigate the impact of GABA on the salt tolerance of eight distinct varieties of bread wheat (Triticum aestivum L.) by examining plant growth rates and physiological and molecular response characteristics. The application of salt stress had a detrimental impact on plant growth markers. Nevertheless, the impact was mitigated by the administration of GABA in comparison to the control treatment. When the cultivars Gemmiza 7, Gemmiza 9, and Gemmiza 12 were exposed to GABA at two distinct salt concentrations, there was a substantial increase in both the leaf chlorophyll content and photosynthetic rate. Both the control wheat cultivars and the plants exposed to salt treatment and GABA treatment showed alterations in stress-related biomarkers and antioxidants. This finding demonstrated that GABA plays a pivotal role in mitigating the impact of salt treatments on wheat cultivars. Among the eight examined kinds of wheat, CV. Gemmiza 7 and CV. Gemmiza 11 exhibited the most significant alterations in the expression of their TaSOS1 genes. CV. Misr 2, CV. Sakha 94, and CV. Sakha 95 exhibited the highest degree of variability in the expression of the NHX1, DHN3, and GR genes, respectively. The application of GABA to wheat plants enhances their ability to cope with salt stress by reducing the presence of reactive oxygen species (ROS) and other stress indicators, regulating stomatal aperture, enhancing photosynthesis, activating antioxidant enzymes, and upregulating genes involved in salt stress tolerance.

Transcriptomic analysis identifies candidate genes for Aphanomyces root rot disease resistance in pea.

BACKGROUND: Aphanomyces euteiches is a soil-borne oomycete that causes root rot in pea and other legume species. Symptoms of Aphanomyces root rot (ARR) include root discoloration and wilting, leading to significant yield losses in pea production. Resistance to ARR is known to be polygenic but the roles of single genes in the pea immune response are still poorly understood. This study uses transcriptomics to elucidate the immune response of two pea genotypes varying in their levels of resistance to A. euteiches. RESULTS: In this study, we inoculated roots of the pea (P. sativum L.) genotypes 'Linnea' (susceptible) and 'PI180693' (resistant) with two different A. euteiches strains varying in levels of virulence. The roots were harvested at 6 h post-inoculation (hpi), 20 hpi and 48 hpi, followed by differential gene expression analysis. Our results showed a time- and genotype-dependent immune response towards A. euteiches infection, involving several WRKY and MYB-like transcription factors, along with genes associated with jasmonic acid (JA) and abscisic acid (ABA) signaling. By cross-referencing with genes segregating with partial resistance to ARR, we identified 39 candidate disease resistance genes at the later stage of infection. Among the genes solely upregulated in the resistant genotype 'PI180693', Psat7g091800.1 was polymorphic between the pea genotypes and encoded a Leucine-rich repeat receptor-like kinase reminiscent of the Arabidopsis thaliana FLAGELLIN-SENSITIVE 2 receptor. CONCLUSIONS: This study provides new insights into the gene expression dynamics controlling the immune response of resistant and susceptible pea genotypes to A. euteiches infection. We present a set of 39 candidate disease resistance genes for ARR in pea, including the putative immune receptor Psat7g091800.1, for future functional validation.

LRcell: detecting the source of differential expression at the sub-cell-type level from bulk RNA-seq data.

Given most tissues are consist of abundant and diverse (sub-)cell types, an important yet unaddressed problem in bulk RNA-seq analysis is to identify at which (sub-)cell type(s) the differential expression occurs. Single-cell RNA-sequencing (scRNA-seq) technologies can answer the question, but they are often labor-intensive and cost-prohibitive. Here, we present LRcell, a computational method aiming to identify specific (sub-)cell type(s) that drives the changes observed in a bulk RNA-seq experiment. In addition, LRcell provides pre-embedded marker genes computed from putative scRNA-seq experiments as options to execute the analyses. We conduct a simulation study to demonstrate the effectiveness and reliability of LRcell. Using three different real datasets, we show that LRcell successfully identifies known cell types involved in psychiatric disorders. Applying LRcell to bulk RNA-seq results can produce a hypothesis on which (sub-)cell type(s) contributes to the differential expression. LRcell is complementary to cell type deconvolution methods.

Contrastive learning-based computational histopathology predict differential expression of cancer driver genes.

MOTIVATION: Digital pathological analysis is run as the main examination used for cancer diagnosis. Recently, deep learning-driven feature extraction from pathology images is able to detect genetic variations and tumor environment, but few studies focus on differential gene expression in tumor cells. RESULTS: In this paper, we propose a self-supervised contrastive learning framework, HistCode, to infer differential gene expression from whole slide images (WSIs). We leveraged contrastive learning on large-scale unannotated WSIs to derive slide-level histopathological features in latent space, and then transfer it to tumor diagnosis and prediction of differentially expressed cancer driver genes. Our experiments showed that our method outperformed other state-of-the-art models in tumor diagnosis tasks, and also effectively predicted differential gene expression. Interestingly, we found the genes with higher fold change can be more precisely predicted. To intuitively illustrate the ability to extract informative features from pathological images, we spatially visualized the WSIs colored by the attention scores of image tiles. We found that the tumor and necrosis areas were highly consistent with the annotations of experienced pathologists. Moreover, the spatial heatmap generated by lymphocyte-specific gene expression patterns was also consistent with the manually labeled WSIs.

Expression patterns of Mal genes and association with differential maltose and maltotriose transport rate of two Saccharomyces pastorianus yeasts.

Beer brewing is a well-known process that still faces great challenges, such as the total consumption of sugars present in the fermentation media. Lager-style beer, a major worldwide beer type, is elaborated by Saccharomyces pastorianus (Sp) yeast, which must ferment high maltotriose content worts, but its consumption represents a notable problem, especially among Sp strains belonging to group I. Factors, such as fermentation conditions, presence of maltotriose transporters, transporter copy number variation, and genetic regulation variations contribute to this issue. We assess the factors affecting fermentation in two Sp yeast strains: SpIB1, with limited maltotriose uptake, and SpIB2, known for efficient maltotriose transport. Here, SpIB2 transported significantly more maltose (28%) and maltotriose (32%) compared with SpIB1. Furthermore, SpIB2 expressed all MAL transporters (ScMALx1, SeMALx1, ScAGT1, SeAGT1, MTT1, and MPHx) on the first day of fermentation, whereas SpIB1 only exhibited ScMalx1, ScAGT1, and MPH2/3 genes. Some SpIB2 transporters had polymorphic transmembrane domains (TMD) resembling MTT1, accompanied by higher expression of these transporters and its positive regulator genes, such as MAL63. These findings suggest that, in addition to the factors mentioned above, positive regulators of Mal transporters contribute significantly to phenotypic diversity in maltose and maltotriose consumption among the studied lager yeast strains.IMPORTANCEBeer, the third most popular beverage globally with a 90% market share in the alcoholic beverage industry, relies on Saccharomyces pastorianus (Sp) strains for lager beer production. These strains exhibit phenotypic diversity in maltotriose consumption, a crucial process for the acceptable organoleptic profile in lager beer. This diversity ranges from Sp group II strains with a notable maltotriose-consuming ability to Sp group I strains with limited capacity. Our study highlights that differential gene expression of maltose and maltotriose transporters and its upstream trans-elements, such as MAL gene-positive regulators, adds complexity to this variation. This insight can contribute to a more comprehensive analysis needed to the development of controlled and efficient biotechnological processes in the beer brewing industry.

Obesity fosters severe disease outcomes in a mouse model of coronavirus infection associated with transcriptomic abnormalities.

Obesity has been identified as an independent risk factor for severe outcomes in humans with coronavirus disease 2019 (COVID-19) and other infectious diseases. Here, we established a mouse model of COVID-19 using the murine betacoronavirus, mouse hepatitis virus 1 (MHV-1). C57BL/6 and C3H/HeJ mice exposed to MHV-1 developed mild and severe disease, respectively. Obese C57BL/6 mice developed clinical manifestations similar to those of lean controls. In contrast, all obese C3H/HeJ mice succumbed by 8 days postinfection, compared to a 50% mortality rate in lean controls. Notably, both lean and obese C3H/HeJ mice exposed to MHV-1 developed lung lesions consistent with severe human COVID-19, with marked evidence of diffuse alveolar damage (DAD). To identify early predictive biomarkers of worsened disease outcomes in obese C3H/HeJ mice, we sequenced RNA from whole blood 2 days postinfection and assessed changes in gene and pathway expression. Many pathways uniquely altered in obese C3H/HeJ mice postinfection aligned with those found in humans with severe COVID-19. Furthermore, we observed altered gene expression related to the unfolded protein response and lipid metabolism in infected obese mice compared to their lean counterparts, suggesting a role in the severity of disease outcomes. This study presents a novel model for studying COVID-19 and elucidating the mechanisms underlying severe disease outcomes in obese and other hosts.

Ecological transcriptomics reveals stress response pathways of a ground-herb species in a waterlogging gradient of Amazonian riparian forests.

Environmental stress is a fundamental facet of life and a significant driver of natural selection in the wild. Gene expression diversity may facilitate adaptation to environmental changes, without necessary genetic change, but its role in adaptive divergence remains largely understudied in Neotropical systems. In Amazonian riparian forests, species distribution is predominantly influenced by species' waterlogging tolerance. The flooding gradient delineates distinct wetland forest types, shaping habitats and species characteristics. Here we investigated the molecular basis of environmental stress response in a tropical ground-herb species (Ischnosiphon puberulus) to environmental variation in Amazonian riparian forests. We compared environmental variables and gene expression profiles from individuals collected in two forest types: Igapó and Terra firme in the Amazonian riparian forests. Predictable seasonal flooding poses a significant challenge in Igapó compared to Terra firme environments, with the former presenting higher water column height and longer flooding duration. Our findings suggest that contrasting environmental conditions related to flooding regimes are important drivers of population genetic differentiation and differential gene expression in I. puberulus. Enriched gene ontology terms highlight associations with environmental stresses, such as defence response, water transport, phosphorylation, root development, response to auxin, salicylic acid and oxidative stress. By uncovering key environmental stress response pathways conserved across populations, I. puberulus offers novel genetic insights into the molecular basis of plant reactions to environmental constraints found in flooded areas of this highly biodiverse neotropical ecosystem.

De novo assembly of transcriptomes and differential gene expression analysis using short-read data from emerging model organisms - a brief guide.

Many questions in biology benefit greatly from the use of a variety of model systems. High-throughput sequencing methods have been a triumph in the democratization of diverse model systems. They allow for the economical sequencing of an entire genome or transcriptome of interest, and with technical variations can even provide insight into genome organization and the expression and regulation of genes. The analysis and biological interpretation of such large datasets can present significant challenges that depend on the 'scientific status' of the model system. While high-quality genome and transcriptome references are readily available for well-established model systems, the establishment of such references for an emerging model system often requires extensive resources such as finances, expertise and computation capabilities. The de novo assembly of a transcriptome represents an excellent entry point for genetic and molecular studies in emerging model systems as it can efficiently assess gene content while also serving as a reference for differential gene expression studies. However, the process of de novo transcriptome assembly is non-trivial, and as a rule must be empirically optimized for every dataset. For the researcher working with an emerging model system, and with little to no experience with assembling and quantifying short-read data from the Illumina platform, these processes can be daunting. In this guide we outline the major challenges faced when establishing a reference transcriptome de novo and we provide advice on how to approach such an endeavor. We describe the major experimental and bioinformatic steps, provide some broad recommendations and cautions for the newcomer to de novo transcriptome assembly and differential gene expression analyses. Moreover, we provide an initial selection of tools that can assist in the journey from raw short-read data to assembled transcriptome and lists of differentially expressed genes.

A rapid return to normal: temporal gene expression patterns following cold exposure in the bumble bee Bombus impatiens.

Bumble bees are common in cooler climates and many species likely experience periodic exposure to very cold temperatures, but little is known about the temporal dynamics of cold response mechanisms following chill exposure, especially how persistent effects of cold exposure may facilitate tolerance of future events. To investigate molecular processes involved in the temporal response by bumble bees to acute cold exposure, we compared mRNA transcript abundance in Bombus impatiens workers exposed to 0°C for 75 min (inducing chill coma) and control bees maintained at a constant ambient temperature (28°C). We sequenced the 3' end of mRNA transcripts (TagSeq) to quantify gene expression in thoracic tissue of bees at several time points (0, 10, 30, 120 and 720 min) following cold exposure. Significant differences from control bees were only detectable within 30 min after the treatment, with most occurring at the 10 min recovery time point. Genes associated with gluconeogenesis and glycolysis were most notably upregulated, while genes related to lipid and purine metabolism were downregulated. The observed patterns of expression indicate a rapid recovery after chill coma, suggesting an acute differential transcriptional response during recovery from chill coma and return to baseline expression levels within an hour, with no long-term gene expression markers of this cold exposure. Our work highlights the functions and pathways important for acute cold recovery, provides an estimated time frame for recovery from cold exposure in bumble bees, and suggests that cold hardening may be less important for these heterothermic insects.