Combining spatial-temporal and phylogenetic analysis approaches for improved understanding on global H5N1 transmission

Background Since late 2003, the highly pathogenic influenza A H5N1 had initiated several outbreak waves that swept across the Eurasia and Africa continents. Getting prepared for reassortment or mutation of H5N1 viruses has become a global priority. Although the spreading mechanism of H5N1 has been studied from different perspectives, its main transmission agents and spread route problems remain unsolved. Methodology/Principal Findings Based on a compilation of the time and location of global H5N1 outbreaks from November 2003 to December 2006, we report an interdisciplinary effort that combines the geospatial informatics approach with a bioinformatics approach to form an improved understanding on the transmission mechanisms of H5N1 virus. Through a spherical coordinate based analysis, which is not conventionally done in geographical analyses, we reveal obvious spatial and temporal clusters of global H5N1 cases on different scales, which we consider to be associated with two different transmission modes of H5N1 viruses. Then through an interdisciplinary study of both geographic and phylogenetic analysis, we obtain a H5N1 spreading route map. Our results provide insight on competing hypotheses as to which avian hosts are responsible for the spread of H5N1. Conclusions/Significance We found that although South China and Southeast Asia may be the virus pool of avian flu, East Siberia may be the source of the H5N1 epidemic. The concentration of migratory birds from different places increases the possibility of gene mutation. Special attention should be paid to East Siberia, Middle Siberia and South China for improved surveillance of H5N1 viruses and monitoring of migratory birds.

Molecular characterization and phylogenetic analysis of small ruminant lentiviruses isolated from Canadian sheep and goats

Background: Small Ruminant Lentiviruses (SRLV) are widespread in Canadian sheep and goats and represent an important health issue in these animals. There is however no data about the genetic diversity of Caprine Arthritis Encephalitis Virus (CAEV) or Maedi Visna Virus (MVV) in this country. Findings: We performed a molecular and phylogenetic analysis of sheep and goat lentiviruses from a small geographic area in Canada using long sequences from the gag region of 30 infected sheep and 36 infected goats originating from 14 different flocks. Pairwise DNA distance and phylogenetic analyses revealed that all SRLV sequences obtained from sheep clustered tightly with prototypical Maedi visna sequences from America. Similarly, all SRLV strains obtained from goats clustered tightly with prototypical US CAEV-Cork strain. Conclusions: The data reported in this study suggests that Canadian and US SRLV strains share common origins. In addition, the molecular data failed to bring to light any evidence of past cross species transmission between sheep and goats, which is consistent with the type of farming practiced in this part of the country where single species flocks predominate and where opportunities of cross species transmissions are proportionately low.

Phylogenetic analysis of the gag region encoding the matrix protein of small ruminant lentiviruses: comparative analysis and molecular epidemiological applications

Little sequence information exists on the matrix-protein (MA) encoding region of small ruminant lentiviruses (SRLV). Fifty-two novel sequences were established and permitted a first phylogenetic analysis of this region of the SRLV genome. The variability of the MA encoding region is higher compared to the gag region encoding the capsid protein and surprisingly close to that reported for the env gene. In contrast to primate lentiviruses, the deduced amino acid sequences of the N- and C-terminal domains of MA are variable. This permitted to pinpoint a basic domain in the N-terminal domain that is conserved in all lentiviruses and likely to play an important functional role. Additionally, a seven amino acid insertion was detected in all MVV strains, which may be used to differentiate CAEV and MVV isolates. A molecular epidemiology analysis based on these sequences indicates that the Italian lentivirus strains are closely related to each other and to the CAEV-CO strain, a prototypic strain isolated three decades ago in the US. This suggests a common origin of the SRLV circulating in the monitored flocks, possibly related to the introduction of infected goats in a negative population. Finally, this study shows that the MA region is suitable for phylogenetic studies and may be applied to monitor SRLV eradication programs.

Phylogenetic analysis of "Candidatus Mycoplasma turicensis" isolates from pet cats in the United Kingdom, Australia, and South Africa, with analysis of risk factors for infection

Two hemotropic mycoplasmas have been recognized in cats, Mycoplasma haemofelis and "Candidatus Mycoplasma haemominutum." We recently described a third feline hemoplasma species, designated "Candidatus Mycoplasma turicensis," in a Swiss cat with hemolytic anemia. This isolate induced anemia after experimental transmission to two specific-pathogen-free cats and analysis of the 16S rRNA gene revealed its close relationship to rodent hemotropic mycoplasmas. The agent was recently shown to be prevalent in Swiss pet cats. We sought to investigate the presence and clinical importance of "Candidatus Mycoplasma turicensis" infection in pet cats outside of Switzerland and to perform the molecular characterization of isolates from different countries. A "Candidatus Mycoplasma turicensis"-specific real-time PCR assay was applied to blood samples from 426 United Kingdom (UK), 147 Australian, and 69 South African pet cats. The 16S rRNA genes of isolates from different countries were sequenced and signalment and laboratory data for the cats were evaluated for associations with "Candidatus Mycoplasma turicensis" infection. Infections were detected in samples from UK, Australian, and South African pet cats. Infection was associated with the male gender, and "Candidatus Mycoplasma haemominutum" and M. haemofelis coinfection. Coinfected cats exhibited significantly lower packed cell volume (PCV) values than uninfected cats. Phylogenetic analyses revealed that some Australian and South African "Candidatus Mycoplasma turicensis" isolates branched away from the remaining isolates. In summary, "Candidatus Mycoplasma turicensis" infection in pet cats exists over a wide geographical area and significantly decreased PCV values are observed in cats coinfected with other feline hemoplasmas.

Prediction of whole-genome DNA-DNA similarity, determination of G+C content and phylogenetic analysis within the family Pasteurellaceae by multilocus sequence analysis (MLSA)

Genome predictions based on selected genes would be a very welcome approach for taxonomic studies, including DNA-DNA similarity, G+C content and representative phylogeny of bacteria. At present, DNA-DNA hybridizations are still considered the gold standard in species descriptions. However, this method is time-consuming and troublesome, and datasets can vary significantly between experiments as well as between laboratories. For the same reasons, full matrix hybridizations are rarely performed, weakening the significance of the results obtained. The authors established a universal sequencing approach for the three genes recN, rpoA and thdF for the Pasteurellaceae, and determined if the sequences could be used for predicting DNA-DNA relatedness within the family. The sequence-based similarity values calculated using a previously published formula proved most useful for species and genus separation, indicating that this method provides better resolution and no experimental variation compared to hybridization. By this method, cross-comparisons within the family over species and genus borders easily become possible. The three genes also serve as an indicator of the genome G+C content of a species. A mean divergence of around 1 % was observed from the classical method, which in itself has poor reproducibility. Finally, the three genes can be used alone or in combination with already-established 16S rRNA, rpoB and infB gene-sequencing strategies in a multisequence-based phylogeny for the family Pasteurellaceae. It is proposed to use the three sequences as a taxonomic tool, replacing DNA-DNA hybridization.

Phylogenetic analysis of Prevotella nigrescens, Prevotella intermedia and Porphyromonas gingivalis clinical strains reveals a clear species clustering

Prevotella nigrescens, Prevotella intermedia and Porphyromonas gingivalis are oral pathogens from the family Bacteroidaceae, regularly isolated from cases of gingivitis and periodontitis. In this study, the phylogenetic variability of these three bacterial species was investigated by means of 16S rRNA (rrs) gene sequence comparisons of a set of epidemiologically and geographically diverse isolates. For each of the three species, the rrs gene sequences of 11 clinical isolates as well as the corresponding type strains was determined. Comparison of all rrs sequences obtained with those of closely related species revealed a clear clustering of species, with only a little intraspecies variability but a clear difference in the rrs gene with respect to the next related taxon. The results indicate that the three species form stable, homogeneous genetic groups, which favours an rrs-based species identification of these oral pathogens. This is especially useful given the 7% sequence divergence between Prevotella intermedia and Prevotella nigrescens, since phenotypic distinction between the two Prevotella species is inconsistent or involves techniques not applicable in routine identification.

Phylogenetic analysis of Pasteurella multocida subspecies and molecular identification of feline P. multocida subsp. septica by 16S rRNA gene sequencing

Pasteurella multocida is commonly found in the oral cavity of cats and dogs. In humans it is known as an opportunistic pathogen after bites from these animals. Phenotypic identification of P. multocida based on biochemical reactions is often limited and usually only done on a species level, even though 3 subspecies are described. For molecular taxonomy and diagnostic purposes a phylogenetic analysis of the three subspecies of P. multocida based on their 16S rRNA (rrs) gene sequence was therefore carried out. We found P. multocida subsp. septica on a distinguished branch on the phylogenetic tree of Pasteurellaceae, due to a 1.5% divergence of its rrs gene compared to the two other, more closely related subspecies multocida and gallicida. This phylogenetic divergence can be used for the identification of P. multocida subsp. septica by rrs gene determination since they form a phylogenetically well isolated and defined group as shown with a set of feline isolates. Comparison to routine phenotypic identification shows the advantage of the sequence-based identification over conventional methods. It is therefore helpful for future unambiguous identification and molecular taxonomy of P. multocida as well as for epidemiological investigations.

Phylogenetic analysis of receptor FgfrL1 shows divergence of the C-terminal end in rodents

FGFRL1 is a member of the fibroblast growth factor receptor (FGFR) family. Similar to the classical receptors FGFR1-FGFR4, it contains three extracellular Ig-like domains and a single transmembrane domain. However, it lacks the intracellular tyrosine kinase domain that would be required for signal transduction, but instead contains a short intracellular tail with a peculiar histidine-rich motif. This motif has been conserved during evolution from mollusks to echinoderms and vertebrates. Only the sequences of FgfrL1 from a few rodents diverge at the C-terminal region from the canonical sequence, as they appear to have suffered a frameshift mutation within the histidine-rich motif. This mutation is observed in mouse, rat and hamster, but not in the closely related rodents mole rat (Nannospalax) and jerboa (Jaculus), suggesting that it has occurred after branching of the Muridae and Cricetidae from the Dipodidae and Spalacidae. The consequence of the frameshift is a deletion of a few histidine residues and an extension of the C-terminus by about 40 unrelated amino acids. A similar frameshift mutation has also been observed in a human patient with a craniosynostosis syndrome as well as in several patients with colorectal cancer and bladder tumors, suggesting that the histidine-rich motif is prone to mutation. The reason why this motif was conserved during evolution in most species, but not in mice, is not clear.

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

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

Phylogenetic relationships of German heavy draught horse breeds inferred from mitochondrial DNA D-loop variation

We analysed a 610-bp mitochondrial (mt)DNA D-loop fragment in a sample of German draught horse breeds and compared the polymorphic sites with sequences from Arabian, Hanoverian, Exmoor, Icelandic, Sorraia and Przewalski's Horses as well as with Suffolk, Shire and Belgian horses. In a total of 65 horses, 70 polymorphic sites representing 47 haplotypes were observed. The average percentage of polymorphic sites was 11.5% for the mtDNA fragment analysed. In the nine different draught horse breeds including South German, Mecklenburg, Saxon Thuringa coldblood, Rhenisch German, Schleswig Draught Horse, Black Forest Horse, Shire, Suffolk and Belgian, 61 polymorphic sites and 24 haplotypes were found. The phylogenetic analysis failed to show monophyletic groups for the draught horses. The analysis indicated that the draught horse populations investigated consist of diverse genetic groups with respect to their maternal lineage.

Phylogeny and prediction of genetic similarity of Cronobacter and related taxa by multilocus sequence analysis (MLSA)

Multilocus sequence analysis (MLSA) based on recN, rpoA and thdF genes was done on more than 30 species of the family Enterobacteriaceae with a focus on Cronobacter and the related genus Enterobacter. The sequences provide valuable data for phylogenetic, taxonomic and diagnostic purposes. Phylogenetic analysis showed that the genus Cronobacter forms a homogenous cluster related to recently described species of Enterobacter, but distant to other species of this genus. Combining sequence information on all three genes is highly representative for the species' %GC-content used as taxonomic marker. Sequence similarity of the three genes and even of recN alone can be used to extrapolate genetic similarities between species of Enterobacteriaceae. Finally, the rpoA gene sequence, which is the easiest one to determine, provides a powerful diagnostic tool to identify and differentiate species of this family. The comparative analysis gives important insights into the phylogeny and genetic relatedness of the family Enterobacteriaceae and will serve as a basis for further studies and clarifications on the taxonomy of this large and heterogeneous family.

Phylogenetic relationship of equine Actinobacillus species and distribution of RTX toxin genes among clusters

Equine Actinobacillus species were analysed phylogenetically by 16S rRNA gene (rrs) sequencing focusing on the species Actinobacillus equuli, which has recently been subdivided into the non-haemolytic A. equuli subsp. equuli and the haemolytic A. equuli subsp. haemolyticus. In parallel we determined the profile for RTX toxin genes of the sample of strains by PCR testing for the presence of the A. equuli haemolysin gene aqx, and the toxin genes apxI, apxII, apxIII and apxIV, which are known in porcine pathogens such as Actinobacillus pleuropneumoniae and Actinobacillus suis. The rrs-based phylogenetic analysis revealed two distinct subclusters containing both A. equuli subsp. equuli and A. equuli subsp. haemolyticus distributed through both subclusters with no correlation to taxonomic classification. Within one of the rrs-based subclusters containing the A. equuli subsp. equuli type strain, clustered as well the porcine Actinobacillus suis strains. This latter is known to be also phenotypically closely related to A. equuli. The toxin gene analysis revealed that all A. equuli subsp. haemolyticus strains from both rrs subclusters specifically contained the aqx gene while the A. suis strains harboured the genes apxI and apxII. The aqx gene was found to be specific for A. equuli subsp. haemolyticus, since A. equuli subsp. equuli contained no aqx nor any of the other RTX genes tested. The specificity of aqx for the haemolytic equine A. equuli and ApxI and ApxII for the porcine A. suis indicates a role of these RTX toxins in host species predilection of the two closely related species of bacterial pathogens and allows PCR based diagnostic differentiation of the two.

Virological and phylogenetic characterization of attenuated small ruminant lentivirus isolates eluding efficient serological detection

Three field isolates of small ruminant lentiviruses (SRLVs) were derived from a mixed flock of goats and sheep certified for many years as free of caprine arthritis encephalitis virus (CAEV). The phylogenetic analysis of pol sequences permitted to classify these isolates as A4 subtype. None of the animals showed clinical signs of SRLV infection, confirming previous observations which had suggested that this particular subtype is highly attenuated, at least for goats. A quantitative real time PCR strategy based on primers and probes derived from a highly variable env region permitted us to classify the animals as uninfected, singly or doubly infected. The performance of different serological tools based on this classification revealed their profound inadequacy in monitoring animals infected with this particular SRLV subtype. In vitro, the isolates showed differences in their cytopathicity and a tendency to replicate more efficiently in goat than sheep cells, especially in goat macrophages. By contrast, in vivo, these viruses reached significantly higher viral loads in sheep than in goats. Both env subtypes infected goats and sheep with equal efficiency. One of these, however, reached significantly higher viral loads in both species. In conclusion, we characterized three isolates of the SRLV subtype A4 that efficiently circulate in a mixed herd of goats and sheep in spite of their apparent attenuation and a strict physical separation between goats and sheep. The poor performance of the serological tools applied indicates that, to support an SRLV eradication campaign, it is imperative to develop novel, subtype specific tools.

Social Meets Molecular: Combining Phylogenetic and Latent Class Analyses to Understand HIV-1 Transmission in Switzerland

Switzerland has a complex human immunodeficiency virus (HIV) epidemic involving several populations. We examined transmission of HIV type 1 (HIV-1) in a national cohort study. Latent class analysis was used to identify socioeconomic and behavioral groups among 6,027 patients enrolled in the Swiss HIV Cohort Study between 2000 and 2011. Phylogenetic analysis of sequence data, available for 4,013 patients, was used to identify transmission clusters. Concordance between sociobehavioral groups and transmission clusters was assessed in correlation and multiple correspondence analyses. A total of 2,696 patients were infected with subtype B, 203 with subtype C, 196 with subtype A, and 733 with recombinant subtypes (mainly CRF02_AG and CRF01_AE). Latent class analysis identified 8 patient groups. Most transmission clusters of subtype B were shared between groups of gay men (groups 1-3) or between the heterosexual groups "heterosexual people of lower socioeconomic position" (group 4) and "injection drug users" (group 8). Clusters linking homosexual and heterosexual groups were associated with "older heterosexual and gay people on welfare" (group 5). "Migrant women in heterosexual partnerships" (group 6) and "heterosexual migrants on welfare" (group 7) shared non-B clusters with groups 4 and 5. Combining approaches from social and molecular epidemiology can provide insights into HIV-1 transmission and inform the design of prevention strategies.