Transcriptional changes in the hookworm, Ancylostoma caninum, during the transition from a free-living to a parasitic larva

Background

Third-stage larvae (L3) of the canine hookworm, Ancylostoma caninum, undergo arrested development preceding transmission to a host. Many of the mRNAs up-regulated at this stage are likely to encode proteins that facilitate the transition from a free-living to a parasitic larva. The initial phase of mammalian host invasion by A. caninum L3 (herein termed “activation”) can be mimicked in vitro by culturing L3 in serum-containing medium.

Methodology/Principal Findings

The mRNAs differentially transcribed between activated and non-activated L3 were identified by suppression subtractive hybridisation (SSH). The analysis of these mRNAs on a custom oligonucleotide microarray printed with the SSH expressed sequence tags (ESTs) and publicly available A. caninum ESTs (non-subtracted) yielded 602 differentially expressed mRNAs, of which the most highly represented sequences encoded members of the pathogenesis-related protein (PRP) superfamily and proteases. Comparison of these A. caninum mRNAs with those of Caenorhabditis elegans larvae exiting from developmental (dauer) arrest demonstrated unexpectedly large differences in gene ontology profiles. C. elegans dauer exiting L3 up-regulated expression of mostly intracellular molecules involved in growth and development. Such mRNAs are virtually absent from activated hookworm larvae, and instead are over-represented by mRNAs encoding extracellular proteins with putative roles in host-parasite interactions.

Conclusions/Significance

Although this should not invalidate C. elegans dauer exit as a model for hookworm activation, it highlights the limitations of this free-living nematode as a model organism for the transition of nematode larvae from a free-living to a parasitic state.     

Core and accessory genome architecture in a group of Pseudomonas aeruginosa Mu-like phages

Background

Bacteriophages that infect the opportunistic pathogen Pseudomonas aeruginosa have been classified into several groups. One of them, which includes temperate phage particles with icosahedral heads and long flexible tails, bears genomes whose architecture and replication mechanism, but not their nucleotide sequences, are like those of coliphage Mu. By comparing the genomic sequences of this group of P. aeruginosa phages one could draw conclusions about their ontogeny and evolution.

Results

Two newly isolated Mu-like phages of P. aeruginosa are described and their genomes sequenced and compared with those available in the public data banks. The genome sequences of the two phages are similar to each other and to those of a group of P. aeruginosa transposable phages. Comparing twelve of these genomes revealed a common genomic architecture in the group. Each phage genome had numerous genes with homologues in all the other genomes and a set of variable genes specific for each genome. The first group, which comprised most of the genes with assigned functions, was named “core genome”, and the second group, containing mostly short ORFs without assigned functions was called “accessory genome”. Like in other phage groups, variable genes are confined to specific regions in the genome.

Conclusion

Based on the known and inferred functions for some of the variable genes of the phages analyzed here, they appear to confer selective advantages for the phage survival under particular host conditions. We speculate that phages have developed a mechanism for horizontally acquiring genes to incorporate them at specific loci in the genome that help phage adaptation to the selective pressures imposed by the host.

Electronic supplementary material

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

Rapid Detection and Identification of Human Hookworm Infections through High Resolution Melting (HRM) Analysis

Background

Hookworm infections are still endemic in low and middle income tropical countries with greater impact on the socioeconomic and public health of the bottom billion of the world''s poorest people. In this study, a real-time polymerase chain reaction (PCR) coupled with high resolution melting-curve (HRM) analysis was evaluated for an accurate, rapid and sensitive tool for species identification focusing on the five human hookworm species.

Methods

Real-time PCR coupled with HRM analysis targeting the second internal transcribed spacer (ITS-2) of nuclear ribosomal DNA as the genetic marker was used to identify and distinguish hookworm species in human samples. Unique and distinct characteristics of HRM patterns were produced for each of the five hookworm species. The melting curves were characterized by peaks of 79.24±0.05°C and 83.00±0.04°C for Necator americanus, 79.12±0.10°C for Ancylostoma duodenale, 79.40±0.10°C for Ancylostoma ceylanicum, 79.63±0.05°C for Ancylostoma caninum and 79.70±0.14°C for Ancylostoma braziliense. An evaluation of the method''s sensitivity and specificity revealed that this assay was able to detect as low as 0.01 ng/µl hookworm DNA and amplification was only recorded for hookworm positive samples.

Conclusion

The HRM assay developed in this study is a rapid and straightforward method for the diagnosis, identification and discrimination of five human hookworms. This assay is simple compared to other probe-based genotyping methods as it does not require multiplexing, DNA sequencing or post-PCR processing. Therefore, this method offers a new alternative for rapid detection of human hookworm species.     

Distance-based assessment of the localization of functional annotations in 3D genome reconstructions

Background

Recent studies used the contact data or three-dimensional (3D) genome reconstructions from Hi-C (chromosome conformation capture with next-generation sequencing) to assess the co-localization of functional genomic annotations in the nucleus. These analyses dichotomized data point pairs belonging to a functional annotation as “close” or “far” based on some threshold and then tested for enrichment of “close” pairs. We propose an alternative approach that avoids dichotomization of the data and instead directly estimates the significance of distances within the 3D reconstruction.

Results

We applied this approach to 3D genome reconstructions for Plasmodium falciparum, the causative agent of malaria, and Saccharomyces cerevisiae and compared the results to previous approaches. We found significant 3D co-localization of centromeres, telomeres, virulence genes, and several sets of genes with developmentally regulated expression in P. falciparum; and significant 3D co-localization of centromeres and long terminal repeats in S. cerevisiae. Additionally, we tested the experimental observation that telomeres form three to seven clusters in P. falciparum and S. cerevisiae. Applying affinity propagation clustering to telomere coordinates in the 3D reconstructions yielded six telomere clusters for both organisms.

Conclusions

Distance-based assessment replicated key findings, while avoiding dichotomization of the data (which previously yielded threshold-sensitive results).

Electronic supplementary material

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

The bandit, a new DNA transposon from a hookworm-possible horizontal genetic transfer between host and parasite

Background

An enhanced understanding of the hookworm genome and its resident mobile genetic elements should facilitate understanding of the genome evolution, genome organization, possibly host-parasite co-evolution and horizontal gene transfer, and from a practical perspective, development of transposon-based transgenesis for hookworms and other parasitic nematodes.

Methodology/Principal Findings

A novel mariner-like element (MLE) was characterized from the genome of the dog hookworm, Ancylostoma caninum, and termed bandit. The consensus sequence of the bandit transposon was 1,285 base pairs (bp) in length. The new transposon was flanked by perfect terminal inverted repeats of 32 nucleotides in length with a common target site duplication TA, and it encoded an open reading frame (ORF) of 342 deduced amino acid residues. Phylogenetic comparisons confirmed that the ORF encoded a mariner-like transposase, which included conserved catalytic domains, and that the bandit transposon belonged to the cecropia subfamily of MLEs. The phylogenetic analysis also indicated that the Hsmar1 transposon from humans was the closest known relative of bandit, and that bandit and Hsmar1 constituted a clade discrete from the Tc1 subfamily of MLEs from the nematode Caenorhabditis elegans. Moreover, homology models based on the crystal structure of Mos1 from Drosophila mauritiana revealed closer identity in active site residues of the catalytic domain including Ser281, Lys289 and Asp293 between bandit and Hsmar1 than between Mos1 and either bandit or Hsmar1. The entire bandit ORF was amplified from genomic DNA and a fragment of the bandit ORF was amplified from RNA, indicating that this transposon is actively transcribed in hookworms.

Conclusions/Significance

A mariner-like transposon termed bandit has colonized the genome of the hookworm A. caninum. Although MLEs exhibit a broad host range, and are identified in other nematodes, the closest phylogenetic relative of bandit is the Hsmar1 element of humans. This surprising finding suggests that bandit was transferred horizontally between hookworm parasites and their mammalian hosts.     

Detection of copy number variations and their effects in Chinese bulls

Background

Copy number variations (CNVs) are a main source of genomic structural variations underlying animal evolution and production traits. Here, with one pure-blooded Angus bull as reference, we describe a genome-wide analysis of CNVs based on comparative genomic hybridization arrays in 29 Chinese domesticated bulls and examined their effects on gene expression and cattle growth traits.

Results

We identified 486 copy number variable regions (CNVRs), covering 2.45% of the bovine genome, in 24 taurine (Bos taurus), together with 161 ones in 2 yaks (Bos grunniens) and 163 ones in 3 buffaloes (Bubalus bubalis). Totally, we discovered 605 integrated CNVRs, with more “loss” events than both “gain” and “both” ones, and clearly clustered them into three cattle groups. Interestingly, we confirmed their uneven distributions across chromosomes, and the differences of mitochondrion DNA copy number (gain: taurine, loss: yak & buffalo). Furthermore, we confirmed approximately 41.8% (253/605) and 70.6% (427/605) CNVRs span cattle genes and quantitative trait loci (QTLs), respectively. Finally, we confirmed 6 CNVRs in 9 chosen ones by using quantitative PCR, and further demonstrated that CNVR22 had significantly negative effects on expression of PLA2G2D gene, and both CNVR22 and CNVR310 were associated with body measurements in Chinese cattle, suggesting their key effects on gene expression and cattle traits.

Conclusions

The results advanced our understanding of CNV as an important genomic structural variation in taurine, yak and buffalo. This study provides a highly valuable resource for Chinese cattle’s evolution and breeding researches.

Electronic supplementary material

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

Cyclic nucleotide-gated ion channel gene family in rice,identification, characterization and experimental analysis of expression response to plant hormones,biotic and abiotic stresses

Background

Cyclic nucleotide-gated channels (CNGCs) are Ca²⁺-permeable cation transport channels, which are present in both animal and plant systems. They have been implicated in the uptake of both essential and toxic cations, Ca²⁺ signaling, pathogen defense, and thermotolerance in plants. To date there has not been a genome-wide overview of the CNGC gene family in any economically important crop, including rice (Oryza sativa L.). There is an urgent need for a thorough genome-wide analysis and experimental verification of this gene family in rice.

Results

In this study, a total of 16 full length rice CNGC genes distributed on chromosomes 1–6, 9 and 12, were identified by employing comprehensive bioinformatics analyses. Based on phylogeny, the family of OsCNGCs was classified into four major groups (I-IV) and two sub-groups (IV-A and IV- B). Likewise, the CNGCs from all plant lineages clustered into four groups (I-IV), where group II was conserved in all land plants. Gene duplication analysis revealed that both chromosomal segmentation (OsCNGC1 and 2, 10 and 11, 15 and 16) and tandem duplications (OsCNGC1 and 2) significantly contributed to the expansion of this gene family. Motif composition and protein sequence analysis revealed that the CNGC specific domain “cyclic nucleotide-binding domain (CNBD)” comprises a “phosphate binding cassette” (PBC) and a “hinge” region that is highly conserved among the OsCNGCs. In addition, OsCNGC proteins also contain various other functional motifs and post-translational modification sites. We successively built a stringent motif: (LI-X(2)-[GS]-X-[FV]-X-G-[1]-ELL-X-W-X(12,22)-SA-X(2)-T-X(7)-[EQ]-AF-X-L) that recognizes the rice CNGCs specifically. Prediction of cis-acting regulatory elements in 5′ upstream sequences and expression analyses through quantitative qPCR demonstrated that OsCNGC genes were highly responsive to multiple stimuli including hormonal (abscisic acid, indoleacetic acid, kinetin and ethylene), biotic (Pseudomonas fuscovaginae and Xanthomonas oryzae pv. oryzae) and abiotic (cold) stress.

Conclusions

There are 16 CNGC genes in rice, which were probably expanded through chromosomal segmentation and tandem duplications and comprise a PBC and a “hinge” region in the CNBD domain, featured by a stringent motif. The various cis-acting regulatory elements in the upstream sequences may be responsible for responding to multiple stimuli, including hormonal, biotic and abiotic stresses.

Electronic supplementary material

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

Identification of a New Target of miR-16, Vacuolar Protein Sorting 4a

Rationale

The rationale was to utilize a bioinformatics approach to identify miRNA binding sites in genes with single nucleotide mutations (SNPs) to discover pathways in heart failure (HF).

Objective

The objective was to focus on the genes containing miRNA binding sites with miRNAs that were significantly altered in end-stage HF and in response to a left ventricular assist device (LVAD).

Methods and Results

BEDTools v2.14.3 was used to discriminate SNPs within predicted 3′UTR miRNA binding sites. A member of the miR-15/107 family, miR-16, was decreased in the circulation of end-stage HF patients and increased in response to a LVAD (p<0.001). MiR-16 decreased Vacuolar Protein Sorting 4a (VPS4a) expression in HEK 293T cells (p<0.01). The SNP rs16958754 was identified in the miR-15/107 family binding site of VPS4a which abolished direct binding of miR-16 to the 3′UTR of VPS4a (p<0.05). VPS4a was increased in the circulation of end-stage HF patients (p<0.001), and led to a decrease in the number of HEK 293T cells in vitro (p<0.001).

Conclusions

We provide evidence that miR-16 decreases in the circulation of end-stage HF patients and increases with a LVAD. Modeling studies suggest that miR-16 binds to and decreases expression of VPS4a. Overexpression of VPS4a decreases cell number. Together, these experiments suggest that miR-16 and VPS4a expression are altered in end-stage HF and in response to unloading with a LVAD. This signaling pathway may lead to reduced circulating cell number in HF.     

A novel model of double replications and random loss accounts for rearrangements in the Mitogenome of Samariscus latus (Teleostei: Pleuronectiformes)

Background

Although more than one thousand complete mitochondrial DNA (mtDNA) sequences have been determined in teleostean fishes, only a few gene rearrangements have been observed, and genome-scale rearrangements are even rarer. However, flatfishes (Pleuronectiformes) have been identified as having diverse types of mitochondrial gene rearrangements. It has been reported that tongue soles and the blue flounder mitogenomes exhibit different types of large-scale gene rearrangements.

Results

In the present study, the complete mitochondrial genome of another flatfish, Samariscus latus, was sequenced, and genome-scale rearrangements were observed. The genomic features of this flounder are different from those of any other studied vertebrates, including flatfish species too. The mitogenome of S. latus is characterized by the duplication and translocation of the control region (CR). The genes located between the two CRs are divided into two clusters in which their relative orders are maintained.

Conclusions

We propose a “Double Replications and Random Loss” model to explain the rearrangement events in S. latus mitogenome. This model consists of the following steps. First, the CR was duplicated and translocated. Subsequently, double replications of the mitogenome were successively initiated from the two CRs, leading to the duplication of the genes between the two CRs. Finally, one of each pair of duplicated genes was lost in a random event.

Electronic supplementary material

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

Bimodal signatures of germline methylation are linked with gene expression plasticity in the coral Acropora millepora

Background

In invertebrates, genes belonging to dynamically regulated functional categories appear to be less methylated than “housekeeping” genes, suggesting that DNA methylation may modulate gene expression plasticity. To date, however, experimental evidence to support this hypothesis across different natural habitats has been lacking.

Results

Gene expression profiles were generated from 30 pairs of genetically identical fragments of coral Acropora millepora reciprocally transplanted between distinct natural habitats for 3 months. Gene expression was analyzed in the context of normalized CpG content, a well-established signature of historical germline DNA methylation. Genes with weak methylation signatures were more likely to demonstrate differential expression based on both transplant environment and population of origin than genes with strong methylation signatures. Moreover, the magnitude of expression differences due to environment and population were greater for genes with weak methylation signatures.

Conclusions

Our results support a connection between differential germline methylation and gene expression flexibility across environments and populations. Studies of phylogenetically basal invertebrates such as corals will further elucidate the fundamental functional aspects of gene body methylation in Metazoa.

Electronic supplementary material

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

Genomic comparative analysis of the environmental Enterococcus mundtii against enterococcal representative species

Background

Enterococcus mundtii is a yellow-pigmented microorganism rarely found in human infections. The draft genome sequence of E. mundtii was recently announced. Its genome encodes at least 2,589 genes and 57 RNAs, and 4 putative genomic islands have been detected. The objective of this study was to compare the genetic content of E. mundtii with respect to other enterococcal species and, more specifically, to identify genes coding for putative virulence traits present in enterococcal opportunistic pathogens.

Results

An in-depth mining of the annotated genome was performed in order to uncover the unique properties of this microorganism, which allowed us to detect a gene encoding the antimicrobial peptide mundticin among other relevant features. Moreover, in this study a comparative genomic analysis against commensal and pathogenic enterococcal species, for which genomic sequences have been released, was conducted for the first time. Furthermore, our study reveals significant similarities in gene content between this environmental isolate and the selected enterococci strains (sharing an “enterococcal gene core” of 805 CDS), which contributes to understand the persistence of this genus in different niches and also improves our knowledge about the genetics of this diverse group of microorganisms that includes environmental, commensal and opportunistic pathogens.

Conclusion

Although E. mundtii CRL1656 is phylogenetically closer to E. faecium, frequently responsible of nosocomial infections, this strain does not encode the most relevant relevant virulence factors found in the enterococcal clinical isolates and bioinformatic predictions indicate that it possesses the lowest number of putative pathogenic genes among the most representative enterococcal species. Accordingly, infection assays using the Galleria mellonella model confirmed its low virulence.

Electronic supplementary material

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

Analysis of whole genome sequencing for the Escherichia coli O157:H7 typing phages

Background

Shiga toxin producing Escherichia coli O157 can cause severe bloody diarrhea and haemolytic uraemic syndrome. Phage typing of E. coli O157 facilitates public health surveillance and outbreak investigations, certain phage types are more likely to occupy specific niches and are associated with specific age groups and disease severity. The aim of this study was to analyse the genome sequences of 16 (fourteen T4 and two T7) E. coli O157 typing phages and to determine the genes responsible for the subtle differences in phage type profiles.

Results

The typing phages were sequenced using paired-end Illumina sequencing at The Genome Analysis Centre and the Animal Health and Veterinary Laboratories Agency and bioinformatics programs including Velvet, Brig and Easyfig were used to analyse them. A two-way Euclidian cluster analysis highlighted the associations between groups of phage types and typing phages. The analysis showed that the T7 typing phages (9 and 10) differed by only three genes and that the T4 typing phages formed three distinct groups of similar genomic sequences: Group 1 (1, 8, 11, 12 and 15, 16), Group 2 (3, 6, 7 and 13) and Group 3 (2, 4, 5 and 14). The E. coli O157 phage typing scheme exhibited a significantly modular network linked to the genetic similarity of each group showing that these groups are specialised to infect a subset of phage types.

Conclusion

Sequencing the typing phage has enabled us to identify the variable genes within each group and to determine how this corresponds to changes in phage type.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1470-z) 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.     

Resources and costs for microbial sequence analysis evaluated using virtual machines and cloud computing

Background

The widespread popularity of genomic applications is threatened by the “bioinformatics bottleneck” resulting from uncertainty about the cost and infrastructure needed to meet increasing demands for next-generation sequence analysis. Cloud computing services have been discussed as potential new bioinformatics support systems but have not been evaluated thoroughly.

Results

We present benchmark costs and runtimes for common microbial genomics applications, including 16S rRNA analysis, microbial whole-genome shotgun (WGS) sequence assembly and annotation, WGS metagenomics and large-scale BLAST. Sequence dataset types and sizes were selected to correspond to outputs typically generated by small- to midsize facilities equipped with 454 and Illumina platforms, except for WGS metagenomics where sampling of Illumina data was used. Automated analysis pipelines, as implemented in the CloVR virtual machine, were used in order to guarantee transparency, reproducibility and portability across different operating systems, including the commercial Amazon Elastic Compute Cloud (EC2), which was used to attach real dollar costs to each analysis type. We found considerable differences in computational requirements, runtimes and costs associated with different microbial genomics applications. While all 16S analyses completed on a single-CPU desktop in under three hours, microbial genome and metagenome analyses utilized multi-CPU support of up to 120 CPUs on Amazon EC2, where each analysis completed in under 24 hours for less than $60. Representative datasets were used to estimate maximum data throughput on different cluster sizes and to compare costs between EC2 and comparable local grid servers.

Conclusions

Although bioinformatics requirements for microbial genomics depend on dataset characteristics and the analysis protocols applied, our results suggests that smaller sequencing facilities (up to three Roche/454 or one Illumina GAIIx sequencer) invested in 16S rRNA amplicon sequencing, microbial single-genome and metagenomics WGS projects can achieve cost-efficient bioinformatics support using CloVR in combination with Amazon EC2 as an alternative to local computing centers.