CRISPR-Cas in Pseudomonas aeruginosa provides transient population-level immunity against high phage exposures.

The prokaryotic adaptive immune system, CRISPR-Cas (clustered regularly interspaced short palindromic repeats; CRISPR-associated), requires the acquisition of spacer sequences that target invading mobile genetic elements such as phages. Previous work has identified ecological variables that drive the evolution of CRISPR-based immunity of the model organism Pseudomonas aeruginosa PA14 against its phage DMS3vir, resulting in rapid phage extinction. However, it is unclear if and how stable such acquired immunity is within bacterial populations, and how this depends on the environment. Here, we examine the dynamics of CRISPR spacer acquisition and loss over a 30-day evolution experiment and identify conditions that tip the balance between long-term maintenance of immunity versus invasion of alternative resistance strategies that support phage persistence. Specifically, we find that both the initial phage dose and reinfection frequencies determine whether or not acquired CRISPR immunity is maintained in the long term, and whether or not phage can coexist with the bacteria. At the population genetics level, emergence and loss of CRISPR immunity are associated with high levels of spacer diversity that subsequently decline due to invasion of bacteria carrying pilus-associated mutations. Together, these results provide high resolution of the dynamics of CRISPR immunity acquisition and loss and demonstrate that the cumulative phage burden determines the effectiveness of CRISPR over ecologically relevant timeframes.

The Interfascicular Matrix of Energy Storing Tendons Houses Heterogenous Cell Populations Disproportionately Affected by Aging.

Energy storing tendons such as the human Achilles and equine superficial digital flexor tendon (SDFT) are prone to injury, with incidence increasing with aging, peaking in the 5th decade of life in the human Achilles tendon. The interfascicular matrix (IFM), which binds tendon fascicles, plays a key role in energy storing tendon mechanics, and aging alterations to the IFM negatively impact tendon function. While the mechanical role of the IFM in tendon function is well-established, the biological role of IFM-resident cell populations remains to be elucidated. Therefore, the aim of this study was to identify IFM-resident cell populations and establish how these populations are affected by aging. Cells from young and old SDFTs were subjected to single cell RNA-sequencing, and immunolabelling for markers of each resulting population used to localise cell clusters. Eleven cell clusters were identified, including tenocytes, endothelial cells, mural cells, and immune cells. One tenocyte cluster localised to the fascicular matrix, whereas nine clusters localised to the IFM. Interfascicular tenocytes and mural cells were preferentially affected by aging, with differential expression of genes related to senescence, dysregulated proteostasis and inflammation. This is the first study to establish heterogeneity in IFM cell populations, and to identify age-related alterations specific to IFM-localised cells.

Further Step in the Transition from Conventional Plasticizers to Versatile Bioplasticizers Obtained by the Valorization of Levulinic Acid and Glycerol.

In the last two decades, the use of phthalates has been restricted worldwide due to their well-known toxicity. Nonetheless, phthalates are still widely used for their versatility, high plasticization effect, low cost, and lack of valuable alternatives. This study presents the fully bio-based and versatile glycerol trilevulinate plasticizer (GT) that was obtained by the valorization of glycerol and levulinic acid. The mild-conditions and solvent-free esterification used to synthesize GT was optimized by investigating the product by Fourier transform infrared and NMR spectroscopy. An increasing content of GT, from 10 to 40 parts by weight per hundred parts of resin (phr), was tested with poly(vinyl chloride), poly(3-hydroxybutyrate), poly(3-hydroxybutyrate-co-3-hydroxyvalerate), poly(lactic acid), and poly(caprolactone), which typically present complicated processability and/or mechanical properties. GT produced a significant plasticization effect on both amorphous and semicrystalline polymers, reducing their glass-transition temperature and stiffness, as observed by differential scanning calorimetry measurements and tensile tests. Remarkably, GT also decreased both the melting temperature and crystallinity degree of semicrystalline polymers. Furthermore, GT underwent enzyme-mediated hydrolysis to its initial constituents, envisioning a promising prospective for environmental safety and upcycling. Furthermore, 50% inhibitory concentration (IC50) tests, using mouse embryo fibroblasts, proved that GT is an unharmful alternative plasticizer, which makes it potentially applicable in the biomedical field.

Impact of Interleukin 10 Deficiency on Intestinal Epithelium Responses to Inflammatory Signals.

Interleukin 10 (IL-10) is a pleiotropic, anti-inflammatory cytokine that has a major protective role in the intestine. Although its production by cells of the innate and adaptive immune system has been extensively studied, its intrinsic role in intestinal epithelial cells is poorly understood. In this study, we utilised both ATAC sequencing and RNA sequencing to define the transcriptional response of murine enteroids to tumour necrosis factor (TNF). We identified that the key early phase drivers of the transcriptional response to TNF within intestinal epithelium were NFκB transcription factor dependent. Using wild-type and Il10-/- enteroid cultures, we showed an intrinsic, intestinal epithelium specific effect of IL-10 deficiency on TNF-induced gene transcription, with significant downregulation of identified NFκB target genes Tnf, Ccl20, and Cxcl10, and delayed overexpression of NFκB inhibitor encoding genes, Nfkbia and Tnfaip3. IL-10 deficiency, or immunoblockade of IL-10 receptor, impacted on TNF-induced endogenous NFκB activity and downstream NFκB target gene transcription. Intestinal epithelium-derived IL-10 appears to play a crucial role as a positive regulator of the canonical NFκB pathway, contributing to maintenance of intestinal homeostasis. This is particularly important in the context of an inflammatory environment and highlights the potential for future tissue-targeted IL-10 therapeutic intervention.

Inter-kingdom relationships in Crohn's disease explored using a multi-omics approach.

The etiology of Crohn's disease (CD) is multifactorial. Bacterial and fungal microbiota are involved in the onset and/or progression of the disease. A bacterial dysbiosis in CD patients is accepted; however, less is known about the mycobiome and the relationships between the two communities. We investigated the interkingdom relationships, their metabolic consequences, and the changes in the fungal community during relapse and remission in CD.Two cohorts were evaluated: a British cohort (n = 63) comprising CD and ulcerative colitis patients, and controls. The fungal and bacterial communities of biopsy and fecal samples were analyzed, with the fecal volatiles; datasets were also integrated; and a Dutch cohort (n = 41) comprising CD patients and healthy controls was analyzed for stability of the gut mycobiome.A dysbiosis of the bacterial community was observed in biopsies and stool. Results suggest Bacteroides is likely key in CD and may modulate Candida colonization. A dysbiosis of the fungal community was observed only in the Dutch cohort; Malassezia and Candida were increased in patients taking immunosuppressants. Longitudinal analysis showed an increase in Cyberlindnera in relapse. Saccharomyces was dominant in all fecal samples, but not in biopsies, some of which did not yield fungal reads; amino acid degradation was the main metabolic change associated with CD and both bacteria and fungi might be implicated.We have shown that Bacteroides and yeasts may play a role in CD; understanding their role and relationship in the disease would shed new light on the development and treatment of CD.

Phage gene expression and host responses lead to infection-dependent costs of CRISPR immunity.

CRISPR-Cas immune systems are widespread in bacteria and archaea, but not ubiquitous. Previous work has demonstrated that CRISPR immunity is associated with an infection-induced fitness cost, which may help explain the patchy distribution observed. However, the mechanistic basis of this cost has remained unclear. Using Pseudomonas aeruginosa PA14 and its phage DMS3vir as a model, we perform a 30-day evolution experiment under phage mediated selection. We demonstrate that although CRISPR is initially selected for, bacteria carrying mutations in the phage receptor rapidly invade the population following subsequent reinfections. We then test three potential mechanisms for the observed cost of CRISPR: (1) autoimmunity from the acquisition of self-targeting spacers, (2) immunopathology or energetic costs from increased cas gene expression and (3) toxicity caused by phage gene expression prior to CRISPR-mediated cleavage. We find that phages can express genes before the immune system clears the infection and that expression of these genes can have a negative effect on host fitness. While infection does not lead to increased expression of cas genes, it does cause differential expression of multiple other host processes that may further contribute to the cost of CRISPR immunity. In contrast, we found little support for infection-induced autoimmunological and immunopathological effects. Phage gene expression prior to cleavage of the genome by the CRISPR-Cas immune system is therefore the most parsimonious explanation for the observed phage-induced fitness cost.

Maternal influences on oral and faecal microbiota maturation in neonatal calves in beef and dairy production systems.

BACKGROUND: The dam is considered an important source of microbes for the calf; consequently, the development of calf microbiota may vary with farming system due to differences between the contact the calf has with the dam. The objective of this study was to characterise the early changes in the composition of oral and faecal microbiota in beef and dairy calves (N = 10) using high-throughput sequencing of the 16S rRNA gene. The microbiota of calves was compared to selected anatomical niches on their dams which were likely to contribute to the vertical transfer of microbes. RESULTS: A total of 14,125 amplicon sequence variants (ASVs) were identified and taxonomically assigned. The oral microbiota of calves and their dams were composed of more similar microbes after the first 4 weeks of life than immediately after calving. The faecal microbiota of four-week old calves was composed of microbes which were more similar to those found in the oral microbiota of calves and adult cows than the faecal microbiota of adult cows. Specific ASVs were identified in the oral microbiota of four-week old calves that were also present in cow niches at calving, whereas very few ASVs were present in the calf faecal microbiota at four-weeks of age were present in any adult cow niche at calving. These results were observed in both beef and dairy calves. CONCLUSIONS: We did not observe any marked differences in the maturation of the oral and faecal microbiota between beef or dairy calves, despite dairy calves having very limited contact with their dam. This suggests the development of gastrointestinal microbiota in calves may not be affected by continued vertical transmission of microbes from the dam. Although the calf faecal microbiota changed over the first four-weeks of life, it was composed of microbes which were phylogenetically closer to those in the oral microbiota of calves and adult cows than the faeces of adult cows. There was little evidence of persistent microbial seeding of the calf faeces from anatomical niches on the cow at calving in either beef or dairy animals.

Characterization of Circulating Clostridium difficile Strains, Host Response and Intestinal Microbiome in Hospitalized Children With Diarrhea.

BACKGROUND: Clostridium difficile is capable of causing severe enterocolitis in adults. The significance of toxin-producing C. difficile in children with diarrhea is unclear and practice differs on whether to institute treatment. We aimed to characterize the microbiome in relation to the presence of C. difficile and co-infection with other pathogens and to describe host response to infection. METHODS: Participants were children with acute diarrhea, 0-16 years of age, from whom stool samples had been submitted to the hospital laboratory for routine microbiology/virology. Convenience sampling was used for 50 prospective and 150 retrospective samples. No participants were treated for C. difficile. Rates of culture positivity for C. difficile, presence of toxin and PCR-ribotype were compared between age groups. Presence of other potential pathogens, comorbidities and complications were recorded. Microbiotal diversity was measured by 16S profiling. RESULTS: Nineteen of 77 (25%) children <2 years of age and 13 of 119 (11%) children >2 years of age were C. difficile positive, of whom 10 (53%) and 9 (69%), respectively, carried toxigenic strains. Increased Shannon diversity was seen in children carrying C. difficile, with altered milieu. Presence of C. difficile was not associated with adverse clinical outcomes. In stools containing both Norovirus and C. difficile, there was increased relative abundance of verrucomicrobia. CONCLUSIONS: Children with diarrhea regularly carried toxigenic and non-toxigenic strains of C. difficile, demonstrating enhanced microbiotal diversity, and change in milieu, without apparent morbidity. This unexpected finding is contrary to that seen in adults with C. difficile disease.

Host genetics and geography influence microbiome composition in the sponge Ircinia campana.

Marine sponges are hosts to large, diverse communities of microorganisms. These microbiomes are distinct among sponge species and from seawater bacterial communities, indicating a key role of host identity in shaping its resident microbial community. However, the factors governing intraspecific microbiome variability are underexplored and may shed light on the evolutionary and ecological relationships between host and microbiome. Here, we examined the influence of genetic variation and geographic location on the composition of the Ircinia campana microbiome. We developed new microsatellite markers to genotype I. campana from two locations in the Florida Keys, USA, and characterized their microbiomes using V4 16S rRNA amplicon sequencing. We show that microbial community composition and diversity is influenced by host genotype, with more genetically similar sponges hosting more similar microbial communities. We also found that although I. campana was not genetically differentiated between sites, microbiome composition differed by location. Our results demonstrate that both host genetics and geography influence the composition of the sponge microbiome. Host genotypic influence on microbiome composition may be due to stable vertical transmission of the microbial community from parent to offspring, making microbiomes more similar by descent. Alternatively, sponge genotypic variation may reflect variation in functional traits that influence the acquisition of environmental microbes. This study reveals drivers of microbiome variation within and among locations, and shows the importance of intraspecific variability in mediating eco-evolutionary dynamics of host-associated microbiomes.

DNA extraction and amplicon production strategies deeply inf luence the outcome of gut mycobiome studies.

Microbial ecology studies are often performed through extraction of metagenomic DNA followed by amplification and sequencing of a marker. It is known that each step may bias the results. These biases have been explored for the study of bacterial communities, but rarely for fungi. Our aim was therefore to evaluate methods for the study of the gut mycobiome. We first evaluated DNA extraction methods in fungal cultures relevant to the gut. Afterwards, to assess how these methods would behave with an actual sample, stool from a donor was spiked with cells from the same cultures. We found that different extraction kits favour some species and bias against others. In terms of amplicon sequencing, we evaluated five primer sets, two for ITS2 and one for ITS1, 18S and 28S rRNA. Results showed that 18S rRNA outperformed the other markers: it was able to amplify all the species in the mock community and to discriminate among them. ITS primers showed both amplification and sequencing biases, the latter related to the variable length of the product. We identified several biases in the characterisation of the gut mycobiome and showed how crucial it is to be aware of these before drawing conclusions from the results of these studies.

Application of long read sequencing to determine expressed antigen diversity in Trypanosoma brucei infections.

Antigenic variation is employed by many pathogens to evade the host immune response, and Trypanosoma brucei has evolved a complex system to achieve this phenotype, involving sequential use of variant surface glycoprotein (VSG) genes encoded from a large repertoire of ~2,000 genes. T. brucei express multiple, sometimes closely related, VSGs in a population at any one time, and the ability to resolve and analyse this diversity has been limited. We applied long read sequencing (PacBio) to VSG amplicons generated from blood extracted from batches of mice sacrificed at time points (days 3, 6, 10 and 12) post-infection with T. brucei TREU927. The data showed that long read sequencing is reliable for resolving variant differences between VSGs, and demonstrated that there is significant expressed diversity (449 VSGs detected across 20 mice) and across the timeframe of study there was a clear semi-reproducible pattern of expressed diversity (median of 27 VSGs per sample at day 3 post infection (p.i.), 82 VSGs at day 6 p.i., 187 VSGs at day 10 p.i. and 132 VSGs by day 12 p.i.). There was also consistent detection of one VSG dominating expression across replicates at days 3 and 6, and emergence of a second dominant VSG across replicates by day 12. The innovative application of ecological diversity analysis to VSG reads enabled characterisation of hierarchical VSG expression in the dataset, and resulted in a novel method for analysing such patterns of variation. Additionally, the long read approach allowed detection of mosaic VSG expression from very few reads-the earliest in infection that such events have been detected. Therefore, our results indicate that long read analysis is a reliable tool for resolving diverse gene expression profiles, and provides novel insights into the complexity and nature of VSG expression in trypanosomes, revealing significantly higher diversity than previously shown and the ability to identify mosaic gene formation early during the infection process.

Temperature Driven Changes in Benthic Bacterial Diversity Influences Biogeochemical Cycling in Coastal Sediments.

Marine sediments are important sites for global biogeochemical cycling, mediated by macrofauna and microalgae. However, it is the microorganisms that drive these key processes. There is strong evidence that coastal benthic habitats will be affected by changing environmental variables (rising temperature, elevated CO2), and research has generally focused on the impact on macrofaunal biodiversity and ecosystem services. Despite their importance, there is less understanding of how microbial community assemblages will respond to environmental changes. In this study, a manipulative mesocosm experiment was employed, using next-generation sequencing to assess changes in microbial communities under future environmental change scenarios. Illumina sequencing generated over 11 million 16S rRNA gene sequences (using a primer set biased toward bacteria) and revealed Bacteroidetes and Proteobacteria dominated the total bacterial community of sediment samples. In this study, the sequencing coverage and depth revealed clear changes in species abundance within some phyla. Bacterial community composition was correlated with simulated environmental conditions, and species level community composition was significantly influenced by the mean temperature of the environmental regime (p = 0.002), but not by variation in CO2 or diurnal temperature variation. Species level changes with increasing mean temperature corresponded with changes in NH4 concentration, suggesting there is no functional redundancy in microbial communities for nitrogen cycling. Marine coastal biogeochemical cycling under future environmental conditions is likely to be driven by changes in nutrient availability as a direct result of microbial activity.

Comparison of the human gastric microbiota in hypochlorhydric states arising as a result of Helicobacter pylori-induced atrophic gastritis, autoimmune atrophic gastritis and proton pump inhibitor use.

Several conditions associated with reduced gastric acid secretion confer an altered risk of developing a gastric malignancy. Helicobacter pylori-induced atrophic gastritis predisposes to gastric adenocarcinoma, autoimmune atrophic gastritis is a precursor of type I gastric neuroendocrine tumours, whereas proton pump inhibitor (PPI) use does not affect stomach cancer risk. We hypothesised that each of these conditions was associated with specific alterations in the gastric microbiota and that this influenced subsequent tumour risk. 95 patients (in groups representing normal stomach, PPI treated, H. pylori gastritis, H. pylori-induced atrophic gastritis and autoimmune atrophic gastritis) were selected from a cohort of 1400. RNA extracted from gastric corpus biopsies was analysed using 16S rRNA sequencing (MiSeq). Samples from normal stomachs and patients treated with PPIs demonstrated similarly high microbial diversity. Patients with autoimmune atrophic gastritis also exhibited relatively high microbial diversity, but with samples dominated by Streptococcus. H. pylori colonisation was associated with decreased microbial diversity and reduced complexity of co-occurrence networks. H. pylori-induced atrophic gastritis resulted in lower bacterial abundances and diversity, whereas autoimmune atrophic gastritis resulted in greater bacterial abundance and equally high diversity compared to normal stomachs. Pathway analysis suggested that glucose-6-phospahte1-dehydrogenase and D-lactate dehydrogenase were over represented in H. pylori-induced atrophic gastritis versus autoimmune atrophic gastritis, and that both these groups showed increases in fumarate reductase. Autoimmune and H. pylori-induced atrophic gastritis were associated with different gastric microbial profiles. PPI treated patients showed relatively few alterations in the gastric microbiota compared to healthy subjects.

Local adaptation of a bacterium is as important as its presence in structuring a natural microbial community.

Local adaptation of a species can affect community composition, yet the importance of local adaptation compared with species presence per se is unknown. Here we determine how a compost bacterial community exposed to elevated temperature changes over 2 months as a result of the presence of a focal bacterium, Pseudomonas fluorescens SBW25, that had been pre-adapted or not to the compost for 48 days. The effect of local adaptation on community composition is as great as the effect of species presence per se, with these results robust to the presence of an additional strong selection pressure: an SBW25-specific virus. These findings suggest that evolution occurring over ecological time scales can be a key driver of the structure of natural microbial communities, particularly in situations where some species have an evolutionary head start following large perturbations, such as exposure to antibiotics or crop planting and harvesting.

The Gut Microbiota of Wild Mice.

The gut microbiota profoundly affects the biology of its host. The composition of the microbiota is dynamic and is affected by both host genetic and many environmental effects. The gut microbiota of laboratory mice has been studied extensively, which has uncovered many of the effects that the microbiota can have. This work has also shown that the environments of different research institutions can affect the mouse microbiota. There has been relatively limited study of the microbiota of wild mice, but this has shown that it typically differs from that of laboratory mice (and that maintaining wild caught mice in the laboratory can quite quickly alter the microbiota). There is also inter-individual variation in the microbiota of wild mice, with this principally explained by geographical location. In this study we have characterised the gut (both the caecum and rectum) microbiota of wild caught Mus musculus domesticus at three UK sites and have investigated how the microbiota varies depending on host location and host characteristics. We find that the microbiota of these mice are generally consistent with those described from other wild mice. The rectal and caecal microbiotas of individual mice are generally more similar to each other, than they are to the microbiota of other individuals. We found significant differences in the diversity of the microbiotas among mice from different sample sites. There were significant correlations of microbiota diversity and body weight, a measure of age, body-mass index, serum concentration of leptin, and virus, nematode and mite infection.

Long non-coding RNAs and enhancer RNAs regulate the lipopolysaccharide-induced inflammatory response in human monocytes.

Early reports indicate that long non-coding RNAs (lncRNAs) are novel regulators of biological responses. However, their role in the human innate immune response, which provides the initial defence against infection, is largely unexplored. To address this issue, here we characterize the long non-coding RNA transcriptome in primary human monocytes using RNA sequencing. We identify 76 enhancer RNAs (eRNAs), 40 canonical lncRNAs, 65 antisense lncRNAs and 35 regions of bidirectional transcription (RBT) that are differentially expressed in response to bacterial lipopolysaccharide (LPS). Crucially, we demonstrate that knockdown of nuclear-localized, NF-κB-regulated, eRNAs (IL1ß-eRNA) and RBT (IL1ß-RBT46) surrounding the IL1ß locus, attenuates LPS-induced messenger RNA transcription and release of the proinflammatory mediators, IL1ß and CXCL8. We predict that lncRNAs can be important regulators of the human innate immune response.

From a single whole exome read to notions of clinical screening: primary ciliary dyskinesia and RSPH9 p.Lys268del in the Arabian Peninsula.

Primary ciliary dyskinesia (PCD) is a genetic disorder, usually autosomal recessive, causing early respiratory disease and later subfertility. Whole exome sequencing may enable efficient analysis for locus heterogeneous disorders such as PCD. We whole-exome-sequenced one consanguineous Saudi Arabian with clinically diagnosed PCD and normal laterality, to attempt ab initio molecular diagnosis. We reviewed 13 known PCD genes and potentially autozygous regions (extended homozygosity) for homozygous exon deletions, non-dbSNP codon, splice-site base variants or small indels. Homozygous non-dbSNP changes were also reviewed exome-wide. One single molecular read representing RSPH9 p.Lys268del was observed, with no wild-type reads, and a notable deficiency of mapped reads at this location. Among all observations, RSPH9 was the strongest candidate for causality. Searching unmapped reads revealed seven more mutant reads. Direct assay for p.Lys268del (MboII digest) confirmed homozygosity in the affected individual, then confirmed homozygosity in three siblings with bronchiectasis. Our finding in southwest Saudi Arabia indicates that p.Lys268del, previously observed in two Bedouin families (Israel, UAE), is geographically widespread in the Arabian Peninsula. Analogous with cystic fibrosis CFTR p.Phe508del, screening for RSPH9 p.Lys268del (which lacks sentinel dextrocardia) in those at risk would help in early diagnosis, tailored clinical management, genetic counselling and primary prevention.

TRAM (Transcriptome Mapper): database-driven creation and analysis of transcriptome maps from multiple sources.

BACKGROUND: Several tools have been developed to perform global gene expression profile data analysis, to search for specific chromosomal regions whose features meet defined criteria as well as to study neighbouring gene expression. However, most of these tools are tailored for a specific use in a particular context (e.g. they are species-specific, or limited to a particular data format) and they typically accept only gene lists as input. RESULTS: TRAM (Transcriptome Mapper) is a new general tool that allows the simple generation and analysis of quantitative transcriptome maps, starting from any source listing gene expression values for a given gene set (e.g. expression microarrays), implemented as a relational database. It includes a parser able to assign univocal and updated gene symbols to gene identifiers from different data sources. Moreover, TRAM is able to perform intra-sample and inter-sample data normalization, including an original variant of quantile normalization (scaled quantile), useful to normalize data from platforms with highly different numbers of investigated genes. When in 'Map' mode, the software generates a quantitative representation of the transcriptome of a sample (or of a pool of samples) and identifies if segments of defined lengths are over/under-expressed compared to the desired threshold. When in 'Cluster' mode, the software searches for a set of over/under-expressed consecutive genes. Statistical significance for all results is calculated with respect to genes localized on the same chromosome or to all genome genes. Transcriptome maps, showing differential expression between two sample groups, relative to two different biological conditions, may be easily generated. We present the results of a biological model test, based on a meta-analysis comparison between a sample pool of human CD34+ hematopoietic progenitor cells and a sample pool of megakaryocytic cells. Biologically relevant chromosomal segments and gene clusters with differential expression during the differentiation toward megakaryocyte were identified. CONCLUSIONS: TRAM is designed to create, and statistically analyze, quantitative transcriptome maps, based on gene expression data from multiple sources. The release includes FileMaker Pro database management runtime application and it is freely available at http://apollo11.isto.unibo.it/software/, along with preconfigured implementations for mapping of human, mouse and zebrafish transcriptomes.