Search

Mass-drug administration (MDA) of human populations using praziquantel monotherapy has become the primary strategy for controlling and potentially eliminating the major neglected tropical disease schistosomiasis. To understand how long-term MDA impacts schistosome populations, we analysed whole-genome sequence data of 570 Schistosoma mansoni samples (and the closely related outgroup species, S. rodhaini) from eight countries incorporating both publicly-available sequence data and new parasite material. This revealed broad-scale genetic structure across countries but with extensive transmission over hundreds of kilometres. We characterised variation across the transient receptor potential melastatin ion channel, TRPM PZQ , a target of praziquantel, which has recently been found to influence praziquantel susceptibility. Functional profiling of TRPM PZQ variants found in endemic populations identified four mutations that reduced channel sensitivity to praziquantel, indicating standing variation for resistance. Analysis of parasite infrapopulations sampled from individuals pre- and post-treatment identified instances of treatment failure, further indicative of potential praziquantel resistance. As schistosomiasis is targeted for elimination as a public health problem by 2030 in all currently endemic countries, and even interruption of transmission in selected African regions, we provide an in-depth genomic characterisation of endemic populations and an approach to identify emerging praziquantel resistance alleles.

Schistosoma mansoni is a parasitic flatworm that causes the major neglected tropical disease schistosomiasis. The miracidium is the first larval stage of the life cycle. It swims and infects a freshwater snail, transforms into a mother sporocyst, where its stem cells generate daughter sporocysts that give rise to human-infective cercariae larvae. To understand the miracidium at cellular and molecular levels, we created a whole-body atlas of its ~365 cells. Single-cell RNA sequencing identified 19 transcriptionally distinct cell clusters. In situ hybridisation of tissue-specific genes revealed that 93% of the cells in the larva are somatic (57% neural, 19% muscle, 13% epidermal or tegument, 2% parenchyma, and 2% protonephridia) and 7% are stem. Whereas neurons represent the most diverse somatic cell types, trajectory analysis of the two main stem cell populations indicates that one of them is the origin of the tegument lineage and the other likely contains pluripotent cells. Furthermore, unlike the somatic cells, each of these stem populations shows sex-biased transcriptional signatures suggesting a cell-type-specific gene dosage compensation for sex chromosome-linked loci. The miracidium represents a simple developmental stage with which to gain a fundamental understanding of the molecular biology and spatial architecture of schistosome cells., Competing Interests: TA, KR, CD, KA, GS, JG, JC, SD, GR, MB No competing interests declared, (© 2024, Attenborough, Rawlinson et al.)

Recently, a S168T variant in the acetylcholine receptor subunit ACR-8 was associated with levamisole resistance in the parasitic helminth Haemonchus contortus. Here, we used the Xenopus laevis oocyte expression system and two-electrode voltage-clamp electrophysiology to measure the functional impact of this S168T variant on the H. contortus levamisole-sensitive acetylcholine receptor, L-AChR-1.1. Expression of the ACR-8 S168T variant significantly reduced the current amplitude elicited by levamisole compared to acetylcholine, with levamisole changing from a full to partial agonist on the recombinant L-AChR. Functional validation of the S168T mutation on modulating levamisole activity at the receptor level highlights its critical importance as both a mechanism and a marker of levamisole resistance., Competing Interests: Declaration of competing interest The authors report no conflict of interest for this work., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Haemonchus contortus is a parasitic haematophagous nematode that primarily affects small ruminants and causes significant economic loss to the global livestock industry. Treatment of haemonchosis typically relies on broad-spectrum anthelmintics, resistance to which is an important cause of treatment failure. Resistance to levamisole remains less widespread than to other major anthelmintic classes, prompting the need for more effective and accurate surveillance to maintain its efficacy. Loop-primer endonuclease cleavage loop-mediated isothermal amplification (LEC-LAMP) is a recently developed diagnostic method that facilitates multiplex target detection with single nucleotide polymorphism (SNP) specificity and portable onsite testing. In this study, we designed a new LEC-LAMP assay and applied it to detect the levamisole resistance marker S168T in H. contortus. We explored multiplexing probes for both the resistant S168T and the susceptible S168 alleles in a single-tube assay. We then included a generic probe to detect the acr-8 gene in the multiplex assay, which could facilitate the quantification of both resistance markers and overall genetic material from H. contortus in a single step. Our results showed promising application of these technologies, demonstrating a proof-of-concept assay which is amenable to detection of resistance alleles within the parasite population, with the potential for multiplex detection, and point-of-care application enabled by lateral flow end-point detection. However, further optimisation and validation is necessary., (Copyright © 2023. Published by Elsevier Ltd.)

The flatworm Schistosoma mansoni is an important but neglected pathogen that causes the disease schistosomiasis in millions of people worldwide. The parasite has a complex life cycle, undergoing sexual reproduction in a mammalian host and asexual replication in a snail host. Understanding the molecular mechanisms that the parasite uses to transition between hosts and develop into dimorphic reproductively competent adults may reveal new strategies for control. We present the first comprehensive transcriptomic analysis of S. mansoni, from eggs to sexually naïve worms. Focusing on eight life stages spanning free-living water-borne and parasitic stages from both intermediate and definitive hosts, we have generated deep RNA-seq data for five replicates per group for a total of 75 data sets. The data were produced using a single approach to increase the accuracy of stage-to-stage comparisons and made accessible via a user-friendly tool to visualise and explore gene expression ( https://lifecycle.schisto.xyz/ ). These data are valuable for understanding the biology and sex-specific development of schistosomes and the interpretation of complementary genomic and functional genetics studies., (© 2023. The Author(s).)

Our knowledge of the diet of wild octopus paralarvae, Octopus vulgaris, is restricted to the first 2 weeks of its planktonic phase when they are selective hunters found near the coastline. These small paralarvae, bearing only three suckers per arm, are transported by oceanic currents from the coast towards offshore waters, where they complete the planktonic phase over 2 months. Here, we have investigated the trophic ecology of O. vulgaris paralarvae in two contrasting upwelling sub-regions of the Iberian Canary current (ICC) eastern boundary upwelling system and have evaluated dietary change as paralarvae develop (inferred by counting the number of suckers per arm, ranging from three to 15) along the coastal-oceanic gradient during their planktonic phase. Using high-throughput amplicon sequencing, we have characterised the diet of 100 paralarvae collected along the Northwest Iberian Peninsula (n = 65, three to five suckers per arm) and off the west coast of Morocco (n = 35, three to 15 suckers per arm), identifying up to 87 different prey species. The diet of paralarvae varied along the ICC, with crabs (53.4%), siphonophores (12.2%), copepods (12.3%), cnidarians (8.4%) and pteropods (3.7%) accounting for 90% of the variability detected off NW Iberian Peninsula, whereas off W Morocco, crabs (46.2%), copepods (23.1%), cnidarians (12.9%), krill (9.3%) and fishes (4.2%) explained 95.6% of the variability observed using frequency of observance (FOO%) data. Ontogenetic changes in the diet based on groups of paralarvae with similar numbers per arm were evidenced by the decreasing contribution of coastal meroplankton and an increase in oceanic holoplankton, including siphonophores, copepods, pteropods and krill. Trophic niche breadth values ranged from 0.06 to 0.67, with averaged values ranging from 0.23 to 0.33 (generalist = 1 and specialist = 0), suggesting that O. vulgaris paralarvae are selective predators through their ontogenetic transition between coastal and oceanic environments., (© 2023. The Author(s).)

Haemonchus contortus is a haematophagous parasitic nematode that infects small ruminants and causes significant animal health concerns and economic losses within the livestock industry on a global scale. Treatment primarily depends on broad-spectrum anthelmintics, however, resistance is established or rapidly emerging against all major drug classes. Levamisole (LEV) remains an important treatment option for parasite control, as resistance to LEV is less prevalent than to members of other major classes of anthelmintics. LEV is an acetylcholine receptor (AChR) agonist that, when bound, results in paralysis of the worm. Numerous studies implicated the AChR sub-unit, ACR-8, in LEV sensitivity and in particular, the presence of a truncated acr-8 transcript or a deletion in the acr-8 locus in some resistant isolates. Recently, a single non-synonymous SNP in acr-8 conferring a serine-to-threonine substitution (S168T) was identified that was strongly associated with LEV resistance. Here, we investigate the role of genetic variation at the acr-8 locus in a controlled genetic cross between the LEV susceptible MHco3(ISE) and LEV resistant MHco18(UGA2004) isolates of H. contortus. Using single worm PCR assays, we found that the presence of S168T was strongly associated with LEV resistance in the parental isolates and F3 progeny of the genetic cross surviving LEV treatment. We developed and optimised an allele-specific PCR assay for the detection of S168T and validated the assay using laboratory isolates and field samples that were phenotyped for LEV resistance. In the LEV-resistant field population, a high proportion (>75%) of L 3 encoded the S168T variant, whereas the variant was absent in the susceptible isolates studied. These data further support the potential role of acr-8 S168T in LEV resistance, with the allele-specific PCR providing an important step towards establishing a sensitive molecular diagnostic test for LEV resistance., Competing Interests: Declaration of competing interest The authors report that they have no conflict of interest., (Copyright © 2022 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Like other pathogens, parasitic helminths can rapidly evolve resistance to drug treatment. Understanding the genetic basis of anthelmintic drug resistance in parasitic nematodes is key to tracking its spread and improving the efficacy and sustainability of parasite control. Here, we use an in vivo genetic cross between drug-susceptible and multi-drug-resistant strains of Haemonchus contortus in a natural host-parasite system to simultaneously map resistance loci for the three major classes of anthelmintics. This approach identifies new alleles for resistance to benzimidazoles and levamisole and implicates the transcription factor cky-1 in ivermectin resistance. This gene is within a locus under selection in ivermectin-resistant populations worldwide; expression analyses and functional validation using knockdown experiments support that cky-1 is associated with ivermectin survival. Our work demonstrates the feasibility of high-resolution forward genetics in a parasitic nematode and identifies variants for the development of molecular diagnostics to combat drug resistance in the field., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

Populations within schistosomiasis control areas, especially those in Africa, are recommended to receive regular mass drug administration (MDA) with praziquantel (PZQ) as the main strategy for controlling the disease. The impact of PZQ treatment on schistosome genetics remains poorly understood, and is limited by a lack of high-resolution genetic data on the population structure of parasites within these control areas. We generated whole-genome sequence data from 174 individual miracidia collected from both children and adults from fishing communities on islands in Lake Victoria in Uganda that had received either annual or quarterly MDA with PZQ over four years, including samples collected immediately before and four weeks after treatment. Genome variation within and between samples was characterised and we investigated genomic signatures of natural selection acting on these populations that could be due to PZQ treatment. The parasite population on these islands was more diverse than found in nearby villages on the lake shore. We saw little or no genetic differentiation between villages, or between the groups of villages with different treatment intensity, but slightly higher genetic diversity within the pre-treatment compared to post-treatment parasite populations. We identified classes of genes significantly enriched within regions of the genome with evidence of recent positive selection among post-treatment and intensively treated parasite populations. The differential selection observed in post-treatment and pre-treatment parasite populations could be linked to any reduced susceptibility of parasites to praziquantel treatment., Competing Interests: The authors have declared that no competing interests exist.

The neglected tropical disease trichuriasis is caused by the whipworm Trichuris trichiura, a soil-transmitted helminth that has infected humans for millennia. Today, T. trichiura infects as many as 500 million people, predominantly in communities with poor sanitary infrastructure enabling sustained faecal-oral transmission. Using whole-genome sequencing of geographically distributed worms collected from human and other primate hosts, together with ancient samples preserved in archaeologically-defined latrines and deposits dated up to one thousand years old, we present the first population genomics study of T. trichiura. We describe the continent-scale genetic structure between whipworms infecting humans and baboons relative to those infecting other primates. Admixture and population demographic analyses support a stepwise distribution of genetic variation that is highest in Uganda, consistent with an African origin and subsequent translocation with human migration. Finally, genome-wide analyses between human samples and between human and non-human primate samples reveal local regions of genetic differentiation between geographically distinct populations. These data provide insight into zoonotic reservoirs of human-infective T. trichiura and will support future efforts toward the implementation of genomic epidemiology of this globally important helminth., (© 2022. The Author(s).)

The antiparasitic drug ivermectin plays an essential role in human and animal health globally. However, ivermectin resistance is widespread in veterinary helminths and there are growing concerns of sub-optimal responses to treatment in related helminths of humans. Despite decades of research, the genetic mechanisms underlying ivermectin resistance are poorly understood in parasitic helminths. This reflects significant uncertainty regarding the mode of action of ivermectin in parasitic helminths, and the genetic complexity of these organisms; parasitic helminths have large, rapidly evolving genomes and differences in evolutionary history and genetic background can confound comparisons between resistant and susceptible populations. We undertook a controlled genetic cross of a multi-drug resistant and a susceptible reference isolate of Haemonchus contortus, an economically important gastrointestinal nematode of sheep, and ivermectin-selected the F2 population for comparison with an untreated F2 control. RNA-seq analyses of male and female adults of all populations identified high transcriptomic differentiation between parental isolates, which was significantly reduced in the F2, allowing differences associated specifically with ivermectin resistance to be identified. In all resistant populations, there was constitutive upregulation of a single gene, HCON_00155390:cky-1, a putative pharyngeal-expressed transcription factor, in a narrow locus on chromosome V previously shown to be under ivermectin selection. In addition, we detected sex-specific differences in gene expression between resistant and susceptible populations, including constitutive upregulation of a P-glycoprotein, HCON_00162780:pgp-11, in resistant males only. After ivermectin selection, we identified differential expression of genes with roles in neuronal function and chloride homeostasis, which is consistent with an adaptive response to ivermectin-induced hyperpolarisation of neuromuscular cells. Overall, we show the utility of a genetic cross to identify differences in gene expression that are specific to ivermectin selection and provide a framework to better understand ivermectin resistance and response to treatment in parasitic helminths., Competing Interests: The authors have declared that no competing interests exist.

Background: Soil-transmitted helminths (STH) are targeted for control through mass drug-administration campaigns to prevent morbidity affecting at-risk groups in endemic regions. Although broadly successful, the use of albendazole and mebendazole achieved variable progress, with deficiencies against Trichuris trichiura and a predictable low efficacy against Strongyloides stercoralis . Novel drug combinations offer a potential solution, providing they can be delivered safely and maintain efficacy against all STH species. Here we present the protocol of a clinical trial to evaluate a fixed-dose combination (FDC) tablet containing albendazole and ivermectin that will be compared against albendazole against STH . Methods: An adaptive phase II/III randomized controlled trial will be undertaken in STH endemic sites in Ethiopia, Kenya and Mozambique to evaluate an oral FDC of 400 mg albendazole and either 9- or 18 mg ivermectin. FDC will be administered as a single dose or single doses over three-consecutive days and assessed against a single dose of 400 mg albendazole. In the phase II trial, 126 T. trichiura -infected children weighting 15 to 45 kg will be treated in a dose-escalation manner to determine safety objectives. In the phase III trial, 1097 participants aged 5 to 18 years old infected with T. trichiura, hookworm and S. stercoralis will be recruited to determine safety and efficacy. The trial will be open-label with blinded outcome assessors. Cure rate measured 21-days after-treatment in duplicate Kato-Katz is the primary efficacy outcome. Secondary objectives include efficacy evaluation by quantitative polymerase chain reaction (PCR) as an outcome measurement, description of pharmacokinetic parameters, palatability and acceptability evaluations, and monitoring of anthelmintic resistance. Conclusions: This trial with registrational goals seeks to evaluate an innovative fixed-dose combination of albendazole and ivermectin co-formulated tablets, with the goal of providing an anthelmintic regimen with improved efficacy and spectrum of coverage against STH. ClinicalTrials.gov registration: NCT05124691 (18/11/2021)., Competing Interests: Competing interests: Jaime Algorta is employee of Laboratorios Liconsa SA, an organization of the STOP consortium., (Copyright: © 2022 Krolewiecki A et al.)

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 F 1 crosses. Egg hatch tests revealed a stark difference between the two F 1 crosses with ED 50 estimates for MHco18/3 being 9 times greater than those for MHco3/18. Resistance factors based on ED 50 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.F 1 and MHco18/3.F 1 both had similar frequencies of the F200Y resistant allele (45.3% and 44.3%, respectively), whereas for F167Y, MHco18/3.F 1 had a two-fold greater frequency of the resistant-allele compared to MHco3/18.F 1 (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 (r c  = 0.97), but were less comparable for codon 167 (r c  = 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., (Copyright © 2021 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Haemonchus contortus is a pathogenic gastrointestinal nematode of small ruminants and, in part due to its capacity to develop resistance to drugs, contributes to significant losses in the animal production sector worldwide. Despite decades of research, comparatively little is known about the specific mechanism(s) driving resistance to drugs such as ivermectin in this species. Here we describe a genome-wide approach to detect evidence of selection by ivermectin treatment in a field population of H. contortus from Sweden, using parasites sampled from the same animals before and seven days after ivermectin exposure followed by whole-genome sequencing. Despite an 89% reduction in parasites recovered after treatment measured by the fecal egg count reduction test, the surviving population was highly genetically similar to the population before treatment, suggesting that resistance has likely evolved over time and that resistance alleles are present on diverse haplotypes. Pairwise gene and SNP frequency comparisons indicated the highest degree of differentiation was found at the terminal end of chromosome 4, whereas the most striking difference in nucleotide diversity was observed in a region on chromosome 5 previously reported to harbor a major quantitative trait locus involved in ivermectin resistance. These data provide novel insight into the genome-wide effect of ivermectin selection in a field population as well as confirm the importance of the previously established quantitative trait locus in the development of resistance to ivermectin., (Copyright © 2021 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Gastrointestinal (GI) helminth infections cause significant morbidity in both humans and animals worldwide. Specific and sensitive diagnosis is central to the surveillance of such infections and to determine the effectiveness of treatment strategies used to control them. In this article, we: (i) assess the strengths and limitations of existing methods applied to the diagnosis of GI helminth infections of humans and livestock; (ii) examine high-throughput sequencing approaches, such as targeted molecular barcoding and shotgun sequencing, as tools to define the taxonomic composition of helminth infections; and (iii) discuss the current understanding of the interactions between helminths and microbiota in the host gut. Stool-based diagnostics are likely to serve as an important tool well into the future; improved diagnostics of helminths and their environment in the gut may assist the identification of biomarkers with the potential to define the health/disease status of individuals and populations, and to identify existing or emerging anthelmintic resistance., (© 2022. The Author(s).)

We present a genome assembly and annotation of an individual female Cercopithifilaria johnstoni , a parasitic filarial nematode that is transmitted by hard ticks (Ixodidae) to infect a broad range of native Australian murid and marsupial hosts. The genome sequence is 76.9 Mbp in length, and although in draft form (N50 = 99 kbp, N50[n] = 232), is largely complete based on universally conserved orthologs (BUSCOs; genome = 94.9%, protein = 96.5%) and relative to other related filarial species. These data represent the first genomic resources for the genus Cercopithifilaria , a group of parasites with a broad host range, and form the basis for comparative analysis with the human-infective parasite, Onchocerca volvulus , both of which are responsible for similar eye and skin pathologies in their respective hosts., Competing Interests: No competing interests were disclosed., (Copyright: © 2021 McCann K et al.)

Over 250 million people suffer from schistosomiasis, a tropical disease caused by parasitic flatworms known as schistosomes. Humans become infected by free-swimming, water-borne larvae, which penetrate the skin. The earliest intra-mammalian stage, called the schistosomulum, undergoes a series of developmental transitions. These changes are critical for the parasite to adapt to its new environment as it navigates through host tissues to reach its niche, where it will grow to reproductive maturity. Unravelling the mechanisms that drive intra-mammalian development requires knowledge of the spatial organisation and transcriptional dynamics of different cell types that comprise the schistomulum body. To fill these important knowledge gaps, we perform single-cell RNA sequencing on two-day old schistosomula of Schistosoma mansoni. We identify likely gene expression profiles for muscle, nervous system, tegument, oesophageal gland, parenchymal/primordial gut cells, and stem cells. In addition, we validate cell markers for all these clusters by in situ hybridisation in schistosomula and adult parasites. Taken together, this study provides a comprehensive cell-type atlas for the early intra-mammalian stage of this devastating metazoan parasite.

Anthelmintic resistance is a significant threat to livestock production systems worldwide and is emerging as an important issue in companion animal parasite management. It is also an emerging concern for the control of human soil-transmitted helminths and filaria. An important aspect of managing anthelmintic resistance is the ability to utilise diagnostic tests to detect its emergence at an early stage. In host-parasite systems where resistance is already widespread, diagnostics have a potentially important role in determining those drugs that remain the most effective. The development of molecular diagnostics for anthelmintic resistance is one focus of the Consortium for Anthelmintic Resistance and Susceptibility (CARS) group. The present paper reflects discussions of this issue that occurred at the most recent meeting of the group in Wisconsin, USA, in July 2019. We compare molecular resistance diagnostics with in vivo and in vitro phenotypic methods, and highlight the advantages and disadvantages of each. We assess whether our knowledge on the identity of molecular markers for resistance towards the different drug classes is sufficient to provide some expectation that molecular tests for field use may be available in the short-to-medium term. We describe some practical aspects of such tests and how our current capabilities compare to the requirements of an 'ideal' test. Finally, we describe examples of drug class/parasite species interactions that provide the best opportunity for commercial use of molecular tests in the near future. We argue that while such prototype tests may not satisfy the requirements of an 'ideal' test, their potential to provide significant advances over currently-used phenotypic methods warrants their development as field diagnostics., (Copyright © 2020 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Background: Guinea worm-Dracunculus medinensis-was historically one of the major parasites of humans and has been known since antiquity. Now, Guinea worm is on the brink of eradication, as efforts to interrupt transmission have reduced the annual burden of disease from millions of infections per year in the 1980s to only 54 human cases reported globally in 2019. Despite the enormous success of eradication efforts to date, one complication has arisen. Over the last few years, hundreds of dogs have been found infected with this previously apparently anthroponotic parasite, almost all in Chad. Moreover, the relative numbers of infections in humans and dogs suggests that dogs are currently the principal reservoir on infection and key to maintaining transmission in that country., Principal Findings: In an effort to shed light on this peculiar epidemiology of Guinea worm in Chad, we have sequenced and compared the genomes of worms from dog, human and other animal infections. Confirming previous work with other molecular markers, we show that all of these worms are D. medinensis, and that the same population of worms are causing both infections, can confirm the suspected transmission between host species and detect signs of a population bottleneck due to the eradication efforts. The diversity of worms in Chad appears to exclude the possibility that there were no, or very few, worms present in the country during a 10-year absence of reported cases., Conclusions: This work reinforces the importance of adequate surveillance of both human and dog populations in the Guinea worm eradication campaign and suggests that control programs aiming to interrupt disease transmission should stay aware of the possible emergence of unusual epidemiology as pathogens approach elimination., Competing Interests: The authors have declared that no competing interests exist. Author Ouakou Tchindebet was unable to confirm their authorship contributions. On their behalf, the corresponding author has reported their contributions to the best of their knowledge.

Haemonchus contortus is a globally distributed and economically important gastrointestinal pathogen of small ruminants and has become a key nematode model for studying anthelmintic resistance and other parasite-specific traits among a wider group of parasites including major human pathogens. Here, we report using PacBio long-read and OpGen and 10X Genomics long-molecule methods to generate a highly contiguous 283.4 Mbp chromosome-scale genome assembly including a resolved sex chromosome for the MHco3(ISE).N1 isolate. We show a remarkable pattern of conservation of chromosome content with Caenorhabditis elegans, but almost no conservation of gene order. Short and long-read transcriptome sequencing allowed us to define coordinated transcriptional regulation throughout the parasite's life cycle and refine our understanding of cis- and trans-splicing. Finally, we provide a comprehensive picture of chromosome-wide genetic diversity both within a single isolate and globally. These data provide a high-quality comparison for understanding the evolution and genomics of Caenorhabditis and other nematodes and extend the experimental tractability of this model parasitic nematode in understanding helminth biology, drug discovery and vaccine development, as well as important adaptive traits such as drug resistance.

Background: Chromosome-level assemblies are indispensable for accurate gene prediction, synteny assessment, and understanding higher-order genome architecture. Reference and draft genomes of key helminth species have been published, but little is yet known about the biology of their chromosomes. Here, we present the complete genome of the tapeworm Hymenolepis microstoma, providing a reference quality, end-to-end assembly that represents the first fully assembled genome of a spiralian/lophotrochozoan, revealing new insights into chromosome evolution., Results: Long-read sequencing and optical mapping data were added to previous short-read data enabling complete re-assembly into six chromosomes, consistent with karyology. Small genome size (169 Mb) and lack of haploid variation (1 SNP/3.2 Mb) contributed to exceptionally high contiguity with only 85 gaps remaining in regions of low complexity sequence. Resolution of repeat regions reveals novel gene expansions, micro-exon genes, and spliced leader trans-splicing, and illuminates the landscape of transposable elements, explaining observed length differences in sister chromatids. Syntenic comparison with other parasitic flatworms shows conserved ancestral linkage groups indicating that the H. microstoma karyotype evolved through fusion events. Strikingly, the assembly reveals that the chromosomes terminate in centromeric arrays, indicating that these motifs play a role not only in segregation, but also in protecting the linear integrity and full lengths of chromosomes., Conclusions: Despite strong conservation of canonical telomeres, our results show that they can be substituted by more complex, species-specific sequences, as represented by centromeres. The assembly provides a robust platform for investigations that require complete genome representation.

Haemonchus contortus is a haematophagous parasitic nematode of veterinary interest. We have performed a survey of its genome-wide diversity using single-worm whole genome sequencing of 223 individuals sampled from 19 isolates spanning five continents. We find an African origin for the species, together with evidence for parasites spreading during the transatlantic slave trade and colonisation of Australia. Strong selective sweeps surrounding the β-tubulin locus, a target of benzimidazole anthelmintic drug, are identified in independent populations. These sweeps are further supported by signals of diversifying selection enriched in genes involved in response to drugs and other anthelmintic-associated biological functions. We also identify some candidate genes that may play a role in ivermectin resistance. Finally, genetic signatures of climate-driven adaptation are described, revealing a gene acting as an epigenetic regulator and components of the dauer pathway. These results begin to define genetic adaptation to climate in a parasitic nematode.

Plant defensins are a large family of proteins, most of which have antifungal activity against a broad spectrum of fungi. However, little is known about how they exert their activity. The mechanisms of action of only a few members of the family have been investigated and, in most cases, there are still a number of unknowns. To gain a better understanding of the antifungal mechanisms of a set of four defensins, NaD1, DmAMP1, NbD6, and SBI6, we screened a pooled collection of the nonessential gene deletion set of Saccharomyces cerevisiae Strains with increased or decreased ability to survive defensin treatment were identified based on the relative abundance of the strain-specific barcode as determined by MiSeq next-generation sequencing. Analysis of the functions of genes that are deleted in strains with differential growth in the presence of defensin provides insight into the mechanism of action. The screen identified a novel role for the vacuole in the mechanisms of action for defensins NbD6 and SBI6. The effect of these defensins on vacuoles was further confirmed by using confocal microscopy in both S. cerevisiae and the cereal pathogen Fusarium graminearum These results demonstrate the utility of this screening method to identify novel mechanisms of action for plant defensins., (Copyright © 2019 American Society for Microbiology.)

Whole-genome sequencing is being rapidly applied to the study of helminth genomes, including de novo genome assembly, population genetics, and diagnostic applications. Although late-stage juvenile and adult parasites typically produce sufficient DNA for molecular analyses, these parasitic stages are almost always inaccessible in the live host; immature life stages found in the environment for which samples can be collected non-invasively offer a potential alternative; however, these samples typically yield very low quantities of DNA, can be environmentally resistant, and are susceptible to contamination, often from bacterial or host DNA. Here, we have tested five low-input DNA extraction protocols together with a low-input sequencing library protocol to assess the feasibility of whole-genome sequencing of individual immature helminth samples. These approaches do not use whole-genome amplification, a common but costly approach to increase the yield of low-input samples. We first tested individual parasites from two species spotted onto FTA cards-egg and L1 stages of Haemonchus contortus and miracidia of Schistosoma mansoni -before further testing on an additional five species- Ancylostoma caninum , Ascaridia dissimilis , Dirofilaria immitis , Strongyloides stercoralis , and Trichuris muris -with an optimal protocol. A sixth species- Dracunculus medinensis -was included for comparison. Whole-genome sequencing followed by analyses to determine the proportion of on- and off-target mapping revealed successful sample preparations for six of the eight species tested with variation both between species and between different life stages from some species described. These results demonstrate the feasibility of whole-genome sequencing of individual parasites, and highlight a new avenue toward generating sensitive, specific, and information-rich data for the diagnosis and surveillance of helminths., (Copyright © 2019 Doyle, Sankaranarayanan, Allan, Berger, Jimenez Castro, Collins, Crellen, Duque-Correa, Ellis, Jaleta, Laing, Maitland, McCarthy, Moundai, Softley, Thiele, Ouakou, Tushabe, Webster, Weiss, Lok, Devaney, Kaplan, Cotton, Berriman and Holroyd.)

Anthelmintic resistance is a threat to global food security. In order to alleviate the selection pressure for resistance and maintain drug efficacy, management strategies increasingly aim to preserve a proportion of the parasite population in 'refugia', unexposed to treatment. While persuasive in its logic, and widely advocated as best practice, evidence for the ability of refugia-based approaches to slow the development of drug resistance in parasitic helminths is currently limited. Moreover, the conditions needed for refugia to work, or how transferable those are between parasite-host systems, are not known. This review, born of an international workshop, seeks to deconstruct the concept of refugia and examine its assumptions and applicability in different situations. We conclude that factors potentially important to refugia, such as the fitness cost of drug resistance, the degree of mixing between parasite sub-populations selected through treatment or not, and the impact of parasite life-history, genetics and environment on the population dynamics of resistance, vary widely between systems. The success of attempts to generate refugia to limit anthelmintic drug resistance are therefore likely to be highly dependent on the system in hand. Additional research is needed on the concept of refugia and the underlying principles for its application across systems, as well as empirical studies within systems that prove and optimise its usefulness., (Copyright © 2019 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Background: Infections with helminths cause an enormous disease burden in billions of animals and plants worldwide. Large scale use of anthelmintics has driven the evolution of resistance in a number of species that infect livestock and companion animals, and there are growing concerns regarding the reduced efficacy in some human-infective helminths. Understanding the mechanisms by which resistance evolves is the focus of increasing interest; robust genetic analysis of helminths is challenging, and although many candidate genes have been proposed, the genetic basis of resistance remains poorly resolved., Results: Here, we present a genome-wide analysis of two genetic crosses between ivermectin resistant and sensitive isolates of the parasitic nematode Haemonchus contortus, an economically important gastrointestinal parasite of small ruminants and a model for anthelmintic research. Whole genome sequencing of parental populations, and key stages throughout the crosses, identified extensive genomic diversity that differentiates populations, but after backcrossing and selection, a single genomic quantitative trait locus (QTL) localised on chromosome V was revealed to be associated with ivermectin resistance. This QTL was common between the two geographically and genetically divergent resistant populations and did not include any leading candidate genes, suggestive of a previously uncharacterised mechanism and/or driver of resistance. Despite limited resolution due to low recombination in this region, population genetic analyses and novel evolutionary models supported strong selection at this QTL, driven by at least partial dominance of the resistant allele, and that large resistance-associated haplotype blocks were enriched in response to selection., Conclusions: We have described the genetic architecture and mode of ivermectin selection, revealing a major genomic locus associated with ivermectin resistance, the most conclusive evidence to date in any parasitic nematode. This study highlights a novel genome-wide approach to the analysis of a genetic cross in non-model organisms with extreme genetic diversity, and the importance of a high-quality reference genome in interpreting the signals of selection so identified.

The replacement of forest areas with human-dominated landscapes usually leads to fragmentation, altering the structure and function of the forest. Here we studied the dynamics of forest patch sizes at a global level, examining signals of a critical transition from an unfragmented to a fragmented state, using the MODIS vegetation continuous field. We defined wide regions of connected forest across continents and big islands, and combined five criteria, including the distribution of patch sizes and the fluctuations of the largest patch over the last sixteen years, to evaluate the closeness of each region to a fragmentation threshold. Regions with the highest deforestation rates-South America, Southeast Asia, Africa-all met these criteria and may thus be near a critical fragmentation threshold. This implies that if current forest loss rates are maintained, wide continental areas could suddenly fragment, triggering extensive species loss and degradation of ecosystems services.

Antarctica is seriously affected by climate change, particularly at the Western Antarctic Peninsula (WAP) where a rapid regional warming is observed. Potter Cove is a WAP fjord at Shetland Islands that constitutes a biodiversity hotspot where over the last years, Potter Cove annual air temperatures averages increased by 0.66 °C, coastal glaciers declined, and suspended particulate matter increased due to ice melting. Macroalgae are the main energy source for all consumers and detritivores of Potter Cove. Some effects of climate change favor pioneer macroalgae species that exploit new ice-free areas and can also decline rates of photosynthesis and intensify competition between species due to the increase of suspended particulate matter. In this study, we evaluated possible consequences of climate change at Potter Cove food web by simulating the extinction of macroalgae and detritus using a topological approach with thresholds of extinction. Thresholds represent the minimum number of incoming links necessary for species' survival. When we simulated the extinctions of macroalgae species at random, a threshold of extinction beyond 50% was necessary to obtain a significant number of secondary extinctions, while with a 75% threshold a real collapse of the food web occurred. Our results indicate that Potter Cove food web is relative robust to macroalgae extinction. This is dramatically different from what has been found in other food webs, where the reduction of 10% in prey intake caused a disproportionate increase of secondary extinctions. Robustness of the Potter Cove food web was mediated by omnivory and redundancy, which had an important relevance in this food web. When we eliminated larger-biomass species more secondary extinctions occurred, a similar response was observed when more connected species were deleted, yet there was no correlation between species of larger-biomass and high-degree. This similarity could be explained because both criteria involved key species that produced an emerging effect on the food web. In this way, large-biomass and high-degree species could be acting as source for species with few trophic interactions or low redundancy. Based on this work, we expect the Potter Cove food web to be robust to changes in macroalgae species caused by climate change until a high threshold of stress is reached, and then negative effects are expected to spread through the entire food web leading to its collapse., Competing Interests: The authors declare there are no competing interests.

A 'sibling' species of the model organism Caenorhabditis elegans has long been sought for use in comparative analyses that would enable deep evolutionary interpretations of biological phenomena. Here, we describe the first sibling species of C. elegans, C. inopinata n. sp., isolated from fig syconia in Okinawa, Japan. We investigate the morphology, developmental processes and behaviour of C. inopinata, which differ significantly from those of C. elegans. The 123-Mb C. inopinata genome was sequenced and assembled into six nuclear chromosomes, allowing delineation of Caenorhabditis genome evolution and revealing unique characteristics, such as highly expanded transposable elements that might have contributed to the genome evolution of C. inopinata. In addition, C. inopinata exhibits massive gene losses in chemoreceptor gene families, which could be correlated with its limited habitat area. We have developed genetic and molecular techniques for C. inopinata; thus C. inopinata provides an exciting new platform for comparative evolutionary studies.

The search for general properties in network structure has been a central issue for food web studies in recent years. One such property is the small-world topology that combines a high clustering and a small distance between nodes of the network. This property may increase food web resilience but make them more sensitive to the extinction of connected species. Food web theory has been developed principally from freshwater and terrestrial ecosystems, largely omitting marine habitats. If theory needs to be modified to accommodate observations from marine ecosystems, based on major differences in several topological characteristics is still on debate. Here we investigated if the small-world topology is a common structural pattern in marine food webs. We developed a novel, simple and statistically rigorous method to examine the largest set of complex marine food webs to date. More than half of the analyzed marine networks exhibited a similar or lower characteristic path length than the random expectation, whereas 39% of the webs presented a significantly higher clustering than its random counterpart. Our method proved that 5 out of 28 networks fulfilled both features of the small-world topology: short path length and high clustering. This work represents the first rigorous analysis of the small-world topology and its associated features in high-quality marine networks. We conclude that such topology is a structural pattern that is not maximized in marine food webs; thus it is probably not an effective model to study robustness, stability and feasibility of marine ecosystems., Competing Interests: The authors have declared that no competing interests exist.

The parasitic nematode Haemonchus contortus is an economically and clinically important pathogen of small ruminants, and a model system for understanding the mechanisms and evolution of traits such as anthelmintic resistance. Anthelmintic resistance is widespread and is a major threat to the sustainability of livestock agriculture globally; however, little is known about the genome architecture and parameters such as recombination that will ultimately influence the rate at which resistance may evolve and spread. Here, we performed a genetic cross between two divergent strains of H. contortus, and subsequently used whole-genome resequencing of a female worm and her brood to identify the distribution of genome-wide variation that characterizes these strains. Using a novel bioinformatic approach to identify variants that segregate as expected in a pseudotestcross, we characterized linkage groups and estimated genetic distances between markers to generate a chromosome-scale F1 genetic map. We exploited this map to reveal the recombination landscape, the first for any helminth species, demonstrating extensive variation in recombination rate within and between chromosomes. Analyses of these data also revealed the extent of polyandry, whereby at least eight males were found to have contributed to the genetic variation of the progeny analyzed. Triploid offspring were also identified, which we hypothesize are the result of nondisjunction during female meiosis or polyspermy. These results expand our knowledge of the genetics of parasitic helminths and the unusual life-history of H. contortus, and enhance ongoing efforts to understand the genetic basis of resistance to the drugs used to control these worms and for related species that infect livestock and humans throughout the world. This study also demonstrates the feasibility of using whole-genome resequencing data to directly construct a genetic map in a single generation cross from a noninbred nonmodel organism with a complex lifecycle., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background: Treatment of onchocerciasis using mass ivermectin administration has reduced morbidity and transmission throughout Africa and Central/South America. Mass drug administration is likely to exert selection pressure on parasites, and phenotypic and genetic changes in several Onchocerca volvulus populations from Cameroon and Ghana-exposed to more than a decade of regular ivermectin treatment-have raised concern that sub-optimal responses to ivermectin's anti-fecundity effect are becoming more frequent and may spread., Methodology/principal Findings: Pooled next generation sequencing (Pool-seq) was used to characterise genetic diversity within and between 108 adult female worms differing in ivermectin treatment history and response. Genome-wide analyses revealed genetic variation that significantly differentiated good responder (GR) and sub-optimal responder (SOR) parasites. These variants were not randomly distributed but clustered in ~31 quantitative trait loci (QTLs), with little overlap in putative QTL position and gene content between the two countries. Published candidate ivermectin SOR genes were largely absent in these regions; QTLs differentiating GR and SOR worms were enriched for genes in molecular pathways associated with neurotransmission, development, and stress responses. Finally, single worm genotyping demonstrated that geographic isolation and genetic change over time (in the presence of drug exposure) had a significantly greater role in shaping genetic diversity than the evolution of SOR., Conclusions/significance: This study is one of the first genome-wide association analyses in a parasitic nematode, and provides insight into the genomics of ivermectin response and population structure of O. volvulus. We argue that ivermectin response is a polygenically-determined quantitative trait (QT) whereby identical or related molecular pathways but not necessarily individual genes are likely to determine the extent of ivermectin response in different parasite populations. Furthermore, we propose that genetic drift rather than genetic selection of SOR is the underlying driver of population differentiation, which has significant implications for the emergence and potential spread of SOR within and between these parasite populations.

Deficiency in pleomorphic adenoma gene 1 (PLAG1) leads to reduced fertility in male mice, but the mechanism by which PLAG1 contributes to reproduction is unknown. To investigate the involvement of PLAG1 in testicular function, we determined (i) the spatial distribution of PLAG1 in the testis using X-gal staining; (ii) transcriptomic consequences of PLAG1 deficiency in knock-out and heterozygous mice compared to wild-type mice using RNA-seq; and (iii) morphological and functional consequences of PLAG1 deficiency by determining testicular histology, daily sperm production and sperm motility in knock-out and wild-type mice. PLAG1 was sparsely expressed in germ cells and in Sertoli cells. Genes known to be involved in spermatogenesis were downregulated in the testes of knock-out mice, as well as Hsd17b3, which encodes a key enzyme in androgen biosynthesis. In the absence of Plag1, a number of genes involved in immune processes and epididymis-specific genes were upregulated in the testes. Finally, loss of PLAG1 resulted in significantly lowered daily sperm production, in reduced sperm motility, and in several animals, in sloughing of the germinal epithelium. Our results demonstrate that the subfertility seen in male PLAG1-deficient mice is, at least in part, the result of significantly reduced sperm output and sperm motility.

Preventive chemotherapy has long been practiced against nematode parasites of livestock, leading to widespread drug resistance, and is increasingly being adopted for eradication of human parasitic nematodes even though it is similarly likely to lead to drug resistance. Given that the genetic architecture of resistance is poorly understood for any nematode, we have analyzed multidrug resistant Teladorsagia circumcincta, a major parasite of sheep, as a model for analysis of resistance selection. We introgressed a field-derived multiresistant genotype into a partially inbred susceptible genetic background (through repeated backcrossing and drug selection) and performed genome-wide scans in the backcross progeny and drug-selected F2 populations to identify the major genes responsible for the multidrug resistance. We identified variation linking candidate resistance genes to each drug class. Putative mechanisms included target site polymorphism, changes in likely regulatory regions and copy number variation in efflux transporters. This work elucidates the genetic architecture of multiple anthelmintic resistance in a parasitic nematode for the first time and establishes a framework for future studies of anthelmintic resistance in nematode parasites of humans.

The common octopus ( Octopus vulgaris ) is an attractive species for aquaculture, however, several challenges inhibit sustainable commercial production. Little is known about the early paralarval stages in the wild, including diet and intestinal microbiota, which likely play a significant role in development and vitality of this important life stage. High throughput sequencing was used to characterize the gastrointestinal microbiome of wild O. vulgaris paralarvae collected from two different upwelling regions off the coast of North West Spain ( n = 41) and Morocco ( n = 35). These were compared to that of paralarvae reared with Artemia for up to 25 days in captivity ( n = 29). In addition, the gastrointestinal microbiome of zooplankton prey (crabs, copepod and krill) was also analyzed to determine if the microbial communities present in wild paralarvae are derived from their diet. Paralarvae reared in captivity with Artemia showed a depletion of bacterial diversity, particularly after day 5, when almost half the bacterial species present on day 0 were lost and two bacterial families (Mycoplasmataceae and Vibrionaceae) dominated the microbial community. In contrast, bacterial diversity increased in wild paralarvae as they developed in the oceanic realm of both upwelling systems, likely due to the exposure of new bacterial communities via ingestion of a wide diversity of prey. Remarkably, the bacterial diversity of recently hatched paralarvae in captivity was similar to that of wild paralarvae and zooplankton, thus suggesting a marked effect of the diet in both the microbial community species diversity and evenness. This study provides a comprehensive overview of the bacterial communities inhabiting the gastrointestinal tract of O. vulgaris paralarvae, and reveals new research lines to challenge the current bottlenecks preventing sustainable octopus aquaculture.

Background: Wolbachia are intracellular bacteria found in arthropods and several filarial nematode species. The filarial Wolbachia have been proposed to be involved in the immunopathology associated with onchocerciasis. Higher Wolbachia-to-nematode ratios have been reported in the savannah-ecotype compared to the forest-ecotype, and have been interpreted as consistent with a correlation between Wolbachia density and disease severity. However, factors such as geographic stratification and ivermectin drug exposure can lead to significant genetic heterogeneity in the nematode host populations, so we investigated whether Wolbachia copy number variation is also associated with these underlying factors., Methods: Genomic DNA was prepared from single adult nematodes representing forest and savannah ecotypes sampled from Togo, Ghana, Côte d'Ivoire and Mali. A qPCR assay was developed to measure the number of Wolbachia genome(s) per nematode genome. Next-generation sequencing (NGS) was also used to measure relative Wolbachia copy number, and independently verify the qPCR assay., Results: Significant variation was observed within the forest (range: 0.02 to 452.99; median: 10.58) and savannah (range: 0.01 to 1106.25; median: 9.10) ecotypes, however, no significant difference between ecotypes (P = 0.645) was observed; rather, strongly significant Wolbachia variation was observed within and between the nine study communities analysed (P = 0.021), independent of ecotype. Analysis of ivermectin-treated and untreated nematodes by qPCR showed no correlation (P = 0.869); however, an additional analysis of a subset of the nematodes by qPCR and NGS revealed a correlation between response to ivermectin treatment and Wolbachia copy number (P = 0.020)., Conclusions: This study demonstrates that extensive within and between population variation exists in the Wolbachia content of individual adult O. volvulus. The origin and functional significance of such variation (up to ~ 100,000-fold between worms; ~10 to 100-fold between communities) in the context of the proposed mutualistic relationship between the worms and the bacteria, and between the presence of Wolbachia and clinical outcome of infection, remains unclear. These data do not support a correlation between Wolbachia copy number and forest or savannah ecotype, and may have implications for the development of anti-Wolbachia drugs as a macrofilaricidal treatment of onchocerciasis. The biological significance of a correlation between variation in Wolbachia copy number and ivermectin response remains unexplained.

Background: Genetic surveillance of the human filarial parasite, Onchocerca volvulus, from onchocerciasis endemic regions will ideally focus on genotyping individual infective larval stages collected from their intermediate host, Simuliid blackflies. However, blackflies also transmit other Onchocerca species, including the cattle parasite O. ochengi, which are difficult to distinguish from the human parasite based on morphological characteristics alone. This study describes a versatile approach to discriminate between O. volvulus and O. ochengi that is demonstrated using parasite infective larvae dissected from blackflies., Results: A speciation assay was designed based on genetic differentiation between O. volvulus and O. ochengi mitochondrial genome sequences that can be performed in high-throughput high-resolution melt (HRM)- or lower throughput conventional restriction fragment length polymorphism (RFLP) analyses. This assay was validated on 185 Onchocerca larvae dissected from blackflies captured from 14 communities in Ghana throughout 2011-2013. The frequency of O. ochengi was approximately 67 % of all larvae analysed, which is significantly higher than previously reported in this region. Furthermore, the species distribution was not uniform throughout the study region, with 25 %, 47 % and 93 % of O. volvulus being found in the western-most (Black Volta, Tain and Tombe), the central (Pru) and eastern-most (Daka) river basins, respectively., Conclusions: This tool provides a simple and cost-effective approach to determine the identity and distribution of two Onchocerca species, and will be valuable for future genetic studies that focus on parasites collected from blackflies. The results presented highlight the need to discriminate Onchocerca species in transmission studies, as the frequency of each species varied significantly between the communities studied.

Little is known about the prevalence, functionality and ecological roles of temperate phages for members of the mycolic acid producing bacteria, the Mycolata. While many lytic phages infective for these organisms have been isolated, and assessed for their suitability for use as biological control agents of activated sludge foaming, no studies have investigated how temperate phages might be induced for this purpose. Bioinformatic analysis using the PHAge Search Tool (PHAST) on Mycolata whole genome sequence data in GenBank for members of the genera Gordonia, Mycobacterium, Nocardia, Rhodococcus, and Tsukamurella revealed 83% contained putative prophage DNA sequences. Subsequent prophage inductions using mitomycin C were conducted on 17 Mycolata strains. This led to the isolation and genome characterization of three novel Caudovirales temperate phages, namely GAL1, GMA1, and TPA4, induced from Gordonia alkanivorans, Gordonia malaquae, and Tsukamurella paurometabola, respectively. All possessed highly distinctive dsDNA genome sequences.

Background: DNA methylation at a gene promoter region has the potential to regulate gene transcription. Patterns of methylation over multiple CpG sites in a region are often complex and cell type specific, with the region showing multiple allelic patterns in a sample. This complexity is commonly obscured when DNA methylation data is summarised as an average percentage value for each CpG site (or aggregated across CpG sites). True representation of methylation patterns can only be fully characterised by clonal analysis. Deep sequencing provides the ability to investigate clonal DNA methylation patterns in unprecedented detail and scale, enabling the proper characterisation of the heterogeneity of methylation patterns. However, the sheer amount and complexity of sequencing data requires new synoptic approaches to visualise the distribution of allelic patterns., Results: We have developed a new analysis and visualisation software tool "Methpat", that extracts and displays clonal DNA methylation patterns from massively parallel sequencing data aligned using Bismark. Methpat was used to analyse multiplex bisulfite amplicon sequencing on a range of CpG island targets across a panel of human cell lines and primary tissues. Methpat was able to represent the clonal diversity of epialleles analysed at specific gene promoter regions. We also used Methpat to describe epiallelic DNA methylation within the mitochondrial genome., Conclusions: Methpat can summarise and visualise epiallelic DNA methylation results from targeted amplicon, massively parallel sequencing of bisulfite converted DNA in a compact and interpretable format. Unlike currently available tools, Methpat can visualise the diversity of epiallelic DNA methylation patterns in a sample.

Background: DNA methylation is a complex epigenetic marker that can be analyzed using a wide variety of methods. Interpretation and visualization of DNA methylation data can mask complexity in terms of methylation status at each CpG site, cellular heterogeneity of samples and allelic DNA methylation patterns within a given DNA strand. Bisulfite sequencing is considered the gold standard, but visualization of massively parallel sequencing results remains a significant challenge., Findings: We created a program called Methpat that facilitates visualization and interpretation of bisulfite sequencing data generated by massively parallel sequencing. To demonstrate this, we performed multiplex PCR that targeted 48 regions of interest across 86 human samples. The regions selected included known gene promoters associated with cancer, repetitive elements, known imprinted regions and mitochondrial genomic sequences. We interrogated a range of samples including human cell lines, primary tumours and primary tissue samples. Methpat generates two forms of output: a tab-delimited text file for each sample that summarizes DNA methylation patterns and their read counts for each amplicon, and a HTML file that summarizes this data visually. Methpat can be used with publicly available whole genome bisulfite sequencing and reduced representation bisulfite sequencing datasets with sufficient read depths., Conclusions: Using Methpat, complex DNA methylation data derived from massively parallel sequencing can be summarized and visualized for biological interpretation. By accounting for allelic DNA methylation states and their abundance in a sample, Methpat can unmask the complexity of DNA methylation and yield further biological insight in existing datasets.

Mitochondrial targeting signals (MTSs) are responsible for trafficking nuclear encoded proteins to their final destination within mitochondria. These sequences are diverse, sharing little amino acid homology and vary significantly in length, and although the formation of a positively-charged amphiphilic alpha helix within the MTS is considered to be necessary and sufficient to mediate import, such a feature does not explain their diversity, nor how such diversity influences target sequence function, nor how such dissimilar signals interact with a single, evolutionarily conserved import mechanism. An in silico analysis of 296 N-terminal, matrix destined MTSs from Homo sapiens, Mus musculus, Saccharomyces cerevisiae, Arabidopsis thaliana, and Oryza sativa was undertaken to investigate relationships between MTSs, and/or, relationships between an individual targeting signal sequence and the protein that it imports. We present evidence that suggests MTS diversity is influenced in part by physiochemical and N-terminal characteristics of their mature sequences, and that some of these correlated characteristics are evolutionarily maintained across a number of taxa. Importantly, some of these associations begin to explain the variation in MTS length and composition.

Background: Many SNP discrimination strategies employ natural restriction endonucleases to discriminate between allelic states. However, SNPs are often not associated with a restriction site and therefore, a number of attempts have been made to generate sequence-adaptable restriction endonucleases. In this study, a simple, sequence-adaptable SNP discrimination mechanism between a 'wild-type' and 'mutant' template is demonstrated. This model differs from other artificial restriction endonuclease models as cis- rather than trans-orientated regions of single stranded DNA were generated and cleaved, and therefore, overcomes potential issues of either inefficient or non-specific binding when only a single variant is targeted., Results: A series of mismatch 'bubbles' that spanned 0-5-bp surrounding a point mutation was generated and analysed for sensitivity to S1 nuclease. In this model, generation of oligonucleotide-mediated ssDNA mismatch 'bubbles' in the presence of S1 nuclease resulted in the selective degradation of the mutant template while maintaining wild-type template integrity. Increasing the size of the mismatch increased the rate of mutant sequence degradation, until a threshold above which discrimination was lost and the wild-type sequence was degraded. This level of fine discrimination was possible due to the development of a novel high-resolution melting curve assay to empirically determine changes in Tm (~5.0°C per base-pair mismatch) and to optimise annealing conditions (~18.38°C below Tm) of the mismatched oligonucleotide sets., Conclusions: The in vitro 'cleavage bubble' model presented is sequence-adaptable as determined by the binding oligonucleotide, and hence, has the potential to be tailored to discriminate between any two or more SNPs. Furthermore, the demonstrated fluorometric assay has broad application potential, offering a rapid, sensitive and high-throughput means to determine Tm and annealing rates as an alternative to conventional hybridisation detection strategies.

Mitochondrial DNA (mtDNA) disorders include a vast range of pathological conditions, despite each sharing a mutual inability to produce ATP efficiently as a result of defective oxidative phosphorylation. There is no clear consensus regarding an effective therapeutic approach, and consequently the current treatment strategies are largely supportive rather than curative. This is almost certainly the result of there being virtually no defined genotype-phenotype relationships among the mtDNA disorders; hence an identical mutation may be responsible for multiple phenotypes, or the same phenotype may be produced by different mutations. In light of this, the development of gene therapy to treat mtDNA disorders offers a promising approach, as it potentially circumvents the complication of the aforementioned genotype-phenotype inconsistency and ultimately the current inability to treat individual disorders with sufficient efficacy. Such an approach will ultimately require the combination of efficient mitochondrial targeting, and an effective therapeutic molecule. Although promising proof-of-principle developments in this field have been demonstrated, the realization of a successful therapeutic mitochondrial gene therapy strategy has not come to fruition. This review critiques the key approaches under development by discussing the theory underlying each strategy, and detailing the current progress made. We also emphasize the potential hurdles that must be acknowledged and overcome if the potential of a therapeutic gene therapy to treat mitochondrial DNA disorders is to be realized.

Background: Polyethylenimine (PEI) is one of the most efficient and versatile non-viral vectors available for gene delivery. Despite many advantages over viral vectors, PEI is still limited by lower transfection efficiency compared to its viral counterparts. Considerable investigation is devoted to the modification of PEI to incorporate virus-like properties to improve its efficacy, including the incorporation of the protein transduction domain (PTD) polyarginine (Arg); itself demonstrated to facilitate membrane translocation of molecular cargo. There is, however, limited understanding of the underlying mechanisms of gene delivery facilitated by both PEI and PEI-bioconjugates such as PEI-polyarginine (PEI-Arg) within live cells, which once elucidated will provide valuable insights into the development of more efficient non-viral gene delivery vectors., Methods: PEI and PEI-Arg were investigated for their ability to facilitate DNA internalization and gene expression within live COS-7 cells, in terms of the percentage of cells transfected and the relative amount of gene expression per cell. Intracellular trafficking of vectors was investigated using fluorescent microscopy during the first 5 h post transfection. Finally, nocodazole and aphidicolin were used to investigate the role of microtubules and mitosis, respectively, and their impact on PEI and PEI-Arg mediated gene delivery and expression., Results: PEI-Arg maintained a high cellular DNA uptake efficiency, and facilitated as much as 2-fold more DNA internalization compared to PEI alone. PEI, but not PEI-Arg, displayed microtubule-facilitated trafficking, and was found to accumulate within close proximity to the nucleus. Only PEI facilitated significant gene expression, whereas PEI-Arg conferred negligible expression. Finally, while not exclusively dependent, microtubule trafficking and, to a greater extent, mitotic events significantly contributed to PEI facilitated gene expression., Conclusion: PEI polyplexes are trafficked by an indirect association with microtubules, following endosomal entrapment. PEI facilitated expression is significantly influenced by a mitotic event, which is increased by microtubule organization center (MTOC)-associated localization of PEI polyplexes. PEI-Arg, although enhancing DNA internalization per cell, did not improve gene expression, highlighting the importance of microtubule trafficking for PEI vectors and the impact of the Arg peptide to intracellular trafficking. This study emphasizes the importance of a holistic approach to investigate the mechanisms of novel gene delivery vectors.