A hybrid zone between hydrothermal vent mussels (Bivalvia: Mytilidae) from the Mid-Atlantic Ridge

This study provides the first example of a hybrid zone between animal taxa distributed along the mid-ocean ridge system. We examined the distribution and genetic structure of deep-sea hydrothermal vent mussels (Bivalvia: Mytilidae) along a 2888-km portion of the Mid-Atlantic Ridge between 37 degrees 50' N and 14 degrees 45' N latitude. Mitochondrial DNA (mtDNA), allozymes and multivariate-morphometric evidence discriminated between individuals of a northern species, Bathymodiolus azoricus, and a southern species, B. puteoserpentis, that were separated by an intermediate ridge segment almost devoid of mussels. A small sample of mussels from Broken Spur, a vent locality along this intermediate zone, revealed a mixed population with gene frequencies and morphology that were broadly intermediate to those of the northern and southern species. Multilocus clines in mtDNA and allozyme frequencies were centred over the intermediate zone. We consider intrinsic and extrinsic processes that might limit genetic exchange across this hybrid zone.     

Distinct patterns of genetic differentiation among annelids of eastern Pacific hydrothermal vents

Population genetic and phylogenetic analyses of mitochondrial COI from five deep-sea hydrothermal vent annelids provided insights into their dispersal modes and barriers to gene flow. These polychaetes inhabit vent fields located along the East Pacific Rise (EPR) and Galapagos Rift (GAR), where hundreds to thousands of kilometers can separate island-like populations. Long-distance dispersal occurs via larval stages, but larval life histories differ among these taxa. Mitochondrial gene flow between populations of Riftia pachyptila, a siboglinid worm with neutrally buoyant lecithothrophic larvae, is diminished across the Easter Microplate region, which lies at the boundary of Indo-Pacific and Antarctic deep-sea provinces. Populations of the siboglinid Tevnia jerichonana are similarly subdivided. Oasisia alvinae is not found on the southern EPR, but northern EPR populations of this siboglinid are subdivided across the Rivera Fracture Zone. Mitochondrial gene flow of Alvinella pompejana, an alvinellid with large negatively buoyant lecithotrophic eggs and arrested embryonic development, is unimpeded across the Easter Microplate region. Gene flow in the polynoid Branchipolynoe symmytilida also is unimpeded across the Easter Microplate region. However, A. pompejana populations are subdivided across the equator, whereas B. symmitilida populations are subdivided between the EPR and GAR axes. The present findings are compared with similar evidence from codistributed species of annelids, molluscs and crustaceans to identify potential dispersal filters in these eastern Pacific ridge systems.     

Directional dispersal between mid-ocean ridges: deep-ocean circulation and gene flow in Ridgeia piscesae

This study examined relationships between bathymetrically induced deep-ocean currents and the dispersal of the hydrothermal vent tubeworm Ridgeia piscesae along the northeast Pacific ridge system. A robust diagnostic model of deep-ocean circulation in this region predicted strong southeasterly currents following contours of the Blanco Transform Fault, a 450-km lateral offset that separates the Gorda and Juan de Fuca ridge systems. Such currents should facilitate the southward dispersal of R. piscesae larvae. Immigration rates for populations north and south of the Blanco Transform Fault were estimated from molecular population genetic data. Mitochondrial DNA evidence revealed population subdivision across the Blanco Transform Fault, and a strong directional bias in gene flow that was consistent with predictions of the circulation model. The distribution of mitochondrial diversity between the northern and southern populations of R. piscesae suggests that the Gorda Ridge tubeworms have maintained larger effective population sizes than the northern populations, a pattern that also exists in co-occurring limpets. Together, these data suggest that the northern vent fields may experience a higher frequency of habitat turnover and consequently more rapid losses of genetic diversity.     

Characterization of a novel spirochete associated with the hydrothermal vent polychaete annelid, Alvinella pompejana

A highly integrated, morphologically diverse bacterial community is associated with the dorsal surface of Alvinella pompejana, a polychaetous annelid that inhabits active high-temperature deep-sea hydrothermal vent sites along the East Pacific Rise (EPR). Analysis of a previously prepared bacterial 16S ribosomal DNA (rDNA) library identified a spirochete most closely related to an endosymbiont of the oligochete Olavius loisae. This spirochete phylotype (spirochete A) comprised only 2.2% of the 16S rDNA clone library but appeared to be much more dominant when the same sample was analyzed by denaturing gradient gel electrophoresis (DGGE) and the terminal restriction fragment length polymorphism procedure (12 to 18%). PCR amplification of the community with spirochete-specific primers used in conjunction with DGGE analysis identified two spirochete phylotypes. The first spirochete was identical to spirochete A but was present in only one A. pompejana specimen. The second spirochete (spirochete B) was 84.5% similar to spirochete A and, more interestingly, was present in the epibiont communities of all of the A. pompejana specimens sampled throughout the geographic range of the worm (13 degrees N to 32 degrees S along the EPR). The sequence variation of the spirochete B phylotype was less than 3% for the range of A. pompejana specimens tested, suggesting that a single spirochete species was present in the A. pompejana epibiotic community. Additional analysis of the environments surrounding the worm revealed that spirochetes are a ubiquitous component of high-temperature vents and may play an important role in this unique ecosystem.     

Diversity among three novel groups of hyperthermophilic deep-sea Thermococcus species from three sites in the northeastern Pacific Ocean

Eight new strains of deep-sea hyperthermophilic sulfur reducers were isolated from hydrothermal vent fields at 9 degrees 50'N East Pacific Rise (EPR) and at the Cleft and CoAxial segments along the Juan de Fuca Ridge (JdFR). 16S rRNA gene sequence analysis showed that each strain belongs to the genus Thermococcus. Restriction fragment length polymorphism patterns of the 16S/23S rRNA intergenic spacer region revealed that these isolates fell into three groups: those from the EPR, those from fluid and rock sources on the JdFR, and those isolated from Paralvinella spp. polychaete vent worms from the JdFR. The optimum-temperature specific growth rates and the temperature ranges for growth were significantly higher and broader for those strains isolated from worms relative to those isolated from low-temperature diffuse hydrothermal fluids. Furthermore, the worm-derived isolates generally produced a larger array of proteases and amylases based on zymogram analyses. The zymogram patterns also changed with growth temperature suggesting that these organisms alter their lytic protein suites in response to changes in temperature. This study suggests that there is significant phenotypic diversity in Thermococcus that is not apparent from their highly conserved 16S rRNA nucleotide sequences.     

Cytonuclear disequilibrium in a hybrid zone involving deep-sea hydrothermal vent mussels of the genus Bathymodiolus

A hybrid zone involving the deep-sea mussels, Bathymodiolus azoricus and B. puteoserpentis, was recently discovered at Broken Spur hydrothermal vent field (29 degrees 10' N, 43 degrees 10' W) along an intermediate segment of the Mid-Atlantic Ridge axis. Examination of nuclear (allozymes) and cytoplasmic (mitochondrial DNA) gene markers in a new sample from Broken Spur revealed significant cytonuclear disequilibrium caused by an excess of the parental types (coupling phase) and a deficiency of recombinants (repulsion phase). An assignment test of individual multilocus genotypes also revealed an excess of parental genotypes in the admixed population. These results support the hypothesis that the Broken Spur mussel population comprises a nonequilibrium mixture of parental immigrants and hybrid individuals.     

Comparative phylogeography among hydrothermal vent species along the East Pacific Rise reveals vicariant processes and population expansion in the South

The use of sequence polymorphism from individual mitochondrial genes to infer past demography has recently proved controversial because of the recurrence of selective sweeps acting over genes and the need for unlinked multilocus data sets. However, comparative analyses using several species for one gene and/or multiple genes for one species can serve as a test for potential selective effects and clarify our understanding of historical demographic effects. This study compares nucleotide polymorphisms in mitochondrial cytochrome oxidase I across seven deep-sea hydrothermal vent species that live along the volcanically active East Pacific Rise. Approximate Bayesian Computation (ABC) method, developed to trace shared vicariant events across species pairs, indicates the occurrence of two across species divergence times, and suggests that the present geographical patterns of genetic differentiation may be explained by two periods of significant population isolation. The oldest period dates back 11.6 Ma and is associated with the vent limpet Lepetodrilus elevatus , while the most recent period of isolation is 1.3 Ma, which apparently affected all species examined and coincides with a transition zone across the equator. Moreover, significant negative Tajima's D and star-like networks were observed for all southern lineages, suggesting that these lineages experienced a concomitant demographic and geographical expansion about 100 000–300 000 generations ago. This expansion may have initiated from a wave of range expansions during the secondary colonization of new sites along the Southern East Pacific Rise (founder effects below the equator) or recurrent bottleneck events because of the increase of eruptive phases associated with the higher spreading rates of the ridge in this region.     

GENETIC STRUCTURE OF GIANT CLAM (TRIDACNA MAXIMA) POPULATIONS IN THE WEST PACIFIC IS NOT CONSISTENT WITH DISPERSAL BY PRESENT-DAY OCEAN CURRENTS

The Pacific marine biota, particularly species with long planktonic larval stages, are thought to disperse widely throughout the Pacific via ocean currents. The little genetic data available to date has supported this view in that little or no significant regional differentiation of populations has been found over large geographical distances. However, recent data from giant clams has demonstrated not only significant regional differentiation of populations, but routes of gene flow that run perpendicular to the main present-day ocean currents. Extensive surveys of genetic variation at eight polymorphic loci in 19 populations of the giant clam Tridacna maxima, sampled throughout the West and Central Pacific, confirmed that the patterns of variation seen so far in T. gigas were not unique to that species, and may reflect a fundamental genetic structuring of shallow-water marine taxa. Populations of T. maxima within highly connected reef systems like the Great Barrier Reef were panmictic (average F_ST < 0.003), but highly significant genetic differences between reef groups on different archipelagos (average F_ST = 0.084) and between West and Central Pacific regions (average F_ST = 0.156) were found. Inferred gene flow was high (N_em usually > 5) between the Philippines and the Great Barrier Reef, between the Philippines and Melanesia (the Solomon Islands and Fiji), and between the Philippines and the Central Pacific island groups (Marshall Islands, Kiribati, Tuvalu and Cook Islands). Gene flow was low between these three sets of island chains (N_em < 2). These routes of gene flow are perpendicular to present-day ocean currents. It is suggested that the spatial patterns of gene frequencies reflect past episodes of dispersal at times of lower sea levels which have not been erased by subsequent dispersal by present-day circulation. The patterns are consistent with extensive dispersal of marine species in the Pacific, and with traditional views of dispersal from the Indo-Malay region. However, they demonstrate that dispersal along present-day ocean surface currents cannot be assumed, that other mechanisms may operate today or that major dispersal events are intermittent (perhaps separated by several thousands of years), and that the nature and timing of dispersal of Pacific marine species is more complex than has been thought.     

Gene flow of Acanthaster planci (L.) in relation to ocean currents revealed by microsatellite analysis

Population outbreaks of the coral-eating starfish, Acanthaster planci , are hypothesized to spread to many localities in the Indo-Pacific Ocean through dispersal of planktonic larvae. To elucidate the gene flow of A. planci across the Indo-Pacific in relation to ocean currents and to test the larval dispersal hypothesis, the genetic structure among 23 samples over the Indo-Pacific was analysed using seven highly polymorphic microsatellite loci. The F -statistics and genetic admixture analysis detected genetically distinct groups in accordance with ocean current systems, that is, the Southeast African group (Kenya and Mayotte), the Northwestern Pacific group (the Philippines and Japan), Palau, the North Central Pacific group (Majuro and Pohnpei), the Great Barrier Reef, Fiji, and French Polynesia, with a large genetic break between the Indian and Pacific Oceans. A pattern of significant isolation by distance was observed among all samples ( P =  0.001, r  = 0.88, n  = 253, Mantel test), indicating restricted gene flow among the samples in accordance with geographical distances. The data also indicated strong gene flow within the Southeast African, Northwestern Pacific, and Great Barrier Reef groups. These results suggest that the western boundary currents have strong influence on gene flow of this species and may trigger secondary outbreaks.     

A signal for independent coastal and continental histories among North American wolves

Relatively little genetic variation has been uncovered in surveys across North American wolf populations. Pacific Northwest coastal wolves, in particular, have never been analysed. With an emphasis on coastal Alaska wolf populations, variation at 11 microsatellite loci was assessed. Coastal wolf populations were distinctive from continental wolves and high levels of diversity were found within this isolated and relatively small geographical region. Significant genetic structure within southeast Alaska relative to other populations in the Pacific Northwest, and lack of significant correlation between genetic and geographical distances suggest that differentiation of southeast Alaska wolves may be caused by barriers to gene flow, rather than isolation by distance. Morphological research also suggests that coastal wolves differ from continental populations. A series of studies of other mammals in the region also has uncovered distinctive evolutionary histories and high levels of endemism along the Pacific coast. Divergence of these coastal wolves is consistent with the unique phylogeographical history of the biota of this region and re-emphasizes the need for continued exploration of this biota to lay a framework for thoughtful management of southeast Alaska.     

Ecological zones rather than molecular forms predict genetic differentiation in the malaria vector Anopheles gambiae s.s. in Ghana

The malaria mosquito Anopheles gambiae s.s. is rapidly becoming a model for studies on the evolution of reproductive isolation. Debate has centered on the taxonomic status of two forms (denoted M and S) within the nominal taxon identified by point mutations in the X-linked rDNA region. Evidence is accumulating that there are significant barriers to gene flow between these forms, but that the barriers are not complete throughout the entire range of their distribution. We sampled populations from across Ghana and southern Burkina Faso, West Africa, from areas where the molecular forms occurred in both sympatry and allopatry. Neither Bayesian clustering methods nor F(ST)-based analysis of microsatellite data found differentiation between the M and S molecular forms, but revealed strong differentiation among different ecological zones, irrespective of M/S status and with no detectable effect of geographical distance. Although no M/S hybrids were found in the samples, admixture analysis detected evidence of contemporary interform gene flow, arguably most pronounced in southern Ghana where forms occur sympatrically. Thus, in the sampled area of West Africa, lack of differentiation between M and S forms likely reflects substantial introgression, and ecological barriers appear to be of greater importance in restricting gene flow.     

Faunal changes and geographic crypticism indicate the occurrence of a biogeographic transition zone along the southern East Pacific Rise

Aim Deep‐sea hydrothermal vents have now been reported along all active mid‐ocean ridges and back‐arc basins, but the boundaries of biogeographic entities remain questionable owing to methodological issues. Here we examine biogeographic patterns of the vent fauna along the East Pacific Rise (EPR) and determine the relative roles of regional and local factors on the distribution of biodiversity associated with mussel beds along a poorly explored zone, the southern EPR (SEPR). Location East Pacific Rise. Methods A species list of macrobenthic invertebrates along the EPR was compiled from the literature and supplemented with data recovered during the French research cruise BIOSPEEDO carried out in 2004 along the SEPR. Biogeographic patterns were assessed by combining the identification of morphological species with a molecular barcoding approach. A multivariate regression tree (MRT) analysis was performed to identify any geographic breaks, and an empirical distribution of species richness was compared with predictions provided by a mid‐domain effect model. Macrofaunal community structure associated with mussel beds along the SEPR was analysed in relation to environmental factors using cluster and canonical redundancy analyses. Results Sequencing of the cytochrome c oxidase subunit I gene revealed the occurrence of several cryptic species complexes along the EPR, with the equator separating the southern and northern clades. Furthermore, during the BIOSPEEDO cruise at least 10 still unnamed species were collected between 7°25′ S and 21°33′ S. The shift in community structure identified by MRT analysis was located south of 17°34′ S or south of 13°59′ S, depending on the data used, suggesting that the southern part of the SEPR (17°25′–21°33′ S) constitutes a biogeographic transition zone in the vent fauna along the EPR. At a regional scale, latitude combined with the type of venting was significantly correlated with the community structure associated with mussel beds. Main conclusions Together, the molecular data, in situ observations, and the distribution of species suggest that the high diversity of vent fauna species presently observed between 17°25′ S and 21°33′ S is probably a result of the overlap of several distinct biogeographic provinces. We argue that this area thus constitutes a biogeographic vent fauna transition zone along the EPR.     

Vertical distribution of bacterial and archaeal communities along discrete layers of a deep-sea cold sediment sample at the East Pacific Rise (∼13°N)

The community structure and vertical distribution of prokaryotes in a deep-sea (ca. 3,191 m) cold sediment sample (ca. 43 cm long) collected at the East Pacific Rise (EPR) approximately 13 degrees N were studied with 16SrDNA-based molecular analyses. Total community DNA was extracted from each of four discrete layers EPRDS-1, -2, -3 and -4 (from top to bottom) and 16S rDNA were amplified by PCR. Cluster analysis of DGGE profiles revealed that the bacterial communities shifted sharply between EPRDS-1 and EPRDS-2 in similarity coefficient at merely 49%. Twenty-three sequences retrieved from DGGE bands fell into 11 groups based on BLAST and bootstrap analysis. The dominant groups in the bacterial communities were Chloroflexi, Gamma proteobacteria, Actinobacterium and unidentified bacteria, with their corresponding percentages varying along discrete layers. Pairwise Fst (F-statistics) values between the archaeal clone libraries indicated that the archaeal communities changed distinctly between EPRDS-2 and EPRDS-3. Sequences from the archaeal libraries were divided to eight groups. Crenarchaea Marine Group I (MGI) was prevalent in EPRDS-1 at 83%, while Uncultured Crenarchaea group II B (UCII B) abounded in EPRDS-4 at 61%. Our results revealed that the vertically stratified distribution of prokaryotic communities might be in response to the geochemical settings and suggested that the sampling area was influenced by hydrothermalism. The copresence of members related to hydrothermalism and cold deep-sea environments in the microbial community indicated that the area might be a transitional region from hydrothermal vents to cold deep-sea sediments.     

Population genetic subdivision in the New Zealand greenshell mussel (Perna canaliculus) inferred from single-strand conformation polymorphism analysis of mitochondrial DNA

Single-strand conformation polymorphism (SSCP) analysis of the NADH IV region of the mitochondrial DNA (mtDNA) molecule in greenshell mussels (Perna canaliculus) indicated strong population genetic structuring in this endemic New Zealand species. A northern and a southern group were differentiated by frequency shifts in common haplotypes and by the occurrence of a unique southern haplotype at approximately 20% frequency. This split occurred south of Cook Strait (the body of water between the North and the South Island) at approximately 42 degrees S latitude. Northern populations were less genetically diverse than southern populations and mussels from the west coast of the South Island were most distinct from northern mussels. We hypothesize that the unique haplotype VIII originated in the lower South Island, and that its spread northwards was obstructed by the opening of Cook Strait approximately 15 000-16 000 years ago and the subsequent establishment of present-day surface water circulation patterns in Greater Cook Strait. We suggest that present-day strong tidal flows and turbulent mixing of water masses in Cook Strait, and intense up-welling on the east and west coasts in this region, represent a barrier to gene flow between mussels located in the North Island and northern South Island vs. mussels in most of the South Island and Stewart Island.     

Geographic structure in the Southern Ocean circumpolar brittle star Ophionotus victoriae (Ophiuridae) revealed from mtDNA and single‐nucleotide polymorphism data

Marine systems have traditionally been thought of as “open” with few barriers to gene flow. In particular, many marine organisms in the Southern Ocean purportedly possess circumpolar distributions that have rarely been well verified. Here, we use the highly abundant and endemic Southern Ocean brittle star Ophionotus victoriae to examine genetic structure and determine whether barriers to gene flow have existed around the Antarctic continent. Ophionotus victoriae possesses feeding planktotrophic larvae with presumed high dispersal capability, but a previous study revealed genetic structure along the Antarctic Peninsula. To test the extent of genetic differentiation within O. victoriae, we sampled from the Ross Sea through the eastern Weddell Sea. Whereas two mitochondrial DNA markers (16S rDNA and COI) were employed to allow comparison to earlier work, a 2b‐RAD single‐nucleotide polymorphism (SNP) approach allowed sampling of loci across the genome. Mitochondrial data from 414 individuals suggested three major lineages, but 2b‐RAD data generated 1,999 biallelic loci that identified four geographically distinct groups from 89 samples. Given the greater resolution by SNP data, O. victoriae can be divided into geographically distinct populations likely representing multiple species. Specific historical scenarios that explain current population structure were examined with approximate Bayesian computation (ABC) analyses. Although the Bransfield Strait region shows high diversity possibly due to mixing, our results suggest that within the recent past, dispersal processes due to strong currents such as the Antarctic Circumpolar Current have not overcome genetic subdivision presumably due to historical isolation, questioning the idea of large open circumpolar populations in the Southern Ocean.     

Ancient trade routes shaped the genetic structure of horses in eastern Eurasia

Animal exchange networks have been shown to play an important role in determining gene flow among domestic animal populations. The Silk Road is one of the oldest continuous exchange networks in human history, yet its effectiveness in facilitating animal exchange across large geographical distances and topographically challenging landscapes has never been explicitly studied. Horses are known to have been traded along the Silk Roads; however, extensive movement of horses in connection with other human activities may have obscured the genetic signature of the Silk Roads. To investigate the role of the Silk Roads in shaping the genetic structure of horses in eastern Eurasia, we analysed microsatellite genotyping data from 455 village horses sampled from 17 locations. Using least‐cost path methods, we compared the performance of models containing the Silk Roads as corridors for gene flow with models containing single landscape features. We also determined whether the recent isolation of former Soviet Union countries from the rest of Eurasia has affected the genetic structure of our samples. The overall level of genetic differentiation was low, consistent with historically high levels of gene flow across the study region. The spatial genetic structure was characterized by a significant, albeit weak, pattern of isolation by distance across the continent with no evidence for the presence of distinct genetic clusters. Incorporating landscape features considerably improved the fit of the data; however, when we controlled for geographical distance, only the correlation between genetic differentiation and the Silk Roads remained significant, supporting the effectiveness of this ancient trade network in facilitating gene flow across large geographical distances in a topographically complex landscape.     

Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone

The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long‐standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M. galloprovincialis, across a well‐studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best‐supported scenario was a long period of allopatric isolation during the first three‐quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.     

Population genetic structure of diploid sexual and polyploid apomictic hawthorns (Crataegus; Rosaceae) in the Pacific Northwest

Polyploidy and gametophytic apomixis are two important and associated processes in plants. Many hawthorn species are polyploids and can reproduce both sexually and apomictically. However, the population genetic structure of these species is poorly understood. Crataegus douglasii is represented exclusively by self-compatible tetraploid pseudogamous apomicts across North America, whereas Crataegus suksdorfii found in the Pacific Northwest is known to include self-incompatible diploid sexuals as well as polyploid apomicts. We compare population structure and genetic variability in these two closely related taxa using microsatellite and chloroplast sequence markers. Using 13 microsatellite loci located on four linkage groups, 251 alleles were detected in 239 individuals sampled from 15 localities. Within-population multilocus genotypic variation and molecular diversity are greatest in diploid sexuals and lowest in triploid apomicts. Apart from the isolation of eastern North American populations of C. douglasii , there is little evidence of isolation by distance in this taxon. Genetic diversity in western populations of C. douglasii suggests that gene flow is frequent, and that colonization and establishment are often successful. In contrast, local populations of C. suksdorfii are more markedly differentiated. Gene flow appears to be limited primarily by distance in diploids and by apomixis and self-compatibility in polyploids. We infer that apomixis and reproductive barriers between cytotypes are factors that reduce the frequency of gene flow among populations, and may ultimately lead to allopatric speciation in C. suksdorfii . Our findings shed light on evolution in woody plants that show heterogeneous ploidy levels and reproductive systems.     

Phylogeography and genetic connectivity of the marine macro‐alga Sargassum ilicifolium (Phaeophyceae,Ochrophyta) in the northwestern Pacific1

The evolutionary influences of historical and contemporary factors on the population connectivity and phylogeographic structure of a brown seaweed, Sargassum ilicifolium, were elucidated using the nuclear ITS2 and mitochondrial COI markers for the collections newly sampled within its distribution range in the northwestern Pacific (NWP). Significant genetic structure at variable levels was identified between populations (pairwise F_ST) and among populations grouped by geographical proximity (Φ_CT among regions). The adjacent groups of populations with moderate structure revealed from AMOVA appeared to have high genetic connectivity. However, a lack of genealogical concordance with the geographic distribution was uncovered for S. ilicifolium from the NWP. Such genetic homogeneity is interpreted as a result of the interaction between postglacial recolonization and dynamic oceanic current regimes in the region. Two separated glacial refugia, the South China Sea and the Okinawa Trough, in the marginal seas of east China were recognized based on the presence of endemic haplotypes and high haplotype diversity in the populations at southern China and northeast of Taiwan. Populations persisting in these refugia may have served as the source for recolonization in the NWP with the rise of sea level during the warmer interglacial periods. The role of oceanic currents in maintaining genetic connectivity of S. ilicifolium in the region was further corroborated by the coherence between the direction of oceanic currents and that of gene flow, especially along the eastern coast of Taiwan. This study underlines the interaction between historical postglacial recolonization and contemporary coastal hydrodynamics in contributing to population connectivity and distribution for this tropical seaweed in the NWP.