Whole blood transcriptome signature predicts severe forms of COVID-19: Results from the COVIDeF cohort study.

COVID-19 is associated with heterogeneous outcome. Early identification of a severe progression of the disease is essential to properly manage the patients and improve their outcome. Biomarkers reflecting an increased inflammatory response, as well as individual features including advanced age, male gender, and pre-existing comorbidities, are risk factors of severe COVID-19. Yet, these features show limited accuracy for outcome prediction. The aim was to evaluate the prognostic value of whole blood transcriptome at an early stage of the disease. Blood transcriptome of patients with mild pneumonia was profiled. Patients with subsequent severe COVID-19 were compared to those with favourable outcome, and a molecular predictor based on gene expression was built. Unsupervised classification discriminated patients who would later develop a COVID-19-related severe pneumonia. The corresponding gene expression signature reflected the immune response to the viral infection dominated by a prominent type I interferon, with IFI27 among the most over-expressed genes. A 48-genes transcriptome signature predicting the risk of severe COVID-19 was built on a training cohort, then validated on an external independent cohort, showing an accuracy of 81% for predicting severe outcome. These results identify an early transcriptome signature of severe COVID-19 pneumonia, with a possible relevance to improve COVID-19 patient management.

Blood transcriptomic signature in type-2 biomarker-low severe asthma and asthma control.

BACKGROUND: Patients with type-2 (T2) cytokine-low severe asthma often have persistent symptoms despite suppression of T2 inflammation with corticosteroids. OBJECTIVES: We sought to analyze whole blood transcriptome from 738 samples in T2-biomarker-high/-low patients with severe asthma to relate transcriptomic signatures to T2 biomarkers and asthma symptom scores. METHODS: Bulk RNA-seq data were generated for blood samples (baseline, week 24, week 48) from 301 participants recruited to a randomized clinical trial of corticosteroid optimization in severe asthma. Unsupervised clustering, differential gene expression analysis, and pathway analysis were performed. Patients were grouped by T2-biomarker status and symptoms. Associations between clinical characteristics and differentially expressed genes (DEGs) associated with biomarker and symptom levels were investigated. RESULTS: Unsupervised clustering identified 2 clusters; cluster 2 patients were blood eosinophil-low/symptom-high and more likely to be receiving oral corticosteroids (OCSs). Differential gene expression analysis of these clusters, with and without stratification for OCSs, identified 2960 and 4162 DEGs, respectively. Six hundred twenty-seven of 2960 genes remained after adjusting for OCSs by subtracting OCS signature genes. Pathway analysis identified dolichyl-diphosphooligosaccharide biosynthesis and assembly of RNA polymerase I complex as significantly enriched pathways. No stable DEGs were associated with high symptoms in T2-biomarker-low patients, but numerous associated with elevated T2 biomarkers, including 15 that were upregulated at all time points irrespective of symptom level. CONCLUSIONS: OCSs have a considerable effect on whole blood transcriptome. Differential gene expression analysis demonstrates a clear T2-biomarker transcriptomic signature, but no signature was found in association with T2-biomarker-low patients, including those with a high symptom burden.

Multi-omics analysis reveals changes in tryptophan and cholesterol metabolism before and after sexual maturation in captive macaques.

Rhesus macaques (Macaca mulatta, RMs) are widely used in sexual maturation studies due to their high genetic and physiological similarity to humans. However, judging sexual maturity in captive RMs based on blood physiological indicators, female menstruation, and male ejaculation behavior can be inaccurate. Here, we explored changes in RMs before and after sexual maturation based on multi-omics analysis and identified markers for determining sexual maturity. We found that differentially expressed microbiota, metabolites, and genes before and after sexual maturation showed many potential correlations. Specifically, genes involved in spermatogenesis (TSSK2, HSP90AA1, SOX5, SPAG16, and SPATC1) were up-regulated in male macaques, and significant changes in gene (CD36), metabolites (cholesterol, 7-ketolithocholic acid, and 12-ketolithocholic acid), and microbiota (Lactobacillus) related to cholesterol metabolism were also found, suggesting the sexually mature males have stronger sperm fertility and cholesterol metabolism compared to sexually immature males. In female macaques, most differences before and after sexual maturity were related to tryptophan metabolism, including changes in IDO1, IDO2, IFNGR2, IL1Β, IL10, L-tryptophan, kynurenic acid (KA), indole-3-acetic acid (IAA), indoleacetaldehyde, and Bifidobacteria, indicating that sexually mature females exhibit stronger neuromodulation and intestinal immunity than sexually immature females. Cholesterol metabolism-related changes (CD36, 7-ketolithocholic acid, 12-ketolithocholic acid) were also observed in female and male macaques. Exploring differences before and after sexual maturation through multi-omics, we identified potential biomarkers of sexual maturity in RMs, including Lactobacillus (for males) and Bifidobacterium (for females) valuable for RM breeding and sexual maturation research.

Comparative analysis of whole blood transcriptomics between European and local Caribbean pigs in response to feed restriction in a tropical climate.

BACKGROUND: Feed restriction occurs frequently during pig growth, either due to economic reasons or stressful environmental conditions. Local breeds are suggested to have better tolerance to periods of feed restriction. However, the mechanisms underlying the response to feed restriction in different breeds is largely unknown. The aims of the present study were (1) to compare the blood transcriptome profile in response to feed restriction and refeeding of two contrasted breeds, Large White (LW), which has been selected for high performance, and Creole (CR), which is adapted to tropical conditions, and (2) to investigate the effect of a moderate feed restriction and refeeding on whole blood transcriptome. Analysis of blood transcriptome allows to study the response to feed restriction and refeeding in a dynamic way. RNAseq was performed on blood samples of growing LW and CR pigs at two time points: after 3 weeks of feed restriction and after 3 weeks of refeeding. The data was compared with samples from control animals offered the same diet on an ad libitum basis throughout the whole experiment. RESULTS: In terms of performance (body weight and feed efficiency), CR pigs were less impacted by feed restriction than LW. The transcriptional response to feed restriction and refeeding between CR and LW was contrasted both in terms of number of DEGs and enriched pathways. CR demonstrated a stronger transcriptional response to feed restriction whereas LW had a stronger response to refeeding. Differences in the transcriptional response to feed restriction between CR and LW were related to cell stress response (Aldosterone Signalling, Protein ubiquitination, Unfolded Protein Signalling) whereas after refeeding, differences were linked to thermogenesis, metabolic pathways and cell proliferation (p38 MAPK, ERK/MAPK pathway). In both breeds, transcriptional changes related to the immune response were found after restriction and refeeding. CONCLUSIONS: Altogether, the present study indicates that blood transcriptomics can be a useful tool to study differential genetic response to feed restriction in a dynamic way. The results indicate a differential response of blood gene expression to feed restriction and refeeding between breeds, affecting biological pathways that are in accordance with performance and thermoregulatory results.

The turning point of COVID-19 severity is associated with a unique circulating neutrophil gene signature.

COVID-19 has a broad spectrum of clinical manifestations associated with the host immune response heterogeneity. Despite the advances in COVID-19 research, it is still crucial to seek a panel of molecular markers that enable accurate stratification of COVID-19 patients. Here, we performed a study that combined analysis of blood transcriptome, demographic data, clinical aspects and laboratory findings from 66 participants classified into different degrees of COVID-19 severity and healthy subjects. We identified a perturbation in blood-leukocyte transcriptional profile associated with COVID-19 aggravation, which was mainly related to processes that disfavoured lymphocyte activation and favoured neutrophil activation. This transcriptional profile stratified patients according to COVID-19 severity. Hence, it enabled identification of a turning point in transcriptional dynamics that distinguished disease outcomes and non-hospitalized from hospitalized moderate patients. Central genes of this unique neutrophil signature were S100A9, ANXA3, CEACAM6, VNN1, OLFM4, IL1R2, TCN1 and CD177. Our study indicates the molecular changes that are linked with the differing clinical aspects presented by humans when suffering from COVID-19, which involve neutrophil activation.

Milk yield variation partially attributed to blood oxygen-mediated neutrophil activation in lactating dairy goats.

Blood oxygen is an essential component for numerous biological processes of mammalian animals. Milk production of ruminants largely relies on the supply of nutrients, such as glucose, amino acids and fatty acids. To define the regulatory role of blood oxygen availability in regard to milk production, seventy-five healthy Guanzhong dairy goats with similar body weight, days in milk and parities were selected. For each animal, milk yield was recorded and milk sample was collected to determine compositions. Milk vein blood was collected to determine parameters including blood gas, physio-biochemistry and haematology. Another blood sample was prepared for transcriptome and RT-qPCR. Results showed that both pressure of oxygen (pO2) in the milk vein (positively) and numbers of neutrophils in mammary vein (negatively) were associated with milk yield of the animals. To learn the role of pO2 in blood cell functionality, twelve animals (six with higher yield (H-group) and six with lower yield (L-group)) from seventy-five goats were selected. Compared with animals in L-group, goats in H-group were higher in pO2 but lower in pCO2, lactate, lactate dehydrogenase activity and neutrophil abundance in milk vein, compared with L-group. The blood transcriptome analysis suggested that compared with L-group, animals in H-group were depressed in functionality including neutrophil activation and metabolic pathways including glycolysis, NF-κB and HIF-1. Our result revealed that lower milk production could be associated with neutrophil activation responding to low pO2 in the mammary vein. In the meantime, we highlighted the potential importance of blood oxygen as a milk yield regulator.

Blood transcriptome profiling identifies 2 candidate endotypes of atopic dermatitis.

BACKGROUND: Few studies have analyzed the blood transcriptome in atopic dermatitis (AD). OBJECTIVE: We explored blood transcriptomic features of moderate to severe AD. METHODS: Blood messenger RNA sequencing on 60 adults from the TREATgermany registry including 49 patients before and after dupilumab treatment, as well as from an independent cohort of 31 patients and 43 controls was performed. Patient clustering, differential expression, correlation and coexpression network analysis, and unsupervised learning were conducted. RESULTS: AD patients showed pronounced inflammatory expression signatures with increased myeloid and IL-5-related patterns, and clearly segregated into 2 distinct clusters, with striking differences in particular for transcripts involved in eosinophil signaling. The eosinophil-high endotype showed a more pronounced global dysregulation, a positive correlation between disease activity and signatures related to IL-5 signaling, and strong correlations with several target proteins of antibodies or small molecules under development for AD. In contrast, the eosinophil-low endotype showed little transcriptomic dysregulation and no association between disease activity and gene expression. Clinical improvement with receipt of dupilumab was accompanied by a decrease of innate immune responses and an increase of lymphocyte signatures including B-cell activation and natural killer cell composition and/or function. The proportion of super responders was higher in the eosinophil-low endotype (32% vs 11%). Continued downregulation of IL18RAP, IFNG, and granzyme A in the eosinophil-high endotype suggests a residual disturbance of natural killer cell function despite clinical improvement. CONCLUSION: AD can be stratified into eosinophilic and noneosinophilic endotypes; such stratification may be useful when assessing stratified trial designs and treatment strategies.

Whole blood transcriptome characterization of young female triathlon athletes following an endurance exercise: a pilot study.

The vast majority of studies focusing on the effects of endurance exercise on hematological parameters and leukocyte gene expression were performed in adult men, so our aim was to investigate these changes in young females. Four young (age 15.3 ± 1.3 yr) elite female athletes completed an exercise session, in which they accomplished the cycling and running disciplines of a junior triathlon race. Blood samples were taken immediately before the exercise, right after the exercise, and then 1, 2, and 7 days later. Analysis of cell counts and routine biochemical parameters were complemented by RNA sequencing (RNA-seq) to whole blood samples. The applied exercise load did not trigger remarkable changes in either cardiovascular or biochemical parameters; however, it caused a significant increase in the percentage of neutrophils and a significant reduction in the ratio of lymphocytes immediately after exercise. Furthermore, endurance exercise induced a characteristic gene expression pattern change in the blood transcriptome. Gene set enrichment analysis (GSEA) using the Reactome database revealed that the expression of genes involved in immune processes and neutrophil granulocyte activation was upregulated, whereas the expression of genes important in translation and rRNA metabolism was downregulated. Comparison of a set of immune cell gene signatures (ImSig) and our transcriptomic data identified 15 overlapping genes related to T-cell functions and involved in podosome formation and adhesion to the vessel wall. Our results suggest that RNA-seq to whole blood together with ImSig analysis are useful tools for the investigation of systemic responses to endurance exercise.

Blood transcriptome analysis revealed the immune changes and immunological adaptation of wildness training giant pandas.

The giant panda (Ailuropoda melanoleuca) is a global flagship species for biodiversity conservation. As the time for captive giant pandas to be released into the wild matures, wildness training is provided to allow adaptation to their natural environment. It is assumed that changes in the immune system would be integral in this adaptation from captive to wild, where many more pathogens would be encountered in their natural habitats. Therefore, this study aims to determine the expression changes of immune-related genes and their potential as immunoassay markers for adaptation monitoring in wildness training giant pandas, and then to understand the adaptation strategy of wildness training giant pandas to the wild environment, thereby improving the success rate of panda reintroduction. We obtained 300 differentially expressed genes (DEGs) by RNA-seq, with 239 up-regulated and 61 down-regulated DEGs in wildness training giant pandas compared to captive pandas. Functional enrichment analysis indicated that up-regulated DEGs were enriched in several immune-related terms and pathways. There were 21 immune-related DEGs, in which most of them were up-regulated in wildness training giant pandas, including several critical innate and cellular immune genes. IL1R2 was the most significantly up-regulated gene and is a signature of homeostasis within the immune system. In the protein-protein interaction (PPI) analysis, CXCL8, CXCL10, and CCL5 were identified as the hub immune genes. Our results suggested that wildness training giant pandas have stronger innate and cellular immunity than captive giant pandas, and we proposed that a gene set of CXCL8, CXCL10, CCL5, CD3D, NFKBIA, TBX21, IL12RB2, and IL1R2 may serve as potential immunoassay markers to monitor and assess the immune status of wildness training giant pandas. Our study offers the first insight into immune alterations of wildness training giant pandas, paving the way for monitoring and evaluating the immune status of giant pandas when reintroducing them into the wild.

Cellular transcriptional alterations of peripheral blood in Alzheimer's disease.

BACKGROUND: Alzheimer's disease (AD), a progressive neurodegenerative disease, is the most common cause of dementia worldwide. Accumulating data support the contributions of the peripheral immune system in AD pathogenesis. However, there is a lack of comprehensive understanding about the molecular characteristics of peripheral immune cells in AD. METHODS: To explore the alterations of cellular composition and the alterations of intrinsic expression of individual cell types in peripheral blood, we performed cellular deconvolution in a large-scale bulk blood expression cohort and identified cell-intrinsic differentially expressed genes in individual cell types with adjusting for cellular proportion. RESULTS: We detected a significant increase and decrease in the proportion of neutrophils and B lymphocytes in AD blood, respectively, which had a robust replicability across other three AD cohorts, as well as using alternative algorithms. The differentially expressed genes in AD neutrophils were enriched for some AD-associated pathways, such as ATP metabolic process and mitochondrion organization. We also found a significant enrichment of protein-protein interaction network modules of leukocyte cell-cell activation, mitochondrion organization, and cytokine-mediated signaling pathway in neutrophils for AD risk genes including CD33 and IL1B. Both changes in cellular composition and expression levels of specific genes were significantly associated with the clinical and pathological alterations. A similar pattern of perturbations on the cellular proportion and gene expression levels of neutrophils could be also observed in mild cognitive impairment (MCI). Moreover, we noticed an elevation of neutrophil abundance in the AD brains. CONCLUSIONS: We revealed the landscape of molecular perturbations at the cellular level for AD. These alterations highlight the putative roles of neutrophils in AD pathobiology.

Postseptic Cognitive Impairment and Expression of APOE in Peripheral Blood: The Cognition After SepsiS (CASS) Observational Pilot Study.

BACKGROUND: Cognitive impairment after sepsis is an important clinical problem. Determinants of postseptic cognitive impairment are not well understood. We thus undertook a systems biology approach to exploring a possible role for apolipoprotein E (APOE) in postseptic cognitive impairment. DESIGN: Prospective, observational cohort. SETTING: Intermountain Medical Center, a tertiary referral center in Utah. PATIENTS/PARTICIPANTS: Patients with sepsis admitted to study intensive care units. INTERVENTIONS: None. METHODS: We obtained peripheral blood for deep sequencing of RNA and followed up survivors at 6 months with a battery of cognitive instruments. We defined cognitive impairment based on the 6-month Hayling test of executive function. In our primary analysis, we employed weighted network analysis. Secondarily, we compared variation in gene expression between patients with normal versus impaired cognition. MEASUREMENTS AND MAIN RESULTS: We enrolled 40 patients, of whom 34 were follow-up eligible and 31 (91%) completed follow-up; 1 patient's RNA sample was degraded-the final analytic cohort was 30 patients. Mean Hayling test score was 5.8 (standard deviation 1.1), which represented 20% with impaired executive function. The network module containing APOE was dominated by low-expression genes, with no association on primary analysis (P = .8). Secondary analyses suggested several potential lines of future investigation, including oxidative stress. CONCLUSIONS: In this prospective pilot cohort, executive dysfunction affected 1 in 5 survivors of sepsis. The APOE gene was sparsely transcribed in peripheral leukocytes and not associated with cognitive impairment. Future lines of research are suggested.

Blood Transcriptomic Stratification of Short-term Risk in Contacts of Tuberculosis.

BACKGROUND: The highest risk of tuberculosis arises in the first few months after exposure. We reasoned that this risk reflects incipient disease among tuberculosis contacts. Blood transcriptional biomarkers of tuberculosis may predate clinical diagnosis, suggesting they offer improved sensitivity to detect subclinical incipient disease. Therefore, we sought to test the hypothesis that refined blood transcriptional biomarkers of active tuberculosis will improve stratification of short-term disease risk in tuberculosis contacts. METHODS: We combined analysis of previously published blood transcriptomic data with new data from a prospective human immunodeficiency virus (HIV)-negative UK cohort of 333 tuberculosis contacts. We used stability selection as an alternative computational approach to identify an optimal signature for short-term risk of active tuberculosis and evaluated its predictive value in independent cohorts. RESULTS: In a previously published HIV-negative South African case-control study of patients with asymptomatic Mycobacterium tuberculosis infection, a novel 3-gene transcriptional signature comprising BATF2, GBP5, and SCARF1 achieved a positive predictive value (PPV) of 23% for progression to active tuberculosis within 90 days. In a new UK cohort of 333 HIV-negative tuberculosis contacts with a median follow-up of 346 days, this signature achieved a PPV of 50% (95% confidence interval [CI], 15.7-84.3) and negative predictive value of 99.3% (95% CI, 97.5-99.9). By comparison, peripheral blood interferon gamma release assays in the same cohort achieved a PPV of 5.6% (95% CI, 2.1-11.8). CONCLUSIONS: This blood transcriptional signature provides unprecedented opportunities to target therapy among tuberculosis contacts with greatest risk of incident disease.

Comparative analysis of lung and blood transcriptomes in mice exposed to multi-walled carbon nanotubes.

Pulmonary exposure to multi-walled carbon nanotubes (MWCNT) causes inflammation, fibroproliferation, immunotoxicity, and systemic responses in rodents. However, the search for representative biomarkers of exposure is an ongoing endeavor. Whole blood gene expression profiling is a promising new approach for the identification of novel disease biomarkers. We asked if the whole blood transcriptome reflects pathology-specific changes in lung gene expression caused by MWCNT. To answer this question, we performed mRNA sequencing analysis of the whole blood and lung in mice administered MWCNT or vehicle solution via pharyngeal aspiration and sacrificed 56 days later. The pattern of lung mRNA expression as determined using Ingenuity Pathway Analysis (IPA) was indicative of continued inflammation, immune cell trafficking, phagocytosis, and adaptive immune responses. Simultaneously, innate immunity-related transcripts (Plunc, Bpifb1, Reg3g) and cancer-related pathways were downregulated. IPA analysis of the differentially expressed genes in the whole blood suggested increased hematopoiesis, predicted activation of cancer/tumor development pathways, and atopy. There were several common upregulated genes between whole blood and lungs, important for adaptive immune responses: Cxcr1, Cd72, Sharpin, and Slc11a1. Trim24, important for TH2 cell effector function, was downregulated in both datasets. Hla-dqa1 mRNA was upregulated in the lungs and downregulated in the blood, as was Lilrb4, which controls the reactivity of immune response. "Cancer" disease category had opposing activation status in the two datasets, while the only commonality was "Hypersensitivity". Transcriptome changes occurring in the lungs did not produce a completely replicable pattern in whole blood; however, specific systemic responses may be shared between transcriptomic profiles.

Transcriptome analysis revealed potential mechanisms of differences in physiological stress responses between caged male and female magpies.

BACKGROUND: Under caged conditions, birds are affected more severely by environmental stressors such as dietary structure, activity space, human disturbances, and pathogens, which may be reflected in the gene expression in peripheral blood or other tissues. Elucidating the molecular mechanism of these stress responses will help improve animal welfare. RESULTS: In the present study, the blood transcriptomes of six male and five female caged magpies (Pica pica) were sequenced, and a total of ~ 100 Gb in clean reads were generated using the Illumina HiSeq 2000 sequencer. A total of 420,291 unigenes were identified after assembly, of which 179,316 were annotated in five databases, 7471 were assigned to 269 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and 566 were assigned to the Clusters of Orthologous Groups (COG) functional classification "defense mechanisms". Analysis of differentially expressed genes (DEGs) showed that 2657 unigenes were differentially expressed between males and females (q < 0.1), and these DEGs were assigned to 45 KEGG pathways involving stress resistance, immunity, energy metabolism, reproduction, lifespan regulation, and diseases. Further analysis revealed that females might be more sensitive to stress through upregulation of c-Jun N-terminal kinases (JNKs) and 5'AMP-activated protein kinase (AMPK), and were also possibly more sensitive to dynamic changes in energy. Females expressed higher major histocompatibility complex (MHC) class II levels than males, enhancing resistance to pathogens, and the DEGs related to reproduction included MAPK, CaMK, CPEB, and Cdc25. The genes related to stress, energy, and immunity were also likely related to the regulation of longevity. The upregulated JNKs in females might prolong lifespan and relieve antioxidant stress. Females may also activate the AMPK pathway and implement dietary restrictions to prolong lifespan, whereas males may upregulate SIRT1 and CRAB to increase lifespan. CONCLUSIONS: Female magpies might be more sensitive to stress and dynamic changes in energy thus enhanced resistance to pathogens, and the genes related to stress, energy, and immunity were also possibly related to the regulation of longevity. Further confirmations with techniques such as RT-qPCR and western blot are necessary to validate the above arguments.

Blood Gene Signatures of Chagas Cardiomyopathy With or Without Ventricular Dysfunction.

Chagas disease, caused by the protozoan parasite Trypanosoma cruzi, affects 7 million people in Latin American areas of endemicity. About 30% of infected patients will develop chronic Chagas cardiomyopathy (CCC), an inflammatory cardiomyopathy characterized by hypertrophy, fibrosis, and myocarditis. Further studies are necessary to understand the molecular mechanisms of disease progression. Transcriptome analysis has been increasingly used to identify molecular changes associated with disease outcomes. We thus assessed the whole-blood transcriptome of patients with Chagas disease. Microarray analysis was performed on blood samples from 150 subjects, of whom 30 were uninfected control patients and 120 had Chagas disease (1 group had asymptomatic disease, and 2 groups had CCC with either a preserved or reduced left ventricular ejection fraction [LVEF]). Each Chagas disease group displayed distinct gene expression and functional pathway profiles. The most different expression patterns were between CCC groups with a preserved or reduced LVEF. A more stringent analysis indicated that 27 differentially expressed genes, particularly those related to natural killer (NK)/CD8+ T-cell cytotoxicity, separated the 2 groups. NK/CD8+ T-cell cytotoxicity could play a role in determining Chagas disease progression. Understanding genes associated with disease may lead to improved insight into CCC pathogenesis and the identification of prognostic factors for CCC progression.

Characterization of the peripheral blood transcriptome and adaptive evolution of the MHC I and TLR gene families in the wolf (Canis lupus).

BACKGROUND: The wolf (Canis lupus) is one of the most widely distributed terrestrial mammals, because it is well adapted to various ecological niches and their corresponding pathogen environments. Immunological competence is a crucial factor involved in adapting to a changing environment and fighting pathogen infection in animals. In this study, the peripheral blood transcriptome of wolves was generated via RNA-seq to advance understanding of the wolf immunome, with a special focus on the major histocompatibility complex class I (MHC I) and toll-like receptor (TLR) gene families, which are involved in pathogen recognition and defense. RESULTS: The blood transcriptomic libraries of eight wolves originating from Tibet and Inner Mongolia were sequenced, and approximately 383 million reads were generated. Using a genome-guided assembly strategy, we obtained 123,851 unigenes, with a mean length of 845 bp and an N50 length of 1121 bp. On the basis of BLAST searches against the NCBI non-redundant protein database (Nr), a total of 36,192 (29.22%) unigenes were annotated. For functional classification, 24,663 unigenes were assigned to 13,016 Gene Ontology (GO) terms belonging to 51 sub-categories of the three main GO categories. Additionally, 7682 unigenes were classified into 6 Kyoto Encyclopedia of Genes and Genomes (KEGG) categories, in which the most represented functional sub-categories were signal transduction and the immune system, and 16,238 unigenes were functionally classified into 25 Eukaryotic Orthologous Groups (KOG) categories. We observed an overall higher ω (d N/d S) value at antigen-binding sites (ABSs) than at non-ABS regions as well as clear evidence of intergenic/intragenic recombination events at wolf MHC I loci. Additionally, our analysis revealed that carnivorous TLRs were dominated by purifying selection, with mean ω values at each TLR locus ranging from 0.173 to 0.527. However, we also found significant instances of positive selection that acted on several codons in pathogen recognition domains and were linked to species-specific differences in pathogen recognition. CONCLUSIONS: This study represents the first attempt to characterize the blood transcriptome of the wolf and to highlight the value of investigating the immune system. Balancing selection and recombination have contributed to the historical evolution of wolf MHC I genes. Moreover, TLRs in carnivores have undergone adaptive evolution against the background of purifying selection, and a high level of adaptive evolution was detected in the wolf TLR system.

Transcriptome profiling of Arabian horse blood during training regimens.

BACKGROUND: Arabian horses are believed to be one of the oldest and most influential horse breeds in the world. Blood is the main tissue involved in maintaining body homeostasis, and it is considered a marker of the processes taking place in the other tissues. Thus, the aim of our study was to identify the genetic basis of changes occurring in the blood of Arabian horses subjected to a training regimen and to compare the global gene expression profiles between different training periods (T1: after a slow canter phase that is considered a conditioning phase, T2: after an intense gallop phase, and T3: at the end of the racing season) and between trained and untrained horses (T0). RNA sequencing was performed on 37 samples with a 75-bp single-end run on a HiScanSQ platform (Illumina), and differentially expressed genes (DEGs) were identified based on DESeq2 (v1.11.25) software. RESULTS: An increase in the number of DEGs between subsequent training periods was observed, and the highest amount of DEGs (440) was detected between untrained horses (T0) and horses at the end of the racing season (T3). The comparisons of the T2 vs. T3 transcriptomes and the T0 vs. T3 transcriptomes showed a significant gain of up-regulated genes during long-term exercise (up-regulation of 266 and 389 DEGs in the T3 period compared to T2 and T0, respectively). Forty differentially expressed genes were detected between the T1 and T2 periods, and 296 between T2 and T3. Functional annotation showed that the most abundant genes up-regulated in exercise were involved in pathways regulating cell cycle (PI3K-Akt signalling pathway), cell communication (cAMP-dependent pathway), proliferation, differentiation and apoptosis, as well as immunity processes (Jak-STAT signalling pathway). CONCLUSIONS: We investigated whether training causes permanent transcriptome changes in horse blood as a reflection of adaptation to conditioning and the maintenance of fitness to compete in flat races. The present study identified the overrepresented molecular pathways and genes that are essential for maintaining body homeostasis during long-term exercise in Arabian horses. Selected DEGs should be further investigated as markers that are potentially associated with racing performance in Arabian horses.

RNA-Seq analysis of seasonal and individual variation in blood transcriptomes of healthy managed bottlenose dolphins.

BACKGROUND: The blood transcriptome can reflect both systemic exposures and pathological changes in other organs of the body because immune cells recirculate through the blood, lymphoid tissues, and affected sites. In human and veterinary medicine, blood transcriptome analysis has been used successfully to identify markers of disease or pathological conditions, but can be confounded by large seasonal changes in expression. In comparison, the use of transcriptomic based analyses in wildlife has been limited. Here we report a longitudinal study of four managed bottlenose dolphins located in Waikoloa, Hawaii, serially sampled (approximately monthly) over the course of 1 year to establish baseline information on the content and variation of the dolphin blood transcriptome. RESULTS: Illumina based RNA-seq analyses were carried out using both the Ensembl dolphin genome and a de novo blood transcriptome as guides. Overall, the blood transcriptome encompassed a wide array of cellular functions and processes and was relatively stable within and between animals over the course of 1 year. Principal components analysis revealed moderate clustering by sex associated with the variation among global gene expression profiles (PC1, 22 % of variance). Limited seasonal change was observed, with < 2.5 % of genes differentially expressed between winter and summer months (FDR < 0.05). Among the differentially expressed genes, cosinor analysis identified seasonal rhythmicity for the observed changes in blood gene expression, consistent with studies in humans. While the proportion of seasonally variant genes in these dolphins is much smaller than that reported in humans, the majority of those identified in dolphins were also shown to vary with season in humans. Gene co-expression network analysis identified several gene modules with significant correlation to age, sex, or hematological parameters. CONCLUSIONS: This longitudinal analysis of healthy managed dolphins establishes a preliminary baseline for blood transcriptome analysis in this species. Correlations with hematological parameters, distinct from muted seasonal effects, suggest that the otherwise relatively stable blood transcriptome may be a useful indicator of health and exposure. A robust database of gene expression in free-ranging and managed dolphins across seasons with known adverse health conditions or contaminant exposures will be needed to establish predictive gene expression profiles suitable for biomonitoring.

Whole blood transcriptomic signature of Cushing's syndrome.

OBJECTIVE: Cushing's syndrome is characterized by high morbidity and mortality with high interindividual variability. Easily measurable biomarkers, in addition to the hormone assays currently used for diagnosis, could reflect the individual biological impact of glucocorticoids. The aim of this study is to identify such biomarkers through the analysis of whole blood transcriptome. DESIGN: Whole blood transcriptome was evaluated in 57 samples from patients with overt Cushing's syndrome, mild Cushing's syndrome, eucortisolism, and adrenal insufficiency. Samples were randomly split into a training cohort to set up a Cushing's transcriptomic signature and a validation cohort to assess this signature. METHODS: Total RNA was obtained from whole blood samples and sequenced on a NovaSeq 6000 System (Illumina). Both unsupervised (principal component analysis) and supervised (Limma) methods were used to explore the transcriptome profile. Ridge regression was used to build a Cushing's transcriptome predictor. RESULTS: The transcriptomic profile discriminated samples with overt Cushing's syndrome. Genes mostly associated with overt Cushing's syndrome were enriched in pathways related to immunity, particularly neutrophil activation. A prediction model of 1500 genes built on the training cohort demonstrated its discriminating value in the validation cohort (accuracy .82) and remained significant in a multivariate model including the neutrophil proportion (P = .002). Expression of FKBP5, a single gene both overexpressed in Cushing's syndrome and implied in the glucocorticoid receptor signaling, could also predict Cushing's syndrome (accuracy .76). CONCLUSIONS: Whole blood transcriptome reflects the circulating levels of glucocorticoids. FKBP5 expression could be a nonhormonal marker of Cushing's syndrome.

Comparative Blood Transcriptome Analysis of Semi-Natural and Controlled Environment Populations of Yangtze Finless Porpoise.

The Yangtze finless porpoises (Neophocaena asiaeorientalis asiaeorientalis) living in different environments display significant differences in behavior and physiology. To compare and analyze gene expression differences between an ex situ population and a controlled environment population of the Yangtze finless porpoise, we sequenced the transcriptome of blood tissues living in a semi-natural reserve and an artificial facility, respectively. We identified 6860 differentially expressed genes (DEGs), of which 6603 were up-regulated and 257 were down-regulated in the controlled environment vs ex situ comparison. GO and KEGG enrichment analysis showed that the up-regulated genes in the controlled environment population were significantly associated with glucose metabolism, amino acid metabolism, and the nervous system, while those up-regulated in the ex situ population were significantly associated with energy supply and biosynthesis. Further analysis showed that metabolic and hearing-related genes were significantly affected by changes in the environment, and key metabolic genes such as HK, PFK, IDH, and GLS and key hearing-related genes such as OTOA, OTOF, SLC38A1, and GABBR2 were identified. These results suggest that the controlled environment population may have enhanced glucose metabolic ability via activation of glycolysis/gluconeogenesis, the TCA cycle, and inositol phosphate metabolism, while the ex situ population may meet higher energy requirements via enhancement of the amino acid metabolism of the liver and muscle and oxidative phosphorylation. Additionally, the acoustic behavior and auditory-related genes of Yangtze finless porpoise may show responsive changes and differential expression under different environment conditions, and thus the auditory sensitivity may also show corresponding adaptive characteristics. This study provides a new perspective for further exploration of the responsive changes of the two populations to various environments and provides a theoretical reference for further improvements in conservation practices for the Yangtze finless porpoise.