Insights from entire mitochondrial genome sequences into the phylogeny of ticks of the genera Haemaphysalis and Archaeocroton with the elevation of the subgenus Alloceraea Schulze, 1919 back to the status of a genus.

We used entire mitochondrial (mt) genome sequences (14.5-15 kbp) to resolve the phylogeny of the four main lineages of the Haematobothrion ticks: Alloceraea, Archaeocroton, Bothriocroton and Haemaphysalis. In our phylogenetic trees, Alloceraea was the sister to Archaeocroton sphenodonti, a tick of an archetypal reptile, the tuatara, from New Zealand, to the exclusion of the rest of the species of Haemaphysalis. The mt genomes of all four of the Alloceraea species that have been sequenced so far had a substantial insert, 132-312 bp, between the tRNA-Glu (E) gene and the nad1 gene in their mt genomes. This insert was not found in any of the other eight subgenera of Haemaphysalis. The mt genomes of 13 species of Haemaphysalis from NCBI GenBank were added to the most recent data set on Haemaphysalis and its close relatives to help resolve the phylogeny of Haemaphysalis, including five new subgenera of Haemaphysalis not previously considered by other authors: Allophysalis (structurally primitive), Aboimisalis (structurally primitive), Herpetobia (structurally intermediate), Ornithophysalis (structurally advanced) and Segalia (structurally advanced). We elevated Alloceraea Schulze, 1919 to the status of genus because Alloceraea Schulze, 1919 is phylogenetically distinct from the other subgenera of Haemaphysalis. Moreover, we propose that the subgenus Allophysalis is the sister to the rest of the Haemaphysalis (14 subgenera) and that the 'structurally primitive' subgenera Hoogstraal and Kim comprise early diverging lineages. Our matrices of the pairwise genetic difference (percent) of mt genomes and partial 16S rRNA sequences indicated that the mt genome sequence of Al. kitaokai (gb# OM368280) may not be Al. kitaokai Hoogstraal, 1969 but rather another species of Alloceraea. In a similar way, the mt genome sequence of H. (Herpetobia) nepalensis Hoogstraal, 1962 (gb# NC_064124) was only 2% genetically different to that of H. (Allophysalis) tibetensis Hoogstraal, 1965 (gb# OM368293): this indicates to us that they are the same species. Alloceraea cretacea may be better placed in a genus other than Alloceraea Schulze, 1919. Reptiles may have been the host to the most recent common ancestor of Archaeocroton and Alloceraea.

A high fidelity approach to assembling the complex Borrelia genome.

BACKGROUND: Bacteria of the Borrelia burgdorferi sensu lato (s.l.) complex can cause Lyme borreliosis. Different B. burgdorferi s.l. genospecies vary in their host and vector associations and human pathogenicity but the genetic basis for these adaptations is unresolved and requires completed and reliable genomes for comparative analyses. The de novo assembly of a complete Borrelia genome is challenging due to the high levels of complexity, represented by a high number of circular and linear plasmids that are dynamic, showing mosaic structure and sequence homology. Previous work demonstrated that even advanced approaches, such as a combination of short-read and long-read data, might lead to incomplete plasmid reconstruction. Here, using recently developed high-fidelity (HiFi) PacBio sequencing, we explored strategies to obtain gap-free, complete and high quality Borrelia genome assemblies. Optimizing genome assembly, quality control and refinement steps, we critically appraised existing techniques to create a workflow that lead to improved genome reconstruction. RESULTS: Despite the latest available technologies, stand-alone sequencing and assembly methods are insufficient for the generation of complete and high quality Borrelia genome assemblies. We developed a workflow pipeline for the de novo genome assembly for Borrelia using several types of sequence data and incorporating multiple assemblers to recover the complete genome including both circular and linear plasmid sequences. CONCLUSION: Our study demonstrates that, with HiFi data and an ensemble reconstruction pipeline with refinement steps, chromosomal and plasmid sequences can be fully resolved, even for complex genomes such as Borrelia. The presented pipeline may be of interest for the assembly of further complex microbial genomes.

An invasive human commensal and a native marsupial maintain tick populations at the urban fringe.

Ticks (Acari: Ixodidae) are major disease vectors globally making it increasingly important to understand how altered vertebrate communities in urban areas shape tick population dynamics. In urban landscapes of Australia, little is known about which native and introduced small mammals maintain tick populations preventing host-targeted tick management and leading to human-wildlife conflict. Here, we determined (1) larval, nymphal, and adult tick burdens on host species and potential drivers, (2) the number of ticks supported by the different host populations, and (3) the proportion of medically significant tick species feeding on the different host species in Northern Sydney. We counted 3551 ticks on 241 mammals at 15 sites and found that long-nosed bandicoots (Perameles nasuta) hosted more ticks of all life stages than other small mammals but introduced black rats (Rattus rattus) were more abundant at most sites (33%-100%) and therefore important in supporting larval and nymphal ticks in our study areas. Black rats and bandicoots hosted a greater proportion of medically significant tick species including Ixodes holocyclus than other hosts. Our results show that an introduced human commensal contributes to maintaining urban tick populations and suggests ticks could be managed by controlling rat populations on urban fringes.

Rejection of the name Borreliella and all proposed species comb. nov. placed therein.

Rejection (nomen rejiciendum) of the name Borreliella and all new combinations therein is being requested on grounds of risk to human health and patient safety (Principle 1, subprinciple 2 and Rule 56a) and violation to aim for stability of names, to avoid useless creation of names (Principle 1, subprinciple 1 and 3) and that names should not be changed without sufficient reason (Principle 9 of the International Code of Nomenclature of Prokaryotes).

Phylogenetic characterisation of two novel Anaplasmataceae from Australian Ixodes holocyclus ticks: 'Candidatus Neoehrlichia australis' and 'Candidatus Neoehrlichia arcana'.

Recently, two novel species of Anaplasmataceae were detected in the Australian paralysis tick, Ixodes holocyclus, by 16S rRNA gene metabarcoding. Analysis of these sequences suggested that these novel organisms are closely related to the genus 'Candidatus Neoehrlichia'. In this study, phylogenetic analysis of 16S rRNA (1264 bp), groESL (1047 bp) and gltA (561 bp) gene sequences, and concatenated (2872 bp) sequences, all concur that these novel species belong in the genus 'Candidatus Neoehrlichia' and are most closely related to, but distinct from the only other recognised members of this genus, 'Candidatus Neoehrlichia mikurensis' and 'Candidatus Neoehrlichia lotoris'. Based on their unique molecular signature, we propose to designate these species 'Candidatus Neoehrlichia australis' (reference strain HT41R) and 'Candidatus Neoehrlichia arcana' (reference strain HT94R). Identical 'Candidatus Neoehrlichia australis' 16S rRNA, groESL and gltA sequences were detected in 34/391 (8.7 %) individual Ixodes holocyclus ticks, and sequences were most similar to 'Candidatus Neoehrlichia lotoris' (96.2 %, 83.1 % and 67.2 %, respectively) and 'Candidatus Neoehrlichia mikurensis' (96.2 %, 84 % and 68.4 % respectively). Likewise, identical 'Candidatus Neoehrlichia arcana' 16S rRNA, groESL and gltA sequences were detected in 12/391 (3.1 %) Ixodes holocyclus ticks, and sequences were most similar to 'Candidatus Neoehrlichia lotoris' (98.5 %, 88.7 % and 79.3 %, respectively) and 'Candidatus Neoehrlichia mikurensis' (96.3 %, 84 % and 67.4 % respectively). These new species are the first Anaplasmataceae (except Wolbachia spp.) to be found to be endemic to Australia. The pathogenic consequences of these organisms are yet to be determined.

Dietary cholesterol modulates pathogen blocking by Wolbachia.

The bacterial endosymbiont Wolbachia pipientis protects its hosts from a range of pathogens by limiting their ability to form infections inside the insect. This "pathogen blocking" could be explained by innate immune priming by the symbiont, competition for host-derived resources between pathogens and Wolbachia, or the direct modification of the cell or cellular environment by Wolbachia. Recent comparative work in Drosophila and the mosquito Aedes aegypti has shown that an immune response is not required for pathogen blocking, implying that there must be an additional component to the mechanism. Here we have examined the involvement of cholesterol in pathogen blocking using a system of dietary manipulation in Drosophila melanogaster in combination with challenge by Drosophila C virus (DCV), a common fly pathogen. We observed that flies reared on cholesterol-enriched diets infected with the Wolbachia strains wMelPop and wMelCS exhibited reduced pathogen blocking, with viral-induced mortality occurring 2-5 days earlier than flies reared on Standard diet. This shift toward greater virulence in the presence of cholesterol also corresponded to higher viral copy numbers in the host. Interestingly, an increase in dietary cholesterol did not have an effect on Wolbachia density except in one case, but this did not directly affect the strength of pathogen blocking. Our results indicate that host cholesterol levels are involved with the ability of Wolbachia-infected flies to resist DCV infections, suggesting that cholesterol contributes to the underlying mechanism of pathogen blocking.

Comparison of Sanger and next generation sequencing performance for genotyping Cryptosporidium isolates at the 18S rRNA and actin loci.

Cryptosporidium is an important enteric pathogen that infects a wide range of humans and animals. Rapid and reliable detection and characterisation methods are essential for understanding the transmission dynamics of the parasite. Sanger sequencing, and high-throughput sequencing (HTS) on an Ion Torrent platform, were compared with each other for their sensitivity and accuracy in detecting and characterising 25 Cryptosporidium-positive human and animal faecal samples. Ion Torrent reads (n = 123,857) were obtained at both 18S rRNA and actin loci for 21 of the 25 samples. Of these, one isolate at the actin locus (Cattle 05) and three at the 18S rRNA locus (HTS 10, HTS 11 and HTS 12), suffered PCR drop-out (i.e. PCR failures) when using fusion-tagged PCR. Sanger sequences were obtained for both loci for 23 of the 25 samples and showed good agreement with Ion Torrent-based genotyping. Two samples both from pythons (SK 02 and SK 05) produced mixed 18S and actin chromatograms by Sanger sequencing but were clearly identified by Ion Torrent sequencing as C. muris. One isolate (SK 03) was typed as C. muris by Sanger sequencing but was identified as a mixed C. muris and C. tyzzeri infection by HTS. 18S rRNA Type B sequences were identified in 4/6 C. parvum isolates when deep sequenced but were undetected in Sanger sequencing. Sanger was cheaper than Ion Torrent when sequencing a small numbers of samples, but when larger numbers of samples are considered (n = 60), the costs were comparative. Fusion-tagged amplicon based approaches are a powerful way of approaching mixtures, the only draw-back being the loss of PCR efficiency on low-template samples when using primers coupled to MID tags and adaptors. Taken together these data show that HTS has excellent potential for revealing the "true" composition of species/types in a Cryptosporidium infection, but that HTS workflows need to be carefully developed to ensure sensitivity, accuracy and contamination are controlled.

The first finding of Dictyocaulus cervi and Dictyocaulus skrjabini (Nematoda) in feral fallow deer (Dama dama) in Australia.

Feral deer are widespread throughout Australia with the capacity to impact livestock production via transmission of parasites. Samples of Dama dama (fallow deer), Rusa unicolor (sambar deer), Cervus elaphus (red deer) and an unidentified deer were sourced from various locations in south-eastern Australia for examination for parasites. Adult nematodes were collected from the lungs of all deer species across four separate geographical locations. The nematodes were identified as species of Dictyocaulus through both morphological and molecular means. Species identification based on morphological features was difficult, with many measurements from described species overlapping. Molecular analyses targeting three markers, namely 18S rRNA, ITS2, and cox1 revealed the presence of two distinct species: Dictyocaulus cervi and Dictyocaulus skrjabini. These are the first genetically confirmed reports of species of Dictyocaulus in feral deer in Australia, and although cross-transmission of species of Dictyocaulus with livestock has not yet been reported, it cannot be completely discounted without further research.

Unraveling the phylogenetics of genetically closely related species, Haemaphysalis japonica and Haemaphysalis megaspinosa, using entire tick mitogenomes and microbiomes.

Ticks have a profound impact on public health. Haemaphysalis is one of the most widespread genera in Asia, including Japan. The taxonomy and genetic differentiation of Haemaphysalis spp. is challenging. For instance, previous studies struggled to distinguish Haemaphysalis japonica and Haemaphysalis megaspinosa due to the dearth of nucleotide sequence polymorphisms in widely used barcoding genes. The classification of H. japonica japonica and its related sub-species Haemaphysalis japonica douglasi or Haemaphysalis jezoensis is also confused due to their high morphological similarity and a lack of molecular data that support the current classification. We used mitogenomes and microbiomes of H. japonica and H. megaspinosa to gain deeper insights into the phylogenetic relationships and genetic divergence between two species. Phylogenetic analyses of concatenated nucleotide sequences of protein-coding genes and ribosomal DNA genes distinguished H. japonica and H. megaspinosa as monophyletic clades, with further subdivision within the H. japonica clade. The 16S rRNA and NAD5 genes were valuable markers for distinguishing H. japonica and H. megaspinosa. Population genetic structure analyses indicated that genetic variation within populations accounted for a large proportion of the total variation compared to variation between populations. Microbiome analyses revealed differences in alpha and beta diversity between H. japonica and H. megaspinosa: H. japonica had the higher diversity. Coxiella sp., a likely endosymbiont, was found in both Haemaphysalis species. The abundance profiles of likely endosymbionts, pathogens, and commensals differed between H. japonica and H. megaspinosa: H. megaspinosa was more diverse.

The first cryptic genus of Ixodida, Cryptocroton n. gen. for Amblyomma papuanum Hirst, 1914: a tick of North Queensland, Australia, and Papua New Guinea.

We describe a new genus Cryptocroton n. gen. for Amblyomma papuanum Hirst, 1914, a tick of North Queensland, Australia, and Papua New Guinea.

New insights into the molecular phylogeny, biogeographical history, and diversification of Amblyomma ticks (Acari: Ixodidae) based on mitogenomes and nuclear sequences.

BACKGROUND: Amblyomma is the third most diversified genus of Ixodidae that is distributed across the Indomalayan, Afrotropical, Australasian (IAA), Nearctic and Neotropical biogeographic ecoregions, reaching in the Neotropic its highest diversity. There have been hints in previously published phylogenetic trees from mitochondrial genome, nuclear rRNA, from combinations of both and morphology that the Australasian Amblyomma or the Australasian Amblyomma plus the Amblyomma species from the southern cone of South America, might be sister-group to the Amblyomma of the rest of the world. However, a stable phylogenetic framework of Amblyomma for a better understanding of the biogeographic patterns underpinning its diversification is lacking. METHODS: We used genomic techniques to sequence complete and nearly complete mitochondrial genomes -ca. 15 kbp- as well as the nuclear ribosomal cluster -ca. 8 kbp- for 17 Amblyomma ticks in order to study the phylogeny and biogeographic pattern of the genus Amblyomma, with particular emphasis on the Neotropical region. The new genomic information generated here together with genomic information available on 43 ticks (22 other Amblyomma species and 21 other hard ticks-as outgroup-) were used to perform probabilistic methods of phylogenetic and biogeographic inferences and time-tree estimation using biogeographic dates. RESULTS: In the present paper, we present the strongest evidence yet that Australasian Amblyomma may indeed be the sister-group to the Amblyomma of the rest of the world (species that occur mainly in the Neotropical and Afrotropical zoogeographic regions). Our results showed that all Amblyomma subgenera (Cernyomma, Anastosiella, Xiphiastor, Adenopleura, Aponomma and Dermiomma) are not monophyletic, except for Walkeriana and Amblyomma. Likewise, our best biogeographic scenario supports the origin of Amblyomma and its posterior diversification in the southern hemisphere at 47.8 and 36.8 Mya, respectively. This diversification could be associated with the end of the connection of Australasia and Neotropical ecoregions by the Antarctic land bridge. Also, the biogeographic analyses let us see the colonization patterns of some neotropical Amblyomma species to the Nearctic. CONCLUSIONS: We found strong evidence that the main theater of diversification of Amblyomma was the southern hemisphere, potentially driven by the Antarctic Bridge's intermittent connection in the late Eocene. In addition, the subgeneric classification of Amblyomma lacks evolutionary support. Future studies using denser taxonomic sampling may lead to new findings on the phylogenetic relationships and biogeographic history of Amblyomma genus.

Characterisation and comparative genomics of three new Varanus-associated Borrelia spp. from Indonesia and Australia.

BACKGROUND: Borrelia are important disease-causing tick- and louse-borne spirochaetes than can infect a wide variety of vertebrates, including humans and reptiles. Reptile-associated (REP) Borrelia, once considered a peculiarity, are now recognised as a distinct and important evolutionary lineage, and are increasingly being discovered worldwide in association with novel hosts. Numerous novel Borrelia spp. associated with monitor lizards (Varanus spp.) have been recently identified throughout the Indo-Pacific region; however, there is a lack of genomic data on these Borrelia. METHODS: We used metagenomic techniques to sequence almost complete genomes of novel Borrelia spp. from Varanus varius and Varanus giganteus from Australia, and used long- and short-read technologies to sequence the complete genomes of two strains of a novel Borrelia sp. previously isolated from ticks infesting Varanus salvator from Indonesia. We investigated intra- and interspecies genomic diversity, including plasmid diversity and relatedness, among Varanus-associated Borrelia and other available REP Borrelia and, based on 712 whole genome orthologues, produced the most complete phylogenetic analysis, to the best of our knowledge, of REP Borrelia to date. RESULTS: The genomic architecture of Varanus-associated Borrelia spp. is similar to that of Borrelia spp. that cause relapsing fever (RF), and includes a highly conserved megaplasmid and numerous smaller linear and circular plasmids that lack structural consistency between species. Analysis of PF32 and PF57/62 plasmid partitioning genes indicated that REP Borrelia plasmids fall into at least six distinct plasmid families, some of which are related to previously defined Borrelia plasmid families, whereas the others appear to be unique. REP Borrelia contain immunogenic variable major proteins that are homologous to those found in Borrelia spp. that cause RF, although they are limited in copy number and variability and have low sequence identities to RF variable major proteins. Phylogenetic analyses based on single marker genes and 712 single copy orthologs also definitively demonstrated the monophyly of REP Borrelia as a unique lineage. CONCLUSIONS: In this work we present four new genomes from three novel Borrelia, and thus double the number of REP Borrelia genomes publicly available. The genomic characterisation of these Borrelia clearly demonstrates their distinctiveness as species, and we propose the names Borrelia salvatorii, 'Candidatus Borrelia undatumii', and 'Candidatus Borrelia rubricentralis' for them.

Seventy-eight entire mitochondrial genomes and nuclear rRNA genes provide insight into the phylogeny of the hard ticks, particularly the Haemaphysalis species, Africaniella transversale and Robertsicus elaphensis.

Hoogstraal and Kim (1985) proposed from morphology, three groups of Haemaphysalis subgenera: (i) the "structurally advanced"; (ii) the "structurally intermediate"; and (iii) the "structurally primitive" subgenera. Nuclear gene phylogenies, however, did not indicate monophyly of these morphological groups but alas, only two mitochondrial (mt) genomes from the "structurally intermediate" subgenera had been sequenced. The phylogeny of Haemaphysalis has not yet been resolved. We aimed to resolve the phylogeny of the genus Haemaphysalis, with respect to the subgenus Alloceraea. We presented 15 newly sequenced and annotated mt genomes from 15 species of ticks, five species of which have not been sequenced before, and four new 18S rRNA and 28S rRNA nuclear gene sequences. Our datasets were constructed from 10 mt protein-coding genes, cox1, and the 18S and 28S nuclear rRNA genes. We found a 132-bp insertion between tRNA-Glu (E) gene and the nad1 gene in the mt genome of Haemaphysalis (Alloceraea) inermis that resembles insertions in H. (Alloceraea) kitaokai and Rhipicephalus (Boophilus) geigyi. Our mt phylogenies had the three species of Amblyomma (Aponomma) we sequenced embedded in the main clade of Amblyomma: Am. (Aponomma) fimbriatum, Am. (Aponomma) gervaisi and Am. (Aponomma) latum. This is further support for the hypothesis that the evolution of eyes appears to have occurred in the most-recent-common-ancestor of Amblyocephalus (i.e. Amblyomminae plus Rhipicephalinae) and that eyes were subsequently lost in the most-recent-common-ancestor of the subgenus Am. (Aponomma). Either Africaniella transversale or Robertsicus elaphensis, or perhaps Af. transversale plus Ro. elaphensis, appear to be the sister-group to the rest of the metastriate Ixodida. Our cox1 phylogenies did not indicate monophyly of the "structurally primitive", "structurally intermediate" nor the "structurally advanced" groups of Haemaphysalis subgenera. Indeed, the subgenus Alloceraea may be the only monophyletic subgenus of the genus Haemaphysalis sequenced thus far. All of our mt genome and cox1 phylogenies had the subgenus Alloceraea in a clade that was separate from the rest of the Haemaphysalis ticks. If Alloceraea is indeed the sister to the rest of the Haemaphysalis subgenera this would resonate with the argument of Hoogstraal and Kim (1985), that Alloceraea was a subgenus of "primitive" Haemaphysalis. Alectorobius capensis from Japan had a higher genetic-identity to A. sawaii, which was also from Japan, than to the A. capensis from South Africa. This indicates that A. capensis from Japan may be a cryptic species with respect to the A. capensis from South Africa.

A new subgenus, Australixodes n. subgen. (Acari: Ixodidae), for the kiwi tick, Ixodes anatis Chilton, 1904, and validation of the subgenus Coxixodes Schulze, 1941 with a phylogeny of 16 of the 22 subgenera of Ixodes Latreille, 1795 from entire mitochondrial genome sequences.

A new subgenus, Australixodes n. subgen., is described for the kiwi tick, Ixodes anatis Chilton, 1904. The subgenus Coxixodes Schulze, 1941 is validated for the platypus tick, Ixodes (Coxixodes) ornithorhynchi Lucas, 1846, sister group of the subgenus Endopalpiger Schulze, 1935. A phylogeny from mitochondrial genomes of 16 of the 22 subgenera of Ixodes (32 spp.) indicates, for the first time, the relationships of the subgenera of Ixodes Latreille, 1795, the largest genus of ticks.

Metatranscriptomic profiling reveals diverse tick-borne bacteria, protozoans and viruses in ticks and wildlife from Australia.

Tick-borne zoonoses are emerging globally due to changes in climate and land use. While the zoonotic threats associated with ticks are well studied elsewhere, in Australia, the diversity of potentially zoonotic agents carried by ticks and their significance to human and animal health is not sufficiently understood. To this end, we used untargeted metatranscriptomics to audit the prokaryotic, eukaryotic and viral biomes of questing ticks and wildlife blood samples from two urban and rural sites in New South Wales, Australia. Ixodes holocyclus and Haemaphysalis bancrofti were the main tick species collected, and blood samples from Rattus rattus, Rattus fuscipes, Perameles nasuta and Trichosurus vulpecula were also collected and screened for tick-borne microorganisms using metatranscriptomics followed by conventional targeted PCR to identify important microbial taxa to the species level. Our analyses identified 32 unique tick-borne taxa, including 10 novel putative species. Overall, a wide range of tick-borne microorganisms were found in questing ticks including haemoprotozoa such as Babesia, Theileria, Hepatozoon and Trypanosoma spp., bacteria such as Borrelia, Rickettsia, Ehrlichia, Neoehrlichia and Anaplasma spp., and numerous viral taxa including Reoviridiae (including two coltiviruses) and a novel Flaviviridae-like jingmenvirus. Of note, a novel hard tick-borne relapsing fever Borrelia sp. was identified in questing H. bancrofti ticks which is closely related to, but distinct from, cervid-associated Borrelia spp. found throughout Asia. Notably, all tick-borne microorganisms were phylogenetically unique compared to their relatives found outside Australia, and no foreign tick-borne human pathogens such as Borrelia burgdorferi s.l. or Babesia microti were found. This work adds to the growing literature demonstrating that Australian ticks harbour a unique and endemic microbial fauna, including potentially zoonotic agents which should be further studied to determine their relative risk to human and animal health.

Controversies in bacterial taxonomy: The example of the genus Borrelia.

In this paper we survey key issues in bacterial taxonomy and review the literature regarding the recent genus separation proposed for the genus Borrelia. We discuss how information on members of the genus Borrelia is increasing but detailed knowledge on the relevant features is available only for a small subset of species. The data accumulated here show that there is considerable overlap in ecology, clinical aspects and molecular features between clades that argue against splitting of the genus Borrelia.

Evaluation of 16S next-generation sequencing of hypervariable region 4 in wastewater samples: An unsuitable approach for bacterial enteric pathogen identification.

Recycled wastewater can carry human-infectious microbial pathogens and therefore wastewater treatment strategies must effectively eliminate pathogens before recycled wastewater is used to supplement drinking and agricultural water supplies. This study characterised the bacterial composition of four wastewater treatment plants (WWTPs) (three waste stabilisation ponds and one oxidation ditch WWTP using activated sludge treatment) in Western Australia. The hypervariable region 4 (V4) of the bacterial 16S rRNA (16S) gene was sequenced using next-generation sequencing (NGS) on the Illumina MiSeq platform. Sequences were pre-processed in USEARCH v10.0 and denoised into zero-radius taxonomic units (ZOTUs) with UNOISE3. Taxonomy was assigned to the ZOTUs using QIIME 2 and the Greengenes database and cross-checked with the NCBI nr/nt database. Bacterial composition of all WWTPs and treatment stages (influent, intermediate and effluent) were dominated by Proteobacteria (29.0-87.4%), particularly Betaproteobacteria (9.0-53.5%) and Gammaproteobacteria (8.6-34.6%). Nitrifying bacteria (Nitrospira spp.) were found only in the intermediate and effluent of the oxidation ditch WWTP, and denitrifying and floc-forming bacteria were detected in all WWTPs, particularly from the families Comamonadaceae and Rhodocyclales. Twelve pathogens were assigned taxonomy by the Greengenes database, but comparison of sequences from genera and families known to contain pathogens to the NCBI nr/nt database showed that only three pathogens (Arcobacter venerupis, Laribacter hongkongensis and Neisseria canis) could be identified in the dataset at the V4 region. Importantly, Enterobacteriaceae genera could not be differentiated. Family level taxa assigned by Greengenes database agreed with NCBI nr/nt in most cases, however, BLAST analyses revealed erroneous taxa in Greengenes database. This study highlights the importance of validating taxonomy of NGS sequences with databases such as NCBI nr/nt, and recommends including the V3 region of 16S in future short amplicon NGS studies that aim to identify bacterial enteric pathogens, as this will improve taxonomic resolution of most, but not all, Enterobacteriaceae species.

Genome-wide analysis of Borrelia turcica and 'Candidatus Borrelia tachyglossi' shows relapsing fever-like genomes with unique genomic links to Lyme disease Borrelia.

Borrelia are tick-borne bacteria that in humans are the aetiological agents of Lyme disease and relapsing fever. Here we present the first genomes of B. turcica and B. tachyglossi, members of a recently described and rapidly expanding Borrelia clade associated with reptile (B. turcica) or echidna (B. tachyglossi) hosts, transmitted by hard ticks, and of unknown pathogenicity. Borrelia tachyglossi and B. turcica genomes are similar to those of relapsing fever Borrelia species, containing a linear ~ 900 kb chromosome, a single long (> 70 kb) linear plasmid, and numerous short (< 40 kb) linear and circular plasmids, as well as a suite of housekeeping and macronutrient biosynthesis genes which are not found in Lyme disease Borrelia. Additionally, both B. tachyglossi and B. turcica contain paralogous vsp and vlp proteins homologous to those used in the multiphasic antigen-switching system used by relapsing fever Borrelia to evade vertebrate immune responses, although their number was greatly reduced compared to human-infectious species. However, B. tachyglossi and B. turcica chromosomes also contain numerous genes orthologous to Lyme disease Borrelia-specific genes, demonstrating a unique evolutionary, and potentially phenotypic link between these groups. Borrelia tachyglossi and B. turcica genomes also have unique genetic features, including degraded and deleted tRNA modification genes, and an expanded range of macronutrient salvage and biosynthesis genes compared to relapsing fever and Lyme disease Borrelia. These genomes and genomic comparisons provide an insight into the biology and evolutionary origin of these Borrelia, and provide a valuable resource for future work.

Recent insights into the tick microbiome gained through next-generation sequencing.

The tick microbiome comprises communities of microorganisms, including viruses, bacteria and eukaryotes, and is being elucidated through modern molecular techniques. The advent of next-generation sequencing (NGS) technologies has enabled the genes and genomes within these microbial communities to be explored in a rapid and cost-effective manner. The advantages of using NGS to investigate microbiomes surpass the traditional non-molecular methods that are limited in their sensitivity, and conventional molecular approaches that are limited in their scalability. In recent years the number of studies using NGS to investigate the microbial diversity and composition of ticks has expanded. Here, we provide a review of NGS strategies for tick microbiome studies and discuss the recent findings from tick NGS investigations, including the bacterial diversity and composition, influential factors, and implications of the tick microbiome.