Identification of constrained sequence elements across 239 primate genomes.

Noncoding DNA is central to our understanding of human gene regulation and complex diseases1,2, and measuring the evolutionary sequence constraint can establish the functional relevance of putative regulatory elements in the human genome3-9. Identifying the genomic elements that have become constrained specifically in primates has been hampered by the faster evolution of noncoding DNA compared to protein-coding DNA10, the relatively short timescales separating primate species11, and the previously limited availability of whole-genome sequences12. Here we construct a whole-genome alignment of 239 species, representing nearly half of all extant species in the primate order. Using this resource, we identified human regulatory elements that are under selective constraint across primates and other mammals at a 5% false discovery rate. We detected 111,318 DNase I hypersensitivity sites and 267,410 transcription factor binding sites that are constrained specifically in primates but not across other placental mammals and validate their cis-regulatory effects on gene expression. These regulatory elements are enriched for human genetic variants that affect gene expression and complex traits and diseases. Our results highlight the important role of recent evolution in regulatory sequence elements differentiating primates, including humans, from other placental mammals.

Total evidence phylogeny of platyrrhine primates and a comparison of undated and tip-dating approaches.

There have been multiple published phylogenetic analyses of platyrrhine primates (New World monkeys) using both morphological and molecular data, but relatively few that have integrated both types of data into a total evidence approach. Here, we present phylogenetic analyses of recent and fossil platyrrhines, based on a total evidence data set of 418 morphological characters and 10.2 kilobases of DNA sequence data from 17 nuclear genes taken from previous studies, using undated and tip-dating approaches in a Bayesian framework. We compare the results of these analyses with molecular scaffold analyses using maximum parsimony and Bayesian approaches, and we use a formal information theoretic approach to identify unstable taxa. After a posteriori pruning of unstable taxa, the undated and tip-dating topologies appear congruent with recent molecular analyses and support largely similar relationships, with strong support for Stirtonia as a stem alouattine, Neosaimiri as a stem saimirine, Cebupithecia as a stem pitheciine, and Lagonimico as a stem callitrichid. Both analyses find three Greater Antillean subfossil platyrrhines (Xenothrix, Antillothrix, and Paralouatta) to form a clade that is related to Callicebus, congruent with a single dispersal event by the ancestor of this clade to the Greater Antilles. They also suggest that the fossil Proteropithecia may not be closely related to pitheciines, and that all known platyrrhines older than the Middle Miocene are stem taxa. Notably, the undated analysis found the Early Miocene Panamacebus (currently recognized as the oldest known cebid) to be unstable, and the tip-dating analysis placed it outside crown Platyrrhini. Our tip-dating analysis supports a late Oligocene or earliest Miocene (20.8-27.0 Ma) age for crown Platyrrhini, congruent with recent molecular clock analyses.

Two hundred and five newly assembled mitogenomes provide mixed evidence for rivers as drivers of speciation for Amazonian primates.

Mitochondrial DNA remains a cornerstone for molecular ecology, especially for study species from which high-quality tissue samples cannot be easily obtained. Methods using mitochondrial markers are usually reliant on reference databases, but these are often incomplete. Furthermore, available mitochondrial genomes often lack crucial metadata, such as sampling location, limiting their utility for many analyses. Here, we assembled 205 new mitochondrial genomes for platyrrhine primates, most from the Amazon and with known sampling locations. We present a dated mitogenomic phylogeny based on these samples along with additional published platyrrhine mitogenomes, and use this to assess support for the long-standing riverine barrier hypothesis (RBH), which proposes that river formation was a major driver of speciation in Amazonian primates. Along the Amazon, Negro, and Madeira rivers, we found mixed support for the RBH. While we identified divergences that coincide with a river barrier, only some occur synchronously and also overlap with the proposed dates of river formation. The most compelling evidence is for the Amazon river potentially driving speciation within bearded saki monkeys (Chiropotes spp.) and within the smallest extant platyrrhines, the marmosets and tamarins. However, we also found that even large rivers do not appear to be barriers for some primates, including howler monkeys (Alouatta spp.), uakaris (Cacajao spp.), sakis (Pithecia spp.), and robust capuchins (Sapajus spp.). Our results support a more nuanced, clade-specific effect of riverine barriers and suggest that other evolutionary mechanisms, besides the RBH and allopatric speciation, may have played an important role in the diversification of platyrrhines.

Antibiotic resistance profiles and population structure of disease-associated Staphylococcus aureus infecting patients in Fort Portal Regional Referral Hospital, Western Uganda.

Tackling antimicrobial resistance (AMR) is particularly challenging in low-resource settings such as Fort Portal Regional Referral Hospital (FPRRH) in Western Uganda. Specific knowledge of local AMR epidemiology is required to inform evidence-based improvement of antibiotic stewardship measures in the hospital. To address this, we combined existing antimicrobial susceptibility testing (AST) from FPRRH, with whole genome sequencing (WGS) of 41 Staphylococcus aureus isolates (2017-2019). AST revealed 73 % (30 of 41) of isolates were resistant to one or more antibiotics and 29 % (12 of 41) were multi-drug resistant (MDR). Resistance phenotypes were largely explained by the presence of antibiotic resistance genes in WGS data. Five isolates were methicillin-resistant S. aureus (MRSA) and MDR. Although all isolates were susceptible to clindamycin, a 24 % carriage of erm genes suggests potential for rapid development of resistance. We inferred a population structure for the S. aureus isolates by comparing their core genomes. Twenty isolates formed a tight cluster corresponding to multilocus sequence typing clonal complex (CC) 152, a CC found to be particularly prevalent in northern Africa. The frequency of genes associated with methicillin, chloramphenicol and ciprofloxacin resistance were significantly lower among CC152 strains than non-CC152 strains; thus, in keeping with previous work, we find that CC152 is almost exclusively methicillin-sensitive S. aureus (MSSA). Also, in agreement with other studies, we observed that the occurrence of Panton-Valentine leukocidin toxin-encoding genes was significantly higher among CC152 strains than non-CC152 strains. However, we also observed that the coagulase gene was over-represented in this CC, further defining the virulence strategy of this important pathogen. By generating detailed information about the epidemiology of circulating S. aureus and their antibiotic susceptibility, our study has provided, for the first time, data on which evidence-based infection and AMR interventions at FPRRH can be based.

Fungal microbiomes are determined by host phylogeny and exhibit widespread associations with the bacterial microbiome.

Interactions between hosts and their resident microbial communities are a fundamental component of fitness for both agents. Though recent research has highlighted the importance of interactions between animals and their bacterial communities, comparative evidence for fungi is lacking, especially in natural populations. Using data from 49 species, we present novel evidence of strong covariation between fungal and bacterial communities across the host phylogeny, indicative of recruitment by hosts for specific suites of microbes. Using co-occurrence networks, we demonstrate marked variation across host taxonomy in patterns of covariation between bacterial and fungal abundances. Host phylogeny drives differences in the overall richness of bacterial and fungal communities, but the effect of diet on richness was only evident in the mammalian gut microbiome. Sample type, tissue storage and DNA extraction method also affected bacterial and fungal community composition, and future studies would benefit from standardized approaches to sample processing. Collectively these data indicate fungal microbiomes may play a key role in host fitness and suggest an urgent need to study multiple agents of the animal microbiome to accurately determine the strength and ecological significance of host-microbe interactions.

Nitric oxide (NO) elicits aminoglycoside tolerance in Escherichia coli but antibiotic resistance gene carriage and NO sensitivity have not co-evolved.

The spread of multidrug-resistance in Gram-negative bacterial pathogens presents a major clinical challenge, and new approaches are required to combat these organisms. Nitric oxide (NO) is a well-known antimicrobial that is produced by the immune system in response to infection, and numerous studies have demonstrated that NO is a respiratory inhibitor with both bacteriostatic and bactericidal properties. However, given that loss of aerobic respiratory complexes is known to diminish antibiotic efficacy, it was hypothesised that the potent respiratory inhibitor NO would elicit similar effects. Indeed, the current work demonstrates that pre-exposure to NO-releasers elicits a > tenfold increase in IC50 for gentamicin against pathogenic E. coli (i.e. a huge decrease in lethality). It was therefore hypothesised that hyper-sensitivity to NO may have arisen in bacterial pathogens and that this trait could promote the acquisition of antibiotic-resistance mechanisms through enabling cells to persist in the presence of toxic levels of antibiotic. To test this hypothesis, genomics and microbiological approaches were used to screen a collection of E. coli clinical isolates for antibiotic susceptibility and NO tolerance, although the data did not support a correlation between increased carriage of antibiotic resistance genes and NO tolerance. However, the current work has important implications for how antibiotic susceptibility might be measured in future (i.e. ± NO) and underlines the evolutionary advantage for bacterial pathogens to maintain tolerance to toxic levels of NO.

Insights into hominid evolution from the gorilla genome sequence.

Gorillas are humans' closest living relatives after chimpanzees, and are of comparable importance for the study of human origins and evolution. Here we present the assembly and analysis of a genome sequence for the western lowland gorilla, and compare the whole genomes of all extant great ape genera. We propose a synthesis of genetic and fossil evidence consistent with placing the human-chimpanzee and human-chimpanzee-gorilla speciation events at approximately 6 and 10 million years ago. In 30% of the genome, gorilla is closer to human or chimpanzee than the latter are to each other; this is rarer around coding genes, indicating pervasive selection throughout great ape evolution, and has functional consequences in gene expression. A comparison of protein coding genes reveals approximately 500 genes showing accelerated evolution on each of the gorilla, human and chimpanzee lineages, and evidence for parallel acceleration, particularly of genes involved in hearing. We also compare the western and eastern gorilla species, estimating an average sequence divergence time 1.75 million years ago, but with evidence for more recent genetic exchange and a population bottleneck in the eastern species. The use of the genome sequence in these and future analyses will promote a deeper understanding of great ape biology and evolution.

A simple procedure for directly obtaining haplotype sequences of diploid genomes.

BACKGROUND: Almost all genome sequencing projects neglect the fact that diploid organisms contain two genome copies and consequently what is published is a composite of the two. This means that the relationship between alternate alleles at two or more linked loci is lost. We have developed a simplified method of directly obtaining the haploid sequences of each genome copy from an individual organism. RESULTS: The diploid sequences of three groups of cattle samples were obtained using a simple sample preparation procedure requiring only a microscope and a haemocytometer. Samples were: 1) lymphocytes from a single Angus steer; 2) sperm cells from an Angus bull; 3) lymphocytes from East African Zebu (EAZ) cattle collected and processed in a field laboratory in Eastern Kenya. Haploid sequence from a fosmid library prepared from lymphocytes of an EAZ cow was used for comparison. Cells were serially diluted to a concentration of one cell per microlitre by counting with a haemocytometer at each dilution. One microlitre samples, each potentially containing a single cell, were lysed and divided into six aliquots (except for the sperm samples which were not divided into aliquots). Each aliquot was amplified with phi29 polymerase and sequenced. Contigs were obtained by mapping to the bovine UMD3.1 reference genome assembly and scaffolds were assembled by joining adjacent contigs that were within a threshold distance of each other. Scaffolds that appeared to contain artefacts of CNV or repeats were filtered out leaving scaffolds with an N50 length of 27-133 kb and a 88-98 % genome coverage. SNP haplotypes were assembled with the Single Individual Haplotyper program to generate an N50 size of 97-201 kb but only ~27-68 % genome coverage. This method can be used in any laboratory with no special equipment at only slightly higher costs than conventional diploid genome sequencing. A substantial body of software for analysis and workflow management was written and is available as supplementary data. CONCLUSIONS: We have developed a set of laboratory protocols and software tools that will enable any laboratory to obtain haplotype sequences at only modestly greater cost than traditional mixed diploid sequences.

Population genomics of domestic and wild yeasts.

Since the completion of the genome sequence of Saccharomyces cerevisiae in 1996 (refs 1, 2), there has been a large increase in complete genome sequences, accompanied by great advances in our understanding of genome evolution. Although little is known about the natural and life histories of yeasts in the wild, there are an increasing number of studies looking at ecological and geographic distributions, population structure and sexual versus asexual reproduction. Less well understood at the whole genome level are the evolutionary processes acting within populations and species that lead to adaptation to different environments, phenotypic differences and reproductive isolation. Here we present one- to fourfold or more coverage of the genome sequences of over seventy isolates of the baker's yeast S. cerevisiae and its closest relative, Saccharomyces paradoxus. We examine variation in gene content, single nucleotide polymorphisms, nucleotide insertions and deletions, copy numbers and transposable elements. We find that phenotypic variation broadly correlates with global genome-wide phylogenetic relationships. S. paradoxus populations are well delineated along geographic boundaries, whereas the variation among worldwide S. cerevisiae isolates shows less differentiation and is comparable to a single S. paradoxus population. Rather than one or two domestication events leading to the extant baker's yeasts, the population structure of S. cerevisiae consists of a few well-defined, geographically isolated lineages and many different mosaics of these lineages, supporting the idea that human influence provided the opportunity for cross-breeding and production of new combinations of pre-existing variations.

Dynamic repertoire of a eukaryotic transcriptome surveyed at single-nucleotide resolution.

Recent data from several organisms indicate that the transcribed portions of genomes are larger and more complex than expected, and that many functional properties of transcripts are based not on coding sequences but on regulatory sequences in untranslated regions or non-coding RNAs. Alternative start and polyadenylation sites and regulation of intron splicing add additional dimensions to the rich transcriptional output. This transcriptional complexity has been sampled mainly using hybridization-based methods under one or few experimental conditions. Here we applied direct high-throughput sequencing of complementary DNAs (RNA-Seq), supplemented with data from high-density tiling arrays, to globally sample transcripts of the fission yeast Schizosaccharomyces pombe, independently from available gene annotations. We interrogated transcriptomes under multiple conditions, including rapid proliferation, meiotic differentiation and environmental stress, as well as in RNA processing mutants to reveal the dynamic plasticity of the transcriptional landscape as a function of environmental, developmental and genetic factors. High-throughput sequencing proved to be a powerful and quantitative method to sample transcriptomes deeply at maximal resolution. In contrast to hybridization, sequencing showed little, if any, background noise and was sensitive enough to detect widespread transcription in >90% of the genome, including traces of RNAs that were not robustly transcribed or rapidly degraded. The combined sequencing and strand-specific array data provide rich condition-specific information on novel, mostly non-coding transcripts, untranslated regions and gene structures, thus improving the existing genome annotation. Sequence reads spanning exon-exon or exon-intron junctions give unique insight into a surprising variability in splicing efficiency across introns, genes and conditions. Splicing efficiency was largely coordinated with transcript levels, and increased transcription led to increased splicing in test genes. Hundreds of introns showed such regulated splicing during cellular proliferation or differentiation.

Identification of a mutation associated with fatal Foal Immunodeficiency Syndrome in the Fell and Dales pony.

The Fell and Dales are rare native UK pony breeds at risk due to falling numbers, in-breeding, and inherited disease. Specifically, the lethal Mendelian recessive disease Foal Immunodeficiency Syndrome (FIS), which manifests as B-lymphocyte immunodeficiency and progressive anemia, is a substantial threat. A significant percentage (∼10%) of the Fell ponies born each year dies from FIS, compromising the long-term survival of this breed. Moreover, the likely spread of FIS into other breeds is of major concern. Indeed, FIS was identified in the Dales pony, a related breed, during the course of this work. Using a stepwise approach comprising linkage and homozygosity mapping followed by haplotype analysis, we mapped the mutation using 14 FIS-affected, 17 obligate carriers, and 10 adults of unknown carrier status to a ∼1 Mb region (29.8 - 30.8 Mb) on chromosome (ECA) 26. A subsequent genome-wide association study identified two SNPs on ECA26 that showed genome-wide significance after Bonferroni correction for multiple testing: BIEC2-692674 at 29.804 Mb and BIEC2-693138 at 32.19 Mb. The associated region spanned 2.6 Mb from ∼29.6 Mb to 32.2 Mb on ECA26. Re-sequencing of this region identified a mutation in the sodium/myo-inositol cotransporter gene (SLC5A3); this causes a P446L substitution in the protein. This gene plays a crucial role in the regulatory response to osmotic stress that is essential in many tissues including lymphoid tissues and during early embryonic development. We propose that the amino acid substitution we identify here alters the function of SLC5A3, leading to erythropoiesis failure and compromise of the immune system. FIS is of significant biological interest as it is unique and is caused by a gene not previously associated with a mammalian disease. Having identified the associated gene, we are now able to eradicate FIS from equine populations by informed selective breeding.

Understanding the Impact of Temperate Bacteriophages on Their Lysogens Through Transcriptomics.

Temperate phages are found integrated as prophages in the majority of bacterial genomes. Some prophages are cryptic and fixed in the bacterial chromosome, but others are active and can be triggered into a replicative form either spontaneously or by exposure to inducing factors. Prophages are commonly associated with the ability to confer toxin production or other virulence-associated traits on their host cell. More recent studies have shown they can play a much bigger role in altering the physiology of their hosts. The technique described here has enabled us to investigate how prophages affect gene expression in the opportunistic bacterium Pseudomonas aeruginosa. In this work, the growth of the wild-type P. aeruginosa strain PAO1 was compared with that of isogenic lysogens carrying different combinations of prophages from the Liverpool Epidemic Strain (LES) LESB58. In a lysogen culture, a proportion of bacterial cells will be supporting lytic bacteriophage replication (spontaneous induction) with a high level of expression per cell of late phage genes, such as those associated with the assembly of phage particles, thus masking the low-level gene expression associated with lysogen-restricted gene expression. The impact of spontaneous induction can thus obscure prophage gene expression across a lysogen population. Growth profiling experiments were used to identify spontaneous induction, which was minimal during the early exponential growth phase. This study reports how to prepare sample cultures during the early exponential growth phase and how to set up adequate controls despite low cell numbers. These protocols ensure the reliable and reproducible comparison of wild-type and lysogenic bacteria under various conditions, thus improving the transcriptomic profiling of prophage genomes and aiding in the identification of previously unrecognized prophage functions.

The landscape of tolerated genetic variation in humans and primates.

Personalized genome sequencing has revealed millions of genetic differences between individuals, but our understanding of their clinical relevance remains largely incomplete. To systematically decipher the effects of human genetic variants, we obtained whole genome sequencing data for 809 individuals from 233 primate species, and identified 4.3 million common protein-altering variants with orthologs in human. We show that these variants can be inferred to have non-deleterious effects in human based on their presence at high allele frequencies in other primate populations. We use this resource to classify 6% of all possible human protein-altering variants as likely benign and impute the pathogenicity of the remaining 94% of variants with deep learning, achieving state-of-the-art accuracy for diagnosing pathogenic variants in patients with genetic diseases. One Sentence Summary: Deep learning classifier trained on 4.3 million common primate missense variants predicts variant pathogenicity in humans.

A global catalog of whole-genome diversity from 233 primate species.

The rich diversity of morphology and behavior displayed across primate species provides an informative context in which to study the impact of genomic diversity on fundamental biological processes. Analysis of that diversity provides insight into long-standing questions in evolutionary and conservation biology and is urgent given severe threats these species are facing. Here, we present high-coverage whole-genome data from 233 primate species representing 86% of genera and all 16 families. This dataset was used, together with fossil calibration, to create a nuclear DNA phylogeny and to reassess evolutionary divergence times among primate clades. We found within-species genetic diversity across families and geographic regions to be associated with climate and sociality, but not with extinction risk. Furthermore, mutation rates differ across species, potentially influenced by effective population sizes. Lastly, we identified extensive recurrence of missense mutations previously thought to be human specific. This study will open a wide range of research avenues for future primate genomic research.

The landscape of tolerated genetic variation in humans and primates.

Personalized genome sequencing has revealed millions of genetic differences between individuals, but our understanding of their clinical relevance remains largely incomplete. To systematically decipher the effects of human genetic variants, we obtained whole-genome sequencing data for 809 individuals from 233 primate species and identified 4.3 million common protein-altering variants with orthologs in humans. We show that these variants can be inferred to have nondeleterious effects in humans based on their presence at high allele frequencies in other primate populations. We use this resource to classify 6% of all possible human protein-altering variants as likely benign and impute the pathogenicity of the remaining 94% of variants with deep learning, achieving state-of-the-art accuracy for diagnosing pathogenic variants in patients with genetic diseases.

Deletion of glyceraldehyde-3-phosphate dehydrogenase (gapN) in Clostridium saccharoperbutylacetonicum N1-4(HMT) using CLEAVE™ increases the ATP pool and accelerates solvent production.

The development and advent of mutagenesis tools for solventogenic clostridial species in recent years has allowed for the increased refinement of industrially relevant strains. In this study we have utilised CLEAVE™, a CRISPR/Cas genome editing system developed by Green Biologics Ltd., to engineer a strain of Clostridium saccharoperbutylacetonicum N1-4(HMT) with potentially useful solvents titres and energy metabolism. As one of two enzymes responsible for the conversion of glyceraldehyde-3-phosphate (GAP) to 3-phosphoglyceric acid in glycolysis, it was hypothesised that deletion of gapN would increase ATP and NADH production that could in turn improve solvent production. Herein, whole genome sequencing has been used to evaluate CLEAVE™ and the successful knockout of gapN, demonstrating a clean knockout with no other detectable variations from the wild type sequence. Elevated solvent levels were detected during the first 24 h of batch fermentation, indicating an earlier shift to solventogenesis. A 2.4-fold increase in ATP concentration was observed, and quantitation of NAD(P)H derivatives revealed a more reducing cytoplasm for the gapN strain. These findings expand our understanding of clostridium carbon metabolism and report a new approach to optimising biofuel production.

Subclonal phylogenetic structures in cancer revealed by ultra-deep sequencing.

During the clonal expansion of cancer from an ancestral cell with an initiating oncogenic mutation to symptomatic neoplasm, the occurrence of somatic mutations (both driver and passenger) can be used to track the on-going evolution of the neoplasm. All subclones within a cancer are phylogenetically related, with the prevalence of each subclone determined by its evolutionary fitness and the timing of its origin relative to other subclones. Recently developed massively parallel sequencing platforms promise the ability to detect rare subclones of genetic variants without a priori knowledge of the mutations involved. We used ultra-deep pyrosequencing to investigate intraclonal diversification at the Ig heavy chain locus in 22 patients with B-cell chronic lymphocytic leukemia. Analysis of a non-polymorphic control locus revealed artifactual insertions and deletions resulting from sequencing errors and base substitutions caused by polymerase misincorporation during PCR amplification. We developed an algorithm to differentiate genuine haplotypes of somatic hypermutations from such artifacts. This proved capable of detecting multiple rare subclones with frequencies as low as 1 in 5000 copies and allowed the characterization of phylogenetic interrelationships among subclones within each patient. This study demonstrates the potential for ultra-deep resequencing to recapitulate the dynamics of clonal evolution in cancer cell populations.

Validation of FASTFISH-ID: A new commercial platform for rapid fish species authentication via universal closed-tube barcoding.

Seafood represents up to 20% of animal protein consumption in global food consumption and is a critical dietary and income resource for the world's population. Currently, over 30% of marine fish stocks are harvested at unsustainable levels, and the industry faces challenges related to Illegal, Unregulated and Unreported (IUU) fishing. Accurate species identification is one critical component of successful stock management and helps combat fraud. Existing DNA-based technologies permit identification of seafood even when morphological features are removed, but are either too time-consuming, too expensive, or too specific for widespread use throughout the seafood supply chain. FASTFISH-ID is an innovative commercial platform for fish species authentication, employing closed-tube barcoding in a portable device. This method begins with asymmetric PCR amplification of the full length DNA barcode sequence and subsequently interrogates the resulting single-stranded DNA with a universal set of Positive/Negative probes labeled in two fluorescent colors. Each closed-tube reaction generates two species-specific fluorescent signatures that are then compared to a cloud-based library of previously validated fluorescent signatures. This novel approach results in rapid, automated species authentication without the need for complex, time consuming, identification by DNA sequencing, or repeated analysis with a panel of species-specific tests. Performance of the FASTFISH-ID platform was assessed in a blinded study carried out in three laboratories located in the UK and North America. The method exhibited a 98% success rate among the participating laboratories when compared to species identification via conventional DNA barcoding by sequencing. Thus, FASTFISH-ID is a promising new platform for combating seafood fraud across the global seafood supply chain.

Functional Genomics of a Symbiotic Community: Shared Traits in the Olive Fruit Fly Gut Microbiota.

The olive fruit fly Bactrocera oleae is a major pest of olives worldwide and houses a specialized gut microbiota dominated by the obligate symbiont "Candidatus Erwinia dacicola." Candidatus Erwinia dacicola is thought to supplement dietary nitrogen to the host, with only indirect evidence for this hypothesis so far. Here, we sought to investigate the contribution of the symbiosis to insect fitness and explore the ecology of the insect gut. For this purpose, we examined the composition of bacterial communities associated with Cretan olive fruit fly populations, and inspected several genomes and one transcriptome assembly. We identified, and reconstructed the genome of, a novel component of the gut microbiota, Tatumella sp. TA1, which is stably associated with Mediterranean olive fruit fly populations. We also reconstructed a number of pathways related to nitrogen assimilation and interactions with the host. The results show that, despite variation in taxa composition of the gut microbial community, core functions related to the symbiosis are maintained. Functional redundancy between different microbial taxa was observed for genes involved in urea hydrolysis. The latter is encoded in the obligate symbiont genome by a conserved urease operon, likely acquired by horizontal gene transfer, based on phylogenetic evidence. A potential underlying mechanism is the action of mobile elements, especially abundant in the Ca. E. dacicola genome. This finding, along with the identification, in the studied genomes, of extracellular surface structure components that may mediate interactions within the gut community, suggest that ongoing and past genetic exchanges between microbes may have shaped the symbiosis.

Large-scale and significant expression from pseudogenes in Sodalis glossinidius - a facultative bacterial endosymbiont.

The majority of bacterial genomes have high coding efficiencies, but there are some genomes of intracellular bacteria that have low gene density. The genome of the endosymbiont Sodalis glossinidius contains almost 50 % pseudogenes containing mutations that putatively silence them at the genomic level. We have applied multiple 'omic' strategies, combining Illumina and Pacific Biosciences Single-Molecule Real-Time DNA sequencing and annotation, stranded RNA sequencing and proteome analysis to better understand the transcriptional and translational landscape of Sodalis pseudogenes, and potential mechanisms for their control. Between 53 and 74 % of the Sodalis transcriptome remains active in cell-free culture. The mean sense transcription from coding domain sequences (CDSs) is four times greater than that from pseudogenes. Comparative genomic analysis of six Illumina-sequenced Sodalis isolates from different host Glossina species shows pseudogenes make up ~40 % of the 2729 genes in the core genome, suggesting that they are stable and/or that Sodalis is a recent introduction across the genus Glossina as a facultative symbiont. These data shed further light on the importance of transcriptional and translational control in deciphering host-microbe interactions. The combination of genomics, transcriptomics and proteomics gives a multidimensional perspective for studying prokaryotic genomes with a view to elucidating evolutionary adaptation to novel environmental niches.