OBJECTIVE: Current data suggests that Bacille Calmette-Guerin (BCG) vaccination contributes to nonspecific enhancement of resistance to various infections. Thus, BCG vaccination induces both specific immunity against mycobacteria and non-specific \"trained immunity\" against various pathogens. To understand the fundamental mechanisms of \"trained\" immunity, studies of transcriptome changes occurring during BCG vaccination in innate immunity cells, as well as in their precursors, are necessary. Furthermore, this data possesses important significance for practical applications associated with the development of recombinant BCG strains aimed to enhance innate immunity against diverse infectious agents.
METHODS: We performed RNA sequencing of innate immune cells derived from murine bone marrow and spleen three days after subcutaneous BCG vaccination. Using fluorescence-activated cell sorting we obtained three cell populations for each mouse from both control and BCG vaccinated groups: bone marrow monocytes and neutrophils and splenic NK-cells. Then double-indexed cDNA libraries for Illumina sequencing from the collected samples were prepared, the resulting cDNA library mix was subjected to NovaSeq 6000 sequencing. This paper describes the collection of 24 RNA sequencing samples comprising 4 sets of immune cell populations obtained from subcutaneously BCG-vaccinated and control mice.

BACKGROUND: Host responses to infection are a major determinant of outcome. However, the existence of different response profiles in patients with endocarditis has not been addressed. Our objective was to apply transcriptomics to identify endotypes in patients with infective endocarditis.
METHODS: Thirty-two patients with infective endocarditis were studied. Clinical data and a blood sample were collected at diagnosis, and RNA sequenced. Gene expression was used to identify two clusters (endocarditis endotypes EE1 and EE2). RNA sequencing was repeated after surgery. Transcriptionally active cell populations were identified by deconvolution. Differences between endotypes in clinical data, survival, gene expression and molecular pathways involved were assessed. Identified endotypes were recapitulated in a cohort of COVID19 patients.
RESULTS: 18 and 14 patients were assigned to EE1 and EE2 respectively, with no differences in clinical data. Patients assigned to EE2 showed an enrichment in genes related to T-cell maturation and a decrease in the activation of the STAT pathway, with higher counts of active T-cells and lower counts of neutrophils. Fourteen patients (9 in EE1 and 5 in EE2) were submitted to surgery. Surgery in EE2 patients shifted gene expression towards a EE1-like profile. In-hospital mortality was higher in EE1 (56% vs 14%, p=0.027) with adjusted hazard ratio of 12.987 (95% confidence interval 3.356 - 50]. Translation of these endotypes to COVID19 and non-COVID septic patients yielded similar results in cell populations and outcome.
CONCLUSIONS: Gene expression reveals two endotypes in patients with acute endocarditis, with different underlying pathogenetic mechanisms, response to surgery and outcome.

Pearl millet is a nutri-cereal that is mostly grown in harsh environments, making it an ideal crop to study heat tolerance mechanisms at the molecular level. Despite having a better-inbuilt tolerance to high temperatures than other crops, heat stress negatively affects the crop, posing a threat to productivity gain. Hence, to understand the heat-responsive genes, the leaf and root samples of two contrasting pearl millet inbreds, EGTB 1034 (heat tolerant) and EGTB 1091 (heat sensitive), were subjected to heat-treated conditions and generated genome-wide transcriptomes. We discovered 13,464 differentially expressed genes (DEGs), of which 6932 were down-regulated and 6532 up-regulated in leaf and root tissues. The pairwise analysis of the tissue-based transcriptome data of the two genotypes demonstrated distinctive genotype and tissue-specific expression of genes. The root exhibited a higher number of DEGs compared to the leaf, emphasizing different adaptive strategies of pearl millet. A large number of genes encoding ROS scavenging enzymes, WRKY, NAC, enzymes involved in nutrient uptake, protein kinases, photosynthetic enzymes, and heat shock proteins (HSPs) and several transcription factors (TFs) involved in cross-talking of temperature stress responsive mechanisms were activated in the stress conditions. Ribosomal proteins emerged as pivotal hub genes, highly interactive with key genes expressed and involved in heat stress response. The synthesis of secondary metabolites and metabolic pathways of pearl millet were significantly enriched under heat stress. Comparative synteny analysis of HSPs and TFs in the foxtail millet genome demonstrated greater collinearity with pearl millet compared to proso millet, rice, sorghum, and maize. In this study, 1906 unannotated DEGs were identified, providing insight into novel participants in the molecular response to heat stress. The identified genes hold promise for expediting varietal development for heat tolerance in pearl millet and similar crops, fostering resilience and enhancing grain yield in heat-prone environments.

Heterocapsa bohaiensis is a newly identified dinoflagellate species that causes harmful blooms in coastal areas in China, Malaysian, and New Caledonian. These blooms have led to substantial economic losses for local aquaculture. Previous studies have mainly focused on understanding the toxicity of H. bohaiensis. However, the causes of H. bohaiensis blooms remain unknown. In this study, we aimed to ascertain nitrogen (N) and phosphorus (P) requirements for the growth and reproduction of H. bohaiensis. Additionally, we sought to understand the functional mechanisms by comparing the transcriptomes of H. bohaiensis under nutrient-limited conditions and control conditions. The results revealed a wide range of acceptable N:P ratios for H. bohainensis, attributed to a mechanism involving nutrient storage, which allowed H. bohainensis to sustain its growth even when either nitrate or phosphate was depleted. Higher N:P ratios (>27.5) were more conducive to the growth of H. bohainensis than f/2 medium or low ratios, which is related to the N:P ratios absorbed by H. bohainensis. The toxicity of H. bohainensis was significantly enhanced in N-limited or P-limited states. These findings underscore the significance of the physiological metabolism of H. bohainensis in adapting to environmental stresses induced by human activities and establishing the dominance of blooms.

Soybean rust (SBR), caused by an obligate biotrophic pathogen Phakopsora pachyrhizi, is a devastating disease of soybean worldwide. However, the mechanisms underlying plant invasion by P. pachyrhizi are poorly understood, which hinders the development of effective control strategies for SBR. Here we performed detailed histological characterization on the infection cycle of P. pachyrhizi in soybean and conducted a high-resolution transcriptional dissection of P. pachyrhizi during infection. This revealed P. pachyrhizi infection leads to significant changes in gene expression with 10 co-expressed gene modules, representing dramatic transcriptional shifts in metabolism and signal transduction during different stages throughout the infection cycle. Numerous genes encoding secreted protein are biphasic expressed, and are capable of inhibiting programmed cell death triggered by microbial effectors. Notably, three co-expressed P. pachyrhizi apoplastic effectors (PpAE1, PpAE2, and PpAE3) were found to suppress plant immune responses and were essential for P. pachyrhizi infection. Double-stranded RNA coupled with nanomaterials significantly inhibited SBR infection by targeting PpAE1, PpAE2, and PpAE3, and provided long-lasting protection to soybean against P. pachyrhizi. Together, this study revealed prominent changes in gene expression associated with SBR and identified P. pachyrhizi virulence effectors as promising targets of RNA interference-based soybean protection strategy against SBR.

Eusocial insects, such as ants and termites, are characterized by high levels of coordinated social organization. This is contrasted by solitary insects that display more limited forms of collective behavior. It has been hypothesized that this gradient in sociobehavioral sophistication is positively correlated with chemical profile complexity, due to a potentially increased demand for diversity in chemical communication mechanisms in insects with higher levels of social complexity. However, this claim has rarely been assessed empirically. Here, we compare different levels of chemical and transcriptomic complexity in selected species of the order Blattodea that represent different levels of social organization, from solitary to eusocial. We primarily focus on cuticular hydrocarbon (CHC) complexity, since it has repeatedly been demonstrated that CHCs are key signaling molecules conveying a wide variety of chemical information in solitary as well as eusocial insects. We assessed CHC complexity and divergence between our studied taxa of different social complexity levels as well as the differentiation of their respective repertoires of CHC biosynthesis gene transcripts. Surprisingly, we did not find any consistent pattern of chemical complexity correlating with social complexity, nor did the overall chemical divergence or transcriptomic repertoire of CHC biosynthesis genes reflect on the levels of social organization. Our results challenge the assumption that increasing social complexity is generally reflected in more complex chemical profiles and point toward the need for a more cautious and differentiated view on correlating complexity on a chemical, genetic, and social level.

For almost 200 years, the taxonomy of cutthroat trout (Oncorhynchus clarkii), a salmonid native to Western North America, has been in flux as ichthyologists and fisheries biologists have tried to describe the diversity within these fishes. Starting in the 1950s, Robert Behnke reexamined the cutthroat trout and identified 14 subspecies based on morphological traits, Pleistocene events, and modern geographic ranges. His designations became instrumental in recognizing and preserving the remaining diversity of cutthroat trout. Over time, molecular techniques (i.e. karyotypes, allozymes, mitochondrial DNA, SNPs, and microsatellite arrays) have largely reinforced Behnke\'s phylogenies, but have also revealed that some relationships are consistently weakly supported. To further resolve these relationships, we generated de novo transcriptomes for nine cutthroat subspecies, as well as a Bear River Bonneville form and two Colorado River lineages (blue and green). We present phylogenies of these subspecies generated from multiple sets of orthologous genes extracted from our transcriptomes. We confirm many of the relationships identified in previous morphological and molecular studies, as well as discuss the importance of significant differences apparent in our phylogenies from these studies within a geological perspective. Specific findings include three distinct clades: (1) Bear River Bonneville form and Yellowstone cutthroat trout; (2) Bonneville cutthroat trout (n = 2); and (3) Greenback and Rio Grande cutthroat trout. We also identify potential gene transfer between Bonneville cutthroat trout and a population of Colorado River green lineage cutthroat trout. Using these findings, it appears that additional groups warrant species-level consideration if other recent species elevations are retained.

Dairy cattle with high (HM) versus low muscle (LM) reserves as determined by longissimus dorsi muscle depth (LDD) in late gestation exhibit differential muscle mobilization related to subsequent milk production. Moreover, branched-chain volatile fatty acid (BCVFA) supplementation increased blood glucose levels. We hypothesized that differences in HM and LM reflect distinct muscle metabolism and that BCVFA supplementation altered metabolic pathways. At 42 days before expected calving (BEC), Holstein dairy cows were enrolled in a 2 × 2 factorial study of diet and muscle reserves, by assignment to control (CON)- or BCVFA-supplemented diets and LDD of HM (>4.6 cm) or LM (≤4.6 cm) groups: HM-CON (n = 13), HM-BCVFA (n = 10), LM-CON (n = 9), and LM-BCVFA (n = 9). Longisumus dorsi muscle was biopsied at 21 days BEC, total RNA was isolated, and protein-coding gene expression was measured with RNA sequencing. Between HM and LM, 713 genes were differentially expressed and 481 between BCVFA and CON (P < 0.05). Transcriptional signatures indicated differential distribution of type II fibers between groups, with MYH1 greater in LM cattle and MYH2 greater in HM cattle (P < 0.05). Signatures of LM cattle relative to HM cattle indicated greater activation of autophagy, ubiquitin-proteasome, and Ca2+-calpain pathways. HM cattle displayed greater expression of genes that encode extracellular matrix proteins and factors that regulate their proteolysis and turnover. BCVFA modified transcriptomes by increasing expression of genes that regulate fatty acid degradation and flux of carbons into the tricarboxylic acid cycle as acetyl CoA. Molecular signatures support distinct metabolic strategies between LM and HM cattle and that BCVFA supplementation increased substrates for energy generation.NEW & NOTEWORTHY Muscle biopsies of the longissimus dorsi of prepartum dairy cattle indicate that molecular signatures support distinct metabolic strategies between low- and high-muscle cattle and that branched-chain volatile fatty acid supplementation increased substrates for energy generation.

Genome size variation in eukaryotes has myriad effects on organismal biology from the genomic to whole-organism level. Large genome size may be associated with lower selection efficiency because lower effective population sizes allow fixation of deleterious mutations via genetic drift, increasing genome size and decreasing selection efficiency. Because of a hypothesized negative relationship between genome size and recombination rate per base pair, increased genome size could also increase the effect of linked selection in the genome, decreasing the efficiency with which natural selection can fix or remove mutations. We used a transcriptomic dataset of 15 and a subset of six Neotropical salamander species ranging in genome size from 12 to 87 pg to study the relationship between genome size and efficiency of selection. We estimated dN/dS of salamanders with small and large genomes and tested for relaxation of selection in the larger genomes. Contrary to our expectations, we did not find a significant relationship between genome size and selection efficiency or strong evidence for higher dN/dS values in species with larger genomes for either species set. We also found little evidence for relaxation of selection in species with larger genomes. A positive correlation between genome size and range size (a proxy of population size) in this group disagrees with predictions of stronger drift in species with larger genomes. Our results highlight the complex interactions between the many forces shaping genomic variation in organisms with genomic gigantism.

Epidemiological evidence suggests existing comorbidity between postmenopausal osteoporosis (OP) and cardiovascular disease (CVD), but identification of possible shared genes is lacking. The skeletal global transcriptomes were analyzed in trans-iliac bone biopsies (n = 84) from clinically well-characterized postmenopausal women (50 to 86 years) without clinical CVD using microchips and RNA sequencing. One thousand transcripts highly correlated with areal bone mineral density (aBMD) were further analyzed using bioinformatics, and common genes overlapping with CVD and associated biological mechanisms, pathways and functions were identified. Fifty genes (45 mRNAs, 5 miRNAs) were discovered with established roles in oxidative stress, inflammatory response, endothelial function, fibrosis, dyslipidemia and osteoblastogenesis/calcification. These pleiotropic genes with possible CVD comorbidity functions were also present in transcriptomes of microvascular endothelial cells and cardiomyocytes and were differentially expressed between healthy and osteoporotic women with fragility fractures. The results were supported by a genetic pleiotropy-informed conditional False Discovery Rate approach identifying any overlap in single nucleotide polymorphisms (SNPs) within several genes encoding aBMD- and CVD-associated transcripts. The study provides transcriptional and genomic evidence for genes of importance for both BMD regulation and CVD risk in a large collection of postmenopausal bone biopsies. Most of the transcripts identified in the CVD risk categories have no previously recognized roles in OP pathogenesis and provide novel avenues for exploring the mechanistic basis for the biological association between CVD and OP.