Morphometric and physical characteristics distinguishing adult Patagonian lamprey, Geotria macrostoma from the pouched lamprey, Geotria australis.

The pouched lamprey, Geotria australis Gray, 1851, has long been considered monotypic in the Geotriidae family with a wide southern temperate distribution across Australasia and South America. Recent studies have provided molecular and morphological evidence for a second Geotria species in South America; Geotria macrostoma (Burmeister, 1868). The aim of this study was to determine morphometric and physical characteristics of adult G. macrostoma that further differentiate this re-instated species of Geotriidae from G. australis. The diagnostic features discriminating immature adult G. macrostoma from G. australis when entering fresh water, are distinct differences in dentition, oral papillae and fimbriae counts and differences in coloration. In addition, G. macrostoma display greater growth of the prebranchial region and oral disc and has a deeper body depth and higher condition factor. All current ecological knowledge of the genus Geotria is based on Australasian populations, which may not be applicable to G. macrostoma. To ensure the conservation and protection of the Patagonian lamprey as a re-identified species, further investigations are needed to understand its life history, biology and ecology throughout its range.

Non-ammocoete larvae of Palaeozoic stem lampreys.

Ammocoetes-the filter-feeding larvae of modern lampreys-have long influenced hypotheses of vertebrate ancestry1-7. The life history of modern lampreys, which develop from a superficially amphioxus-like ammocoete to a specialized predatory adult, appears to recapitulate widely accepted scenarios of vertebrate origin. However, no direct evidence has validated the evolutionary antiquity of ammocoetes, and their status as models of primitive vertebrate anatomy is uncertain. Here we report larval and juvenile forms of four stem lampreys from the Palaeozoic era (Hardistiella, Mayomyzon, Pipiscius, and Priscomyzon), including a hatchling-to-adult growth series of the genus Priscomyzon from Late Devonian Gondwana. Larvae of all four genera lack the defining traits of ammocoetes. They instead display features that are otherwise unique to adult modern lampreys, including prominent eyes, a cusped feeding apparatus, and posteriorly united branchial baskets. Notably, phylogenetic analyses find that these non-ammocoete larvae occur in at least three independent lineages of stem lamprey. This distribution strongly implies that ammocoetes are specializations of modern-lamprey life history rather than relics of vertebrate ancestry. These phylogenetic insights also suggest that the last common ancestor of hagfishes and lampreys was a macrophagous predator that did not have a filter-feeding larval phase. Thus, the armoured 'ostracoderms' that populate the cyclostome and gnathostome stems might serve as better proxies than living cyclostomes for the last common ancestor of all living vertebrates.

Revalidation of the Argentinian pouched lamprey Geotria macrostoma (Burmeister, 1868) with molecular and morphological evidence.

BACKGROUND: The Argentinian pouched lamprey, classified as Petromyzon macrostomus Burmeister, 1868 was first described in 1867 in De La Plata River, in Buenos Aires, Argentina, and subsequently recorded in several rivers from Patagonia. Since its original description, the validity of P. macrostomus was questioned by several ichthyologists and 36 years after its original discovery it was considered a junior synonym of Geotria australis Gray, 1851. For a long time, the taxonomic status of G. australis has been uncertain, largely due to the misinterpretations of the morphological alterations that occur during sexual maturation, including the arrangement of teeth, size and position of fins and cloaca, and the development of an exceptionally large gular pouch in males. In this study, the taxonomic status of Geotria from across the "species" range was evaluated using both molecular analysis and examination of morphological characteristics. METHODOLOGY/PRINCIPAL FINDINGS: Phylogenetic and species delimitation analyses based on mitochondrial DNA sequences of Cytochrome b (Cyt b) and Cytochrome C Oxidase Subunit 1 (COI) genes, along with morphological analysis of diagnostic characters reported in the original descriptions of the species were used to assess genetic and morphological variation within Geotria and to determine the specific status of the Argentinian lamprey. These analyses revealed that Geotria from Argentina constitutes a well differentiated lineage from Chilean and Australasian populations. The position of the cloaca and the distance between the second dorsal and caudal fins in sub-adult individuals, and at previous life stages, can be used to distinguish between the two species. In addition, the genetic distance between G. macrostoma and G. australis for the COI and Cyt b mitochondrial genes is higher than both intra- and inter-specific distances reported for other Petromyzontiformes. CONCLUSIONS/SIGNIFICANCE: Our results indicate that the Argentinian pouched lamprey, found along a broad latitudinal gradient on the south-west Atlantic coast of South America, should be named as Geotria macrostoma (Burmeister, 1868) and not as G. australis Gray 1851, returning to its earliest valid designation in Argentina. Geotria macrostoma can now be considered as the single lamprey species inhabiting Argentinian Patagonia, with distinct local adaptations and evolutionary potential. It is essential that this distinctiveness is recognized in order to guide future conservation and management actions against imminent threats posed by human actions in the major basins of Patagonia.

Phylogeography of the European brook lamprey (Lampetra planeri) and the European river lamprey (Lampetra fluviatilis) species pair based on mitochondrial data.

The European river lamprey Lampetra fluviatilis and the European brook lamprey Lampetra planeri (Block 1784) are classified as a paired species, characterized by notably different life histories but morphological similarities. Previous work has further shown limited genetic differentiation between these two species at the mitochondrial DNA level. Here, we expand on this previous work, which focused on lamprey species from the Iberian Peninsula in the south and mainland Europe in the north, by sequencing three mitochondrial marker regions of Lampetra individuals from five river systems in Ireland and five in southern Italy. Our results corroborate the previously identified pattern of genetic diversity for the species pair. We also show significant genetic differentiation between Irish and mainland European lamprey populations, suggesting another ichthyogeographic district distinct from those previously defined. Finally, our results stress the importance of southern Italian L. planeri populations, which maintain several private alleles and notable genetic diversity.

Identification and characterization of the lamprey cathepsin genes.

Cathepsins are key mammalian proteases that play an important role in the immune response. Several studies have revealed the versatile and critical functions of cathepsins. Here, we obtained ten kinds of cathepsin homologs and identified seven homologs with complete coding sequences. Phylogenetic analysis verified their identities and supported the classification of cathepsins into seven families, which is similar to other vertebrates. Tissue-specific expression analysis showed that all lamprey cathepsins (L-cathepsins) are present in the supraneural body (SB), kidney, gill, intestine, brain, heart, and liver, but their relative abundance varied among tissues. Additionally, we focused on the lamprey cathepsin L (L-cathepsin L) and used recombinant L-cathepsin L protein (rL-cathepsin L) to prepare anti rL-cathepsin L polyclonal antibodies, which were used to detect its distribution in lamprey tissues. The L-cathepsin L protein was primarily detected in the SB, kidney, gill, intestine, brain, and liver via western blot and immunohistochemistry assays. Importantly, quantitative real-time PCR (RT-PCR) revealed that the expression level of L-cathepsins mRNA significantly increased after exposure to three different stimuli (poly I:C, Staphylococcus aureus (S.a) and Vibro anguilarum (V.an)). This suggested that L-cathepsins may participate in defense processes. These results revealed that L-cathepsins may play key roles in the immune response to exogenous stimuli. The findings provide important information for future studies aiming to understand the molecular mechanisms underlying the immune response to pathogen invasion in lamprey.

Parasitism offers large rewards but carries high risks: Predicting parasitic strategies under different life history conditions in lampreys.

The loss of parasitism in metazoan lineages is often seen as unlikely, but it has occurred in some lineages (e.g., leeches, lampreys). How and why parasitism is lost is aptly addressed by studying lampreys, because extant species include a range of feeding modes and parasitism has been lost repeatedly. An individual-based model was developed to determine whether variations in survival and growth rates in the larval and juvenile stages could favour parasitic or nonparasitic strategies. A realization of the model for a Lampetra spp. population, a genus which includes parasitic and nonparasitic animals, indicated that both strategies could be successful. A different model realization of the nonparasitic species Lethenteron appendix also agreed with expectations, and only nonparasitic strategies were successful. Modelling anadromous Petromyzon marinus produced only parasitic animals, as expected, but suggested two different adult sizes should appear in the population, which has not been reported in the literature. Finally, a realization of an Ichthyomyzon castaneus population, known to be parasitic only, rarely selected for parasitism (c. 7% of model iterations), possibly because the population used to parameterize the model was unusual for the species. The results suggest that nonparasitic lineages in lampreys are common because parasitism, while offering better growth, also has lower survival. Additionally, nonparasitic species may be generated at different rates because growth and survival thresholds in the model favouring parasitism are close to observed estimates in some populations. Loss of parasitism can occur when life stages have different trade-offs in growth and survivability.

A lamprey view on the origins of neuroendocrine regulation of the thyroid axis.

This mini review summarizes the current knowledge of the hypothalamic-pituitary-thyroid (HPT) endocrine system in lampreys, jawless vertebrates. Lampreys and hagfish are the only two extant members of the class of agnathans, the oldest lineage of vertebrates. The high conservation of the hypothalamic-pituitary-gonadal (HPG) axis in lampreys makes the lamprey model highly appropriate for comparative and evolutionary analyses. However, there are still many unknown questions concerning the hypothalamic-pituitary (HP) axis in its regulation of thyroid activities in lampreys. As an example, the hypothalamic and pituitary hormone(s) that regulate the HPT axis have not been confirmed and/or characterized. Similar to gnathostomes (jawed vertebrates), lampreys produce thyroxine (T4) and triiodothyronine (T3) from thyroid follicles that are suggested to be involved in larval development, metamorphosis, and reproduction. The existing data provide evidence of a primitive, overlapping yet functional HPG and HPT endocrine system in lamprey. We hypothesize that lampreys are in an evolutionary intermediate stage of hypothalamic-pituitary development, leading to the emergence of the highly specialized HPG and HPT endocrine axes in jawed vertebrates. Study of the ancient lineage of jawless vertebrates, the agnathans, is key to understanding the origins of the neuroendocrine system in vertebrates.

A Noninvasive Tool to Assess the Distribution of Pacific Lamprey (Entosphenus tridentatus) in the Columbia River Basin.

The Pacific lamprey (Entosphenus tridentatus) is an anadromous fish once abundant throughout coastal basins of western North America that has suffered dramatic declines in the last century due primarily to human activities. Here, we describe the development of an environmental DNA (eDNA) assay to detect Pacific lamprey in the Columbia River basin. The eDNA assay successfully amplified tissue derived DNA of Pacific lamprey collected from 12 locations throughout the Columbia River basin. The assay amplifies DNA from other Entosphenus species found outside of the Columbia River basin, but is species-specific within this basin. As a result, the assay presented here may be useful for detecting Entosphenus spp. in geographic range beyond the Columbia River Basin. The assay did not amplify tissue or synthetically derived DNA of 14 commonly sympatric non-target species, including lampreys of the genus Lampetra, which are morphologically similar to Pacific lamprey in the freshwater larval stage.

Molecular cytogenetic characterization of chromosome site-specific repetitive sequences in the Arctic lamprey (Lethenteron camtschaticum, Petromyzontidae).

All extant lamprey karyotypes are characterized by almost all dot-shaped microchromosomes. To understand the molecular basis of chromosome structure in lampreys, we performed chromosome C-banding and silver staining and chromosome mapping of the 18S-28S and 5S ribosomal RNA (rRNA) genes and telomeric TTAGGG repeats in the Arctic lamprey (Lethenteron camtschaticum). In addition, we cloned chromosome site-specific repetitive DNA sequences and characterized them by nucleotide sequencing, chromosome in situ hybridization, and filter hybridization. Three types of repetitive sequences were detected; a 200-bp AT-rich repetitive sequence, LCA-EcoRIa that co-localized with the 18S-28S rRNA gene clusters of 3 chromosomal pairs; a 364-bp AT-rich LCA-EcoRIb sequence that showed homology to the EcoRI sequence family from the sea lamprey (Petromyzon marinus), which contains short repeats as centromeric motifs; and a GC-rich 702-bp LCA-ApaI sequence that was distributed on nearly all chromosomes and showed significant homology with the integrase-coding region of a Ty3/Gypsy family long terminal repeat (LTR) retrotransposon. All three repetitive sequences are highly conserved within the Petromyzontidae or within Petromyzontidae and Mordaciidae. Molecular cytogenetic characterization of these site-specific repeats showed that they may be correlated with programed genome rearrangement (LCA-EcoRIa), centromere structure and function (LCA-EcoRIb), and site-specific amplification of LTR retroelements through homogenization between non-homologous chromosomes (LCA-ApaI).

Inferring the demographic history underlying parallel genomic divergence among pairs of parasitic and nonparasitic lamprey ecotypes.

Understanding the evolutionary mechanisms generating parallel genomic divergence patterns among replicate ecotype pairs remains an important challenge in speciation research. We investigated the genomic divergence between the anadromous parasitic river lamprey (Lampetra fluviatilis) and the freshwater-resident nonparasitic brook lamprey (Lampetra planeri) in nine population pairs displaying variable levels of geographic connectivity. We genotyped 338 individuals with RAD sequencing and inferred the demographic divergence history of each population pair using a diffusion approximation method. Divergence patterns in geographically connected population pairs were better explained by introgression after secondary contact, whereas disconnected population pairs have retained a signal of ancient migration. In all ecotype pairs, models accounting for differential introgression among loci outperformed homogeneous migration models. Generating neutral predictions from the inferred divergence scenarios to detect highly differentiated markers identified greater proportions of outliers in disconnected population pairs than in connected pairs. However, increased similarity in the most divergent genomic regions was found among connected ecotype pairs, indicating that gene flow was instrumental in generating parallelism at the molecular level. These results suggest that heterogeneous genomic differentiation and parallelism among replicate ecotype pairs have partly emerged through restricted introgression in genomic islands.

Regulatory evolution of Tbx5 and the origin of paired appendages.

The diversification of paired appendages has been a major factor in the evolutionary radiation of vertebrates. Despite its importance, an understanding of the origin of paired appendages has remained elusive. To address this problem, we focused on T-box transcription factor 5 (Tbx5), a gene indispensable for pectoral appendage initiation and development. Comparison of gene expression in jawless and jawed vertebrates reveals that the Tbx5 expression in jawed vertebrates is derived in having an expression domain that extends caudal to the heart and gills. Chromatin profiling, phylogenetic footprinting, and functional assays enabled the identification of a Tbx5 fin enhancer associated with this apomorphic pattern of expression. Comparative functional analysis of reporter constructs reveals that this enhancer activity is evolutionarily conserved among jawed vertebrates and is able to rescue the finless phenotype of tbx5a mutant zebrafish. Taking paleontological evidence of early vertebrates into account, our results suggest that the gain of apomorphic patterns of Tbx5 expression and regulation likely contributed to the morphological transition from a finless to finned condition at the base of the vertebrate lineage.

Genesis of the vertebrate FoxP subfamily member genes occurred during two ancestral whole genome duplication events.

The vertebrate FoxP subfamily genes play important roles in the construction of essential functional modules involved in physiological and developmental processes. To explore the adaptive evolution of functional modules associated with the FoxP subfamily member genes, it is necessary to study the gene duplication process. We detected four member genes of the FoxP subfamily in sea lampreys (a representative species of jawless vertebrates) through genome screenings and phylogenetic analyses. Reliable paralogons (i.e. paralogous chromosome segments) have rarely been detected in scaffolds of FoxP subfamily member genes in sea lampreys due to the considerable existence of HTH_Tnp_Tc3_2 transposases. However, these transposases did not alter gene numbers of the FoxP subfamily in sea lampreys. The coincidence between the "1-4" gene duplication pattern of FoxP subfamily genes from invertebrates to vertebrates and two rounds of ancestral whole genome duplication (1R- and 2R-WGD) events reveal that the FoxP subfamily of vertebrates was quadruplicated in the 1R- and 2R-WGD events. Furthermore, we deduced that a synchronous gene duplication process occurred for the FoxP subfamily and for three linked gene families/subfamilies (i.e. MIT family, mGluR group III and PLXNA subfamily) in the 1R- and 2R-WGD events using phylogenetic analyses and mirror-dendrogram methods (i.e. algorithms to test protein-protein interactions). Specifically, the ancestor of FoxP1 and FoxP3 and the ancestor of FoxP2 and FoxP4 were generated in 1R-WGD event. In the subsequent 2R-WGD event, these two ancestral genes were changed into FoxP1, FoxP2, FoxP3 and FoxP4. The elucidation of these gene duplication processes shed light on the phylogenetic relationships between functional modules of the FoxP subfamily member genes.

The 'Tully monster' is a vertebrate.

Problematic fossils, extinct taxa of enigmatic morphology that cannot be assigned to a known major group, were once a major issue in palaeontology. A long-favoured solution to the 'problem of the problematica', particularly the 'weird wonders' of the Cambrian Burgess Shale, was to consider them representatives of extinct phyla. A combination of new evidence and modern approaches to phylogenetic analysis has now resolved the affinities of most of these forms. Perhaps the most notable exception is Tullimonstrum gregarium, popularly known as the Tully monster, a large soft-bodied organism from the late Carboniferous Mazon Creek biota (approximately 309-307 million years ago) of Illinois, USA, which was designated the official state fossil of Illinois in 1989. Its phylogenetic position has remained uncertain and it has been compared with nemerteans, polychaetes, gastropods, conodonts, and the stem arthropod Opabinia. Here we review the morphology of Tullimonstrum based on an analysis of more than 1,200 specimens. We find that the anterior proboscis ends in a buccal apparatus containing teeth, the eyes project laterally on a long rigid bar, and the elongate segmented body bears a caudal fin with dorsal and ventral lobes. We describe new evidence for a notochord, cartilaginous arcualia, gill pouches, articulations within the proboscis, and multiple tooth rows adjacent to the mouth. This combination of characters, supported by phylogenetic analysis, identifies Tullimonstrum as a vertebrate, and places it on the stem lineage to lampreys (Petromyzontida). In addition to increasing the known morphological disparity of extinct lampreys, a chordate affinity for T. gregarium resolves the nature of a soft-bodied fossil which has been debated for more than 50 years.

European Lampreys: New Insights on Postglacial Colonization, Gene Flow and Speciation.

Ice ages are known to be the most dominant palaeoclimatic feature occurring on Earth, producing severe climatic oscillations and consequently shaping the distribution and the population structure of several species. Lampreys constitute excellent models to study the colonization of freshwater systems, as they commonly appear in pairs of closely related species of anadromous versus freshwater resident adults, thus having the ability to colonize new habitats, through the anadromous species, and establish freshwater resident derivates. We used 10 microsatellite loci to investigate the spatial structure, patterns of gene flow and migration routes of Lampetra populations in Europe. We sampled 11 populations including the migratory L. fluviatilis and four resident species, L. planeri, L. alavariensis, L. auremensis and L. lusitanica, the last three endemic to the Iberian Peninsula. In this southern glacial refugium almost all sampled populations represent a distinct genetic cluster, showing high levels of allopatric differentiation, reflecting long periods of isolation. As result of their more recent common ancestor, populations from northern Europe are less divergent among them, they are represented by fewer genetic clusters, and there is evidence of strong recent gene flow among populations. These previously glaciated areas from northern Europe may have been colonized from lampreys expanding out of the Iberian refugia. The pair L. fluviatilis/L. planeri is apparently at different stages of speciation in different locations, showing evidences of high reproductive isolation in the southern refugium, and low differentiation in the north.

The p53-Mdm2 interaction and the E3 ligase activity of Mdm2/Mdm4 are conserved from lampreys to humans.

The extant jawless vertebrates, represented by lampreys and hagfish, are the oldest group of vertebrates and provide an interesting genomic evolutionary pivot point between invertebrates and jawed vertebrates. Through genome analysis of one of these jawless vertebrates, the Japanese lamprey (Lethenteron japonicum), we identified all three members of the important p53 transcription factor family--Tp53, Tp63, and Tp73--as well as the Mdm2 and Mdm4 genes. These genes and their products are significant cellular regulators in human cancer, and further examination of their roles in this most distant vertebrate relative sheds light on their origin and coevolution. Their important role in response to DNA damage has been highlighted by the discovery of multiple copies of the Tp53 gene in elephants. Expression of lamprey p53, Mdm2, and Mdm4 proteins in mammalian cells reveals that the p53-Mdm2 interaction and the Mdm2/Mdm4 E3 ligase activity existed in the common ancestor of vertebrates and have been conserved for >500 million years of vertebrate evolution. Lamprey Mdm2 degrades human p53 with great efficiency, but this interaction is not blocked by currently available small molecule inhibitors of the human HDM2 protein, suggesting utility of lamprey Mdm2 in the study of the human p53 signaling pathway.

The mitogenomes of the pouched lamprey (Geotria australis) and least brook lamprey (Lampetra aepyptera) with phylogenetic considerations.

We report the mitogenomes of the pouched lamprey (Geotria australis) and least brook lamprey (Lampetra aepyptera) in the families Geotriidae and Petromyzontidae, respectively. Both of the mitogenomes contain the 37 typical vertebrate genes. Their gene order and contents are identical to those of previously described lamprey mitogenomes. The mitogenome of G. australis (17 080 bp) is the largest among the 10 reported lamprey mitogenomes, owed to two long noncoding regions. The mitogenome of L. aepyptera is 77 bp longer (16 236 bp) than that of the congeneric European river lamprey L. fluviatilis, a size difference mostly due to different copy numbers of tandem repeats in the noncoding regions. The phylogenetic analysis supports that the pouched lamprey (Geotriidae) diverged earlier from the common ancestor of lampreys than the Petromyzonids, and the placement of the least brook lamprey in the genus Lampetra.

Using stable isotopes and C:N ratios to examine the life-history strategies and nutritional sources of larval lampreys.

Natural abundance stable-isotope analysis (δ(13)C and δ(15)N) and C:N ratios were used to study the ammocoete phase of two common non-parasitic lamprey species (least brook lamprey Lampetra aepyptera and American brook lamprey Lethenteron appendix) in two tributaries of the Ohio River (U.S.A.). The C:N ratios suggest that each species employs different lipid accumulation strategies to support its metamorphosis and recruitment into an adult animal. Ammocoete δ(13)C values generally increased with increasing C:N values. In contrast to δ(13)C, ammocoete δ(15)N values were weakly related to the total length (LT) in L. aepyptera, but positively correlated to both LT and C:N ratios in L. appendix. In L. appendix, C:N also correlated positively with LT, and presumably age. A Bayesian mixing model using δ(13)C and δ(15)N was used to estimate nutritional subsidies of different potential food resources to ammocoetes at each site. The models suggested that although nutritional subsidies to ammocoetes varied as a function of site, ammocoetes were generally reliant on large contributions (42-62% at three sites) from aquatic plants. Contributions from aquatic sediment organic matter were also important at all sites (32-63%) for ammocoetes, with terrestrially derived plant materials contributing smaller amounts (4-33%). These findings provide important insights into the feeding ecology and nutrition of two species of lampreys. They also suggest that similar and other quantitative approaches are required to (1) fully understand how the observed stable-isotopes ratios are established in ammocoetes and (2) better assess ammocoete nutritional subsidies in different natal streams.

Low reproductive isolation and highly variable levels of gene flow reveal limited progress towards speciation between European river and brook lampreys.

Ecologically based divergent selection is a factor that could drive reproductive isolation even in the presence of gene flow. Population pairs arrayed along a continuum of divergence provide a good opportunity to address this issue. Here, we used a combination of mating trials, experimental crosses and population genetic analyses to investigate the evolution of reproductive isolation between two closely related species of lampreys with distinct life histories. We used microsatellite markers to genotype over 1000 individuals of the migratory parasitic river lamprey (Lampetra fluviatilis) and freshwater-resident nonparasitic brook lamprey (Lampetra planeri) distributed in 10 sympatric and parapatric population pairs in France. Mating trials, parentage analyses and artificial fertilizations demonstrated a low level of reproductive isolation between species even though size-assortative mating may contribute to isolation. Most parapatric population pairs were strongly differentiated due to the joint effects of geographic distance and barriers to migration. In contrast, we found variable levels of gene flow between sympatric populations ranging from panmixia to moderate differentiation, which indicates a gradient of divergence with some population pairs that may correspond to alternative morphs or ecotypes of a single species and others that remain partially isolated. Ecologically based divergent selection may explain these variable levels of divergence among sympatric population pairs, but incomplete genome swamping following secondary contact could have also played a role. Overall, this study illustrates how highly differentiated phenotypes can be maintained despite high levels of gene flow that limit the progress towards speciation.

Genome-wide analysis of the ATP-binding cassette (ABC) transporter gene family in sea lamprey and Japanese lamprey.

BACKGROUND: Lampreys are extant representatives of the jawless vertebrate lineage that diverged from jawed vertebrates around 500 million years ago. Lamprey genomes contain information crucial for understanding the evolution of gene families in vertebrates. The ATP-binding cassette (ABC) gene family is found from prokaryotes to eukaryotes. The recent availability of two lamprey draft genomes from sea lamprey Petromyzon marinus and Japanese lamprey Lethenteron japonicum presents an opportunity to infer early evolutionary events of ABC genes in vertebrates. RESULTS: We conducted a genome-wide survey of the ABC gene family in two lamprey draft genomes. A total of 37 ABC transporters were identified and classified into seven subfamilies; namely seven ABCA genes, 10 ABCB genes, 10 ABCC genes, three ABCD genes, one ABCE gene, three ABCF genes, and three ABCG genes. The ABCA subfamily has expanded from three genes in sea squirts, seven and nine in lampreys and zebrafish, to 13 and 16 in human and mouse. Conversely, the multiple copies of ABCB1-, ABCG1-, and ABCG2-like genes found in sea squirts have contracted in the other species examined. ABCB2 and ABCB3 seem to be new additions in gnathostomes (not in sea squirts or lampreys), which coincides with the emergence of the gnathostome-specific adaptive immune system. All the genes in the ABCD, ABCE and ABCF subfamilies were conserved and had undergone limited duplication and loss events. In the sea lamprey transcriptomes, the ABCE and ABCF gene subfamilies were ubiquitously and highly expressed in all tissues while the members in other gene subfamilies were differentially expressed. CONCLUSIONS: Thirteen more lamprey ABC transporter genes were identified in this study compared with a previous study. By concatenating the same gene sequences from the two lampreys, more full length sequences were obtained, which significantly improved both the assignment of gene names and the phylogenetic trees compared with a previous analysis using partial sequences. The ABC gene subfamilies in chordates have undergone obvious expansion or contraction. The ABCA subfamily showed the highest gene expansion rate during chordate evolution. The evolution of ABC transporters in lampreys requires further evaluation because the present results are based on a draft genome.