Molecular Epidemiology of Shiga Toxin-Producing Escherichia coli (STEC) on New Zealand Dairy Farms: Application of a Culture-Independent Assay and Whole-Genome Sequencing

New Zealand has a relatively high incidence of human cases of Shiga toxin-producing Escherichia coli (STEC), with 8.9 STEC cases per 100,000 people reported in 2016. Previous research showed living near cattle and contact with cattle feces as significant risk factors for STEC infections in humans in New Zealand, but infection was not linked to food-associated factors. During the 2014 spring calving season, a random, stratified, cross-sectional study of dairy farms (n = 102) in six regions across New Zealand assessed the prevalence of the "Top 7" STEC bacteria (serogroups O157, O26, O45, O103, O111, O121, and O145) in young calves (n = 1,508), using a culture-independent diagnostic test (PCR/MALDI-TOF). Twenty percent (306/1,508) of calves on 75% (76/102) of dairy farms were positive for at least one of the "Top 7" STEC bacteria. STEC carriage by calves was associated with environmental factors, increased calf age, region, and increased number of calves in a shared calf pen. The intraclass correlation coefficient (ρ) indicated strong clustering of "Top 7" STEC-positive calves for O157, O26, and O45 serogroups within the same pens and farms, indicating that if one calf was positive, others in the same environment were likely to be positive as well. This finding was further evaluated with whole-genome sequencing, which indicated that a single E. coli O26 clonal strain could be found in calves in the same pen or farm, but different strains existed on different farms. This study provides evidence that would be useful for designing on-farm interventions to reduce direct and indirect human exposure to STEC bacteria.IMPORTANCE Cattle are asymptomatic carriers of Shiga toxin-producing E. coli (STEC) bacteria that can cause bloody diarrhea and kidney failure in humans if ingested. New Zealand has relatively high numbers of STEC cases, and contact with cattle feces and living near cattle are risk factors for human infection. This study assessed the national prevalence of STEC in young dairy cattle by randomly selecting 102 farms throughout New Zealand. The study used a molecular laboratory method that has relatively high sensitivity and specificity compared to traditional methods. "Top 7" STEC was found in 20% of calves on 75% of the farms studied, indicating widespread prevalence across the country. By examining the risk factors associated with calf carriage, potential interventions that could decrease the prevalence of "Top 7" STEC bacteria at the farm level were identified, which could benefit both public health and food safety.

Genomic Analysis of Salmonella enterica Serovar Typhimurium DT160 Associated with a 14-Year Outbreak, New Zealand, 1998-2012

During 1998-2012, an extended outbreak of Salmonella enterica serovar Typhimurium definitive type 160 (DT160) affected >3,000 humans and killed wild birds in New Zealand. However, the relationship between DT160 within these 2 host groups and the origin of the outbreak are unknown. Whole-genome sequencing was used to compare 109 Salmonella Typhimurium DT160 isolates from sources throughout New Zealand. We provide evidence that DT160 was introduced into New Zealand around 1997 and rapidly propagated throughout the country, becoming more genetically diverse over time. The genetic heterogeneity was evenly distributed across multiple predicted functional protein groups, and we found no evidence of host group differentiation between isolates collected from human, poultry, bovid, and wild bird sources, indicating ongoing transmission between these host groups. Our findings demonstrate how a comparative genomic approach can be used to gain insight into outbreaks, disease transmission, and the evolution of a multihost pathogen after a probable point-source introduction.

Genomic Epidemiology and Management of Salmonella in Island Ecosystems Used for Takahe Conservation

Translocation and isolation of threatened wildlife in new environments may have unforeseen consequences on pathogen transmission and evolution in host populations. Disease threats associated with intensive conservation management of wildlife remain speculative without gaining an understanding of pathogen dynamics in meta-populations and how location attributes may determine pathogen prevalence. We determined the prevalence and population structure of an opportunistic pathogen, Salmonella, in geographically isolated translocated sub-populations of an endangered New Zealand flightless bird, the takahe (Porphyrio hochstetteri). Out of the nine sub-populations tested, Salmonella was only isolated from takahe living on one private island. The apparent prevalence of Salmonella in takahe on the private island was 32% (95% CI 13-57%), with two serotypes, Salmonella Mississippi and Salmonella houtenae 40:gt-, identified. Epidemiological investigation of reservoirs on the private island and another island occupied by takahe identified environmental and reptile sources of S. Mississippi and S. houtenae 40:gt- on the private island. Single nucleotide polymorphism analysis of core genomes revealed low-level diversity among isolates belonging to the same serotype and little differentiation according to host and environmental source. The pattern observed may be representative of transmission between sympatric hosts and environmental sources, the presence of a common unsampled source, and/or evidence of a recent introduction into the ecosystem. This study highlights how genomic epidemiology can be used to ascertain and understand disease dynamics to inform the management of disease threats in endangered wildlife populations.

Distinct gut microbiome patterns associate with consensus molecular subtypes of colorectal cancer

Colorectal cancer (CRC) is a heterogeneous disease and recent advances in subtype classification have successfully stratified the disease using molecular profiling. The contribution of bacterial species to CRC development is increasingly acknowledged, and here, we sought to analyse CRC microbiomes and relate them to tumour consensus molecular subtypes (CMS), in order to better understand the relationship between bacterial species and the molecular mechanisms associated with CRC subtypes. We classified 34 tumours into CRC subtypes using RNA-sequencing derived gene expression and determined relative abundances of bacterial taxonomic groups using 16S rRNA amplicon metabarcoding. 16S rRNA analysis showed enrichment of Fusobacteria and Bacteroidetes, and decreased levels of Firmicutes and Proteobacteria in CMS1. A more detailed analysis of bacterial taxa using non-human RNA-sequencing reads uncovered distinct bacterial communities associated with each molecular subtype. The most highly enriched species associated with CMS1 included Fusobacterium hwasookii and Porphyromonas gingivalis. CMS2 was enriched for Selenomas and Prevotella species, while CMS3 had few significant associations. Targeted quantitative PCR validated these findings and also showed an enrichment of Fusobacterium nucleatum, Parvimonas micra and Peptostreptococcus stomatis in CMS1. In this study, we have successfully associated individual bacterial species to CRC subtypes for the first time.

Genomic epidemiology of methicillin-susceptible Staphylococcus aureus across colonisation and skin and soft tissue infection

OBJECTIVES

Staphylococcus aureus skin and soft tissue infection (Sa-SSTI) places a significant burden on healthcare systems. New Zealand has a high incidence of Sa-SSTI, and here most morbidity is caused by a polyclonal methicillin-susceptible (MSSA) bacterial population. However, MSSA also colonise asymptomatically the cornified epithelia of approximately 20% of the population, and their divide between commensalism and pathogenicity is poorly understood. We aimed to see whether MSSA are genetically differentiated across colonisation and SSTI; and given the close interactions between people and pets, whether strains isolated from pets differ from human strains.

METHODS

We compared the genomes of contemporaneous colonisation and clinical MSSA isolates obtained in New Zealand from humans and pets.

RESULTS

Core and accessory genome comparisons revealed a homogeneous bacterial population across colonisation, disease, humans, and pets. The rate of MSSA colonisation in dogs was comparatively low (5.4%).

CONCLUSIONS

In New Zealand, most Sa-SSTI morbidity is caused by a random sample of the colonising MSSA population, consistent with the opportunistic infection model rather than the paradigm distinguishing strains according to their pathogenicity. Thus, studies of the factors determining colonisation and immune-escape may be more beneficial than comparative virulence studies. Contact with house-hold pets may pose low zoonotic risk.

Isolation and characterization of microsatellite loci from Arthropodium cirratum (Asparagaceae)

PREMISE OF THE STUDY

Microsatellite markers were developed for Arthropodium cirratum (Asparagaceae) to study population genetic structure and translocation of this species. These markers were tested for cross-amplification in two other Arthropodium species.

METHODS AND RESULTS

Sixteen microsatellite markers were developed from a genomic library and tested in three populations of A. cirratum. The loci exhibited one to five alleles per locus, with private alleles present in each of the populations. Cross-amplification tests in the two other New Zealand Arthropodium species revealed that many of the loci amplify and demonstrate polymorphism in A. bifurcatum.

CONCLUSIONS

These markers will be useful for determining genetic structure in A. cirratum and for determining the origins of translocated populations of this species.

Insights into the Geobacillus stearothermophilus species based on phylogenomic principles

BACKGROUND

The genus Geobacillus comprises bacteria that are Gram positive, thermophilic spore-formers, which are found in a variety of environments from hot-springs, cool soils, to food manufacturing plants, including dairy manufacturing plants. Despite considerable interest in the use of Geobacillus spp. for biotechnological applications, the taxonomy of this genus is unclear, in part because of differences in DNA-DNA hybridization (DDH) similarity values between studies. In addition, it is also difficult to use phenotypic characteristics to define a bacterial species. For example, G. stearothermophilus was traditionally defined as a species that does not utilise lactose, but the ability of dairy strains of G. stearothermophilus to use lactose has now been well established.

RESULTS

This study compared the genome sequences of 63 Geobacillus isolates and showed that based on two different genomic approaches (core genome comparisons and average nucleotide identity) the Geobacillus genus could be divided into sixteen taxa for those Geobacillus strains that have genome sequences available thus far. In addition, using Geobacillus stearothermophilus as an example, we show that inclusion of the accessory genome, as well as phenotypic characteristics, is not suitable for defining this species. For example, this is the first study to provide evidence of dairy adaptation in G. stearothermophilus - a phenotypic feature not typically considered standard in this species - by identifying the presence of a putative lac operon in four dairy strains.

CONCLUSIONS

The traditional polyphasic approach of combining both genotypic and phenotypic characteristics to define a bacterial species could not be used for G. stearothermophilus where many phenotypic characteristics vary within this taxon. Further evidence of this discordant use of phenotypic traits was provided by analysis of the accessory genome, where the dairy strains contained a putative lac operon. Based on the findings from this study, we recommend that novel bacterial species should be defined using a core genome approach.

Culture independent analysis using gnd as a target gene to assess Escherichia coli diversity and community structure

Current culture methods to investigate changes in Escherichia coli community structure are often slow and laborious. Genes such as gnd (6-phosphogluconate dehydrogenase) have a highly variable nucleotide sequence and may provide a target for E. coli microbiome analysis using culture-independent methods. Metabarcoded PCR primers were used to generate separate libraries from calf faecal samples for high throughput sequencing. Although a total of 348 separate gnd sequence types (gSTs) were identified, 188 were likely to be due to sequencing errors. Of the remaining 160 gSTs, 92 did not match those in a database of 319 separate gnd sequences. ‘Animal’ was the main determinant of E. coli diversity with limited impact of sample type or DNA extraction method on intra-host E. coli community variation from faeces and recto-anal mucosal swab samples. This culture-independent study has addressed the difficulties of quantifying bacterial intra-species diversity and revealed that, whilst individual animals may harbour >50 separate E. coli strains, communities are dominated by <10 strains alongside a large pool of subdominant strains present at low abundances. This method will be useful for characterising the diversity and population structure of E. coli in experimental studies designed to assess the impact of interventions on the gut microbiome.

Inferring Ancestral Recombination Graphs from Bacterial Genomic Data

Homologous recombination is a central feature of bacterial evolution, yet it confounds traditional phylogenetic methods. While a number of methods specific to bacterial evolution have been developed, none of these permit joint inference of a bacterial recombination graph and associated parameters. In this article, we present a new method which addresses this shortcoming. Our method uses a novel Markov chain Monte Carlo algorithm to perform phylogenetic inference under the ClonalOrigin model. We demonstrate the utility of our method by applying it to ribosomal multilocus sequence typing data sequenced from pathogenic and nonpathogenic Escherichia coli serotype O157 and O26 isolates collected in rural New Zealand. The method is implemented as an open source BEAST 2 package, Bacter, which is available via the project web page at http://tgvaughan.github.io/bacter.

Sticky Genomes: Using NGS Evidence to Test Hybrid Speciation Hypotheses

Hypotheses of hybrid origin are common. Here we use next generation sequencing to test a hybrid hypothesis for a non-model insect with a large genome. We compared a putative hybrid triploid stick insect species (Acanthoxyla geisovii) with its putative paternal diploid taxon (Clitarchus hookeri), a relationship that provides clear predictions for the relative genetic diversity within each genome. The parental taxon is expected to have comparatively low allelic diversity that is nested within the diversity of the hybrid daughter genome. The scale of genome sequencing required was conveniently achieved by extracting mRNA and sequencing cDNA to examine expressed allelic diversity. This allowed us to test hybrid-progenitor relationships among non-model organisms with large genomes and different ploidy levels. Examination of thousands of independent loci avoids potential problems produced by the silencing of parts of one or other of the parental genomes, a phenomenon sometimes associated with the process of stabilisation of a hybrid genome. Transcript assembles were assessed for evidence of paralogs and/or alternative splice variants before proceeding. Comparison of transcript assemblies was not an appropriate measure of genetic variability, but by mapping reads back to clusters derived from each species we determined levels of allelic diversity. We found greater cDNA sequence diversity among alleles in the putative hybrid species (Acanthoxyla geisovii) than the non-hybrid. The allelic diversity within the putative paternal species (Clitachus hookeri) nested within the hybrid-daughter genome, supports the current view of a hybrid-progenitor relationship for these stick insect species. Next generation sequencing technology provides opportunities for testing evolutionary hypotheses with non-model organisms, including, as here, genomes that are large due to polyploidy.

Microbial Genomics of a Host-Associated Commensal Bacterium in Fragmented Populations of Endangered Takahe

Isolation of wildlife into fragmented populations as a consequence of anthropogenic-mediated environmental change may alter host-pathogen relationships. Our understanding of some of the epidemiological features of infectious disease in vulnerable populations can be enhanced by the use of commensal bacteria as a proxy for invasive pathogens in natural ecosystems. The distinctive population structure of a well-described meta-population of a New Zealand endangered flightless bird, the takahe (Porphyrio hochstetteri), provided a unique opportunity to investigate the influence of host isolation on enteric microbial diversity. The genomic epidemiology of a prevalent rail-associated endemic commensal bacterium was explored using core genome and ribosomal multilocus sequence typing (rMLST) of 70 Campylobacter sp. nova 1 isolated from one third of the takahe population resident in multiple locations. While there was evidence of recombination between lineages, bacterial divergence appears to have occurred and multivariate analysis of 52 rMLST genes revealed location-associated differentiation of C. sp. nova 1 sequence types. Our results indicate that fragmentation and anthropogenic manipulation of populations can influence host-microbial relationships, with potential implications for niche adaptation and the evolution of micro-organisms in remote environments. This study provides a novel framework in which to explore the complex genomic epidemiology of micro-organisms in wildlife populations.

Using a common commensal bacterium in endangered Takahe as a model to explore pathogen dynamics in isolated wildlife populations

Predicting and preventing outbreaks of infectious disease in endangered wildlife is problematic without an understanding of the biotic and abiotic factors that influence pathogen transmission and the genetic variation of microorganisms within and between these highly modified host communities. We used a common commensal bacterium, Campylobacter spp., in endangered Takahe (Porphyrio hochstetteri) populations to develop a model with which to study pathogen dynamics in isolated wildlife populations connected through ongoing translocations. Takahe are endemic to New Zealand, where their total population is approximately 230 individuals. Takahe were translocated from a single remnant wild population to multiple offshore and mainland reserves. Several fragmented subpopulations are maintained and connected through regular translocations. We tested 118 Takahe from 8 locations for fecal Campylobacter spp. via culture and DNA extraction and used PCR for species assignment. Factors relating to population connectivity and host life history were explored using multivariate analytical methods to determine associations between host variables and bacterial prevalence. The apparent prevalence of Campylobacter spp. in Takahe was 99%, one of the highest reported in avian populations. Variation in prevalence was evident among Campylobacter species identified. C. sp. nova 1 (90%) colonized the majority of Takahe tested. Prevalence of C. jejuni (38%) and C. coli (24%) was different between Takahe subpopulations, and this difference was associated with factors related to population management, captivity, rearing environment, and the presence of agricultural practices in the location in which birds were sampled. Modeling results of Campylobacter spp. in Takahe metapopulations suggest that anthropogenic management of endangered species within altered environments may have unforeseen effects on microbial exposure, carriage, and disease risk. Translocation of wildlife between locations could have unpredictable consequences including the spread of novel microbes between isolated populations.

Biochar in co-contaminated soil manipulates arsenic solubility and microbiological community structure, and promotes organochlorine degradation

We examined the effect of biochar on the water-soluble arsenic (As) concentration and the extent of organochlorine degradation in a co-contaminated historic sheep-dip soil during a 180-d glasshouse incubation experiment. Soil microbial activity, bacterial community and structure diversity were also investigated. Biochar made from willow feedstock (Salix sp) was pyrolysed at 350 or 550°C and added to soil at rates of 10 g kg^-1 and 20 g kg^-1 (representing 30 t ha^-1 and 60 t ha^-1). The isomers of hexachlorocyclohexane (HCH) alpha-HCH and gamma-HCH (lindane), underwent 10-fold and 4-fold reductions in concentration as a function of biochar treatment. Biochar also resulted in a significant reduction in soil DDT levels (P < 0.01), and increased the DDE:DDT ratio. Soil microbial activity was significantly increased (P < 0.01) under all biochar treatments after 60 days of treatment compared to the control. 16S amplicon sequencing revealed that biochar-amended soil contained more members of the Chryseobacterium, Flavobacterium, Dyadobacter and Pseudomonadaceae which are known bioremediators of hydrocarbons. We hypothesise that a recorded short-term reduction in the soluble As concentration due to biochar amendment allowed native soil microbial communities to overcome As-related stress. We propose that increased microbiological activity (dehydrogenase activity) due to biochar amendment was responsible for enhanced degradation of organochlorines in the soil. Biochar therefore partially overcame the co-contaminant effect of As, allowing for enhanced natural attenuation of organochlorines in soil.

Robust identification of noncoding RNA from transcriptomes requires phylogenetically-informed sampling

Noncoding RNAs are integral to a wide range of biological processes, including translation, gene regulation, host-pathogen interactions and environmental sensing. While genomics is now a mature field, our capacity to identify noncoding RNA elements in bacterial and archaeal genomes is hampered by the difficulty of de novo identification. The emergence of new technologies for characterizing transcriptome outputs, notably RNA-seq, are improving noncoding RNA identification and expression quantification. However, a major challenge is to robustly distinguish functional outputs from transcriptional noise. To establish whether annotation of existing transcriptome data has effectively captured all functional outputs, we analysed over 400 publicly available RNA-seq datasets spanning 37 different Archaea and Bacteria. Using comparative tools, we identify close to a thousand highly-expressed candidate noncoding RNAs. However, our analyses reveal that capacity to identify noncoding RNA outputs is strongly dependent on phylogenetic sampling. Surprisingly, and in stark contrast to protein-coding genes, the phylogenetic window for effective use of comparative methods is perversely narrow: aggregating public datasets only produced one phylogenetic cluster where these tools could be used to robustly separate unannotated noncoding RNAs from a null hypothesis of transcriptional noise. Our results show that for the full potential of transcriptomics data to be realized, a change in experimental design is paramount: effective transcriptomics requires phylogeny-aware sampling.

We have analysed more than 400 public transcriptomes, generated using RNA-seq, from almost 40 strains of Bacteria and Archaea. We discovered that the capacity to identify noncoding RNA outputs from this data is strongly dependent on phylogenetic sampling. Our results show that, for the full potential of transcriptomics data as a discovery tool to be realized, a change in experimental design is critical: effective comparative transcriptomics requires phylogeny-aware sampling. We also examined how comparative transcriptomics experiments can be used to effectively identify RNA elements. We find that, for RNA element discovery, a phylogeny-informed sampling approach is more effective than analyses of individual species. Phylogeny-informed sampling reveals a narrow ‘Goldilocks Zone’ (where species are not too similar and not too divergent) for RNA identification using clusters of related species. In stark contrast to protein-coding genes, not only is the phylogenetic window for the effective use of comparative methods for noncoding RNA identification perversely narrow, but few existing datasets sit within this Goldilocks Zone: by aggregating public datasets, we were only able to create one phylogenetic cluster where comparative tools could be used to confidently separate unannotated noncoding RNAs from transcriptional noise.

Quantification of historical livestock importation into New Zealand 1860-1979

AIMS

To quantify the numbers of live cattle, sheep and poultry imported into New Zealand and, where possible, their country of origin from 1860 to 1979.

METHODS

Information on the origin and number of live animal importations into New Zealand was collected for cattle, sheep and poultry for the period 1868-1979 from the annual reports compiled by the New Zealand Registrar General's Office, Government Statistician's Office, Census and Statistics Office, Census and Statistics Department, Customs Department and Department of Statistics. Census data from 1851 to 1871 were also used to estimate the livestock population during this period. The number of animals imported and the mean population for each species in a decade were determined, and the major countries of origin were identified.

RESULTS

A large number of cattle (53,384) and sheep (604,525) were imported in the 1860s, and then there was a marked reduction in importations. Live poultry were imported in relatively small numbers (20,701) from 1880 to 1939, then 1,564,330 live poultry were imported between 1960 and 1979. Australia was the predominant country of origin for sheep between 1868 and 1959 (51,347/60,918; 84.3%) and of cattle between 1868 and 1979 (10,080/15,157; 66.5%). Only 6,712 (11.0%) sheep and 3,909 (25.8%) cattle were imported from the United Kingdom over the same periods, and even fewer from other countries.

CONCLUSIONS

The collated data and historical reports show that from 1860 to 1979 Australia has been the main source of livestock introduced into New Zealand. The pattern of importation showed that large numbers of cattle and sheep were initially imported in the 1860s, probably in response to rapid agricultural expansion. Thereafter importations continued at much reduced numbers. In contrast, relatively small numbers of poultry were introduced until the 1960s when large numbers were imported as part of the development of a modern high-production industry. The overall pattern for both cattle and sheep was of a bottleneck event, as initially a relatively limited number of animals arrived from outside populations, followed by population expansion with ongoing but limited immigration (admixture). Investigation into the genetic population structure of New Zealand's cattle and sheep, as well as their host-associated microorganisms, could reflect the impact of these early historical events.

Mutational dynamics of aroid chloroplast genomes

A characteristic feature of eukaryote and prokaryote genomes is the co-occurrence of nucleotide substitution and insertion/deletion (indel) mutations. Although similar observations have also been made for chloroplast DNA, genome-wide associations have not been reported. We determined the chloroplast genome sequences for two morphotypes of taro (Colocasia esculenta; family Araceae) and compared these with four publicly available aroid chloroplast genomes. Here, we report the extent of genome-wide association between direct and inverted repeats, indels, and substitutions in these aroid chloroplast genomes. We suggest that alternative but not mutually exclusive hypotheses explain the mutational dynamics of chloroplast genome evolution.

De novo sequence assembly and characterisation of a partial transcriptome for an evolutionarily distinct reptile, the tuatara (Sphenodon punctatus)

Background

The tuatara (Sphenodon punctatus) is a species of extraordinary zoological interest, being the only surviving member of an entire order of reptiles which diverged early in amniote evolution. In addition to their unique phylogenetic placement, many aspects of tuatara biology, including temperature-dependent sex determination, cold adaptation and extreme longevity have the potential to inform studies of genome evolution and development. Despite increasing interest in the tuatara genome, genomic resources for the species are still very limited. We aimed to address this by assembling a transcriptome for tuatara from an early-stage embryo, which will provide a resource for genome annotation, molecular marker development and studies of development and adaptation in tuatara.

Results

We obtained 30 million paired-end 50 bp reads from an Illumina Genome Analyzer and assembled them with Velvet and Oases using a range of kmers. After removing redundancy and filtering out low quality transcripts, our transcriptome dataset contained 32911 transcripts, with an N50 of 675 and a mean length of 451 bp. Almost 50% (15965) of these transcripts could be annotated by comparison with the NCBI non-redundant (NR) protein database or the chicken, green anole and zebrafish UniGene sets. A scan of candidate genes and repetitive elements revealed genes involved in immune function, sex differentiation and temperature-sensitivity, as well as over 200 microsatellite markers.

Conclusions

This dataset represents a major increase in genomic resources for the tuatara, increasing the number of annotated gene sequences from just 60 to almost 16,000. This will facilitate future research in sex determination, genome evolution, local adaptation and population genetics of tuatara, as well as inform studies on amniote evolution.

The evolutionary root of flowering plants

Correct rooting of the angiosperm radiation is both challenging and necessary for understanding the origins and evolution of physiological and phenotypic traits in flowering plants. The problem is known to be difficult due to the large genetic distance separating flowering plants from other seed plants and the sparse taxon sampling among basal angiosperms. Here, we provide further evidence for concern over substitution model misspecification in analyses of chloroplast DNA sequences. We show that support for Amborella as the sole representative of the most basal angiosperm lineage is founded on sequence site patterns poorly described by time-reversible substitution models. Improving the fit between sequence data and substitution model identifies Trithuria, Nymphaeaceae, and Amborella as surviving relatives of the most basal lineage of flowering plants. This finding indicates that aquatic and herbaceous species dominate the earliest extant lineage of flowering plants. [; ; ; ; ; .].

Identification of a novel nidovirus associated with a neurological disease of the Australian brushtail possum (Trichosurus vulpecula)

A novel, fatal neurological disease of the Australian brushtail possum (Trichosurus vulpecula) was first identified in 1995 in a research facility and subsequently in free-living possums in New Zealand and termed wobbly possum disease (WPD). The results of previous transmission studies suggested that the aetiological agent of WPD is most likely a virus. However, the identity of the presumed viral agent had not been elucidated. In the current report, we describe identification of a novel virus from tissues of WPD-affected possums using a combination of next generation sequencing and traditional molecular methods. The proportion of possums positive for the novel virus by PCR was significantly higher (p<0.0001) among animals with WPD than clinically healthy possums, strongly suggesting an aetiological involvement of the virus in WPD. Analysis of the partial genomic sequence of the putative WPD virus indicated that it is a novel nidovirus, most closely related to the current members of the family Arteriviridae.

Whole-genome comparison of two Campylobacter jejuni isolates of the same sequence type reveals multiple loci of different ancestral lineage

Campylobacter jejuni ST-474 is the most important human enteric pathogen in New Zealand, and yet this genotype is rarely found elsewhere in the world. Insight into the evolution of this organism was gained by a whole genome comparison of two ST-474, flaA SVR-14 isolates and other available C. jejuni isolates and genomes. The two isolates were collected from different sources, human (H22082) and retail poultry (P110b), at the same time and from the same geographical location. Solexa sequencing of each isolate resulted in 1.659 Mb (H22082) and 1.656 Mb (P110b) of assembled sequences within 28 (H22082) and 29 (P110b) contigs. We analysed 1502 genes for which we had sequences within both ST-474 isolates and within at least one of 11 C. jejuni reference genomes. Although 94.5% of genes were identical between the two ST-474 isolates, we identified 83 genes that differed by at least one nucleotide, including 55 genes with non-synonymous substitutions. These covered 101 kb and contained 672 point differences. We inferred that 22 (3.3%) of these differences were due to mutation and 650 (96.7%) were imported via recombination. Our analysis estimated 38 recombinant breakpoints within these 83 genes, which correspond to recombination events affecting at least 19 loci regions and gives a tract length estimate of 2 kb. This includes a 12 kb region displaying non-homologous recombination in one of the ST-474 genomes, with the insertion of two genes, including ykgC, a putative oxidoreductase, and a conserved hypothetical protein of unknown function. Furthermore, our analysis indicates that the source of this recombined DNA is more likely to have come from C. jejuni strains that are more closely related to ST-474. This suggests that the rates of recombination and mutation are similar in order of magnitude, but that recombination has been much more important for generating divergence between the two ST-474 isolates.

Systematic error in seed plant phylogenomics

Resolving the closest relatives of Gnetales has been an enigmatic problem in seed plant phylogeny. The problem is known to be difficult because of the extent of divergence between this diverse group of gymnosperms and their closest phylogenetic relatives. Here, we investigate the evolutionary properties of conifer chloroplast DNA sequences. To improve taxon sampling of Cupressophyta (non-Pinaceae conifers), we report sequences from three new chloroplast (cp) genomes of Southern Hemisphere conifers. We have applied a site pattern sorting criterion to study compositional heterogeneity, heterotachy, and the fit of conifer chloroplast genome sequences to a general time reversible + G substitution model. We show that non-time reversible properties of aligned sequence positions in the chloroplast genomes of Gnetales mislead phylogenetic reconstruction of these seed plants. When 2,250 of the most varied sites in our concatenated alignment are excluded, phylogenetic analyses favor a close evolutionary relationship between the Gnetales and Pinaceae—the Gnepine hypothesis. Our analytical protocol provides a useful approach for evaluating the robustness of phylogenomic inferences. Our findings highlight the importance of goodness of fit between substitution model and data for understanding seed plant phylogeny.

RNA networks in prokaryotes I: CRISPRs and riboswitches

As with eukaryotes, prokaryotes employ a variety of mechanisms to allow the various types of RNA to interact and perform complex functions as a network. This chapter will detail prokaryotic molecular systems, such as riboswitches and CRISPRs, to show how they perform unique functions within the cell. These systems can interact with each other to gain a higher level of control and here we highlight some examples of such interactions including the cleavage of certain riboswitches by RNaseP, and endoribonuclease cleavage of pre-crRNAs in the CRISPR system. Thanks to such insights, we are beginning to get a glimpse of the prokaryotic RNA infrastructure, just as we have done with eukaryotes.

Index-free de novo assembly and deconvolution of mixed mitochondrial genomes

Second-generation sequencing technology has allowed a very large increase in sequencing throughput. In order to make use of this high throughput, we have developed a pipeline for sequencing and de novo assembly of multiple mitochondrial genomes without the costs of indexing. Simulation studies on a mixture of diverse animal mitochondrial genomes showed that mitochondrial genomes could be reassembled from a high coverage of short (35 nt) reads, such as those generated by a second-generation Illumina Genome Analyzer. We then assessed this experimentally with long-range polymerase chain reaction products from mitochondria of a human, a rat, a bird, a frog, an insect, and a mollusc. Comparison with reference genomes was used for deconvolution of the assembled contigs rather than for mapping of sequence reads. As proof of concept, we report the complete mollusc mitochondrial genome of an olive shell (Amalda northlandica). It has a very unusual putative control region, which contains a structure that would probably only be detectable by next-generation sequencing. The general approach has considerable potential, especially when combined with indexed sequencing of different groups of genomes.