Elongation factor-1 alpha occurs as two copies in bees: implications for phylogenetic analysis of EF-1 alpha sequences in insects

Co-phylogeny of a hyper-symbiotic system: Endosymbiotic bacteria (Gammaproteobacteria), chewing lice (Insecta: Phthiraptera) and birds (Passeriformes)

Chewing lice are hosts to endosymbiotic bacteria as well as themselves being permanent parasites. This offers a unique opportunity to examine the cophylogenetic relationships between three ecologically interconnected organismal groups: birds, chewing lice, and bacteria. Here, we examine the cophylogenetic relationships between lice in the genus Guimaraesiella Eichler, 1949, their endosymbiotic Sodalis-allied bacteria, and a range of bird species from across South China. Both event and distance-based cophylogenetic analyses were explored to compare phylogenies of the three organismal groups. Pair-wise comparisons between lice-endosymbionts and bird-endosymbionts indicated that their evolutionary histories are not independent. However, comparisons between lice and birds, showed mixed results; the distance-based method of ParaFit indicated that their evolutionary histories are not independent, while the event-based method of Jane indicated that their phylogenies were no more congruent than expected by chance. Notably, louse host-switching does not seem to have affected bacterial strains, as conspecific lice sampled from distantly related hosts share bacteria belonging to the same clade.

Prevalence and diversity of parasitic bird lice (Insecta: Psocodea) in northeast Arkansas

Many groups of parasites lack basic information on biodiversity and host associations, which poses challenges for conservation and understanding the ecological relationships between hosts and their parasites. This gap in knowledge is particularly relevant for parasitic species with obscure lifestyles. Ectoparasitc bird lice (Insecta: Psocodea: Phthiraptera) are a group of parasites that has received a relatively substantial research focus, yet patterns of bird-louse relationships and louse diversity remain understudied in many geographic regions, including in parts of the southeastern United States. In this study, we assessed the diversity, prevalence, abundance, and intensity of lice from live and salvaged birds in northeastern Arkansas. We also focused on the frequency of co-occurrence of lice and symbiotic feather mites. Finally, we used nuclear and mitochondrial genes to assess the phylogenic relationships among the most common genera of lice in our sample. We found a total louse prevalence of 10.57% with the highest prevalence on the Passeriformes families Turdidae, Passerellidae, and Parulidae. We also found the louse genera Myrsidea and Brueelia to be the most prevalent and abundant in our sample. Additionally, we reported several novel associations among well-studied bird species. We also found that louse phylogenic patterns tend to reflect host taxonomy and/or ecology. Overall, our results provide important insight into the biodiversity, community structure, and host interactions of parasitic lice from North American birds.

Molecular phylogenetics of the sucking louse genus Lemurpediculus (Insecta: Phthiraptera), ectoparasites of lemurs, with descriptions of three new species

Sucking lice live in intimate association with their hosts and often display a high degree of host specificity. The present study investigated sucking lice of the genus Lemurpediculus from six mouse lemur (Microcebus) and two dwarf lemur (Cheirogaleus) species endemic to the island of Madagascar, considered a biodiversity hotspot. Louse phylogenetic trees were created based on cytochrome C oxidase subunit I (COI), elongation factor 1α (EF1α) and internal transcribed spacer 1 (ITS1) sequences. While clustering according to host species was generally observed for COI and ITS1, suggesting high host specificity of the examined lice, EF1α sequences alone did not distinguish between lice of different Microcebus species, possibly due to rather recent divergence. As bootstrap support for basal tree structure was rather low, further data are necessary to resolve the evolutionary history of louse-mouse lemur associations. Three new species of sucking lice are described: Lemurpediculus zimmermanni sp. Nov. From Microcebus ravelobensis, Lemurpediculus gerpi sp.nov. from Microcebus gerpi, and Lemurpediculus tsimanampesotsae sp. nov. from Microcebus griseorufus. These new species are compared with all known congeneric species and identifying features are illustrated for all known species of Lemurpediculus.

Detailed morphological structure and phylogenetic relationships of Degeeriella punctifer (Phthiraptera: Philopteridae), a parasite of the bearded vulture Gypaetus barbatus (Accipitriformes: Accipitridae)

Habitat loss is one of the main threats to species survival and, in the case of parasites, it is their hosts that provide their habitat. Therefore, extinction even at local scale of host taxa also implies the extinction of their parasites in a process known as co-extinction. This is the case of the bearded vulture (Gypaetus barbatus), which almost became extinct at the beginning of the twentieth century. After several attempts, this species was successfully reintroduced into the Alps at the end of the twentieth century. We collected 25 lice specimens from an electrocuted bearded vulture from Susa (Italian Alps) that were morphologically identified as Degeeriella punctifer. Six individuals were studied by scanning electron microscopy, with particular emphasis on their cephalic sensorial structures, while four further specimens were characterized at molecular level by amplifying partial regions of the 12SrRNA, COX1 and elongation factor 1 alpha (EF-1) genes. From a morphological perspective, the number, type and arrangement of the sensillae on the two distal antennal segments is quite similar to that of other species of the family Philopteridae (Phthiraptera: Ischnocera). The mandibles and tarsal claws allow lice to cling firmly to their host's feathers. Phylogenetic analyses help unravel the paraphyletic nature of the genus Degeeriella and demonstrate the clear differentiation between lice parasitizing Accipitriformes and Falconiformes, as well as the close relationship between D. punctifer, D. fulva, D. nisus and Capraiella sp. that, along with other genera, parasitize rollers (Aves: Coraciiformes).

Molecular Identification of Newly Recorded Louse Columbicola tschulyschman Eichler (Phthiraptera: Philopteridae) From Saudi Arabia With a Note on Genus Columbicola spp. Phylogeny

Chewing lice comprise a large group of ectoparasites that colonize and adversely affect several domestic and wild birds including pigeons. In Saudi Arabia, there is a lack of studies describing such ectoparasites and their infestation rates. Through this work, a new record, Columbicola, tschulyschman Eichler (C. tschulyschman Eichler) was collected from domestic pigeons (Columba livia domestica, Linnaeus). The collected C. tschulyschman Eichler was morphologically identified based on specific taxonomic keys. Mitochondrial (COI) and nuclear (EF-1α) gene fragments were used for molecular identification and phylogenetic reconstruction. In this study, the C. tschulyschman Eichler accounted for around 69.40%. To our knowledge, this is the first report of C. tschulyschman Eichler in Riyadh, Saudi Arabia. To improve the tree topology and differentiate between genera, further studies should utilize the 16s rRNA.

Molecular phylogenetics of the avian feather louse Philopterus-complex (Phthiraptera: Philopteridae)

The avian feather louse Philopterus-complex (Phthiraptera: Ischnocera: Philopteridae) currently contains 12 genera that have been grouped together because of shared morphological characteristics. Although previously lumped into a single genus (Philopterus), more recent morphological treatments have separated the group into several different genera. Here we evaluate the status of these genera using DNA sequence data from 118 ingroup specimens belonging to ten genera in the Philopterus-complex: Australophilopterus Mey, 2004, Cinclosomicola Mey 2004, Clayiella Eichler, 1940, Corcorides Mey, 2004, Mayriphilopterus Mey, 2004, Paraphilopterus Mey 2004, Philopteroides Mey 2004, Philopterus Nitzsch, 1818, Tyranniphilopterus Mey, 2004, and Vinceopterus Gustafsson, Lei, Chu, Zou, and Bush, 2019. Our sampling includes 97 new louse-host association records. Our analyses suggest that the genus Debeauxoecus Conci, 1941, parasitic on pittas (Aves: Pittidae), is outside of the Philopterus-complex, and that there is strong support for the monophyly of a group containing the remaining genera from the complex. Some diverse genera, such as Philopterus (sensu stricto) and Mayriphilopterus are supported as monophyletic, whereas the genera Australophilopterus, Philopteroides, and Tyranniphilopterus are not. The present study is the largest phylogenetic reconstruction of avian lice belonging to the Philopterus-complex to date and suggests that further generic revision is needed in the group to integrate molecular and morphological information.

Using Proteomic Approaches to Unravel the Response of Ctenocephalides felis felis to Blood Feeding and Infection With Bartonella henselae

Among the Ctenocephalides felis felis-borne pathogens, Bartonella henselae, the main aetiological agent of cat scratch disease (CSD), is of increasing comparative biomedical importance. Despite the importance of B. henselae as an emergent pathogen, prevention of the diseases caused by this agent in cats, dogs and humans mostly relies on the use of ectoparasiticides. A vaccine targeting both flea fitness and pathogen competence is an attractive choice requiring the identification of flea proteins/metabolites with a dual effect. Even though recent developments in vector and pathogen -omics have advanced the understanding of the genetic factors and molecular pathways involved at the tick-pathogen interface, leading to discovery of candidate protective antigens, only a few studies have focused on the interaction between fleas and flea-borne pathogens. Taking into account the period of time needed for B. henselae replication in flea digestive tract, the present study investigated flea-differentially abundant proteins (FDAP) in unfed fleas, fleas fed on uninfected cats, and fleas fed on B. henselae-infected cats at 24 hours and 9 days after the beginning of blood feeding. Proteomics approaches were designed and implemented to interrogate differentially expressed proteins, so as to gain a better understanding of proteomic changes associated with the initial B. henselae transmission period (24 hour timepoint) and a subsequent time point 9 days after blood ingestion and flea infection. As a result, serine proteases, ribosomal proteins, proteasome subunit α-type, juvenile hormone epoxide hydrolase 1, vitellogenin C, allantoinase, phosphoenolpyruvate carboxykinase, succinic semialdehyde dehydrogenase, glycinamide ribotide transformylase, secreted salivary acid phosphatase had high abundance in response of C. felis blood feeding and/or infection by B. henselae. In contrast, high abundance of serpin-1, arginine kinase, ribosomal proteins, peritrophin-like protein, and FS-H/FSI antigen family member 3 was strongly associated with unfed cat fleas. Findings from this study provide insights into proteomic response of cat fleas to B. henselae infected and uninfected blood meal, as well as C. felis response to invading B. henselae over an infection time course, thus helping understand the complex interactions between cat fleas and B. henselae at protein levels.

Transcript level is a key factor affecting RNAi efficiency

Efficiency is the basis for the application of RNA interference (RNAi) technology. Actually, RNAi efficiency varies greatly among insect species, tissues and genes. Previous efforts have revealed the mechanisms for variation among insect species and tissues. Here, we investigated the reason for variable efficiency among the target genes in the same insect. First, we tested the genes sampled randomly from Tribolium castaneum, Locusta migratoria and Drosophila S2 cells for both their expression levels and sensitivity to RNAi. The results indicated that the genes with higher expression levels were more sensitive to RNAi. Statistical analysis showed that the correlation coefficients between transcript levels and knockdown efficiencies were 0.8036 (n = 90), 0.7255 (n = 18) and 0.9505 (n = 13), respectively in T. castaneum, L. migratoria and Drosophila S2 cells. Subsequently, ten genes with varied expression level in different tissues (midgut and carcass without midgut) of T. castaneum were tested. The results indicated that the higher knockdown efficiency was always obtained in the tissue where the target gene expressed higher. In addition, three genes were tested in different developmental stages, larvae and pupae of T. castaneum. The results found that when the expression level increased after insect pupation, these genes became more sensitive to RNAi. Thus, all the proofs support unanimously that transcript level is a key factor affecting RNAi sensitivity. This finding allows for a better understanding of the RNAi efficiency variation and lead to effective or efficient use of RNAi technology.

Chewing Lice of Fregata magnificens with First Record of Fregatiella aurifasciata (Phthiraptera: Amblycera) in Brazil

The genus Fregata includes 5 species, with 3 recorded in Brazil, with Fregata magnificens being the most abundant. However, its ectoparasitic fauna is still little known. This study aimed to evaluate the incidence of ectoparasites of F. magnificens residing along the coast of Rio de Janeiro and São Paulo collected by 2 animal rehabilitation centers. Samples were collected from 5 frigatebirds of the Instituto Argonauta in São Paulo and 10 frigatebirds of the Centro de Recuperação de Animais Selvagens (CRAS) in Rio de Janeiro. Species of lice were identified using both morphological and molecular methods. Scanning electron microscopy was also used for identification. Colpocephalum spineum, Fregatiella aurifasciata, and Pectinopygus fregatiphagus were identified. All 3 louse species have previously been recorded from this host outside Brazil, but only P. fregatiphagus has been recorded from Brazil. This paper reports the first occurrence of F. aurifasciata and C. spineum in Brazil. It is also the first record of P. fregatiphagus in the state of Rio de Janeiro.

Combining Nuclear and Mitochondrial Loci Provides Phylogenetic Information in the Philopterus Complex of Lice (Psocodea: Ischnocera: Philopteridae)

The Philopterus Complex includes several lineages of lice that occur on birds. The complex includes the genera Philopterus (Nitzsch, 1818; Psocodea: Philopteridae), Philopteroides (Mey, 2004; Psocodea: Philopteridae), and many other lineages that have sometimes been regarded as separate genera. Only a few studies have investigated the phylogeny of this complex, all of which are based on morphological data. Here we evaluate the utility of nuclear and mitochondrial loci for recovering the phylogeny within this group. We obtained phylogenetic trees from 39 samples of the Philopterus Complex (Psocodea: Philopteridae), using sequences of two nuclear (hyp and TMEDE6) and one mitochondrial (COI) marker. We evaluated trees derived from these genes individually as well as from concatenated sequences. All trees show 20 clearly demarcated taxa (i.e., putative species) divided into five well-supported clades. Percent sequence divergence between putative species (~5-30%) for the COI gene tended to be much higher than those for the nuclear genes (~1-15%), as expected. In cases where species are described, the lineages identified based on molecular divergence correspond to morphologically defined species. In some cases, species that are host generalists exhibit additional underlying genetic variation and such cases need to be explored by further future taxonomic revisions of the Philopterus Complex.

A New Species of Sucking Louse from the Mandrill from Gabon with a Review of Host Associations and Geographical Distributions, and Identification Keys to Members of the Genus pedicinus (Phthiraptera: Anoplura: Pedicinidae)

Members of the sucking louse genus Pedicinus are ectoparasites of cercopithecid primates in Africa, Asia, and Gibraltar. Pedicinus gabonensis n. sp. is described on the basis of adult male and female specimens collected from the mandrill (Mandrillus sphinx) in Gabon. The new species is compared morphologically with other members of the genus Pedicinus, and a nuclear elongation factor 1 alpha gene sequence is provided. Host associations and geographical distributions of the 18 previously recognized species of the genus and of P. gabonensis n. sp. are reviewed. Updated identification keys are provided for males and females of all known valid species of Pedicinus.

Patterns of host-parasite associations in tropical lice and their passerine hosts in Cameroon

Coevolutionary processes that drive the patterns of host-parasite associations can be deduced through congruence analysis of their phylogenies. Feather lice and their avian hosts have previously been used as typical model systems for congruence analysis; however, such analyses are strongly biased toward nonpasserine hosts in the temperate zone. Further, in the Afrotropical region especially, cospeciation studies of lice and birds are entirely missing. This work supplements knowledge of host-parasite associations in lice using cospeciation analysis of feather lice (genus Myrsidea and the Brueelia complex) and their avian hosts in the tropical rainforests of Cameroon. Our analysis revealed a limited number of cospeciation events in both parasite groups. The parasite-host associations in both louse groups were predominantly shaped by host switching. Despite a general dissimilarity in phylogeny for the parasites and hosts, we found significant congruence in host-parasite distance matrices, mainly driven by associations between Brueelia lice and passerine species of the Waxbill (Estrildidae) family, and Myrsidea lice and their Bulbul (Pycnonotidae) host species. As such, our study supports the importance of complex biotic interactions in tropical environments.

A method to generate multilocus barcodes of pinned insect specimens using MiSeq

For molecular insect identification, amplicon sequencing methods are recommended because they offer a cost-effective approach for targeting small sets of informative genes from multiple samples. In this context, high-throughput multilocus amplicon sequencing has been achieved using the MiSeq Illumina sequencing platform. However, this approach generates short gene fragments of <500 bp, which then have to be overlapped using bioinformatics to achieve longer sequence lengths. This increases the risk of generating chimeric sequences or leads to the formation of incomplete loci. Here, we propose a modified nested amplicon sequencing method for targeting multiple loci from pinned insect specimens using the MiSeq Illumina platform. The modification exists in using a three-step nested PCR approach targeting near full-length loci in the initial PCR and subsequently amplifying short fragments of between 300 and 350 bp for high-throughput sequencing using Illumina chemistry. Using this method, we generated 407 sequences of three loci from 86% of all the specimens sequenced. Out of 103 pinned bee specimens of replicated species, 71% passed the 95% sequence similarity threshold between species replicates. This method worked best for pinned specimens aged between 0 and 5 years, with a limit of 10 years for pinned and 14 years for ethanol-preserved specimens. Hence, our method overcomes some of the challenges of amplicon sequencing using short read next generation sequencing and improves the possibility of creating high-quality multilocus barcodes from insect collections.

Two lineages of kingfisher feather lice exhibit differing degrees of cospeciation with their hosts

Unlike most bird species, individual kingfisher species (Aves: Alcedinidae) are typically parasitized by only a single genus of louse (Alcedoffula, Alcedoecus, or Emersoniella). These louse genera are typically specific to a particular kingfisher subfamily. Specifically, Alcedoecus and Emersoniella parasitize Halcyoninae, whereas Alcedoffula parasitizes Alcedininae and Cerylinae. Although Emersoniella is geographically restricted to the Indo-Pacific region, Alcedoecus and Alcedoffula are geographically widespread. We used DNA sequences from two genes, the mitochondrial COI and nuclear EF-1α genes, to infer phylogenies for the two geographically widespread genera of kingfisher lice, Alcedoffula and Alcedoecus. These phylogenies included 47 kingfisher lice sampled from 11 of the 19 currently recognized genera of kingfishers. We compared louse phylogenies to host phylogenies to reconstruct their cophylogenetic history. Two distinct clades occur within Alcedoffula, one that infests Alcedininae and a second that infests Cerylinae. All species of Alcedoecus were found only on host species of the subfamily Halcyoninae. Cophylogenetic analysis indicated that Alcedoecus, as well as the clade of Alcedoffula occurring on Alcedininae, do not show evidence of cospeciation. In contrast, the clade of Alcedoffula occurring on Cerylinae showed strong evidence of cospeciation.

Host and parasite morphology influence congruence between host and parasite phylogenies

Comparisons of host and parasite phylogenies often show varying degrees of phylogenetic congruence. However, few studies have rigorously explored the factors driving this variation. Multiple factors such as host or parasite morphology may govern the degree of phylogenetic congruence. An ideal analysis for understanding the factors correlated with congruence would focus on a diverse host-parasite system for increased variation and statistical power. In this study, we focused on the Brueelia-complex, a diverse and widespread group of feather lice that primarily parasitise songbirds. We generated a molecular phylogeny of the lice and compared this tree with a phylogeny of their avian hosts. We also tested for the contribution of each host-parasite association to the overall congruence. The two trees overall were significantly congruent, but the contribution of individual associations to this congruence varied. To understand this variation, we developed a novel approach to test whether host, parasite or biogeographic factors were statistically associated with patterns of congruence. Both host plumage dimorphism and parasite ecomorphology were associated with patterns of congruence, whereas host body size, other plumage traits and biogeography were not. Our results lay the framework for future studies to further elucidate how these factors influence the process of host-parasite coevolution.

Elongation Factor-1α Accurately Reconstructs Relationships Amongst Psyllid Families (Hemiptera: Psylloidea), with Possible Diagnostic Implications

The superfamily Psylloidea (Hemiptera: Sternorrhyncha) lacks a robust multigene phylogeny. This impedes our understanding of the evolution of this group of insects and, consequently, an accurate identification of individuals, of their plant host associations, and their roles as vectors of economically important plant pathogens. The conserved nuclear gene elongation factor-1 alpha (EF-1α) has been valuable as a higher-level phylogenetic marker in insects and it has also been widely used to investigate the evolution of intron/exon structure. To explore evolutionary relationships among Psylloidea, polymerase chain reaction amplification and nucleotide sequencing of a 250-bp EF-1α gene fragment was applied to psyllids belonging to five different families. Introns were detected in three individuals belonging to two families. The nine genera belonging to the family Aphalaridae all lacked introns, highlighting the possibility of using intron presence/absence as a diagnostic tool at a family level. When paired with cytochrome oxidase I gene sequences, the 250 bp EF-1α sequence appeared to be a very promising higher-level phylogenetic marker for psyllids.

The roles of barriers, refugia, and chromosomal clines underlying diversification in Atlantic Forest social wasps

Phylogeographic studies have sought to explain the genetic imprints of historical climatic changes and geographic barriers within the Brazilian Atlantic Forest (AF) biota, and consequently two processes of diversification (refugia and barriers) have been proposed. Additionally, there is evidence that eustatic changes influenced the biogeographic history of the AF. Here we evaluate these contrasting diversification processes using two AF social wasp species – the mid-montane Synoeca cyanea and the lowland Synoeca aff. septentrionalis. We analyzed several sources of data including multilocus DNA sequence, climatic niche models and chromosomal features. We find support for idiosyncratic phylogeographic patterns between these wasps, involving different levels of population structure and genetic diversity, contrary suitable climatic conditions during the last glaciation, and contrasting historical movements along the AF. Our data indicate that neotectonics and refugia played distinct roles in shaping the genetic structure of these wasps. However, we argue that eustatic changes influenced the demographic expansion but not population structure in AF biota. Notably, these wasps exhibited chromosomal clines, involving chromosome number and decreasing of GC content, latitudinally oriented along the AF. Together, these results reinforce the need to consider individual organismal histories and indicate that barriers and refugia are significant factors in understanding AF evolution.

Molecular characterization, expression analysis and RNAi knock-down of elongation factor 1α and 1γ from Nilaparvata lugens and its yeast-like symbiont

In the present paper, four cDNAs encoding the alpha and gamma subunits of elongation factor 1 (EF-1) were cloned and sequenced from Nilaparvata lugens, named NlEF-1α, NlEF-1γ, and its yeast-like symbiont (YLS), named YsEF-1α and YsEF-1γ, respectively. Comparisons with sequences from other species indicated a greater conservation for EF-1α than for EF-1γ. NlEF-1α has two identical copies. The deduced amino acid sequence homology of NlEF-1α and NlEF-1γ is 96 and 64%, respectively, compared with Homalodisca vitripennis and Locusta migratoria. The deduced amino acid sequence homology of YsEF-1α and YsEF-1γ is 96 and 74%, respectively, compared with Metarhizium anisopliae and Ophiocordyceps sinensis. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) analysis revealed that the expression level of NlEF-1α and NlEF-1γ mRNA in hemolymph, ovary, fat body and salivary glands were higher than the midgut and leg tissue. YsEF-1α and YsEF-1γ was highly expressed in fat body. The expression level of NlEF-1α was higher than that of NlEF-1γ. Through RNA interference (RNAi) of the two genes, the mortality of nymph reached 92.2% at the 11th day after treatment and the ovarian development was severely hindered. The RT-qPCR analysis verified the correlation between mortality, sterility and the down-regulation of the target genes. The expression and synthesis of vitellogenin (Vg) protein in insects injected with NlEF-1α and NlEF-1γ double-stranded RNA (dsRNA) was significantly lower than control groups. Attempts to knockdown the YsEF-1 genes in the YLS was unsuccessful. However, the phenotype of N. lugens injected with YsEF-1α dsRNA was the same as that injected with NlEF-1α dsRNA, possibly due to the high similarity (up to 71.9%) in the nucleotide sequences between NlEF-1α and YsEF-1α. We demonstrated that partial silencing of NlEF-1α and NlEF-1γ genes caused lethal and sterility effect on N. lugens. NlEF-1γ shares low identity with that of other insects and therefore it could be a potential target for RNAi-based pest management.

EVOLUTION IN A PUTATTVELY ANCIENT ASEXUAL APHID LINEAGE: RECOMBINATION AND RAPID KARYOTYPE CHANGE

Ancient asexual lineages are of great potential significance for understanding the evolutionary biology of sex, but their existence is controversial. In part, this is because claims of ancient asexuality have rested on negative evidence-a mere absence of evidence for sexuality in a taxon. M. Meselson has suggested a method, discussed by Judson and Normark (1996) and by Birky (1996), that has the potential to uncover positive evidence of ancient asexuality. Phylogenetic relationships between alleles and interallelic divergences are predicted to be very different in diploid lineages that lack recombination from those in diploid lineages that undergo recombination. I have applied Meselson's method to the putatively ancient asexual aphid tribe Tramini (Homoptera: Aphidoidea: Lachnidae), using the intron-bearing nuclear protein-coding gene elongation factor 1α (EF-1α). I found heterozygosities much lower than intraspecific divergences, indicating that some recombination has occurred, but not discriminating between recombination within an asexual lineage (automixis or mitotic recombination) and outcrossing sex. Species of Tramini (especially in the genus Trama) typically have highly structurally heterozygous karyotypes that appear to be incompatible with regular successful meiosis, and have very high levels of karyotype variability within species. I found very high levels of karyotype variability within lineages with identical EF-1α and mitochondrial (cytochrome oxidase 1 and 2) genotypes, indicating a high rate of karyotype evolution compared to the rate of nucleotide substitution.

Heterozygous mutation of eEF1A1b resulted in spermatogenesis arrest and infertility in male tilapia, Oreochromis niloticus

Eukaryotic elongation factor 1 alpha (eEF1A) is an essential component of the translational apparatus. In the present study, eEF1A1b was isolated from the Nile tilapia. Real-time PCR and Western blot revealed that eEF1A1b was expressed highly in the testis from 90 dah (days after hatching) onwards. In situ hybridization and immunohistochemistry analyses showed that eEF1A1b was highly expressed in the spermatogonia of the testis. CRISPR/Cas9 mediated mutation of eEF1A1b resulted in spermatogenesis arrest and infertility in the F0 XY fish. Consistently, heterozygous mutation of eEF1A1b (eEF1A1b+/-) resulted in an absence of spermatocytes at 90 dah, very few spermatocytes, spermatids and spermatozoa at 180 dah, and decreased Cyp11b2 and serum 11-ketotestosterone level at both stages. Further examination of the fertilization capacity of the sperm indicated that the eEF1A1b+/- XY fish were infertile due to abnormal spermiogenesis. Transcriptomic analyses of the eEF1A1b+/- testis from 180 dah XY fish revealed that key elements involved in spermatogenesis, steroidogenesis and sperm motility were significantly down-regulated compared with the control XY. Transgenic overexpression of eEF1A1b rescued the spermatogenesis arrest phenotype of the eEF1A1b+/- testis. Taken together, our data suggested that eEF1A1b is crucial for spermatogenesis and male fertility in the Nile tilapia.

High-throughput sequencing of multiple amplicons for barcoding and integrative taxonomy

Until now, the potential of NGS for the construction of barcode libraries or integrative taxonomy has been seldom realised. Here, we amplified (two-step PCR) and simultaneously sequenced (MiSeq) multiple markers from hundreds of fig wasp specimens. We also developed a workflow for quality control of the data. Illumina and Sanger sequences accumulated in the past years were compared. Interestingly, primers and PCR conditions used for the Sanger approach did not require optimisation to construct the MiSeq library. After quality controls, 87% of the species (76% of the specimens) had a valid MiSeq sequence for each marker. Importantly, major clusters did not always correspond to the targeted loci. Nine specimens exhibited two divergent sequences (up to 10%). In 95% of the species, MiSeq and Sanger sequences obtained from the same sampling were similar. For the remaining 5%, species were paraphyletic or the sequences clustered into divergent groups on the Sanger + MiSeq trees (>7%). These problematic cases may represent coding NUMTS or heteroplasms. Our results illustrate that Illumina approaches are not artefact-free and confirm that Sanger databases can contain non-target genes. This highlights the importance of quality controls, working with taxonomists and using multiple markers for DNA-taxonomy or species diversity assessment.

An Assessment of Host Associations, Geographic Distributions, and Genetic Diversity of Avian Chewing Lice (Insecta: Phthiraptera) from Benin

Host associations of highly host-specific chewing lice (Insecta: Phthiraptera) across multiple avian species remains fairly undocumented in the West African country of Benin. Two hundred and seventeen bird specimens collected from multiple localities across Benin and housed at the Texas A&M University Biodiversity Research and Teaching Collections were examined for lice. Lice were identified and genetic data (mitochondrial COI and nuclear EF-1α genes) were obtained and phylogenetically analyzed. In total, we found 15 host associations, 7 of which were new to science. Genetically, most lice from Benin were unique and could represent new species. Based on host associations and unique genetic lineages, we estimate we discovered a minimum of 4 and possibly as many as 8 new chewing louse species. Given the lack of current data on chewing louse species distributions in Benin, this study adds to the knowledge of host associations, geographic distribution, and genetic variability of avian chewing louse species in West Africa.

Genomic Mining of Phylogenetically Informative Nuclear Markers in Bark and Ambrosia Beetles

Deep level insect relationships are generally difficult to resolve, especially within taxa of the most diverse and species rich holometabolous orders. In beetles, the major diversity occurs in the Phytophaga, including charismatic groups such as leaf beetles, longhorn beetles and weevils. Bark and ambrosia beetles are wood boring weevils that contribute 12 percent of the diversity encountered in Curculionidae, one of the largest families of beetles with more than 50000 described species. Phylogenetic resolution in groups of Cretaceous age has proven particularly difficult and requires large quantity of data. In this study, we investigated 100 nuclear genes in order to select a number of markers with low evolutionary rates and high phylogenetic signal. A PCR screening using degenerate primers was applied to 26 different weevil species. We obtained sequences from 57 of the 100 targeted genes. Sequences from each nuclear marker were aligned and examined for detecting multiple copies, pseudogenes and introns. Phylogenetic informativeness (PI) and the capacity for reconstruction of previously established phylogenetic relationships were used as proxies for selecting a subset of the 57 amplified genes. Finally, we selected 16 markers suitable for large-scale phylogenetics of Scolytinae and related weevil taxa.

Patterns of cryptic host specificity in duck lice based on molecular data

Documenting patterns of host specificity in parasites relies on the adequate definition of parasite species. In many cases, parasites have simplified morphology, making species delimitation based on traditional morphological characters difficult. Molecular data can help in assessing whether widespread parasites harbour cryptic species and, alternatively, in guiding further taxonomic revision in cases in which there is morphological variation. The duck louse genus Anaticola (Phthiraptera: Philopteridae), based on current taxonomy, contains both host-specific and widespread species. Mitochondrial and nuclear DNA sequences of samples from this genus were used to document patterns of host specificity. The comparison of these patterns with morphological variations in Anaticola revealed a general correspondence between the groups identified by DNA sequences and morphology, respectively. These results suggest that a more thorough taxonomic review of this genus is needed. In general, the groups identified on the basis of molecular data were associated with particular groups of waterfowl (e.g. dabbling ducks, sea ducks, geese) or specific biogeographic regions (e.g. North America, South America, Australia, Eurasia).

Lineage sorting in multihost parasites: Eidmanniella albescens and Fregatiella aurifasciata on seabirds from the Galapagos Islands

Parasites comprise a significant percentage of the biodiversity of the planet and are useful systems to test evolutionary and ecological hypotheses. In this study, we analyze the effect of host species identity and the immediate local species assemblage within mixed species colonies of nesting seabirds on patterns of genetic clustering within two species of multihost ectoparasitic lice. We use three genetic markers (one mitochondrial, COI, and two nuclear, EF1-α and wingless) and maximum likelihood phylogenetic trees to test whether (1) parasites show lineage sorting based on their host species; and (2) switching of lineages to the alternate host species depends on the immediate local species assemblage of individual hosts within a colony. Specifically, we examine the genetic structure of two louse species: Eidmanniella albescens, infecting both Nazca (Sula granti) and blue-footed boobies (Sula nebouxii), and Fregatiella aurifasciata, infecting both great (Fregata minor) and magnificent frigatebirds (Fregata magnificens). We found that host species identity was the only factor explaining the patterns of genetic structure in both parasites. In both cases, there is evident genetic differentiation depending on the host species. Thus, a revision of the taxonomy of these louse species is needed. One possible explanation of this pattern is extremely low louse migration rates between host species, perhaps influenced by fine-scale spatial separation of host species within mixed colonies, and low parasite infrapopulation numbers.

Changes in the Gene Expression Profiles of the Hypopharyngeal Gland of Worker Honeybees in Association with Worker Behavior and Hormonal Factors

The hypopharyngeal glands (HPGs) of worker honeybees undergo physiological changes along with the age-dependent role change from nursing to foraging: nurse bee HPGs secrete mainly major royal jelly proteins, whereas forager HPGs secrete mainly α-glucosidase III, which converts the sucrose in the nectar into glucose and fructose. We previously identified two other genes, Apis mellifera buffy (Ambuffy) and Apis mellifera matrix metalloproteinase 1 (AmMMP1), with enriched expression in nurse bee and forager HPGs, respectively. In the present study, to clarify the molecular mechanisms that coordinate HPG physiology with worker behavior, we first analyzed whether Ambuffy, AmMMP1, mrjp2 (a gene encoding one of major royal jelly protein isoforms), and Hbg3 (a gene encoding α-glucosidase III) expression, is associated with worker behavior in 'single-cohort colonies' where workers of almost the same age perform different tasks. Expression of these genes correlated with the worker’s role, while controlling for age, indicating their regulation associated with the worker’s behavior. Associated gene expression suggested the possible involvement of some hormonal factors in its regulation. We therefore examined the relationship between ecdysone- and juvenile hormone (JH)-signaling, and the expression profiles of these ‘indicator’ genes (nurse bee HPG-selective genes: mrjp2 and Ambuffy, and forager HPG-selective genes: Hbg3 and AmMMP1). Expression of both ecdysone-regulated genes (ecdysone receptor, mushroom body large type Kenyon cell specific protein-1, and E74) and JH-regulated genes (Methoprene tolerant and Krüppel homolog 1) was higher in the forager HPGs than in the nurse bee HPGs, suggesting the possible roles of ecdysone- and JH-regulated genes in worker HPGs. Furthermore, 20-hydroxyecdysone-treatment repressed both nurse bee- and forager-selective gene expression, whereas methoprene-treatment enhanced the expression of forager-selective genes and repressed nurse bee-selective genes in the HPGs. Our findings suggest that both ecdysone- and JH-signaling cooperatively regulate the physiological state of the HPGs in association with the worker’s behavior.

Polyploidy versus endosymbionts in obligately thelytokous thrips

Background

Thelytoky, the parthenogenetic development of females, has independently evolved in several insect orders yet the study of its mechanisms has so far mostly focussed on haplodiploid Hymenoptera, while alternative mechanisms of thelytoky such as polyploidy are far less understood. In haplodiploid insects, thelytoky can be encoded in their genomes, or induced by maternally inherited bacteria such as Wolbachia or Cardinium. Microbially facilitated thelytoky usually results in complete homozygosity due to gamete duplication and can be reverted into arrhenotoky, the parthenogenetic development of males, through treatment with antibiotics. In contrast, genetically encoded thelytoky cannot be removed and may result in conservation of heterozygosity due to gamete fusion. We have probed the obligate thelytoky of the greenhouse thrips, Heliothrips haemorrhoidalis (Bouché), a significant cosmopolitan pest and a model species of thelytoky in the haplodiploid insect order Thysanoptera. Earlier studies suggested terminal fusion as a mechanism for thelytoky in this species, while another study reported presence of Wolbachia; later it was speculated that Wolbachia plays a role in this thrips’ thelytokous reproduction.

Results

By using PCR and sequence analysis, we demonstrated that global population samples of H. haemorrhoidalis were not infected with Wolbachia, Cardinium or any other known bacterial reproductive manipulators. Antibiotic treatment of this thrips did also not result in male production. Some individuals carried two different alleles in two nuclear loci, histone 3 and elongation factor 1 alpha, suggesting heterozygosity. However, the majority of individuals had three different alleles suggesting that they were polyploid. Genetic diversity across both nuclear loci was low in all populations, and absent from mitochondrial cytochrome oxidase I, indicating that this species had experienced genetic bottlenecks, perhaps due to its invasion biology or a switch to thelytoky.

Conclusions

Geographically broad sampling and experimental manipulation revealed low genetic diversity, absence of Wolbachia but presence of three different alleles of nuclear loci in most analysed individuals of obligately thelytokous H. haemorrhoidalis. This suggests that polyploidy may be involved in the thelytokous reproduction of this thrips species, and polyploidy may be a contributing factor in the reproduction of Thysanoptera and other haplodiploid insect orders.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0304-6) contains supplementary material, which is available to authorized users.

Host generalists and specialists emerging side by side: an analysis of evolutionary patterns in the cosmopolitan chewing louse genus Menacanthus

Parasites with wide host spectra provide opportunities to study the ecological parameters of speciation, as well as the process of the evolution of host specificity. The speciose and cosmopolitan louse genus Menacanthus comprises both multi-host and specialised species, allowing exploration of the ecological and historical factors affecting the evolution of parasites using a comparative approach. We used phylogenetic analysis to reconstruct evolutionary relationships in 14 species of Menacanthus based on the sequences of one mitochondrial and one nuclear gene. The results allowed us to validate species identification based on morphology, as well as to explore host distribution by assumed generalist and specialist species. Our analyses confirmed a narrow host use for several species, however in some cases, the supposed host specialists had a wider host spectrum than anticipated. In one case a host generalist (Menacanthus eurysternus) was clustered terminally on a clade almost exclusively containing host specialists. Such a clade topology indicates that the process of host specialisation may not be irreversible in parasite evolution. Finally, we compared patterns of population genetic structure, geographic distribution and host spectra between two selected species, M. eurysternus and Menacanthus camelinus, using haplotype networks. Menacanthus camelinus showed limited geographical distribution in combination with monoxenous host use, whereas M. eurysternus showed a global distribution and lack of host specificity. It is suggested that frequent host switching maintains gene flow between M. eurysternus populations on unrelated hosts in local populations. However, gene flow between geographically distant localities was restricted, suggesting that geography rather than host-specificity is the main factor defining the global genetic diversity of M. eurysternus.

Phylogenetic position of the bee genera Ancyla and Tarsalia (Hymenoptera: Apidae): a remarkable base compositional bias and an early Paleogene geodispersal from North America to the Old World

We address the phylogenetic position of the bee genera Tarsalia and Ancyla (currently forming the tribe Ancylaini) on the basis of morphological, molecular and combined data. We assembled a matrix of 309 morphological characters and 5246 aligned nucleotide positions from six nuclear genes (28S, EF-1a, wingless, POL2, LW-Rhodopsin, NAK). In addition to both constituent genera of Ancylaini, we include all three subtribes of the Eucerini as well as a large number of other tribes from the "eucerine line". The morphological data suggest Ancyla to be sister to Tarsalia+Eucerini and analyses of the entire molecular dataset suggest Tarsalia to be sister to Ancyla+Eucerini. However, analyses of the combined dataset suggests the Ancylaini to be monophyletic. We address possible bias within the molecular data and show that the base composition of two markers (EF-1a and NAK) is significantly heterogeneous among taxa and that this heterogeneity is strong enough to overcome the phylogenetic signal from the other markers. Analyses of a molecular matrix where the heterogeneous partitions have been RY-recoded yield trees that are better resolved and have higher nodal support values than those recovered in analyses of the non-recoded matrix, and strongly suggest the Ancylaini to be a monophyletic sister group to the Eucerini. A dated phylogeny and ancestral range reconstructions suggest that the common ancestor of the Ancylaini reached the Old World from the New World most probably via the Thulean Land Bridge in a time window between 69 and 47 mya, a period that includes the Early Eocene Climatic Optimum. No further exchanges between the New World and the Old World are implied by our data until the period between 22 mya and 13.9 mya. These more recent faunal exchanges probably involved geodispersal over the Bering Land Bridge by less thermophilic lineages.

Climate, physiological tolerance and sex-biased dispersal shape genetic structure of Neotropical orchid bees

Understanding the impact of past climatic events on the demographic history of extant species is critical for predicting species' responses to future climate change. Palaeoclimatic instability is a major mechanism of lineage diversification in taxa with low dispersal and small geographical ranges in tropical ecosystems. However, the impact of these climatic events remains questionable for the diversification of species with high levels of gene flow and large geographical distributions. In this study, we investigate the impact of Pleistocene climate change on three Neotropical orchid bee species (Eulaema bombiformis, E. meriana and E. cingulata) with transcontinental distributions and different physiological tolerances. We first generated ecological niche models to identify species-specific climatically stable areas during Pleistocene climatic oscillations. Using a combination of mitochondrial and nuclear markers, we inferred calibrated phylogenies and estimated historical demographic parameters to reconstruct the phylogeographical history of each species. Our results indicate species with narrower physiological tolerance experienced less suitable habitat during glaciations and currently exhibit strong population structure in the mitochondrial genome. However, nuclear markers with low and high mutation rates show lack of association with geography. These results combined with lower migration rate estimates from the mitochondrial than the nuclear genome suggest male-biased dispersal. We conclude that despite large effective population sizes and capacity for long-distance dispersal, climatic instability is an important mechanism of maternal lineage diversification in orchid bees. Thus, these Neotropical pollinators are susceptible to disruption of genetic connectivity in the event of large-scale climatic changes.

Molecular phylogeny, biogeography, and host plant shifts in the bee genus Melitta (Hymenoptera: Anthophila)

New molecular studies suggested that the family Melittidae is either a paraphyletic group from which all the other bees are derived, or the sister clade to all other existing bees. Studying the historical biogeography and evolution of each major lineage within this group is a key step to understand the origin and early radiation of bees. Melitta is the largest genus of melittid bees, for which a robust molecular phylogeny and a biogeographic analysis are still lacking. Here, we derive a phylogenetic hypothesis from the sequences of seven independent DNA fragments of mitochondrial and nuclear origin. This phylogenetic hypothesis is then used to infer the evolution of the species range and of the host-plant shifts in Melitta. Our results confirmed the monophyly of Melitta, but did not recover all previously defined clades within the genus. We propose new taxa by splitting the genus in three subgenera (including two new subgenera described in the Appendix: Afromelitta subgen. nov., Plesiomelitta subgen. nov.) and describe two new species: Melitta avontuurensis sp. n. and M. richtersveldensis sp. n. Regarding the evolution of host-plant use, our analysis suggests that all species currently specialized on one plant family originated from an ancestor that was specialized on Fabaceae plants. The inferred biogeographic history for the genus supported an African origin. In concordance with previous studies identifying Africa as the geographic origin for many clades of bees, our data bring new evidence for an African origin of melittid bees.

Elongation factor-1α, a putative single-copy nuclear gene, has divergent sets of paralogs in an arachnid

Identification of paralogy in candidate nuclear loci is an important prerequisite in phylogenetics and statistical phylogeography, but one that is often overlooked. One marker commonly assumed to be a single-copy gene and claimed to harbor great utility for inferring recent divergences is elongation factor-1α (EF-1α). To test this hypothesis, we systematically cloned EF-1α in three disjunct populations of the harvestman Metasiro americanus. Here we show that EF-1α has a large number of paralogs in this species. The paralogs do not evolve in a concerted manner, and the paralogs diverged prior to the population divergence. Moreover, the paralogs of M. americanus are not comparable to the highly divergent EF-1α paralogs found in bees and spiders, which are easily recognized and separated through the use of specific primers. We demonstrate statistically that our detection of paralogs cannot be attributed to amplification error. The presence of EF-1α paralogs in M. americanus prevents its use in statistical phylogeography, and the presence of out-paralogs argues against its use in phylogenetic inference among recently diverged clades. These data contradict the common assumption that EF-1α is for most or all taxa a single-copy gene, or that it has a small number of paralogs that are homogenized through gene conversion, unequal crossing over, or other processes.

Convergent intron gains in hymenopteran elongation factor-1α

The eukaryotic translation elongation factor-1α gene (eEF1A) has been used extensively in higher level phylogenetics of insects and other groups, despite being present in two or more copies in several taxa. Orthology assessment has relied heavily on the position of introns, but the basic assumption of low rates of intron loss and absence of convergent intron gains has not been tested thoroughly. Here, we study the evolution of eEF1A based on a broad sample of taxa in the insect order Hymenoptera. The gene is universally present in two copies - F1 and F2 - both of which apparently originated before the emergence of the order. An elevated ratio of non-synonymous versus synonymous substitutions and differences in rates of amino acid replacements between the copies suggest that they evolve independently, and phylogenetic methods clearly cluster the copies separately. The F2 copy appears to be ancient; it is orthologous with the copy known as F1 in Diptera, and is likely present in most insect orders. The hymenopteran F1 copy, which may or may not be unique to this order, apparently originated through retroposition and was originally intron free. During the evolution of the Hymenoptera, it has successively accumulated introns, at least three of which have appeared at the same position as introns in the F2 copy or in eEF1A copies in other insects. The sites of convergent intron gain are characterized by highly conserved nucleotides that strongly resemble specific intron-associated sequence motifs, so-called proto-splice sites. The significant rate of convergent intron gain renders intron-exon structure unreliable as an indicator of orthology in eEF1A, and probably also in other protein-coding genes.

A total-evidence approach to dating with fossils, applied to the early radiation of the hymenoptera

Phylogenies are usually dated by calibrating interior nodes against the fossil record. This relies on indirect methods that, in the worst case, misrepresent the fossil information. Here, we contrast such node dating with an approach that includes fossils along with the extant taxa in a Bayesian total-evidence analysis. As a test case, we focus on the early radiation of the Hymenoptera, mostly documented by poorly preserved impression fossils that are difficult to place phylogenetically. Specifically, we compare node dating using nine calibration points derived from the fossil record with total-evidence dating based on 343 morphological characters scored for 45 fossil (4--20 complete) and 68 extant taxa. In both cases we use molecular data from seven markers (∼5 kb) for the extant taxa. Because it is difficult to model speciation, extinction, sampling, and fossil preservation realistically, we develop a simple uniform prior for clock trees with fossils, and we use relaxed clock models to accommodate rate variation across the tree. Despite considerable uncertainty in the placement of most fossils, we find that they contribute significantly to the estimation of divergence times in the total-evidence analysis. In particular, the posterior distributions on divergence times are less sensitive to prior assumptions and tend to be more precise than in node dating. The total-evidence analysis also shows that four of the seven Hymenoptera calibration points used in node dating are likely to be based on erroneous or doubtful assumptions about the fossil placement. With respect to the early radiation of Hymenoptera, our results suggest that the crown group dates back to the Carboniferous, ∼309 Ma (95% interval: 291--347 Ma), and diversified into major extant lineages much earlier than previously thought, well before the Triassic. [Bayesian inference; fossil dating; morphological evolution; relaxed clock; statistical phylogenetics.].

Multilocus genetic characterization of two ant vectors (Group II "Dirty 22" species) known to contaminate food and food products and spread foodborne pathogens

The U.S. Food and Drug Administration utilizes the presence of filth and extraneous materials as one of the criteria for implementing regulatory actions and assessing adulteration of food products of public health importance. Twenty-two prevalent pest species (also known as the ''Dirty 22'' species) have been considered by this agency as possible vehicles for the spread of foodborne diseases, and the presence of these species is considered an indicator of unsanitary conditions in food processing and storage facilities. In a previous study, we further categorized the Dirty 22 species into four groups: group I includes four cockroach species, group II includes two ant species, group III includes 12 fly species, and group IV includes four rodent species. Here, we describe the development of three nested PCR primer sets and multilocus genetic characterization by amplifying the small subunit rRNA, elongation factor 1-alpha, and wingless (WNT-1) genes of group II Dirty 22 ant species Monomorium pharaonis and Solenopsis molesta. These novel group II Dirty 22 species-specific nested PCR primer sets can be used when the specimens cannot be identified using conventional microscopic methods. These newly developed assays will provide correct identification of group II Dirty 22 ant species, and the information can be used in the control of foodborne pathogens.

Molecular discrimination of phytoseiids associated with the red palm mite Raoiella indica (Acari: Tenuipalpidae) from Mauritius and South Florida

Phytoseiid populations imported from Mauritius for evaluation for a classical biological control program in Florida, USA, were morphologically identified as Amblyseius largoensis Muma, a species associated with the red palm mite in south Florida and the Caribbean. Bayesian analysis and sequence divergences of the mitochondrial 12S rRNA and nuclear Elongation factor--I alpha (EF-Iα) genes and Neighbor-Joining analysis of High-fidelity-RAPD-PCR markers were used to discriminate between the south Florida and Mauritius populations. High-fidelity-RAPD-PCR markers in addition to Bayesian and sequence divergence analyses of the 12S rRNA sequences suggest that the Mauritius and south Florida populations are genetically different but whether these are species or population differences is unknown. The degenerate EF-Iα primers used to survey the phytoseiids amplified two different elongation factor sequences with distinct amino acid translations, the putative EF-Iα and an unknown elongation factor. Variability within the 12S gene was used to develop population-specific primers for identifying the Mauritius phytoseiids in the event they are released in south Florida.

Repeated adaptive divergence of microhabitat specialization in avian feather lice

Background

Repeated adaptive radiations are evident when phenotypic divergence occurs within lineages, but this divergence into different forms is convergent when compared across lineages. Classic examples of such repeated adaptive divergence occur in island (for example, Caribbean Anolis lizards) and lake systems (for example, African cichlids). Host-parasite systems in many respects are analogous to island systems, where host species represent isolated islands for parasites whose life cycle is highly tied to that of their hosts. Thus, host-parasite systems might exhibit interesting cases of repeated adaptive divergence as seen in island and lake systems.

The feather lice of birds spend their entire life cycle on the body of the host and occupy distinct microhabitats on the host: head, wing, body and generalist. These microhabitat specialists show pronounced morphological differences corresponding to how they escape from host preening. We tested whether these different microhabitat specialists were a case of repeated adaptive divergence by constructing both morphological and molecular phylogenies for a diversity of avian feather lice, including many examples of head, wing, body and generalist forms.

Results

Morphological and molecular based phylogenies were highly incongruent, which could be explained by rampant convergence in morphology related to microhabitat specialization on the host. In many cases lice from different microhabitat specializations, but from the same group of birds, were sister taxa.

Conclusions

This pattern indicates a process of repeated adaptive divergence of these parasites within host group, but convergence when comparing parasites across host groups. These results suggest that host-parasite systems might be another case in which repeated adaptive radiations could be relatively common, but potentially overlooked, because morphological convergence can obscure evolutionary relationships.

Cryptic variation in an ecological indicator organism: mitochondrial and nuclear DNA sequence data confirm distinct lineages of Baetis harrisoni Barnard (Ephemeroptera: Baetidae) in southern Africa

BACKGROUND

Baetis harrisoni Barnard is a mayfly frequently encountered in river studies across Africa, but the external morphological features used for identifying nymphs have been observed to vary subtly between different geographic locations. It has been associated with a wide range of ecological conditions, including pH extremes of pH 2.9-10.0 in polluted waters. We present a molecular study of the genetic variation within B. harrisoni across 21 rivers in its distribution range in southern Africa.

RESULTS

Four gene regions were examined, two mitochondrial (cytochrome c oxidase subunit I [COI] and small subunit ribosomal 16S rDNA [16S]) and two nuclear (elongation factor 1 alpha [EF1α] and phosphoenolpyruvate carboxykinase [PEPCK]). Bayesian and parsimony approaches to phylogeny reconstruction resulted in five well-supported major lineages, which were confirmed using a general mixed Yule-coalescent (GMYC) model. Results from the EF1α gene were significantly incongruent with both mitochondrial and nuclear (PEPCK) results, possibly due to incomplete lineage sorting of the EF1α gene. Mean between-clade distance estimated using the COI and PEPCK data was found to be an order of magnitude greater than the within-clade distance and comparable to that previously reported for other recognised Baetis species. Analysis of the Isolation by Distance (IBD) between all samples showed a small but significant effect of IBD. Within each lineage the contribution of IBD was minimal. Tentative dating analyses using an uncorrelated log-normal relaxed clock and two published estimates of COI mutation rates suggest that diversification within the group occurred throughout the Pliocene and mid-Miocene (~2.4-11.5 mya).

CONCLUSIONS

The distinct lineages of B. harrisoni correspond to categorical environmental variation, with two lineages comprising samples from streams that flow through acidic Table Mountain Sandstone and three lineages with samples from neutral-to-alkaline streams found within eastern South Africa, Malawi and Zambia. The results of this study suggest that B. harrisoni as it is currently recognised is not a single species with a wide geographic range and pH-tolerance, but may comprise up to five species under the phylogenetic species concept, each with limited pH-tolerances, and that the B. harrisoni species group is thus in need of taxonomic review.

Over-expression of a cytochrome P450 is associated with resistance to pyriproxyfen in the greenhouse whitefly Trialeurodes vaporariorum

Background

The juvenile hormone mimic, pyriproxyfen is a suppressor of insect embryogenesis and development, and is effective at controlling pests such as the greenhouse whitefly Trialeurodes vaporariorum (Westwood) which are resistant to other chemical classes of insecticides. Although there are reports of insects evolving resistance to pyriproxyfen, the underlying resistance mechanism(s) are poorly understood.

Results

Bioassays against eggs of a German (TV8) population of T. vaporariorum revealed a moderate level (21-fold) of resistance to pyriproxyfen. This is the first time that pyriproxyfen resistance has been confirmed in this species. Sequential selection of TV8 rapidly generated a strain (TV8pyrsel) displaying a much higher resistance ratio (>4000-fold). The enzyme inhibitor piperonyl butoxide (PBO) suppressed this increased resistance, indicating that it was primarily mediated via metabolic detoxification. Microarray analysis identified a number of significantly over-expressed genes in TV8pyrsel as candidates for a role in resistance including cytochrome-P450 dependent monooxygenases (P450s). Quantitative PCR highlighted a single P450 gene (CYP4G61) that was highly over-expressed (81.7-fold) in TV8pyrsel.

Conclusion

Over-expression of a single cytochrome P450 gene (CYP4G61) has emerged as a strong candidate for causing the enhanced resistance phenotype. Further work is needed to confirm the role of the encoded P450 enzyme CYP4G61 in detoxifying pyriproxyfen.

Cophylogeny on a fine scale: Geomydoecus chewing lice and their pocket gopher hosts, Pappogeomys bulleri

Many species of pocket gophers and their ectoparasitic chewing lice have broadly congruent phylogenies, indicating a history of frequent codivergence. For a variety of reasons, phylogenies of codiverging hosts and parasites are expected to be less congruent for more recently diverged taxa. This study is the first of its scale in the pocket gopher and chewing louse system, with its focus entirely on comparisons among populations within a single species of host and 3 chewing louse species in the Geomydoecus bulleri species complex. We examined mitochondrial DNA from a total of 46 specimens of Geomydoecus lice collected from 11 populations of the pocket gopher host, Pappogeomys bulleri. We also examined nuclear DNA from a subset of these chewing lice. Louse phylogenies were compared with a published pocket gopher phylogeny. Contrary to expectations, we observed a statistically significant degree of parallel cladogenesis in these closely related hosts and their parasites. We also observed a higher rate of evolution in chewing louse lineages than in their corresponding pocket gopher hosts. In addition, we found that 1 louse species (Geomydoecus burti) may not be a valid species, that subspecies within G. bulleri are not reciprocally monophyletic, and that morphological and genetic evidence support recognition of a new species of louse, Geomydoecus pricei.