Pyrosequencing analysis identifies discrete populations of Haemonchus contortus from small ruminants

The genus Haemonchus consists of blood-sucking parasitic nematodes in the abomasum of ruminants. Members of this genus are responsible for extensive production losses, particularly of small ruminants in the tropics but are also found in temperate regions. In this study, we examined the internal transcribed spacers-1 and -2 of rRNA in Haemonchus spp. The rRNA region spanning the internal transcribed spacers-1, -2 and the 5.8S rRNA gene was amplified by PCR from each of 10 worms from Swedish sheep, a Swedish goat and Kenyan sheep. The fragments were sequenced and examined with respect to genetic differences fixed among the three isolates. These and additional worms were further analysed with Pyrosequencing technology. This technique allowed us to rapidly analyse 110 individuals in three putative polymorphic nucleotide positions that were initially identified with dideoxy sequencing. The geographical isolates could to a large extent be genetically distinguished, but none of the polymorphic positions were consistent among all individuals within each isolate. Furthermore an alignment of our sequences and a consensus sequence published for Haemonchus contortus revealed two differences in positions 123 and 196 in internal transcribed spacers-2. Although these positions were previously reported as heterogenic, no polymorphism was detected among the 30 worms sequenced in the present study. Modelling of the internal transcribed spacers-2 secondary structure based on our data also identified a new putative long-range interaction. The isolates are best described as populations. In conclusion, consistent differences were not identified and the studied isolates are therefore best described as discrete populations. This study also reveals for the first time the potential of Pyrosequencing technology as a tool in the analysis of nematode population genetics.     

Global patterns reveal strong population structure in Haemonchus contortus, a nematode parasite of domesticated ruminants

We have examined the global population genetic structure of Haemonchus contortus. The genetic variability was studied using both amplified fragment length polymorphism (AFLP) and nad4 sequences of the mitochondrial genome. To examine the performance and information content of the two different marker systems, comparative assessment of population genetic diversity was undertaken in 19 isolates of H. contortus, a parasitic nematode of small ruminants. A total of 150 individual adult worms representing 14 countries from all inhabited continents were analysed. Altogether 1,429 informative AFLP markers were generated using four different primer combinations. Also, the genetic variation was high, which agrees with results from previous AFLP studies of nematode parasites of livestock. The genetic structure was high, indicating limited gene flow between the different isolates and populations from each continent mostly formed monophyletic groups in the phylogenetic analysis. However, for isolates representing Australia, Greece and one laboratory strain that originated from South Africa (WRS), there was no clear genetic relationship between the isolates and the distance between their geographical origins. Basically the same pattern was observed for the mitochondrial marker, although the phylogenetic analysis was less resolved than for AFLP. In contrast with previous findings on the population genetic structure of H. contortus, the calculation of population structure gave high values (Nst=0.59). The strong structure was present also for the four Swedish isolates (Nst=0.16) representing a small geographical area.     

Genetic analysis of inbreeding of two strains of the parasitic nematode Haemonchus contortus

Haemonchus contortus is a sheep parasitic nematode that causes severe economic losses. Previous studies have indicated a high degree of genetic heterogeneity, which is hardly affected by selection for drug resistance. As a tool for the analysis of the population dynamics of H. contortus and its response to drug resistance, we designed a strategy to study the inbreeding of a benzimidazole-sensitive and a benzimidazole-resistant strain. After 15 generations, a theoretical inbreeding coefficient of 0.87 was achieved. The different stages of inbreeding were analysed using restriction fragment polymorphism, microsatellite variability and amplified fragment length polymorphism. Model-based clustering of the amplified fragment length polymorphism genotypes showed that the allele frequencies of the benzimidazole-resistant strain were stable during the last eight generations. In the sensitive strain a gradual shift of allele frequencies was observed, which led to a temporary increase of the genetic diversity around the eight generations.     

Microsatellite variability and heterozygote excess in Elymus trachycaulus populations from British Columbia in Canada

Microsatellite markers were used to assess the genetic diversity and population structure in four populations of Elymus trachycaulus from British Columbia and one population of Elymus alaskanus from Northwest Territories. Fourteen microsatellite loci were used in this study. Our results indicated that E. trachycaulus is highly polymorphic, with an average percentage of polymorphic loci of 96.5% over the four populations. Average expected heterozygosity values (H_E or gene diversity) varied from 0.418 to 0.585 with a mean of 0.497. Most of the genetic variation was found within populations (85%) and the differentiation among populations was found to be 15% (F_st = 0.15). Interpopulation genetic distances corresponded well with the geographic distance between the population sites of origin, as well as morphological characteristics. Tests for Hardy–Weinberg equilibrium (HWE) for all loci and all populations revealed that all loci significantly differ from HWE. Subsequent analysis indicated that departure from HWE at some loci was due to an excess of heterozygotes. Possible explanations for heterozygote excess are discussed. The most likely reason for observed heterozygote excess could be due to the polyploidy nature of the species.     

Multilocus population genetic analysis of the Southwest Pacific malaria vector Anopheles punctulatus

The population structure and history of the cryptic malaria vector species, Anopheles punctulatus (Doenitz), was investigated throughout Papua New Guinea and the Solomon Islands with the aim of detailing genetic subdivisions and the potential for movement through this biogeographically complex region. We obtained larval collections from over 80 sites and utilised a diverse array of molecular markers that evolve through different processes. Individuals were initially identified to species and genotyped using the ribosomal DNA second internal transcribed spacer. DNA sequencing of a single copy nuclear ribosomal protein S9 and the mitochondrial cytochrome oxidase I loci were then investigated and 12 nuclear microsatellite markers were developed and analysed. Our data revealed three genetically distinct populations – one in Papua New Guinea, the second on Buka Island (Bougainville Province, Papua New Guinea), and the third on Guadalcanal Island (Solomon Islands). Genetic differentiation within Papua New Guinea was much lower than that found in studies of other closely related species in the region. The data does suggest that A. punctulatus has undergone a population bottleneck followed by a recent population and range expansion in Papua New Guinea. Humans and regional economic growth may be facilitating this population expansion, as A. punctulatus is able to rapidly occupy human modified landscapes and traverse unsealed roads. We therefore anticipate extensive movement of this species through New Guinea – particularly into the highlands, with a potential increase in malaria frequency in a warming climate – as well as relatively unrestricted gene flow of advantageous alleles that may confound vector control efforts.     

Low genetic differentiation among widely separated populations of the pearl oyster Pinctada fucata as revealed by AFLP

Genetic variation within and among five populations of the pearl oyster Pinctada fucata, from China (Daya Bay, Sanya Bay and Beibu Bay), Japan (Mie Prefecture) and Australia (Port Stephens) was studied using AFLP. Three primer pairs generated 184 loci among which 91.8-97.3% is polymorphic. An overall genetic diversity of 0.38 among populations and an average of 0.37 within populations (ranging from 0.35 in Japanese population to 0.39 in Beibu Bay population) were observed. Genetic differentiation among the five populations is low but significant as indicated by pairwise G_ST (0.0079-0.0404). AMOVA further shows that differentiation is significant among the five populations but is not significant at a broader geographical scale, among the three groups of Chinese, Japanese and Australian populations or among the two groups of Australian and north Pacific populations. The low level of genetic differentiation indicated that P. fucata populations in the west Pacific are genetically linked. Among the five populations, the Australian one is more differentiated from the others, based on both pairwise AMOVA and G_ST analyses, and is genetically isolated by distance as indicated by Mantel test. However, genetic differences among the three Chinese populations are not correlated with the geographic distances, suggesting that Hainan Island and Leizhou Peninsula may act as barriers blocking gene flow.     

Expression and purification of an active cysteine protease of Haemonchus contortus using Caenorhabditis elegans

Many proteolytic enzymes of parasitic nematodes have been identified as possible targets of control. Testing these as vaccine or drug targets is often difficult due to the problems of expressing proteases in a correctly folded, active form in standard expression systems. In an effort to overcome these difficulties we have tested Caenorhabditis elegans as an expression system for a Haemonchus contortus cathepsin L cysteine protease, Hc-CPL-1. Recombinant Hc-CPL-1 with a polyhistidine tag added to the C-terminal was expressed in an active and glycosylated form in C. elegans. Optimal expression was obtained expressing Hc-cpl-1 under control of the promoter of the homologous C. elegans cpl-1 gene. The recombinant protein was purified from liquid cultures by nickel chelation chromatography in sufficient amounts for vaccination studies to be carried out. This study provides proof of principle that active, post-translationally modified parasitic nematode proteases can be expressed in C. elegans and this approach can be extended for expression of known protective antigens.     

Microsatellite diversity and population genetic structure of yellowcheek, Elopichthys bambusa (Cyprinidae) in the Yangtze River

Yellowcheek carp (Elopichthys bambusa) is the only species of genus Elopichthys. It is widely distributed in Chinese freshwaters but currently its populations have declined to threatening level. We examined the genetic diversity and population structure of E. bambusa in the Yangtze River basin. A total of nine polymorphic microsatellite markers were employed to study five populations occurring in middle and lower reaches of the river. The results revealed low-to-moderate genetic diversity. The number of alleles per locus varied between 3 and 8 with an average of 4.8. Observed heterozygosity ranged from 0.15 to a maximum of 1.00. Significant deviations (P < 0.01) from Hardy–Weinberg equilibrium were observed for all the tested locus-population combinations with clear heterozygosity deficits. AMOVA indicated that majority of the variance lies within populations (93.81%) than among the populations (7.05%). Pairwise F_ST and unbiased genetic distance pointed out significant differentiation among the samples from populations with different connections to the Yangtze River. In the UPGMA dendrogram, clustering pattern of populations indicated that most of the populations are reproductively isolated due to anthropogenic interventions. Clustering of PYL and DTL populations shows ongoing gene flow through the mainstream. The recent hydrological alterations and overfishing are major factors shaping the current genetic structure. These results can be helpful for effective management and sustainable conservation of E. bambusa populations.     

Microsatellite analysis reveals the genetic structure and gene flow of the aquatic quillwort Isoetes sinensis, a critically endangered species in China

Isoetes sinensis is a critically endangered aquatic quillwort, occurring in two fragmented sites of China (Xiuning county of Anhui; Jiande county of Zhejiang). During a five-year period (2004-2009), the areas and sizes of the two populations diminished dramatically due to intensive human activities. Genetic structure of the species was investigated using simple sequence repeat makers (SSRs). For seven nuclear microsatellites, high levels of genetic diversity were found within populations (H_E = 0.324-0.447). Strong genetic differentiation was detected between populations (G_ST = 0.376), while weak genetic differentiation was found within populations (G_ST = 0.026-0.080). Notably, in contrast to the source-sink model suggested by previous study (Chen et al., 2009), the migration pattern of I. sinensis along the Xin’an River is best explained by the linear symmetrical non-adjacent flow model (LSNF), which indicates that intensive human activities of recent years have greatly affected the gene exchange pattern among I. sinensis subpopulations.     

SNP genotyping reveals genetic diversity between cultivated landraces and contemporary varieties of tomato

Background

The tomato (Solanum lycopersium L.) is the most widely grown vegetable in the world. It was domesticated in Latin America and Italy and Spain are considered secondary centers of diversification. This food crop has experienced severe genetic bottlenecks and modern breeding activities have been characterized by trait introgression from wild species and divergence in different market classes.

Results

With the aim to examine patterns of polymorphism, characterize population structure and identify putative loci under positive selection, we genotyped 214 tomato accessions (which include cultivated landraces, commercial varieties and wild relatives) using a custom-made Illumina SNP-panel. Most of the 175 successfully scored SNP loci were found to be polymorphic. Population structure analysis and estimates of genetic differentiation indicated that landraces constitute distinct sub-populations. Furthermore, contemporary varieties could be separated in groups (processing, fresh and cherry) that are consistent with the recent breeding aimed at market-class specialization. In addition, at the 95% confidence level, we identified 30, 34 and 37 loci under positive selection between landraces and each of the groups of commercial variety (cherry, processing and fresh market, respectively). Their number and genomic locations imply the presence of some extended regions with high genetic variation between landraces and contemporary varieties.

Conclusions

Our work provides knowledge concerning the level and distribution of genetic variation within cultivated tomato landraces and increases our understanding of the genetic subdivision of contemporary varieties. The data indicate that adaptation and selection have led to a genomic signature in cultivated landraces and that the subpopulation structure of contemporary varieties is shaped by directed breeding and largely of recent origin. The genomic characterization presented here is an essential step towards a future exploitation of the available tomato genetic resources in research and breeding programs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-835) contains supplementary material, which is available to authorized users.     

Evidence for genetic differentiation of Octopus vulgaris (Mollusca, Cephalopoda) fishery populations from the southern coast of Brazil as revealed by microsatellites

Octopus vulgaris is a cephalopod species in several oceans and commonly caught by artisanal and industrial fisheries. In Brazil, O. vulgaris populations are mainly distributed along the southern coast and have been subjected to intensive fishing during recent years. Despite the importance of this marine resource, no genetic study has been carried out to examine genetic differences among populations along the coast of Brazil. In this study, 343 individuals collected by commercial vessels were genotyped at six microsatellite loci to investigate the genetic differences in O. vulgaris populations along the southern coast of Brazil. Genetic structure and levels of differentiation among sampling sites were estimated via a genotype assignment test and F-statistics. Our results indicate that the O. vulgaris stock consists of four genetic populations with an overall significant analogous F_ST (Φ_CT = 0.10710, P < 0.05) value. The genetic diversity was high with an observed heterozygosity of Ho = 0.987. The negative values of F_IS found for most of the loci examined suggested a possible bottleneck process. These findings are important for further steps toward more sustainable octopus fisheries, so that this marine resource can be preserved for long-term utilization.     

Microsatellite analysis of the genetic structure of captive forest musk deer populations and its implication for conservation

Forest musk deer (Moschus berezovskii) are rare as a result of poaching for musk and habitat loss. Some captive populations of forest musk deer have been established for decades in China. However, little genetic information is available for conservation management. In this paper, genetic variations, population structures, and the genetic bottleneck hypothesis were examined using 11 microsatellite loci from captive populations in Miyalo, Jinfeng and Maerkang in Sichuan Province, China. Estimates of genetic variability revealed substantial genetic variation in the three populations. A total of 142 different alleles were observed in 121 forest musk deer and the effective number of alleles per locus varied from 6.76 to 12.95. The average values of observed heterozygosity, expected heterozygosity, and Nei's expected heterozygosity were 0.552, 0.899 and 0.894 respectively. The overall significant (P < 0.001) deficit of heterozygotes because of inbreeding within breeds amounted to 34.5%. The mean F_ST (P < 0.001) showed that approximately 90.2% of the genetic variation was within populations and 9.8% was across populations. The UPGMA diagram, based on Nei's unbiased genetic distance, indicated that the three populations were differentiated into two different groups and it agreed with their origin and history. Bottleneck tests indicated that all three populations have undergone a population bottleneck, suggesting a small effective population size. Acknowledging that the genetic structure of populations has crucial conservation implications, the present genetic information should be taken into account in management plans for the conservation of captive forest musk deer.     

Haemonchus contortus: evaluation of two signal sequence trapping systems for detection of secreted molecules

Given that signal sequences between secreted proteins of different species can be interchanged, it is reasonable to expect that both mammalian and yeast signal sequence trapping (SST) systems would secrete Haemonchus contortus proteins with similar efficiency and quality. To determine if H. contortus cDNAs that contain a signal sequence could re-establish secretion of a reporter protein, mammalian and yeast SST vectors were designed, 10 H. contortus genes selected, and their respective cDNAs cloned into these two SST vectors. The selected molecules included genes known to code for excretory/secretory or membrane-bound proteins as potential test 'positives', and genes known to code for non-secreted proteins as test 'negatives'. While differentiation between secretion and non-secretion was evident in both systems, the results indicated greater efficiency was achieved when the mammalian system was used. Therefore, mammalian SST using COS cells would be a more useful tool to screen H. contortus cDNA libraries for potential secreted and type-1 integral membrane proteins than yeast SST.     

Trypanosoma vivax displays a clonal population structure

African animal trypanosomiasis, or Nagana, is a debilitating and economically costly disease with a major impact on animal health in sub-Saharan Africa. Trypanosoma vivax, one of the principal trypanosome species responsible for the disease, infects a wide host range including cattle, goats, horses and donkeys and is transmitted both cyclically by tsetse flies and mechanically by other biting flies, resulting in a distribution covering large swathes of South America and much of sub-Saharan Africa. While there is evidence for mating in some of the related trypanosome species, Trypanosoma brucei, Trypanosoma congolense and Trypanosoma cruzi, very little work has been carried out to examine this question in T. vivax. Understanding whether mating occurs in T. vivax will provide insight into the dynamics of trait inheritance, for example the spread of drug resistance, as well as examining the origins of meiosis in the order Kinetoplastida. With this in mind we have identified orthologues of eight core meiotic genes within the genome, the presence of which imply that the potential for mating exists in this species. In order to address whether mating occurs, we have investigated a sympatric field population of T. vivax collected from livestock in The Gambia, using microsatellite markers developed for this species. Our analysis has identified a clonal population structure showing significant linkage disequilibrium, homozygote deficits and disagreement with Hardy-Weinberg predictions at six microsatellite loci, indicative of a lack of mating in this population of T. vivax.     

High genetic differentiation and low genetic diversity in Incarvillea younghusbandii, an endemic plant of Qinghai-Tibetan Plateau,revealed by AFLP markers

Incarvillea younghusbandii Sprague (Bignoniaceae) is a perennial herbaceous plant endemic to Qinghai-Tibetan Plateau. As a species of medical and horticulture importance, I. younghusbandii is threatened by over exploitation and habitat fragmentation. In this study, we analyze the genetic diversity and population structure of I. younghusbandii using amplified fragment length polymorphism (AFLP) markers. Our data reveal very low levels of genetic diversity in seven natural populations across Tibet. Specifically, at population level, the average Nei's genetic diversity index (H_E) and Shannon's diversity index (I) were 0.063 and 0.096, respectively. In contrast, high genetic differentiation among populations (Gst = 0.6238, Φ_ST = 0.614) is detected. The results of Neighbor-joining cluster, PCO, and STRUCTURE assignment reveal consistent pattern, suggesting seven well-defined genetic groups that are concordant with their geographical origins. The possible mechanisms and implications of these findings for conservation are discussed.     

Population genetic differentiation of Schisandra chinensis and Schisandra sphenanthera as revealed by ISSR analysis

Schisandra chinensis (Turcz.) Baill. and Schisandra sphenanthera Rehd. et Wils. are well-known Chinese medicinal plants. The population genetic variation of the two species was studied using inter simple sequence repeat (ISSR) molecular markers. High levels of genetic diversity are revealed in both S. chinensis (P = 88.36%, h = 0.2894, I = 0.4396) and S. sphenanthera (P = 84.09%, H = 0.2782, I = 0.4280). However, the population genetic differentiation is significantly different between the two species. The S. sphenanthera harbors as high as 27% of the genetic variation among populations but 73% within populations, whereas in S. chinensis 17% of the genetic variation occurs among populations and 83% within populations. Both significant (P < 0.05) heterozygosity excess and shifted mode of allele frequency distribution are detected in four out of six populations of S. chinensis and one out of five populations of S. sphenanthera, suggesting the occurrence of recent population bottlenecks in the two species. The different patterns of genetic variation in S. chinensis and S. sphenanthera are discussed in relation to their differences in pollination mechanism, geographic distribution and historical events, and the level of gene flow and genetic drift.     

Haemonchus contortus: in vitro and in vivo anthelmintic activity of aqueous and hydro-alcoholic extracts of Hedera helix

In vitro anthelmintic activity of crude extracts of the ripe fruits of Hedera helix was investigated on eggs and adult nematode parasites Haemonchus contortus. Aqueous extract of H. helix was also evaluated for in vivo anthelmintic activity at dose of 1.13 and 2.25 g/kg in sheep artificially infected with H. contortus. ED(50) for egg hatch inhibition was 0.12 and 0.17 mg/ml for aqueous and hydro-alcoholic extracts, respectively. There was no statistically significant difference in the activity of the two extract types (p>0.05). Hydro-alcoholic extract showed better in vitro activity against adult parasites compared to the aqueous extract. Significant faecal egg count reduction (FECR) was detected in groups treated with both doses of H. helix (p<0.05) on day 2 post-treatment. On day 7 post-treatment significant reduction was detected only for higher dose of H. helix (p<0.05) while on day 14 post-treatment there was no significant FECR in both groups treated with H. helix. The percentage of larvae recovered from culturing faeces obtained from groups of sheep treated with lower and higher doses of H. helix was 47.52% and 36.07%, respectively, which was significantly lower than (p<0.05) that recovered from the control group (60%). Significant (p<0.05), dose dependent total worm count reduction (WCR) was observed for groups of sheep treated with H. helix. Increasing the dose of H. helix improved the efficacy against the male than the female parasites. Treatment with both doses of H. helix helped the animals maintain their packed cell volume (PCV) unlike the untreated control group. The overall findings of the current study indicated that H. helix has a potential anthelmintic benefit and further in vitro and in vivo evaluation of the different parts and fractions is needed to make use of this plant for therapeutic purposes.     

Non-destructive genotyping and genetic variation of fanning in a honey bee colony

The relationship between workers from different patrilines in a naturally mated queen honey bee colony is very complex due to queen polyandry, and still poorly characterized. Here, we report a means of determining the genotype of living workers in a natural honey bee colony by a new non-destructive method, which makes it possible to observe the relationship between behaviours and genotypes. DNA was extracted from the exuvia, found at the bottom of each brood cell, and confirmed to be identical to the DNA extracted from the thorax muscle of the bee emerging from that particular cell. The genotypes were thus determined using DNA from the exuviae without having to hurt or kill the organisms. The emerging workers were marked with coloured, numbered tags to enable behavioural observations over their entire life. Using this new method, we determined 20 patrilines in a naturally mated queen colony, and discovered that the patriline composition of bees exhibiting fanning behaviour was significantly different from the patriline composition of the whole colony. Our results confirm that the genetic structure of a natural insect society plays a fundamental role in the division of labour. The new non-destructive method reveals a novel avenue for the determination of relationships between the behaviours and genes of social insects.