Diversity of benzimidazole-resistance alleles in populations of small ruminant parasites

The resistance of gastro-intestinal nematodes of small ruminants (sheep and goat) to benzimidazole anthelmintic drugs seems to be linked primarily to a single mutation in the isotype 1 beta-tubulin gene. This study was carried out to investigate the origin and diversity of benzimidazole-resistance alleles in trichostrongylid nematodes. We sequenced a 550 bp fragment of the isotype 1 beta-tubulin gene from several benzimidazole-resistant Teladorsagia circumcincta populations isolated from dairy goat farms in the central and south-western France. We also sequenced the same beta-tubulin fragment from Trichostrongylus colubriformis and Haemonchus contortus populations in south-western France. We found eight benzimidazole-resistance alleles in all T. circumcincta populations studied, six in H. contortus populations, and only one in T. colubriformis populations. In most cases, only one benzimidazole-resistance allele was present in T. circumcincta and H. contortus populations, but two alleles were found in a fewer number of them. Some T. circumcincta populations shared the same benzimidazole-resistance allele whereas some others had a specific benzimidazole-resistance allele. Similar findings were obtained for H. contortus. As no parasites are introduced once the flock of dairy goat farms has been constituted, these data indicate for the three studied species that rare pre-existing benzimidazole-resistance alleles already present before the isolation of populations had been selected. On the other hand, the fact that some benzimidazole-resistance alleles were specific to one population of T. circumcincta or H. contortus, seems to be in agreement with the hypothesis of the selection of spontaneous mutations. Thus, the origin of benzimidazole-resistance alleles in trichostrongylid nematodes seems to involve primarily the selection of rare alleles and possibly of spontaneous mutations.     

Genetics of the Tubulin Gene Families of Physarum

The organization of the alpha- and beta-tubulin gene families in Physarum was investigated by Mendelian analysis. Restriction endonuclease-generated DNA fragments homologous to alpha- and beta-tubulin show length polymorphisms that can be used as markers for genetic mapping. Analysis of meiotic assortment among progeny of heterozygotes allowed alpha- and beta-tubulin sequence loci to be defined. There are four unlinked alpha-tubulin sequence loci (altA, altB, altC and altD) and at least three unlinked beta-tubulin sequence loci (betA, betB and betC). The alpha-tubulin loci are not linked to the beta-tubulin loci. --Segregation of tubulin sequence loci with respect to ben mutations that confer resistance to antitubulin benzimidazole drugs was used to investigate whether any members of the alpha- or beta-tubulin gene families are allelic to ben loci. The beta-tubulin sequence locus betB is allelic to the resistance locus benD, the betA locus is probably allelic to benA and the alpha-tubulin sequence locus altC may be allelic to benC. The molecular implications of benzimidazole resistance phenotypes when only one of the expressed beta-tubulin gene family members mutates to drug resistance are discussed in relation to tubulin function.     

Identification of an amino acid substitution in the benA, beta-tubulin gene of Aspergillus nidulans that confers thiabendazole resistance and benomyl supersensitivity

We are using molecular genetic techniques to identify sites of interaction of beta-tubulin with benzimidizole anti-microtubule agents. We have developed a marker-rescue technique for cloning mutant alleles of the benA, beta-tubulin gene of Aspergillus nidulans and have used the technique to clone two mutant benA alleles, benA16 and benA19. These are the only A. nidulans alleles known to confer resistance to the benzimidazole antimicrotubule agent thiabendazole and supersensitivity to other benzimidazole antimicrotubule agents including benomyl and its active breakdown product, carbendazim. benA16 has been shown, moreover, to reduce thiabendazole binding to beta-tubulin. We have sequenced the two mutant alleles and have found that they carry different nucleotide changes that cause the same single amino acid substitution, valine for alanine at amino acid 165. Since thiabendazole and carbendazim differ at only one side chain, the R2 group, we conclude that the region around amino acid 165 is involved in the binding of the R2 group of benzimidazole antimicrotubule agents to beta-tubulin.     

The Emergence of Resistance to the Benzimidazole Anthlemintics in Parasitic Nematodes of Livestock Is Characterised by Multiple Independent Hard and Soft Selective Sweeps

Anthelmintic resistance is a major problem for the control of parasitic nematodes of livestock and of growing concern for human parasite control. However, there is little understanding of how resistance arises and spreads or of the “genetic signature” of selection for this group of important pathogens. We have investigated these questions in the system for which anthelmintic resistance is most advanced; benzimidazole resistance in the sheep parasites Haemonchus contortus and Teladorsagia circumcincta. Population genetic analysis with neutral microsatellite markers reveals that T. circumcincta has higher genetic diversity but lower genetic differentiation between farms than H. contortus in the UK. We propose that this is due to epidemiological differences between the two parasites resulting in greater seasonal bottlenecking of H. contortus. There is a remarkably high level of resistance haplotype diversity in both parasites compared with drug resistance studies in other eukaryotic systems. Our analysis suggests a minimum of four independent origins of resistance mutations on just seven farms for H. contortus, and even more for T. circumincta. Both hard and soft selective sweeps have occurred with striking differences between individual farms. The sweeps are generally softer for T. circumcincta than H. contortus, consistent with its higher level of genetic diversity and consequent greater availability of new mutations. We propose a model in which multiple independent resistance mutations recurrently arise and spread by migration to explain the widespread occurrence of resistance in these parasites. Finally, in spite of the complex haplotypic diversity, we show that selection can be detected at the target locus using simple measures of genetic diversity and departures from neutrality. This work has important implications for the application of genome-wide approaches to identify new anthelmintic resistance loci and the likelihood of anthelmintic resistance emerging as selection pressure is increased in human soil-transmitted nematodes by community wide treatment programs.     

Deep amplicon sequencing as a powerful new tool to screen for sequence polymorphisms associated with anthelmintic resistance in parasitic nematode populations

Parasitic gastrointestinal nematodes contribute to significant human morbidity and cause billions of dollars per year in lost agricultural production. Control is dependent on the use of anthelmintic drugs which, in the case of livestock parasites, is severely compromised by the widespread development of drug resistance. There are now concerns regarding the emergence of anthelmintic resistance in parasitic nematodes of humans in response to the selection pressure resulting from mass drug administration programs. Consequently, there is an urgent need for sensitive, scalable and accurate diagnostic tools to detect the emergence of anthelmintic resistance. Detecting and measuring the frequency of resistance-associated mutations in parasite populations has the potential to provide sensitive and quantitative assessment of resistance emergence from an early stage. Here, we describe the development and validation of deep amplicon sequencing as a powerful new approach to detect and quantify the frequency of single nucleotide polymorphisms associated with benzimidazole resistance. We have used parasite communities in sheep to undertake a proof-of-concept study of this approach. Sheep provide an excellent host system, as there are multiple co-infecting trichostrongylid nematode species, each likely with a varying prevalence of benzimidazole resistance. We demonstrate that the approach provides an accurate measure of resistance allele frequencies, and can reliably detect resistance alleles down to a frequency of 0.1%, making it particularly valuable for screening mutations in the early stages of resistance. We illustrate the power of the technique by screening UK sheep flocks for benzimidazole resistance-associated single nucleotide polymorphisms at three different codons of the β-tubulin gene in seven different parasite species from 164 populations (95 from ewes and 69 from lambs) in a single MiSeq sequencing run. This approach provides a powerful new tool to screen for the emergence of anthelmintic resistance mutations in parasitic nematode populations of both animals and humans.     

野生太行猕猴肠道寄生虫感染的初步调查

为了解野生太行猕猴(Macaca mulatta tcheliensis)肠道寄生虫感染情况,采用饱和蔗糖溶液漂浮法、改良抗酸染色法和卢戈碘液染色法对194份粪便样品进行检查,结果共检查出阿米巴原虫(Amoebasp.)、等孢球虫(Isospora sp.)、鞭虫(Trichuris sp.)、圆线虫(Strongylid sp.)、泡翼线虫(Physaloptera sp.)、蛲虫(Enterobius叩.)、蛔虫(Ascarid sp.)7种寄生虫.寄生虫总感染率为97.4%,线虫感染率为93.8%.其中,鞭虫感染率最高,为91.2%;阿米巴原虫和圆线虫次之,分别为67.5%和60.8%.寄生虫混合感染率为87.1%;鞭虫、泡翼线虫、圆线虫3种线虫混合感染率为71%,同时感染其中2种的占42%,同时感染3种的占29%.对鞭虫卵、圆线虫卵和泡翼线虫卵进行感染强度测定,92%以上阳性样品线虫卵的EPG值在2 000以下.等孢球虫为非人灵长类国内新发现种;所检出寄生虫中,阿米巴原虫、鞭虫、泡翼线虫、蛲虫、蛔虫均可感染人,因此具有重要的公共卫生意义.     

Genetic Variability of the β-Tubulin Genes in Benzimidazole-Susceptible and -Resistant Strains of Haemonchus Contortus

Benzimidazole anthelmintics are the most common chemotherapeutic agents used to remove intestinal helminths from farm animals. The development of drug resistance within helminth populations is wide-spread and can render these drugs essentially useless. The mechanism of benzimidazole resistance appears to be common to many species ranging from fungi to nematodes and involves alterations in the genes encoding β-tubulin. During the selection process resulting in resistance, there must be quantitative changes in the population gene pool. Knowledge of these changes would indicate the mechanisms underlying the spread of resistance in the population, which in turn could be used to design more effective drug administration strategies. To this end we have identified allelic variation at two β-tubulin genes in Haemonchus contortus using restriction map analysis of individual adults. Extremely high levels of variation were identified at both loci within a susceptible strain. In two independently derived benzimidazole resistant strains, allele frequencies at both loci were significantly different from the susceptible strain but not from each other. The same alleles at both loci, in both resistant strains, were favored by selection with benzimidazoles, suggesting that both loci are involved in determining benzimidazole resistance. These data confirm that changes in allele frequency, rather than novel genetic rearrangements induced by exposure to the drug, explain the changes associated with benzimidazole resistance. These results also show that any DNA based test for the development of benzimidazole resistance must take into account the frequency of alleles present in the population and not simply test for the presence or absence of specific allelic types.     

Real-time PCR assays for monitoring benzimidazole resistance-associated mutations in Ancylostoma caninum

Frequent and broad application of anthelmintic drugs for treatment of intestinal parasite infection has led to drug resistance that often renders whole populations of livestock unresponsive to treatment. Therefore, it is important to detect mutations associated with drug resistance before it becomes clinically manifest. To monitor developing drug resistance against benzimidazoles (BZ), we developed real-time PCR assays and applied them to analyse the beta-tubulin isotype-1 gene of the hookworm Ancylostoma caninum, an important parasite of dogs. Previously, we developed PCR assays to monitor codon positions 167 and 200. Here, we describe an assay which is able to detect resistance alleles in codon 198. These real-time PCR assays were subsequently applied to screen hookworm specimens recovered from dogs in Georgia. No elevated levels of polymorphisms at the investigated loci were found, suggesting that selection for resistance in the tested samples did not occur.     

Variation in resistance to benzimidizole in different biocontrol agents based on protein sequence homology

Benzimidazole resistance in thirty isolates of the biocontrol agents was studied with reference to specific mutations in the beta-tubulin genes. Our results suggest that apart from the correlation between specific mutations in the beta-tubulin gene and variation in resistance, the overall ratio of polar and non-polar amino-acids may also play a vital role in conferring benzimidazole resistance to these fungi.     

Beta-tubulin complementary DNA sequence variations observed between cyathostomins from benzimidazole-susceptible and -resistant populations

The molecular mechanism of benzimidazole (BZ) resistance in cyathostomins of horses is still unclear. Previous studies revealed that the TTC or TAC polymorphism in codon 200 of the beta-tubulin isotype 1 gene is not as strictly correlated with BZ resistance as in trichostrongyles in sheep. To identify further sites of polymorphism within the beta-tubulin gene related to BZ resistance, complete complementary DNAs (cDNAs) encoding beta-tubulin of adult worms of Cylicocyclus nassatus, Cyathostomum pateratum, Cyathostomum coronatum, Cyathostomum catinatum, Cylicostephanus longibursatus, and Cylicostephanus goldi of a BZ-resistant cyathostomin population were characterized using specific primers. The cDNA sequence of each species spans 1,429 bp, encoding a protein of 448 amino acids. The interspecific identities are 95.2-99.6% at the nucleotide and 98.7-100.0% at the peptide level. The comparison of the amino acid sequences of individuals isolated from the BZ-resistant cyathostomin population with those from individuals of Cc. nassatus, Cy. coronatum, Cy. pateratum, and Cy. catinatum of a BZ-susceptible one showed differing amino acids in 11 positions. The commonness of a phenylalanine to tyrosine mutation at position 167 in all the 6 cyathostomin species isolated from a BZ-resistant population suggests its involvement in the molecular mechanism in BZ resistance.     

Multi-sited mutations of beta-tubulin are involved in benzimidazole resistance and thermotolerance of fungal biocontrol agent Beauveria bassiana

Fungicide resistance and thermotolerance of biocontrol agents in mitosporic fungi are of merits for enhancing fungal formulations against insect pests in the field. Among 20 wild strains of Beauveria bassiana (a well-known fungal biocontrol agent) tested in this study, 19 were sensitive or highly sensitive to carbendazim (methyl 2-benzimidazole carbamate), a typical benzimidazole fungicide, despite low resistance found in one strain. Sequential mutagenesis of a carbendazim-sensitive wild strain [minimal inhibitory concentration (MIC) = 1.32 microg ml(-1)] under artificial selection pressure generated 11 mutants sharing a common MIC of > 1000 microg ml(-1) without visible variation in colony growth and conidiation capacity. This represents at least 758-fold enhancement of the resistance among the mutants. However, accompanied with the enhanced resistance, all the mutants became less thermotolerable. Stressed at 48 degrees C, conidial LT(50)s of the mutants varied from 1.8 to 9.6 min and were lower than the parental LT(50) (36 min). Moreover, the contents of hydrophobin-like proteins in conidial walls declined significantly among the mutants compared with that of the wild parent. Mutations commonly relating to benzimidazole resistance in fungi were located at Q134, F167 and/or E198 around the taxol-binding site of beta-tubulin by sequencing the beta-tubulin of the mutants. Also, mutations of other 37 amino acid residues in the sequences (each having one to five residues mutated) were found for the first time and they were diverse in spatial structure. All mutations restricted to the half of beta-tubulin close to alpha-tubulin were likely involved in variation in each of the traits concerned but their interactions were complicated.     

Systematic mutational analysis of the yeast beta-tubulin gene.

A systematic strategy was used to create a synoptic set of mutations that are distributed throughout the single beta-tubulin gene of Saccharomyces cerevisiae. Clusters of charged amino acids were targeted for mutagenesis and converted to alanine to maximize alterations on the protein's surface and minimize alterations that affect protein folding. Of the 55 mutations we constructed, three confer dominant-lethality, 11 confer recessive-lethality, 10 confer cold-sensitivity, one confers heat-sensitivity, and 27 confer altered resistance to benomyl. Only 11 alleles give no discernible phenotype. In spite of the fact that beta-tubulin is a highly conserved protein, three-fourths of the mutations do not destroy the ability of the protein to support the growth of yeast at 30 degrees C. The lethal substitutions are primarily located in three regions of the protein and presumably identify domains most critical for beta-tubulin function. Interestingly, most of the conditional-lethal alleles produce specific defects in spindle assembly at their restrictive temperature; cytoplasmic microtubules are relatively unaffected. The exceptions are two mutants that contain abnormally long cytoplasmic microtubules. Mutants with specific spindle defects were not observed in our previous collection of beta-tubulin mutants and should be valuable in dissecting spindle function.     

Analysis of the beta-tubulin genes from Enterocytozoon bieneusi isolates from a human and rhesus macaque

Enterocytozoon bieneusi is the most common and clinically significant microsporidium associated with chronic diarrhea and wasting in immunocompromised humans. Albendazole, which is effective against several helminths, protozoa, and microsporidia, is relatively ineffective against infections due to E. bieneusi. A likely explanation for the observed clinical resistance to albendazole was discovered from sequence analysis of the E. bieneusibeta-tubulin from isolates from an infected human and a naturally infected rhesus macaque. The beta-tubulin of E. bieneusi has a substitution at Glu(198), which is one of six amino acids reported to be associated with benzimidazole sensitivity.     

Identification and characterization of benzimidazole resistance in Monilinia fructicola from stone fruit orchards in California

Low and high levels of resistance to the benzimidazole fungicides benomyl and thiophanate-methyl were observed in field isolates of Monilinia fructicola, which is the causative agent of brown rot of stone fruit. Isolates that had low levels of resistance (hereafter referred to as LR isolates) and high levels of resistance (hereafter referred to as HR isolates) were also cold and heat sensitive, respectively. Results from microsatellite DNA fingerprints showed that genetic identities among the populations of sensitive (S), LR, and HR isolates were very high (>0.96). Analysis of DNA sequences of the beta-tubulin gene showed that the LR isolates had a point mutation at codon 6, causing a replacement of the amino acid histidine by tyrosine. Codon 198, which encodes a glutamic acid in S and LR isolates, was converted to a codon for alanine in HR isolates. Based on these point mutations in the beta-tubulin gene, allele-specific PCR assays were developed for rapid detection of benzimidazole-resistant isolates of M. fructicola from stone fruit.     

Isolation of mip (microtubule-interacting protein) mutations of Aspergillus nidulans.

We identified four mutations in two previously undescribed loci involved in microtubule function in Aspergillus nidulans as extragenic suppressors of benA33, a heat-sensitive beta-tubulin mutation. Three of the four mutations map to a locus closely linked to riboB on linkage group VIII; we designated this locus mipA (for microtubule-interacting protein). We were not able to map the remaining suppressor because of chromosomal rearrangements. However, since it recombines with riboB at a significantly higher frequency than the mipA alleles, it is unlikely to be in mipA; thus, we designated it mipB1. The mip mutations are not allelic to the previously identified loci that encode alpha- and beta-tubulin, and it is likely that mipA and mipB encode previously unidentified nontubulin proteins involved in microtubule function. Each of the mip mutations suppresses the heat sensitivity conferred by benA33 and suppresses the blockage of nuclear division and movement conferred by this mutation at high temperatures. Interactions between mipA and benA are allele specific. All of the mipA mutations are cryptic in a wild-type benA background but cause cold sensitivity in combination with benA33. These mutations also confer cold sensitivity in combination with benA31 and benA32 and reduce the resistance conferred by these mutations to the antimicrotubule agent benomyl but do not suppress the heat sensitivity conferred by these alleles. Finally, the mipA alleles suppress the heat sensitivity conferred by benA11, benA17, and benA21 but do not confer cold sensitivity in combination with these alleles.     

Allele-specific PCR for the beta-tubulin codon 200 TTC/TAC polymorphism using single adult and larval small strongyle (Cyathostominae) stages

It has been shown that benzimidazole (BZ) resistance in sheep gastrointestinal nematodes is linked with an increase in beta-tubulin codon 200 tyrosine-expressing alleles in the resistant parasite populations. Here, an allele-specific PCR has been developed for the discrimination of the TAC/TTC polymorphism in the beta-tubulin 200 codon of small strongyles. One reverse primer was used in 2 separate amplifications with 1 of 2 forward primers that differed only in their final 3' nucleotide. The primers flank a facultative intron/exon. Therefore, the amplified fragments are either 251 or 308 bp in size, depending on the presence or absence of the intron in individual worms. Amplification of genomic DNA isolated from single adult small strongyles from a set of 7 species consistently generated allele-specific products. Three worms each of the following species were used: Cylicocyclus nassatus, Cylicocyclus insigne, Cylicocyclus elongatus, Cylicocyclus radiatus, Cyathostomum pateratum, Cyathostomum catinatum, and Cyathostomum coronatum. PCR with DNA isolated from single larvae also reproducibly generated specific fragments. This method might be applied for the future assessment of allele frequencies in susceptible and resistant populations to further investigate the mechanism of BZ-resistance in small strongyles.     

Characterisation of a beta-tubulin gene from the liver fluke, Fasciola hepatica

This study represents the first beta-tubulin sequence from a trematode parasite, namely, the liver fluke, Fasciola hepatica. PCR of genomic DNA showed that at least one beta-tubulin gene from F. hepatica contains no introns. A number of amino acids in the primary sequence of fluke tubulin are different from those described previously in various nematode species and the cestode, Echinococcus multilocularis. beta-Tubulin is an important target for benzimidazole anthelmintics, although (with the exception of triclabendazole) they show limited activity against F. hepatica. The amino acid differences in fluke beta-tubulin are discussed in relation to the selective toxicity of benzimidazoles against helminths and the mechanism of drug resistance.     

A molecular tool for species identification and benzimidazole resistance diagnosis in larval communities of small ruminant parasites

This report describes a molecular method for determining in a first step the generic composition of a nematode community and in a second step, the resistance of each species to benzimidazole (BZ). We first established a polymerase chain reaction (PCR) linked to a restriction fragment length polymorphism strategy using the isotype 1 beta-tubulin gene. This method overcame the limitations of morphological identification of larval stages of trichostrongylid nematode species. Geographically distant isolates from the three main gastrointestinal species in temperate zones, Teladorsagia circumcincta, Haemonchus contortus, and Trichostrongylus colubriformis, were distinguished using this method. We then used an allele-specific PCR (AS-PCR) to detect mutations of residue 200 of the beta-tubulin, which is implicated in BZ resistance. The sequences of several samples confirmed the BZ-resistance genotype determined by AS-PCR. The ability to process large numbers of samples simultaneously makes this PCR-based strategy particularly suitable for epidemiological studies. It may also be useful for monitoring the emergence of resistant alleles in nematode communities.