GUT MICROBIOME DIVERSITY OF THREE RHINOCEROS SPECIES IN EUROPEAN ZOOS.

The wild rhinoceros populations have declined drastically in the past decades because the rhinoceros are heavily hunted for their horns. Zoological institutions aim to conserve rhinoceros populations in captivity, but one of the challenges of ex situ conservation is to provide food sources that resemble those available in the wild. Considering that the mammalian gut microbiota is a pivotal player in their host's health, the gut microbiota of rhinoceros may also play a role in the bioavailability of nutrients. Therefore, this study aims to characterize the fecal microbiome composition of grazing white rhinoceros (WR; Ceratotherium simum) and greater one-horned rhinoceros (GOHR; Rhinoceros unicornis) as well as the browsing black rhinoceros (BR; Diceros bicornis) kept in European zoos. Over the course of 1 yr, 166 fecal samples in total were collected from 9 BR (n = 39), 10 GOHR (n = 56), and 14 WR (n = 71) from 23 zoological institutions. The bacterial composition in the samples was determined using 16S rRNA gene Illumina sequencing. The fecal microbiomes of rhinoceros clustered by species, with BR clustering more distantly from GOHR and WR. Furthermore, the data report clustering of rhinoceros microbiota according to individual rhinoceros and institutional origin, showing that zoological institutions play a significant role in shaping the gut microbiome of rhinoceros species. In addition, BR exhibit a relatively higher microbial diversity than GOHR and WR. BR seem more susceptible to microbial gut changes and appear to have a more diverse microbiome composition among individuals than GOHR and WR. These data expand on the role of gut microbes and can provide baseline data for continued efforts in rhinoceros conservation and health status.

Absence/presence calling in microarray-based CGH experiments with non-model organisms.

Structural variations in genomes are commonly studied by (micro)array-based comparative genomic hybridization. The data analysis methods to infer copy number variation in model organisms (human, mouse) are established. In principle, the procedures are based on signal ratios between test and reference samples and the order of the probe targets in the genome. These procedures are less applicable to experiments with non-model organisms, which frequently comprise non-sequenced genomes with an unknown order of probe targets. We therefore present an additional analysis approach, which does not depend on the structural information of a reference genome, and quantifies the presence or absence of a probe target in an unknown genome. The principle is that intensity values of target probes are compared with the intensities of negative-control probes and positive-control probes from a control hybridization, to determine if a probe target is absent or present. In a test, analyzing the genome content of a known bacterial strain: Staphylococcus aureus MRSA252, this approach proved to be successful, demonstrated by receiver operating characteristic area under the curve values larger than 0.9995. We show its usability in various applications, such as comparing genome content and validating next-generation sequencing reads from eukaryotic non-model organisms.

Staphylococcus aureus ST398 gene expression profiling during ex vivo colonization of porcine nasal epithelium.

BACKGROUND: Staphylococcus aureus is a common human and animal opportunistic pathogen. In humans nasal carriage of S. aureus is a risk factor for various infections. Methicillin-resistant S. aureus ST398 is highly prevalent in pigs in Europe and North America. The mechanism of successful pig colonization by MRSA ST398 is poorly understood. Previously, we developed a nasal colonization model of porcine nasal mucosa explants to identify molecular traits involved in nasal MRSA colonization of pigs. RESULTS: We report the analysis of changes in the transcription of MRSA ST398 strain S0462 during colonization on the explant epithelium. Major regulated genes were encoding metabolic processes and regulation of these genes may represent metabolic adaptation to nasal mucosa explants. Colonization was not accompanied by significant changes in transcripts of the main virulence associated genes or known human colonization factors. Here, we documented regulation of two genes which have potential influence on S. aureus colonization; cysteine extracellular proteinase (scpA) and von Willebrand factor-binding protein (vWbp, encoded on SaPIbov5). Colonization with isogenic-deletion strains (Δvwbp and ΔscpA) did not alter the ex vivo nasal S. aureus colonization compared to wild type. CONCLUSIONS: Our results suggest that nasal colonization with MRSA ST398 is a complex event that is accompanied with changes in bacterial gene expression regulation and metabolic adaptation.

ONT read assembly of the black rhino genome.

OBJECTIVES: The black rhinoceros (Diceros bicornis) is an endangered mammal for which a captive breeding program is part of the conservation effort. Black rhinos in zoo's often suffer from chronic infections and heamochromatosis. Furthermore, breeding is hampered by low male fertility. To aid a research project studying these topics, we sequenced and assembled the genome of a captive male black rhino using ONT sequencing data only. DATA DESCRIPTION: This work produced over 100 Gb whole genome sequencing reads from whole blood. These were assembled into a 2.47 Gb draft genome consisting of 834 contigs with an N50 of 29.53 Mb. The genome annotation was lifted over from an available genome annotation for black rhino, which resulted in the retrieval of over 99% of gene features. This new genome assembly will be a valuable resource in for conservation genetic research in this species.

Interorgan coordination of the murine adaptive response to fasting.

Starvation elicits a complex adaptive response in an organism. No information on transcriptional regulation of metabolic adaptations is available. We, therefore, studied the gene expression profiles of brain, small intestine, kidney, liver, and skeletal muscle in mice that were subjected to 0-72 h of fasting. Functional-category enrichment, text mining, and network analyses were employed to scrutinize the overall adaptation, aiming to identify responsive pathways, processes, and networks, and their regulation. The observed transcriptomics response did not follow the accepted "carbohydrate-lipid-protein" succession of expenditure of energy substrates. Instead, these processes were activated simultaneously in different organs during the entire period. The most prominent changes occurred in lipid and steroid metabolism, especially in the liver and kidney. They were accompanied by suppression of the immune response and cell turnover, particularly in the small intestine, and by increased proteolysis in the muscle. The brain was extremely well protected from the sequels of starvation. 60% of the identified overconnected transcription factors were organ-specific, 6% were common for 4 organs, with nuclear receptors as protagonists, accounting for almost 40% of all transcriptional regulators during fasting. The common transcription factors were PPARα, HNF4α, GCRα, AR (androgen receptor), SREBP1 and -2, FOXOs, EGR1, c-JUN, c-MYC, SP1, YY1, and ETS1. Our data strongly suggest that the control of metabolism in four metabolically active organs is exerted by transcription factors that are activated by nutrient signals and serves, at least partly, to prevent irreversible brain damage.

Low-carbon acclimation in carboxysome-less and photorespiratory mutants of the cyanobacterium Synechocystis sp. strain PCC 6803.

Using metabolic and transcriptomic phenotyping, we studied acclimation of cyanobacteria to low inorganic carbon (LC) conditions and the requirements for coordinated alteration of metabolism and gene expression. To analyse possible metabolic signals for LC sensing and compensating reactions, the carboxysome-less mutant ΔccmM and the photorespiratory mutant ΔglcD1/D2 were compared with wild-type (WT) Synechocystis. Metabolic phenotyping revealed accumulation of 2-phosphoglycolate (2PG) in ΔccmM and of glycolate in ΔglcD1/D2 in LC- but also in high inorganic carbon (HC)-grown mutant cells. The accumulation of photorespiratory metabolites provided evidence for the oxygenase activity of RubisCO at HC. The global gene expression patterns of HC-grown ΔccmM and ΔglcD1/D2 showed differential expression of many genes involved in photosynthesis, high-light stress and N assimilation. In contrast, the transcripts of LC-specific genes, such as those for inorganic carbon transporters and components of the carbon-concentrating mechanism (CCM), remained unchanged in HC cells. After a shift to LC, ΔglcD1/D2 and WT cells displayed induction of many of the LC-inducible genes, whereas ΔccmM lacked similar changes in expression. From the coincidence of the presence of 2PG in ΔccmM without CCM induction and of glycolate in ΔglcD1/D2 with CCM induction, we regard a direct role for 2PG as a metabolic signal for the induction of CCM during LC acclimation as less likely. Instead, our data suggest a potential role for glycolate as a signal molecule for enhanced expression of CCM genes.

Operon structure of Staphylococcus aureus.

In bacteria, gene regulation is one of the fundamental characteristics of survival, colonization and pathogenesis. Operons play a key role in regulating expression of diverse genes involved in metabolism and virulence. However, operon structures in pathogenic bacteria have been determined only by in silico approaches that are dependent on factors such as intergenic distances and terminator/promoter sequences. Knowledge of operon structures is crucial to fully understand the pathophysiology of infections. Presently, transcriptome data obtained from growth curves in a defined medium were used to predict operons in Staphylococcus aureus. This unbiased approach and the use of five highly reproducible biological replicates resulted in 93.5% significantly regulated genes. These data, combined with Pearson's correlation coefficients of the transcriptional profiles, enabled us to accurately compile 93% of the genome in operon structures. A total of 1640 genes of different functional classes were identified in operons. Interestingly, we found several operons containing virulence genes and showed synergistic effects for two complement convertase inhibitors transcribed in one operon. This is the first experimental approach to fully identify operon structures in S. aureus. It forms the basis for further in vitro regulation studies that will profoundly advance the understanding of bacterial pathophysiology in vivo.

Finding maximal transcriptome differences between reprotoxic and non-reprotoxic phthalate responses in rat testis.

The chemical legislation of the EU, Registration, Evaluation, and Authorization of Chemicals (REACH), stipulates that about 30 000 chemical substances are to be assessed on their possible risks. Toxicological evaluation of these compounds will at least partly be based on animal testing. In particular, the assessment of reproductive toxicity is a very complicated, time-consuming and animal-demanding process. Introducing microarray-based technologies can potentially refine in vivo toxicity testing. If compounds of a distinct chemical class induce reproducible gene-expression responses with a recognizable overlap, these gene-expression signatures may indicate intrinsic features of certain compounds, including specific toxicity. In the present study, we have set out the first steps towards this approach for the reproductive toxicity of phthalates. Male rats were treated with a single dose of either reprotoxic or non-reprotoxic phthalates, and were analyzed 24 h afterwards. Subsequently, histopathological and gene-expression profiling analyses were performed. Despite ambiguous histopathological observations, we were able to identify genes with differential expression profiles between the reprotoxic phthalates and the non-reprotoxic counterparts. This shows that differences in gene-expression profiles, indicative of the type of exposure, may be detected earlier, or at lower doses, than classical pathological endpoints. These findings are promising for 'early warning' biomarker analyses and for using toxicogenomics in a category approach. Ultimately, this could lead to a more cost-effective approach for prioritizing the toxicity testing of large numbers of chemicals in a short period of time in hazard assessment of chemicals, which is one of the objectives of the REACH chemical legislation.

Combination of IL-21 and IL-15 enhances tumour-specific cytotoxicity and cytokine production of TCR-transduced primary T cells.

IL-21, and to a lesser extent IL-15, inhibits differentiation of antigen-primed CD8 T cells and promotes their homeostasis and anti-tumour activity. Here, we investigated molecular mechanisms behind tumour-specific responses of primary murine T lymphocytes engineered to express a TCR directed against human gp100/HLA-A2 following short-term exposure to IL-15 and/or IL-21. We demonstrated that IL-15 + IL-21, and to a lesser extent IL-21, enhanced antigen-specific T-cell cytotoxicity, which was related to enhanced expression of granzymes A and B, and perforin 1. Furthermore, IL-15 + IL-21 synergistically enhanced release levels and kinetics of T-cell IFNgamma and IL-2, but not IL-10. Enhanced secretion of IFNgamma was accompanied by increased gene expression and cytosolic protein content, and was restricted to effector memory T cells. To summarize, we show that IL-15 + IL-21 improves antigen-specific responses of TCR-transduced effector T cells at multiple levels, which provides a rationale to treat T cells with a combination of these cytokines prior to their use in adoptive TCR gene therapy.

The core genome of the anaerobic oral pathogenic bacterium Porphyromonas gingivalis.

BACKGROUND: The Gram negative anaerobic bacterium Porphyromonas gingivalis has long been recognized as a causative agent of periodontitis. Periodontitis is a chronic infectious disease of the tooth supporting tissues eventually leading to tooth-loss. Capsular polysaccharide (CPS) of P. gingivalis has been shown to be an important virulence determinant. Seven capsular serotypes have been described. Here, we used micro-array based comparative genomic hybridization analysis (CGH) to analyze a representative of each of the capsular serotypes and a non-encapsulated strain against the highly virulent and sequenced W83 strain. We defined absent calls using Arabidopsis thaliana negative control probes, with the aim to distinguish between aberrations due to mutations and gene gain/loss. RESULTS: Our analyses allowed us to call aberrant genes, absent genes and divergent regions in each of the test strains. A conserved core P. gingivalis genome was described, which consists of 80% of the analyzed genes from the sequenced W83 strain. The percentage of aberrant genes between the test strains and control strain W83 was 8.2% to 13.7%. Among the aberrant genes many CPS biosynthesis genes were found. Most other virulence related genes could be found in the conserved core genome. Comparing highly virulent strains with less virulent strains indicates that hmuS, a putative CobN/Mg chelatase involved in heme uptake, may be a more relevant virulence determinant than previously expected. Furthermore, the description of the 39 W83-specific genes could give more insight in why this strain is more virulent than others. CONCLUSION: Analyses of the genetic content of the P. gingivalis capsular serotypes allowed the description of a P. gingivalis core genome. The high resolution data from three types of analysis of triplicate hybridization experiments may explain the higher divergence between P. gingivalis strains than previously recognized.

Photorespiratory 2-phosphoglycolate metabolism and photoreduction of O2 cooperate in high-light acclimation of Synechocystis sp. strain PCC 6803.

In cyanobacteria, photorespiratory 2-phosphoglycolate (2PG) metabolism is mediated by three different routes, including one route involving the glycine decarboxylase complex (Gcv). It has been suggested that, in addition to conversion of 2PG into non-toxic intermediates, this pathway is important for acclimation to high-light. The photoreduction of O(2) (Mehler reaction), which is mediated by two flavoproteins Flv1 and Flv3 in cyanobacteria, dissipates excess reductants under high-light by the four electron-reduction of oxygen to water. Single and double mutants defective in these processes were constructed to investigate the relation between photorespiratory 2PG-metabolism and the photoreduction of O(2) in the cyanobacterium Synechocystis sp. PCC 6803. The single mutants Deltaflv1, Deltaflv3, and DeltagcvT, as well as the double mutant Deltaflv1/DeltagcvT, were completely segregated but not the double mutant Deltaflv3/DeltagcvT, suggesting that the T-protein subunit of the Gcv (GcvT) and Flv3 proteins cooperate in an essential process. This assumption is supported by the following results: (1) The mutant Deltaflv3/DeltagcvT showed a considerable longer lag phase and sometimes bleached after shifts from slow (low light, air CO(2)) to rapid (standard light, 5% CO(2)) growing conditions. (2) Photoinhibition experiments indicated a decreased ability of the mutant Deltaflv3/DeltagcvT to cope with high-light. (3) Fluorescence measurements showed that the photosynthetic electron chain is reduced in this mutant. Our data suggest that the photorespiratory 2PG-metabolism and the photoreduction of O(2), particularly that catalyzed by Flv3, cooperate during acclimation to high-light stress in cyanobacteria.

Transcriptome analysis of Traf6 function in the innate immune response of zebrafish embryos.

TRAF6 is a key player at the cross-roads of development and immunity. The analysis of its in vivo molecular function is a great challenge since severe developmental defects and early lethality caused by Traf6 deficiency in knock-out mice interfere with analyses of the immune response. In this study we have used a new strategy to analyze the function of Traf6 in a zebrafish-Salmonella infectious disease model. In our approach the effect of a Traf6 translation-blocking morpholino was titrated down to avoid developmental defects and the response to infection under these conditions was studied using the combination of microarray analysis and whole transcriptome deep sequencing. Transcriptome profiling of the traf6 knock-down allowed the identification of a gene set whose responsiveness during infection is highly dependent on Traf6. Expression trend analysis based on the resulting datasets identified nine clusters of genes with characteristic transcription response profiles, demonstrating Traf6 has a dynamic role as a positive and negative regulator. Among the Traf6-dependent genes was a large set of well known anti-microbial and inflammatory genes. Additionally, we identified several genes which were not previously linked to a response to microbial infection, such as the fertility hormone gene gnrh2 and the DNA-damage regulated autophagy modulator 1 gene dram1. With the use of the zebrafish embryo model we have now analyzed the in vivo function of Traf6 in the innate immune response without interference of adaptive immunity.

Serious complications in gene-expression studies with stress perturbation: An example of UV-exposed p53-mutant mouse embryonic fibroblasts.

Reanalysis of our UV study of p53-mutant mouse embryonic fibroblasts revealed an intriguing orchestration of massive transcriptome responses. However, close scrutiny of the data uncovered an affected mRNA/rRNA ratio, effectively inhibiting valid data analysis. UV-dose range-finding showed low-dose UV specific- and high-dose stress-related responses, which represent a plea for UV dose range-finding in experimental design.

Gene expression analysis of tuberous sclerosis complex cortical tubers reveals increased expression of adhesion and inflammatory factors.

Cortical tubers in patients with tuberous sclerosis complex are associated with disabling neurological manifestations, including intractable epilepsy. While these malformations are believed to result from the effects of TSC1 or TSC2 gene mutations, the molecular mechanisms leading to tuber formation, as well as the onset of seizures, remain largely unknown. We used the Affymetrix Gene Chip platform to provide the first genome-wide investigation of gene expression in surgically resected tubers, compared with histological normal perituberal tissue from the same patients or autopsy control tissue. We identified 2501 differentially expressed genes in cortical tubers compared with autopsy controls. Expression of genes associated with cell adhesion, for example, VCAM1, integrins and CD44, or with the inflammatory response, including complement factors, serpinA3, CCL2 and several cytokines, was increased in cortical tubers, whereas genes related to synaptic transmission, for example, the glial glutamate transporter GLT-1, and voltage-gated channel activity, exhibited lower expression. Gene expression in perituberal cortex was distinct from autopsy control cortex suggesting that even in the absence of tissue pathology the transcriptome is altered in TSC. Changes in gene expression yield insights into new candidate genes that may contribute to tuber formation or seizure onset, representing new targets for potential therapeutic development.

Membrane-associated transcripts in Arabidopsis; their isolation and characterization by DNA microarray analysis and bioinformatics.

Membrane-associated, integral membrane and secreted proteins are of key importance in many cellular processes. For most of the 28,952 predicted proteins in Arabidopsis, the actual subcellular localization has not been demonstrated experimentally. So far, their potential membrane-association has been deduced from algorithms that predict transmembrane domains and signal peptides. However, the comprehensiveness and accuracy of these algorithms is still limited. The majority of membrane-associated and secreted proteins is synthesized on membrane-bound polysomes. Therefore, the isolation and characterization of mRNA associated with membrane-bound polysomes offers an experimental tool for the genome-wide identification of these proteins. Here we describe an efficient method to isolate mRNA from membrane-bound polysomes and report on the validation of the method to enrich for transcripts encoding membrane-associated and secreted proteins. The sensitivity and reproducibility of the isolation method was investigated by DNA microarray analysis. Pearson correlations between transcript levels obtained from three replicate isolations showed that the method is highly reproducible. A significant enrichment for mRNAs encoding proteins containing predicted transmembrane domains and signal peptides was observed in the membrane-bound polysomal fraction. In this fraction, 301 transcripts were classified by gene ontologies as 'cellular component unknown', and potentially encode previously unrecognized secreted or membrane-associated proteins.

Silencing of the pollen-specific gene NTP303 and its family members in tobacco affects in vivo pollen tube growth and results in male sterile plants.

In seed plants, successful fertilization requires correct regulation of pollen tube growth. At germination and during growth, the pollen tube interacts with tissues from the pistil while the pollen tube extends via tip growth. Despite the fact that much research has been devoted to the mechanisms regulating pollen tube growth, many aspects are currently unknown. Previously, we have isolated a pollen-specific gene from tobacco--NTP303--that probably functions during pollen tube growth. NTP303 is part of a family of five members. Its expression is regulated both at the transcriptional and at the translational level. While NTP303 transcripts accumulate to high levels between early bi-cellular and mature pollen stages, NTP303 protein is hardly detectable until germination and pollen tube growth. In order to elucidate the role and function of NTP303 in the pollen tube, we studied the effect of NTP303 gene silencing on pollen function. Therefore, we have transformed tobacco plants with NTP303 co-suppression and anti-sense gene constructs. In these plants, the kanamycin resistance trait--which was linked to the NTP303-silencing gene--was not transmitted through the male gametophyte. This indicated that lowering the transcript level of NTP303 and/or its family members interferes with pollen function. Because we could not find a readily distinguishable phenotype in pollen from the hemizygous anti-sense and co-suppression plants, we rescued the defective pollen to produce doubled haploid plants that were homozygous for the NTP303 anti-sense gene. We found that in pollen from these plants the transcript levels of all NTP303 family members were reduced. Although pollen and pollen tubes from these plants appeared completely normal in vitro, the pollen tubes showed slower growth rates in vivo and arrested in the style before they reached the ovary, so that fertilization failed. These data demonstrate that NTP303 and its family members are essential for normal pollen tube growth and indicate several possible functions.

Versatile gene-specific sequence tags for Arabidopsis functional genomics: transcript profiling and reverse genetics applications.

Microarray transcript profiling and RNA interference are two new technologies crucial for large-scale gene function studies in multicellular eukaryotes. Both rely on sequence-specific hybridization between complementary nucleic acid strands, inciting us to create a collection of gene-specific sequence tags (GSTs) representing at least 21,500 Arabidopsis genes and which are compatible with both approaches. The GSTs were carefully selected to ensure that each of them shared no significant similarity with any other region in the Arabidopsis genome. They were synthesized by PCR amplification from genomic DNA. Spotted microarrays fabricated from the GSTs show good dynamic range, specificity, and sensitivity in transcript profiling experiments. The GSTs have also been transferred to bacterial plasmid vectors via recombinational cloning protocols. These cloned GSTs constitute the ideal starting point for a variety of functional approaches, including reverse genetics. We have subcloned GSTs on a large scale into vectors designed for gene silencing in plant cells. We show that in planta expression of GST hairpin RNA results in the expected phenotypes in silenced Arabidopsis lines. These versatile GST resources provide novel and powerful tools for functional genomics.