Mitochondrial Non-coding Region Research Articles

Molecular Phylogeny and Taxonomy of the Coral Genus Cyphastrea (Cnidaria, Scleractinia, Merulinidae) in Japan, With the First Records of Two Species.

The scleractinian coral genus Cyphastrea is widely distributed in the Indo-Pacific region and is common from the subtropical to the warm-temperate regions in Japan. Three new species in this genus have recently been reported from south-eastern Australia or the Red Sea. However, taxonomic and species diversity have been little studied so far in Japan. In this study, we analyzed 112 specimens of Cyphastrea collected from the subtropical to the warm-temperate regions in Japan to clarify the species diversity in the country. This analysis was based on skeletal morphological and molecular analyses using three genetic markers of the nuclear 28S rDNA, histone H3 gene, and the mitochondrial noncoding intergenic region between COI and tRNAmet. The molecular phylogenetic trees showed that our specimens are separated mainly into four clades. Considering the morphological data with the molecular phylogenetic relationships, we confirmed a total of nine species, including two species, C. magna and C. salae, recorded for the first time in Japan. Although eight out of nine species were genetically included within Cyphastrea, one species, C. agassizi, was genetically distant from all other species and was closely related to the genus Leptastrea, suggesting the return of this species to the genus to which it was originally ascribed. Two newly recorded species were reciprocally monophyletic, while the other six species (excluding C. agassizi) clustered in two clades without forming species-specific lineages, including three polyphyletic species. Thus, the species boundary between species in Cyphastrea remains unclear in most species using these three sequenced loci.

The signals of selective constraints on the mitochondrial non-coding control region: insights from comparative mitogenomics of Clupeoid fishes.

The vertebrate mitochondrial genome is characterized by an exceptional organization evolving towards a reduced size. However, the persistence of a non-coding and highly variable control region is against this evolutionary trend that is explained by the presence of conserved sequence motifs or binding sites for nuclear-organized proteins that regulate mtDNA maintenance and expression. We performed a comparative mitogenomic investigation of the non-coding control region to understand its evolutionary patterns in Clupeoid fishes which are widely distributed across oceans of the world, exhibiting exemplary evolutionary potential. We confirmed the ability of sequence flanking the conserved sequence motifs in the control region to form stable secondary structures. The existence of evolutionarily conserved secondary structures without primary structure conservation suggested the action of selective constraints towards maintaining the secondary structure. The functional secondary structure is maintained by retaining the frequency of discontinuous AT and TG repeats along with compensatory base substitutions in the stem forming regions which can be considered as a selective constraint. The nucleotide polymorphism along the flanking regions of conserved sequence motifs can be explained as errors during the enzymatic replication of secondary structure-forming repeat elements. The evidence for selective constraints on secondary structures emphasizes the role of the control region in mitogenome function. Maintenance of high frequency of discontinuous repeats can be proposed as a model of adaptive evolution against the mutations that break the secondary structure involved in the efficient regulation of mtDNA functions substantiating the efficient functioning of the control region even in a high nucleotide polymorphism environment.

Identification of Main Oyster Species and Comparison of Their Genetic Diversity in Zhejiang Coast, South of Yangtze River Estuary

Oysters are an important aquaculture species distributed worldwide, including in Zhejiang Province, located on the east coast of China. Because of the high diversity and complicated introduction history of oysters and their seedlings, there has been much disagreement regarding the origin of each species, and the dominant and indigenous species remain unclear. We sampled 16 batches of oysters from seven sites in three aquaculture bays and found two main oyster species, Crassostrea sikamea and Crassostrea angulata. The former occupied the higher intertidal zone and comprised more than 70% of the cultured oysters. Based on the cytochrome oxidase C subunit I (COI) and mitochondrial noncoding region (MNR), C. sikamea showed higher genetic diversity than C. angulata. The analysis of molecular variance among COI sequences of these species from the Xiangshan Bay populations were comparable to those of other populations and showed that most of the molecular variance was within groups, which was consistent with the low pairwise fixation index FST values. The neutrality test revealed that C. sikamea experienced population expansion events, whereas for C. angulata, the significant Fu’s Fs and non-significant Tajima’s D test results may indicate a possible population expansion event, implying that C. sikamea is likely an indigenous species. The method established based on internal transcribed spacer 1 digestion by the HindIII restriction enzyme is useful for identifying C. sikamea and C. angulata in the local region. The specific primers on the MNR sequence show potential for distinguishing C. sikamea from four other important Crassostrea oysters. These results highlight the abundance of C. sikamea on the Zhejiang coast and lay the foundation for protecting and utilizing the local oyster germplasm resources and for the sustainable development of the oyster industry.

Evolutionary biogeography of the reef-building coral genus Galaxea across the Indo-Pacific ocean

Stony corals (Scleractinia) form the basis for some of the most diverse ecosytems on Earth, but we have much to learn about their evolutionary history and systematic relationships. In order to improve our understanding of species in corals we here investigated phylogenetic relationships between morphologically defined species and genetic lineages in the genus Galaxea (Euphyllidae) using a combined phylogenomic and phylogeographic approach. Previous studies revealed the nominal species G. fascicularis included three genetically well-differentiated lineages (L, S & L+) in the western Pacific, but their distribution and relationship to other species in the genus was unknown. Based on genomic (RAD-seq) and mitochondrial sequence data (non-coding region between cytb and ND2) we investigated whether the morphological taxa represent genetically coherent entities and what is the phylogenetic relationship and spatial distribution of the three lineages of G. fascicularis throughout the observed species range. Using the RAD-seq data, we find that the genus Galaxea is monophyletic and contains three distinct clades: an Indo-Pacific, a Pacific, and a small clade restricted to the Chagos Archipelago. The three lineages of G. fascicularis were associated with different RAD-seq clades, with the ‘L’ lineage showing some morphological distinction from the other two lineages (larger more asymmetrical polyps). In addition to these, three more genetic lineages in G. fascicularis may be distinguished – a Chagossian, an Ogasawaran, and one from the Indian-Red Sea. Among nominal taxa for which we have multiple samples, G. horrescens was the only monophyletic species. The mitochondrial non-coding region is highly conserved apart of the length polymorphism used to define L, S & L+ lineages and lacks the power to distinguish morphological and genetic groups resolved with genomic RAD-sequencing. The polyphyletic nature of most species warrants a careful examination of the accepted taxonomy of this group with voucher collections and their comparison to type specimens to resolve species boundaries. Further insight to the speciation process in corals will require international cooperation for the sharing of specimens to facilitate scientific discovery.

Maternal inheritance of F1 hybrid morphology and colony shape in the coral genus Acropora

BackgroundThe coral genus Acropora contains more than 150 species with very high morphological diversity. This high diversity may have been caused by repeated hybridization via mass spawning. However, we have little information whether hybrids are formed in these corals. Identifying morphological differences between hybrids and their parental species would provide an opportunity to find wild hybrids in the field and to understand how colony shapes of Acropora have become highly diversified throughout evolutionary history. In the two morphologically distinctive coral species Acropora florida and A. intermedia in the Indo-Pacific, their gametes show high rates of bi-directional intercrossing in vitro, and thus these two species are ideal species to investigate the morphological traits of the hybrids.MethodsWe examined morphological characters of F1 hybrids from A. florida to A. intermedia, which were produced from in vitro crossing experiments. To compare morphological differences, we grew juveniles and mature colonies of reciprocal F1 hybrids (FLOint: A. florida eggs × A. intermedia sperm, and INTflo: A. intermedia eggs × A. florida sperm) and of the parental species (purebreds of A. intermedia and A. florida). We analyzed skeletal morphology such as colony size, branch length, and branching number, and compared them with those of a putative F1 hybrid between A. florida and A. intermedia found in the field. We also confirmed the molecular phylogenetic position of F1 hybrids, parental species, and a putative F1 hybrid using the mitochondrial non-coding region.ResultsOur morphological analysis revealed that branching number of the F1 hybrids was intermediate relative to the parental species. Moreover, the FLOint hybrids were morphologically more closely related to the maternal species A. florida, and the INTflo hybrids were to A. intermedia. Molecular data showed that A. florida and A. intermedia were clearly divided into two clades, and that F1 hybrids grouped in the clade based on their maternal parent. A very similar pattern to the INTflo hybrids was obtained for the putative F1 hybrid in nature.DiscussionOur results revealed that F1 hybrids between two Indo-Pacific species A. florida and A. intermedia had intermediate morphology relative to their parent species but reflected the maternal parent more. Similarity to maternal species in hybrids is opposite to the Caribbean Acropora species that had more paternal morphological characters in hybrids. These results further suggest that some genetic factor in eggs is likely to affect determination of colony shape in the Indo-Pacific. At present, we have considered colonies with intermediate morphs between different species to be intra-specific morphological variation, but they may be real F1 hybrids. Indeed, a putative F1 hybrid represented similar morphological and molecular features to the F1 hybrids, and thus it is plausible to be attributed as a “real” F1 hybrid in nature.

Novel Polymorphic Microsatellite Markers Reveal Genetic Differentiation between Two Sympatric Types of Galaxea fascicularis.

The reef-building, scleractinian coral, Galaxea fascicularis, is classified into soft and hard types, based on nematocyst morphology. This character is correlated with the length of the mitochondrial non-coding region (mt-Long: soft colony type, and nematocysts with wide capsules and long shafts; mt-Short: hard colony type, and nematocysts with thin capsules and short shafts). We isolated and characterized novel polymorphic microsatellite markers for G. fascicularis using next-generation sequencing. Based upon the mitochondrial non-coding region, 53 of the 97 colonies collected were mt-Long (mt-L) and 44 were mt-Short (mt-S). Among the 53 mt-L colonies, 27 loci were identified as amplifiable, polymorphic microsatellite loci, devoid of somatic mutations and free of scoring errors. Eleven of those 27 loci were also amplifiable and polymorphic in the 44 mt-S colonies; these 11 are cross-type microsatellite loci. The other 16 loci were considered useful only for mt-L colonies. These 27 loci identified 10 multilocus lineages (MLLs) among the 53 mt-L colonies (N MLL/N = 0.189), and the 11 cross-type loci identified 7 MLLs in 44 mt-S colonies (N MLL/N = 0.159). Significant genetic differentiation between the two types was detected based on the genetic differentiation index (F ST = 0.080, P = 0.001). Bayesian clustering also indicated that these two types are genetically isolated. While nuclear microsatellite genotypes also showed genetic differentiation between mitochondrial types, the mechanism of divergence is not yet clear. These markers will be useful to estimate genetic diversity, differentiation, and connectivity among populations, and to understand evolutionary processes, including divergence of types in G. fascicularis.

Different stocks of brook lamprey in Spain and their origin from Lampetra fluviatilis at two distinct times and places.

Using the mitochondrial non-coding region I, it was shown that the two Spanish Lampetra planeri populations (Cares-Deva and Olabidea-Ugarana) correspond to different genetic units. The Cares-Deva population is probably a recent offshoot of Lampetra fluviatilis, being the Olabidea-Ugarana population less diverse and of older origin.

Haplotypes and variable position detection in the mitochondrial DNA coding region encompassing nucleotide positions 10,716–11,184

This study evaluates the mitochondrial noncoding regions by using the Sanger sequencing method for application in Forensic Science. FTA® Technology (FTA™ paper DNA extraction) was utilized to extract DNA. Portion of coding region encompassing positions from (10,716 to 11,184) amplified in accordance with the Anderson reference sequence. PCR products purified by EZ-10 spin column were then sequenced and detected using the ABI 3730 × L DNA Analyzer. A new polymorphic positions 10,750 and 10,790 that are described may be suitable sources in future for identification purpose. The data obtained can be used to identify variable nucleotide positions characterized by frequent occurrence, most promising for identification variants.

Species Boundaries ofAstreoporaCorals (Scleractinia, Acroporidae) Inferred by Mitochondrial and Nuclear Molecular Markers

The genus Astreopora is a small but ancestral group in Acroporidae, which is one of the most diverse and dominant families of scleractinian coral in Indo-Pacific reefs. We estimated the species boundaries of Astreopora corals using two molecular markers: a mitochondrial non-coding region and a nuclear ribosomal 5.8S region. Seven species (59 specimens) commonly observed around Japan (Astreopora expansa, A. gracilis, A. incrustans, A. listeri, A. myriophthalma, A. cf. suggesta, and Astreopora sp.1) were investigated, and we observed no genetic divergence in the mitochondrial marker, suggesting that these species are closely related, consistent with a species complex or recent divergence, although genotyping by the marker is not so sensitive. In the nuclear 5.8S region, 121 clones consisted of six species were divided into the four major genetic groups. Although there were no monophyletic clades, the two dominant species A. myriophthalma and A. gracilis rarely shared the same haplotypes, suggesting that gene flow is limited between them. However, A. incrustans frequently shared the same haplotypes with A. gracilis although the distributions do not overlap. We found that the ancestral genus Astreopora in Acroporidae shows less genetic variation than traditionally identified morphospecies. Although further research on fertilization rate among these species is required to determine if there are reproductive barriers, the low level of genetic diversification in this genus hints that some ecological differences among acroporid corals play a role in the evolution of scleractinian corals, considering that the other members of this family, Acropora and Montipora, are highly diversified.

Contrasts in the phylogeography of two migratory lampreys in western Europe

The ways that organisms respond to climatic oscillations depend on a wide range of factors, including life-history, behaviour, thermal preferences and physiology, and ecology. To investigate these processes, we compared the phylogeographic patterns found in Lampetra fluviatilis and Petromyzon marinus in Europe. We reanalysed all molecular data available for the mitochondrial non-coding region subunit I for both species. For L. fluviatilis, we also analyzed new and existing data for ATPase subunits 6 and 8, which cover a wider geographic range. In L. fluviatilis, both gene diversity and nucleotide diversity are at least three times higher than in P. marinus in Europe. L. fluviatilis shows population differentiation in Europe and displays a deeper haplotype network, with no predominance of an ancestral haplotype, which contrasts with a star-like pattern for P. marinus. Bayesian skyline plots for the two species fit exponential models and, with estimates of the times to the most recent common ancestor in each species, indicate that P. marinus has much younger populations in Europe, supporting the hypothesis of its relatively recent migration from North America. The differences in phylogeographic structures of these two species are discussed considering the likely effects of differences in their thermal preferences, migration abilities, and times available for diversification.

Mitochondrial genome rearrangement in a hydrothermal vent-endemic lineage of provannid gastropods provides a new DNA marker for phylogeographical studies

The nucleotide sequences of mitochondrial DNA fragments containing genes for cytochromecoxidase subunits I and II, ATPase subunits 6 and 8, large and small subunits of rRNA, and 9 tRNAs were determined using 2 genetically distinct representative specimens of the hydrothermal vent-endemic genus,Alviniconcha, and for the hydrothermal vent-endemic species,Ifremeria nautilei(Provannidae: Caenogastropoda: Mollusca). All 3 specimens showed tRNATrpgene translocation. To determine the time point at which this translocation occurred, we determined the nucleotide sequences of fragments containing 7 tRNA genes for another 3Alviniconchaand 2Ifremeriaspecimens, as well as for 7 related species. Our findings indicated that the translocation occurred in a common ancestor of theAlviniconchaandIfremeriagenera. Relatively long non-coding regions were observed on both sides of the tRNATrpgene and at the original position of this gene in all theAlviniconchaandIfremeriagastropod specimens. Overall, the mitochondrial non-coding regions appear to provide much more information about genetic structure than the mitochondrial COI gene does and could therefore be a valuable marker for future studies.

Genetic relationships of brook lamprey of the genus Lampetra in a Pyrenean stream in Spain

The genetic relationships of a population of brook lamprey Lampetra planeri from Spain were analysed using a fragment of the mitochondrial non-coding region in order to describe its relationships with other European Lampetra populations. DNA sequences were obtained from ten larvae, all corresponding to a single private haplotype. Its closest haplotype differed by five mutations and was found in several samples of Lampetra fluviatlis, both from the North Sea and the Baltic. These results argue for the great conservation value of this brook lamprey population, likely originated from L. fluviatilis ancestors.

Population genetic analysis of a recent range expansion: mechanisms regulating the poleward range limit in the volcano barnacleTetraclita rubescens

As range shifts coincident with climate change have become increasingly well documented, efforts to describe the causes of range boundaries have increased. Three mechanisms-genetic impoverishment, migration load, or a physical barrier to dispersal-are well described theoretically, but the data needed to distinguish among them have rarely been collected. We describe the distribution, abundance, genetic variation, and environment of Tetraclita rubescens, an intertidal barnacle that expanded its northern range limit by several hundreds of kilometres from San Francisco, CA, USA, since the 1970s. We compare geographic variation in abundance with abiotic and biotic patterns, including sea surface temperatures and the distributions of 387 co-occurring species, and describe genetic variation in cytochrome c oxidase subunit I, mitochondrial noncoding region, and nine microsatellite loci from 27 locations between Bahia Magdalena (California Baja Sur, Mexico) and Cape Mendocino (CA, USA). We find very high gene flow, high genetic diversity, and a gradient in physical environmental variation coincident with the range limit. We infer that the primary cause of the northern range boundary in T. rubescens is migration load arising from flow of maladapted alleles into peripheral locations and that environmental change, which could have reduced selection against genotypes immigrating into the newly colonized portion of the range, is the most likely cause of the observed range expansion. Because environmental change could similarly affect all taxa in a region whose distributional limits are established by migration load, these mechanisms may be common causes of range boundaries and largely synchronous multi-species range expansions.

Polymorphism and structural features of two noncoding regions of the liver fluke Fasciola hepatica (Plathelminthes: Trematoda) mitochondrial genome

Structural characteristics and polymorphism of the long (LNR) and short (SNR) mitochondrial noncoding regions were studied in the liver fluke Fasciola hepatica. Flukes were sampled from several populations of Russia and Belarus. LNR amplification yielded a set of nine fragments, neighboring ones differing in length by one tandem repeat (85 bp), published for Australian flukes. The LNR amplification products of different lengths were cloned and sequenced. A comparison of the LNR sequences of Australian and Belarussian flukes revealed three nucleotide substitutions and one point heteroplasmy in the first positions of the imperfect repeat and four adjacent perfect repeats. The positions of the three mutations coincided in the perfect and imperfect repeats. The frequency of mutations was 4.0–4.7 %, while the frequency of heteroplasmic sites varied from 0.1 to 1.2%. It was shown that the mutations and the heteroplasmy of one site could change the structure and stability of the putative secondary structures of the perfect and imperfect repeats. SNR amplification in F. hepatica from several populations yielded fragments that differed from the published SNR sequence of Australian F. hepatica by one transversion (T → G in position 21). Both noncoding regions had several conserved and potential regulatory sequences. The possible causes of heteroplasmy and a concerted origin of substitutions in different repeats are discussed.

The population history of the deep-sea demersal fish Bothrocara hollandi in the Japan Sea revealed by tandem repeat units in the mitochondrial non-coding region

The population history of the deep-sea demersal fish Bothrocara hollandi in the Japan Sea revealed by tandem repeat units in the mitochondrial non-coding region

Parallel Rate Heterogeneity in Chloroplast and Mitochondrial Genomes of Brazil Nut Trees (Lecythidaceae) Is Consistent with Lineage Effects

We investigated whether relative rates of divergence were correlated between the mitochondrial and chloroplast genomes as expected under lineage effects or were genome specific as expected with locus-specific effects. Five mitochondrial noncoding regions (nad1B_C, nad4exon1_2, nad7exon2_3, nad7exon3_4, and rps14-cob) for 21 samples from Lecythidaceae were sequenced. Three chloroplast regions (rpl20-5'rps12, trnS-trnG, and psbA-trnH) were sequenced to expand the taxa in an existing data set. Absolute rates of nucleotide and insertion and deletion (indel) changes were 13 times faster in the chloroplast genome than in the mitochondrial genome. Similar indel length frequency distributions for both organelles suggested that common mechanisms were responsible for generating indels. Molecular clock tests applied to phylogenetic trees estimated from mitochondrial and chloroplast sequences revealed global rate heterogeneity of nucleotide substitution. Maximum likelihood and Tajima's 1D relative rate tests show that Lecythis zabucajo exhibited a rate acceleration for both the mitochondrial and chloroplast sequences. Whereas Eschweilera romeu-cardosoi showed a significant rate slowdown for chloroplast sequences, the mitochondrial sequences for 3 Eschweilera taxa showed evidence for a rate slowdown only when compared with L. zabucajo. Significant rate heterogeneity was also observed for indel changes in the mitochondrial genome but not for the chloroplast. The lack of mitochondrial nucleotide changes for some taxa as well as chloroplast indel homoplasy may have limited the power of relative rate tests to detect rate variation. Relative ratio tests consistently indicated rate proportionality among branch lengths between the mitochondrial and chloroplast phylogenetic trees. The relative ratio tests showed that taxa possessing rate heterogeneity had parallel relative divergence rates in both mitochondrial and chloroplast sequences as expected under lineage effects. A neutral replication-dependent model of rate heterogeneity for both nucleotide and indel changes provides a simple explanation for common patterns of rate heterogeneity across the 2 organelle genomes in Lecythidaceae. The lineage effects observed here were uncoupled from annual/perennial habit because all the species from this study are perennial.

Maternal Inheritance of Mitochondrial DNA (mtDNA) in the Pacific oyster (Crassostrea gigas): a Preliminary Study Using mtDNA Sequence Analysis with Evidence of Random Distribution of MitoTracker-Stained Sperm Mitochondria in Fertilized Eggs

In many bivalve species, paternal and maternal mitochondrial DNA (mtDNA) from sperm and eggs is transmitted to the offspring. This phenomenon is known as doubly uniparental inheritance (DUI). In these species, sperm mtDNA (M type) is inherited by the male gonad of the offspring. Egg mtDNA (F type) is inherited by both male and female somatic cells and female gonadal cells. In Mytilidae, sperm mitochondria are distributed in the cytoplasm of differentiating male germ cells because they are transmitted to the male gonad. In the present study, we investigated maternal inheritance of mtDNA in the Pacific oyster, Crassostrea gigas. Sequence analysis of two mitochondrial non-coding regions revealed an identical sequence pattern in the gametes and adductor muscle samples taken from six males and five females. To observe whether sperm mitochondria were specifically located in the cytoplasm of differentiating germ cells, their distribution was recorded in C. gigas fertilized eggs by vital staining with MitoTracker Green. Although the 1D blastomere was identified in the cytoplasm of differentiating germ cells, sperm mitochondria were located at the 1D blastomere in only 32% of eggs during the 8-cell stage. Thus, in C. gigas, sperm mitochondria do not specifically locate in the germ cell region at the 1D blastomere. We suggest that the distribution of sperm mitochondria is not associated with germ cell formation in C. gigas. Furthermore, as evidenced by the mtDNA sequences of two non-coding regions, we conclude that mitochondrial DNA is maternally inherited in this species.

PCR-RFLP genotyping for Japanese and Korean populations of Pacific oyster using mitochondrial DNA noncoding region

Production of Pacific oyster Crassostrea gigas in Miyagi and Hiroshima (Japan) has gradually increased, with a marked increase in imported oysters from Goseong (Korea), and then cultured oysters of the Miyagi, Hiroshima, and Goseong populations accounting for most oyster consumption in Japan. In this study, we developed a simple PCR-RFLP analysis of a mitochondrial DNA noncoding region (mtDNA-NCR) for differentiating cultured oysters of these populations. PCR amplification yielded a 818 bp fragment comprising the entire mtDNA-NCR and parts of adjacent tRNACys and tRNAAsn genes from all the specimens. By use of a wide restriction test using 30 different restriction enzymes, only a single enzyme only, AluI, produced a unique RFLP pattern enabling us to discriminate Miyagi and Hiroshima oysters from Goseong oysters. This difference is probably because of nucleotide alteration at the presumptive AluI recognizing site on position 439 of the mtDNA-NCR. Our simple, robust and cost-effective PCR-RFLP analysis is potentially useful for population genetic investigation of cultured Pacific oyster, particularly when large numbers of specimens must be analyzed.

GEOGRAPHIC DIFFERENCES IN SPECIES BOUNDARIES AMONG MEMBERS OF THE MONTASTRAEA ANNULARIS COMPLEX BASED ON MOLECULAR AND MORPHOLOGICAL MARKERS

The three members of the Montastraea annularis complex (M. annularis, M. franksi, and M. faveolata) are dominant reef builders in the western Atlantic whose species status has been controversial for over a decade. Although differences in colony morphology and reproductive characteristics exist, interspecific fertilizations are possible in the laboratory and genetic differentiation is slight. Here we compare the three taxa genetically and morphologically in Panama and the Bahamas, widely separated locations spanning most of their geographic ranges. In Panama, analyses of three AFLP loci, a noncoding region of the mitochondrial genome, and ITS sequences reveal that M. faveolata is strongly differentiated genetically. Discriminant function analysis also indicates no overlap with the other two species in the fine structure of the corallites that comprise the colony. Genetic analyses of larvae from interspecific crosses between M. faveolata and the other two taxa confirmed the hybrid status of the larvae, but no examples of the most probable F1 genotype were observed in the field. Although M. annularis and M. franksi were more similar, they also exhibited strong frequency differences at two AFLP loci and in the mitochondrial noncoding region, as well as distinct corallite structure. In the Bahamas, in contrast, the three taxa exhibited overlapping morphologies. Montastraeafranksi and M. annularis were indistinguishable genetically, and M. faveolata was distinct at fewer genetic loci. Once again, however, the most probable F1 genotype involving M. faveolata was not observed. Geographic differences between Panama and the Bahamas explain why past studies have come to different conclusions concerning the status of the three species. In general, the genetic and morphological data suggest a north to south hybridization gradient, with evidence for introgression strongest in the north. However, reproductive data show no such trend, with intrinsic barriers to gene flow comparable or stronger in the north.

GEOGRAPHIC DIFFERENCES IN SPECIES BOUNDARIES AMONG MEMBERS OF THE MONTASTRAEA ANNULARIS COMPLEX BASED ON MOLECULAR AND MORPHOLOGICAL MARKERS

The three members of the Montastraea annularis complex (M. annularis, M. franksi, and M. faveolata) are dominant reef builders in the western Atlantic whose species status has been controversial for over a decade. Although differences in colony morphology and reproductive characteristics exist, interspecific fertilizations are possible in the laboratory and genetic differentiation is slight. Here we compare the three taxa genetically and morphologically in Panama and the Bahamas, widely separated locations spanning most of their geographic ranges. In Panama, analyses of three AFLP loci, a noncoding region of the mitochondrial genome, and ITS sequences reveal that M. faveolata is strongly differentiated genetically. Discriminant function analysis also indicates no overlap with the other two species in the fine structure of the corallites that comprise the colony. Genetic analyses of larvae from interspecific crosses between M. faveolata and the other two taxa confirmed the hybrid status of the larvae, but no examples of the most probable F1 genotype were observed in the field. Although M. annularis and M. franksi were more similar, they also exhibited strong frequency differences at two AFLP loci and in the mitochondrial noncoding region, as well as distinct corallite structure. In the Bahamas, in contrast, the three taxa exhibited overlapping morphologies. Montastraea franksi and M. annularis were indistinguishable genetically, and M. faveolata was distinct at fewer genetic loci. Once again, however, the most probable F1 genotype involving M. faveolata was not observed. Geographic differences between Panama and the Bahamas explain why past studies have come to different conclusions concerning the status of the three species. In general, the genetic and morphological data suggest a north to south hybridization gradient, with evidence for introgression strongest in the north. However, reproductive data show no such trend, with intrinsic barriers to gene flow comparable or stronger in the north.