Exon-primed Intron-crossing Research Articles

The genus Coenagrion Kirby, 1890 (Odonata: Coenagrionidae) in the Russian part of the Caucasus

Distribution, characters and habitats of Coenagrion spp. in the Russian part of the Caucasus is considered, chiefly by the authors’ own data. There are seven species, of which C. australocaspicum Dumont & Heidari, 1996 is for the first time recorded in Russia (Dagestan), C. ornatum (Selys, 1850) and C. scitulum (Rambur, 1842) in Karachay-Cherkess and Dagestan Republics and C. pulchellum (Vander Linden, 1825) in Dagestan. Coenagrion lunulatum is still known from the territory only by a century-old record by A.N. Bartenev from the highlands of West Caucasus. Keys for identification of the Caucasian species of the puella-group in a broad sense—C. australocaspicum, C. puella (Linnaeus, 1758), C. pulchellum and C. ponticum (Bartenef, 1929)—are provided. Coenagrion ponticum is supposed to be a hybridogenic species resulting from hybridiisation between C. puella and C. pulchellum in the past. Three of the EPIC (exon-primed intron-crossing) molecular markers, AgT, MLC and PRMT, were sequenced from the above mentioned puella-group spp. and C. ornatum, from the territory considered and elsewhere, plus from C. persicum (Lohmann, 1993) from Iran. Coenagrion ponticum thus entered molecular analysis for the first time. Coenagrion australocaspicum appeared divergent from the rest of the group. Coenagrion ponticum appeared inside the cluster of C. pulchellum while C. persicum outside it, so the species status of the former is challenged while that of the latter confirmed. The joint tree showed two clades in C. puella (with no difference in outer characters), one composed of specimens originating from the regions to the south of the Main Caucasian Range and from the Caspian areas of Dagestan including foothills, and the other with specimens from Inner Dagestan (plus one locality in the Caspian-faced foothills), West Caucasus, Western Europe, Ural and West Siberia.

Genetic Variation and Population Structure of the Old World Bollworm Helicoverpa armigera (Hübner, 1808) (Lepidoptera: Noctuidae) in Ethiopia

Helicoverpa armigera is one of the most destructive insect pests of economically valuable crops in the world. Despite its economic importance, the population genetic structure of this insect remains unexplored in Ethiopia. To investigate the genetic diversity and population structure of H. armigera, we sampled 170 individuals from 15 populations throughout Ethiopia. We sequenced a fragment of the mitochondrial cytochrome b (cyt b) gene and five exon-primed intron-crossing (EPIC) markers. Twenty cyt b haplotypes with low-to-moderate haplotype diversity (mean Hd = 0.537) and high nucleotide diversity (mean Pi = 0.00339) were identified. The most frequently observed and widely distributed cyt b haplotype was designated as Hap_1 (67.058%), which is identical to sequences found across the globe. Tajima's D and Fu's F for the cyt b data were negative, supporting a model of population expansion. Within populations, a mean of 2.493 alleles/locus was recorded across the five EPIC loci, ranging from 1.200 to 3.600 alleles/locus. The highest mean effective number of alleles/population was 2.369 and the lowest was 1.178. The mean observed heterozygosity (HO) of the five loci (0-0.289; mean 0.104 ± 0.020) was lower than the expected heterozygosity (HE) (0.095-0.523; mean 0.258 ± 0.028). AMOVA detected significant genetic structure with 61% of the total molecular genetic variation of EPIC genotypes occurring between populations, suggesting a considerable degree of differentiation among populations. STRUCTURE analyses clustered the H. armigera populations into three distinct population groups but very low isolation by distance (R2 = 0.0132, P < 0.05).

Comparative phylogeography of Aedes mosquitoes and the role of past climatic change for evolution within Africa.

The study of demographic processes involved in species diversification and evolution ultimately provides explanations for the complex distribution of biodiversity on earth, indicates regions important for the maintenance and generation of biodiversity, and identifies biological units important for conservation or medical consequence. African and forest biota have both received relatively little attention with regard to understanding their diversification, although one possible mechanism is that this has been driven by historical climate change. To investigate this, we implemented a standard population genetics approach along with Approximate Bayesian Computation, using sequence data from two exon‐primed intron‐crossing (EPIC) nuclear loci and mitochondrial cytochrome oxidase subunit I, to investigate the evolutionary history of five medically important and inherently forest dependent mosquito species of the genus Aedes. By testing different demographic hypotheses, we show that Aedes bromeliae and Aedes lilii fit the same model of lineage diversification, admixture, expansion, and recent population structure previously inferred for Aedes aegypti. In addition, analyses of population structure show that Aedes africanus has undergone lineage diversification and expansion while Aedes hansfordi has been impacted by population expansion within Uganda. This congruence in evolutionary history is likely to relate to historical climate‐driven habitat change within Africa during the late Pleistocene and Holocene epoch. We find differences in the population structure of mosquitoes from Tanzania and Uganda compared to Benin and Uganda which could relate to differences in the historical connectivity of forests across the continent. Our findings emphasize the importance of recent climate change in the evolution of African forest biota.

Bottomless barrel-sponge species in the Indo-Pacific?

The use of nuclear markers, in addition to traditional mitochondrial markers, helps to clarify hidden patterns of genetic structure in natural populations (Palumbi & Baker, 1994). This is particularly evident among demosponges that possess slow mitochondrial evolutionary rates compared to Bilateria, where nuclear intron markers can aid in the understanding of shallow level phylogenetic relationships (Shearer et al., 2002). Ideally, these nuclear markers (i) are evolutionary well-conserved across different lineages, (ii) produce amplicons holding a number of sites with sufficient variability to answer the relevant phylogenetic question, (iii) derive from single copy genes (see review in Zhang & Hewitt, 2003). A popular method to amplify intron markers uses EPIC (Exon-Primed, Intron-Crossing) primers that anneal to the more conserved flanking exon regions and subsequently bridge the intron during amplification (Palumbi & Baker, 1994).

Universal PCR primers for ribosomal protein gene introns of fish

Human ribosomal protein (RP) gene sequences with respect to intron/exon structures and corresponding cDNA or genomic data of fish species were obtained from the GenBank database. Based on conserved exon sequences, 128 primer pairs for 41 genes were designed for exon-primed intron-crossing (EPIC) polymerase chain reaction (PCR). In reference to the draft genome sequences of the Pacific bluefin tuna (Thunnus orientalis), 12 primer pairs expected to amplify introns of the bluefin tuna with lengths of 500–1000 bp were selected and applied to six distantly related fish species belonging to the Orders Clupeiformes, Tetraodontiformes, Pleuronectiformes, Perciformes, Scorpaeniformes, and Anguilliformes. PCR amplification was observed for at least four species in each primer pair, and all fragments were larger than those expected for intronless amplification. Single fragment amplification was observed for at least seven primer pairs per species. Fragment sizes of the bluefin tuna for nine primer pairs corresponded to those expected from the genomic data. Thus, our primer pairs are potentially applicable to a wide variety of fish species and serve as an initial step for isolating single-copy nuclear DNA sequences.

Insect Resistance to Bacillus thuringiensis Toxin Cry2Ab Is Conferred by Mutations in an ABC Transporter Subfamily A Protein.

The use of conventional chemical insecticides and bacterial toxins to control lepidopteran pests of global agriculture has imposed significant selection pressure leading to the rapid evolution of insecticide resistance. Transgenic crops (e.g., cotton) expressing the Bt Cry toxins are now used world wide to control these pests, including the highly polyphagous and invasive cotton bollworm Helicoverpa armigera. Since 2004, the Cry2Ab toxin has become widely used for controlling H. armigera, often used in combination with Cry1Ac to delay resistance evolution. Isolation of H. armigera and H. punctigera individuals heterozygous for Cry2Ab resistance in 2002 and 2004, respectively, allowed aspects of Cry2Ab resistance (level, fitness costs, genetic dominance, complementation tests) to be characterised in both species. However, the gene identity and genetic changes conferring this resistance were unknown, as was the detailed Cry2Ab mode of action. No cross-resistance to Cry1Ac was observed in mutant lines. Biphasic linkage analysis of a Cry2Ab-resistant H. armigera family followed by exon-primed intron-crossing (EPIC) marker mapping and candidate gene sequencing identified three independent resistance-associated INDEL mutations in an ATP-Binding Cassette (ABC) transporter gene we named HaABCA2. A deletion mutation was also identified in the H. punctigera homolog from the resistant line. All mutations truncate the ABCA2 protein. Isolation of further Cry2Ab resistance alleles in the same gene from field H. armigera populations indicates unequal resistance allele frequencies and the potential for Bt resistance evolution. Identification of the gene involved in resistance as an ABC transporter of the A subfamily adds to the body of evidence on the crucial role this gene family plays in the mode of action of the Bt Cry toxins. The structural differences between the ABCA2, and that of the C subfamily required for Cry1Ac toxicity, indicate differences in the detailed mode-of-action of the two Bt Cry toxins.

Exon primed intron-crossing (EPIC) markers reveal natural hybridization and introgression in Actinidia (Actinidiaceae) with sympatric distribution

Interspecific hybridization is widespread among plants. Understanding the phylogenetic relationships among species is necessary for revealing the potential hybridization events. Actinidia, best known as kiwifruit genus found throughout a wide range in eastern Asian from Indonesia to Siberia. In this study, phylogenetic relationships of Actinidia species with sympatric distributions were investigated using three chloroplast introns (trnL-F, atpB-rbcL and rpl32-trnL) and three Exon primed intron-crossing (EPIC) markers. Chloroplast phylogeny supports non-monophyly of the five species studied excluding Actinidia fulvicoma var. lanata. The non-monophyly was also revealed by EPIC markers. Our results showed EPIC markers are more variable and informative for phylogenetic inference than that of chloroplast markers. The incongruences between loci from the plastid and nuclear DNA phylogenic trees may stem from incomplete lineage sorting or historical introgression hybridization. Incomplete lineage sorting may explain the non-monophyly between Actinidia chrysantha (section Maculatae) and other four species (section Stellatae), and introgression hybridization and high level of interspecific gene flow may explain the non-monophyly among the species of sect. Stellatae. Thus, natural hybridization and introgression may be common in Actinidia with sympatric distribution.

Primers for low-copy nuclear genes in Metrosideros and cross-amplification in Myrtaceae.

• Premise of the study: Primers were developed to amplify low-copy nuclear genes in Hawaiian Metrosideros (Myrtaceae).• Methods and Results: Data from a pooled 454 Titanium run of the partial transcriptomes of four Metrosideros taxa were used to identify the loci of interest. Ten exon-primed intron-crossing (EPIC) markers were amplified and sequenced directly with success in Metrosideros, as well as in a representative selection of Myrtaceae, including Syzygium, Psidium, and Melaleuca for most of the markers. The loci amplified ranged between 500 and 1100 bp, and up to 117 polymorphic sites were observed within an individual gene alignment. Two introns contained microsatellites in some of the species.• Conclusions: These novel primer pairs should be useful for phylogenetic analysis and population genetics of a broad range of Myrtaceae, particularly the diverse fleshy-fruited tribes Syzygieae and Myrteae.

Developing Exon-Primed Intron-Crossing (EPIC) markers for population genetic studies in three Aedes disease vectors.

Aedes aegypti, Aedes notoscriptus, and Aedes albopictus are important vectors of many arboviruses implicated in human disease such as dengue fever. Genetic markers applied across vector species can provide important information on population structure, gene flow, insecticide resistance, and taxonomy, however, robust microsatellite markers have proven difficult to develop in these species and mosquitoes generally. Here we consider the utility and transferability of 15 Ribosome protein (Rp) Exon-Primed Intron-Crossing (EPIC) markers for population genetic studies in these 3 Aedes species. Rp EPIC markers designed for Ae. aegypti also successfully amplified populations of the sister species, Ae. albopictus, as well as the distantly related species, Ae. notoscriptus. High SNP and good indel diversity in sequenced alleles plus support for amplification of the same regions across populations and species were additional benefits of these markers. These findings point to the general value of EPIC markers in mosquito population studies.

New EPIC nuclear DNA sequence markers to improve the resolution of phylogeographic studies of coenagrionids and other odonates

While phylogeographic data provide valuable information to inform conservation plans, there are comparatively few Odonata phylogeographic studies. This lack of research is partially due to a lack of independent DNA markers with appropriate levels of polymorphism that PCR-amplify in a range of species. We followed an exon-primed, intron-crossing (EPIC) PCR strategy to develop five new, polymorphic nuclear DNA sequence loci (six distinct DNA fragments) for the southern damselfly Coenagrion mercuriale. These markers were: cell division cycle 5 protein (CDC5), arginine methyltransferase (PRMT), acetylglucosaminyl-transferase (AgT), myosin light chain (MLC) and phosphoglucose isomerase (PGI). Between three and five of these new markers could be PCR-amplified in five other species from the genus Coenagrion; one locus (PRMT) can be used in 26 other species of odonates that we examined, including three species of Anisoptera belonging to the genus Onychogomphus. These new nuclear genetic markers will be useful for...

Development of novel, exon-primed intron-crossing (EPIC) markers from EST databases and evaluation of their phylogenetic utility in Commiphora (Burseraceae).

• Premise of the study: Novel nuclear exon-primed intron-crossing (EPIC) markers were developed to increase phylogenetic resolution among recently diverged lineages in the frankincense and myrrh family, Burseraceae, using Citrus, Arabidopsis, and Oryza genome resources.• Methods and Results: Primer pairs for 48 nuclear introns were developed using the genome resource IntrEST and were screened using species of Commiphora and other Burseraceae taxa. Four putative intron regions (RPT6A, BXL2, mtATP Synthase D, and Rab6) sequenced successfully for multiple taxa and recovered phylogenies consistent with those of existing studies. In some cases, these regions yielded informative sequence variation on par with that of the nuclear ribosomal DNA internal transcribed spacer.• Conclusions: The combination of freely available genome resources and our design criteria have uncovered four single-copy nuclear intron regions that are useful for phylogenetic reconstruction of Burseraceae taxa. Because our EPIC primers also amplify Arabidopsis, we recommend their trial in other rosid and eudicot lineages.

Molecular genetic analysis of cereal β-amylase genes using exon-primed intron-crossing (EPIC) PCR

The proteins encoded by cereal β-amylase genes Bamy1 and Bamy2 genes play an important role in seedling germination and in the brewing process. Here, we use exon-primed intron-crossing (EPIC) to analyse Bamy1 and Bamy2 genetic diversity among 38 accessions belonging to six Poaceae tribes. DNA sequence alignment of multiple Poaceae species β-amylase sequences allowed design of EPIC primers that simultaneously amplify Bamy1 and Bamy2 in all the cereal species investigated. The genetic variation observed in the samples investigated is analysed and discussed, and illustrates the effectiveness of this approach for intra- and interspecific analysis in plant species.

Exon-primed intron-crossing (EPIC) markers as a tool for ant phylogeography

Exon-primed intron-crossing (EPIC) markers as a tool for ant phylogeography. Due to their local abundance, diversity of adaptations and worldwide distribution, ants are a classic example of adaptive radiation. Despite this evolutionary and ecological importance, phylogeographical studies on ants have relied largely on mitochondrial markers. In this study we design and test exon-primed intron-crossing (EPIC) markers, which can be widely used to uncover ant intraspecific variation. Candidate markers were obtained through screening the available ant genomes for unlinked conserved exonic regions interspersed with introns. A subset of 15 markers was tested in vitro and showed successful amplification in several phylogenetically distant ant species. These markers represent an important step forward in ant phylogeography and population genetics, allowing for more extensive characterization of variation in ant nuclear DNA without the need to develop species-specific markers.

A Brave New World for an Old World Pest: Helicoverpa armigera (Lepidoptera: Noctuidae) in Brazil

The highly polyphagous Old World cotton bollworm Helicoverpa armigera is a quarantine agricultural pest for the American continents. Historically H. armigera is thought to have colonised the American continents around 1.5 to 2 million years ago, leading to the current H. zea populations on the American continents. The relatively recent species divergence history is evident in mating compatibility between H. zea and H. armigera under laboratory conditions. Despite periodic interceptions of H. armigera into North America, this pest species is not believed to have successfully established significant populations on either continent. In this study, we provide molecular evidence via mitochondrial DNA (mtDNA) cytochrome oxidase I (COI) and cytochrome b (Cyt b) partial gene sequences for the successful recent incursion of H. armigera into the New World, with individuals being detected at two sites (Primavera do Leste, Pedra Preta) within the State of Mato Grosso in Brazil. The mtDNA COI and Cyt b haplotypes detected in the Brazilian H. armigera individuals are common throughout the Old World, thus precluding identification of the founder populations. Combining the two partial mtDNA gene sequences showed that at least two matrilines are present in Brazil, while the inclusion of three nuclear DNA Exon-Primed Intron-Crossing (EPIC) markers identified a further two possible matrilines in our samples. The economic, biosecurity, resistance management, ecological and evolutionary implications of this incursion are discussed in relation to the current agricultural practices in the Americas.

PCR survey of 50 introns in animals: Cross-amplification of homologous EPIC loci in eight non-bilaterian, protostome and deuterostome phyla

Exon Primed Intron Crossing (EPIC) markers provide molecular tools that are susceptible to be variable within species while remaining amplifiable by PCR using potentially universal primers. In this study we tested the possibility of obtaining PCR products from 50 EPIC markers on 23 species belonging to seven different phyla (Porifera, Cnidaria, Arthropoda, Nematoda, Mollusca, Annelida, Echinodermata) using 70 new primer pairs. A previous study had identified and tested those loci in a dozen species, including another phylum, Urochordata (Chenuil et al., 2010). Results were contrasted among species. The best results were achieved with the oyster (Mollusca) where 28 loci provided amplicons susceptible to contain an intron according to their size. This was however not the case with the other mollusk Crepidula fornicata, which seems to have undergone a reduction in intron number or intron size. In the Porifera, 13 loci appeared susceptible to contain an intron, a surprisingly high number for this phylum considering its phylogenetic distance with genomic data used to design the primers. For two cnidarian species, numerous loci (24) were obtained. Ecdysozoan phyla (arthropods and nematodes) proved less successful than others as expected considering reports of their rapid rate of genome evolution and the worst results were obtained for several arthropods. Some general patterns among phyla arose, and we discuss how the results of this EPIC survey may give new insights into genome evolution of the study species. This work confirms that this set of EPIC loci provides an easy-to-use toolbox to identify genetic markers potentially useful for population genetics, phylogeography or phylogenetic studies for a large panel of metazoan species. We then argue that obtaining diploid sequence genotypes for these loci became simple and affordable owing to Next-Generation Sequencing development. Species surveyed in this study belong to several genera (Acanthaster, Alvinocaris, Aplysina, Aurelia, Crepidula, Eunicella, Hediste, Hemimysis, Litoditis, Lophelia, Mesopodopsis, Mya, Ophiocten, Ophioderma, Ostrea, Pelagia, Platynereis, Rhizostoma, Rimicaris), two of them, belonging to the family Vesicomydae and Eunicidae, could not be determined at the genus level.

Exon-primed intron-crossing (EPIC) markers for evolutionary studies of Ficus and other taxa in the fig family (Moraceae).

• Premise of the study: The genus Ficus (fig trees) comprises ca. 750 species of trees, vines, and stranglers found in tropical forests throughout the world. Fig trees are keystone species in many tropical forests, and their relationship with host-specific wasp pollinators has received much attention, although many questions remain unresolved regarding the levels of host specificity, cospeciation, and the role of hybridization in fig and wasp speciation. We developed exon-primed intron-crossing (EPIC) markers to obtain phylogenetic resolution needed to address these questions.• Methods and Results: Expressed sequence tags (ESTs) from F. elastica were compared to Arabidopsis and Populus genomes to locate introns and to design primers in flanking exons. Primer pairs for 80 EPIC markers were tested in samples from divergent clades within Ficus and the outgroup Poulsenia (Moraceae).• Conclusions: Thirty-one EPIC markers were successfully sequenced across Ficus, and 29 of the markers also amplified in Poulsenia, indicating broad transferability within Moraceae. All of the EPIC markers were polymorphic and showed levels of polymorphism similar to that of the widely used internal transcribed spacer (ITS).

Population Genetic Structure of the Cotton Bollworm Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) in India as Inferred from EPIC-PCR DNA Markers

Helicoverpa armigera is an important pest of cotton and other agricultural crops in the Old World. Its wide host range, high mobility and fecundity, and the ability to adapt and develop resistance against all common groups of insecticides used for its management have exacerbated its pest status. An understanding of the population genetic structure in H. armigera under Indian agricultural conditions will help ascertain gene flow patterns across different agricultural zones. This study inferred the population genetic structure of Indian H. armigera using five Exon-Primed Intron-Crossing (EPIC)-PCR markers. Nested alternative EPIC markers detected moderate null allele frequencies (4.3% to 9.4%) in loci used to infer population genetic structure but the apparently genome-wide heterozygote deficit suggests in-breeding or a Wahlund effect rather than a null allele effect. Population genetic analysis of the 26 populations suggested significant genetic differentiation within India but especially in cotton-feeding populations in the 2006–07 cropping season. In contrast, overall pair-wise F ST estimates from populations feeding on food crops indicated no significant population substructure irrespective of cropping seasons. A Baysian cluster analysis was used to assign the genetic make-up of individuals to likely membership of population clusters. Some evidence was found for four major clusters with individuals in two populations from cotton in one year (from two populations in northern India) showing especially high homogeneity. Taken as a whole, this study found evidence of population substructure at host crop, temporal and spatial levels in Indian H. armigera, without, however, a clear biological rationale for these structures being evident.

Population Genetic Structure of Aedes (Stegomyia) aegypti (L.) at a Micro-Spatial Scale in Thailand: Implications for a Dengue Suppression Strategy

BackgroundThe genetic population structure of Aedes (Stegomyia) aegypti (L.), the main vector of dengue virus, is being investigated in areas where a novel dengue suppression program is to be implemented. The aim of the program is to release and establish mosquito populations with impaired virus transmission capabilities. To model effects of the release and devise protocols for its implementation, information about the genetic structure of populations at a range of spatial scales is required.Methodology/Principal FindingsThis study investigates a potential release site in the Hua Sam Rong Subdistrict of Plaeng Yao District, Chachoengsao Province, in eastern Thailand which comprises a complex of five villages within a 10 km radius. Aedes aegypti resting indoors was sampled at four different times of year from houses within the five villages. Genetic markers were used to screen the mosquitoes: two Exon Primed Intron Crossing (EPIC) markers and five microsatellite markers. The raw allele size was determined using several statistical software packages to analyze the population structure of the mosquito. Estimates of effective population size for each village were low, but there was no evidence of genetic isolation by geographic distance.ConclusionsThe presence of temporary genetic structure is possibly caused by genetic drift due to large contributions of adults from a few breeding containers. This suggests that the introduction of mosquitoes into an area needs to proceed through multiple releases and targeting of sites where mosquitoes are emerging in large numbers.

Combining multiple autosomal introns for studying shallow phylogeny and taxonomy of Laurasiatherian mammals: Application to the tribe Bovini (Cetartiodactyla, Bovidae)

Mitochondrial sequences are widely used for species identification and for studying phylogenetic relationships among closely related species or populations of the same species. However, many studies of mammals have shown that the maternal history of the mitochondrial genome can be discordant with the true evolutionary history of the taxa. In such cases, the analyses of multiple nuclear genes can be more powerful for deciphering interspecific relationships.Here, we designed primers for amplifying 13 new exon-primed intron-crossing (EPIC) autosomal loci for studying shallow phylogeny and taxonomy of Laurasiatherian mammals. Three criteria were used for the selection of the markers: gene orthology, a PCR product length between 600 and 1200 nucleotides, and different chromosomal locations in the bovine genome. Positive PCRs were obtained from different species representing the orders Carnivora, Cetartiodactyla, Chiroptera, Perissodactyla and Pholidota.The newly developed markers were analyzed in a phylogenetic study of the tribe Bovini (the group containing domestic and wild cattle, bison, yak, African buffalo, Asian buffalo, and saola) based on 17 taxa and 18 nuclear genes, representing a total alignment of 13,095 nucleotides. The phylogenetic results were compared to those obtained from analyses of the complete mitochondrial genome and Y chromosomal genes. Our analyses support a basal divergence of the saola (Pseudoryx) and a sister-group relationship between yak and bison. These results contrast with recent molecular studies but are in better agreement with morphology. The comparison of pairwise nucleotide distances shows that our nuDNA dataset provides a good signal for identifying taxonomic levels, such as species, genera, subtribes, tribes and subfamilies, whereas the mtDNA genome fails because of mtDNA introgression and higher levels of homoplasy. Accordingly, we conclude that the genus Bison should be regarded as a synonym of Bos, with the European bison relegated to a subspecies rank within Bos bison. We compared our molecular dating estimates to the fossil record in order to propose a biogeographic scenario for the evolution of Bovini during the Neogene.

EvolMarkers: a database for mining exon and intron markers for evolution, ecology and conservation studies

Recent innovations in next-generation sequencing have lowered the cost of genome projects. Nevertheless, sequencing entire genomes for all representatives in a study remains expensive and unnecessary for most studies in ecology, evolution and conservation. It is still more cost-effective and efficient to target and sequence single-copy nuclear gene markers for such studies. Many tools have been developed for identifying nuclear markers, but most of these have focused on particular taxonomic groups. We have built a searchable database, EvolMarkers, for developing single-copy coding sequence (CDS) and exon-primed-intron-crossing (EPIC) markers that is designed to work across a broad range of phylogenetic divergences. The database is made up of single-copy CDS derived from BLAST searches of a variety of metazoan genomes. Users can search the database for different types of markers (CDS or EPIC) that are common to different sets of input species with different divergence characteristics. EvolMarkers can be applied to any taxonomic group for which genome data are available for two or more species. We included 82 genomes in the first version of EvolMarkers and have found the methods to be effective across Placozoa, Cnidaria, Arthropod, Nematoda, Annelida, Mollusca, Echinodermata, Hemichordata, Chordata and plants. We demonstrate the effectiveness of searching for CDS markers within annelids and show how to find potentially useful intronic markers within the lizard Anolis.