Evolution of homospermidine synthase in the convolvulaceae: a story of gene duplication, gene loss, and periods of various selection pressures.

Homospermidine synthase (HSS), the first pathway-specific enzyme of pyrrolizidine alkaloid biosynthesis, is known to have its origin in the duplication of a gene encoding deoxyhypusine synthase. To study the processes that followed this gene duplication event and gave rise to HSS, we identified sequences encoding HSS and deoxyhypusine synthase from various species of the Convolvulaceae. We show that HSS evolved only once in this lineage. This duplication event was followed by several losses of a functional gene copy attributable to gene loss or pseudogenization. Statistical analyses of sequence data suggest that, in those lineages in which the gene copy was successfully recruited as HSS, the gene duplication event was followed by phases of various selection pressures, including purifying selection, relaxed functional constraints, and possibly positive Darwinian selection. Site-specific mutagenesis experiments have confirmed that the substitution of sites predicted to be under positive Darwinian selection is sufficient to convert a deoxyhypusine synthase into a HSS. In addition, analyses of transcript levels have shown that HSS and deoxyhypusine synthase have also diverged with respect to their regulation. The impact of protein-protein interaction on the evolution of HSS is discussed with respect to current models of enzyme evolution.

A transcriptomic approach for exploring the molecular basis for dosha-balancing property-based classification of plants in Ayurveda.

In Ayurveda, a healthy body is defined by a balance among the three doshas (Vata, Pitta, Kapha) and ailments result due to imbalances among them. It prescribes specific plant parts/tissues collected in a season-specific manner for curing dosha-related imbalances but the plants prescribed for treating a particular dosha imbalance belong to taxonomically diverse families and often contain similar classes of phytomolecules, making it difficult to provide a phytochemical validation for any similarity that might exist among them. This exploratory study hypothesised that plants of the same dosha-curing group may have similarity at the transcript level. For proving/disproving the hypothesis, cDNA-AFLP and specific expression subset analysis (SESA) were carried out on the Ayurveda-defined active tissues of four representative plants each of the three dosha-balancing groups. cDNA-AFLP analyses indicated that even though the plants belonging to a particular dosha-group may widely differ at the transcript level, there is a small fraction of transcripts that is monomorphic among their active tissues. SESA (Tester-active tissue cDNA; Driver-cDNA from other major tissue[s]) generated 803 subtractive ESTs from the twelve plants that yielded 150 unigenes upon assembly (of ESTs from each plant separately). Cross-plant EST assembly for plants in the same dosha group also corroborated the results. Although a distinct pattern of transcripts was not observed across all the plants in a particular dosha group, some commonalities were obtained that need further characterization towards searching for the hitherto elusive similarity among plants of the same group.

Agave tequilana MADS genes show novel expression patterns in meristems, developing bulbils and floral organs.

Agave tequilana is a monocarpic perennial species that flowers after 5-8 years of vegetative growth signaling the end of the plant's life cycle. When fertilization is unsuccessful, vegetative bulbils are induced on the umbels of the inflorescence near the bracteoles from newly formed meristems. Although the regulation of inflorescence and flower development has been described in detail for monocarpic annuals and polycarpic species, little is known at the molecular level for these processes in monocarpic perennials, and few studies have been carried out on bulbils. Histological samples revealed the early induction of umbel meristems soon after the initiation of the vegetative to inflorescence transition in A. tequilana. To identify candidate genes involved in the regulation of floral induction, a search for MADS-box transcription factor ESTs was conducted using an A. tequilana transcriptome database. Seven different MIKC MADS genes classified into 6 different types were identified based on previously characterized A. thaliana and O. sativa MADS genes and sequences from non-grass monocotyledons. Quantitative real-time PCR analysis of the seven candidate MADS genes in vegetative, inflorescence, bulbil and floral tissues uncovered novel patterns of expression for some of the genes in comparison with orthologous genes characterized in other species. In situ hybridization studies using two different genes showed expression in specific tissues of vegetative meristems and floral buds. Distinct MADS gene regulatory patterns in A. tequilana may be related to the specific reproductive strategies employed by this species.

[Analysis on correlation between 3-hydroxy-3-methylglutary-coenzyme A reductase gene polymorphism of Glycyrrhiza uralensis and content of glycyrrhizic acid].

OBJECTIVE: To reveal the 3-hydroxy-3-methylglutary-coenzyme A reductase (HMGR) gene polymorphism of Glycyrrhiza uralensis, and the correlation between HMGR gene polymorphism and the content of glycyrrhizic acid. METHOD: Liquorice plants containing different content of glycyrrhizic acid were used as materials. RT-PCR was used to amplify their HMGR gene sequences, which were connected with vector pMD19-T for clone sequencing. Multiple alignments were performed to analyse HMGR gene polymorphism of G. uralensis. Then the correlation between HMGR gene polymorphism and the content of glycyrrhizic acid was revealed. RESULT: HMGR gene sequences polymorphism included codon mutation, base substitution mutation, copy number polymorphism and allele heterozygosity. There were 4 types of mutations in HMGR gene coding amino acid sequences, namely -HSL, -HSV, GALLV, GALSV. Among them, -HSV type was common in liquorice plants, -HSL type only existed in liquorice plants with low content of glycyrrhizic acid, and GALSV type only existed in liquorice plants with high content of glycyrrhizic acid. CONCLUSION: HMGR gene sequences of G. uralensis are highly polymorphic and related to the content of glycyrrhizic acid.

Identification of salt-induced genes from Salicornia brachiata, an extreme halophyte through expressed sequence tags analysis.

Salinity severely affects plant growth and development causing crop loss worldwide. We have isolated a large number of salt-induced genes as well as unknown and hypothetical genes from Salicornia brachiata Roxb. (Amaranthaceae). This is the first description of identification of genes in response to salinity stress in this extreme halophyte plant. Salicornia accumulates salt in its pith and survives even at 2 M NaCl under field conditions. For isolating salt responsive genes, cDNA subtractive hybridization was performed between control and 500 mM NaCl treated plants. Out of the 1200 recombinant clones, 930 sequences were submitted to the NCBI database (GenBank accession: EB484528 to EB485289 and EC906125 to EC906292). 789 ESTs showed matching with different genes in NCBI database. 4.8% ESTs belonged to stress-tolerant gene category and approximately 29% ESTs showed no homology with known functional gene sequences, thus classified as unknown or hypothetical. The detection of a large number of ESTs with unknown putative function in this species makes it an interesting contribution. The 90 unknown and hypothetical genes were selected to study their differential regulation by reverse Northern analysis for identifying their role in salinity tolerance. Interestingly, both up and down regulation at 500 mM NaCl were observed (21 and 10 genes, respectively). Northern analysis of two important salt tolerant genes, ASR1 (Abscisic acid stress ripening gene) and plasma membrane H+ATPase, showed the basal level of transcripts in control condition and an increase with NaCl treatment. ASR1 gene is made full length using 5' RACE and its potential role in imparting salt tolerance is being studied.

[Neprilysin--structure of the gene and protein product and the localization of expression]. / Neprylizyna--budowa genu i produktu bialkowego oraz lokalizacja ekspresji.

Neprilysin (NEP, CD10, CALLA - common acute lymphoblastic leukaemia antigen, neutral endopeptidase, enkephalinase) is a zinc-dependent metallopeptidase, which is involved in the metabolism of a number of regulatory peptides and plays an important role in turning off peptide signalling at the cell surface. NEP gene is located on chromosome 3q 25.1-q25.2 and is composed of 24 exons. Four types of NEP cDNAs have been identified resulting from alternative splicing of exons 1, 1bis, 2a or 2b to the common exon 3. Neprilysin is expressed in normal and malignant hematopoietic cells and in epithelial cells of many organs. In kidneys, it is expressed in podocytes, renal proximal tubular epithelium and in smooth muscles of the vessels.

Functional annotation of 19,841 Populus nigra full-length enriched cDNA clones.

BACKGROUND: Populus is one of favorable model plants because of its small genome. Structural genomics of Populus has reached a breakpoint as nucleotides of the entire genome have been determined. Reaching the post genome era, functional genomics of Populus is getting more important for well-comprehended plant science. Development of bioresorce serving functional genomics is making rapid progress. Huge efforts have achieved deposits of expressed sequence tags (ESTs) in various plant species consequently accelerating functional analysis of genes. ESTs from full-length cDNA clones are especially powerful for accurate molecular annotation. We promoted collection and annotation of the ESTs from Populus full-length enriched cDNA clones as part of functional genomics of tree species. RESULTS: We have been collecting the full-length enriched cDNA of the female poplar (Populus nigra var. italica) for years. By sequencing P. nigra full-length (PnFL) cDNA libraries, we generated about 116,000 5'-end or 3'-end ESTs corresponding to 19,841 nonredundant PnFL clones. Population of PnFL cDNA clones represents 44% of the predicted genes in the Populus genome. CONCLUSION: Our resource of P. nigra full-length enriched clones is expected to provide valuable tools to gain further insight into genome annotation and functional genomics in Populus.

Application of SSH and a macroarray to investigate altered gene expression in Mytilus edulis in response to exposure to benzo[a]pyrene.

The lack of genomic resources for aquatic invertebrates restricts their use as sentinel species in coastal environments. It is known that where genomic data are not available, suppression subtractive hybridisation (SSH) can generate cDNA libraries representative of pollutant-responsive gene transcription in aquatic vertebrates. To assess whether the approach was equally suited to aquatic invertebrates, altered gene expression in digestive gland of the mussel, Mytilus edulis, in response to exposure to benzo[a]pyrene (BaP) (1 mg/l) was investigated with SSH and a nylon macroarray. Screening of the subtracted libraries showed 112/250 up-regulated and 25/55 down-regulated clones were positive for differential expression and characterisation of these identified 87 with unique sequence suitable for array on a nylon membrane. The transcripts isolated were from a diverse range of genes involved in general stress, oxidative stress, cell adhesion, transcriptional and translational regulation, transport mechanisms, energy metabolism, cell metabolism, lipid metabolism, protein turnover and activation, lysosomal activity and 22 cryptic clones. Subsequent use of the clones in macroarray format to analyse expression of BaP-responsive genes (0 vs 4 day exposed) showed 0-100-fold increased levels of the forward-subtracted probes and between 0 and 0.1-fold down-regulation of the reverse-subtracted probes. Only 15% of the clones showed less than 2-fold change in expression. The gene ontology of the transcripts isolated demonstrates that BaP elicits a multitude of responses with a major feature being disruption of cellular redox status. The results indicate that the use of SSH and a macroarray is a robust method to discover novel pollutant-responsive genes in aquatic invertebrates.

Evaluation of normalization methods for cDNA microarray data by k-NN classification.

BACKGROUND: Non-biological factors give rise to unwanted variations in cDNA microarray data. There are many normalization methods designed to remove such variations. However, to date there have been few published systematic evaluations of these techniques for removing variations arising from dye biases in the context of downstream, higher-order analytical tasks such as classification. RESULTS: Ten location normalization methods that adjust spatial- and/or intensity-dependent dye biases, and three scale methods that adjust scale differences were applied, individually and in combination, to five distinct, published, cancer biology-related cDNA microarray data sets. Leave-one-out cross-validation (LOOCV) classification error was employed as the quantitative end-point for assessing the effectiveness of a normalization method. In particular, a known classifier, k-nearest neighbor (k-NN), was estimated from data normalized using a given technique, and the LOOCV error rate of the ensuing model was computed. We found that k-NN classifiers are sensitive to dye biases in the data. Using NONRM and GMEDIAN as baseline methods, our results show that single-bias-removal techniques which remove either spatial-dependent dye bias (referred later as spatial effect) or intensity-dependent dye bias (referred later as intensity effect) moderately reduce LOOCV classification errors; whereas double-bias-removal techniques which remove both spatial- and intensity effect reduce LOOCV classification errors even further. Of the 41 different strategies examined, three two-step processes, IGLOESS-SLFILTERW7, ISTSPLINE-SLLOESS and IGLOESS-SLLOESS, all of which removed intensity effect globally and spatial effect locally, appear to reduce LOOCV classification errors most consistently and effectively across all data sets. We also found that the investigated scale normalization methods do not reduce LOOCV classification error. CONCLUSION: Using LOOCV error of k-NNs as the evaluation criterion, three double-bias-removal normalization strategies, IGLOESS-SLFILTERW7, ISTSPLINE-SLLOESS and IGLOESS-SLLOESS, outperform other strategies for removing spatial effect, intensity effect and scale differences from cDNA microarray data. The apparent sensitivity of k-NN LOOCV classification error to dye biases suggests that this criterion provides an informative measure for evaluating normalization methods. All the computational tools used in this study were implemented using the R language for statistical computing and graphics.

Classification of rice allergenic protein cDNAs belonging to the alpha-amylase/trypsin inhibitor gene family.

Seven cDNA clones encoding rice allergenic proteins were newly isolated. Comparison of the sequences of ten cDNA clones, including the previously isolated three clones results in their classification into four subfamilies. Homologies in the nucleotide sequences among and within subfamilies are 70-85% and above 95%, respectively. A sequence of twenty five amino-acid residues at the C-terminal proximal region is highly conserved among all clones and resembles that of plant lipid transfer proteins.

Rat homologues of yeast sec7p.

Mutations in the Saccharomyces cerevisiae sec7 locus lead to a pleiotropic secretory phenotype that is characterized by an accumulation of Golgi cisternae and a loss of secretory granules. This indicates that the corresponding gene product sec7p is involved in the budding of secretory granules from the Golgi apparatus. Here we report the primary structure of three rat homologues of sec7p, called msec7-1, -2, and -3. The mRNAs of these genes are expressed in all tissues tested. All msec7s share the same domain structure in which an N-terminal coiled-coil domain is followed by a sec7-homology domain and a pleckstrin-homology domain. On the protein level, msec7s are present in all rat tissues tested, with highest protein levels in brain and adrenal. In the adult rat brain, they are present in soluble and membrane-associated pools.

Negative subtraction hybridization: an efficient method to isolate large numbers of condition-specific cDNAs.

BACKGROUND: The construction of cDNA libraries is a useful tool to understand gene expression in organisms under different conditions, but random sequencing of unbiased cDNA collections is laborious and can give rise to redundant EST collections. We aimed to isolate cDNAs of messages induced by switching Aspergillus nidulans from growth on glucose to growth on selected polysaccharides. Approximately 4,700 contigs from 12,320 ESTs were already available from a cDNA library representing transcripts isolated from glucose-grown A. nidulans during asexual development. Our goals were to expand the cDNA collection without repeated sequencing of previously identified ESTs and to find as many transcripts as possible that are specifically induced in complex polysaccharide metabolism. RESULTS: We have devised a Negative Subtraction Hybridization (NSH) method and tested it in A. nidulans. NSH entails screening a plasmid library made from cDNAs prepared from cells grown under a selected physiological condition with labeled cDNA probes prepared from another physiological condition. Plasmids with inserts that failed to hybridize to cDNA probes through two rounds of screening (i.e. negatives) indicate that they are transcripts present at low concentration in the labeled probe pool. Thus, these transcripts will be predominantly condition-specific, along with some rare transcripts. In a screen for transcripts induced by switching the carbon source from glucose to 12 selected polysaccharides, 3,532 negatives were isolated from approximately 100,000 surveyed colonies using this method. Negative clones were end-sequenced and assembled into 2,039 contigs, of which 1,722 were not present in the previously characterized glucose-grown cDNA library. Single-channel microarray hybridization experiments confirmed that the majority of the negatives represented genes that were differentially induced by a switch from growth in glucose to one or more of the polysaccharides. CONCLUSIONS: The Negative Subtraction Hybridization method described here has several practical benefits. This method can be used to screen any existing cDNA library, including full-length and pooled libraries, and does not rely on PCR or sequence information. In addition, NSH is a cost-effective method for the isolation of novel, full-length cDNAs for differentially expressed transcripts or enrichment of rare transcripts.

Molecular cloning of newt sex pheromone precursor cDNAs: evidence for the existence of species-specific forms of pheromones.

Cloning of cDNA encoding a decapeptide pheromone (sodefrin) that attracts conspecific female newts was attempted. A cDNA clone encoding a protein consisting of 189 amino acid residues including a sodefrin sequence was isolated from a Cynops pyrrhogaster abdominal gland cDNA library. Likewise, a cDNA clone encoding a molecule comparable to the sodefrin precursor was obtained from a Cynops ensicauda abdominal gland cDNA library. This clone encoded a precursor protein of 192 amino acid residues, including a sodefrin-like peptide sequence with substitutions of two amino acid residues. This is the first report of a peptide pheromone precursor in vertebrates.

Mammalian HSP40/DNAJ homologs: cloning of novel cDNAs and a proposal for their classification and nomenclature.

We have cloned 10 novel full-length cDNAs of mouse and human HSP40/DNAJ homologs using expressed sequence tag (EST) clones found in the DDBJ/GenBank/EMBL DNA database. In this report, we tentatively designated them mHsp40, mDj3, mDj4, mDj5, mDj6, mDj7, mDj8, hDj9, mDj10, and mDj11. Based on the identity of the deduced amino acid sequences, mHsp40, mDj3, and mDj11 are orthologs of human Hsp40, rat Rdj2, and human Tpr2, respectively. We determined that mDj4 is identical with the recently isolated mouse Mrj (mammalian relative of DnaJ). PSORT analysis (a program that predicts the subcellular localization site of a given protein from its amino acid sequences) revealed that hDj9 has an N-terminal signal peptide; hence, its localization might be extracellular, suggesting that there may be a partner Hsp70 protein that acts together with the hDj9 outside of the cell. The same analysis indicated that mDj7 and mDj10 may have transmembrane domains. In order to simplify the complicated and confusing nomenclature of recently identified mammalian HSP40/DNAJ homologs, we propose here some new rules for their nomenclature. This proposed nomenclature includes the name of species with 2 lowercase letters such as hs (Homo sapiens), mm (Mus musculus) and rn (Rattus norvegicus); Dj standing for DnaJ; the name of types with A, B, and C, which were previously classified as type I, II, and III according to the domain structure of the homologs; and finally Arabic numerals according to the chronological order of registration of the sequence data into the database.

An insight into the salivary transcriptome and proteome of the adult female mosquito Culex pipiens quinquefasciatus.

To obtain an insight into the salivary transcriptome and proteome (sialome) of the adult female mosquito Culex quinquefasciatus, a cDNA library was randomly sequenced, and aminoterminal information for selected proteins and peptides was obtained. cDNA sequence clusters coding for secreted proteins were further analyzed. The transcriptome revealed messages coding for several proteins of known families previously reported in the salivary glands of other blood-feeding insects as well as immune-related products such as C-type lectin, gambicin, and members of the prophenol oxidase cascade. Additionally, several transcripts coding for low-complexity proteins were found, some clearly coding for mucins. Many novel transcripts were found, including a novel endonuclease previously described in crabs and shrimps but not in insects; a hyaluronidase, not described before in mosquito salivary glands but found in venom glands and in salivary glands of sand flies and black flies; several cysteine-rich peptides with possible anticlotting function, including one similar to a previously described nematode family of anti-proteases; and a completely novel family of cysteine- and tryptophane-rich proteins (CWRC family) for which 12 full-length sequences are described. Also described are 14 additional novel proteins and peptides whose function and/or family affiliation are unknown. In total, 54 transcripts coding for full-length proteins are described. That several of these are translated into proteins was confirmed by finding the corresponding aminoterminal sequences in the SDS-PAGE/Edman degradation experiments. Electronic versions of all tables and sequences can be found at http://www.ncbi.nlm.nih.gov/projects/Mosquito/C_quinquefasciatus_sialome.

Age-associated changes in hypothalamic and pituitary neuroendocrine gene expression in the rat.

A cDNA membrane array displaying 1183 probes was used to detect hypothalamic and pituitary changes in gene expression accompanying ageing and age-associated pituitary macroadenomas. Four groups of male Sprague-Dawley rats (3-, 15-, 24-month-old and 24-month-old with prolactinoma) were compared in two independent hybridizations. cDNA array data were confirmed and completed by comparative reverse transcriptase-polymerase chain reaction on selected genes. The expression of 454 and 116 mRNAs was detected in hypothalamus and pituitary, respectively. Growth hormone (GH) mRNA alone represented 85% of total gene expression in the gland of young rats, and other pituitary hormone transcripts 2.8%, while melanin-concentrating hormone (MCH) mRNA, the most expressed neuropeptide transcript involved in neuroendocrine regulation, accounted for only 0.8% of total hypothalamic transcripts. The proportion of genes modified in the hypothalamus and pituitary was rather modest: 1.5% and 5.2%, respectively, for ageing per se, and 1.1% and 5.2% for age-associated macroprolactinomas. Among pituitary specific RNAs, GH mRNA expression was notably decreased with age. At the hypothalamic level, expression of genes directly involved in GH regulation, such as somatostatin and growth hormone-releasing hormone, was not altered, while neuropeptide transcripts involved in feeding behaviour [orexin/hypocretin, MCH, pro-opiomelanocortin (POMC), cocaine- and amphetamine-regulated transcript (CART)] were significantly altered. In addition, a few ubiquitous transcripts (hnRNP-K, PFKm, CCND 2, calponin and set) were differently affected in both tissues. Modifications in hypothalamic orexigenic (orexin, MCH) and anorexigenic (POMC, CART) gene expression are in keeping with an age-associated decrease in energy consumption but a higher one in the presence of macroprolactinomas.

Isolation of differential genes in suspension cultures of Taxus cuspidata induced by additional taxol.

Addition of taxol into suspension cultures of Taxus cuspidata induced cell apoptosis, which was confirmed by gel electrophoresis of the DNA ladders indicating the progressive delineation of fragmented nuclear DNA (nDNA) into distinct bodies. The additional taxol not only changed the microtubule assembly of cells, but also affected the gene expression. Fourteen cDNA fragments, named as TIGT9-22, were isolated after addition of taxol and their GenBank accession numbers were given as BF704560-BF704573, respectively. Among them, TIGT13 and TIGT21 were apparently homogeneous with apbE and carbamoylphosphate synthetase, respectively. Other cDNA fragments showed no significant analogy with the known sequences in GenBank.

Consistency checks for characterizing protein forms.

Proteomics enforces the reverse chronological order on the gene to protein dogma and imposes amino acid sequences as a starting point of an investigation relative to function. By this approach, proteomics data can confirm the presence of multiple forms of a protein. Notwithstanding variations attributed specific individual features of organisms and tissues, from two to over ten protein forms can be identified in a given sample. The present work describes some guidelines for tracking the origin of alternative protein forms and attempts to tag the details of sequence data in the literature. Working via these guidelines we have uncovered a third alternative form of the Pim subfamily of oncogenes. The term form is here combined with the qualification alternative to describe any product of a given gene including closely related paralogs. This paper also emphasizes the need for consistency checks in annotation processes, such as gene clustering, to avoid losing important details describing protein alternative forms. By identifying alternative protein forms, we illustrate the fact that rationalizing of protein function via the identification of protein-protein interactions should in reality be that of identifying (alternative) form-form interactions.

Cataloging altered gene expression during rat hippocampal long-term potentiation by means of differential display.

We have employed mRNA differential display (DD) to generate a catalog of cDNAs whose expression in the hippocampus was regulated during long-term potentiation (LTP) in dentate gyrus of anesthetized rats. DD with 459 combinations of primer pairs revealed that 80 out of approximately 70000 bands displayed showed a reproducible change in their expression level. These cDNAs were categorized into seven groups according to changes in their temporal expression pattern. Some of these cDNAs were induced rapidly, but transiently, after the LTP induction, some induced rapidly and persistently, some induced slowly, and some down-regulated following LTP. This suggests that a complex molecular hierarchy underlies the maintenance of hippocampal LTP.

Intestinal microbiota composition of interleukin-10 deficient C57BL/6J mice and susceptibility to Helicobacter hepaticus-induced colitis.

The mouse pathobiont Helicobacter hepaticus can induce typhlocolitis in interleukin-10-deficient mice, and H. hepaticus infection of immunodeficient mice is widely used as a model to study the role of pathogens and commensal bacteria in the pathogenesis of inflammatory bowel disease. C57BL/6J Il10(-/-) mice kept under specific pathogen-free conditions in two different facilities (MHH and MIT), displayed strong differences with respect to their susceptibilities to H. hepaticus-induced intestinal pathology. Mice at MIT developed robust typhlocolitis after infection with H. hepaticus, while mice at MHH developed no significant pathology after infection with the same H. hepaticus strain. We hypothesized that the intestinal microbiota might be responsible for these differences and therefore performed high resolution analysis of the intestinal microbiota composition in uninfected mice from the two facilities by deep sequencing of partial 16S rRNA amplicons. The microbiota composition differed markedly between mice from both facilities. Significant differences were also detected between two groups of MHH mice born in different years. Of the 119 operational taxonomic units (OTUs) that occurred in at least half the cecum or colon samples of at least one mouse group, 24 were only found in MIT mice, and another 13 OTUs could only be found in MHH samples. While most of the MHH-specific OTUs could only be identified to class or family level, the MIT-specific set contained OTUs identified to genus or species level, including the opportunistic pathogen, Bilophila wadsworthia. The susceptibility to H. hepaticus-induced colitis differed considerably between Il10(-/-) mice originating from the two institutions. This was associated with significant differences in microbiota composition, highlighting the importance of characterizing the intestinal microbiome when studying murine models of IBD.