Margaret Hughes

The Taxonomic and Functional Diversity of Microbes at a Temperate Coastal Site: A ‘Multi-Omic’ Study of Seasonal and Diel Temporal Variation

How microbial communities change over time in response to the environment is poorly understood. Previously a six-year time series of 16S rRNA V6 data from the Western English Channel demonstrated robust seasonal structure within the bacterial community, with diversity negatively correlated with day-length. Here we determine whether metagenomes and metatranscriptomes follow similar patterns. We generated 16S rRNA datasets, metagenomes (1.2 GB) and metatranscriptomes (157 MB) for eight additional time points sampled in 2008, representing three seasons (Winter, Spring, Summer) and including day and night samples. This is the first microbial ‘multi-omic’ study to combine 16S rRNA amplicon sequencing with metagenomic and metatranscriptomic profiling. Five main conclusions can be drawn from analysis of these data: 1) Archaea follow the same seasonal patterns as Bacteria, but show lower relative diversity; 2) Higher 16S rRNA diversity also reflects a higher diversity of transcripts; 3) Diversity is highest in winter and at night; 4) Community-level changes in 16S-based diversity and metagenomic profiles are better explained by seasonal patterns (with samples closest in time being most similar), while metatranscriptomic profiles are better explained by diel patterns and shifts in particular categories (i.e., functional groups) of genes; 5) Changes in key genes occur among seasons and between day and night (i.e., photosynthesis); but these samples contain large numbers of orphan genes without known homologues and it is these unknown gene sets that appear to contribute most towards defining the differences observed between times. Despite the huge diversity of these microbial communities, there are clear signs of predictable patterns and detectable stability over time. Renewed and intensified efforts are required to reveal fundamental deterministic patterns in the most complex microbial communities. Further, the presence of a substantial proportion of orphan sequences underscores the need to determine the gene products of sequences with currently unknown function.

Analysis of gene expression from the Wolbachia genome of a filarial nematode supports both metabolic and defensive roles within the symbiosis

The α-proteobacterium Wolbachia is probably the most prevalent, vertically transmitted symbiont on Earth. In contrast with its wide distribution in arthropods, Wolbachia is restricted to one family of animal-parasitic nematodes, the Onchocercidae. This includes filarial pathogens such as Onchocerca volvulus, the cause of human onchocerciasis, or river blindness. The symbiosis between filariae and Wolbachia is obligate, although the basis of this dependency is not fully understood. Previous studies suggested that Wolbachia may provision metabolites (e.g., haem, riboflavin, and nucleotides) and/or contribute to immune defense. Importantly, Wolbachia is restricted to somatic tissues in adult male worms, whereas females also harbor bacteria in the germline. We sought to characterize the nature of the symbiosis between Wolbachia and O. ochengi, a bovine parasite representing the closest relative of O. volvulus. First, we sequenced the complete genome of Wolbachia strain wOo, which revealed an inability to synthesize riboflavin de novo. Using RNA-seq, we also generated endobacterial transcriptomes from male soma and female germline. In the soma, transcripts for membrane transport and respiration were up-regulated, while the gonad exhibited enrichment for DNA replication and translation. The most abundant Wolbachia proteins, as determined by geLC-MS, included ligands for mammalian Toll-like receptors. Enzymes involved in nucleotide synthesis were dominant among metabolism-related proteins, whereas the haem biosynthetic pathway was poorly represented. We conclude that Wolbachia may have a mitochondrion-like function in the soma, generating ATP for its host. Moreover, the abundance of immunogenic proteins in wOo suggests a role in diverting the immune system toward an ineffective antibacterial response.

Post-genomic approaches to understanding the mechanisms of environmentally induced phenotypic plasticity

SUMMARY Post-genomic techniques offer new and detailed insights into the mechanisms underpinning all biological processes, including phenotypic plasticity and environmentally relevant phenotypes. Although they require access to genomic resources it is now possible to create these for species of comparative or environmental interest even within a modest research project. Here we describe an open transcript screen for genes responding to environmental cold that might account for the acquired cold-specific phenotype in all its complex manifestations. Construction of a cDNA microarray led to a survey of transcript expression levels in seven tissues of carp, as a function of time, and three different extents of cooling. The resulting data delineated a common stress response found in all tissues that comprises genes involved in cellular homeostasis, including energy charge, ATP turnover, protein turnover and stress protein production. These genes respond to kinds of perturbation other than cold and probably form part of a more general stress response common to other species. We also defined tissue-specific response patterns of transcript regulation whose main characteristics were investigated by a profiling technique based on categorisation of gene function. These genes underpin the highly tissue-specific pattern of physiological adaptations observed in the cold-acclimated fish. As a result we have identified a large number of candidate gene targets with which to investigate adaptive responses to environmental challenge.

Characteristics of the genome of Arsenophonus nasoniae, son-killer bacterium of the wasp Nasonia

We report the properties of a draft genome sequence of the bacterium Arsenophonus nasoniae, son‐killer bacterium of Nasonia vitripennis. The genome sequence data from this study are the first for a male‐killing bacterium, and represent a microorganism that is unusual compared with other sequenced symbionts, in having routine vertical and horizontal transmission, two alternating hosts, and being culturable on cell‐free media. The resulting sequence totals c. 3.5 Mbp and is annotated to contain 3332 predicted open reading frames (ORFs). Therefore, Arsenophonus represents a relatively large genome for an insect symbiont. The annotated ORF set suggests that the microbe is capable of a broad array of metabolic functions, well beyond those found for reproductive parasite genomes sequenced to date and more akin to horizontally transmitted and secondary symbionts. We also find evidence of genetic transfer from Wolbachia symbionts, and phage exchange with other gammaproteobacterial symbionts. These findings reflect the complex biology of a bacterium that is able to live, invade and survive multiple host environments while resisting immune responses.

Identification of Three Novel Superantigen-Encoding Genes in Streptococcus equi subsp. zooepidemicus, szeF, szeN, and szeP

ABSTRACT The acquisition of superantigen-encoding genes by Streptococcus pyogenes has been associated with increased morbidity and mortality in humans, and the gain of four superantigens by Streptococcus equi is linked to the evolution of this host-restricted pathogen from an ancestral strain of the opportunistic pathogen Streptococcus equi subsp. zooepidemicus. A recent study determined that the culture supernatants of several S. equi subsp. zooepidemicus strains possessed mitogenic activity but lacked known superantigen-encoding genes. Here, we report the identification and activities of three novel superantigen-encoding genes. The products of szeF, szeN, and szeP share 59%, 49%, and 34% amino acid sequence identity with SPEH, SPEM, and SPEL, respectively. Recombinant SzeF, SzeN, and SzeP stimulated the proliferation of equine peripheral blood mononuclear cells, and tumor necrosis factor alpha (TNF-α) and gamma interferon (IFN-γ) production, in vitro. Although none of these superantigen genes were encoded within functional prophage elements, szeN and szeP were located next to a prophage remnant, suggesting that they were acquired by horizontal transfer. Eighty-one of 165 diverse S. equi subsp. zooepidemicus strains screened, including 7 out of 15 isolates from cases of disease in humans, contained at least one of these new superantigen-encoding genes. The presence of szeN or szeP, but not szeF, was significantly associated with mitogenic activity in the S. equi subsp. zooepidemicus population (P < 0.000001, P < 0.000001, and P = 0.104, respectively). We conclude that horizontal transfer of these novel superantigens from and within the diverse S. equi subsp. zooepidemicus population is likely to have implications for veterinary and human disease.

A De Novo Expression Profiling of Anopheles funestus, Malaria Vector in Africa, Using 454 Pyrosequencing

Background Anopheles funestus is one of the major malaria vectors in Africa and yet there are few genomic tools available for this species compared to An. gambiae. To start to close this knowledge gap, we sequenced the An. funestus transcriptome using cDNA libraries developed from a pyrethroid resistant laboratory strain and a pyrethroid susceptible field strain from Mali. Results Using a pool of life stages (pupae, larvae, adults: females and males) for each strain, 454 sequencing generated 375,619 reads (average length of 182 bp). De novo assembly generated 18,103 contigs with average length of 253 bp. The average depth of coverage of these contigs was 8.3. In total 20.8% of all reads were novel when compared to reference databases. The sequencing of the field strain generated 204,758 reads compared to 170,861 from the insecticide resistant laboratory strain. The contigs most differentially represented in the resistant strain belong to the P450 gene family and cuticular genes which correlates with previous studies implicating both of these gene families in pyrethroid resistance. qPCR carried out on six contigs indicates that these ESTs could be suitable for gene expression studies such as microarray. 31,000 sites were estimated to contain Single Nucleotide Polymorphisms (SNPs) and analysis of SNPs from 20 contigs suggested that most of these SNPs are likely to be true SNPs. Gene conservation analysis confirmed the close phylogenetic relationship between An. funestus and An. gambiae. Conclusion This study represents a significant advance for the genetics and genomics of An. funestus since it provides an extensive set of both Expressed Sequence Tags (ESTs) and SNPs which can be readily adopted for the design of new genomic tools such as microarray or SNP platforms.

Gene expression profiling: metatranscriptomics.

Metatranscriptomics has been developed to help understand how communities respond to changes in their environment. Metagenomic studies provided a snapshot of the genetic composition of the community at any given time. However, short-timescale studies investigating the response of communities to rapid environmental changes (e.g. pollution events or diurnal light availability) require analysis of changes in the abundance and composition of the active fraction of the community. Metatranscriptomics enables researchers to investigate the actively transcribed ribosomal and messenger RNA from a community. It has been applied to environments as diverse as soil and seawater. This chapter outlines sampling protocols and RNA extraction techniques from these two ecosystems, as well as details a method to enrich mRNA in the extracted nucleic acid. Also, a section is dedicated for outlining a bioinformatic procedure for the analysis of metatranscriptomic datasets.

Population genomics of the killer whale indicates ecotype evolution in sympatry involving both selection and drift

The evolution of diversity in the marine ecosystem is poorly understood, given the relatively high potential for connectivity, especially for highly mobile species such as whales and dolphins. The killer whale (Orcinus orca) has a worldwide distribution, and individual social groups travel over a wide geographic range. Even so, regional populations have been shown to be genetically differentiated, including among different foraging specialists (ecotypes) in sympatry. Given the strong matrifocal social structure of this species together with strong resource specializations, understanding the process of differentiation will require an understanding of the relative importance of both genetic drift and local adaptation. Here we provide a high‐resolution analysis based on nuclear single‐nucleotide polymorphic markers and inference about differentiation at both neutral loci and those potentially under selection. We find that all population comparisons, within or among foraging ecotypes, show significant differentiation, including populations in parapatry and sympatry. Loci putatively under selection show a different pattern of structure compared to neutral loci and are associated with gene ontology terms reflecting physiologically relevant functions (e.g. related to digestion). The pattern of differentiation for one ecotype in the North Pacific suggests local adaptation and shows some fixed differences among sympatric ecotypes. We suggest that differential habitat use and resource specializations have promoted sufficient isolation to allow differential evolution at neutral and functional loci, but that the process is recent and dependent on both selection and drift.

Phylogenomics of the killer whale indicates ecotype divergence in sympatry

For many highly mobile species, the marine environment presents few obvious barriers to gene flow. Even so, there is considerable diversity within and among species, referred to by some as the ‘marine speciation paradox’. The recent and diverse radiation of delphinid cetaceans (dolphins) represents a good example of this. Delphinids are capable of extensive dispersion and yet many show fine-scale genetic differentiation among populations. Proposed mechanisms include the division and isolation of populations based on habitat dependence and resource specializations, and habitat release or changing dispersal corridors during glacial cycles. Here we use a phylogenomic approach to investigate the origin of differentiated sympatric populations of killer whales (Orcinus orca). Killer whales show strong specialization on prey choice in populations of stable matrifocal social groups (ecotypes), associated with genetic and phenotypic differentiation. Our data suggest evolution in sympatry among populations of resource specialists.

Characterisation of the genomes of four putative vesiculoviruses: tench rhabdovirus, grass carp rhabdovirus, perch rhabdovirus and eel rhabdovirus European X.

The complete coding sequences were determined for four putative vesiculoviruses isolated from fish. Sequence alignment and phylogenetic analysis based on the predicted amino acid sequences of the five main proteins assigned tench rhabdovirus and grass carp rhabdovirus together with spring viraemia of carp and pike fry rhabdovirus to a lineage that was distinct from the mammalian vesiculoviruses. Perch rhabdovirus, eel virus European X, lake trout rhabdovirus 903/87 and sea trout virus were placed in a second lineage that was also distinct from the recognised genera in the family Rhabdoviridae. Establishment of two new rhabdovirus genera, “Perhabdovirus” and “Sprivivirus”, is discussed.