Multiple-geographic-scale genetic structure of two mangrove tree species: the roles of mating system, hybridization, limited dispersal and extrinsic factors

Chromosome-Level Genome Assembly and Population Genomic Analyses Reveal Geographic Variation and Population Genetic Structure of Prunus tenella

Prunus tenella is a rare and precious relict plant in China. It is an important genetic resource for almond improvement and an indispensable material in ecological protection and landscaping. However, the research into molecular breeding and genetic evolution has been severely restricted due to the lack of genome information. In this investigation, we created a chromosome-level genomic pattern of P. tenella, 231 Mb in length with a contig N50 of 18.1 Mb by Hi-C techniques and high-accuracy PacBio HiFi sequencing. The present assembly predicted 32,088 protein-coding genes, and an examination of the genome assembly indicated that 94.7% among all assembled transcripts were alignable to the genome assembly; most (97.24%) were functionally annotated. By phylogenomic genome comparison, we found that P. tenella is an ancient group that diverged approximately 13.4 million years ago (mya) from 13 additional closely related species and about 6.5 Mya from the cultivated almond. Collinearity analysis revealed that P. tenella is highly syntenic and has high sequence conservation with almond and peach. However, this species also exhibits many presence/absence variants. Moreover, a large inversion at the 7588 kb position of chromosome 5 was observed, which may have a significant association with phenotypic traits. Lastly, population genetic structure analysis in eight different populations indicated a high genetic differentiation among the natural distribution of P. tenella. This high-quality genome assembly provides critical clues and comprehensive information for the systematic evolution, genetic characteristics, and functional gene research of P. tenella. Moreover, it provides a valuable genomic resource for in-depth study in protection, developing, and utilizing P. tenella germplasm resources.

Out to sea: ocean currents and patterns of asymmetric gene flow in an intertidal fish species

Passive dispersal via wind or ocean currents can drive asymmetric gene flow, which influences patterns of genetic variation and the capacity of populations to evolve in response to environmental change. The mangrove rivulus fish (Kryptolebias marmoratus), hereafter "rivulus," is an intertidal fish species restricted to the highly fragmented New World mangrove forests of Central America, the Caribbean, the Bahamas, and Florida. Mangrove patches are biological islands with dramatic differences in both abiotic and biotic conditions compared to adjacent habitat. Over 1,000 individual rivulus across 17 populations throughout its range were genotyped at 32 highly polymorphic microsatellites. Range-wide population genetic structure was evaluated with five complementary approaches that found eight distinct population clusters. However, an analysis of molecular variance indicated significant population genetic structure among regions, populations within regions, sampling locations within populations, and individuals within sampling locations, indicating that rivulus has both broad- and fine-scale genetic differentiation. Integrating range-wide genetic data with biophysical modeling based on 10 years of ocean current data showed that ocean currents and the distance between populations over water drive gene flow patterns on broad scales. Directional migration estimates suggested some significant asymmetries in gene flow that also were mediated by ocean currents and distance. Specifically, populations in the center of the range (Florida Keys) were identified as sinks that received migrants (and alleles) from other populations but failed to export individuals. These populations thus harbor genetic variation, perhaps even from extirpated populations across the range, but ocean currents and complex arrangements of landmasses might prevent the distribution of that genetic variation elsewhere. Hence, the inherent asymmetry of ocean currents shown to impact both genetic differentiation and directional migration rates may be responsible for the complex distribution of genetic variation across the range and observed patterns of metapopulation structure.

The role of oceanic currents in the dispersal and connectivity of the mangrove Rhizophora mangle on the Southwest Atlantic region

Dispersal is a crucial mechanism to living beings, allowing them to reach new resources such that populations and species can occupy new environments. However, directly observing the dispersal mechanisms of widespread species can be costly or even impractical, which is the case for mangrove trees. The influence of ocean currents on mangrove dispersal is increasingly evident; however, few studies mechanistically relate the patterns of population distribution with the dispersal by oceanic currents under an integrated framework. Here, we evaluate the role of oceanic currents on connectivity of Rhizophora mangle along the Southwest Atlantic. We inferred population genetic structure and migration rates, simulated the displacement of propagules and tested our hypotheses with Mantel tests and redundancy analysis. We observed populations structured in two major groups, north and south, which is corroborated by other studies with Rhizophora and other coastal plants. Inferred recent migration rates do not indicate ongoing gene flow between sites. Conversely, long-term migration rates were low across groups and contrasting dispersal patterns within each one, which is consistent with long-distance dispersal events. Our hypothesis tests suggest that both isolation by distance and isolation by oceanography (derived from the oceanic currents) can explain the neutral genetic variation of R. mangle in the region. Our findings expand current knowledge of mangrove connectivity and highlight how the association of molecular methods with oceanographic simulations improve the interpretation of the dispersal process. This integrative approach is a cost- and time-efficient strategy to include dispersal and connectivity data into marine protected areas planning and management.

Oceanographic connectivity explains the intra-specific diversity of mangrove forests at global scales

The distribution of mangrove intra-specific biodiversity can be structured by historical demographic processes that enhance or limit effective population sizes. Oceanographic connectivity (OC) may further structure intra-specific biodiversity by preserving or diluting the genetic signatures of historical changes. Despite its relevance for biogeography and evolution, the role of oceanographic connectivity in structuring the distribution of mangrove's genetic diversity has not been addressed at global scale. Here we ask whether connectivity mediated by ocean currents explains the intra-specific diversity of mangroves. A comprehensive dataset of population genetic differentiation was compiled from the literature. Multigenerational connectivity and population centrality indices were estimated with biophysical modeling coupled with network analyses. The variability explained in genetic differentiation was tested with competitive regression models built upon classical isolation-by-distance (IBD) models considering geographic distance. We show that oceanographic connectivity can explain the genetic differentiation of mangrove populations regardless of the species, region, and genetic marker (significant regression models in 95% of cases, with an average R-square of 0.44 ± 0.23 and Person's correlation of 0.65 ± 0.17), systematically improving IBD models. Centrality indices, providing information on important stepping-stone sites between biogeographic regions, were also important in explaining differentiation (R-square improvement of 0.06 ± 0.07, up to 0.42). We further show that ocean currents produce skewed dispersal kernels for mangroves, highlighting the role of rare long-distance dispersal events responsible for historical settlements. Overall, we demonstrate the role of oceanographic connectivity in structuring mangrove intra-specific diversity. Our findings are critical for mangroves' biogeography and evolution, but also for management strategies considering climate change and genetic biodiversity conservation.

Assessment of the Genetic Diversity and Population Structure of Rhizophora mucronata along Coastal Areas in Thailand

Unique and biodiverse, mangrove ecosystems provide humans with benefits and contribute to coastal protection. Rhizophora mucronata, a member of the Rhizophoraceae family, is prevalent in the mangrove forests of Thailand. R. mucronata's population structure and genetic diversity have received scant attention. Here, we sequenced the entire genome of R. mucronata using 10× Genomics technology and obtained an assembly size of 219 Mb with the N50 length of 542,540 bases. Using 2857 single nucleotide polymorphism (SNP) markers, this study investigated the genetic diversity and population structure of 80 R. mucronata accessions obtained from the mangrove forests in Thailand. The genetic diversity of R. mucronata was moderate (I = 0.573, Ho = 0.619, He = 0.391). Two subpopulations were observed and confirmed from both population structure and principal component analysis (PCA). Analysis of molecular variance (AMOVA) showed that there was more variation within populations than between them. Mean pairwise genetic differentiation (FST = 0.09) showed that there was not much genetic difference between populations. Intriguingly, the predominant clustering pattern in the R. mucronata population did not correspond to the Gulf of Thailand and the Andaman Sea, which are separated by the Malay Peninsula. Several factors could have influenced the R. mucronata genetic pattern, such as hybridization and anthropogenic factors. This research will provide important information for the future conservation and management of R. mucronata in Thailand.

Differential Adaptive Potential and Vulnerability to Climate-Driven Habitat Loss in Brazilian Mangroves

Geographic and environmental differences have been identified as factors influencing Brazilian mangrove trees' genetic diversity. Geographically, distinct species have convergent spatial genetic structures, indicating a limited gene flow between northern and southern populations. Environmentally, genomic studies and common garden experiments have found evidence of local adaptations along the latitudinal gradient of the Brazilian coast. However, little is known about how such adaptive heterogeneity could be affected by a rapidly changing climate in the coming decades, and the combination of deforestation and climate-induced habitat loss may affect these forests and their genetic diversity. Here, we applied two genomic-environmental association methods to model the turnover of potentially adaptive alleles for two dominant mangrove trees: Avicennia germinans and A. schaueriana. We analyzed a total of 134 individuals from six populations of A. germinans and 10 populations of A. schaueriana spanning the Brazilian coast from 1 °S to 28 °S. Gradient forest models identified temperature-related variables as the most important predictors for A. germinans outlier loci, whereas both temperature and precipitation were important for A. schaueriana. We modeled allele frequencies and projected them for future climatic scenarios to estimate adaptively driven vulnerability. We assessed climate-driven habitat loss through climate-only distribution models and calculated annual deforestation rates for each sampled region. Finally, to assess the vulnerability of individual populations, we combined the environmental suitability, deforestation data, and adaptive vulnerability projections. For both species, subtropical populations presented a higher vulnerability than equatorial populations to climate-driven habitat loss. We also identified deforestation rates at the sampled sites that were alarmingly higher than the global average mangrove deforestation rate. Our results provide improved estimates of the impacts of ongoing climate change and human-caused habitat loss on the distribution of mangroves and highlight the importance of site-based conservation strategies that consider individual subtropical and equatorial mangrove forests.

Reduction of Genetic Variation When Far From the Niche Centroid: Prediction for Mangrove Species

The niche-centroid hypothesis states that populations that are distributed near the centroid of the species' ecological niche will have higher fitness-related attributes, such as population abundance and genetic diversity than populations near the edges of the niche. Empirical evidence based on abundance and, more recently, genetic diversity data support this hypothesis. However, there are few studies that test this hypothesis in coastal species, such as mangroves. Here, we focused on the black mangrove Avicennia germinans. We combined ecological, heterozygosity, and allelic richness information from 1,419 individuals distributed in 40 populations with three main goals: (1) test the relationship between distance to the niche centroid and genetic diversity, (2) determine the set of environmental variables that best explain heterozygosity and allelic richness, and (3) predict the spatial variation in genetic diversity throughout most of the species' natural geographic range. We found a strong correlation between the distance to the niche centroid and both observed heterozygosity (Ho; ρ2 = 0.67 P < 0.05) and expected heterozygosity (He; ρ2 = 0.65, P < 0.05). The niche variables that best explained geographic variation in genetic diversity were soil type and precipitation seasonality. This suggests that these environmental variables influence mangrove growth and establishment, indirectly impacting standing genetic variation. We also predicted the spatial heterozygosity of A. germinans across its natural geographic range in the Americas using regression model coefficients. They showed significant power in predicting the observed data (R2 = 0.65 for Ho; R2 = 0.60 for He), even when we considered independent data sets (R2= 0.28 for Ho; R2 = 0.25 for He). Using this approach, several genetic diversity estimates can be implemented and may take advantage of population genomics to improve genetic diversity predictions. We conclude that the level of genetic diversity in A. germinans is in agreement with expectations of the niche-centroid hypothesis, namely that the highest heterozygosity and allelic richness (the basic genetic units for adaptation) are higher at locations of high environmental suitability. This shows that this approach is a potentially powerful tool in the conservation and management of this species, including for modelling changes in the face of climate change.

Population subdivision promoted by a sea-level-change-driven bottleneck: A glimpse from the evolutionary history of the mangrove plant Aegiceras corniculatum

Historic climate changes drive geographical populations of coastal plants to contract and recover dynamically, even die out completely. Species suffering from such bottlenecks usually lose intraspecific genetic diversity, but how do these events influence population subdivision patterns of coastal plants? Here, we investigated this question in the typical coastal plant: mangrove species Aegiceras corniculatum. Inhabiting the intertidal zone of the tropical and subtropical coast of the Indo-West Pacific oceans, its populations are deemed to be greatly shaped by historic sea-level fluctuations. Using dual methods of Sanger and Illumina sequencing, we found that the 18 sampled populations were structured into two groups, namely, the "Indo-Malayan" group, comprising three subgroups (the northern South China Sea, Gulf of Bengal, and Bali), and the "Pan-Australasia" group, comprising the subgroups of the southern South China Sea and Australasia. Based on the approximate Bayesian computations and Stairway Plot, we inferred that the southern South China Sea subgroup, which penetrates the interior of the "Indo-Malayan" group, originated from the Australasia subgroup, accompanied by a severe bottleneck event, with a spot of gene flow from both the Australasia and "Indo-Malayan" groups. Geographical barriers such as the Sundaland underlie the genetic break between Indian and Pacific Oceans, but the discontinuity between southern and northern South China Sea was originated from genetic drift in the bottleneck event. Hence, we revealed a case evidencing that the bottleneck event promoted population subdivision. This conclusion may be applicable in other taxa beyond coastal plants.

Genotypes of Rhizophora Propagules From a Non-mangrove Beach Provide Evidence of Recent Long-Distance Dispersal

Dispersal plays a crucial role in the connectivity of established mangrove populations and in species range dynamics. As species ranges shift in response to climate change, range expansions can occur from incremental short-distance dispersal events and from stochastic long-distance dispersal events. Most population genetic research dealt with historically accumulated events though evidence of actual propagule dispersal allows to estimate genotypic features and origin of founders. In this study, we aim to disentangle a contemporary dispersal event. Using microsatellite markers, we genotyped 60 Rhizophora racemosa drift propagules obtained on a bare unforested coastal area in southern Cameroon, estimated their relationship to 109 adult trees from most proximate sites (which were 3–85 km away), and assessed their relative difference with 873 trees of major mangrove areas (> 300 km) along the Cameroonian coastline. Proximate mangrove populations were considered as potential source populations in assignment tests. However, drift propagules could not be assigned to any of the Cameroonian mangrove sites and were genetically isolated from Cameroonian populations. Drift propagules showed higher levels of genetic diversity and private alleles giving a higher relatedness to each other than to any putative source population. Chloroplast sequences were used to confirm the identity of drift propagules as R. racemosa. We postulate that a complex interaction of ocean currents, estuarine geomorphology, and tidal patterns explain drift propagule dispersal to an area. Most likely the investigated cohort of propagules originated from more southern mangrove areas of the West African range beyond the Cameroonian border. This study unraveled the allelic, genetic, and genotypic features of stranded propagules following a stochastic long-distance dispersal. Transboundary dispersal of these propagules highlights the need for intergovernmental efforts in the management of biodiversity.

Complete Chloroplast Genome Variants Reveal Discrete Long-Distance Dispersal Routes of Rhizophora in the Western Indian Ocean

Historical processes of long-distance migration and ocean-wide expansion feature the global biogeographic pattern of Rhizophora species. Throughout the Indian Ocean, Rhizophora stylosa and Rhizophora mucronata seem to be a young phylogenetic group with an expansion of R. mucronata toward the Western Indian Ocean (WIO) driven by the South Equatorial Current (SEC). Nuclear microsatellites revealed genetic patterns and breaks; however, the estimation of propagule dispersal routes requires maternally inherited cytoplasmic markers. Here, we examine the phylogeography of 21 R. mucronata provenances across a >4,200 km coastal stretch in the WIO using R. stylosa as an outgroup. Full-length chloroplast genome (164,474 bp) and nuclear ribosomal RNA cistron (8,033 bp) sequences were assembled. The boundaries, junction point, sequence orientation, and stretch between LSC/IRb/SSC/IRa/LSC showed no differences with R. stylosa chloroplast genome. A total of 58 mutations in R. mucronata encompassing transitions/transversions, insertions-deletions, and mononucleotide repeats revealed three major haplogroups. Haplonetwork, Bayesian maximum likelihood (ML), and approximate Bayesian computation (ABC) analyses supported discrete historical migration events. An ancient haplogroup A in the Seychelles and eastern Madagascar was as different from other haplogroups as from R. stylosa. A star-like haplonetwork referred as the recent range expansion of haplogroup B from northern Madagascar toward the African mainland coastline, including a single variant spanning >1,800 km across the Mozambique Channel area (MCA). Populations in the south of Delagoa Bight contained haplogroup C and was originated from a unique bottleneck dispersal event. Divergence estimates of pre- and post-Last Glacial Maximum (LGM) illustrated the recent emergence of Rhizophora mangroves in the WIO compared to other oceans. Connectivity patterns could be aligned with the directionality of major ocean currents. Madagascar and the Seychelles each harbored haplogroups A and B, albeit among spatially separated populations, explained from a different migration era. Likewise, the Aldabra Atoll harbored spatially distinct haplotypes. Nuclear ribosomal cistron (8,033 bp) variants corresponded to haplogroups and confirmed admixtures in the Seychelles and Aldabra. These findings shed new light on the origins and dispersal routes of R. mucronata lineages that have shaped their contemporary populations in large regions of the WIO, which may be the important information for defining marine conservation units both at ocean scale and at the level of small islands.

Barrier to Gene Flow of Grey Mangrove Avicennia marina Populations in the Malay Peninsula as Revealed From Nuclear Microsatellites and Chloroplast Haplotypes

Contemporary mangrove forest areas took shape historically and their genetic connectivity depends on sea-faring propagules, subsequent settlement, and persistence in suitable environments. Mangrove species world-wide may experience genetic breaks caused by major land barriers or opposing ocean currents influencing their population genetic structure. For Malay Peninsula, several aquatic species showed strong genetic differentiation between East and West coast regions due to the Sunda shelf flooding since the Last Glacial Maximum. In this study genetic diversity and structure of Avicennia marina populations in Malay Peninsula were assessed using nuclear microsatellite markers and chloroplast sequences. Even though all populations showed identical morphological features of A. marina, three evolutionary significant units were obtained with nuclear and cytoplasmic markers. Avicennia marina along a 586 km stretch of the West coast differed strongly from populations along an 80 km stretch of the East coast featuring chloroplast capture of Avicennia alba in an introgressive A. marina. Over and above this expected East-West division, an intra-regional subdivision was detected among A. marina populations in the narrowest region of the Strait of Malacca. The latter genetic break was supported by an amova, structure, and barrier analysis whereas RST > FST indicated an evolutionary signal of long-lasting divergence. Two different haplotypes along the Western coast showed phylogeographic relationship with either a northern or a putative southern lineage, thereby assuming two Avicennia sources facing each other during Holocene occupation with prolonged separation in the Strait of Malacca. Migrate-n model testing supported a northward unidirectional stepping-stone migration route, although with an unclear directionality at the genetic break position, most likely due to weak oceanic currents. Low levels of genetic diversity and southward connectivity was detected for East coast Avicennia populations. We compared the fine-scale spatial genetic structure (FSGS) of Avicennia populations along the exposed coast in the East vs. the sheltered coast in the West. A majority of transects from both coastlines revealed no within-site kinship-based FSGS, although the remoteness of the open sea is important for Avicennia patches to maintain a neighborhood. The results provide new insights for mangrove researchers and managers for future in-depth ecological-genetic-based species conservation efforts in Malay Peninsula.

Expansion of the mangrove species Rhizophora mucronata in the Western Indian Ocean launched contrasting genetic patterns

Estimates of population structure and gene flow allow exploring the historical and contemporary processes that determine a species’ biogeographic pattern. In mangroves, large-scale genetic studies to estimate gene flow have been conducted predominantly in the Indo-Pacific and Atlantic region. Here we examine the genetic diversity and connectivity of Rhizophora mucronata across a > 3,000 km coastal stretch in the Western Indian Ocean (WIO) including WIO islands. Based on 359 trees from 13 populations and using 17 polymorphic microsatellite loci we detected genetic breaks between populations of the (1) East African coastline, (2) Mozambique Channel Area (3) granitic Seychelles, and (4) Aldabra and northern Madagascar. Genetic structure, diversity levels, and patterns of inferred connectivity, aligned with the directionality of major ocean currents, driven by bifurcation of the South Equatorial Current, northward into the East African Coastal Current and southward into the Mozambique Channel Area. A secondary genetic break between nearby populations in the Delagoa Bight coincided with high inbreeding levels and fixed loci. Results illustrate how oceanographic processes can connect and separate mangrove populations regardless of geographic distance.

Geographical and environmental contributions to genomic divergence in mangrove forests

Assessing the relative importance of geographical and environmental factors to the spatial distribution of genetic variation can provide information about the processes that maintain genetic variation in natural populations. With a globally wide but very restricted habitat distribution, mangrove trees are a useful model for studies aiming to understand the contributions of these factors. Mangroves occur along the continent–ocean interface of tropical and subtropical latitudes, regions considered inhospitable to many other types of plants. Here, we used landscape genomics approaches to investigate the relative contributions of geographical and environmental variables to the genetic variation of two black mangrove species, Avicennia schaueriana and Avicennia germinans, along the South American coast. Using single nucleotide polymorphisms, our results revealed an important role of ocean currents and geographical distance in the gene flow of A. schaueriana and an isolation-by-environment pattern in the organization of the genetic diversity of A. germinans. Additionally, for A. germinans, we observed significant correlations between genetic variation with evidence of selection and the influence of precipitation regimens, solar radiation and temperature patterns. These discoveries expand our knowledge about the evolution of mangrove trees and provide important information to predict future responses of coastal species to the expected global changes during this century.

Channel network structure determines genetic connectivity of landward-seaward Avicennia marina populations in a tropical bay

Mangrove ecosystems along the East African coast are often characterized by a disjunct zonation pattern of seaward and landward Avicennia marina trees. This disjunct zonation may be maintained through different positions in the tidal frame, yielding different dispersal settings. The spatial configuration of the landscape and coastal processes such as tides and waves is expected to largely influence the extent of propagule transport and subsequent regeneration. We hypothesized that landward sites would keep a stronger genetic structure over short distances in comparison with enhanced gene flow among regularly flooded seaward fringes. We tested this hypothesis from densely vegetated A. marina transects of a well-documented mangrove system (Gazi Bay, Kenya) and estimated local gene flow and kinship-based fine-scale genetic structure. Ten polymorphic microsatellite markers in 457 A. marina trees revealed no overall significant difference in levels of allele or gene diversities between sites that differ in hydrological proximity. Genetic structure and connectivity of A. marina populations however indicated an overall effect of geographic distance and revealed a pronounced distinction between channels and topographic setting. Migration models allowed to infer gene flow directionality among channels, and indicated a bidirectional steppingstone between seaward and nearest located landward stands. Admixed gene pools without any fine-scale structure were found within the wider and more exposed Kidogoweni channel, suggesting open systems. Elevated kinship values and structure over 5 to 20 m distance were only detected in two distant landward and seaward transects near the mouth of the Mkurumuji River, indicating local retention and establishment. Overall, our findings show that patterns of A. marina connectivity are explained by hydrological proximity, channel network structure, and hydrokinetic energy, rather than just their positioning as disjunct landward or seaward zones.

More than meets the eye: syntopic and morphologically similar mangrove killifish species show different mating systems and patterns of genetic structure along the Brazilian coast

Different mating systems can strongly affect the extent of genetic diversity and population structure among species. Given the increased effects of genetic drift on reduced population size, theory predicts that species undergoing self-fertilisation should have greater population structure than outcrossed species; however, demographic dynamics may affect this scenario. The mangrove killifish clade is composed of the two only known examples of self-fertilising species among vertebrates (Kryptolebias marmoratus and Kryptolebias hermaphroditus). A third species in this clade, Kryptolebias ocellatus, inhabits mangrove forests in southeast Brazil; however, its mating system and patterns of genetic structure have been rarely explored. Here, we examined the genetic structure and phylogeographic patterns of K. ocellatus along its distribution, using mitochondrial DNA and microsatellites to compare its patterns of genetic structure with the predominantly selfing and often-syntopic, K. hermaphroditus. Our results indicate that K. ocellatus reproduces mainly by outcrossing, with no current evidence of selfing, despite being an androdioecious species. Our results also reveal a stronger population subdivision in K. ocellatus compared to K. hermaphroditus, contrary to the theoretical predictions based on reproductive biology of the two species. Our findings indicate that, although morphologically similar, K. ocellatus and K. hermaphroditus had remarkably different evolutionary histories when colonising the same mangrove areas in southeastern Brazil, with other factors (e.g., time of colonisation, dispersal/establishment capacity) having more profound effects on the current population structuring of those species than differences in mating systems.

Genetic Structure and Connectivity of the Red Mangrove at Different Geographic Scales through a Complex Transverse Hydrological System from Freshwater to Marine Ecosystems

Mangrove forests are ecologically and economically valuable resources composed of trees morphologically and physiologically adapted to thrive across a range of habitats. Although, mangrove trees have high dispersion capacity, complexity of hydrological systems may lead to a fine-scale genetic structure (FSGS). The Transverse Coastal Corridor (TCC) is an interesting case of hydrological systems from fresh to marine waters where mangrove forests dominate. We evaluated genetic diversity and structure of Rhizophora mangle across a range of hydrological conditions within the TCC using inter-simple sequence repeat molecular markers. Sampling included four hydrological systems, two localities inside each system, and fringe and dwarf trees. Genetic differentiation was evaluated at local (<100 km) and fine (<10 km) scales through a set of analyses, and genetic diversity was evaluated at all scale levels and between fringe and dwarf physiognomic types. Rhizophora mangle exhibited a high genetic structure at both scales with high genetic diversity. The genetic structure observed among hydrological systems likely reflects the historical dispersion of mangroves, whereas the FSGS reflect contemporary processes such as seed dispersal restriction, habitat fragmentation, and local water flow regimes. A higher genetic diversity for dwarf than for fringe trees and differentiation between both physiognomic types at a fine-scale were observed and discussed.

Local adaptation of a dominant coastal tree to freshwater availability and solar radiation suggested by genomic and ecophysiological approaches

Local adaptation is often a product of environmental variations in geographical space and has implications for biodiversity conservation. We investigated the role of latitudinal heterogeneity in climate on the organization of genetic and phenotypic variation in the dominant coastal tree Avicennia schaueriana. In a common garden experiment, samples from an equatorial region, with pronounced seasonality in precipitation, accumulated less biomass, and showed lower stomatal conductance and transpiration, narrower xylem vessels, smaller leaves and higher reflectance of long wavelengths by the stem epidermis than samples from a subtropical region, with seasonality in temperature and no dry season. Transcriptomic differences identified between trees sampled under field conditions at equatorial and subtropical sites, were enriched in functional categories such as responses to temperature, solar radiation, water deficit, photosynthesis and cell wall biosynthesis. Remarkably, the diversity based on genome-wide SNPs revealed a north-south genetic structure and signatures of selection were identified for loci associated with photosynthesis, anthocyanin accumulation and the responses to osmotic and hypoxia stresses. Our results suggest the existence of divergence in key resource-use characteristics, likely driven by seasonality in water deficit and solar radiation. These findings provide a basis for conservation plans and for predicting coastal plants responses to climate change.

Molecular responses to freshwater limitation in the mangrove tree Avicennia germinans (Acanthaceae)

Environmental variation along the geographical space can shape populations by natural selection. In the context of global warming and changing precipitation regimes, it is crucial to understand the role of environmental heterogeneity in tropical trees adaptation, given their disproportional contribution to water and carbon biogeochemical cycles. Here, we investigated how heterogeneity in freshwater availability along tropical wetlands has influenced molecular variations of the black mangrove (Avicennia germinans). A total of 57 trees were sampled at seven sites differing markedly in precipitation regime and riverine freshwater inputs. Using 2,297 genome-wide single nucleotide polymorphic markers, we found signatures of natural selection by the association between variations in allele frequencies and environmental variables, including the precipitation of the warmest quarter and the annual precipitation. Additionally, we found candidate loci for selection based on statistical deviations from neutral expectations of interpopulation differentiation. Most candidate loci within transcribed sequences were functionally associated with central aspects of drought tolerance or plant response to drought. Moreover, our results suggest the occurrence of the rapid evolution of a population, probably in response to sudden and persistent limitations in plant access to soil water, following a road construction in 1974. Observations supporting rapid evolution included the reduction in tree size and changes in allele frequencies and in transcript expression associated with increased drought tolerance through the accumulation of osmoprotectants and antioxidants, biosynthesis of cuticles, protection against protein degradation, stomatal closure, photorespiration and photosynthesis. We describe a major role of spatial heterogeneity in freshwater availability in the specialization of this typically tropical tree.

Hybridization and polyploidy in the weeping lizard Liolaemus chiliensis (Squamata: Liolaemidae)

In reptiles, polyploidy is an unusual phenomenon that can originate from interspecific hybridization. In Chile, the lizard Liolaemus chiliensis provides a unique model with which to study the origin of polyploidy because it occurs in populations of diploid individuals and a few populations that also contain triploid and diploid–triploid mosaic lizards. To test whether L. chiliensis polyploids were hybrids between lineages within the species, we analysed the karyotype, the mitochondrial gene Cytb, seven microsatellite loci, and the linear and geometric morphometry of individuals throughout the species distribution. We found polyploidy at different localities throughout the distribution. Hybrids were detected with the microsatellite loci and morphological data. A clear relationship between hybridization and polyploidy could not be established because triploid and mosaic lizards were hybrids and purebreds. However, our results support the occurrence of both phenomena, and further research is needed to clarify how hybridization could be involved in the origin of polyploidy.

The presence of a cryptic barrier in the West Pacific Ocean suggests the effect of glacial climate changes on a widespread sea-dispersed plant, Vigna marina (Fabaceae)

Ocean currents are an important driver of evolution for sea-dispersed plants, enabling them to maintain reciprocal gene flow via sea-dispersed diaspores and obtain wide distribution ranges. Although geographic barriers are known to be the primary factors shaping present genetic structure of sea-dispersed plants, cryptic barriers which form clear genetic structure within oceanic regions are poorly understood. To test the presence of a cryptic barrier, we conducted a phylogeographic study together with past demographic inference for a widespread sea-dispersed plant, Vigna marina, using 308 individuals collected from the entire Indo-West Pacific (IWP) region. Chloroplast DNA variation showed strong genetic structure that separated populations into three groups: North Pacific (NP), South Pacific (SP) and Indian Ocean (IN) (F'CT among groups = 0.954-1.000). According to the Approximate Bayesian computation inference, splitting time between NP and SP was approximately 20,200 years (95%HPD, 4,530-95,400) before present. Moreover, a signal of recent population expansion was detected in the NP group. This study clearly showed the presence of a cryptic barrier in the West Pacific region of the distributional range of V. marina. The locations of the cryptic barrier observed in V. marina corresponded to the genetic breaks found in other plants, suggesting the presence of a common cryptic barrier for sea-dispersed plants. Demographic inference suggested that genetic structure related to this cryptic barrier has been present since the last glacial maximum and may reflect patterns of past population expansion from refugia.

The integration and application of genomic information in mangrove conservation

Article impact statement: Genomic information can help prevent species loss, facilitate reserve design and maximize adaptive potential of natural mangrove populations.

Geographical range size and latitude predict population genetic structure in a global survey

While genetic diversity within species is influenced by both geographical distance and environmental gradients, it is unclear what other factors are likely to promote population genetic structure. Using a machine learning framework and georeferenced DNA sequences from more than 8000 species, we demonstrate that geographical attributes of the species range, including total size, latitude and elevation, are the most important predictors of which species are likely to contain structured genetic variation. While latitude is well known as an important predictor of biodiversity, our work suggests that it also plays a key role in shaping diversity within species.

Short-distance barriers affect genetic variability of Rhizophora mangle L. in the Yucatan Peninsula

The environmental variability at local scale results in different physiognomic types of mangrove forest. However, this variability has never been considered in studies of mangrove genetic variability. This study analyzed the genetic and morphological variability and structure of Rhizophora mangle at regional and local scales in the Yucatan Peninsula. Thirteen mangrove populations (eight scrub and five tall), located in seven sites, were sampled, and their morphological variability and relationship with the availability of phosphorus and salinity were analyzed. The diversity and genetic structure were estimated at different hierarchical levels with nine microsatellites, also Bayesian inference and Principal Coordinates Analysis were used. We found a great morphological variability of R. mangle that responded to local environmental variability and not to the precipitation gradient of the peninsula. The genetic diversity found in the peninsula was greater than that reported for other populations in Mexico and was grouped into two regions: the Gulf of Mexico and the Caribbean Sea. At a local scale, tall and scrub mangroves had significant genetic differentiation suggesting that ecological barriers promote genetic differentiation within sites. These results need to be considered in future population genetic studies and for mangrove management and conservation.

Extremophiles as a Model of a Natural Ecosystem: Transcriptional Coordination of Genes Reveals Distinct Selective Responses of Plants Under Climate Change Scenarios

The goal of this research was to generate networks of co-expressed genes to explore the genomic responses of Rhizophora mangle L. populations to contrasting environments and to use gene network analysis to investigate their capacity for adaptation in the face of historical and future perturbations and climatic changes. RNA sequencing data were generated for R. mangle samples collected under field conditions from contrasting climate zones in the equatorial and subtropical regions of Brazil. A gene co-expression network was constructed using Pearson's correlation coefficient, showing correlations among 78,364 transcriptionally coordinated genes. Each region exhibited two distinct network profiles; genes correlated with the oxidative stress response showed higher relative expression levels in subtropical samples than in equatorial samples, whereas genes correlated with the hyperosmotic salinity response, heat response and UV response had higher expression levels in the equatorial samples than in the subtropical samples. In total, 992 clusters had enriched ontology terms, which suggests that R. mangle is under higher stress in the equatorial region than in the subtropical region. Increased heat may thus pose a substantial risk to species diversity at the center of its distribution range in the Americas. This study, which was performed using trees in natural field conditions, allowed us to associate the specific responses of genes previously described in controlled environments with their responses to the local habitat where the species occurs. The study reveals the effects of contrasting environments on gene expression in R. mangle, shedding light on the different abiotic variables that may contribute to the genetic divergence previously described for the species through the use of simple sequence repeats (SSRs). These effects may result from two fundamental processes in evolution, namely, phenotypic plasticity and natural selection.

Population genetic structure, introgression, and hybridization in the genus Rhizophora along the Brazilian coast

Mangrove plants comprise plants with similar ecological features that have enabled them to adapt to life between the sea and the land. Within a geographic region, different mangrove species share not only similar adaptations but also similar genetic structure patterns. Along the eastern coast of South America, there is a subdivision between the populations north and south of the continent's northeastern extremity. Here, we aimed to test for this north‐south genetic structure in Rhizophora mangle, a dominant mangrove plant in the Western Hemisphere. Additionally, we aimed to study the relationships between R. mangle, R. racemosa, and R. × harrisonii and to test for evidence of hybridization and introgression. Our results confirmed the north‐south genetic structure pattern in R. mangle and revealed a less abrupt genetic break in the northern population than those observed in Avicennia species, another dominant and widespread mangrove genus in the Western Hemisphere. These results are consistent with the role of oceanic currents influencing sea‐dispersed plants and differences between Avicennia and Rhizophora propagules in longevity and establishment time. We also observed that introgression and hybridization are relevant biological processes in the northeastern coast of South America and that they are likely asymmetric toward R. mangle, suggesting that adaptation might be a process maintaining this hybrid zone.

Contrasting Effects of Historical Sea Level Rise and Contemporary Ocean Currents on Regional Gene Flow of Rhizophora racemosa in Eastern Atlantic Mangroves

Mangroves are seafaring taxa through their hydrochorous propagules that have the potential to disperse over long distances. Therefore, investigating their patterns of gene flow provides insights on the processes involved in the spatial genetic structuring of populations. The coastline of Cameroon has a particular geomorphological history and coastal hydrology with complex contemporary patterns of ocean currents, which we hypothesize to have effects on the spatial configuration and composition of present-day mangroves within its spans. A total of 982 trees were sampled from 33 transects (11 sites) in 4 estuaries. Using 11 polymorphic SSR markers, we investigated genetic diversity and structure of Rhizophora racemosa, a widespread species in the region. Genetic diversity was low to moderate and genetic differentiation between nearly all population pairs was significant. Bayesian clustering analysis, PCoA, estimates of contemporary migration rates and identification of barriers to gene flow were used and complemented with estimated dispersal trajectories of hourly released virtual propagules, using high-resolution surface current from a mesoscale and tide-resolving ocean simulation. These indicate that the Cameroon Volcanic Line (CVL) is not a present-day barrier to gene flow. Rather, the Inter-Bioko-Cameroon (IBC) corridor, formed due to sea level rise, allows for connectivity between two mangrove areas that were isolated during glacial times by the CVL. Genetic data and numerical ocean simulations indicated that an oceanic convergence zone near the Cameroon Estuary complex (CEC) presents a strong barrier to gene flow, resulting in genetic discontinuities between the mangrove areas on either side. This convergence did not result in higher genetic diversity at the CEC as we had hypothesized. In conclusion, the genetic structure of Rhizophora racemosa is maintained by the contrasting effects of the contemporary oceanic convergence and historical climate change-induced sea level rise.

Species distribution and introgressive hybridization of two Avicennia species from the Western Hemisphere unveiled by phylogeographic patterns

Background

Mangrove plants grow in the intertidal zone in tropical and subtropical regions worldwide. The global latitudinal distribution of the mangrove is mainly influenced by climatic and oceanographic features. Because of current climate changes, poleward range expansions have been reported for the major biogeographic regions of mangrove forests in the Western and Eastern Hemispheres. There is evidence that mangrove forests also responded similarly after the last glaciation by expanding their ranges. In this context, the use of genetic tools is an informative approach for understanding how historical processes and factors impact the distribution of mangrove species. We investigated the phylogeographic patterns of two Avicennia species, A. germinans and A. schaueriana, from the Western Hemisphere using nuclear and chloroplast DNA markers.

Results

Our results indicate that, although Avicennia bicolor, A. germinans and A. schaueriana are independent lineages, hybridization between A. schaueriana and A. germinans is a relevant evolutionary process. Our findings also reinforce the role of long-distance dispersal in widespread mangrove species such as A. germinans, for which we observed signs of transatlantic dispersal, a process that has, most likely, contributed to the breadth of the distribution of A. germinans. However, along the southern coast of South America, A. schaueriana is the only representative of the genus. The distribution patterns of A. germinans and A. schaueriana are explained by their different responses to past climate changes and by the unequal historical effectiveness of relative gene flow by propagules and pollen.

Conclusions

We observed that A. bicolor, A. germinans and A. schaueriana are three evolutionary lineages that present historical and ongoing hybridization on the American continent. We also inferred a new evidence of transatlantic dispersal for A. germinans, which may have contributed to its widespread distribution. Despite the generally wider distribution of A. germinans, only A. schaueriana is found in southern South America, which may be explained by the different demographic histories of these two species and the larger proportion of gene flow produced by propagules rather than pollen in A. schaueriana. These results highlight that these species responded in different ways to past events, indicating that such differences may also occur in the currently changing world.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0343-z) contains supplementary material, which is available to authorized users.