The draft genome of the pest tephritid fruit fly Bactrocera tryoni: resources for the genomic analysis of hybridising species

Background

The tephritid fruit flies include a number of economically important pests of horticulture, with a large accumulated body of research on their biology and control. Amongst the Tephritidae, the genus Bactrocera, containing over 400 species, presents various species groups of potential utility for genetic studies of speciation, behaviour or pest control. In Australia, there exists a triad of closely-related, sympatric Bactrocera species which do not mate in the wild but which, despite distinct morphologies and behaviours, can be force-mated in the laboratory to produce fertile hybrid offspring. To exploit the opportunities offered by genomics, such as the efficient identification of genetic loci central to pest behaviour and to the earliest stages of speciation, investigators require genomic resources for future investigations.

Results

We produced a draft de novo genome assembly of Australia’s major tephritid pest species, Bactrocera tryoni. The male genome (650 -700 Mbp) includes approximately 150Mb of interspersed repetitive DNA sequences and 60Mb of satellite DNA. Assessment using conserved core eukaryotic sequences indicated 98% completeness. Over 16,000 MAKER-derived gene models showed a large degree of overlap with other Dipteran reference genomes. The sequence of the ribosomal RNA transcribed unit was also determined. Unscaffolded assemblies of B. neohumeralis and B. jarvisi were then produced; comparison with B. tryoni showed that the species are more closely related than any Drosophila species pair. The similarity of the genomes was exploited to identify 4924 potentially diagnostic indels between the species, all of which occur in non-coding regions.

Conclusions

This first draft B. tryoni genome resembles other dipteran genomes in terms of size and putative coding sequences. For all three species included in this study, we have identified a comprehensive set of non-redundant repetitive sequences, including the ribosomal RNA unit, and have quantified the major satellite DNA families. These genetic resources will facilitate the further investigations of genetic mechanisms responsible for the behavioural and morphological differences between these three species and other tephritids. We have also shown how whole genome sequence data can be used to generate simple diagnostic tests between very closely-related species where only one of the species is scaffolded.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1153) contains supplementary material, which is available to authorized users.     

Enhanced methods for unbiased deep sequencing of Lassa and Ebola RNA viruses from clinical and biological samples

We have developed a robust RNA sequencing method for generating complete de novo assemblies with intra-host variant calls of Lassa and Ebola virus genomes in clinical and biological samples. Our method uses targeted RNase H-based digestion to remove contaminating poly(rA) carrier and ribosomal RNA. This depletion step improves both the quality of data and quantity of informative reads in unbiased total RNA sequencing libraries. We have also developed a hybrid-selection protocol to further enrich the viral content of sequencing libraries. These protocols have enabled rapid deep sequencing of both Lassa and Ebola virus and are broadly applicable to other viral genomics studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0519-7) contains supplementary material, which is available to authorized users.     

Acidobacteria appear to dominate the microbiome of two sympatric Caribbean Sponges and one Zoanthid

Background

Marine invertebrate-associated microbial communities are interesting examples of complex symbiotic systems and are a potential source of biotechnological products.

Results

In this work, pyrosequencing-based assessment from bacterial community structures of sediments, two sponges, and one zoanthid collected in the Mexican Caribbean was performed. The results suggest that the bacterial diversity at the species level is higher in the sediments than in the animal samples. Analysis of bacterial communities’ structure showed that about two thirds of the bacterial diversity in all the samples belongs to the phyla Acidobacteria and Proteobacteria. The genus Acidobacterium appears to dominate the bacterial community in all the samples, reaching almost 80% in the sponge Hyrtios.

Conclusions

Our evidence suggests that the sympatric location of these benthonic species may lead to common bacterial structure features among their bacterial communities. The results may serve as a first insight to formulate hypotheses that lead to more extensive studies of sessile marine organisms’ microbiomes from the Mexican Caribbean.

Electronic supplementary material

The online version of this article (doi:10.1186/0717-6287-47-67) contains supplementary material, which is available to authorized users.     

Concerted Evolution in the Ribosomal RNA Cistron

Background

Gene conversion is the mechanism proposed to be responsible for the homogenization of multigene families such as the nuclear ribosomal gene clusters. This concerted evolutionary process prevents individual genes in gene clusters from accumulating mutations. The mechanism responsible for concerted evolution is not well understood but recombination during meiosis has been hypothesized to play a significant role in this homogenization. In this study we tested the hypothesis of unequal crossing over playing a significant role in gene conversion events within the ribosomal RNA cistron during meiosis, mitosis or both life stages in the fungal tree pathogen Ceratocystis manginecans.

Methods

Ceratocystis manginecans, a haploid ascomycete, reproduces homothallically and was found to have two distinct sequences within the internally transcribed spacer (ITS) region of the ribosomal RNA cistron. The different ITS types were scored using PCR-RFLP assays and chi-square analyses to determine the level of significance of the changes in the ratios of the ITS types.

Results

The relative ratios of the two ITS sequence types changed when the fungal isolates were cultured vegetatively or allowed to produced sexual structures and spores. These active changes were shown to occur more frequently during meiosis than mitosis.

Conclusion

The evidence presented provides concrete support for homogenization in the rRNA gene clusters found in this fungus and that the most reasonable explanation for this process is unequal crossing over.     

Intergenomic single nucleotide polymorphisms as a tool for bacterial artificial chromosome contig building of homoeologous Brassica napus regions

Background

Homoeologous sequences pose a particular challenge if bacterial artificial chromosome (BAC) contigs shall be established for specific regions of an allopolyploid genome. Single nucleotide polymorphisms (SNPs) differentiating between homoeologous genomes (intergenomic SNPs) may represent a suitable screening tool for such purposes, since they do not only identify homoeologous sequences but also differentiate between them.

Results

Sequence alignments between Brassica rapa (AA) and Brassica oleracea (CC) sequences mapping to corresponding regions on chromosomes A1 and C1, respectively were used to identify single nucleotide polymorphisms between the A and C genomes. A large fraction of these polymorphisms was also present in Brassica napus (AACC), an allopolyploid species that originated from hybridisation of A and C genome species. Intergenomic SNPs mapping throughout homoeologous chromosome segments spanning approximately one Mbp each were included in Illumina’s GoldenGate® Genotyping Assay and used to screen multidimensional pools of a Brassica napus bacterial artificial chromosome library with tenfold genome coverage. Based on the results of 50 SNP assays, a BAC contig for the Brassica napus A subgenome was established that spanned the entire region of interest. The C subgenome region was represented in three BAC contigs.

Conclusions

This proof-of-concept study shows that sequence resources of diploid progenitor genomes can be used to deduce intergenomic SNPs suitable for multiplex polymerase chain reaction (PCR)-based screening of multidimensional BAC pools of a polyploid organism. Owing to their high abundance and ease of identification, intergenomic SNPs represent a versatile tool to establish BAC contigs for homoeologous regions of a polyploid genome.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-560) contains supplementary material, which is available to authorized users.     

Comparative genomic analysis of Ralstonia solanacearum reveals candidate genes for host specificity

Background

Ralstonia solanacearum is a vascular soil-borne plant pathogen with an unusually broad host range. This economically destructive and globally distributed bacterium has thousands of distinct lineages within a heterogeneous and taxonomically disputed species complex. Some lineages include highly host-adapted strains (ecotypes), such as the banana Moko disease-causing strains, the cold-tolerant potato brown rot strains (also known as R3bv2) and the recently emerged Not Pathogenic to Banana (NPB) strains.

Results

These distinct ecotypes offer a robust model to study host adaptation and the emergence of ecotypes because the polyphyletic Moko strains include lineages that are phylogenetically close to the monophyletic brown rot and NPB strains. Draft genomes of eight new strains belonging to these three model ecotypes were produced to complement the eleven publicly available R. solanacearum genomes. Using a suite of bioinformatics methods, we searched for genetic and evolutionary features that distinguish ecotypes and propose specific hypotheses concerning mechanisms of host adaptation in the R. solanacearum species complex. Genome-wide, few differences were identified, but gene loss events, non-synonymous polymorphisms, and horizontal gene transfer were identified among type III effectors and were associated with host range differences.

Conclusions

This extensive comparative genomics analysis uncovered relatively few divergent features among closely related strains with contrasting biological characteristics; however, several virulence factors were associated with the emergence of Moko, NPB and brown rot and could explain host adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1474-8) contains supplementary material, which is available to authorized users.     

Analysis of the Pantoea ananatis pan-genome reveals factors underlying its ability to colonize and interact with plant,insect and vertebrate hosts

Background

Pantoea ananatis is found in a wide range of natural environments, including water, soil, as part of the epi- and endophytic flora of various plant hosts, and in the insect gut. Some strains have proven effective as biological control agents and plant-growth promoters, while other strains have been implicated in diseases of a broad range of plant hosts and humans. By analysing the pan-genome of eight sequenced P. ananatis strains isolated from different sources we identified factors potentially underlying its ability to colonize and interact with hosts in both the plant and animal Kingdoms.

Results

The pan-genome of the eight compared P. ananatis strains consisted of a core genome comprised of 3,876 protein coding sequences (CDSs) and a sizeable accessory genome consisting of 1,690 CDSs. We estimate that ~106 unique CDSs would be added to the pan-genome with each additional P. ananatis genome sequenced in the future. The accessory fraction is derived mainly from integrated prophages and codes mostly for proteins of unknown function. Comparison of the translated CDSs on the P. ananatis pan-genome with the proteins encoded on all sequenced bacterial genomes currently available revealed that P. ananatis carries a number of CDSs with orthologs restricted to bacteria associated with distinct hosts, namely plant-, animal- and insect-associated bacteria. These CDSs encode proteins with putative roles in transport and metabolism of carbohydrate and amino acid substrates, adherence to host tissues, protection against plant and animal defense mechanisms and the biosynthesis of potential pathogenicity determinants including insecticidal peptides, phytotoxins and type VI secretion system effectors.

Conclusions

P. ananatis has an ‘open’ pan-genome typical of bacterial species that colonize several different environments. The pan-genome incorporates a large number of genes encoding proteins that may enable P. ananatis to colonize, persist in and potentially cause disease symptoms in a wide range of plant and animal hosts.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-404) contains supplementary material, which is available to authorized users.     

Small RNAs from plants,bacteria and fungi within the order Hypocreales are ubiquitous in human plasma

Background

The human microbiome plays a significant role in maintaining normal physiology. Changes in its composition have been associated with bowel disease, metabolic disorders and atherosclerosis. Sequences of microbial origin have been observed within small RNA sequencing data obtained from blood samples. The aim of this study was to characterise the microbiome from which these sequences are derived.

Results

Abundant non-human small RNA sequences were identified in plasma and plasma exosomal samples. Assembly of these short sequences into longer contigs was the pivotal novel step in ascertaining their origin by BLAST searches. Most reads mapped to rRNA sequences. The taxonomic profiles of the microbes detected were very consistent between individuals but distinct from microbiomes reported at other sites. The majority of bacterial reads were from the phylum Proteobacteria, whilst for 5 of 6 individuals over 90% of the more abundant fungal reads were from the phylum Ascomycota; of these over 90% were from the order Hypocreales. Many contigs were from plants, presumably of dietary origin. In addition, extremely abundant small RNAs derived from human Y RNAs were detected.

Conclusions

A characteristic profile of a subset of the human microbiome can be obtained by sequencing small RNAs present in the blood. The source and functions of these molecules remain to be determined, but the specific profiles are likely to reflect health status. The potential to provide biomarkers of diet and for the diagnosis and prognosis of human disease is immense.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-933) contains supplementary material, which is available to authorized users.     

Stool microbiota composition is associated with the prospective risk of Plasmodium falciparum infection

Background

In humans it is unknown if the composition of the gut microbiota alters the risk of Plasmodium falciparum infection or the risk of developing febrile malaria once P. falciparum infection is established. Here we collected stool samples from a cohort composed of 195 Malian children and adults just prior to an intense P. falciparum transmission season. We assayed these samples using massively parallel sequencing of the 16S ribosomal RNA gene to identify the composition of the gut bacterial communities in these individuals. During the ensuing 6-month P. falciparum transmission season we examined the relationship between the stool microbiota composition of individuals in this cohort and their prospective risk of both P. falciparum infection and febrile malaria.

Results

Consistent with prior studies, stool microbial diversity in the present cohort increased with age, although the overall microbiota profile was distinct from cohorts in other regions of Africa, Asia and North America. Age-adjusted Cox regression analysis revealed a significant association between microbiota composition and the prospective risk of P. falciparum infection; however, no relationship was observed between microbiota composition and the risk of developing febrile malaria once P. falciparum infection was established.

Conclusions

These findings underscore the diversity of gut microbiota across geographic regions, and suggest that strategic modulation of gut microbiota composition could decrease the risk of P. falciparum infection in malaria-endemic areas, potentially as an adjunct to partially effective malaria vaccines.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1819-3) contains supplementary material, which is available to authorized users.     

Non-coding RNA gene families in the genomes of anopheline mosquitoes

Background

Only a small fraction of the mosquito species of the genus Anopheles are able to transmit malaria, one of the biggest killer diseases of poverty, which is mostly prevalent in the tropics. This diversity has genetic, yet unknown, causes. In a further attempt to contribute to the elucidation of these variances, the international “Anopheles Genomes Cluster Consortium” project (a.k.a. “16 Anopheles genomes project”) was established, aiming at a comprehensive genomic analysis of several anopheline species, most of which are malaria vectors. In the frame of the international consortium carrying out this project our team studied the genes encoding families of non-coding RNAs (ncRNAs), concentrating on four classes: microRNA (miRNA), ribosomal RNA (rRNA), small nuclear RNA (snRNA), and in particular small nucleolar RNA (snoRNA) and, finally, transfer RNA (tRNA).

Results

Our analysis was carried out using, exclusively, computational approaches, and evaluating both the primary NGS reads as well as the respective genome assemblies produced by the consortium and stored in VectorBase; moreover, the results of RNAseq surveys in cases in which these were available and meaningful were also accessed in order to obtain supplementary data, as were “pre-genomic era” sequence data stored in nucleic acid databases. The investigation included the identification and analysis, in most species studied, of ncRNA genes belonging to several families, as well as the analysis of the evolutionary relations of some of those genes in cross-comparisons to other members of the genus Anopheles.

Conclusions

Our study led to the identification of members of these gene families in the majority of twenty different anopheline taxa. A set of tools for the study of the evolution and molecular biology of important disease vectors has, thus, been obtained.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1038) contains supplementary material, which is available to authorized users.     

Safety in numbers: multiple occurrences of highly similar homologs among Azotobacter vinelandii carbohydrate metabolism proteins probably confer adaptive benefits

Background

Gene duplication and horizontal gene transfer are common processes in bacterial and archaeal genomes, and are generally assumed to result in either diversification or loss of the redundant gene copies. However, a recent analysis of the genome of the soil bacterium Azotobacter vinelandii DJ revealed an abundance of highly similar homologs among carbohydrate metabolism genes. In many cases these multiple genes did not appear to be the result of recent duplications, or to function only as a means of stimulating expression by increasing gene dosage, as the homologs were located in varying functional genetic contexts. Based on these initial findings we here report in-depth bioinformatic analyses focusing specifically on highly similar intra-genome homologs, or synologs, among carbohydrate metabolism genes, as well as an analysis of the general occurrence of very similar synologs in prokaryotes.

Results

Approximately 900 bacterial and archaeal genomes were analysed for the occurrence of synologs, both in general and among carbohydrate metabolism genes specifically. This showed that large numbers of highly similar synologs among carbohydrate metabolism genes are very rare in bacterial and archaeal genomes, and that the A. vinelandii DJ genome contains an unusually large amount of such synologs. The majority of these synologs were found to be non-tandemly organized and localized in varying but metabolically relevant genomic contexts. The same observation was made for other genomes harbouring high levels of such synologs. It was also shown that highly similar synologs generally constitute a very small fraction of the protein-coding genes in prokaryotic genomes. The overall synolog fraction of the A. vinelandii DJ genome was well above the data set average, but not nearly as remarkable as the levels observed when only carbohydrate metabolism synologs were considered.

Conclusions

Large numbers of highly similar synologs are rare in bacterial and archaeal genomes, both in general and among carbohydrate metabolism genes. However, A. vinelandii and several other soil bacteria harbour large numbers of highly similar carbohydrate metabolism synologs which seem not to result from recent duplication or transfer events. These genes may confer adaptive benefits with respect to certain lifestyles and environmental factors, most likely due to increased regulatory flexibility and/or increased gene dosage.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-192) contains supplementary material, which is available to authorized users.     

Horizontal transfer of a non-autonomous Helitron among insect and viral genomes

Background

The movement of mobile elements among species by horizontal transposon transfer (HTT) influences the evolution of genomes through the modification of structure and function. Helitrons are a relatively new lineage of DNA-based (class II) transposable elements (TEs) that propagate by rolling-circle replication, and are capable of acquiring host DNA. The rapid spread of Helitrons among animal lineages by HTT is facilitated by shuttling in viral particles or by unknown mechanisms mediated by close organism associations (e.g. between hosts and parasites).

Results

A non-autonomous Helitron independently annotated as BmHel-2 from Bombyx mori and the MITE01 element from Ostrinia nubilalis was predicted in the genomes of 24 species in the insect Order Lepidoptera. Integrated Helitrons retained ≥ 65% sequence identity over a 250 bp consensus, and were predicted to retain secondary structures inclusive of a 3′-hairpin and a 5′-subterminal inverted repeat. Highly similar Hel-2 copies were predicted in the genomes of insects and associated viruses, which along with a previous documented case of real-time virus-insect cell line transposition suggests that this Helitron has likely propagated by HTT.

Conclusions

These findings provide evidence that insect virus may mediate the HTT of Helitron-like TEs. This movement may facilitate the shuttling of DNA elements among insect genomes. Further sampling is required to determine the putative role of HTT in insect genome evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1318-6) contains supplementary material, which is available to authorized users.