Alternative splicing of pre-mRNA modulates the immune response in Holstein cattle naturally infected with Mycobacterium avium subsp. paratuberculosis

Little is known about the role of alternative splicing (AS) in regulating gene expression in Mycobacteria-infected individuals in distinct stages of infection. Pre-mRNA AS consists of the removal of introns and the assembly of exons contained in eukaryotic genes. AS events can influence transcript stability or structure with important physiological consequences. Using RNA-Seq data from peripheral blood (PB) and ileocecal valve (ICV) samples collected from Holstein cattle with focal and diffuse paratuberculosis (PTB)-associated histopathological lesions in gut tissues and without lesions (controls), we detected differential AS profiles between the infected and control groups. Four of the identified AS events were experimentally validated by reverse transcription-digital droplet PCR (RT-ddPCR). AS events in several genes correlated with changes in gene expression. In the ICV of animals with diffuse lesions, for instance, alternatively spliced genes correlated with changes in the expression of genes involved in endocytosis, antigen processing and presentation, complement activation, and several inflammatory and autoimmune diseases in humans. Taken together, our results identified common mechanisms of AS involvement in the pathogenesis of PTB and human diseases and shed light on novel diagnostic and therapeutic interventions to control these diseases.

Whole-genome long-read sequencing to unveil Enterococcus antimicrobial resistance in dairy cattle farms exposed a widespread occurrence of Enterococcus lactis

Enterococcus faecalis (Efs) and Enterococcus faecium (Efm) are major causes of multiresistant healthcare-associated or nosocomial infections. Efm has been traditionally divided into clades A (healthcare associated) and B (community associated) but clade B has been recently reassigned to Enterococcus lactis (Elc). However, identification techniques do not routinely differentiate Elc from Efm. As part of a longitudinal study to investigate the antimicrobial resistance of Enterococcus in dairy cattle, isolates initially identified as Efm were confirmed as Elc after Oxford-Nanopore long-fragment whole-genome sequencing and genome comparisons. An Efm-specific PCR assay was developed and used to identify isolates recovered from animal feces on five farms, resulting in 44 Efs, 23 Efm, and 59 Elc. Resistance, determined by broth microdilution, was more frequent in Efs than in Efm and Elc but all isolates were susceptible to ampicillin, daptomycin, teicoplanin, tigecycline, and vancomycin. Genome sequencing analysis of 32 isolates identified 23 antimicrobial resistance genes (ARGs, mostly plasmid-located) and 2 single nucleotide polymorphisms associated with resistance to 10 antimicrobial classes, showing high concordance with phenotypic resistance. Notably, linezolid resistance in Efm was encoded by the optrA gene, located in plasmids downstream of the fexA gene. Although most Elc lacked virulence factors and genetic determinants of resistance, one isolate carried a plasmid with eight ARGs. This study showed that Elc is more prevalent than Efm in dairy cattle but carries fewer ARGs and virulence genes. However, Elc can carry multi-drug-resistant plasmids like those harbored by Efm and could act as a donor of ARGs for other pathogenic enterococcal species.IMPORTANCEEnterococcus species identification is crucial due to differences in pathogenicity and antibiotic resistance profiles. The failure of traditional methods or whole-genome sequencing-based taxonomic classifiers to distinguish Enterococcus lactis (Elc) from Enterococcus faecium (Efm) results in a biased interpretation of Efm epidemiology. The Efm species-specific real-time PCR assay developed here will help to properly identify Efm (only the formerly known clade A) in future studies. Here, we showed that Elc is prevalent in dairy cattle, and although this species carries fewer genetic determinants of resistance (GDRs) than Enterococcus faecalis (Efs) and Efm, it can carry multi-drug-resistant (MDR) plasmids and could act as a donor of resistance genes for other pathogenic enterococcal species. Although all isolates (Efs, Efm, and Elc) were susceptible to critically or highly important antibiotics like daptomycin, teicoplanin, tigecycline, and vancomycin, the presence of GDRs in MDR-plasmids is a concern since antimicrobials commonly used in livestock could co-select and confer resistance to critically important antimicrobials not used in food-producing animals.

Summary-data based Mendelian randomization identifies gene expression regulatory polymorphisms associated with bovine paratuberculosis by modulation of the nuclear factor Kappa β (NF-κß)-mediated inflammatory response

Genome-wide association studies (GWAS) have identified host genetic variants associated with paratuberculosis (PTB) susceptibility. Most of the GWAS-identified SNPs are in non-coding regions. Connecting these non-coding variants and downstream affected genes is a challenge and, up to date, only a few functional mutations or expression quantitative loci (cis-eQTLs) associated with PTB susceptibility have been identified. In the current study, the associations between imputed whole-genome sequence genotypes and whole RNA-Sequencing data from peripheral blood (PB) and ileocecal valve (ICV) samples of Spanish Holstein cows (N = 16) were analyzed with TensorQTL. This approach allowed the identification of 88 and 37 cis-eQTLs regulating the expression levels of 90 and 37 genes in PB and ICV samples, respectively (False discorey rate, FDR ≤ 0.05). Next, we applied summary-based data Mendelian randomization (SMR) to integrate the cis-eQTL dataset with GWAS data obtained from a cohort of 813 culled cattle that were classified according to the presence or absence of PTB-associated histopathological lesions in gut tissues. After multiple testing corrections (FDR ≤ 0.05), we identified two novel cis-eQTLs affecting the expression of the early growth response factor 4 (EGR4) and the bovine neuroblastoma breakpoint family member 6-like protein isoform 2 (MGC134040) that showed pleiotropic associations with the presence of multifocal and diffuse lesions in gut tissues; P = 0.002 and P = 0.017, respectively. While EGR4 acts as a brake on T-cell proliferation and cytokine production through interaction with the nuclear factor Kappa β (NF-κß), MGC134040 is a target gene of NF-κß. Our findings provide a better understanding of the genetic factors influencing PTB outcomes, confirm that the multifocal lesions are localized/confined lesions that have different underlying host genetics than the diffuse lesions, and highlight regulatory SNPs and regulated-gene targets to design future functional studies.

Monitoring within-farm transmission dynamics of antimicrobial-resistant Campylobacter in dairy cattle using broth microdilution and long-read whole genome sequencing

Campylobacter jejuni and Campylobacter coli are important foodborne zoonotic pathogens and cause for concern due to the increasing trend in antimicrobial resistance. A long-run surveillance study was conducted in animals from different age groups in five dairy cattle farms to investigate the within-farm diversity and transmission dynamics of resistant Campylobacter throughout time. The resistance phenotype of the circulating isolates (170 C. jejuni and 37 C. coli) was determined by broth microdilution and a selection of 56 isolates were whole genome sequenced using the Oxford-Nanopore long-fragment sequencing technology resulting in completely resolved and circularized genomes (both chromosomes and plasmids). C. jejuni was isolated from all farms while C. coli was isolated from only two farms, but resistance rates were higher in C. coli than in C. jejuni and in calves than in adult animals. Some genotypes (e.g. ST-48, gyrA_T86I/tet(O)/blaOXA-61 in farm F1; ST-12000, aadE-Cc/tet(O)/blaOXA-489 in F4) persisted throughout the study while others were only sporadically detected. Acquisition of extracellular genes from other isolates and intracellular mutational events were identified as the processes that led to the emergence of the resistant genotypes that spread within the herds. Monitoring with Oxford Nanopore Technologies sequencing helped to decipher the complex molecular epidemiology underlying the within-farm dissemination of resistant Campylobacter.

Within-farm dynamics of ESBL-producing Escherichia coli in dairy cattle: Resistance profiles and molecular characterization by long-read whole-genome sequencing

A longitudinal study was designed in five dairy cattle farms to assess the within-farm dynamics of ESBL-/AmpC-/carbapenemase-producing E. coli and their resistance profiles, along with the genes conferring the resistance phenotypes. Twelve samplings were performed over a period of 16 months, collecting rectal feces from apparently healthy animals in three age groups (calves, heifers, and lactating cows) that were subjected to selective isolation in cefotaxime-containing media. Minimum inhibitory concentrations were determined by broth microdilution for 197 cefotaxime-resistant E. coli (1–3 isolates per age group and sampling date), and 41 of them were selected for long-read whole-genome sequencing. Cefotaxime-resistant E. coli were detected in the five farms, but isolation frequency and resistance profiles varied among farms and age groups. The genetic profiling of a selection of isolates recovered in two of the farms was described in full detail, showing the predominance of a few genomic subtypes of E. coli in one farm (F1) and great variability of strains in another one (F4). Two predominant distinct strains carrying the blaCTX-M-1 gene in IncX1 plasmids successively spread and persisted in F1 over a prolonged period. In F4, 13 different MLST types carrying a high diversity of ESBL-encoding genes in 6 different plasmid types were observed, probably as the result of multiple source contamination events. In both farms, the presence of certain plasmid types with the same repertoire of ARGs in different E. coli STs strongly suggested the occurrence of horizontal transfer of such plasmids among strains circulating within the farms. Considering the public health importance of ESBL-producing E. coli both as pathogens and as vectors for resistance mechanisms, the presence of β-lactamase- and other AMR-encoding genes in plasmids that can be readily transferred between bacteria is a concern that highlights the need for One Health surveillance.

Reduction of the resistome risk from cow slurry and manure microbiomes to soil and vegetable microbiomes

In cow farms, the interaction between animal and environmental microbiomes creates hotspots for antibiotic resistance dissemination. A shotgun metagenomic approach was used to survey the resistome risk in five dairy cow farms. To this purpose, 10 environmental compartments were sampled: 3 of them linked to productive cows (fresh slurry, stored slurry, slurry-amended pasture soil); 6 of them to non-productive heifers and dry cows (faeces, fresh manure, aged manure, aged manure-amended orchard soil, vegetables-lettuces and grazed soil); and, finally, unamended control soil. The resistome risk was assessed using MetaCompare, a computational pipeline which scores the resistome risk according to possible links between antibiotic resistance genes (ARGs), mobile genetic elements (MGEs) and human pathogens. The resistome risk decreased from slurry and manure microbiomes to soil and vegetable microbiomes. In total (sum of all the compartments), 18,157 ARGs were detected: 24% related to ansamycins, 21% to multidrugs, 14% to aminoglycosides, 12% to tetracyclines, 9% to β-lactams, and 9% to macrolide-lincosamide-streptogramin B. All but two of the MGE-associated ARGs were only found in the animal dejections (not in soil or vegetable samples). Several ARGs with potential as resistome risk markers (based on their presence in hubs of co-occurrence networks and high dissemination potential) were identified. As a precautionary principle, improved management of livestock dejections is necessary to minimize the risk of antibiotic resistance.

Whole genome-based characterisation of antimicrobial resistance and genetic diversity in Campylobacter jejuni and Campylobacter coli from ruminants

Campylobacter, a leading cause of gastroenteritis in humans, asymptomatically colonises the intestinal tract of a wide range of animals.Although antimicrobial treatment is restricted to severe cases, the increase of antimicrobial resistance (AMR) is a concern. Considering the significant contribution of ruminants as reservoirs of resistant Campylobacter, Illumina whole-genome sequencing was used to characterise the mechanisms of AMR in Campylobacter jejuni and Campylobacter coli recovered from beef cattle, dairy cattle, and sheep in northern Spain. Genome analysis showed extensive genetic diversity that clearly separated both species. Resistance genotypes were identified by screening assembled sequences with BLASTn and ABRicate, and additional sequence alignments were performed to search for frameshift mutations and gene modifications. A high correlation was observed between phenotypic resistance to a given antimicrobial and the presence of the corresponding known resistance genes. Detailed sequence analysis allowed us to detect the recently described mosaic tet(O/M/O) gene in one C. coli, describe possible new alleles of blaOXA-61-like genes, and decipher the genetic context of aminoglycoside resistance genes, as well as the plasmid/chromosomal location of the different AMR genes and their implication for resistance spread. Updated resistance gene databases and detailed analysis of the matched open reading frames are needed to avoid errors when using WGS-based analysis pipelines for AMR detection in the absence of phenotypic data.

Peripheral blood mononuclear cells (PBMC) microbiome is not affected by colon microbiota in healthy goats

Background

The knowledge about blood circulating microbiome and its functional relevance in healthy individuals remains limited. An assessment of changes in the circulating microbiome was performed by sequencing peripheral blood mononuclear cells (PBMC) bacterial DNA from goats supplemented or not in early life with rumen liquid transplantation.

Results

Most of the bacterial DNA associated to PBMC was identified predominantly as Proteobacteria (55%) followed by Firmicutes (24%), Bacteroidetes (11%) and Actinobacteria (8%). The predominant genera found in PBMC samples were Pseudomonas, Prevotella, Sphingomonas, Acinetobacter, Corynebacterium and Ruminococcus. Other genera such as Butyrivibrivio, Bifidobacterium, Dorea and Coprococcus were also present in lower proportions. Several species known as blood pathogens or others involved in gut homeostasis such as Faecalibacterium prausnitzii were also identified. However, the PBMC microbiome phylum composition differed from that in the colon of goats (P ≤ 0.001), where Firmicutes was the predominant phylum (83%). Although, rumen liquid administration in early-life altered bacterial community structure and increased Tlr5 expression (P = 0.020) in colon pointing to higher bacterial translocation, less than 8% of OTUs in colon were also observed in PBMCs.

Conclusions

Data suggest that in physiological conditions, PBMC microbiome differs from and is not affected by colon gut microbiota in small ruminants. Although, further studies with larger number of animals and covering other animal tissues are required, results point to a common circulating bacterial profile on mammals being phylum Proteobacteria, and genera Pseudomonas and Prevotella the most abundants. All suggest that PBMC microbiome in healthy ruminants could be implicated in homeostatic condition. This study expands our knowledge about PBMC microbiome contribution to health in farm animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-021-00091-7.

Borrelia burgdorferi infection induces long-term memory-like responses in macrophages with tissue-wide consequences in the heart

Lyme carditis is an extracutaneous manifestation of Lyme disease characterized by episodes of atrioventricular block of varying degrees and additional, less reported cardiomyopathies. The molecular changes associated with the response to Borrelia burgdorferi over the course of infection are poorly understood. Here, we identify broad transcriptomic and proteomic changes in the heart during infection that reveal a profound down-regulation of mitochondrial components. We also describe the long-term functional modulation of macrophages exposed to live bacteria, characterized by an augmented glycolytic output, increased spirochetal binding and internalization, and reduced inflammatory responses. In vitro, glycolysis inhibition reduces the production of tumor necrosis factor (TNF) by memory macrophages, whereas in vivo, it produces the reversion of the memory phenotype, the recovery of tissue mitochondrial components, and decreased inflammation and spirochetal burdens. These results show that B. burgdorferi induces long-term, memory-like responses in macrophages with tissue-wide consequences that are amenable to be manipulated in vivo.

Lyme carditis is a manifestation of Lyme disease characterized by episodes of atrioventricular block and additional cardiomyopathies. This study describes the proteomic and transcriptomic changes in the heart upon infection with Borrelia burgdorferi, and identifies innate immune memory hallmarks specific to the response to the spirochete that are amenable to therapeutic manipulation.

Intestinal epithelial deletion of the glucocorticoid receptor NR3C1 alters expression of inflammatory mediators and barrier function

Glucocorticoids (GCs) are important hormones involved in the regulation of multiple physiologic functions. GCs are also widely used in anti-inflammatory/immunosuppressant drugs. GCs are synthesized by the adrenal cortex as part of the hypothalamus-pituitary-adrenal axis and also by intestinal epithelial cells, among other peripheral sites. GCs are one of the main therapy choices for the exacerbations of inflammatory bowel disease, but they are not useful to prolong remission, and development of tolerance with secondary treatment failure is frequent. Thus, GC actions at the intestinal epithelial level are of great importance, both physiologically and pharmacologically. We generated a tamoxifen-inducible nuclear receptor subfamily 3 group C member 1 (NR3C1)^ΔIEC mouse model to study the effects of GCs on epithelial cells in vivo. Nr3c1 deletion in epithelial cells of the small intestine and colon was associated with limited colonic inflammation at 1 wk postdeletion, involving augmented epithelial proliferation and mucus production, plus local and systemic immune/inflammatory changes. This phenotype regressed substantially, but not completely, after 2 wk. The mechanism may involve augmented inflammatory signaling by epithelial cells or defective barrier function. We conclude that the epithelial GC receptor plays a significant role in colonic homeostasis in basal conditions, but its deficiency can be compensated in the short term. Future studies are required to assess the impact of Nr3c1 deletion in other conditions such as experimental colitis.-Aranda, C. J., Arredondo-Amador, M., Ocón, B., Lavín, J. L., Aransay, A. M., Martínez-Augustin, O., Sánchez de Medina, F. Intestinal epithelial deletion of the glucocorticoid receptor NR3C1 alters expression of inflammatory mediators and barrier function.

MiR-873-5p acts as an epigenetic regulator in early stages of liver fibrosis and cirrhosis

Glycine N-methyltransferase (GNMT) is the most abundant methyltransferase in the liver and a master regulator of the transmethylation flux. GNMT downregulation leads to loss of liver function progressing to fibrosis, cirrhosis, and hepatocellular carcinoma. Moreover, GNMT deficiency aggravates cholestasis-induced fibrogenesis. To date, little is known about the mechanisms underlying downregulation of GNMT levels in hepatic fibrosis and cirrhosis. On this basis, microRNAs are epigenetic regulatory elements that play important roles in liver pathology. In this work, we aim to study the regulation of GNMT by microRNAs during liver fibrosis and cirrhosis. Luciferase assay on the 3ʹUTR-Gnmt was used to confirm in silico analysis showing that GNMT is potentially targeted by the microRNA miR-873-5p. Correlation between GNMT and miR-873-5p in human cholestasis and cirrhosis together with miR-873-5p inhibition in vivo in different mouse models of liver cholestasis and fibrosis [bile duct ligation and Mdr2 (Abcb4)^-/- mouse] were then assessed. The analysis of liver tissue from cirrhotic and cholestatic patients, as well as from the animal models, showed that miR-873-5p inversely correlated with the expression of GNMT. Importantly, high circulating miR-873-5p was also detected in cholestastic and cirrhotic patients. Preclinical studies with anti-miR-873-5p treatment in bile duct ligation and Mdr2^-/- mice recovered GNMT levels in association with ameliorated inflammation and fibrosis mainly by counteracting hepatocyte apoptosis and cholangiocyte proliferation. In conclusion, miR-873-5p emerges as a novel marker for liver fibrosis, cholestasis, and cirrhosis and therapeutic approaches based on anti-miR-873-5p may be effective treatments for liver fibrosis and cholestatic liver disease.

Identification of a highly active tannase enzyme from the oral pathogen Fusobacterium nucleatum subsp. polymorphum

BACKGROUND

Tannases are tannin-degrading enzymes that have been described in fungi and bacteria as an adaptative mechanism to overcome the stress conditions associated with the presence of these phenolic compounds.

RESULTS

We have identified and expressed in E. coli a tannase from the oral microbiota member Fusobacterium nucleatum subs. polymorphum (TanB_Fnp). TanB_Fnp is the first tannase identified in an oral pathogen. Sequence analyses revealed that it is closely related to other bacterial tannases. The enzyme exhibits biochemical properties that make it an interesting target for industrial use. TanB_Fnp has one of the highest specific activities of all bacterial tannases described to date and shows optimal biochemical properties such as a high thermal stability: the enzyme keeps 100% of its activity after prolonged incubations at different temperatures up to 45 °C. TanB_Fnp also shows a wide temperature range of activity, maintaining above 80% of its maximum activity between 22 and 55 °C. The use of a panel of 27 esters of phenolic acids demonstrated activity of TanB_Fnp only against esters of gallic and protocatechuic acid, including tannic acid, gallocatechin gallate and epigallocatechin gallate. Overall, TanB_Fnp possesses biochemical properties that make the enzyme potentially useful in biotechnological applications.

CONCLUSIONS

We have identified and characterized a metabolic enzyme from the oral pathogen Fusobacterium nucleatum subsp. polymorphum. The biochemical properties of TanB_Fnp suggest that it has a major role in the breakdown of complex food tannins during oral processing. Our results also provide some clues regarding its possible participation on bacterial survival in the oral cavity. Furthermore, the characteristics of this enzyme make it of potential interest for industrial use.

SOX17 regulates cholangiocyte differentiation and acts as a tumor suppressor in cholangiocarcinoma

BACKGROUND & AIMS

Cholangiocarcinoma (CCA) is a biliary malignancy linked to genetic and epigenetic abnormalities, such as hypermethylation of SOX17 promoter. Here, the role of SOX17 in cholangiocyte differentiation and cholangiocarcinogenesis was studied.

METHODS

SOX17 expression/function was evaluated along the differentiation of human induced pluripotent stem cells (iPSC) into cholangiocytes, in the dedifferentiation process of normal human cholangiocytes (NHC) in culture and in cholangiocarcinogenesis. Lentiviruses for SOX17 overexpression or knockdown were used. Gene expression and DNA methylation profiling were performed.

RESULTS

SOX17 expression is induced in the last stage of cholangiocyte differentiation from iPSC and regulates the acquisition of biliary markers. SOX17 becomes downregulated in NHC undergoing dedifferentiation; experimental SOX17 knockdown in differentiated NHC downregulated biliary markers and promoted baseline and Wnt-dependent proliferation. SOX17 expression is lower in human CCA than in healthy tissue, which correlates with worse survival after tumor resection. In CCA cells, SOX17 overexpression decreased their tumorigenic capacity in murine xenograft models, which was related to increased oxidative stress and apoptosis. In contrast, SOX17 overexpression in NHC did not affect their survival but inhibited their baseline proliferation. In CCA cells, SOX17 inhibited migration, anchorage-independent growth and Wnt/β-catenin-dependent proliferation, and restored the expression of biliary markers and primary cilium length. In human CCA, SOX17 promoter was found hypermethylated and its expression inversely correlates with the methylation grade. In NHC, Wnt3a decreased SOX17 expression in a DNMT-dependent manner, whereas in CCA, DNMT1 inhibition or silencing upregulated SOX17.

CONCLUSIONS

SOX17 regulates the differentiation and maintenance of the biliary phenotype and functions as a tumor suppressor for CCA, being a potential prognostic marker and a promising therapeutic target.

LAY SUMMARY

Understanding the molecular mechanisms involved in the pathogenesis of CCA is key in finding new valuable diagnostic and prognostic biomarkers, as well as therapeutic targets. This study provides evidence that SOX17 regulates the differentiation and maintenance of the biliary phenotype, and its downregulation promotes their tumorigenic transformation. SOX17 acts as a tumor suppressor in CCA and its genetic, molecular and/or pharmacological restoration may represent a new promising therapeutic strategy. Moreover, SOX17 expression correlates with the outcome of patients after tumor resection, being a potential prognostic biomarker.

A synbiotic composed of Lactobacillus fermentum CECT5716 and FOS prevents the development of fatty acid liver and glycemic alterations in rats fed a high fructose diet associated with changes in the microbiota

We investigated the effect of a high fructose diet (HFD) on Sprague Dawley rats and the impact of a synbiotic composed of Lactobacillus fermentum CECT5716 and fructooligosaccharides. Feeding the HFD for 5 weeks resulted in liver steatosis and insulin resistance but not obesity. These changes were associated with increased production of short-chain fatty acids and increased Bacteroidetes in feces, with an augmented Bacteroidetes/Firmicutes ratio, among other changes in the microbiota. In addition, barrier function was weakened, with increased LPS plasma levels. These data are consistent with increased fructose availability in the distal gut due to saturation of absorptive mechanisms, leading to dysbiosis, endotoxemia, hepatic steatosis, and insulin resistance. Treatment with the synbiotic prevented some of the pathological effects, so that treated rats did not develop steatosis or systemic inflammation, while dysbiosis and barrier function were greatly ameliorated. In addition, the synbiotic had hypolipidemic effects. The synbiotic composed by L. fermentum CECT5716 and fructooligosaccharides has beneficial effects in a model of metabolic syndrome induced by a HFD, suggesting it might be clinically useful in this type of condition, particularly considering that high fructose intake has been related to metabolic syndrome in humans.

Genomics and transcriptomics characterization of genes expressed during postharvest at 4°C by the edible basidiomycete Pleurotus ostreatus

Pleurotus ostreatus is an industrially cultivated basidiomycete with nutritional and environmental applications. Its genome, which was sequenced by the Joint Genome Institute, has become a model for lignin degradation and for fungal genomics and transcriptomics studies. The complete P. ostreatus genome contains 35 Mbp organized in 11 chromosomes, and two different haploid genomes have been individually sequenced. In this work, genomics and transcriptomics approaches were employed in the study of P. ostreatus under different physiological conditions. Specifically, we analyzed a collection of expressed sequence tags (EST) obtained from cut fruit bodies that had been stored at 4°C for 7 days (postharvest conditions). Studies of the 253 expressed clones that had been automatically and manually annotated provided a detailed picture of the life characteristics of the self-sustained fruit bodies. The results suggested a complex metabolism in which autophagy, RNA metabolism, and protein and carbohydrate turnover are increased. Genes involved in environment sensing and morphogenesis were expressed under these conditions. The data improve our understanding of the decay process in postharvest mushrooms and highlight the use of high-throughput techniques to construct models of living organisms subjected to different environmental conditions.