Reconstructing species phylogeny of the carabid beetles Ohomopterus using multiple nuclear DNA sequences: heterogeneous information content and the performance of simultaneous analyses

We attempted a phylogenetic reconstruction for the carabid subgenus Ohomopterus (genus Carabus), a notable case of radiation with mitochondrial introgression across species. Sequence data from five nuclear single copy loci were used, including wingless (Wg), phosphoenolpyruvate carboxykinase (PepCK), cytochrome c (Cytc), elongation factor-1alpha (EF-1alpha), and an anonymous single copy locus (Carab1). Sequences of Cytc, EF-1alpha, and Carab1 included intron or intron-like parts with length variation. The analysis of individual loci resulted in low resolution of the phylogenetic relationships, and the monophyly of several morphologically recognized species for which multiple specimens were analyzed was not revealed. Several specimens were heterozygous, with non-monophyletic alleles observed in three of the five loci at which alleles in heterozygotes were separated. In a simultaneous analysis of the five loci with ambiguously aligned parts eliminated and heterozygotic sites treated as missing, the resulting tree was well resolved, but the branch support was generally weak because of conflicting phylogenetic signals from different loci. We also attempted to incorporate allelic sequence data plus the ambiguously aligned parts in the analysis, by using all possible combinations of alleles from different loci in heterozygotic individuals, but the resultant tree was not supported more strongly. Nonetheless, these simultaneous analyses provided support for the monophyly of several species and species groups, and revealed the basic evolutionary trend of OHOMOPTERUS: initial widespread groups with simpler genitalia and the origination of exaggerated genitalia in a derived clade. This study exemplifies problems inherent in the phylogenetic reconstruction of closely related organisms where low levels of variation limit the information content from each locus, while heterozygosity, different phylogenetic history of multiple loci, and alignment ambiguity further hamper phylogenetic reconstruction unless several loci converge on a uniform signal.     

Phylogenetic relationships of the carabid subfamily Harpalinae (Coleoptera) based on molecular sequence data

The carabid subfamily Harpalinae contains most of the species of carabid beetles. This subfamily, with over 19,000 species, radiated in the Cretaceous to yield a large clade that is diverse in morphological form and ecological habit. While there are several morphological, cytological, and chemical characters that unite most harpalines, the placement of some tribes within the subfamily remains controversial, as does the sister group relationships to this large group. In this study, DNA sequences from the 28S rDNA gene and the wingless nuclear protein-coding gene were collected from 52 carabid genera representing 31 harpaline tribes in addition to more than 21 carabid outgroup taxa to reconstruct the phylogeny of this group. Molecular sequence data from these genes, along with additional data from the 18S rDNA gene, were analyzed with a variety of phylogenetic analysis methods, separately for each gene and in a combined data approach. Results indicated that the subfamily Harpalinae is monophyletic with the enigmatic tribes of Morionini, Peleciini, and Pseudomorphini included within it. Brachinine bombardier beetles are closely related to Harpalinae as they form the sister group to harpalines or, in some analyses, are included within it or with austral psydrines. The austral psydrines are the sister group to Harpalinae+Brachinini clade in most analyses and austral psydrines+Brachinini+Harpalinae clade is strongly supported.     

Widespread genealogical nonmonophyly in species of Pinus subgenus Strobus

Phylogenetic relationships among Pinus species from subgenus Strobus remain unresolved despite combined efforts based on nrITS and cpDNA. To provide greater resolution among these taxa, a 900-bp intron from a late embryogenesis abundant (LEA)-like gene (IFG8612)was sequenced from 39 pine species, with two or more alleles representing 33 species. Nineteen of 33 species exhibited allelic nonmonphyly in the strict consensus tree, and 10 deviated significantly from allelic monophyly based on topology incongruence tests. Intraspecific nucleotide diversity ranged from 0.0 to 0.0211, and analysis of variance shows that nucleotide diversity was strongly associated (P < 0.0001)with the degree of species monophyly. Although species nonmonophyly complicates phylogenetic interpretations, this nuclear locus offers greater topological support than previously observed for cpDNA or nrITS. Lacking evidence for hybridization, recombination, or imperfect taxonomy, we feel that incomplete lineage sorting remains the best explanation for the polymorphisms shared among species. Depending on the species, coalescent expectations indicate that reciprocal monophyly will be more likely than paraphyly in 1.71 to 24.0 x 10(6) years, and that complete genome-wide coalescence in these species may require up to 76.3 x 10(6) years. The absence of allelic coalescence is a severe constraint in the application of phylogenetic methods in Pinus, and taxa sharing similar life history traits with Pinus are likely to show species nonmonophyly using nuclear markers.     

Incorporating allelic variation for reconstructing the evolutionary history of organisms from multiple genes: An example from Rosa in North America

Allelic variation within individuals holds information regarding the relationships of organisms, which is expected to be particularly important for reconstructing the evolutionary history of closely related taxa. However, little effort has been committed to incorporate such information for reconstructing the phylogeny of organisms. Haplotype trees represent a solution when one nonrecombinant marker is considered, but there is no satisfying method when multiple genes are to be combined. In this paper, we propose an algorithm that converts a distance matrix of alleles to a distance matrix among organisms. This algorithm allows the incorporation of allelic variation for reconstructing the phylogeny of organisms from one or more genes. The method is applied to reconstruct the phylogeny of the seven native diploid species of Rosa sect. Cinnamomeae in North America. The glyceralgehyde 3-phosphate dehydrogenase (GAPDH), the triose phosphate isomerase (TPI), and the malate synthase (MS) genes were sequenced for 40 individuals from these species. The three genes had little genetic variation, and most species showed incomplete lineage sorting, suggesting these species have a recent origin. Despite these difficulties, the networks (NeighborNet) of organisms reconstructed from the matrix obtained with the algorithm recovered groups that more closely match taxonomic boundaries than did the haplotype trees. The combined network of individuals shows that species west of the Rocky Mountains, Rosa gymnocarpa and R. pisocarpa, form exclusive groups and that together they are distinct from eastern species. In the east, three groups were found to be exclusive: R. nitida-R. palustris, R. foliolosa, and R. blanda-R. woodsii. These groups are congruent with the morphology and the ecology of species. The method is also useful for representing hybrid individuals when the relationships are reconstructed using a phylogenetic network.     

Relationships among pest flour beetles of the genus Tribolium (Tenebrionidae) inferred from multiple molecular markers

Model species often provide initial hypotheses and tools for studies of development, genetics and molecular evolution in closely related species. Flour beetles of the genus Tribolium Macleay (1825) are one group with potential for such comparative studies. Tribolium castaneum (Herbst 1797) is an increasingly useful developmental genetic system. The convenience with which congeneric and other species of tenebrionid flour beetles can be reared in the laboratory makes this group attractive for comparative studies on a small phylogenetic scale. Here we present the results of phylogenetic analyses of relationships among the major pest species of Tribolium based on two mitochondrial and three nuclear markers (cytochrome oxidase 1, 16S ribosomal DNA, wingless, 28S ribosomal DNA and histone H3). The utility of partitioning the dataset in a manner informed by biological structure and function is demonstrated by comparing various partitioning strategies. In parsimony and partitioned Bayesian analyses of the combined dataset, the castaneum and confusum species groups are supported as monophyletic and as each other's closest relatives. However, a sister group relationship between this clade and Tribolium brevicornis (Leconte 1859) is not supported. The inferred phylogeny provides an evolutionary framework for comparative studies using flour beetles.     

Phylogenetic discordance at the species boundary: comparative gene genealogies among rapidly radiating Heliconius butterflies

Recent adaptive radiations provide excellent model systems for understanding speciation, but rapid diversification can cause problems for phylogenetic inference. Here we use gene genealogies to investigate the phylogeny of recent speciation in the heliconiine butterflies. We sequenced three gene regions, intron 3 ( approximately 550 bp) of sex-linked triose-phosphate isomerase (Tpi), intron 3 ( approximately 450 bp) of autosomal mannose-phosphate isomerase (Mpi), and 1,603 bp of mitochondrial cytochrome oxidase subunits I and II (COI and COII), for 37 individuals from 25 species of Heliconius and related genera. The nuclear intron sequences evolved at rates similar to those of mitochondrial coding sequences, but the phylogenetic utility of introns was restricted to closely related geographic populations and species due to high levels of indel variation. For two sister species pairs, Heliconius erato-Heliconius himera and Heliconius melpomene-Heliconius cydno, there was highly significant discordance between the three genes. At mtDNA and Tpi, the hypotheses of reciprocal monophyly and paraphyly of at least one species with respect to its sister could not be distinguished. In contrast alleles sampled from the third locus, Mpi, showed polyphyletic relationships between both species pairs. In all cases, recent coalescence of mtDNA lineages within species suggests that polyphyly of nuclear genes is not unexpected. In addition, very similar alleles were shared between melpomene and cydno, implying recent gene flow. Our finding of discordant genealogies between genes is consistent with models of adaptive speciation with ongoing gene flow and highlights the need for multiple locus comparisons to resolve phylogeny among closely related species.     

Phylogenetic relationships among holarctic populations of Chironomus entis and Chironomus plumosus in view of possible horizontal transfer of mitochondrial genes

In eight Holarctic populations of two typical chironomid sibling species of the plumosus group, Chrionomus entis and Chironomus plumosus, nucleotides sequences of mitochondrial (cytb) and nuclear (gb2b) gene regions were examined. The phylogenetic trees reflecting the evolutionary histories of the nuclear and mitochondrial markers exhibited significant differences. On the tree based on the nuclear gene sequences the populations clustered according to their species affiliation, whereas on the tree based on the mitochondrial gene sequences the populations were grouped according to their geographic position. This discrepancy is probably explained by mitochondrial gene flow between sympatric species with incomplete reproductive isolation (sibling species). Based on our results together with the earlier data on nuclear and mitochondrial gene sequences of some other species from the phylogenetic group plumosus, a scheme of phylogenetic relationships within this group is proposed. This scheme is in many ways different from the traditional view on the evolutionary relationships among species of the plumosus group.     

Mechanical barriers to introgressive hybridization revealed by mitochondrial introgression patterns in Ohomopterus ground beetle assemblages

To reveal the role of diverged body size and genital morphology in reproductive isolation among closely related species, we examined patterns of, and factors limiting, introgressive hybridization between sympatric Ohomopterus ground beetles in central Japan using mitochondrial NADH dehydrogenase subunit 5 (ND5) gene sequences. We sampled 17 local assemblages that consisted of two to five species and estimated levels of interspecific gene flow using the genetic distance, D(A), and maximum-likelihood estimates of gene flow. Sharing of haplotypes or haplotype lineages was detected between six of seven species that occurred in the study areas, indicating mitochondrial introgression. The intensity and direction of mitochondrial gene flow were variable among species pairs. To determine the factors affecting introgression patterns, we tested the relationships between interspecific D(A) and five independent variables: difference in body size, difference in genital size, phylogenetic relatedness (nuclear gene sequence divergence), habitat difference, and species richness of the assemblage. Body and genital size differences contributed significantly to preventing gene flow. Thus, mechanical isolation mechanisms reduce the chance of introgressive hybridization between closely related species. Our results highlight the role of morphological divergence in speciation and assemblage formation processes through mechanical isolation.     

Interspecific divergence, intrachromosomal recombination, and phylogenetic utility of Y-chromosomal genes in Drosophila

Reconstruction of phylogenetic relationships among recently diverged species is complicated by three general problems: segregation of polymorphisms that pre-date species divergence, gene flow during and after speciation, and intra-locus recombination. In light of these difficulties, the Y chromosome offers several important advantages over other genomic regions as a source of phylogenetic information. These advantages include the absence of recombination, rapid coalescence, and reduced opportunity for interspecific introgression due to hybrid male sterility. In this report, we test the phylogenetic utility of Y-chromosomal sequences in two groups of closely related and partially inter-fertile Drosophila species. In the D. bipectinata species complex, Y-chromosomal loci unambiguously recover the phylogeny most consistent with previous multi-locus analysis and with reproductive relationships, and show no evidence of either post-speciation gene flow or persisting ancestral polymorphisms. In the D. simulans species complex, the situation is complicated by the duplication of at least one Y-linked gene region, followed by intrachromosomal recombination between the duplicate genes that scrambles their genealogy. We suggest that Y-chromosomal sequences are a useful tool for resolving phylogenetic relationships among recently diverged species, especially in male-heterogametic organisms that conform to Haldane's rule. However, duplication of Y-linked genes may not be uncommon, and special care should be taken to distinguish between orthologous and paralogous sequences.     

Genomes,multiple origins,and lineage recombination in the Glycine tomentella (Leguminosae) polyploid complex: histone H3-D gene sequences

Relationships among the various diploid and polyploid taxa that comprise Glycine tomentella have been hypothesized from crossing studies, isozyme data, and repeat length variation for the 5S nuclear ribosomal gene loci. However, several key questions have persisted, and detailed phylogenetic evidence from homoeologous nuclear genes has been lacking. The histone H3-D locus is single copy in diploid Glycine species and has been used to elucidate relationships among diploid races of G. tomentella, providing a framework for testing genome origins in the polyploid complex. For all six G. tomentella polyploid races (T1-T6), alleles at two homoeologous histone H3-D loci were isolated and analyzed phylogenetically with alleles from diploid Glycine species, permitting the identification of all of the homoeologous genomes of the complex. Allele networks were constructed to subdivide groups of homoeologous alleles further, and two-locus genotypes were constructed using these allele classes. Results suggest that some races have more than one origin and that interfertility within races has led to lineage recombination. Most alleles in polyploids are identical or closely related to alleles in diploids, suggesting recency of polyploid origins and spread beyond Australia. These features parallel the other component of the Glycine subgenus Glycine polyploid complex, G. tabacina, one of whose races shares a diploid genome with a G. tomentella polyploid race.     

Trapped in freshwater: the internal anatomy of the entoproct Loxosomatoides sirindhornae

Background

The identification of vast numbers of unknown organisms using DNA sequences becomes more and more important in ecological and biodiversity studies. In this context, a fragment of the mitochondrial cytochrome c oxidase I (COI) gene has been proposed as standard DNA barcoding marker for the identification of organisms. Limitations of the COI barcoding approach can arise from its single-locus identification system, the effect of introgression events, incomplete lineage sorting, numts, heteroplasmy and maternal inheritance of intracellular endosymbionts. Consequently, the analysis of a supplementary nuclear marker system could be advantageous.

Results

We tested the effectiveness of the COI barcoding region and of three nuclear ribosomal expansion segments in discriminating ground beetles of Central Europe, a diverse and well-studied invertebrate taxon. As nuclear markers we determined the 18S rDNA: V4, 18S rDNA: V7 and 28S rDNA: D3 expansion segments for 344 specimens of 75 species. Seventy-three species (97%) of the analysed species could be accurately identified using COI, while the combined approach of all three nuclear markers provided resolution among 71 (95%) of the studied Carabidae.

Conclusion

Our results confirm that the analysed nuclear ribosomal expansion segments in combination constitute a valuable and efficient supplement for classical DNA barcoding to avoid potential pitfalls when only mitochondrial data are being used. We also demonstrate the high potential of COI barcodes for the identification of even closely related carabid species.     

Positive selection and interallelic recombination at the merozoite surface antigen-1 (MSA-1) locus of Plasmodium falciparum.

DNA sequences of alleles at the merozoite surface antigen-1 (MSA-1) gene locus of the malaria parasite Plasmodium falciparum show evidence of repeated past recombination events between alleles. These include both (1) nonreciprocal recombination events that have homogenized certain gene regions among alleles and (2) reciprocal recombination events that have combined allelic segments with divergent evolutionary histories, thereby enhancing allelic diversity. In three different gene regions, the rate of nonsynonymous nucleotide substitution significantly exceeds that of synonymous nucleotide substitution, implying that positive Darwinian selection has acted to diversify alleles at the amino acid level. The MSA-1 polymorphism seems to be quite ancient; the two major allelic types have been maintained for approximately 35 Myr.     

Phylogeny vs genome reshuffling: horizontal gene transfer

The evolutionary events in organisms can be tracked to the transfer of genetic material. The inheritance of genetic material among closely related organisms is a slow evolutionary process. On the other hand, the movement of genes among distantly related species can account for rapid evolution. The later process has been quite evident in the appearance of antibiotic resistance genes among human and animal pathogens. Phylogenetic trees based on such genes and those involved in metabolic activities reflect the incongruencies in comparison to the 16S rDNA gene, generally used for taxonomic relationships. Such discrepancies in gene inheritance have been termed as horizontal gene transfer (HGT) events. In the post-genomic era, the explosion of known sequences through large-scale sequencing projects has unraveled the weakness of traditional 16S rDNA gene tree based evolutionary model. Various methods to scrutinize HGT events include atypical composition, abnormal sequence similarity, anomalous phylogenetic distribution, unusual phyletic patterns, etc. Since HGT generates greater genetic diversity, it is likely to increase resource use and ecosystem resilience.     

Phylogeny of the superfamily Gelechioidea (Lepidoptera: Obtectomera), with an exploratory application on geometric morphometrics

The superfamily Gelechioidea (Lepidoptera: Obtectomera) has a high species diversity. It consists of more than 18,400 described species and has a global distribution. Among it, large numbers of species were reported to be economically important to people's production and life. However, relationships among families or subfamilies in Gelechioidea have been exceptionally difficult to resolve using morphology or single gene genealogies. Multiple gene genealogies had been used in the molecular phylogenetic studies on Gelechioidea during the past years, but their phylogenetic relationships remain to be controversial mainly due to their limited taxa sampling relative to such high species diversity. In this paper, 89 ingroup species representing 55 genera are sequenced and added to the data downloaded from GenBank, and six species representing four closely related superfamilies are chosen as outgroup. The molecular phylogeny of Gelechioidea is reconstructed based on the concatenated data set composed of one mitochondrial marker (COI) and seven nuclear markers (CAD, EF-1ɑ, GAPDH, IDH, MDH, RpS5, wingless). The phylogenetic results, taking into consideration of the comparative morphological study, show that the clade of Gelechioidea is strongly supported and separated from other superfamilies, which further proves its monophyly. Five families are newly defined: Autostichidae sensu nov., Depressariidae sensu nov., Peleopodidae sensu nov., Ashinagidae sensu nov. and Epimarptidae sensu nov. Meanwhile, a monophyletic “SSABM” clade considered to be closely related is proposed for the first time, consisting of Stathmopodidae, Scythrididae, Ashinagidae, Blastobasidae and Momphidae. Moreover, geometric morphometric analyses using merged landmark data set from fore and hind wings of 118 representative species are conducted. The phenetic tree shows that the monophyly and phylogenetic relationships correspond with the results of molecular phylogeny largely, which well proves its importance and potential application in both phylogenetic reconstruction and species identification.     

Reevaluation of phylogeny in the Drosophila obscura species group based on combined analysis of nucleotide sequences.

The Drosophila obscura species group has served as an important model system in many evolutionary and population genetic studies. Despite the amount of study this group has received, some phylogenetic relationships remain unclear. While individual analysis of different nuclear, mitochondrial, allozyme, restriction fragment, and morphological data partitions are able to discern relationships among closely related species, they are unable to resolve relationships among the five obscura species subgroups. A combined analysis of several nucleotide data sets is able to provide resolution and support for some nodes not seen or well supported in analyses of individual loci. A phylogeny of the obscura species group based on combined analysis of nucleotide sequences from six mitochondrial and five nuclear loci is presented here. The results of several different combined analyses indicate that the Old World obscura and subobscura subgroups form a monophyletic clade, although they are unable to resolve the relationships among the major lineages within the obscura species group.     

Nuclear gene sequences resolve species phylogeny and mitochondrial introgression in Leptocarabus beetles showing trans-species polymorphisms

We studied the phylogenetic relationships among Japanese Leptocarabus ground beetles, which show extensive trans-species polymorphisms in mitochondrial gene genealogies. Simultaneous analysis of combined nuclear data with partial sequences from the long-wavelength rhodopsin, wingless, phosphoenolpyruvate carboxykinase, and 28S rRNA genes resolved the relationships among the five species, although separate analyses of these genes provided topologies with low resolution. For both the nuclear gene tree resulting from the combined data from four genes and a mitochondrial cytochrome oxidase subunit I (COI) gene tree, we applied a Bayesian divergence time estimation using a common calibration method to identify mitochondrial introgression events that occurred after speciation. Three mitochondrial lineages shared by two or three species were likely subject to introgression due to interspecific hybridization because the coalescent times for these lineages were much shorter than the corresponding speciation times estimated from nuclear gene sequences. We demonstrated that when species phylogeny is fully resolved with nuclear gene sequence data, comparative analysis of nuclear and mitochondrial gene trees can be used to infer introgressive hybridization events that might cause trans-species polymorphisms in mitochondrial gene trees.     

Phylogenetic Relationships among Holarctic Populations of Chironomus entis and Chironomus plumosus in View of Possible Horisontal Transfer of Mitochondrial Genes

In eight Holarctic populations of two typical chironomid sibling species of the plumosus group, Chironomus entisandChironomus plumosus, nucleotide sequences of mitochondrial (cytb) and nuclear (gb2b) gene regions were examined. The phylogenetic trees reflecting the evolutionary histories of the nuclear and mitochondrial markers exhibited significant differences. On the tree based on the nuclear gene sequences the populations clustered according to their species affiliation, whereas on the tree based on the mitochondrial gene sequences the populations were grouped according to their geographic position. This discrepancy is probably explained by mitochondrial gene flow between sympatric species with incomplete reproductive isolation (sibling species). Based on our results together with the earlier data on nuclear and mitochondrial gene sequences of some other species from the phylogenetic group plumosus, a scheme of phylogenetic relationships within this group is proposed. This scheme is in many ways different from the traditional view on the evolutionary relationships among species of the plumosus group.     

Influence of tree species and soil properties on ground beetle (Coleoptera: Carabidae) communities

Although Carabidae is among the best-studied families of beetles in Europe from the faunistic point of view, there is still a lack of available information on the ecological requirements of the particular carabid species. The habitat preferences that determine the distribution of species are largely influenced by habitat structure and microclimate. In addition to other factors, these habitat parameters are influenced by the nature of the vegetation. Therefore, our study investigated the influence of tree species on carabid beetle communities. We conducted the research at 9 stands in the Borová Hora Arboretum (Zvolen, Central Slovakia). Each studied site represents a monoculture of one of nine tree species. At each site, some soil and leaf litter attributes (pH, conductivity, and content of H, C, N and P) were evaluated. Ground beetles were collected by pitfall trapping during the vegetation periods in 2008–2011. In total, 3012 individuals of 29 species were obtained. Significant differences in the total dynamic activity and species richness of the carabid beetle communities among the compared forest stands were revealed. The results of the research confirmed statistically significant relationships among 1) the soil conductivity and both the richness and Shannon diversity of the ground beetle communities, 2) the litter and soil N content and richness, the Shannon diversity and the species composition of the ground beetle communities. The Shannon diversity and richness were negatively related to the soil conductivity and positively related with the N content. Our research showed that dominant tree species indirectly influence diversity and composition of carabid communities via the soil properties.     

An empirical demonstration of using pentatricopeptide repeat (PPR) genes as plant phylogenetic tools: Phylogeny of Verbenaceae and the Verbena complex

The pentatricopeptide repeat (PPR) gene family, with hundreds of members in land plant genomes, has been recognized as a tremendous resource for plant phylogenetic studies based on publicly available genomic data from model organisms. However, whether this appealing nuclear gene marker system can be readily applied to non-model organisms remains questionable, particularly given the potential uncertainties in designing specific primers to only amplify the locus of interest from the sea of PPR genes. Here we demonstrate empirically the use of PPR genes in the family Verbenaceae and the Verbena complex. We also lay out a general scheme to design locus-specific primers to amplify and sequence PPR genes in non-model organisms. Intergeneric relationships within the family Verbenaceae were fully resolved with strong support. Relationships among the closely related genera within the Verbena complex and among some species groups within each genus were also well resolved, but resolution among very closely related species was limited. Our results suggest that PPR genes can be readily employed in non-model organisms. They may be best used to resolve relationships in a spectrum from among distantly related genera to among not-so-closely related congeneric species, but may have limited use among very closely related species.