Molecular characterisation of environmental enterococci derived from water samples and assessment of associated public health hazards

Enterococci are versatile Gram-positive bacteria that can survive under extreme conditions. Most enterococci are non-virulent and found in the gastrointestinal tract of humans and animals. Other strains are opportunistic pathogens that contribute to a large number of nosocomial infections globally. Epidemiological studies demonstrated a direct relationship between the density of enterococci in surface waters and the risk of swimmer-associated gastroenteritis. The distribution of infectious enterococcal strains from the hospital environment or other sources to environmental water bodies through sewage discharge or other means, could increase the prevalence of these strains in the human population. Environmental water quality studies may benefit from focusing on a subset of Enterococcus spp. that are consistently associated with sources of faecal pollution such as domestic sewage, rather than testing for the entire genus. E. faecalis and E. faecium are potentially good focal species for such studies, as they have been consistently identified as the dominant Enterococcus spp. in human faeces and sewage. On the other hand enterococcal infections are predominantly caused by E. faecalis and E. faecium. The characterisation of E. faecalis and E. faecium is important in studying their population structures, particularly in environmental samples. In developing and implementing rapid, robust molecular genotyping techniques, it is possible to more accurately establish the relationship between human and environmental enterococci. Of particular importance, is to determine the distribution of high risk enterococcal clonal complexes, such as E. faecium clonal complex 17 and E. faecalis clonal complexes 2 and 9 in recreational waters. These clonal complexes are recognized as particularly pathogenic enterococcal genotypes that cause severe disease in humans globally. The Pimpama-Coomera watershed is located in South East Queensland, Australia and was investigated in this study mainly because it is used intensively for agriculture and recreational purposes and has a strong anthropogenic impact. The primary aim of this study was to develop novel, universally applicable, robust, rapid and cost effective genotyping methods which are likely to yield more definitive results for the routine monitoring of E. faecalis and E. faecium, particularly in environmental water sources. To fullfill this aim, new genotyping methods were developed based on the interrogation of highly informative single nucleotide polymorphisms (SNPs) located in housekeeping genes of both E. faecalis and E. faecium. SNP genotyping was successfully applied in field investigations of the Coomera watershed, South-East Queensland, Australia. E. faecalis and E. faecium isolates were grouped into 29 and 23 SNP profiles respectively. This study showed the high longitudinal diversity of E. faecalis and E. faecium over a period of two years, and both human-related and human-specific SNP profiles were identified. Furthermore, 4.25% of E. faecium strains isolated from water was found to correspond to the important clonal complex-17 (CC17). Strains that belong to CC17 cause the majority of hospital outbreaks and clinical infections globally. Of the six sampling sites of the Coomera River, Paradise Point had the highest number of human-related and human-specific E. faecalis and E. faecium SNP profiles. The secondary aim of this study was to determine the antibiotic-resistance profiles and virulence traits associated with environmental E. faecalis and E. faecium isolates compared to human pathogenic E. faecalis and E. faecium isolates. This was performed to predict the potential health risks associated with coming into contact with these strains in the Coomera watershed. In general, clinical isolates were found to be more resistant to all the antibiotics tested compared to water isolates and they harbored more virulence traits. Multi-drug resistance was more prevalent in clinical isolates (71.18% of E. faecalis and 70.3 % of E. faecium) compared to water isolates (only 5.66 % E. faecium). However, tetracycline, gentamicin, ciprofloxacin and ampicillin resistance was observed in water isolates. The virulence gene esp was the most prevalent virulence determinant observed in clinical isolates (67.79% of E. faecalis and 70.37 % of E. faecium), and this gene has been described as a human-specific marker used for microbial source tracking (MST). The presence of esp in water isolates (16.36% of E. faecalis and 19.14% of E. faecium) could be indicative of human faecal contamination in these waterways. Finally, in order to compare overall gene expression between environmental and clinical strains of E. faecalis, a comparative gene hybridization study was performed. The results of this investigation clearly demonstrated the up-regulation of genes associated with pathogenicity in E. faecalis isolated from water. The expression study was performed at physiological temperatures relative to ambient temperatures. The up-regulation of virulence genes demonstrates that environmental strains of E. faecalis can pose an increased health risk which can lead to serious disease, particularly if these strains belong to the virulent CC17 group. The genotyping techniques developed in this study not only provide a rapid, robust and highly discriminatory tool to characterize E. faecalis and E. faecium, but also enables the efficient identification of virulent enterococci that are distributed in environmental water sources.

A molecular phylogeny for the Tribe Dacini (Diptera: Tephritidae) : systematic and biogeographic implications

With well over 700 species, the Tribe Dacini is one of the most species-rich clades within the dipteran family Tephritidae, the true fruit flies. Nearly all Dacini belong to one of two very large genera, Dacus Fabricius and Bactrocera Macquart. The distribution of the genera overlap in or around the Indian subcontinent, but the greatest diversity of Dacus is in Africa and the greatest diversity of Bactrocera is in south-east Asia and the Pacific. The monophyly of these two genera has not been rigorously established, with previous phylogenies only including a small number of species and always heavily biased to one genus over the other. Moreover, the subgeneric taxonomy within both genera is complex and the monophyly of many subgenera has not been explicitly tested. Previous hypotheses about the biogeography of the Dacini based on morphological reviews and current distributions of taxa have invoked an out-of-India hypothesis; however this has not been tested in a phylogenetic framework. We attempted to resolve these issues with a dated, molecular phylogeny of 125 Dacini species generated using 16S, COI, COII and white eye genes. The phylogeny shows that Bactrocera is not monophyletic, but rather consists of two major clades: Bactrocera s.s. and the ‘Zeugodacus group of subgenera’ (a recognised, but informal taxonomic grouping of 15 Bactrocera subgenera). This ‘Zeugodacus’ clade is the sister group to Dacus, not Bactrocera and, based on current distributions, split from Dacus before that genus moved into Africa. We recommend that taxonomic consideration be given to raising Zeugodacus to genus level. Supportive of predictions following from the out-of-India hypothesis, the first common ancestor of the Dacini arose in the mid-Cretaceous approximately 80 mya. Major divergence events occurred during the Indian rafting period and diversification of Bactrocera apparently did not begin until after India docked with Eurasia (50–35 mya). In contrast, diversification in Dacus, at approximately 65 mya, apparently began much earlier than predicted by the out-of-India hypothesis, suggesting that, if the Dacini arose on the Indian plate, then ancestral Dacus may have left the plate in the mid to late Cretaceous via the well documented India–Madagascar–Africa migration route. We conclude that the phylogeny does not disprove the predictions of an out-of-India hypothesis for the Dacini, although modification of the original hypothesis is required.

A mitochondrial genome phylogeny of termites (Blattodea: Termitoidae) : robust support for interfamilial relationships and molecular synapomorphies define major clades

Despite their ecological significance as decomposers and their evolutionary significance as the most speciose eusocial insect group outside the Hymenoptera, termite (Blattodea: Termitoidae or Isoptera) evolutionary relationships have yet to be well resolved. Previous morphological and molecular analyses strongly conflict at the family level and are marked by poor support for backbone nodes. A mitochondrial (mt) genome phylogeny of termites was produced to test relationships between the recognised termite families, improve nodal support and test the phylogenetic utility of rare genomic changes found in the termite mt genome. Complete mt genomes were sequenced for 7 of the 9 extant termite families with additional representatives of each of the two most speciose families Rhinotermitidae (3 of 7 subfamilies) and Termitidae (3 of 8 subfamilies). The mt genome of the well supported sister group of termites, the subsocial cockroach Cryptocercus, was also sequenced. A highly supported tree of termite relationships was produced by all analytical methods and data treatment approaches, however the relationship of the termites + Cryptocercus clade to other cockroach lineages was highly affected by the strong nucleotide compositional bias found in termites relative to other dictyopterans. The phylogeny supports previously proposed suprafamilial termite lineages, the Euisoptera and Neoisoptera, a later derived Kalotermitidae as sister group of the Neoisoptera and a monophyletic clade of dampwood (Stolotermitidae, Archotermopsidae) and harvester termites (Hodotermitidae). In contrast to previous termite phylogenetic studies, nodal supports were very high for family-level relationships within termites. Two rare genomic changes in the mt genome control region were found to be molecular synapomorphies for major clades. An elongated stem-loop structure defined the clade Polyphagidae + (Cryptocercus + termites), and a further series of compensatory base changes in this stem loop is synapomorphic for the Neoisoptera. The complicated repeat structures first identified in Reticulitermes, composed of short (A-type) and long (B-type repeats) defines the clade Heterotermitinae + Termitidae, while the secondary loss of A-type repeats is synapomorphic for the non-macrotermitine Termitidae.

Beat phenomena in metal nanowires, and their implications for resonance-based elastic property measurements

The elastic properties of 1D nanostructures such as nanowires are often measured experimentally through actuation of the nanowire at its resonance frequency, and then relating the resonance frequency to the elastic stiffness using elementary beam theory. In the present work, we utilize large scale molecular dynamics simulations to report a novel beat phenomenon in [110]oriented Ag nanowires. The beat phenomenon is found to arise from the asymmetry of the lattice spacing in the orthogonal elementary directions of the [110] nanowire, i.e. the [-110] and [001] directions, which results in two different principal moments of inertia. Because of this, actuations imposed along any other direction are found to decompose into two orthogonal vibrational components based on the actuation angle relative to these two elementary directions, with this phenomenon being generalizable to <110> FCC nanowires of different materials (Cu, Au, Ni, Pd and Pt). The beat phenomenon is explained using a discrete moment of inertia model based on the hard sphere assumption, the model is utilized to show that surface effects enhance the beat phenomenon, while the effect is reduced with increasing nanowires cross-sectional size or aspect ratio. Most importantly, due to the existence of the beat phenomena, we demonstrate that in resonance experiments only a single frequency component is expected to be observed, particularly when the damping ratio is relatively large or very small. Furthermore, for a large range of actuation angles, the lower frequency is more likely to be detected than the higher one, which implies that experimental predictions of Young’s modulus obtained from resonance may in fact be under predictions. The present study therefore has significant implications for experimental interpretations of Young’s modulus as obtained via resonance testing.

Supramolecular selection in molecular alloys

Complexes of the type \[M(phen)3](PF6)2 (M = Ni(II), Fe(II), Ru(II) and phen = 1,10-phenanthroline) were found to co-crystallize to form molecular alloys (solid solutions of molecules) with general formula \[MAxMB1–x(phen)3](PF6)2·0.5H2O in which the relative concentrations of the metal complexes in the crystals closely match those in the crystallizing solution. Consequently, the composition of the co-crystals can be accurately predicted and controlled by modulating the relative concentrations of the metal complexes in the crystallizing solution. Although they are chemically and structurally similar, complexes of the type \[M(bipy)3](PF6)2 (M = Ni(II), Fe(II), Ru(II) and bipy = 2,2′-bipyridine) display markedly different behavior upon co-crystallization. In this case, the resulting co-crystals of general formula \[MAxMB1–x(bipy)3](PF6)2 have relative concentrations of the constituent complexes that are markedly different from the relative concentrations of the complexes initially present in the crystallizing solution. For example, when the nickel and iron complexes are co-crystallized from a solution containing a 50:50 ratio of each, the result is the formation of some crystals with a higher proportion of iron and others with a higher proportion of nickel. The relative concentrations of the metal complexes in the crystals can vary from those in the crystallizing solutions by as much as 15%. This result was observed for a range of combinations of metal complexes (Ni/Fe, Ni/Ru, and Fe/Ru) and a range of starting concentrations in the crystallizing solutions (90:10 through to 10:90 in 10% increments). To explain this remarkable result, we introduce the concept of “supramolecular selection”, which is a process driven by molecular recognition that leads to the partially selective aggregation of like molecules during crystallization.

Surface effects on the dual-mode vibration of <110> silver nanowires with different cross-sections

Dual-mode vibration of nanowires has been reported experimentally through actuation of the nanowire at its resonance frequency, which is expected to open up a variety of new modalities for the NEMS that could operate in the nonlinear regime. In the present work, we utilize large scale molecular dynamics simulations to investigate the dual-mode vibration of <110> Ag nanowires with triangular, rhombic and truncated rhombic cross-sections. By incorporating the generalized Young-Laplace equation into Euler-Bernoulli beam theory, the influence of surface effects on the dual-mode vibration is studied. Due to the different lattice spacing in principal axes of inertia of the {110} atomic layers, the NW is also modeled as a discrete system to reveal the influence from such specific atomic arrangement. It is found that the <110> Ag NW will under a dual-mode vibration if the actuation direction is deviated from the two principal axes of inertia. The predictions of the two first mode natural frequencies by the classical beam model appear underestimated comparing with the MD results, which are found to be enhanced by the discrete model. Particularly, the predictions by the beam theory with the contribution of surface effects are uniformly larger than the classical beam model, which exhibit better agreement with MD results for larger cross-sectional size. However, for ultrathin NWs, current consideration of surface effects is still experiencing certain inaccuracy. In all, for all different cross-sections, the inclusion of surface effects is found to reduce the difference between the two first mode natural frequencies. This trend is observed consistent with MD results. This study provides a first comprehensive investigation on the dual-mode vibration of <110> oriented Ag NWs, which is supposed to benefit the applications of NWs that acting as a resonating beam.

Beyond barcoding : a mitochondrial genomics approach to molecular phylogenetics and diagnostics of blowflies (Diptera: Calliphoridae)

Members of the Calliphoridae (blowflies) are significant for medical and veterinary management, due to the ability of some species to consume living flesh as larvae, and for forensic investigations due to the ability of others to develop in corpses. Due to the difficulty of accurately identifying larval blowflies to species there is a need for DNA-based diagnostics for this family, however the widely used DNA-barcoding marker, cox1, has been shown to fail for several groups within this family. Additionally, many phylogenetic relationships within the Calliphoridae are still unresolved, particularly deeper level relationships. Sequencing whole mt genomes has been demonstrated both as an effective method for identifying the most informative diagnostic markers and for resolving phylogenetic relationships. Twenty-seven complete, or nearly so, mt genomes were sequenced representing 13 species, seven genera and four calliphorid subfamilies and a member of the related family Tachinidae. PCR and sequencing primers developed for sequencing one calliphorid species could be reused to sequence related species within the same superfamily with success rates ranging from 61% to 100%, demonstrating the speed and efficiency with which an mt genome dataset can be assembled. Comparison of molecular divergences for each of the 13 protein-coding genes and 2 ribosomal RNA genes, at a range of taxonomic scales identified novel targets for developing as diagnostic markers which were 117–200% more variable than the markers which have been used previously in calliphorids. Phylogenetic analysis of whole mt genome sequences resulted in much stronger support for family and subfamily-level relationships. The Calliphoridae are polyphyletic, with the Polleninae more closely related to the Tachinidae, and the Sarcophagidae are the sister group of the remaining calliphorids. Within the Calliphoridae, there was strong support for the monophyly of the Chrysomyinae and Luciliinae and for the sister-grouping of Luciliinae with Calliphorinae. Relationships within Chrysomya were not well resolved. Whole mt genome data, supported the previously demonstrated paraphyly of Lucilia cuprina with respect to L. sericata and allowed us to conclude that it is due to hybrid introgression prior to the last common ancestor of modern sericata populations, rather than due to recent hybridisation, nuclear pseudogenes or incomplete lineage sorting.

Analysis of Cone-Beam CT using prior information

Treatment plans for conformal radiotherapy are based on an initial CT scan. The aim is to deliver the prescribed dose to the tumour, while minimising exposure to nearby organs. Recent advances make it possible to also obtain a Cone-Beam CT (CBCT) scan, once the patient has been positioned for treatment. A statistical model will be developed to compare these CBCT scans with the initial CT scan. Changes in the size, shape and position of the tumour and organs will be detected and quantified. Some progress has already been made in segmentation of prostate CBCT scans [1],[2],[3]. However, none of the existing approaches have taken full advantage of the prior information that is available. The planning CT scan is expertly annotated with contours of the tumour and nearby sensitive objects. This data is specific to the individual patient and can be viewed as a snapshot of spatial information at a point in time. There is an abundance of studies in the radiotherapy literature that describe the amount of variation in the relevant organs between treatments. The findings from these studies can form a basis for estimating the degree of uncertainty. All of this information can be incorporated as an informative prior into a Bayesian statistical model. This model will be developed using scans of CT phantoms, which are objects with known geometry. Thus, the accuracy of the model can be evaluated objectively. This will also enable comparison between alternative models.

An investigation into factors affecting electron density calibration for a megavoltage cone-beam CT system

There is a growing interest in the use of megavoltage cone-beam computed tomography (MV CBCT) data for radiotherapy treatment planning. To calculate accurate dose distributions, knowledge of the electron density (ED) of the tissues being irradiated is required. In the case of MV CBCT, it is necessary to determine a calibration-relating CT number to ED, utilizing the photon beam produced for MV CBCT. A number of different parameters can affect this calibration. This study was undertaken on the Siemens MV CBCT system, MVision, to evaluate the effect of the following parameters on the reconstructed CT pixel value to ED calibration: the number of monitor units (MUs) used (5, 8, 15 and 60 MUs), the image reconstruction filter (head and neck, and pelvis), reconstruction matrix size (256 by 256 and 512 by 512), and the addition of extra solid water surrounding the ED phantom. A Gammex electron density CT phantom containing EDs from 0.292 to 1.707 was imaged under each of these conditions. The linear relationship between MV CBCT pixel value and ED was demonstrated for all MU settings and over the range of EDs. Changes in MU number did not dramatically alter the MV CBCT ED calibration. The use of different reconstruction filters was found to affect the MV CBCT ED calibration, as was the addition of solid water surrounding the phantom. Dose distributions from treatment plans calculated with simulated image data from a 15 MU head and neck reconstruction filter MV CBCT image and a MV CBCT ED calibration curve from the image data parameters and a 15 MU pelvis reconstruction filter showed small and clinically insignificant differences. Thus, the use of a single MV CBCT ED calibration curve is unlikely to result in any clinical differences. However, to ensure minimal uncertainties in dose reporting, MV CBCT ED calibration measurements could be carried out using parameter-specific calibration measurements.

Molecular dynamics studies of the effects of branching characteristics on the crystalline structure of polyethylene

Molecular dynamics simulations were carried out on single chain models of linear low-density polyethylene in vacuum to study the effects of branch length, branch content, and branch distribution on the polymer’s crystalline structure at 300 K. The trans/gauche (t/g) ratios of the backbones of the modeled molecules were calculated and utilized to characterize their degree of crystallinity. The results show that the t/g ratio decreases with increasing branch content regardless of branch length and branch distribution, indicating that branch content is the key molecular parameter that controls the degree of crystallinity. Although t/g ratios of the models with the same branch content vary, they are of secondary importance. However, our data suggests that branch distribution (regular or random) has a significant effect on the degree of crystallinity for models containing 10 hexyl branches/1,000 backbone carbons. The fractions of branches that resided in the equilibrium crystalline structures of the models were also calculated. On average, 9.8% and 2.5% of the branches were found in the crystallites of the molecules with ethyl and hexyl branches while C13 NMR experiments showed that the respective probabilities of branch inclusion for ethyl and hexyl branches are 10% and 6% [Hosoda et al., Polymer 1990, 31, 1999–2005]. However, the degree of branch inclusion seems to be insensitive to the branch content and branch distribution.