An inhibitor of collagen-stimulated platelet activation from the salivary glands of the Haementeria officinalis leech. II. Cloning of the cDNA and expression.

Salivary glands of the leech Haementeria officinalis contain a protein, leech antiplatelet protein (LAPP), that specifically blocks collagen-mediated platelet aggregation (Connolly, T. M., Jacobs, J. W., and Condra, C. (1992) J. Biol. Chem. 267, 6893-6898). Degenerate oligonucleotides whose sequences were derived from two short peptides from V8 digests of the native LAPP were used as primers to generate a polymerase chain reaction (PCR) product which contains the cDNA region coding for the sequence between these two peptides. Using this PCR product as a hybridization probe, phage containing cDNA clones were isolated containing the entire deduced amino acid sequence for LAPP. Computer analysis of the amino acid sequence predicts a peptidase cleavage site between a 21-residue pre-peptide and a mature protein of 126 amino acids. A DNA insert to express the predicted mature LAPP protein was generated by PCR amplification using phage-derived cDNA clones as a substrate. This insert encoded a fusion protein with the leader sequence of the yeast alpha mating factor and the mature LAPP cDNA. These PCR products were cloned into the yeast expression vector pKH4 alpha 2. A KEX 2 Lys-Arg endopeptidase cleavage site was placed NH2-terminal to the predicted mature protein. This vector transfected into the yeast Saccharomyces cerevisiae directs expression of a secreted mature protein at levels up to 200 mg of LAPP/liter of culture medium. The recombinant protein was comparable to native LAPP in its electrophoretic mobility, its reactivity with anti-LAPP antisera, and its biological activity including inhibition of collagen-stimulated platelet aggregation and the adhesion of platelets to collagen. Availability of significant quantities of recombinant LAPP opens the way to further biochemical structure/function studies and to studies on the effects of an inhibitor of collagen-stimulated platelet aggregation in vivo.     

Molecular cloning and characterization of LKv1, a novel voltage‐gated potassium channel in leech

We have cloned a novel voltage‐gated K channel, LKv1, in two species of leech. The properties of LKv1 expressed in transiently transfected HEK293 cells is that of a delayed rectifier current. LKv1 may be a major modulator of excitability in leech neurons, since antibody localization studies show that LKv1 is expressed in the soma and axons of all neurons in both the central and peripheral nervous systems. Comparison of the biophysical and pharmacological properties of LKv1 with native voltage‐gated conductances in leech neurons suggests that LKv1 may correspond to the previously characterized delayed rectifier current, I_K. Phylogenetic analysis of LKv1 shows that it is related to the Shaker subfamily of voltage‐gated K channels although it occupies a separate branch from that of the monophyletic Shaker clade composed of the flatworm, Aplysia, Drosophila, and mammalian Shaker homologs as well as from that of two recently identified Shaker‐related K channels in jellyfish. Thus, this analysis indicates that this group of voltage‐gated K channels contains several evolutionarily divergent lineages. © 1999 John Wiley & Sons, Inc. J Neurobiol 38: 287–299, 1999     

Conserved genetic regions across angiosperms as tools to develop single‐copy nuclear markers in gymnosperms: an example using cycads

Several individuals of the Caribbean Zamia clade and other cycad genera were used to identify single‐copy nuclear genes for phylogeographic and phylogenetic studies in Cycadales. Two strategies were employed to select target loci: (i) a tblastX search of Arabidopsis conserved ortholog sequence (COS) set and (ii) a tblastX search of Arabidopsis‐Populus‐Vitis‐Oryza Shared Single‐Copy genes (APVO SSC) against the EST Zamia databases in GenBank. From the first strategy, 30 loci were selected, and from the second, 16 loci. In both cases, the matching GenBank accessions of Zamia were used as a query for retrieving highly similar sequences from Cycas, Picea, Pinus species or Ginkgo biloba. After retrieving and aligning all the sequences in each locus, intron predictions were completed to assist in primer design. PCR was carried out in three rounds to detect paralogous loci. A total of 29 loci were successfully amplified as a single band of which 20 were likely single‐copy loci. These loci showed different diversity and divergence levels. A preliminary screening allowed us to select 8 promising loci (40S, ATG2, BG, GroES, GTP, LiSH, PEX4 and TR) for the Zamia pumila complex and 4 loci (COS26, GroES, GTP and HTS) for all other cycad genera.     

The molecular phylogenetics of Trachymyrmex Forel ants and their fungal cultivars provide insights into the origin and coevolutionary history of ‘higher‐attine’ ant agriculture

The fungus‐growing ants and their fungal cultivars constitute a classic example of a mutualism that has led to complex coevolutionary dynamics spanning c. 55–65 Ma. Of the five agricultural systems practised by fungus‐growing ants, higher‐attine agriculture, of which leaf‐cutter agriculture is a derived subset, remains poorly understood despite its relevance to ecosystem function and human agriculture across the Neotropics and parts of North America. Among the ants practising higher‐attine agriculture, the genus Trachymyrmex Forel, as currently defined, shares most‐recent common ancestors with both the leaf‐cutter ants and the higher‐attine genera Sericomyrmex Mayr and Xerolitor Sosa‐Calvo et al. Although previous molecular‐phylogenetic studies have suggested that Trachymyrmex is a paraphyletic grade, until now insufficient taxon sampling has prevented a full investigation of the evolutionary history of this group and limited the possibility of resolving its taxonomy. Here we describe the results of phylogenetic analyses of 38 Trachymyrmex species, including 27 of the 49 described species and at least 11 new species, using four nuclear markers, as well as phylogenetic analyses of the fungi cultivated by 23 species of Trachymyrmex using two markers. We generated new genetic data for 112 ants (402 new gene sequences) and 95 fungi (153 new gene sequences). Our results corroborate previous findings that Trachymyrmex, as currently defined, is paraphyletic. We propose recognizing two new genera, Mycetomoellerius gen.n. and Paratrachymyrmex gen.n. , and restricting the continued use of Trachymyrmex to the clade of nine largely North American species that contains the type species [Trachymyrmex septentrionalis (McCook)] and that is the sister group of the leaf‐cutting ants. Our fungal cultivar phylogeny generally corroborates previously observed broad patterns of ant–fungus association, but it also reveals further violations of those patterns. Higher‐attine fungi are divided into two groups: (i) the single species Leucoagaricus gongylophorus (Möller); and (ii) its sister clade, consisting of multiple species, recently referred to as Leucoagaricus Singer ‘clade B’. Our phylogeny indicates that, although most non‐leaf‐cutting higher‐attine ants typically cultivate species in clade B, some species cultivate L. gongylophorus, whereas still others cultivate fungi typically associated with lower‐attine agriculture. This indicates that the attine agricultural systems, which are currently defined by associations between ants and fungi, are not entirely congruent with ant and fungal phylogenies. They may, however, be correlated with as yet poorly understood biological traits of the ants and/or of their microbiomes.     

Phylogeny and species delimitation of the C‐genome diploid species in Oryza

Abstract The diploid Oryza species with C‐genome type possesses abundant genes useful for rice improvement and provides parental donors of many tetraploid species with the C‐genome (BBCC, CCDD). Despite extensive studies, the phylogenetic relationship among the C‐genome species and the taxonomic status of some taxa remain controversial. In this study, we reconstructed the phylogeny of three diploid species with C‐genome (Oryza officinalis, O. rhizomatis, and O. eichingeri) based on sequences of 68 nuclear single‐copy genes. We obtained a fully resolved phylogenetic tree, clearly indicating the sister relationship of O. officinalis and O. rhizomatis, with O. eichingeri being the more divergent lineage. Incongruent phylogenies of the C‐genome species found in previous studies might result from lineage sorting, introgression/hybridization and limited number of genetic markers used. We further applied a recently developed Bayesian species delimitation method to investigate the species status of the Sri Lankan and African O. eichingeri. Analyses of two datasets (68 genes with a single sample, and 10 genes with multiple samples) support the distinct species status of the Sri Lankan and African O. eichingeri. In addition, we evaluated the impact of the number of sampled individuals and loci on species delimitation. Our simulation suggests that sampling multiple individuals is critically important for species delimitation, particularly for closely related species.     

A complete molecular phylogeny of Claravis confirms its paraphyly within small New World ground‐doves (Aves: Peristerinae) and implies multiple plumage state transitions

The three species in the genus Claravis (Aves: Peristerinae) are unique among members of the small New World ground‐dove clade. All three species inhabit forested areas rather than open scrubby habitat, and exhibit obvious sexual dichromatism. However, the phylogenetic relationships within Claravis remain unknown. The only molecular phylogenetic study to include more than one species of Claravis indicated the genus is paraphyletic. Here we include molecular data from all three Claravis species, including sequences from a museum skin of the previously unsampled Claravis geoffroyi (purple‐winged ground‐dove). Using both mitochondrial and nuclear loci, we estimate phylogenies and divergence times for the small New World ground‐dove clade. We also use ancestral state reconstruction methods to infer the evolution of male blue plumage (and thus sexual dimorphism) in the clade. As in the previous study we recover Claravis as a paraphyletic group, but with Claravis geoffroyi as the sister species to Claravis mondetoura (maroon‐chested ground‐dove). This result has important implications for the evolutionary history of the small New World ground‐dove clade. In particular, we recover multiple independent transitions between the monomorphic and dimorphic plumage states, which perhaps indicates sexual dimorphism arose twice in the group.     

Microglia of medicinal leech (Hirudo medicinalis) express a specific activation marker homologous to vertebrate ionized calcium‐binding adapter molecule 1 (Iba1/alias aif‐1)

The Ionized calcium‐Binding Adapter molecule 1 (Iba1), also known as Allograft Inflammatory Factor 1 (AIF‐1), is a 17 kDa cytokine‐inducible protein, produced by activated macrophages during chronic transplant rejection and inflammatory reactions in Vertebrates. In mammalian central nervous system (CNS), Iba1 is a sensitive marker associated with activated macrophages/microglia and is upregulated following neuronal death or brain lesions. The medicinal leech Hirudo medicinalis is able to regenerate its CNS after injury, leading to a complete functional repair. Similar to Vertebrates, leech neuroinflammatory processes are linked to microglia activation and recruitment at the lesion site. We identified a gene, named Hmiba1, coding a 17.8 kDa protein showing high similarity with Vertebrate AIF‐1. The present work constitutes the first report on an Iba1 protein in the nervous system of an invertebrate. Immunochemistry and gene expression analyses showed that HmIba1, like its mammalian counterpart, is modulated in leech CNS by mechanical injury or chemical stimuli (ATP). We presently demonstrate that most of leech microglial cells migrating and accumulating at the lesion site specifically expressed the activation marker HmIba1. While the functional role of Iba1, whatever species, is still unclear in reactive microglia, this molecule appeared as a good selective marker of activated cells in leech and presents an interesting tool to investigate the functions of these cells during nerve repair events. © 2014 Wiley Periodicals, Inc. Develop Neurobiol 74: 987–1001, 2014     

Identification of ‘Candidatus Phytoplasma asteris' (16SrI‐B) Causing Witches' Broom Disease of Mallotus japonicus in Korea

Mallotus japonicus with witches' broom disease were observed in Jeollabuk‐do, Korea. A phytoplasma from the infected leaves was identified, based on the 16S rDNA, 16S‐23S intergenic spacer region, and fragment of rp operon and tuf gene sequences. The 16S rDNA sequences exhibited maximum (99.7%) similarity with Iranian lettuce phytoplasma, the rp operon sequences exhibited 100% similarity with Goldenrain stunt phytoplasma, and the tuf gene sequences exhibited 99.8% similarity with Japanese spurge yellows phytoplasma. Results of the sequence analysis and phylogenetic studies confirmed that the phytoplasma associated with M. japonicus in Korea was an isolate of Aster Yellows group (subgroup16SrI‐B).     

The complete amino acid sequences of four globins from the land leech Haemadipsa zeylanica var. japonica

The amino acid sequences of four globins from the land leech, Haemadipsa zeylanica var. japonica, were determined using nucleotide sequencing and protein sequencing. The mature globin-molecules were composed of 146 amino acid residues for M-1 globin, 156 for M-2 globin, 143 for D-1 globin, and 149 for D-2 globin. Alignment of the four kinds of globins by Clustal X revealed 22 invariant amino acids. The four globins were 26–33% identical. A striking feature of amino acid alteration was: the replacement of the E7 distal-His of D-1 globin by phenylalanine because histidine is conserved among the rest of the globins of H. zeylanica, those of other representative species (Lumbricus and Tylorrhynchus) of Annelida and most other hemoglobins. A phylogenetic tree constructed of 18 globin structures including two species of leeches, H. zeylanica (a land leech) and Macrobdella decora (a freshwater leech), T. heterochaetus (a representative species of polychaetes), L. terrestris (a representative species of oligochaetes), and human α and β globins strongly indicated that the leech globins first separated from globin lineage of annelids.     

Solenicola setigera is the first characterized member of the abundant and cosmopolitan uncultured marine stramenopile group MAST‐3

Culture‐independent molecular methods based on the amplification, cloning and sequencing of small‐subunit (SSU) rRNA genes are a powerful tool to study the diversity of prokaryotic and eukaryotic microorganisms for which morphological features are not conspicuous. In recent years, molecular data from environmental surveys have revealed several clades of protists lacking cultured and/or described members. Among them are various clades of marine stramenopiles (heterokonts), which are thought to play an essential ecological role as grazers, being abundant and distributed in oceans worldwide. In this work, we show that Solenicola setigera, a distinctive widespread colonial marine protist, is a member of the environmental clade MArine STramenopile 3 (MAST‐3). Solenicola is generally considered as a parasite or an epiphyte of the diatom Leptocylindrus mediterraneus. So far, the ultrastructural, morphological and ecological data available were insufficient to elucidate its phylogenetic position, even at the division or class level. We determined SSU rRNA gene sequences of S. setigera specimens sampled from different locations and seasons in the type locality, the Gulf of Lions, France. They were closely related, though not identical, which, together with morphological differences under electron microscopy, suggest the occurrence of several species. Solenicola sequences were well nested within the MAST‐3 clade in phylogenetic trees. Since Solenicola is the first identified member of this abundant marine clade, we propose the name Solenicolida for the MAST‐3 phylogenetic group.     

Characterization of 13 polymorphic microsatellite loci in the Japanese land leech

We developed 13 polymorphic microsatellite loci of the Japanese land leech (Haemadipsa japonica; Haemadipsidea) using an Illumina MiSeq sequencing approach. A total of 42,064 nuclear DNA contigs were filtered for microsatellite motifs, among which 30,873 simple sequence repeat loci were identified. From these sequences, we selected 30 primer sets, and 13 of these loci were successfully amplified. Polymorphism of the 13 loci was tested using 16 individuals sampled from sixteen populations across Japan. The number of alleles and polymorphism information content varied from 5 to 17 and 0.335 to 0.883, respectively, and observed and expected heterozygosity values ranged from 0.143 to 0.875 and 0.349 to 0.893, respectively, indicating that these loci are polymorphic. Furthermore, we established useful multiplex PCR using these loci. The 13 microsatellite loci described in this paper are the first nuclear microsatellite markers for a land leech species.     

Stratified bacterial community in the bladder of the medicinal leech, Hirudo verbana

Most animals harbour symbiotic microorganisms inside their body, where intimate interactions occur between the partners. The medicinal leech, Hirudo verbana, possesses 17 pairs of excretory bladders that harbour a large number of intracellular and extracellular symbiotic bacteria. In this study, we characterized the bladder symbionts using molecular phylogenetic analyses, transmission electron microscopy (TEM) and fluorescence in situ hybridization (FISH). Restriction fragment length polymorphism (RFLP) and sequence analyses of 16S rRNA gene clone libraries suggested that six bacterial species co‐colonize the leech bladders. Phylogenetic analyses revealed that these species belong to the α‐Proteobacteria (Ochrobactrum symbiont), β‐Proteobacteria (Beta‐1 and Beta‐2 symbionts), δ‐Proteobacteria (Bdellovibrio symbiont) and Bacteroidetes (Niabella and Sphingobacterium symbionts). Species‐specific PCR detection and FISH confirmed the localization of the symbiotic bacteria in the bladders. The Ochrobactrum, Beta‐1, Bdellovibrio and Sphingobacterium symbionts were consistently detected in 13 leeches from two populations, while infection rate of the other symbionts ranged between 20% and 100% in the two leech populations. Transmission electron microscopy observations of the bladders revealed epithelial cells harbouring a number of intracellular bacilli and an additional type of extracellular, rod‐shaped bacteria in the luminal region. Fluorescence in situ hybridization with group‐specific oligonucleotide probes revealed the spatial organization of the bacterial species in the bladder: the Ochrobactrum symbiont was located intracellularly inside epithelial cells; the Bacteroidetes were localized close to the epithelium in the lumen of the bladder; and the Bacteroidetes layer was covered with dense β‐proteobacterial cells. These results clearly demonstrate that a simple but organized microbial community exists in the bladder of the medicinal leech.     

The impact of anchored phylogenomics and taxon sampling on phylogenetic inference in narrow‐mouthed frogs (Anura,Microhylidae)

Despite considerable progress in unravelling the phylogenetic relationships of microhylid frogs, relationships among subfamilies remain largely unstable and many genera are not demonstrably monophyletic. Here, we used five alternative combinations of DNA sequence data (ranging from seven loci for 48 taxa to up to 73 loci for as many as 142 taxa) generated using the anchored phylogenomics sequencing method (66 loci, derived from conserved genome regions, for 48 taxa) and Sanger sequencing (seven loci for up to 142 taxa) to tackle this problem. We assess the effects of character sampling, taxon sampling, analytical methods and assumptions in phylogenetic inference of microhylid frogs. The phylogeny of microhylids shows high susceptibility to different analytical methods and datasets used for the analyses. Clades inferred from maximum‐likelihood are generally more stable across datasets than those inferred from parsimony. Parsimony trees inferred within a tree‐alignment framework are generally better resolved and better supported than those inferred within a similarity‐alignment framework, even under the same cost matrix (equally weighted) and same treatment of gaps (as a fifth nucleotide state). We discuss potential causes for these differences in resolution and clade stability among discovery operations. We also highlight the problem that commonly used algorithms for model‐based analyses do not explicitly model insertion and deletion events (i.e. gaps are treated as missing data). Our results corroborate the monophyly of Microhylidae and most currently recognized subfamilies but fail to provide support for relationships among subfamilies. Several taxonomic updates are provided, including naming of two new subfamilies, both monotypic.     

Using metabarcoding to compare the suitability of two blood‐feeding leech species for sampling mammalian diversity in North Borneo

The application of high‐throughput sequencing (HTS) for metabarcoding of mixed samples offers new opportunities in conservation biology. Recently, the successful detection of prey DNA from the guts of leeches has raised the possibility that these, and other blood‐feeding invertebrates, might serve as useful samplers of mammals. Yet little is known about whether sympatric leech species differ in their feeding preferences, and whether this has a bearing on their relative suitability for monitoring local mammalian diversity. To address these questions, we collected spatially matched samples of two congeneric leech species Haemadipsa picta and Haemadipsa sumatrana from lowland rainforest in Borneo. For each species, we pooled ~500 leeches into batches of 10 individuals, performed PCR to target a section of the mammalian 16S rRNA locus and undertook sequencing of amplicon libraries using an Illumina MiSeq. In total, we identified sequences from 14 mammalian genera, spanning nine families and five orders. We found greater numbers of detections, and higher diversity of OTUs, in H. picta compared with H. sumatrana, with rodents only present in the former leech species. However, comparison of samples from across the landscape revealed no significant difference in mammal community composition between the leech species. We therefore suggest that H. picta is the more suitable iDNA sampler in this degraded Bornean forest. We conclude that the choice of invertebrate sampler can influence the detectability of different mammal groups and that this should be accounted for when designing iDNA studies.     

A β1‐6/β1‐3 galactosidase from Bifidobacterium animalis subsp. lactis Bl‐04 gives insight into sub‐specificities of β‐galactoside catabolism within Bifidobacterium

The Bifidobacterium genus harbours several health promoting members of the gut microbiota. Bifidobacteria display metabolic specialization by preferentially utilizing dietary or host‐derived β‐galactosides. This study investigates the biochemistry and structure of a glycoside hydrolase family 42 (GH42) β‐galactosidase from the probiotic Bifidobacterium animalis subsp. lactis Bl‐04 (BlGal42A). BlGal42A displays a preference for undecorated β1‐6 and β1‐3 linked galactosides and populates a phylogenetic cluster with close bifidobacterial homologues implicated in the utilization of N‐acetyl substituted β1‐3 galactosides from human milk and mucin. A long loop containing an invariant tryptophan in GH42, proposed to bind substrate at subsite + 1, is identified here as specificity signature within this clade of bifidobacterial enzymes. Galactose binding at the subsite ? 1 of the active site induced conformational changes resulting in an extra polar interaction and the ordering of a flexible loop that narrows the active site. The amino acid sequence of this loop provides an additional specificity signature within this GH42 clade. The phylogenetic relatedness of enzymes targeting β1‐6 and β1‐3 galactosides likely reflects structural differences between these substrates and β1‐4 galactosides, containing an axial galactosidic bond. These data advance our molecular understanding of the evolution of sub‐specificities that support metabolic specialization in the gut niche.     

Different filtering strategies of genotyping‐by‐sequencing data provide complementary resolutions of species boundaries and relationships in a clade of sexually deceptive orchids

Ongoing hybridization and retained ancestral polymorphism in rapidly radiating lineages could mask recent cladogenetic events. This presents a challenge for the application of molecular phylogenetic methods to resolve differences between closely related taxa. We reanalyzed published genotyping‐by‐sequencing (GBS) data to infer the phylogeny of four species within the Ophrys sphegodes complex, a recently radiated clade of orchids. We used different data filtering approaches to detect different signals contained in the dataset generated by GBS and estimated their effects on maximum likelihood trees, global F_ST and bootstrap support values. We obtained a maximum likelihood tree with high bootstrap support, separating the species by using a large dataset based on loci shared by at least 30% of accessions. Bootstrap and F_ST values progressively decreased when filtering for loci shared by a higher number of accessions. However, when filtering more stringently to retain homozygous and organellar loci, we identified two main clades. These clades group individuals independently from their a priori species assignment, but were associated with two organellar haplotype clusters. We infer that a less stringent filtering preferentially selects for rapidly evolving lineage‐specific loci, which might better delimit lineages. In contrast, when using homozygous/organellar DNA loci the signature of a putative hybridization event in the lineage prevails over the most recent phylogenetic signal. These results show that using differing filtering strategies on GBS data could dissect the organellar and nuclear DNA phylogenetic signal and yield novel insights into relationships between closely related species.     

COMPARATIVE MOLECULAR EVOLUTION OF NEWLY DISCOVERED PICOCYANOBACTERIAL STRAINS REVEALS A PHYLOGENETICALLY INFORMATIVE VARIABLE REGION OF β‐PHYCOERYTHRIN1

The genetic diversity and phylogenetic position of 10 strains of picocyanobacteria from the Arabian Sea were examined using partial sequences from three loci: 16S rDNA, RNA polymerase rpoC1, and two elements of the phycoerythrin (PE) locus, cpeA and cpeB which encode for the α and β subunit of PE. Nine of the strains showed nearly identical spectral phenotypes based on the in vivo excitation spectrum for PE fluorescence emission and appear to be strains synthesizing a phycourobilin (PUB)–lacking PE. These strains include one, Synechococcus sp. G2.1, already known to be closely related to filamentous cyanobacteria and not to the commonly studied 5.1 subcluster of marine Synechococcus. The 10th strain was a PE‐lacking strain that was of interest because it was isolated from open‐ocean conditions where picocyanobacteria with this phenotype are relatively uncommon. Phylogenetic analysis of the concatenated 16S rDNA and rpoC1 data sets showed that none of the previously described strains were members of the 5.1 subcluster of marine Synechococcus, nor were they closely related to strain G2.1. Instead, they form a well‐supported and previously undescribed clade of cyanobacteria that is sister to Cyanobium. Thus, these strains represent the first PE‐containing Cyanobium from oceanic waters, and the lineage they define includes a strain with a PE‐lacking phenotype from the same environment. Analysis of the PE sequence data showed the PE apoprotein has evolved independently in the G2.1 lineage and the Cyanobium‐like lineage represented by the study strains. It also revealed a hypervariable region of the β‐subunit not described previously; variation in this region shows a pattern among a wide range of PE‐containing organisms congruent with the phylogenetic relationships inferred from other genes. This suggests that the PUB‐lacking spectral phenotype is more likely to have evolved in distantly related phylogenetic lineages by either divergent or convergent evolution than by lateral gene transfer. Both the conserved PE gene sequences and the inferred amino acid sequences for the hypervariable region show considerable divergence among Prochlorococcus PEs, red algal PEs, PUB‐containing PEs from the marine Synechococcus 5.1 subcluster, PEs from the Cyanobium‐like strains, and PEs from other cyanobacteria (including strain G2.1). Thus, it appears that the hypervariable region of the PE gene can be used as a taxon‐specific marker.     

Diversity and selective pressures of anticoagulants in three medicinal leeches (Hirudinida: Hirudinidae,Macrobdellidae)

Although medicinal leeches have long been used as treatment for various ailments because of their potent anticoagulation factors, neither the full diversity of salivary components that inhibit coagulation, nor the evolutionary selection acting on them has been thoroughly investigated. Here, we constructed expressed sequence tag libraries from salivary glands of two species of medicinal hirudinoid leeches, Hirudo verbana and Aliolimnatis fenestrata, and identified anticoagulant‐orthologs through BLASTx searches. The data set then was augmented by the addition of a previously constructed EST library from the macrobdelloid leech Macrobdella decora. The identified orthologs then were compared and contrasted with well‐characterized anticoagulants from a variety of leeches with different feeding habits, including non‐sanguivorous species. Moreover, four different statistical methods for predicting signatures of positive and negative evolutionary pressures were used for 10 rounds each to assess the level and type of selection acting on the molecules as a whole and on specific sites. In total, sequences showing putative BLASTx‐orthology with five and three anticoagulant‐families were recovered in the A. fenestrata and H. verbana EST libraries respectively. Selection pressure analyses predicted high levels of purifying selection across the anticoagulant diversity, although a few isolated sites showed signatures of positive selection. This study represents a first attempt at mapping the anticoagulant repertoires in a comparative fashion across several leech families.