Technical strategy for monozygotic twin discrimination by single-nucleotide variants.

Monozygotic (MZ) twins are theoretically genetically identical. Although they are revealed to accumulate mutations after the zygote splits, discriminating between twin genomes remains a formidable challenge in the field of forensic genetics. Single-nucleotide variants (SNVs) are responsible for a substantial portion of genetic variation, thus potentially serving as promising biomarkers for the identification of MZ twins. In this study, we sequenced the whole genome of a pair of female MZ twins when they were 27 and 33 years old to approximately 30 × coverage using peripheral blood on an Illumina NovaSeq 6000 Sequencing System. Potentially discordant SNVs supported by whole-genome sequencing were validated extensively by amplicon-based targeted deep sequencing and Sanger sequencing. In total, we found nine bona fide post-twinning SNVs, all of which were identified in the younger genomes and found in the older genomes. None of the SNVs occurred within coding exons, three of which were observed in introns, supported by whole-exome sequencing results. A double-blind test was employed, and the reliability of MZ twin discrimination by discordant SNVs was endorsed. All SNVs were successfully detected when input DNA amounts decreased to 0.25 ng, and reliable detection was limited to seven SNVs below 0.075 ng input. This comprehensive analysis confirms that SNVs could serve as cost-effective biomarkers for MZ twin discrimination.

Metabarcoding study to reveal the structural community of strongylid nematodes in domesticated horses in Thailand.

BACKGROUND: Mixed strongylid infections significantly impact equine health and performance. Traditional microscopy-based methods exhibit limitations in accurately identifying strongylid species. Nemabiome deep amplicon sequencing approach previously succeeded in describing the strongylid communities in livestock including equids. However, there are no available studies that describe the structural communities of strongylid parasites in horses in Thailand. Therefore, this study was undertaken encompassing the ITS-2 rDNA metabarcoding assay to characterize strongylid species within horse fecal samples collected from a cohort of yearlings at the largest domesticated stud farm in Thailand. In addition, to investigate the capability of ITS-2 rDNA in assessing the phylogenetic relationships among the identified strongylid species. RESULTS: The study identified 14 strongylid species in the examined equine populations, each with varying prevalence. Notably, Cylicocyclus nassatus and Cylicostephanus longibursatus were identified as the predominant species, with Strongylus spp. conspicuously absent. The phylogenetic analysis of 207 amplicon sequence variants (ASVs) displayed a complex relationship among the investigated cyathostomin species, with some species are positioned across multiple clades, demonstrating close associations with various species and genera. CONCLUSION: The ITS-2 nemabiome sequencing technique provided a detailed picture of horse strongylid parasite species in the studied population. This establishes a foundation for future investigations into the resistance status of these parasites and enables efforts to mitigate their impact.

Establishment of a PEG-mediated protoplast transformation system based on DNA and CRISPR/Cas9 ribonucleoprotein complexes for banana.

BACKGROUND: To date, CRISPR/Cas9 RNP editing tools have not been applied to the genetic modification of banana. Here, the establishment of a PEG-mediated banana protoplast transformation system makes it possible to build an efficient DNA-free method for a site-directed mutagenesis system. RESULTS: Protoplasts constitute a versatile platform for transient expression in plant science. In this study, we established a PEG-mediated banana protoplast transformation system. This system was further optimized for successfully delivering CRISPR/Cas9 and CRISPR/Cas12a plasmids and CRISPR/Cas9 ribonucleoproteins (RNPs) for targeted delivery of the PDS gene into banana protoplasts. Specific bands were observed in PCR-Restriction Enzyme Digestion (PCR-RE) assays, and Sanger sequencing of single clones further confirmed the occurrence of indels at target sites. Deep amplicon sequencing results showed that the editing efficiency of the CRISPR/Cas9 system was higher than that of the other two systems. CONCLUSIONS: The PEG-mediated banana protoplast transformation system can serve as a rapid and effective tool for transient expression assays and sgRNA validation in banana. The application of the CRISPR/Cas9 RNP system enables the generation of banana plants engineered by DNA-free gene editing.

Assessing anthelmintic resistance risk in the post-genomic era: a proof-of-concept study assessing the potential for widespread benzimidazole-resistant gastrointestinal nematodes in North American cattle and bison.

As genomic research continues to improve our understanding of the genetics of anthelmintic drug resistance, the revolution in DNA sequencing technologies will provide increasing opportunities for large-scale surveillance for the emergence of drug resistance. In most countries, parasite control in cattle and bison has mainly depended on pour-on macrocyclic lactone formulations resulting in widespread ivermectin resistance. Consequently, there is an increased interest in using benzimidazole drugs which have been used comparatively little in cattle and bison in recent years. This situation, together with our understanding of benzimidazole resistance genetics, provides a practical opportunity to use deep-amplicon sequencing to assess the risk of drug resistance emergence. In this paper, we use deep-amplicon sequencing to scan for those mutations in the isotype-1 ß-tubulin gene previously associated with benzimidazole resistance in many trichostrongylid nematode species. We found that several of these mutations occur at low frequency in many cattle and bison parasite populations in North America, suggesting increased use of benzimidazole drugs in cattle has the potential to result in widespread emergence of resistance in multiple parasite species. This work illustrates a post-genomic approach to large-scale surveillance of early emergence of anthelmintic resistance in the field.

Ultra-deep amplicon sequencing indicates absence of low-grade mosaicism with normal cells in patients with type-1 NF1 deletions.

Different types of large NF1 deletion are distinguishable by breakpoint location and potentially also by the frequency of mosaicism with normal cells lacking the deletion. However, low-grade mosaicism with fewer than 10% normal cells has not yet been excluded for all NF1 deletion types since it is impossible to assess by the standard techniques used to identify such deletions, including MLPA and array analysis. Here, we used ultra-deep amplicon sequencing to investigate the presence of normal cells in the blood of 20 patients with type-1 NF1 deletions lacking mosaicism according to MLPA. The ultra-deep sequencing entailed the screening of 96 amplicons for heterozygous SNVs located within the NF1 deletion region. DNA samples from three previously identified patients with type-2 NF1 deletions and low-grade mosaicism with normal cells as determined by FISH or microsatellite marker analysis were used to validate our methodology. In these type-2 NF1 deletion samples, proportions of 5.3%, 6.6% and 15.0% normal cells, respectively, were detected by ultra-deep amplicon sequencing. However, using this highly sensitive method, none of the 20 patients with type-1 NF1 deletions included in our analysis exhibited low-grade mosaicism with normal cells in blood, thereby supporting the view that the vast majority of type-1 deletions are germline deletions.

The immune response of the scallop Argopecten purpuratus is associated with changes in the host microbiota structure and diversity.

All organisms live in close association with a variety of microorganisms called microbiota. Furthermore, several studies support a fundamental role of the microbiota on the host health and homeostasis. In this context, the aim of this work was to determine the structure and diversity of the microbiota associated with the scallop Argopecten purpuratus, and to assess changes in community composition and diversity during the host immune response. To do this, adult scallops were immune challenged and sampled after 24 and 48 h. Activation of the immune response was established by transcript overexpression of several scallop immune response genes in hemocytes and gills, and confirmed by protein detection of the antimicrobial peptide big defensin in gills of Vibrio-injected scallops at 24 h post-challenge. Then, the major bacterial community profile present in individual scallops was assessed by denaturing gradient gel electrophoresis (DGGE) of 16S rDNA genes and dendrogram analyses, which indicated a clear clade differentiation of the bacterial communities noticeable at 48 h post-challenge. Finally, the microbiota structure and diversity from pools of scallops were characterized using 16S deep amplicon sequencing. The results revealed an overall modulation of the microbiota abundance and diversity according to scallop immune status, allowing for prediction of some changes in the functional potential of the microbial community. Overall, the present study showed that changes in the structure and diversity of bacterial communities associated with the scallop A. purpuratus are detected after the activation of the host immune response. Now, the relevance of microbial balance disruption in the immune capacity of the scallop remains to be elucidated.

Microbiota Analysis of an Environmental Slurry and Its Potential Role as a Reservoir of Bovine Digital Dermatitis Pathogens.

At present, very little information exists regarding what role the environmental slurry may play as an infection reservoir and/or route of transmission for bovine digital dermatitis (DD), a disease which is a global problem in dairy herds. To investigate whether DD-related bacteria belong to the indigenous microbiota of the dairy herd environment, we used deep amplicon sequencing of the 16S rRNA gene in 135 slurry samples collected from different sites in 22 dairy farms, with and without DD-infected cows. Both the general bacterial populations and digital dermatitis-associated Treponema were targeted in this study. The results revealed significant differences in the bacterial communities between the herds, with only 12 bacterial taxa shared across at least 80% of all the individual samples. These differences in the herd microbiota appeared to reflect mainly between-herd variation. Not surprisingly, the slurry was dominated by ubiquitous gastrointestinal bacteria, such as Ruminococcaceae and Lachnospiraceae Despite the low relative abundance of spirochetes, which ranged from 0 to 0.6%, we were able to detect small amounts of bacterial DNA from DD-associated treponemes in the slurry. However, the DD-associated Treponema spp. were detected only in samples from herds with reported DD problems. These data indicate that treponemes involved in the pathogenesis of DD are not part of the normal environmental microflora in dairy herds without clinical DD and, consequently, that slurry is not a primary reservoir of infection.IMPORTANCE Bovine digital dermatitis (DD), a dermal disease which causes lameness in dairy cattle, is a serious problem worldwide. To control this disease, the infection reservoirs and transmission routes of DD pathogens need to be clarified. The dairy herd slurry may be a pathogen reservoir of DD-associated bacteria. The rationale for the present study was, therefore, to examine whether DD-associated bacteria are always present in slurry or if they are found only in DD-afflicted herds. The results strongly indicated that DD Treponema spp. are not part of the indigenous slurry and, therefore, do not comprise an infection reservoir in healthy herds. This study applied next-generation sequencing technology to decipher the microbial compositions of environmental slurry of dairy herds with and without digital dermatitis.

Metagenomic analysis of oomycete communities from the rhizosphere of field pea on the Canadian prairies.

Oomycetes are a diverse group of microorganisms; however, little is known about their composition and biodiversity in agroecosystems. Illumina MiSeq was used to determine the type and abundance of oomycetes associated with pea root rot in the Canadian prairies. Additional objectives of the study were to identify differences in oomycete communities associated with pea root health and compare oomycete communities among the 3 prairie provinces, where field peas are commonly cultivated. Samples of soil from the rhizosphere of field pea (Pisum sativum L.) were collected from patches of asymptomatic or diseased plants from 26 commercial fields in 2013 and 2014. Oomycete communities were characterized using metagenomic analysis of the ITS1 region on Illumina MiSeq. From 105 identified operational taxonomic units (OTUs), 45 and 16 oomycete OTUs were identified at species and genus levels, respectively. Pythium was the most prevalent genus and Pythium heterothallicum the most prevalent species in all 3 provinces in both 2013 and 2014. Aphanomyces euteiches, a very important pea root rot pathogen in regions of the prairies, was detected in 57% of sites but at very low abundance (<0.2%). Multivariate analysis revealed differences in the relative abundance of species in oomycete communities between asymptomatic and diseased sites, and among years and provinces. This study demonstrated that deep amplicon sequencing can provide information on the composition and diversity of oomycete communities in agricultural soils.

The molecular spectrum and clinical impact of DIS3 mutations in multiple myeloma.

Multiple myeloma (MM) is a plasma cell neoplasm that presents with a major biological and clinical heterogeneity. We here investigated the spectrum of clonal and subclonal mutations of DIS3, an active part of the exosome complex, that may play a role in the development or progression of MM. The whole coding sequence of DIS3 was subjected to deep sequencing in 81 uniformly-treated MM patients and 12 MM cell lines and the overall occurrence of DIS3 mutations as well as the presence of DIS3 mutations in minor and major subclones were correlated with cytogenetic alterations and clinical parameters. Our study identified DIS3 mutations in 9/81 patients that were associated with 13q14 deletions and IGH translocations on the cytogenetic level. Specifically, we detected seven novel somatic DIS3 single nucleotide variants (SNVs) and defined three hot spot mutations within the RNB domain. Lastly, we found a trend towards a shorter median overall survival for patients with DIS3 mutations, and patients carrying DIS3 mutations in minor subclones of their tumours showed a significantly worse response to therapy compared to patients with DIS3 mutations in the major subclone.

Gastrointestinal nematodes in German outdoor-reared pigs based on faecal egg count and next-generation sequencing nemabiome data.

BACKGROUND: There is a higher risk for nematode infections associated with outdoor-reared pigs. Next to Ascaris suum, Oesophagostomum dentatum and Trichuris suis, there is the potential of infections with other nodular worm species, Hyostrongylus rubidus, Stongyloides ransomi and Metastrongylus spp. lungworms. Next-generation sequencing methods describing the nemabiome have not yet been established for porcine nematodes. METHODS: FLOTAC was used for faecal egg counts of porcine gastrointestinal nematodes and lungworms in piglets, fatteners and adults individually. A nemabiome analyses based on ITS-2 gene region metabarcoding was used to differentiate strongyle species. Additionally, questionnaire data was analysed using mixed-effect regression to identify potential risk factors associated with parasite occurrences and egg shedding intensity. RESULTS: On 15 of 17 farms nematode eggs were detected. Ascaris suum, strongyles and T. suis were detected on 82%, 70% and 35% of the 17 farms, respectively. Lungworms were detected on one out of four farms with access to pasture. Strongyloides ransomi was not detected. 32% (CI 28-36%), 27% (24-31%), 5% (4-7%) and 3% (0.9-8%) of the samples where tested positive for strongyles, A. suum, T. suis and lungworms, respectively. The nemabiome analysis revealed three different strongyle species, with O. dentatum being the most common (mean 93.9%), followed by O. quadrispinulatum (5.9%) and the hookworm Globocephalus urosubulatus (0.1%). The bivariate and multivariate risk factor analyses showed among others that cleaning once a week compared to twice a week increased the odds significantly for being infected with A. suum (OR 78.60) and strongyles (2077.59). Access to pasture was associated with higher odds for A. suum (43.83) and strongyles (14.21). Compared to shallow litter systems, deep litter and free range systems resulted in significant higher odds for strongyles (85.74, 215.59, respectively) and T. suis (200.33, 623.08). CONCLUSIONS: Infections with A. suum, O. dentatum, O. quadrispinulatum, T. suis, Metastrongylus spp. and G. urosubulatus are present in German outdoor-reared pigs. This is the first report of G. urosubulatus in domestic pigs in Europe. Metabarcoding based on the ITS-2 region is a suitable tool to analyse the porcine nemabiome. Furthermore, management practices have the potential of reducing the risk of parasite infections.

Deep Amplicon Sequencing Reveals Culture Selection of Mycobacterium Tuberculosis by Clinical Samples.

The commonly-used drug susceptibility testing (DST) relies on bacterial culture and faces shortcomings such as long turnaround time and clone/subclone selection. We developed a targeted deep amplification sequencing (DAS) method directly applied to clinical specimens. In this DAS panel, we examined 941 drug-resistant mutations associated with 20 anti-tuberculosis drugs with an initial amount of 4 pg DNA and reduced clinical testing time from 20 days to two days. A prospective study was conducted using 115 clinical specimens mainly with Xpert® Mycobacterium tuberculosis/rifampicin (Xpert MTB/RIF) assay positive to evaluate drug-resistant mutation detection. DAS was performed on culture-free specimens, while culture-dependent isolates were used for phenotypic DST, DAS, and whole-genome sequencing (WGS). For in silico molecular DST, our result based on DAS panel revealed the similar accuracy to three published reports based on WGS. For 82 isolates, application of DAS showed better sensitivity (93.03% vs. 92.16%), specificity (96.10% vs. 95.02%), and accuracy (91.33% vs. 90.62%) than Mykrobe software using WGS. Compared to culture-dependent WGS, culture-free DAS provides a full picture of sequence variation at population level, exhibiting in detail the gain-and-loss variants caused by bacterial culture. Our study performs a systematic verification of the advantages of DAS in clinical applications and comprehensively illustrates the discrepancy in Mycobacterium tuberculosis before and after culture.

Use of Deep-Amplicon Sequencing (DAS), Real-Time PCR and In Situ Hybridization to Detect H. pylori and Other Pathogenic Helicobacter Species in Feces from Children.

BACKGROUND: Detecting Helicobacter pylori in fecal samples is easier and more comfortable than invasive techniques, especially in children. Thus, the objective of the present work was to detect H. pylori in feces from children by molecular methods as an alternative for diagnostic and epidemiological studies. METHODS: Forty-five fecal samples were taken from pediatric patients who presented symptoms compatible with H. pylori infection. HpSA test, culture, real-time quantitative PCR (qPCR), fluorescence in situ hybridization (FISH), direct viable count associated with FISH (DVC-FISH), and Illumina-based deep-amplicon sequencing (DAS) were applied. RESULTS: No H. pylori colonies were isolated from the samples. qPCR analysis detected H. pylori in the feces of 24.4% of the patients. In comparison, DVC-FISH analysis showed the presence of viable H. pylori cells in 53.3% of the samples, 37% of which carried 23S rRNA mutations that confer resistance to clarithromycin. After DAS, H. pylori-specific 16S rDNA sequences were detected in 26 samples. In addition, DNA from H. hepaticus was identified in 10 samples, and H. pullorum DNA was detected in one sample. CONCLUSION: The results of this study show the presence of H. pylori, H. hepaticus, and H. pullorum in children's stools, demonstrating the coexistence of more than one Helicobacter species in the same patient. The DVC-FISH method showed the presence of viable, potentially infective H. pylori cells in a high percentage of the children's stools. These results support the idea that fecal-oral transmission is probably a common route for H. pylori and suggest possible fecal-oral transmission of other pathogenic Helicobacter species.

A novel metabarcoded deep amplicon sequencing tool for disease surveillance and determining the species composition of Trypanosoma in cattle and other farm animals.

The World Health Organization (WHO) and the Food and Agriculture Organization (FAO) have developed strategies to control trypanosomiasis in humans and livestock in endemic areas. These require a better understanding of the distribution of different Trypanosoma species and improved predictions of where they might appear in the future, based on accurate diagnosis and robust surveillance systems. Here, we describe a metabarcoding deep amplicon sequencing method to identify and determine the Trypanosoma species in co-infecting communities. First, four morphological verified Trypanosoma species (T. brucei, T. congolense, T. vivax and T. theileri) were used to prepare test DNA pools derived from different numbers of parasites to evaluate the method's detection threshold for each of the four species and to assess the accuracy of their proportional quantification. Having demonstrated the accurate determination of species composition in Trypanosoma communities, the method was applied to determine its detection threshold using blood samples collected from cattle with confirmed Trypanosoma infections based on a PCR assay. Each sample showed a different Trypanosoma species composition based on the proportion of MiSeq reads. Finally, we applied the assay to field samples to develop new insight into the species composition of Trypanosoma communities in cattle, camels, buffalo, horses, sheep, and goat in endemically infected regions of Pakistan. We confirmed that Trypanosoma evansi is the major species in Pakistan and for the first time showed the presence of Trypanosoma theileri. The metabarcoding deep amplicon sequencing method and bioinformatics pathway have several potential applications in animal and human research, including evaluation of drug treatment responses, understanding of the emergence and spread of drug resistance, and description of species interactions during co-infections and determination of host and geographic distribution of trypanosomiasis in humans and livestock.

Phenotypic and genotypic analysis of benzimidazole resistance in reciprocal genetic crosses of Haemonchus contortus.

Haemonchus contortus is arguably one of the most economically important and ubiquitous parasites of livestock globally and commonly involved in cases of anthelmintic resistance. Here, we performed reciprocal genetic crosses using susceptible (MHco3(ISE)) and multiple anthelmintic resistant (MHco18(UGA2004)) H. contortus isolates. Resultant admixed populations were designated MHco3/18 or MHco18/3, where the lead isolate reflects the origin of the females. Three independent filial generations were generated for each cross, which were subjected to bioassays, molecular approaches and population genetic analyses to investigate the phenotypic and genotypic inheritance of benzimidazole (BZ) resistance at each stage. A panel of microsatellite markers confirmed the success of the genetic cross as markers from both parents were seen in the F1 crosses. Egg hatch tests revealed a stark difference between the two F1 crosses with ED50 estimates for MHco18/3 being 9 times greater than those for MHco3/18. Resistance factors based on ED50 estimates ranged from 6 to 57 fold in the filial progeny compared to MHco3(ISE) parents. Molecular analysis of the F167Y and F200Y SNP markers associated with BZ resistance were analysed by pyrosequencing and MiSeq deep amplicon sequencing, which showed that MHco3/18.F1 and MHco18/3.F1 both had similar frequencies of the F200Y resistant allele (45.3% and 44.3%, respectively), whereas for F167Y, MHco18/3.F1 had a two-fold greater frequency of the resistant-allele compared to MHco3/18.F1 (18.2% and 8.8%, respectively). Comparison between pyrosequencing and MiSeq amplicon sequencing revealed that the allele frequencies derived from both methods were concordant at codon 200 (rc = 0.97), but were less comparable for codon 167 (rc = 0.55). The use of controlled reciprocal genetic crosses have revealed a potential difference in BZ resistance phenotype dependent on whether the resistant allele is paternally or maternally inherited. These findings provide new insight and prompt further investigation into the inheritance of BZ resistance in H. contortus.

Phylogenetic analysis suggests single and multiple origins of dihydrofolate reductase mutations in Plasmodium vivax.

Pyrimethamine was first introduced for the treatment of malaria in Asia and Africa during the early 1980s, replacing chloroquine, and has become the first line of drugs in many countries. In recent years, development of pyrimethamine resistance in Plasmodium vivax has become a barrier to effective malaria control strategies. Here, we describe the use of meta-barcoded deep amplicon sequencing technology to assess the evolutionary origin of pyrimethamine resistance by analysing the flanking region of dihydrofolate reductase (dhfr) locus. The genetic modelling suggests that 58R and 173L single mutants and 58R/117N double mutants are present on a single lineage; suggesting a single origin of these mutations. The triple mutants (57L/58R/117N, 58R/61M/117N and 58R/117N/173L) share the lineage of 58R/117N, suggesting a common origin. In contrast, the 117N mutant is present on two separate lineages suggesting that there are multiple origins of this mutation. We characterised the allele frequency of the P. vivax dhfr locus. Our results support the view that the single mutation of 117N and double mutations of 58R/117N arise commonly, whereas the single mutation of 173L and triple mutations of 57L/58R/117N, 58R/61M/117N and 58R/117N/173L are less common. Our work will help to inform mitigation strategies for pyrimethamine resistance in P. vivax.

Early detection of metastatic uveal melanoma by the analysis of tumor-specific mutations in cell-free plasma DNA.

BACKGROUND: Eye salvaging therapy of malignant melanomas of the uvea can preserve the eye in most cases, but still about half of patients die from metastatic disease. Previous analyses of cell-free DNA from plasma had shown detectable levels of tumor-specific GNAQ/GNA11 mutations in patients with the clinical diagnosis of progressive disease. However, data on the time span that elapses from the detection of ctDNA in plasma to the clinical detection of metastases (diagnostic lead time) are missing. METHODS: We examined 135 patients with uveal melanoma. Cell-free DNA was isolated from a total of 807 blood samples which were taken over a period of up to 41 months and analyzed for the presence of GNAQ/GNA11 mutations by deep amplicon sequencing. RESULTS: Twenty-one of the 135 patients developed metastases or recurrence. A ctDNA signal was identified in the plasma of 17 of the 21 patients. In 10 patients, this ctDNA signal preceded the clinical diagnosis of metastasis by 2-10 months. In 10 other patients, a ctDNA signal was only detected in samples obtained shortly before or after radiotherapy. The presence of a ctDNA signal in 16 of the remaining 125 patients was linked to clinical manifestation of metastases (n = 14) or tumor recurrence (n = 2) with a sensitivity and specificity of 80% and 96%, respectively. CONCLUSION: Detection of ctDNA in plasma can provide a diagnostic lead time over the clinical diagnosis of metastases or tumor recurrence. Longer lead times are to be expected if intervals between sampling are shortened.

Wild ruminants as reservoirs of domestic livestock gastrointestinal nematodes.

Gastrointestinal nematode (GIN) infections in cattle cause appetite suppression which leads to poor feed conversion, reduced weight gain and reduced milk production. Overuse and exclusive reliance on anthelmintic drugs has resulted in widespread resistance in many parasitic nematode species infecting livestock making control increasingly difficult. Wild ruminants are competent hosts of a number of nematode species that typically infect and are best adapted for cattle, sheep, and goats. Thus, the potential exists for wild ruminants to act as reservoirs in the translocation of domestic GIN, including those carrying anthelmintic resistance mutations as well as susceptible genotypes. The potential for parasite exchange is heightened by interfaces or ecotones between managed and wild rangelands, and by perturbations linked to climate warming that can increasingly alter the distributions of wild ungulates and their interactions with domestic and free-ranging ruminants. To investigate the extent to which wild ruminants harbour parasites capable of infecting domestic ruminants we first performed an epidemiological study of feces from wildlife hosts that spanned 16 states and included white-tailed deer (85 % of the samples), pronghorn, elk, mule deer, bighorn sheep, moose, cattle, and caribou across the United States. All samples were cultured to third stage larvae and nematode DNA was isolated and PCR amplified. Among the 548 wild ruminant samples received, 33 % (181 samples) were positive for nematode DNA, among which half (84 samples) contained DNA from GIN species commonly found in cattle. DNA from cattle GIN species was detected in 46 % of samples from the Northeast, 42 % from the Southeast, 10 % from the Midwest, 0 % from the Southwest and 11 % from the West. Deep amplicon sequencing of the ITS-2 rDNA indicated that Ostertagia and Trichostrongylus were present in 90 % and 69 % of the nematode DNA positive samples, respectively, whereas Haemonchus, Cooperia and Oesophagostomum were present in 26 %, 2 % and 10 % of the samples, respectively. These data clearly show that wild ruminants commonly harbour multiple parasite species whose primary hosts are domestic cattle, and suggest that further work is warranted to investigate their specific roles in the management of anthelmintic resistance.

A novel metabarcoded 18S ribosomal DNA sequencing tool for the detection of Plasmodium species in malaria positive patients.

Various PCR based methods have been described for the diagnosis of malaria, but most depend on the use of Plasmodium species-specific probes and primers; hence only the tested species are identified and there is limited available data on the true circulating species diversity. Sensitive diagnostic tools and platforms for their use are needed to detect Plasmodium species in both clinical cases and asymptomatic infections that contribute to disease transmission. We have recently developed for the first time a novel high throughput 'haemoprotobiome' metabarcoded DNA sequencing method and applied it for the quantification of haemoprotozoan parasites (Theleria and Babesia) of livestock. Here, we describe a novel, high throughput method using an Illumina MiSeq platform to demonstrate the proportions of Plasmodium species in metabarcoded DNA samples derived from human malaria patients. Plasmodium falciparum and Plasmodium vivax positive control gDNA was used to prepare mock DNA pools of parasites to evaluate the detection threshold of the assay for each of the two species. The different mock pools demonstrate the accurate detection ability and to show the proportions of each of the species being present. We then applied the assay to malaria-positive human samples to show the species composition of Plasmodium communities in the Punjab province of Pakistan and in the Afghanistan-Pakistan tribal areas. The diagnostic performance of the deep amplicon sequencing method was compared to an immunochromatographic assay that is widely used in the region. The deep amplicon sequencing showed that P. vivax was present in 69.8%, P. falciparum in 29.5% and mixed infection in 0.7% patients examined. The immunochromatographic assay showed that P. vivax was present in 65.6%, P. falciparum in 27.4%, mixed infection 0.7% patients and 6.32% malaria-positive cases were negative in immunochromatographic assay, but positive in the deep amplicon sequencing. Overall, metabarcoded DNA sequencing demonstrates better diagnostic performance, greatly increasing the estimated prevalence of Plasmodium infection. The next-generation sequencing method using metabarcoded DNA has potential applications in the diagnosis, surveillance, treatment, and control of Plasmodium infections, as well as to study the parasite biology.

Development of a deep amplicon sequencing method to determine the species composition of piroplasm haemoprotozoa.

Piroplasmosis is caused by tick-borne haemoprotozoa of the genera Theileria and Babesia. These parasitic infections can seriously impact on the health of livestock and production. Piroplasms of multiple species can be present in a single host, but reliable molecular diagnostic tools are needed in order to understand the composition of these complex parasite communities. Theileria and Babesia vary in their epidemiology, drug sensitivity, pathogenicity and interaction with co-infecting species, but are similar in that infected animals become persistent carriers after recovery from primary infection, acting as reservoir hosts. Here, we describe for the first time the use of a deep amplicon sequencing platform to identify proportions of piroplasm species in co-infecting communities and develop the concept of a "haemoprotobiome". First, four phenotypically-verified species of Theileria and Babesia were used to prepare mock DNA pools with random numbers of the parasites amplified by four different numbers of PCR cycles to assess sequence representation for each species. Second, we evaluated the detection threshold of the deep amplicon sequencing assay for each of the four species and to assess the accuracy of proportional quantification of all four species. Finally, we applied the assay to the field samples to afford insight of the species composition of piroplasm communities in small and large ruminants in the Punjab province of Pakistan. The "haemoprotobiome" concept has several potential applications in veterinary and human research, including understanding of responses to drug treatment; parasite epidemiology and ecology; species interactions during mixed infections; and parasite control strategies.

A high throughput deep amplicon sequencing method to show the emergence and spread of Calicophoron daubneyi rumen fluke infection in United Kingdom cattle herds.

The prevalence of C. daubneyi infection in the United Kingdom has increased, but despite the potential for rumen flukes to cause production loss in ruminant livestock, understanding of their emergence and spread is poor. Here we describe the development of a method to explore the multiplicity of C. daubneyi infection and patterns of the parasite's emergence and spread, based on Illumina MiSeq deep sequencing of meta barcoded amplicons of a fragment of the cytochrome c oxidase subunit I (mt-COX-1) locus. Our results show high levels of genetic diversity in 32 C. daubneyi populations derived from finished prime cattle consigned to slaughter from northern United Kingdom. The results are consistent with a single introduction of C. daubneyi infection to some of the farms where the cattle had been grazed during their lifetime and multiple introductions to most. The results illustrate the impact of high levels of animal movements in the United Kingdom, whereby multiple common mt-COX-1 haplotypes were identified in 26 populations in the absence of geographical clustering of clades. This has implications for the adaptability of environmental and intermediate host stages of the parasite to changing climatic and animal management conditions, or of parasitic stages to exposure to anthelmintic drugs; potentially allowing for greater pathogenicity, or the development of anthelmintic resistance, respectively.