Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.

Presence of potentially pathogenic Babesia sp. for human in Ixodes ricinus in Switzerland

We have designed and performed a new PCR method based on the 18S rRNA in order to individuate the presence and the identity of Babesia parasites. Out of 1159 Ixodes ricinus (Acari: Ixodidae) ticks collected in four areas of Switzerland, nine were found to contain Babesia DNA. Sequencing of the short amplicon obtained (411-452 bp) allowed the identification of three human pathogenic species: Babesia microti, B. divergens, for the first time in Switzerland, Babesia sp. EU1. We also report coinfections with B. sp. EU1-Borrelia burgdorferi sensu stricto and Babesia sp. EU1-B. afzelii.

Reliable gene expression measurements from degraded RNA by quantitative real-time PCR depend on short amplicons and a proper normalization

Quantitative reverse transcriptase real-time PCR (QRT-PCR) is a robust method to quantitate RNA abundance. The procedure is highly sensitive and reproducible as long as the initial RNA is intact. However, breaks in the RNA due to chemical or enzymatic cleavage may reduce the number of RNA molecules that contain intact amplicons. As a consequence, the number of molecules available for amplification decreases. We determined the relation between RNA fragmentation and threshold values (Ct values) in subsequent QRT-PCR for four genes in an experimental model of intact and partially hydrolyzed RNA derived from a cell line and we describe the relation between RNA integrity, amplicon size and Ct values in this biologically homogenous system. We demonstrate that degradation-related shifts of Ct values can be compensated by calculating delta Ct values between test genes and the mean values of several control genes. These delta Ct values are less sensitive to fragmentation of the RNA and are unaffected by varying amounts of input RNA. The feasibility of the procedure was demonstrated by comparing Ct values from a larger panel of genes in intact and in partially degraded RNA. We compared Ct values from intact RNA derived from well-preserved tumor material and from fragmented RNA derived from formalin-fixed, paraffin-embedded (FFPE) samples of the same tumors. We demonstrate that the relative abundance of gene expression can be based on FFPE material even when the amount of RNA in the sample and the extent of fragmentation are not known.

Intestinal Tritrichomonas foetus infection in cats in Switzerland detected by in vitro cultivation and PCR

Tritrichomonas foetus, a parasite well known for its significance as venereally transmitted pathogen in cattle, has recently been identified as a cause of chronic large-bowel diarrhea in domestic cats in the US, UK, and, more recently, also in Norway. In a period of 3 months (October to December 2007), 45 cats of Switzerland suffering from chronic diarrhea were investigated for intestinal infections, including a search for trichomonads. A commercially available in vitro culture system was used to screen for infection, complemented with a PCR and subsequent amplicon sequencing to support speciation. The PCR is based upon amplification of a sequence derived from the internal transcribed spacer region 1 (ITS1) on the ribosomal RNA gene (rRNA) using primers designed to detect a broad range of genera and species belonging to the family of Trichomonadidae. The method was furthermore adapted to the uracil DNA glycosylase (UDG) system in order to prevent carry-over contamination and it included a recombinant internal control to track for inhibitory reactions. Eleven out of the 45 cats were culture-positive, as revealed by microscopic identification of trichomonadid organisms. One of the isolates was subjected to scanning electron microscopy and findings revealed the presence of three flagella, thus placing the isolate into the gender Tritrichomonas sp. PCR and subsequent amplicon sequencing were carried out with ten of the 11 isolates. A total homology with published T. foetus sequences was confirmed in all of the cases. T. foetus therefore appears to range among those organisms that can cause chronic diarrhea in cats in Switzerland.

The Role of Acidic Fibroblast Growth Factor and Fibroblast Growth Factor Receptor 4 in High Grade Serous Carcinoma

Background: High grade serous carcinoma whether ovarian, tubal or primary peritoneal, continues to be the most lethal gynecologic malignancy in the USA. Although combination chemotherapy and aggressive surgical resection has improved survival in the past decade the majority of patients still succumb to chemo-resistant disease recurrence. It has recently been reported that amplification of 5q31-5q35.3 is associated with poor prognosis in patients with high grade serous ovarian carcinoma. Although the amplicon contains over 50 genes, it is notable for the presence of several members of the fibroblast growth factor signaling axis. In particular acidic fibroblast growth factor (FGF1) has been demonstrated to be one of the driving genes in mediating the observed prognostic effect of the amplicon in ovarian cancer patients. This study seeks to further validate the prognostic value of fibroblast growth receptor 4 (FGFR4), another candidate gene of the FGF/FGFR axis located in the same amplicon. The emphasis will be delineating the role the FGF1/FGFR4 signaling axis plays in high grade serous ovarian carcinoma; and test the feasibility of targeting the FGF1/FGFR4 axis therapeutically. Materials and Methods: Spearman and Pearson correlation studies on data generated from array CGH and transcriptome profiling analyses on 51 microdissected tumor samples were used to identify genes located on chromosome 5q31-35.3 that showed significant correlation between DNA and mRNA copy numbers. Significant correlation between FGF1 and FGFR4 DNA copy numbers was further validated by qPCR analysis on DNA isolated from 51 microdissected tumor samples. Immunolocalization and quantification of FGFR4 expression were performed on paraffin embedded tissue samples from 183 cases of high-grade serous ovarian carcinoma. The expression was then correlated with clinical data to assess impact on survival. The expression of FGF1 and FGFR4 in vitro was quantified by real-time PCR and western blotting in six high-grade serous ovarian carcinoma cell lines and compared to those in human ovarian surface epithelial cells to identify overexpression. The effect of FGF1 on these cell lines after serum starvation was quantified for in vitro cellular proliferation, migration/invasion, chemoresistance and survival utilizing a combination of commercially available colorimetric, fluorometric and electrical impedance assays. FGFR4 expression was then transiently silenced via siRNA transfection and the effects on response to FGF1, cellular proliferation, and migration were quantified. To identify relevant cellular pathways involved, responsive cell lines were transduced with different transcription response elements using the Cignal-Lenti reporter system and treated with FGF1 with and without transient FGFR4 knock down. This was followed by western blot confirmation for the relevant phosphoproteins. Anti-FGF1 antibodies and FGFR trap proteins were used to attempt inhibition of FGF mediated phenotypic changes and relevant signaling in vitro. Orthotopic intraperitoneal tumors were established in nude mice using serous cell lines that have been previously transfected with luciferase expressing constructs. The mice were then treated with FGFR trap protein. Tumor progression was then followed via bioluminescent imaging. The FGFR4 gene from 52 clinical samples was sequenced to screen for mutations. Results: FGFR4 DNA and mRNA copy numbers were significantly correlated and FGFR4 DNA copy number was significantly correlated with that of FGF1. Survival of patients with high FGFR4 expressing tumors was significantly shorter that those with low expression(median survival 28 vs 55 month p< 0.001) In a multivariate cox regression model FGFR expression significantly increased risk of death (HR 2.1, p<0.001). FGFR4 expression was significantly higher in all cell lines tested compared to HOSE, OVCA432 cell line in particular had very high expression suggesting amplification. FGF1 was also particularly overexpressed in OVCA432. FGF1 significantly increased cell survival after serum deprivation in all cell lines. Transient knock down of FGFR4 caused significant reduction in cell migration and proliferation in vitro and significantly decreased the proliferative effects of FGF1 in vitro. FGFR1, FGFR4 traps and anti-FGF1 antibodies did not show activity in vitro. OVCA432 transfected with the cignal lenti reporter system revealed significant activation of MAPK, NFkB and WNT pathways, western blotting confirmed the results. Reverse phase protein array (RPPA) analysis also showed activation of MAPK, AKT, WNT pathways and down regulation of E Cadherin. FGFR trap protein significantly reduced tumor growth in vivo in an orthotopic mouse model. Conclusions: Overexpression and amplification of several members of the FGF signaling axis present on the amplicon 5q31-35.3 is a negative prognostic indicator in high grade serous ovarian carcinoma and may drive poor survival associated with that amplicon. Activation of The FGF signaling pathway leads to downstream activation of MAPK, AKT, WNT and NFkB pathways leading to a more aggressive cancer phenotype with increased tumor growth, evasion of apoptosis and increased migration and invasion. Inhibition of FGF pathway in vivo via FGFR trap protein leads to significantly decreased tumor growth in an orthotopic mouse model.

Neurotropic astrovirus in cattle with nonsuppurative encephalitis in Europe

Encephalitis is a frequently diagnosed condition in cattle with neurological diseases. Many affected animals present with a nonsuppurative inflammatory reaction pattern in the brain. While this pattern supports a viral etiology, the causative pathogen remains unknown in a large proportion of cases. Using viral metagenomics, we identified an astrovirus (bovine astrovirus [BoAstV]-CH13) in the brain of a cow with nonsuppurative encephalitis. Additionally, BoAstV RNA was detected with reverse transcription-PCR and in situ hybridization in about one fourth (5/22 animals) of cattle with nonsuppurative encephalitis of unknown etiology. Viral RNA was found primarily in neurons and at the site of pathology. These findings support the notion that BoAstV infection is a common cause of encephalitis in cattle. Phylogenetically, BoAstV-CH13 was closely related to rare astrovirus isolates from encephalitis cases in animals and a human patient. Future research needs to be directed toward the pathogenic mechanisms, epidemiology, and potential cross-species transmission of these neurotropic astroviruses.

High-resolution microarray analysis of chromosome 20q in human colon cancer metastasis model systems

Amplification of human chromosome 20q DNA is the most frequently occurring chromosomal abnormality detected in sporadic colorectal carcinomas and shows significant correlation with liver metastases. Through comprehensive high-resolution microarray comparative genomic hybridization and microarray gene expression profiling, we have characterized chromosome 20q amplicon genes associated with human colorectal cancer metastasis in two in vitro metastasis model systems. The results revealed increasing complexity of the 20q genomic profile from the primary tumor-derived cell lines to the lymph node and liver metastasis derived cell lines. Expression analysis of chromosome 20q revealed a subset of over expressed genes residing within the regions of genomic copy number gain in all the tumor cell lines, suggesting these are Chromosome 20q copy number responsive genes. Bases on their preferential expression levels in the model system cell lines and known biological function, four of the over expressed genes mapping to the common intervals of genomic copy gain were considered the most promising candidate colorectal metastasis-associated genes. Validation of genomic copy number and expression array data was carried out on these genes, with one gene, DNMT3B, standing out as expressed at a relatively higher levels in the metastasis-derived cell lines compared with their primary-derived counterparts in both the models systems analyzed. The data provide evidence for the role of chromosome 20q genes with low copy gain and elevated expression in the clonal evolution of metastatic cells and suggests that such genes may serve as early biomarkers of metastatic potential. The data also support the utility of the combined microarray comparative genomic hybridization and expression array analysis for identifying copy number responsive genes in areas of low DNA copy gain in cancer cells. ^

Microbiota and Microsatellite genotypes of Pacific oysters stemming from three oyster bed in the Northern Wadden Sea

Background: Studies of oyster microbiomes have revealed that a limited number of microbes, including pathogens, can dominate microbial communities in host tissues such as gills and gut. Much of the bacterial diversity however remains underexplored and unexplained, although environmental conditions and host genetics have been implicated. We used 454 next generation 16S rRNA amplicon sequencing of individually tagged PCR reactions to explore the diversity of bacterial communities in gill tissue of the invasive Pacific oyster Crassostrea gigas stemming from genetically differentiated beds under ambient outdoor conditions and after a multifaceted disturbance treatment imposing stress on the host. Results: While the gill associated microbial communities in oysters were dominated by few abundant taxa (i.e. Sphingomonas, Mycoplasma) the distribution of rare bacterial groups correlated to relatedness between the hosts under ambient conditions. Exposing the host to disturbance broke apart this relationship by removing rare phylotypes thereby reducing overall microbial diversity. Shifts in the microbiome composition in response to stress did not result in a net increase in genera known to contain potentially pathogenic strains. Conclusion: The decrease in microbial diversity and the disassociation between population genetic structure of the hosts and their associated microbiome suggest that disturbance (i.e. stress) may play a significant role for the assembly of the natural microbiome. Such community shifts may in turn also feed back on the course of disease and the occurrence of mass mortality events in oyster populations.

Molecular differentiation by PCR-RFLP technique shows that Antarctic fishes of the genus Notothenia, claimed to be two species, are phenotypic varieties with the same genetic pattern (Table 1)

The taxonomy of Antarctic fishes has been predominantly based on morphological characteristics rather than on genetic criteria. A typical example is the Notothenia group, which includes N. coriiceps Richardson, 1844, N. neglecta Nybelin, 1951 and N. rossii Richardson, 1844. The Polymerase Chain Reaction and Restriction Fragment Length Polymorphism (PCR-RFLP) technique was used to determine whether N. coriiceps Richardson, 1844 and N. neglecta Nybelin, 1951 are different or whether they are the same species with morphological, physiological and behavioural variability. N. rossii was used as control. Mitochondrial DNA (mtDNA) was isolated from muscle specimens of N. coriiceps Richardson, 1844, N. neglecta Nybelin, 1951 and N. rossii, which were collected in Admiralty Bay, King George Island. The DNA was used to amplify a fragment (690 base pairs) of the mitochondrial gene coding region of NADH dehydrogenase subunit 2. Further, the amplicon was digested with the following restriction enzymes: DdeI, HindIII and RsaI. The results showed a variation of the digestion pattern of the fragment amplified between N. rossii, and N. coriiceps Richardson, 1844 or N. neglecta Nybelin, 1951. However, no differences were found between N. coriiceps Richardson, 1844 and N. neglecta Nybelin, 1951, on the grounds of the same genetic pattern shown by the two fish.

Geochemistry of pockmark sediments from the inner Oslofjord

Pockmarks are geological features that are found on the bottom of lakes and oceans all over the globe. Some are active, seeping oil or methane, while others are inactive. Active pockmarks are well studied since they harbor specialized microbial communities that proliferate on the seeping compounds. Such communities are not found in inactive pockmarks. Interestingly, inactive pockmarks are known to have different macrofaunal communities compared to the surrounding sediments. It is undetermined what the microbial composition of inactive pockmarks is and if it shows a similar pattern as the macrofauna. The Norwegian Oslo Fjord contains many inactive pockmarks and they are well suited to study the influence of these geological features on the microbial community in the sediment. Here we present a detailed analysis of the microbial communities found in three inactive pockmarks and two control samples at two core depth intervals. The communities were analyzed using high-throughput amplicon sequencing of the 16S rRNA V3 region. Microbial communities of surface pockmark sediments were indistinguishable from communities found in the surrounding seabed. In contrast, pockmark communities at 40 cm sediment depth had a significantly different community structure from normal sediments at the same depth. Statistical analysis of chemical variables indicated significant differences in the concentrations of total carbon and non-particulate organic carbon between 40 cm pockmark and reference sample sediments. We discuss these results in comparison with the taxonomic classification of the OTUs identified in our samples. Our results indicate that microbial surface sediment communities are affect by the water column, while the 40 cm communities are affect by local conditions within the sediment.

Metadata und statistic analysis of archael and bacterial sequences originating from sediments of the Håkon Mosby mud volcano (Z1-Z3, all habitats)

DNA extraction was carried out as described on the MICROBIS project pages (http://icomm.mbl.edu/microbis ) using a commercially available extraction kit. We amplified the hypervariable regions V4-V6 of archaeal and bacterial 16S rRNA genes using PCR and several sets of forward and reverse primers (http://vamps.mbl.edu/resources/primers.php). Massively parallel tag sequencing of the PCR products was carried out on a 454 Life Sciences GS FLX sequencer at Marine Biological Laboratory, Woods Hole, MA, following the same experimental conditions for all samples. Sequence reads were submitted to a rigorous quality control procedure based on mothur v30 (doi:10.1128/AEM.01541-09) including denoising of the flow grams using an algorithm based on PyroNoise (doi:10.1038/nmeth.1361), removal of PCR errors and a chimera check using uchime (doi:10.1093/bioinformatics/btr381). The reads were taxonomically assigned according to the SILVA taxonomy (SSURef v119, 07-2014; doi:10.1093/nar/gks1219) implemented in mothur and clustered at 98% ribosomal RNA gene V4-V6 sequence identity. V4-V6 amplicon sequence abundance tables were standardized to account for unequal sampling effort using 1000 (Archaea) and 2300 (Bacteria) randomly chosen sequences without replacement using mothur and then used to calculate inverse Simpson diversity indices and Chao1 richness (doi:10.2307/4615964). Bray-Curtis dissimilarities (doi:10.2307/1942268) between all samples were calculated and used for 2-dimensional non metric multidimensional scaling (NMDS) ordinations with 20 random starts (doi:10.1007/BF02289694). Stress values below 0.2 indicated that the multidimensional dataset was well represented by the 2D ordination. NMDS ordinations were compared and tested using Procrustes correlation analysis (doi:10.1007/BF02291478). All analyses were carried out with the R statistical environment and the packages vegan (available at: http://cran.r-project.org/package=vegan), labdsv (available at: http://cran.r-project.org/package=labdsv), as well as with custom R scripts. Operational taxonomic units at 98% sequence identity (OTU0.03) that occurred only once in the whole dataset were termed absolute single sequence OTUs (SSOabs; doi:10.1038/ismej.2011.132). OTU0.03 sequences that occurred only once in at least one sample, but may occur more often in other samples were termed relative single sequence OTUs (SSOrel). SSOrel are particularly interesting for community ecology, since they comprise rare organisms that might become abundant when conditions change.16S rRNA amplicons and metagenomic reads have been stored in the sequence read archive under SRA project accession number SRP042162.

Small Changes in pH Have Direct Effects on Marine Bacterial Community Composition: A Microcosm Approach

As the atmospheric CO2 concentration rises, more CO2 will dissolve in the oceans, leading to a reduction in pH. Effects of ocean acidification on bacterial communities have mainly been studied in biologically complex systems, in which indirect effects, mediated through food web interactions, come into play. These approaches come close to nature but suffer from low replication and neglect seasonality. To comprehensively investigate direct pH effects, we conducted highly-replicated laboratory acidification experiments with the natural bacterial community from Helgoland Roads (North Sea). Seasonal variability was accounted for by repeating the experiment four times (spring, summer, autumn, winter). Three dilution approaches were used to select for different ecological strategies, i.e. fast-growing or low-nutrient adapted bacteria. The pH levels investigated were in situ seawater pH (8.15-8.22), pH 7.82 and pH 7.67, representing the present-day situation and two acidification scenarios projected for the North Sea for the year 2100. In all seasons, both automated ribosomal intergenic spacer analysis and 16S ribosomal amplicon pyrosequencing revealed pH-dependent community shifts for two of the dilution approaches. Bacteria susceptible to changes in pH were different members of Gammaproteobacteria, Flavobacteriaceae, Rhodobacteraceae, Campylobacteraceae and further less abundant groups. Their specific response to reduced pH was often context-dependent. Bacterial abundance was not influenced by pH. Our findings suggest that already moderate changes in pH have the potential to cause compositional shifts, depending on the community assembly and environmental factors. By identifying pH-susceptible groups, this study provides insights for more directed, in-depth community analyses in large-scale and long-term experiments.

Attachment to the nuclear matrix mediates specific alterations in chromatin structure

The DNA in eukaryotic chromosomes is organized into a series of loops that are permanently attached at their bases to the nuclear scaffold or matrix at sequences known as scaffold-attachment or matrix-attachment regions. At present, it is not clear what effect affixation to the nuclear matrix has on chromatin architecture in important regulatory regions such as origins of replication or the promoter regions of genes. In the present study, we have investigated cell-cycle-dependent changes in the chromatin structure of a well characterized replication initiation zone in the amplified dihydrofolate reductase domain of the methotrexate-resistant Chinese hamster ovary cell line CHOC 400. Replication can initiate at any of multiple potential sites scattered throughout the 55-kilobase intergenic region in this domain, with two subregions (termed ori-β and ori-γ) being somewhat preferred. We show here that the chromatin in the ori-β and ori-γ regions undergoes dramatic alterations in micrococcal nuclease hypersensitivity as cells cross the G1/S boundary, but only in those copies of the amplicon that are affixed to the nuclear matrix. In contrast, the fine structure of chromatin in the promoter of the dihydrofolate reductase gene does not change detectably as a function of matrix attachment or cell-cycle position. We suggest that attachment of DNA to the nuclear matrix plays an important role in modulating chromatin architecture, and this could facilitate the activity of origins of replication.

Multiplex detection of four pathogenic retroviruses using molecular beacons

We describe a multiplex nucleic acid assay that identifies and determines the abundance of four different pathogenic retroviruses (HIV-1, HIV-2, and human T-lymphotrophic virus types I and II). Retroviral DNA sequences are amplified in a single, sealed tube by simultaneous PCR assays, and the resulting amplicons are detected in real time by the hybridization of four differently colored, amplicon-specific molecular beacons. The color of the fluorescence generated in the course of amplification identifies which retroviruses are present, and the number of thermal cycles required for the intensity of each color to rise significantly above background provides an accurate measure of the number of copies of each retroviral sequence that were present originally in the sample. Fewer than 10 retroviral genomes can be detected. Moreover, 10 copies of a rare retrovirus can be detected in the presence of 100,000 copies of an abundant retrovirus. Ninety-six samples can be analyzed in 3 hr on a single plate, and the use of a closed-tube format eliminates crossover contamination. Utilizing previously well characterized clinical samples, we demonstrate that each of the pathogenic retroviruses can be identified correctly and no false positives occur. This assay enables the rapid and reliable screening of donated blood and transplantable tissues.

FGF-2 regulation of neurogenesis in adult hippocampus after brain injury

Fibroblast growth factor-2 (FGF-2) promotes proliferation of neuroprogenitor cells in culture and is up-regulated within brain after injury. Using mice genetically deficient in FGF-2 (FGF-2−/− mice), we addressed the importance of endogenously generated FGF-2 on neurogenesis within the hippocampus, a structure involved in spatial, declarative, and contextual memory, after seizures or ischemic injury. BrdUrd incorporation was used to mark dividing neuroprogenitor cells and NeuN expression to monitor their differentiation into neurons. In the wild-type strain, hippocampal FGF-2 increased after either kainic acid injection or middle cerebral artery occlusion, and the numbers of BrdUrd/NeuN-positive cells significantly increased on days 9 and 16 as compared with the controls. In FGF-2−/− mice, BrdUrd labeling was attenuated after kainic acid or middle cerebral artery occlusion, as was the number of neural cells colabeled with both BrdUrd and NeuN. After FGF-2−/− mice were injected intraventricularly with a herpes simplex virus-1 amplicon vector carrying FGF-2 gene, the number of BrdUrd-labeled cells increased significantly to values equivalent to wild-type littermates after kainate seizures. These results indicate that endogenously synthesized FGF-2 is necessary and sufficient to stimulate proliferation and differentiation of neuroprogenitor cells in the adult hippocampus after brain insult.