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

Seabird nutrient subsidies enrich mangrove ecosystems and are exported to nearby coastal habitats

Eutrophication by human-derived nutrient enrichment is a major threat to mangroves, impacting productivity, ecological functions, resilience, and ecosystem services. Natural mangrove nutrient enrichment processes, however, remain largely uninvestigated. Mobile consumers such as seabirds are important vectors of cross-ecosystem nutrient subsidies to islands but how they influence mangrove ecosystems is poorly known. We assessed the contribution, uptake, cycling, and transfer of nutrients from seabird colonies in remote mangrove systems free of human stressors. We found that nutrients from seabird guano enrich mangrove plants, reduce nutrient limitations, enhance mangrove invertebrate food webs, and are exported to nearby coastal habitats through tidal flow. We show that seabird nutrient subsidies in mangroves can be substantial, improving the nutrient status and health of mangroves and adjacent coastal habitats. Conserving mobile consumers, such as seabirds, is therefore vital to preserve and enhance their role in mangrove productivity, resilience, and provision of diverse functions and services.

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.

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.

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.

Inheritance of DNA methylation differences in the mangrove Rhizophora mangle

The capacity to respond to environmental challenges ultimately relies on phenotypic variation which manifests from complex interactions of genetic and nongenetic mechanisms through development. While we know something about genetic variation and structure of many species of conservation importance, we know very little about the nongenetic contributions to variation. Rhizophora mangle is a foundation species that occurs in coastal estuarine habitats throughout the neotropics where it provides critical ecosystem functions and is potentially threatened by anthropogenic environmental changes. Several studies have documented landscape-level patterns of genetic variation in this species, but we know virtually nothing about the inheritance of nongenetic variation. To assess one type of nongenetic variation, we examined the patterns of DNA sequence and DNA methylation in maternal plants and offspring from natural populations of R. mangle from the Gulf Coast of Florida. We used a reduced representation bisulfite sequencing approach (epi-genotyping by sequencing; epiGBS) to address the following questions: (a) What are the levels of genetic and epigenetic diversity in natural populations of R. mangle? (b) How are genetic and epigenetic variation structured within and among populations? (c) How faithfully is epigenetic variation inherited? We found low genetic diversity but high epigenetic diversity from natural populations of maternal plants in the field. In addition, a large portion (up to ~25%) of epigenetic differences among offspring grown in common garden was explained by maternal family. Therefore, epigenetic variation could be an important source of response to challenging environments in the genetically depauperate populations of this foundation species.