Phylogeographic reconstruction of a bacterial species with high levels of lateral gene transfer

Background Phylogeographic reconstruction of some bacterial populations is hindered by low diversity coupled with high levels of lateral gene transfer. A comparison of recombination levels and diversity at seven housekeeping genes for eleven bacterial species, most of which are commonly cited as having high levels of lateral gene transfer shows that the relative contributions of homologous recombination versus mutation for Burkholderia pseudomallei is over two times higher than for Streptococcus pneumoniae and is thus the highest value yet reported in bacteria. Despite the potential for homologous recombination to increase diversity, B. pseudomallei exhibits a relative lack of diversity at these loci. In these situations, whole genome genotyping of orthologous shared single nucleotide polymorphism loci, discovered using next generation sequencing technologies, can provide very large data sets capable of estimating core phylogenetic relationships. We compared and searched 43 whole genome sequences of B. pseudomallei and its closest relatives for single nucleotide polymorphisms in orthologous shared regions to use in phylogenetic reconstruction. Results Bayesian phylogenetic analyses of >14,000 single nucleotide polymorphisms yielded completely resolved trees for these 43 strains with high levels of statistical support. These results enable a better understanding of a separate analysis of population differentiation among >1,700 B. pseudomallei isolates as defined by sequence data from seven housekeeping genes. We analyzed this larger data set for population structure and allele sharing that can be attributed to lateral gene transfer. Our results suggest that despite an almost panmictic population, we can detect two distinct populations of B. pseudomallei that conform to biogeographic patterns found in many plant and animal species. That is, separation along Wallace's Line, a biogeographic boundary between Southeast Asia and Australia. Conclusion We describe an Australian origin for B. pseudomallei, characterized by a single introduction event into Southeast Asia during a recent glacial period, and variable levels of lateral gene transfer within populations. These patterns provide insights into mechanisms of genetic diversification in B. pseudomallei and its closest relatives, and provide a framework for integrating the traditionally separate fields of population genetics and phylogenetics for other bacterial species with high levels of lateral gene transfer.

A Python package for Bayesian estimation using Markov Chain Monte Carlo

Markov chain Monte Carlo (MCMC) estimation provides a solution to the complex integration problems that are faced in the Bayesian analysis of statistical problems. The implementation of MCMC algorithms is, however, code intensive and time consuming. We have developed a Python package, which is called PyMCMC, that aids in the construction of MCMC samplers and helps to substantially reduce the likelihood of coding error, as well as aid in the minimisation of repetitive code. PyMCMC contains classes for Gibbs, Metropolis Hastings, independent Metropolis Hastings, random walk Metropolis Hastings, orientational bias Monte Carlo and slice samplers as well as specific modules for common models such as a module for Bayesian regression analysis. PyMCMC is straightforward to optimise, taking advantage of the Python libraries Numpy and Scipy, as well as being readily extensible with C or Fortran.

Using a Bayesian network to model colonisation with Vancomycin Resistant Enterococcus (VRE)

Objective: Effective management of multi-resistant organisms is an important issue for hospitals both in Australia and overseas. This study investigates the utility of using Bayesian Network (BN) analysis to examine relationships between risk factors and colonization with Vancomycin Resistant Enterococcus (VRE). Design: Bayesian Network Analysis was performed using infection control data collected over a period of 36 months (2008-2010). Setting: Princess Alexandra Hospital (PAH), Brisbane. Outcome of interest: Number of new VRE Isolates Methods: A BN is a probabilistic graphical model that represents a set of random variables and their conditional dependencies via a directed acyclic graph (DAG). BN enables multiple interacting agents to be studied simultaneously. The initial BN model was constructed based on the infectious disease physician‟s expert knowledge and current literature. Continuous variables were dichotomised by using third quartile values of year 2008 data. BN was used to examine the probabilistic relationships between VRE isolates and risk factors; and to establish which factors were associated with an increased probability of a high number of VRE isolates. Software: Netica (version 4.16). Results: Preliminary analysis revealed that VRE transmission and VRE prevalence were the most influential factors in predicting a high number of VRE isolates. Interestingly, several factors (hand hygiene and cleaning) known through literature to be associated with VRE prevalence, did not appear to be as influential as expected in this BN model. Conclusions: This preliminary work has shown that Bayesian Network Analysis is a useful tool in examining clinical infection prevention issues, where there is often a web of factors that influence outcomes. This BN model can be restructured easily enabling various combinations of agents to be studied.

Targeted control of feral pigs in far north Queensland : defining management units using molecular techniques

The feral pig, Sus scrofa, is a widespread and abundant invasive species in Australia. Feral pigs pose a significant threat to the environment, agricultural industry, and human health, and in far north Queensland they endanger World Heritage values of the Wet Tropics. Historical records document the first introduction of domestic pigs into Australia via European settlers in 1788 and subsequent introductions from Asia from 1827 onwards. Since this time, domestic pigs have been accidentally and deliberately released into the wild and significant feral pig populations have become established, resulting in the declaration of this species as a class 2 pest in Queensland. The overall objective of this study was to assess the population genetic structure of feral pigs in far north Queensland, in particular to enable delineation of demographically independent management units. The identification of ecologically meaningful management units using molecular techniques can assist in targeting feral pig control to bring about effective long-term management. Molecular genetic analysis was undertaken on 434 feral pigs from 35 localities between Tully and Innisfail. Seven polymorphic and unlinked microsatellite loci were screened and fixation indices (FST and analogues) and Bayesian clustering methods were used to identify population structure and management units in the study area. Sequencing of the hyper-variable mitochondrial control region (D-loop) of 35 feral pigs was also examined to identify pig ancestry. Three management units were identified in the study at a scale of 25 to 35 km. Even with the strong pattern of genetic structure identified in the study area, some evidence of long distance dispersal and/or translocation was found as a small number of individuals exhibited ancestry from a management unit outside of which they were sampled. Overall, gene flow in the study area was found to be influenced by environmental features such as topography and land use, but no distinct or obvious natural or anthropogenic geographic barriers were identified. Furthermore, strong evidence was found for non-random mating between pigs of European and Asian breeds indicating that feral pig ancestry influences their population genetic structure. Phylogenetic analysis revealed two distinct mitochondrial DNA clades, representing Asian domestic pig breeds and European breeds. A significant finding was that pigs of Asian origin living in Innisfail and south Tully were not mating randomly with European breed pigs populating the nearby Mission Beach area. Feral pig control should be implemented in each of the management units identified in this study. The control should be coordinated across properties within each management unit to prevent re-colonisation from adjacent localities. The adjacent rainforest and National Park Estates, as well as the rainforest-crop boundary should be included in a simultaneous control operation for greater success.

Molecular characterisation and expression analysis of interferon gamma in response to natural chlamydia infection in the koala, phascolarctos cinereus

Interferon gamma (IFNγ) is a key Th1 cytokine, with a principal role in the immune response against intracellular organisms such as Chlamydia. Along with being responsible for significant morbidity in human populations, Chlamydia is also responsible for wide spread infection and disease in many animal hosts, with reports that many Australian koala subpopulations are endemically infected. An understanding of the role played by IFNγ in koala chlamydial diseases is important for the establishment of better prophylactic and therapeutic approaches against chlamydial infection in this host. A limited number of IFNγ sequences have been published from marsupials and no immune reagents to measure expression have been developed. Through preliminary analysis of the koala transcriptome, we have identified the full coding sequence of the koala IFNγ gene. Transcripts were identified in spleen and lymph node tissue samples. Phylogenetic analysis demonstrated that koala IFNγ is closely related to other marsupial IFNγ sequences and more distantly related to eutherian mammals. To begin to characterise the role of this important cytokine in the koala's response to chlamydial infection, we developed a quantitative real time PCR assay and applied it to a small cohort of koalas with and without active chlamydial disease, revealing significant differences in expression patterns between the groups. Description of the IFNγ sequence from the koala will not only assist in understanding this species' response to its most important pathogen but will also provide further insight into the evolution of the marsupial immune system

Molecular analysis of the acinetobacter baumannii biofilm-associated protein

Acinetobacter baumannii is a multidrug-resistant pathogen associated with hospital outbreaks of infection across the globe, particularly in the intensive care unit. The ability of A. baumannii to survive in the hospital environment for long periods is linked to antibiotic resistance and its capacity to form biofilms. Here we studied the prevalence, expression, and function of the A. baumannii biofilm-associated protein (Bap) in 24 carbapenem-resistant A. baumannii ST92 strains isolated from a single institution over a 10-year period. The bap gene was highly prevalent, with 22/24 strains being positive for bap by PCR. Partial sequencing of bap was performed on the index case strain MS1968 and revealed it to be a large and highly repetitive gene approximately 16 kb in size. Phylogenetic analysis employing a 1,948-amino-acid region corresponding to the C terminus of Bap showed that BapMS1968 clusters with Bap sequences from clonal complex 2 (CC2) strains ACICU, TCDC-AB0715, and 1656-2 and is distinct from Bap in CC1 strains. By using overlapping PCR, the bapMS1968 gene was cloned, and its expression in a recombinant Escherichia coli strain resulted in increased biofilm formation. A Bap-specific antibody was generated, and Western blot analysis showed that the majority of A. baumannii strains expressed an ∼200-kDa Bap protein. Further analysis of three Bap-positive A. baumannii strains demonstrated that Bap is expressed at the cell surface and is associated with biofilm formation. Finally, biofilm formation by these Bap-positive strains could be inhibited by affinity-purified Bap antibodies, demonstrating the direct contribution of Bap to biofilm growth by A. baumannii clinical isolates.

Molecular analysis of type 3 fimbrial genes from Escherichia coli, Klebsiella and Citrobacter species

Background Catheter-associated urinary tract infection (CAUTI) is the most common nosocomial infection in the United States and is caused by a range of uropathogens. Biofilm formation by uropathogens that cause CAUTI is often mediated by cell surface structures such as fimbriae. In this study, we characterised the genes encoding type 3 fimbriae from CAUTI strains of Escherichia coli, Klebsiella pneumoniae, Klebsiella oxytoca, Citrobacter koseri and Citrobacter freundii. Results Phylogenetic analysis of the type 3 fimbrial genes (mrkABCD) from 39 strains revealed they clustered into five distinct clades (A-E) ranging from one to twenty-three members. The majority of sequences grouped in clade A, which was represented by the mrk gene cluster from the genome sequenced K. pneumoniae MGH78578. The E. coli and K. pneumoniae mrkABCD gene sequences clustered together in two distinct clades, supporting previous evidence for the occurrence of inter-genera lateral gene transfer. All of the strains examined caused type 3 fimbriae mediated agglutination of tannic acid treated human erythrocytes despite sequence variation in the mrkD-encoding adhesin gene. Type 3 fimbriae deletion mutants were constructed in 13 representative strains and were used to demonstrate a direct role for type 3 fimbriae in biofilm formation. Conclusions The expression of functional type 3 fimbriae is common to many Gram-negative pathogens that cause CAUTI and is strongly associated with biofilm growth. Our data provides additional evidence for the spread of type 3 fimbrial genes by lateral gene transfer. Further work is now required to substantiate the clade structure reported here by examining more strains as well as other bacterial genera that make type 3 fimbriae and cause CAUTI.

Echolocation calls, wing shape, diet and phylogenetic diagnosis of the endemic Chinese bat Myotis pequinius

We describe the echolocation calls, flight morphology and diet of the endemic Chinese bat Myotis pequinius Thomas, 1908. Orientation calls are broadband, and reach low terminal frequencies. Diet comprised 80% beetles by volume. Wing shape and call design suggest that the bats fly in cluttered habitats, and the possession of moderately long ears and the dietary composition imply they forage at least sometimes by gleaning. Myotis pequinius resembles a larger Oriental version of the western Palaearctic species M. nattereri. Phylogenetic analysis based on sequences of the cytochrome b gene of mitochondrial DNA (1,140 base pairs) from a range of Palaearctic Myotis species confirmed that M. pequinius is close to the nattereri group, and is a sister-species to the eastern Palaearctic M. bombinus. One bat sequenced from China could not be identified from available species descriptions. It was smaller than M. pequinius, and also differed from it in sequence divergence by 6.7%, suggesting the existence of additional, cryptic taxonomic diversity in this group. Our phylogenetic analysis also supports the recognition of M. schaubi as a species distinct from M. nattereri in Transcaucasia and south-western Asia. Myotis nattereri tschuliensis is more closely related to M. schaubi than to M. nattereri, and is best considered either as a subspecies of M. schaubi, or possibly as a distinct species.

Comparative analysis and distribution of Omega-3 lcPUFA Biosynthesis genes in marine molluscs

Recent research has identified marine molluscs as an excellent source of omega-3 long-chain polyunsaturated fatty acids (lcPUFAs), based on their potential for endogenous synthesis of lcPUFAs. In this study we generated a representative list of fatty acyl desaturase (Fad) and elongation of very long-chain fatty acid (Elovl) genes from major orders of Phylum Mollusca, through the interrogation of transcriptome and genome sequences, and various publicly available databases. We have identified novel and uncharacterised Fad and Elovl sequences in the following species: Anadara trapezia, Nerita albicilla, Nerita melanotragus, Crassostrea gigas, Lottia gigantea, Aplysia californica, Loligo pealeii and Chlamys farreri. Based on alignments of translated protein sequences of Fad and Elovl genes, the haeme binding motif and histidine boxes of Fad proteins, and the histidine box and seventeen important amino acids in Elovl proteins, were highly conserved. Phylogenetic analysis of aligned reference sequences was used to reconstruct the evolutionary relationships for Fad and Elovl genes separately. Multiple, well resolved clades for both the Fad and Elovl sequences were observed, suggesting that repeated rounds of gene duplication best explain the distribution of Fad and Elovl proteins across the major orders of molluscs. For Elovl sequences, one clade contained the functionally characterised Elovl5 proteins, while another clade contained proteins hypothesised to have Elovl4 function. Additional well resolved clades consisted only of uncharacterised Elovl sequences. One clade from the Fad phylogeny contained only uncharacterised proteins, while the other clade contained functionally characterised delta-5 desaturase proteins. The discovery of an uncharacterised Fad clade is particularly interesting as these divergent proteins may have novel functions. Overall, this paper presents a number of novel Fad and Elovl genes suggesting that many mollusc groups possess most of the required enzymes for the synthesis of lcPUFAs.

Phylogenetic relationships of unclassified, satellitic Pasteurellaceae obtained from different species of birds as demonstrated by 16S rRNA gene sequence comparison

Avian haemophili demonstrating in vitro satellitic growth, also referred to as the V-factor or NAD requirement, have mainly been classified with Avibacterium paragallinarum (Haemophilus paragallinarum), Avibacterium avium (Pasteurella avium), Avibacterium volantium (Pasteurella volantium) and Avibacterium sp. A (Pasteurella species A). The aim of the present study was to assess the taxonomic position of 18 V-factor-requiring isolates of unclassified Haemophilus-like organisms isolated from galliforme, anseriforme, columbiforme and gruiforme birds as well as kestrels and psittacine birds including budgerigars by conventional phenotypic tests and 16S rRNA gene sequencing. All isolates shared phenotypical characteristics which allowed classification with Pasteurellaceae. Haemolysis of bovine red blood cells was negative. Haemin (X-factor) was not required for growth. Maximum-likelihood phylogenetic analysis including bootstrap analysis showed that six isolates were related to the avian 16S rRNA group and were classified as Avibacterium according to 16S rRNA sequence analysis. Surprisingly, the other 12 isolates were unrelated to Avibacterium. Two isolates were unrelated to any of the known 16S rRNA groups of Pasteurellaceae. Two isolates were related to Volucribacter of the avian 16S rRNA group. Seven isolates belonged to the Testudinis 16S rRNA group and out of these, two isolates were closely related to taxa 14 and 32 of Bisgaard, whereas four other isolates were found to form a genus-like group distantly related to taxon 40 and one isolated remained distantly related to other members of the Testudinis group. One isolate was closely related to taxon 26 (a member of Actinobacillus sensu stricto). The study documented major genetic diversity among V-factor-requiring avian isolates beyond the traditional interpretation that they only belong to Avibacterium, underlining the limited value of satellitic growth for identification of avian members of Pasteurellaceae. Our study also emphasized that these organisms will never be isolated without the use of special media satisfying the V-factor requirement.

An inordinate fondness for Fusarium: Phylogenetic diversity of fusaria cultivated by ambrosia beetles in the genus Euwallacea on avocado and other plant hosts

Ambrosia beetle fungiculture represents one of the most ecologically and evolutionarily successful symbioses, as evidenced by the 11 independent origins and 3500 species of ambrosia beetles. Here we document the evolution of a clade within Fusarium associated with ambrosia beetles in the genus Euwallacea (Coleoptera: Scolytinae). Ambrosia Fusarium Clade (AFC) symbionts are unusual in that some are plant pathogens that cause significant damage in naive natural and cultivated ecosystems, and currently threaten avocado production in the United States, Israel and Australia. Most AFC fusaria produce unusual clavate macroconidia that serve as a putative food source for their insect mutualists. AFC symbionts were abundant in the heads of four Euwallacea spp., which suggests that they are transported within and from the natal gallery in mandibular mycangia. In a four-locus phylogenetic analysis, the AFC was resolved in a strongly supported monophyletic group within the previously described Cade 3 of the Fusarium solani species complex (FSSC). Divergence-time estimates place the origin of the AFC in the early Miocene similar to 21.2 Mya, which coincides with the hypothesized adaptive radiation of the Xyleborini. Two strongly supported clades within the AFC (Clades A and B) were identified that include nine species lineages associated with ambrosia beetles, eight with Euwallacea spp. and one reportedly with Xyleborus ferrugineus, and two lineages with no known beetle association. More derived lineages within the AFC showed fixation of the clavate (club-shaped) macroconidial trait, while basal lineages showed a mix of clavate and more typical fusiform macroconidia. AFC lineages consisted mostly of genetically identical individuals associated with specific insect hosts in defined geographic locations, with at least three interspecific hybridization events inferred based on discordant placement in individual gene genealogies and detection of recombinant loci. Overall, these data are consistent with a strong evolutionary trend toward obligate symbiosis coupled with secondary contact and interspecific hybridization. (C) 2013 Elsevier Inc. All rights reserved.

Analysis of Puumala hantavirus in a bank vole population in northern Finland: evidence for co-circulation of two genetic lineages and frequent reassortment between strains

"In this study, for the first time, two distinct genetic lineages of Puumala virus (PUUV) were found within a small sampling area and within a single host genetic lineage (Ural mtDNA) at Pallasjarvi, northern Finland. Lung tissue samples of 171 bank voles (Myodes glareolus) trapped in September 1998 were screened for the presence of PUUV nucleocapsid antigen and 25 were found to be positive. Partial sequences of the PUUV small (S), medium (M) and large (L) genome segments were recovered from these samples using RT-PCR. Phylogenetic analysis revealed two genetic groups of PUUV sequences that belonged to the Finnish and north Scandinavian lineages. This presented a unique opportunity to study inter-lineage reassortment in PUUV; indeed, 32% of the studied bank voles appeared to carry reassortant virus genomes. Thus, the frequency of inter-lineage reassortment in PUUV was comparable to that of intra-lineage reassortment observed previously (Razzauti, M., Plyusnina, A., Henttonen, H. & Plyusnin, A. (2008). J Gen Virol 89, 1649-1660). Of six possible reassortant S/M/L combinations, only two were found at Pallasjarvi and, notably, in all reassortants, both S and L segments originated from the same genetic lineage, suggesting a non-random pattern for the reassortment. These findings are discussed in connection to PUUV evolution in Fermoscandia."

Comparative Nucleotide and Amino Acid Sequence Analysis of the Sequence-Specific RNA-Binding Rotavirus Nonstructural Protein NSP3

NSP3, an acidic nonstructural protein, encoded by gene 7 has been implicated as the key player in the assembly of the 11 viral plus-strand RNAs into the early replication intermediates during rotavirus morphogenesis. To date, the sequence or NSP3 from only three animal rotaviruses (SA11, SA114F, and bovine UK) has been determined and that from a human strain has not been reported. To determine the genetic diversity among gene 7 alleles from group A rotaviruses, the nucleotide sequence of the NSP3 gene from 13 strains belonging to nine different G serotypes, from both humans and animals, has been determined. Based on the amino acid sequence identity as well as phylogenetic analysis, NSP3 from group A rotaviruses falls into three evolutionarily related groups, i.e., the SA11 group, the Wa group, and the S2 group. The SA 11/SA114F gene appears to have a distant ancestral origin from that of the others and codes for a polypeptide of 315 amino acids (aa) in length. NSP3 from all other group A rotaviruses is only 313 aa in length because of a 2-amino-acid deletion near the carboxy-terminus, While the SA114F gene has the longest 3' untranslated region (UTR) of 132 nucleotides, that from other strains suffered deletions of varying lengths at two positions downstream of the translational termination codon. In spite of the divergence of the nucleotide (nt) sequence in the protein coding region, a stretch of about 80 nt in the 3' UTR is highly conserved in the NSP3 gene from all the strains. This conserved sequence in the 3' UTR might play an important role in the regulation of expression of the NSP3 gene. (C) 1995 Academic Press, Inc.

Molecular Identification and Phylogeny of Bactrocera Species (Diptera: Tephritidae)

Fruit flies that belong to the genus Bactrocera (Diptera: Tephritidae) are major invasive pests of agricultural crops in Asia and Australia. Increased transboundary movement of agricultural produce has resulted in the chance introduction of many invasive species that include Bactrocera mainly as immature stages. Therefore quick and accurate species diagnosis is important at the port of entry, where morphological identification has a limited role, as it requires the presence of adult specimens and the availability of a specialist. Unfortunately when only immature stages are present, a lacunae in their taxonomy impedes accurate species diagnosis. At this juncture, molecular species diagnostics based on COX-I have become handy, because diagnosis is not limited by developmental stages. Yet another method of quick and accurate species diagnosis for Bactrocera spp. is based on the development of species-specific markers. This study evaluated the utility of COX-I for the quick and accurate species diagnosis of eggs, larvae, pupae and adults of B. zonata Saunders, B. tau Walker, and B. dorsalis Hendel. Furthermore the utility of species-specific markers in differentiating B. zonata (500bp) and B. tau (220bp) was shown. Phylogenetic relationships among five subgenera, viz., Austrodacus, Bactrocera, Daculus, Notodacus and Zeugodacus have been resolved employing the 5' region of COX-I (1490-2198); where COX-I sequences for B. dorsalis Hendel, B. tau Walker, B. correcta Bezzi and B. zonata Saunders from India were compared with other NCBI-GenBank accessions. Phylogenetic analysis employing Maximum Parsimony (MP) and Bayesian phylogenetic approach (BP) showed that the subgenus Bactrocera is monophyletic.

A comprehensive analysis of the chorion locus in silkmoth

Despite more than 40 years of intense study, essential features of the silkmoth chorion (eggshell) are still not fully understood. To determine the precise structure of the chorion locus, we performed extensive EST analysis, constructed a bacterial artificial chromosome (BAC) contig, and obtained a continuous genomic sequence of 871,711 base pairs. We annotated 127 chorion genes in two segments interrupted by a 164 kb region with 5 non-chorion genes, orthologs of which were on chorion bearing scaffolds in 4 ditrysian families. Detailed transcriptome analysis revealed expression throughout choriogenesis of most chorion genes originally categorized as ``middle'', and evidence for diverse regulatory mechanisms including cis-elements, alternative splicing and promoter utilization, and antisense RNA. Phylogenetic analysis revealed multigene family associations and faster evolution of early chorion genes and transcriptionally active pseudogenes. Proteomics analysis identified 99 chorion proteins in the eggshell and micropyle localization of 1 early and 6 Hc chorion proteins.