"Escalibur"-A practical pipeline for the de novo analysis of nucleotide variation in nonmodel eukaryotes.

The revolution in genomics has enabled large-scale population genetic investigations of a wide range of organisms, but there has been a relatively limited focus on improving analytical pipelines. To efficiently analyse large data sets, highly integrated and automated software pipelines, which are easy to use, efficient, reliable, reproducible and run in multiple computational environments, are required. A number of software workflows have been developed to handle and process such data sets for population genetic analyses, but effective, specialized pipelines for genetic and statistical analyses of nonmodel organisms are lacking. For most species, resources for variomes (sets of genetic variations found in populations of species) are not available, and/or genome assemblies are often incomplete and fragmented, complicating the selection of the most suitable reference genome when multiple assemblies are available. Additionally, the biological samples used often contain extraneous DNA from sources other than the species under investigation (e.g., microbial contamination), which needs to be removed prior to genetic analyses. For these reasons, we established a new pipeline, called Escalibur, which includes: functionalities, such as data trimming and mapping; selection of a suitable reference genome; removal of contaminating read data; recalibration of base calls; and variant-calling. Escalibur uses a proven gatk variant caller and workflow description language (WDL), and is, therefore, a highly efficient and scalable pipeline for the genome-wide identification of nucleotide variation in eukaryotes. This pipeline is available at https://gitlab.unimelb.edu.au/bioscience/escalibur (version 0.3-beta) and is essentially applicable to any prokaryote or eukaryote.

Metaphor-A workflow for streamlined assembly and binning of metagenomes.

Recent advances in bioinformatics and high-throughput sequencing have enabled the large-scale recovery of genomes from metagenomes. This has the potential to bring important insights as researchers can bypass cultivation and analyze genomes sourced directly from environmental samples. There are, however, technical challenges associated with this process, most notably the complexity of computational workflows required to process metagenomic data, which include dozens of bioinformatics software tools, each with their own set of customizable parameters that affect the final output of the workflow. At the core of these workflows are the processes of assembly-combining the short-input reads into longer, contiguous fragments (contigs)-and binning, clustering these contigs into individual genome bins. The limitations of assembly and binning algorithms also pose different challenges depending on the selected strategy to execute them. Both of these processes can be done for each sample separately or by pooling together multiple samples to leverage information from a combination of samples. Here we present Metaphor, a fully automated workflow for genome-resolved metagenomics (GRM). Metaphor differs from existing GRM workflows by offering flexible approaches for the assembly and binning of the input data and by combining multiple binning algorithms with a bin refinement step to achieve high-quality genome bins. Moreover, Metaphor generates reports to evaluate the performance of the workflow. We showcase the functionality of Metaphor on different synthetic datasets and the impact of available assembly and binning strategies on the final results.

Genetic, Morphological and Antigenic Relationships between Mesonivirus Isolates from Australian Mosquitoes and Evidence for Their Horizontal Transmission.

The Mesoniviridae are a newly assigned family of viruses in the order Nidovirales. Unlike other nidoviruses, which include the Coronaviridae, mesoniviruses are restricted to mosquito hosts and do not infect vertebrate cells. To date there is little information on the morphological and antigenic characteristics of this new group of viruses and a dearth of mesonivirus-specific research tools. In this study we determined the genetic relationships of recent Australian isolates of Alphamesonivirus 4 (Casuarina virus-CASV) and Alphamesonivirus 1 (Nam Dinh virus-NDiV), obtained from multiple mosquito species. Australian isolates of NDiV showed high-level similarity to the prototype NDiV isolate from Vietnam (99% nucleotide (nt) and amino acid (aa) identity). Isolates of CASV from Central Queensland were genetically very similar to the prototype virus from Darwin (95-96% nt and 91-92% aa identity). Electron microscopy studies demonstrated that virion diameter (≈80 nm) and spike length (≈10 nm) were similar for both viruses. Monoclonal antibodies specific to CASV and NDiV revealed a close antigenic relationship between the two viruses with 13/34 mAbs recognising both viruses. We also detected NDiV RNA on honey-soaked nucleic acid preservation cards fed on by wild mosquitoes supporting a possible mechanism of horizontal transmission between insects in nature.

Murine and related chapparvoviruses are nephro-tropic and produce novel accessory proteins in infected kidneys.

Mouse kidney parvovirus (MKPV) is a member of the provisional genus Chapparvovirus that causes renal disease in immune-compromised mice, with a disease course reminiscent of polyomavirus-associated nephropathy in immune-suppressed kidney transplant patients. Here we map four major MKPV transcripts, created by alternative splicing, to a common initiator region, and use mass spectrometry to identify "p10" and "p15" as novel chapparvovirus accessory proteins produced in MKPV-infected kidneys. p15 and the splicing-dependent putative accessory protein NS2 are conserved in all near-complete amniote chapparvovirus genomes currently available (from mammals, birds and a reptile). In contrast, p10 may be encoded only by viruses with >60% amino acid identity to MKPV. We show that MKPV is kidney-tropic and that the bat chapparvovirus DrPV-1 and a non-human primate chapparvovirus, CKPV, are also found in the kidneys of their hosts. We propose, therefore, that many mammal chapparvoviruses are likely to be nephrotropic.

An Atypical Parvovirus Drives Chronic Tubulointerstitial Nephropathy and Kidney Fibrosis.

The occurrence of a spontaneous nephropathy with intranuclear inclusions in laboratory mice has puzzled pathologists for over 4 decades, because its etiology remains elusive. The condition is more severe in immunodeficient animals, suggesting an infectious cause. Using metagenomics, we identify the causative agent as an atypical virus, termed "mouse kidney parvovirus" (MKPV), belonging to a divergent genus of Parvoviridae. MKPV was identified in animal facilities in Australia and North America, is transmitted via a fecal-oral or urinary-oral route, and is controlled by the adaptive immune system. Detailed analysis of the clinical course and histopathological features demonstrated a stepwise progression of pathology ranging from sporadic tubular inclusions to tubular degeneration and interstitial fibrosis and culminating in renal failure. In summary, we identify a widely distributed pathogen in laboratory mice and establish MKPV-induced nephropathy as a new tool for elucidating mechanisms of tubulointerstitial fibrosis that shares molecular features with chronic kidney disease in humans.

Detection of Toscana virus from an adult traveler returning to Australia with encephalitis.

Toscana virus (TOSV) is identified in sandflies, animals, and humans around the Mediterranean Sea. TOSV has not been reported in Australia. During investigations of cerebrospinal fluid samples from patients with encephalitis, TOSV genetic sequences were identified in a traveler returning to Australia from Europe. TOSV should be considered, especially during May to October, in travelers to Australia who embarked in countries in and around the Mediterranean Sea and who subsequently present for medical care because of neurological symptoms.

Assessing the genetic diversity of Cu resistance in mine tailings through high-throughput recovery of full-length copA genes.

Characterizing the genetic diversity of microbial copper (Cu) resistance at the community level remains challenging, mainly due to the polymorphism of the core functional gene copA. In this study, a local BLASTN method using a copA database built in this study was developed to recover full-length putative copA sequences from an assembled tailings metagenome; these sequences were then screened for potentially functioning CopA using conserved metal-binding motifs, inferred by evolutionary trace analysis of CopA sequences from known Cu resistant microorganisms. In total, 99 putative copA sequences were recovered from the tailings metagenome, out of which 70 were found with high potential to be functioning in Cu resistance. Phylogenetic analysis of selected copA sequences detected in the tailings metagenome showed that topology of the copA phylogeny is largely congruent with that of the 16S-based phylogeny of the tailings microbial community obtained in our previous study, indicating that the development of copA diversity in the tailings might be mainly through vertical descent with few lateral gene transfer events. The method established here can be used to explore copA (and potentially other metal resistance genes) diversity in any metagenome and has the potential to exhaust the full-length gene sequences for downstream analyses.

Genome sequence of the pathogenic intestinal spirochete brachyspira hyodysenteriae reveals adaptations to its lifestyle in the porcine large intestine.

Brachyspira hyodysenteriae is an anaerobic intestinal spirochete that colonizes the large intestine of pigs and causes swine dysentery, a disease of significant economic importance. The genome sequence of B. hyodysenteriae strain WA1 was determined, making it the first representative of the genus Brachyspira to be sequenced, and the seventeenth spirochete genome to be reported. The genome consisted of a circular 3,000,694 base pair (bp) chromosome, and a 35,940 bp circular plasmid that has not previously been described. The spirochete had 2,122 protein-coding sequences. Of the predicted proteins, more had similarities to proteins of the enteric Escherichia coli and Clostridium species than they did to proteins of other spirochetes. Many of these genes were associated with transport and metabolism, and they may have been gradually acquired through horizontal gene transfer in the environment of the large intestine. A reconstruction of central metabolic pathways identified a complete set of coding sequences for glycolysis, gluconeogenesis, a non-oxidative pentose phosphate pathway, nucleotide metabolism, lipooligosaccharide biosynthesis, and a respiratory electron transport chain. A notable finding was the presence on the plasmid of the genes involved in rhamnose biosynthesis. Potential virulence genes included those for 15 proteases and six hemolysins. Other adaptations to an enteric lifestyle included the presence of large numbers of genes associated with chemotaxis and motility. B. hyodysenteriae has diverged from other spirochetes in the process of accommodating to its habitat in the porcine large intestine.