Exon-primed intron-crossing (EPIC) PCR markers of Helicoverpa armigera (Lepidoptera: Noctuidae)

The effect of missing data on evolutionary analysis of sequence capture bycatch, with application to an agricultural pest

Sequence capture is a genomic technique that selectively enriches target sequences before high throughput next-generation sequencing, to generate specific sequences of interest. Off-target or 'bycatch' data are often discarded from capture experiments, but can be leveraged to address evolutionary questions under some circumstances. Here, we investigated the effects of missing data on a variety of evolutionary analyses using bycatch from an exon capture experiment on the global pest moth, Helicoverpa armigera. We added > 200 new samples from across Australia in the form of mitogenomes obtained as bycatch from targeted sequence capture, and combined these into an additional larger dataset to total > 1000 mitochondrial cytochrome c oxidase subunit I (COI) sequences across the species' global distribution. Using discriminant analysis of principal components and Bayesian coalescent analyses, we showed that mitogenomes assembled from bycatch with up to 75% missing data were able to return evolutionary inferences consistent with higher coverage datasets and the broader literature surrounding H. armigera. For example, low-coverage sequences broadly supported the delineation of two H. armigera subspecies and also provided new insights into the potential for geographic turnover among these subspecies. However, we also identified key effects of dataset coverage and composition on our results. Thus, low-coverage bycatch data can offer valuable information for population genetic and phylodynamic analyses, but caution is required to ensure the reduced information does not introduce confounding factors, such as sampling biases, that drive inference. We encourage more researchers to consider maximizing the potential of the targeted sequence approach by examining evolutionary questions with their off-target bycatch where possible-especially in cases where no previous mitochondrial data exists-but recommend stratifying data at different genome coverage thresholds to separate sampling effects from genuine genomic signals, and to understand their implications for evolutionary research.

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).

Invasion origin, rapid population expansion, and the lack of genetic structure of cotton bollworm (Helicoverpa armigera) in the Americas

In 2013, Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) was officially declared as present in Brazil and, after two years, the species was detected in the Caribbean and North America. Information on genetic features and accurate distribution of pests is the basis for agricultural protection policies. Furthermore, such knowledge is imperative to develop control strategies, understand the geographical range, and genetic patterns of this species in the Americas. Here, we carried out the widest sampling of H. armigera in the South American continent and Puerto Rico, after we estimated the diversity, demographic parameters, and genetic structure. The Internal Transcribed Spacer 1 (ITS1) nuclear marker was used to investigate the presence of putative hybrids between H. armigera and H. zea, and they were observed at a frequency of 1.5%. An ABC analysis, based in COI gene fragment, suggested Europe as the origin of South America specimens of H. armigeraand following a movement northward through the Caribbean. Three mtDNA genes and three nDNA markers revealed high genetic diversity distributed without the defined population structure of H. armigera in South America. Most of the genetic variation is within populations with a multidirectional expansion of H. armigera among morphoclimatic regions. High genetic diversity, rapid population expansion, and hybridization have implications for pest management since they suggest that adaptive alleles are spread through wide areas in South America that favor rapid local adaptation of H. armigera to new and disturbed environments (e.g., in agricultural areas).

Mitochondrial DNA genomes of five major Helicoverpa pest species from the Old and New Worlds (Lepidoptera: Noctuidae)

Five species of noctuid moths, Helicoverpa armigera, H. punctigera, H. assulta, H. zea, and H. gelotopoeon, are major agricultural pests inhabiting various and often overlapping global distributions. Visual identification of these species requires a great deal of expertise and misidentification can have repercussions for pest management and agricultural biosecurity. Here, we report on the complete mitochondrial genomes of H. assulta assulta and H. assulta afra, H. gelotopoeon, H. punctigera, H. zea, and H. armigera armigera and H. armigera conferta' assembled from high-throughput sequencing data. This study significantly increases the mitogenome resources for these five agricultural pests with sequences assembled from across different continents, including an H. armigera individual collected from an invasive population in Brazil. We infer the phylogenetic relationships of these five Helicoverpa species based on the 13 mitochondrial DNA protein-coding genes (PCG's) and show that two publicly available mitogenomes of H. assulta (KP015198 and KR149448) have been misidentified or incorrectly assembled. We further consolidate existing PCR-RFLP methods to cover all five Helicoverpa pest species, providing an updated method that will contribute to species differentiation and to future monitoring efforts of Helicoverpa pest species across different continents. We discuss the value of Helicoverpa mitogenomes to assist with species identification in view of the context of the rapid spread of H. armigera in the New World. With this work, we provide the molecular resources necessary for future studies of the evolutionary history and ecology of these species.

Movement Ecology of Pest Helicoverpa: Implications for Ongoing Spread

The recent introduction and spread of Helicoverpa armigera throughout South America highlight the invasiveness and adaptability of moths in the Helicoverpa genus. Long-range movement in three key members, H. armigera, H. zea, and H. punctigera, occurs by migration and international trade. These movements facilitate high population admixture and genetic diversity, with important economic, biosecurity, and control implications in today's agricultural landscape. This is particularly true for the spread of resistance alleles to transgenic crops expressing Bacillus thuringiensis (Bt) toxins that are planted over vast areas to suppress Helicoverpa spp. The ability to track long-distance movement through radar technology, population genetic markers, and/or long-distance dispersal modeling has advanced in recent years, yet we still know relatively little about the population trajectories or migratory routes in Helicoverpa spp. Here, we consider how experimental and theoretical approaches can be integrated to fill key knowledge gaps and assist management practices.

Detection of sister-species in invasive populations of the fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae) from Uganda

The fall armyworm (FAW) Spodoptera frugiperda (J. E. Smith) is a species native to the Americas. This polyphagous lepidopteran pest was first reported in Nigeria and the Democratic Republic of São Tomé and Principe in 2016, but its presence in eastern Africa has not been confirmed via molecular characterisation. In this study, FAW specimens from western and central Uganda were identified based on the partial mtDNA COI gene sequences, with mtDNA COI haplotypes matching those identified in Nigeria and São Tomé. In this study, we sequence an additional partial mtDNA Cyt b gene and also the partial mtDNA COIII gene in Ugandan FAW samples. We detected identical mitochondrial DNA haplotypes for both the mtDNA Cyt b and COI partial genes, while combining the mtDNA COI/Cyt b haplotypes and mtDNA COIII haplotypes enabled a new maternal lineage in the Ugandan corn-preferred FAW samples to be identified. Our results suggested that the African incursions of S. frugiperda involved at least three maternal lineages. Recent full genome, phylogenetic and microsatellite analyses provided evidence to support S. frugiperda as likely consisted of two sympatric sister species known as the corn-preferred and rice-preferred strains. In our Ugandan FAW populations, we identified the presence of mtDNA haplotypes representative of both sister species. It is not known if both FAW sister species were originally introduced together or separately, and whether they have since spread as a single population. Further analyses of additional specimens originally collected from São Tomé, Nigeria and throughout Africa would be required to clarify this issue. Importantly, our finding showed that the genetic diversity of the African corn-preferred FAW species is higher than previously reported. This potentially contributed to the success of FAW establishment in Africa. Furthermore, with the additional maternal lineages detected, there is likely an increase in paternal lineages, thereby increasing the diversity of the African FAW population. Knowledge of the FAW genetic diversity will be needed to assess the risks of introducing Bt-resistance traits and to understand the FAW incursion pathways into the Old World and its potential onward spread. The agricultural implications of the presence of two evolutionary divergent FAW lineages (the corn and the rice lineage) in the African continent are further considered and discussed.

Population Genetic Structure and Potential Incursion Pathways of the Bluetongue Virus Vector Culicoides brevitarsis (Diptera: Ceratopogonidae) in Australia

Culicoides brevitarsis is a vector of the bluetongue virus (BTV), which infects sheep and cattle. It is an invasive species in Australia with an assumed Asian/South East Asian origin. Using one mitochondrial marker (i.e., part of the cytochrome oxidase subunit I gene) and six nuclear markers, we inferred population genetic structure and possible incursion pathways for Australian C. brevitarsis. Nine mitochondrial haplotypes, with low nucleotide sequence diversity (0.0-0.7%) among these, were identified in a sample of 70 individuals from seven sites. Both sets of markers revealed a homogeneous population structure, albeit with evidence of isolation by distance and two genetically distinct clusters distributed along a north-to-south cline. No evidence of a cryptic species complex was found. The geographical distribution of the mitochondrial haplotypes is consistent with at least two incursion pathways into Australia since the arrival of suitable livestock hosts. By contrast, 15 mitochondrial haplotypes, with up to four times greater nucleotide sequence diversity (0.0-2.9%) among these, were identified in a sample of 16 individuals of the endemic C. marksi (sampled from a site in South Australia and another in New South Wales). A phylogenetic tree inferred using the mitochondrial marker revealed that the Australian and Japanese samples of C. brevitarsis are as evolutionarily different from one another as some of the other Australian species (e.g., C. marksi, C. henryi, C. pallidothorax) are. The phylogenetic tree placed four of the species endemic to Australia (C. pallidothorax, C. bundyensis, C. marksi, C. henryi) in a clade, with a fifth such species (C. bunrooensis) sharing a common ancestor with that clade and a clade comprising two Japanese species (C. verbosus, C. kibunensis).

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.

Transgenic crops expressing the insecticidal protein Cry2Ab from Bacillus thuringiensis (Bt) are used worldwide to suppress damage by lepidopteran pests, often used in combination with Cry1Ac toxin to delay resistance evolution. Until now, the Cry2Ab mode of action and the mechanism of resistance were unknown, with field-isolated Cry2Ab resistant Helicoverpa armigera showing no cross-resistance to Cry1Ac. In this study, biphasic linkage analysis of a Cry2Ab-resistant H. armigera family followed by EPIC marker mapping and candidate gene sequencing identified three independent INDEL mutations in an ATP-Binding Cassette transporter subfamily A gene (ABCA2). A deletion mutation was identified in the same gene of resistant H. punctigera. All four mutations are predicted to truncate the ABCA2 protein. This is the first molecular genetic characterization of insect resistance to the Cry2Ab toxin, and detection of diverse Cry2Ab resistance alleles will contribute to understanding the micro-evolutionary processes that underpinned lepidopteran Bt-resistance.

High nucleotide diversity and limited linkage disequilibrium in Helicoverpa armigera facilitates the detection of a selective sweep

Insecticides impose extreme selective pressures on populations of target pests and so insecticide resistance loci of these species may provide the footprints of 'selective sweeps'. To lay the foundation for future genome-wide scans for selective sweeps and inform genome-wide association study designs, we set out to characterize some of the baseline population genomic parameters of one of the most damaging insect pests in agriculture worldwide, Helicoverpa armigera. To this end, we surveyed nine Z-linked loci in three Australian H. armigera populations. We find that estimates of π are in the higher range among other insects and linkage disequilibrium decays over short distances. One of the surveyed loci, a cytochrome P450, shows an unusual haplotype configuration with a divergent allele at high frequency that led us to investigate the possibility of an adaptive introgression around this locus.

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.

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.

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.

Will the real Bemisia tabaci please stand up?

Since Panayiotis Gennadius first identified the whitefly, Aleyrodes tabaci in 1889, there have been numerous revisions of the taxonomy of what has since become one of the world's most damaging insect pests. Most of the taxonomic revisions have been based on synonymising different species under the name Bemisia tabaci. It is now considered that there is sufficient biological, behavioural and molecular genetic data to support its being a cryptic species complex composed of at least 34 morphologically indistinguishable species. The first step in revising the taxonomy of this complex involves matching the A. tabaci collected in 1889 to one of the members of the species complex using molecular genetic data. To do this we extracted and then amplified a 496 bp fragment from the 3' end of the mitochondrial DNA cytochrome oxidase one (mtCOI) gene belonging to a single whitefly taken from Gennadius' original 1889 collection. The sequence identity of this 123 year-old specimen enabled unambiguous assignment to a single haplotype known from 13 Mediterranean locations across Greece and Tunisia. This enabled us to unambiguously assign the Gennadius A. tabaci to the member of the B. tabaci cryptic species complex known as Mediterranean or as it is commonly, but erroneously referred to, as the 'Q-biotype'. Mediterranean is therefore the real B. tabaci. This study demonstrates the importance of matching museum syntypes with known species to assist in the delimitation of cryptic species based on the organism's biology and molecular genetic data. This study is the first step towards the reclassification of B. tabaci which is central to an improved understanding how best to manage this globally important agricultural and horticultural insect pest complex.

A novel primer set for multilocus phylogenetic inference in East African cichlid fishes

The cichlid fishes in the East African Great Lakes are a prime model system for the study of adaptive radiation. Therefore, the availability of an elaborate phylogenetic framework is an important prerequisite. Previous phylogenetic hypotheses on East African cichlids are mainly based on mitochondrial and/or fragment-based markers, and, to date, no taxon-rich phylogeny exists that is based on multilocus DNA sequence data. Here, we present the design of an extensive new primer set (24 nuclear makers) for East African cichlids that will be used for multilocus phylogenetic analyses in the future. The primers are designed to work for both Sanger sequencing and next-generation sequencing with the 454 technology. As a proof of principle, we validate these primers in a phylogenetically representative set of 16 cichlid species from Lake Tanganyika and main river systems in the area and provide a basic evaluation of the markers with respect to marker length and diversity indices.

Development of Exon-Primed Intron-Crossing (EPIC) PCR primers for the malaria vector Anopheles pseudopunctipennis (Diptera: Culicidae)

Using the Anopheles gambiae Giles genome as a template, we designed, screened and identified 14 novel Exon-Primed Intron-Crossing (EPIC) PCR primer pairs for Anopheles pseudopunctipennis Theobald 1901, a major vector of human Plasmodium sp. in South America. These primers were designed to target the conserved regions flanking consecutive exons of different genes and enabled the amplification of 17 loci of which nine were polymorphic. Polymorphisms at these loci ranged from two to four alleles. Intron length polymorphism analysis is a useful tool, which will allow the study of the population structure of this mosquito species, which remains poorly understood.

Quantitative real-time PCR with SYBR Green detection to assess gene duplication in insects: study of gene dosage in Drosophila melanogaster (Diptera) and in Ostrinia nubilalis (Lepidoptera)

BACKGROUND

The accurate determination of the number of copies of a gene in the genome (gene dosage) is essential for a number of genetic analyses. Quantitative real time PCR (qPCR) with TaqMan detection has shown advantages over traditional Southern-blot and FISH techniques, however the high costs of the required labeled probes is an important limitation of this method. qPCR with SYBR Green I detection is a simple and inexpensive alternative, but it has never been applied to the determination of the copy number of low copy number genes in organisms with high allelic variability (as some insects), where a very small margin of error is essential.

FINDINGS

We have tested the suitability of the qPCR with SYBR Green I detection methodology for the detection of low copy number genes in two insects: the genetically well characterized Drosophila melanogaster (Diptera) and the poor genetically characterized Ostrinia nubilalis (Lepidoptera). The system was applied to determine the copy number of: (1) the O. nubilalis cadherin gene, involved in the mode of action of Bacillus thuringiensis toxins, which showed indirect evidence of duplication, and (2) the D. melanogaster BarH1 and BarH2 genes, located within the Bar region of the X chromosome, to clearly determine whether they both are covered by the tandem duplication in the classical Bar (B1) mutant. Our results showed that the O. nubilalis cadherin gene is an autosomal single copy gene and that BarH1, but not BarH2, is duplicated in the Drosophila B1 mutant.

CONCLUSIONS

This work shows that qPCR with SYBR Green I detection can be specific and accurate enough to distinguish between one and two gene copies per haploid genome of genes with high allelic variability. The technique is sensitive enough to give reliable results with a minimum amount of sample (DNA from individual thoraxes) and to detect gene duplications in tandem.

Genome-wide identification of intron fragment insertion mutations and their potential use as SCAR molecular markers in the soybean

Introns often have a high probability of mutation as a result of DNA insertions and deletions (indels). In this study, 503 introns with exon-derived insertions were identified using a comprehensive search of the soybean genome. Of the 375 pairs of PCR primer sets designed for the loci in question, 161 primer sets amplified length polymorphism among nine soybean varieties and were identified as soybean gene-intron-driven functional sequence characterized amplified region (SCAR) markers. These SCAR markers are distributed among all 20 of the soybean chromosomes, and they developed from numerous genes involved in various physiological and biochemical processes that influence important agronomic traits of the soybean. The development of these novel gene-driven functional SCAR markers was fast and cost effective, and their use will facilitate molecular-assisted breeding of the soybean.

Generation of microsatellite repeat families by RTE retrotransposons in lepidopteran genomes

Background

Developing lepidopteran microsatellite DNA markers can be problematical, as markers often exhibit multiple banding patterns and high frequencies of non-amplifying "null" alleles. Previous studies identified sequences flanking simple sequence repeat (SSR) units that are shared among many lepidopteran species and can be grouped into microsatellite-associated DNA families. These families are thought to be associated with unequal crossing-over during DNA recombination or with transposable elements (TEs).

Results

We identified full-length lepidopteran non-LTR retrotransposable elements of the RTE clade in Heliconius melpomene and Bombyx mori. These retroelements possess a single open reading frame encoding the Exonuclease/Endonuclease/Phosphatase and the Reverse Transcriptase/nLTR domains, a 5' UTR (untranslated region), and an extremely short 3' UTR that regularly consists of SSR units. Phylogenetic analysis supported previous suggestions of horizontal transfer among unrelated groups of organisms, but the diversity of lepidopteran RTE elements appears due to ancient divergence of ancestral elements rather than introgression by horizontal transfer. Similarity searches of lepidopteran genomic sequences in GenBank identified partial RTE elements, usually consisting of the 3' terminal region, in 29 species. Furthermore, we identified the C-terminal end of the Reverse Transcriptase/nLTR domain and the associated 3' UTR in over 190 microsatellite markers from 22 lepidopteran species, accounting for 10% of the lepidopteran microsatellites in GenBank. Occasional retrotransposition of autonomous elements, frequent retrotransposition of 3' partial elements, and DNA replication slippage during retrotransposition offers a mechanistic explanation for the association of SSRs with RTE elements in lepidopteran genomes.

Conclusions

Non-LTR retrotransposable elements of the RTE clade therefore join a diverse group of TEs as progenitors of SSR units in various organisms. When microsatellites are isolated using standard SSR enrichment protocols and primers designed at complementary repeated regions, amplification from multiple genomic sites can cause scoring difficulties that compromise their utility as markers. Screening against RTE elements in the isolation procedure provides one strategy for minimizing this problem.