High throughput whole rumen metagenome profiling using untargeted massively parallel sequencing.

BACKGROUND: Variation of microorganism communities in the rumen of cattle (Bos taurus) is of great interest because of possible links to economically or environmentally important traits, such as feed conversion efficiency or methane emission levels. The resolution of studies investigating this variation may be improved by utilizing untargeted massively parallel sequencing (MPS), that is, sequencing without targeted amplification of genes. The objective of this study was to develop a method which used MPS to generate "rumen metagenome profiles", and to investigate if these profiles were repeatable among samples taken from the same cow. Given faecal samples are much easier to obtain than rumen fluid samples; we also investigated whether rumen metagenome profiles were predictive of faecal metagenome profiles. RESULTS: Rather than focusing on individual organisms within the rumen, our method used MPS data to generate quantitative rumen micro-biome profiles, regardless of taxonomic classifications. The method requires a previously assembled reference metagenome. A number of such reference metagenomes were considered, including two rumen derived metagenomes, a human faecal microflora metagenome and a reference metagenome made up of publically available prokaryote sequences. Sequence reads from each test sample were aligned to these references. The "rumen metagenome profile" was generated from the number of the reads that aligned to each contig in the database. We used this method to test the hypothesis that rumen fluid microbial community profiles vary more between cows than within multiple samples from the same cow. Rumen fluid samples were taken from three cows, at three locations within the rumen. DNA from the samples was sequenced on the Illumina GAIIx. When the reads were aligned to a rumen metagenome reference, the rumen metagenome profiles were repeatable (P < 0.00001) by cow regardless of location of sampling rumen fluid. The repeatability was estimated at 9%, albeit with a high standard error, reflecting the small number of animals in the study. Finally, we compared rumen microbial profiles to faecal microbial profiles. Our hypothesis, that there would be a stronger correlation between faeces and rumen fluid from the same cow than between faeces and rumen fluid from different cows, was not supported by our data (with much greater significance of rumen versus faeces effect than animal effect in mixed linear model). CONCLUSIONS: We have presented a simple and high throughput method of metagenome profiling to assess the similarity of whole metagenomes, and illustrated its use on two novel datasets. This method utilises widely used freeware. The method should be useful in the exploration and comparison of metagenomes.

Enhanced Arbovirus Surveillance with High-Throughput Metatranscriptomic Processing of Field-Collected Mosquitoes.

Surveillance programs are essential for the prevention and control of mosquito-borne arboviruses that cause serious human and animal diseases. Viral metatranscriptomic sequencing can enhance surveillance by enabling untargeted, high-throughput arbovirus detection. We used metatranscriptomic sequencing to screen field-collected mosquitoes for arboviruses to better understand how metatranscriptomics can be utilised in routine surveillance. Following a significant flood event in 2016, more than 56,000 mosquitoes were collected over seven weeks from field traps set up in Victoria, Australia. The traps were split into samples of 1000 mosquitoes or less and sequenced on the Illumina HiSeq. Five arboviruses relevant to public health (Ross River virus, Sindbis virus, Trubanaman virus, Umatilla virus, and Wongorr virus) were detected a total of 33 times in the metatranscriptomic data, with 94% confirmed using reverse transcription quantitative PCR (RT-qPCR). Analysis of metatranscriptomic cytochrome oxidase I (COI) sequences enabled the detection of 12 mosquito and two biting midge species. Screening of the same traps by an established public health arbovirus surveillance program corroborated the metatranscriptomic arbovirus and mosquito species detections. Assembly of genome sequences from the metatranscriptomic data also led to the detection of 51 insect-specific viruses, both known and previously undescribed, and allowed phylogenetic comparison to past strains. We have demonstrated how metatranscriptomics can enhance surveillance by enabling untargeted arbovirus detection, providing genomic epidemiological data, and simultaneously identifying vector species from large, unsorted mosquito traps.

De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea.

Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga's genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host-pathogen interactions.

Re-Evaluation of the Podosphaera tridactyla Species Complex in Australia.

The Podosphaera tridactyla species complex is highly variable morphologically and causes powdery mildew on a wide range of Prunus species, including stone fruit. A taxonomic revision of the Po. tridactyla species complex in 2020 identified 12 species, seven of which were newly characterised. In order to clarify which species of this complex are present in Australia, next generation sequencing was used to isolate the fungal ITS+28S and host matK chloroplast gene regions from 56 powdery mildew specimens of stone fruit and ornamental Prunus species accessioned as Po. tridactyla or Oidium sp. in Australian reference collections. The specimens were collected in Australia, Switzerland, Italy and Korea and were collected from 1953 to 2018. Host species were confirmed using matK phylogenetic analysis, which identified that four had been misidentified as Prunus but were actually Malusprunifolia. Podosphaera species were identified using ITS+28S phylogenetic analysis, recognising three Podosphaera species on stone fruit and related ornamental Prunus hosts in Australia. These were Po.pannosa, the rose powdery mildew, and two species in the Po. tridactyla species complex: Po. ampla, which was the predominant species, and a previously unidentified species from peach, which we describe here as Po. cunningtonii.

Sensitivity and specificity of metatranscriptomics as an arbovirus surveillance tool.

The ability to identify all the viruses within a sample makes metatranscriptomic sequencing an attractive tool to screen mosquitoes for arboviruses. Practical application of this technique, however, requires a clear understanding of its analytical sensitivity and specificity. To assess this, five dilutions (1:1, 1:20, 1:400, 1:8,000 and 1:160,000) of Ross River virus (RRV) and Umatilla virus (UMAV) isolates were spiked into subsamples of a pool of 100 Culex australicus mosquitoes. The 1:1 dilution represented the viral load of one RRV-infected mosquito in a pool of 100 mosquitoes. The subsamples underwent nucleic acid extraction, mosquito-specific ribosomal RNA depletion, and Illumina HiSeq sequencing. The viral load of the subsamples was also measured using reverse transcription droplet digital PCR (RT-ddPCR) and quantitative PCR (RT-qPCR). Metatranscriptomic sequencing detected both RRV and UMAV in the 1:1, 1:20 and 1:400 subsamples. A high specificity was achieved, with 100% of RRV and 99.6% of UMAV assembled contigs correctly identified. Metatranscriptomic sequencing was not as sensitive as RT-qPCR or RT-ddPCR; however, it recovered whole genome information and detected 19 other viruses, including four first detections for Australia. These findings will assist arbovirus surveillance programs in utilising metatranscriptomics in routine surveillance activities to enhance arbovirus detection.

Analysis of the Indole Diterpene Gene Cluster for Biosynthesis of the Epoxy-Janthitrems in Epichloë Endophytes.

Epoxy-janthitrems are a class of indole diterpenes with structural similarity to lolitrem B. Two taxa of asexual Epichloë endophytes have been reported to produce epoxy-janthitrems, LpTG-3 (Lolium perenne Taxonomic Group 3; e.g., NEA12) and LpTG-4 (e.g., E1). Epichloë epoxy-janthitrems are not well understood, the biosynthetic pathway and associated gene complement have not been described and while the literature suggests they are associated with superior protection against pasture insect pests and are tremorgenic in grazing mammals, these properties have not been confirmed using isolated and purified compounds. Whole genome sequence analysis was used to identify candidate genes for epoxy-janthitrem biosynthesis that are unique to epoxy-janthitrem producing strains of Epichloë. A gene, jtmD, was identified with homology to aromatic prenyl transferases involved in synthesis of indole diterpenes. The location of the epoxy-janthitrem biosynthesis gene cluster (JTM locus) was determined in the assembled nuclear genomes of NEA12 and E1. The JTM locus contains cluster 1 and cluster 2 of the lolitrem B biosynthesis gene cluster (LTM locus), as well as four genes jtmD, jtmO, jtm01, and jtm02 that are unique to Epichloë spp. that produce epoxy-janthitrems. Expression of each of the genes identified was confirmed using transcriptome analysis of perennial ryegrass-NEA12 and perennial ryegrass-E1 symbiota. Sequence analysis confirmed the genes are functionally similar to those involved in biosynthesis of related indole diterpene compounds. RNAi silencing of jtmD and in planta assessment in host-endophyte associations confirms the role of jtmD in epoxy-janthitrem production. Using LCMS/MS technologies, a biosynthetic pathway for the production of epoxy-janthitrems I-IV in Epichloë endophytes is proposed.

Metagenomic arbovirus detection using MinION nanopore sequencing.

With its small size and low cost, the hand-held MinION sequencer is a powerful tool for in-field surveillance. Using a metagenomic approach, it allows non-targeted detection of viruses in a sample within a few hours. This study aimed to determine the ability of the MinION to metagenomically detect and characterise a virus from an infected mosquito. RNA was extracted from an Aedes notoscriptus mosquito infected with Ross River virus (RRV), converted into cDNA and sequenced on the MinION. Bioinformatic analysis of the MinION reads led to detection of full-length RRV, with reads of up to 2.5kb contributing to the assembly. The cDNA was also sequenced on the MiSeq sequencer, and both platforms recovered the RRV genome with >98% accuracy. This proof of concept study demonstrates the metagenomic detection of an arbovirus, using the MinION, directly from a mosquito with minimal sample purification.