The integration and application of genomic information in mangrove conservation

Untapped rich microbiota of mangroves of Pakistan: diversity and community compositions

The mangrove ecosystem is the world's fourth most productive ecosystem in terms of service value and offering rich biological resources. Microorganisms play vital roles in these ecological processes, thus researching the mangroves-microbiota is crucial for a deeper comprehension of mangroves dynamics. Amplicon sequencing that targeted V4 region of 16S rRNA gene was employed to profile the microbial diversities and community compositions of 19 soil samples, which were collected from the rhizosphere of 3 plant species (i.e., Avicennia marina, Ceriops tagal, and Rhizophora mucronata) in the mangrove forests of Lasbela coast, Pakistan. A total of 67 bacterial phyla were observed from three mangroves species, and these taxa were classified into 188 classes, 453 orders, 759 families, and 1327 genera. We found that Proteobacteria (34.9-38.4%) and Desulfobacteria (7.6-10.0%) were the dominant phyla followed by Chloroflexi (6.6-7.3%), Gemmatimonadota (5.4-6.8%), Bacteroidota (4.3-5.5%), Planctomycetota (4.4-4.9%) and Acidobacteriota (2.7-3.4%), Actinobacteriota (2.5-3.3%), and Crenarchaeota (2.5-3.3%). After considering the distribution of taxonomic groups, we prescribe that the distinctions in bacterial community composition and diversity are ascribed to the changes in physicochemical attributes of the soil samples (i.e., electrical conductivity (ECe), pH, total organic matter (OM), total organic carbon (OC), available phosphorus (P), and extractable potassium (CaCO3). The findings of this study indicated a high-level species diversity in Pakistani mangroves. The outcomes may also aid in the development of effective conservation policies for mangrove ecosystems, which have been hotspots for anthropogenic impacts in Pakistan. To our knowledge, this is the first microbial research from a Pakistani mangrove forest.

The core mangrove microbiome reveals shared taxa potentially involved in nutrient cycling and promoting host survival

BACKGROUND

Microbes have fundamental roles underpinning the functioning of our planet, they are involved in global carbon and nutrient cycling, and support the existence of multicellular life. The mangrove ecosystem is nutrient limited and if not for microbial cycling of nutrients, life in this harsh environment would likely not exist. The mangroves of Southeast Asia are the oldest and most biodiverse on the planet, and serve vital roles helping to prevent shoreline erosion, act as nursery grounds for many marine species and sequester carbon. Despite these recognised benefits and the importance of microbes in these ecosystems, studies examining the mangrove microbiome in Southeast Asia are scarce.cxs RESULTS: Here we examine the microbiome of Avicenia alba and Sonneratia alba and identify a core microbiome of 81 taxa. A further eight taxa (Pleurocapsa, Tunicatimonas, Halomonas, Marinomonas, Rubrivirga, Altererythrobacte, Lewinella, and Erythrobacter) were found to be significantly enriched in mangrove tree compartments suggesting key roles in this microbiome. The majority of those identified are involved in nutrient cycling or have roles in the production of compounds that promote host survival.

CONCLUSION

The identification of a core microbiome furthers our understanding of mangrove microbial biodiversity, particularly in Southeast Asia where studies such as this are rare. The identification of significantly different microbial communities between sampling sites suggests environmental filtering is occurring, with hosts selecting for a microbial consortia most suitable for survival in their immediate environment. As climate change advances, many of these microbial communities are predicted to change, however, without knowing what is currently there, it is impossible to determine the magnitude of any deviations. This work provides an important baseline against which change in microbial community can be measured.

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.

Comparative Analysis and Phylogenetic Relationships of Ceriops Species (Rhizophoraceae) and Avicennia lanata (Acanthaceae): Insight into the Chloroplast Genome Evolution between Middle and Seaward Zones of Mangrove Forests

Ceriops and Avicennia are true mangroves in the middle and seaward zones of mangrove forests, respectively. The chloroplast genomes of Ceriops decandra, Ceriops zippeliana, and Ceriops tagal were assembled into lengths of 166,650, 166,083 and 164,432 bp, respectively, whereas Avicennia lanata was 148,264 bp in length. The gene content and gene order are highly conserved among these species. The chloroplast genome contains 125 genes in A. lanata and 129 genes in Ceriops species. Three duplicate genes (rpl2, rpl23, and trnM-CAU) were found in the IR regions of the three Ceriops species, resulting in expansion of the IR regions. The rpl32 gene was lost in C. zippeliana, whereas the infA gene was present in A. lanata. Short repeats (<40 bp) and a lower number of SSRs were found in A. lanata but not in Ceriops species. The phylogenetic analysis supports that all Ceriops species are clustered in Rhizophoraceae and A. lanata is in Acanthaceae. In a search for genes under selective pressures of coastal environments, the rps7 gene was under positive selection compared with non-mangrove species. Finally, two specific primer sets were developed for species identification of the three Ceriops species. Thus, this finding provides insightful genetic information for evolutionary relationships and molecular markers in Ceriops and Avicennia species.

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 &lt; 0.05) and expected heterozygosity (He; ρ2 = 0.65, P &lt; 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.

Phylogeographic pattern of a cryptoviviparous mangrove, Aegiceras corniculatum, in the Indo-West Pacific, provides insights for conservation actions

This study identified the historical geoclimatic factors which caused low genetic diversity and strong phylogeographic structure in a cryptoviviparous mangrove. The phylogeographic pattern was used to suggest conservation actions. Phylogeographic studies are used to understand the spatial distribution and evolution of genetic diversity, and have major conservation implications, especially for threatened taxa like the mangroves. This study aimed to assess the phylogeographic pattern of Aegiceras corniculatum, a cryptoviviparous mangrove, across its distribution range in the Indo-West Pacific (IWP) region. We genotyped 398 samples, collected from 37 populations, at four chloroplast DNA (cpDNA) loci, and identified the influence of historical processes on the contemporary population structure of the species. Low genetic diversity at the population level was observed. The evolutionary relationship between 12 cpDNA haplotypes suggested a strong phylogeographic structure, which was further validated by the clustering algorithms and proportioning of maximum variation among hierarchical population groups. The magnitude and direction of historical gene flow indicated that the species attained its wide distribution from its likely ancestral area of the Malay Archipelago. The divergence time estimates of the haplotypes indicated that the geoclimatic changes during the Pleistocene, especially the glacial sea-level changes and emergence of landmasses, hindered genetic exchange and created genetic differentiation between the phylogenetic groups. The species overwintered the last glacial maxima in multiple refugia in the IWP, as identified by the environmental niche modelling. Overall, our findings indicated that ancient glacial vicariance had influenced the present genetic composition of A. corniculatum, which was maintained by the current demographic features of this region. We discussed how these findings can be used to prioritize areas for conservation actions, restore disturbed habitats and prevent further genetic erosion.

Global phylogeography of a pantropical mangrove genus Rhizophora

Rhizophora is a key genus for revealing the formation process of the pantropical distribution of mangroves. In this study, in order to fully understand the historical scenario of Rhizophora that achieved pantropical distribution, we conducted phylogeographic analyses based on nucleotide sequences of chloroplast and nuclear DNA as well as microsatellites for samples collected worldwide. Phylogenetic trees suggested the monophyly of each AEP and IWP lineages respectively except for R. samoensis and R. × selala. The divergence time between the two lineages was 10.6 million years ago on a dated phylogeny, and biogeographic stochastic mapping analyses supported these lineages separated following a vicariant event. These data suggested that the closure of the Tethys Seaway and the reduction in mangrove distribution followed by Mid-Miocene cooling were key factors that caused the linage diversification. Phylogeographic analyses also suggested the formation of the distinctive genetic structure at the AEP region across the American continents around Pliocene. Furthermore, long-distance trans-pacific dispersal occurred from the Pacific coast of American continents to the South Pacific and formed F1 hybrid, resulting in gene exchange between the IWP and AEP lineages after 11 million years of isolation. Considering the phylogeny and phylogeography with divergence time, a comprehensive picture of the historical scenario behind the pantropical distribution of Rhizophora is updated.

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.

Distinct fungal communities associated with different organs of the mangrove Sonneratia alba in the Malay Peninsula

Mangrove forests are key tropical marine ecosystems that are rich in fungi, but our understanding of fungal communities associated with mangrove trees and their various organs remains limited because much of the diversity lies within the microbiome. In this study, we investigated the fungal communities associated with the mangrove tree Sonneratia alba throughout Peninsular Malaysia and Singapore. At each sampling location, we collected leaves, fruits, pneumatophores and sediment samples and performed amplicon sequencing of the ribosomal internal transcribed spacer 1 to characterise the associated communities. Results show distinct fungal communities at each sampled location with further differentiation according to the plant part. We find a significant distance decay of similarity, particularly for sediment samples due to the greater variability of sediment environments relative to the more stable fungal habitats provided by living plant organs. We are able to assign taxonomy to the majority of sequences from leaves and fruits, but a much larger portion of the sequences recovered from pneumatophores and sediment samples could not be identified. This pattern underscores the limited mycological research performed in marine environments and demonstrates the need for a concerted research effort on multiple species to fully characterise the coastal microbiome and its role in the functioning of marine ecosystems.

Land masses and oceanic currents drive population structure of Heritiera littoralis, a widespread mangrove in the Indo-West Pacific

Phylogeographic forces driving evolution of sea-dispersed plants are often influenced by regional and species characteristics, although not yet deciphered at a large spatial scale for many taxa like the mangrove species Heritiera littoralis. This study aimed to assess geographic distribution of genetic variation of this widespread mangrove in the Indo-West Pacific region and identify the phylogeographic factors influencing its present-day distribution. Analysis of five chloroplast DNA fragments' sequences from 37 populations revealed low genetic diversity at the population level and strong genetic structure of H. littoralis in this region. The estimated divergence times between the major genetic lineages indicated that glacial level changes during the Pleistocene epoch induced strong genetic differentiation across the Indian and Pacific Oceans. In comparison to the strong genetic break imposed by the Sunda Shelf toward splitting the lineages of the Indian and Pacific Oceans, the genetic differentiation between Indo-Malesia and Australasia was not so prominent. Long-distance dispersal ability of H. littoralis propagules helped the species to attain transoceanic distribution not only across South East Asia and Australia, but also across the Indian Ocean to East Africa. However, oceanic circulation pattern in the South China Sea was found to act as a barrier creating further intraoceanic genetic differentiation. Overall, phylogeographic analysis in this study revealed that glacial vicariance had profound influence on population differentiation in H. littoralis and caused low genetic diversity except for the refugia populations near the equator which might have persisted through glacial maxima. With increasing loss of suitable habitats due to anthropogenic activities, these findings therefore emphasize the urgent need for conservation actions for all populations throughout the distribution range of H. littoralis.

Multiplex microsatellite PCR panels for the neotropical red mangrove, Rhizophora mangle: combining efforts towards a cost-effective and modifiable tool to better inform conservation and management

Better-informed mangrove conservation and management practices are needed as the ecosystem services provided by these intertidal forests continue to be threatened by increasing anthropogenic pressures and climate change. Multiple layers of knowledge are required to achieve this goal, including insights into population genetics of mangrove species. Understanding the importance of population-genetic insights to conservation, multiple research groups have developed microsatellite loci for the widespread, neotropical red mangrove, Rhizophora mangle. However, although a wealth of genetic markers exist, empirical research is limited in the number of these loci employed. Here, we designed two multiplex PCR panels that combine seven novel loci developed for this work and eight previously-developed loci from three research groups to generate 15-locus genotypes, more than twice the average number of loci used in previous research, in only two PCR. We demonstrated utility in R. mangle from four sites across ~ 2500 km near this species’ northern latitudinal limits, and that these multiplex panels were better able to delineate populations than data subsets with numbers of loci comparable to previous research. We focus our discussion on how this tool is a more-informative, efficient (both in terms of time and resources), and easily-modifiable alternative to address many pressing conservation and management issues, such as the generation of baseline genetic data for areas not yet studied, better defining management units, and monitoring genetic effects of restoration projects. We also provide a quick protocol that outlines each step in this procedure to facilitate the use of this tool by others.

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.

Random PCR-based genotyping by sequencing technology GRAS-Di (genotyping by random amplicon sequencing, direct) reveals genetic structure of mangrove fishes

While various technologies for high-throughput genotyping have been developed for ecological studies, simple methods tolerant to low-quality DNA samples are still limited. In this study, we tested the availability of a random PCR-based genotyping-by-sequencing technology, genotyping by random amplicon sequencing, direct (GRAS-Di). We focused on population genetic analysis of estuarine mangrove fishes, including two resident species, the Amboina cardinalfish (Fibramia amboinensis, Bleeker, 1853) and the Duncker's river garfish (Zenarchopterus dunckeri, Mohr, 1926), and a marine migrant, the blacktail snapper (Lutjanus fulvus, Forster, 1801). Collections were from the Ryukyu Islands, southern Japan. PCR amplicons derived from ~130 individuals were pooled and sequenced in a single lane on a HiSeq2500 platform, and an average of three million reads was obtained per individual. Consensus contigs were assembled for each species and used for genotyping of single nucleotide polymorphisms by mapping trimmed reads onto the contigs. After quality filtering steps, 4,000-9,000 putative single nucleotide polymorphisms were detected for each species. Although DNA fragmentation can diminish genotyping performance when analysed on next-generation sequencing technology, the effect was small. Genetic differentiation and a clear pattern of isolation-by-distance was observed in F. amboinensis and Z. dunckeri by means of principal component analysis, F_ST and the admixture analysis. By contrast, L. fulvus comprised a genetically homogeneous population with directional recent gene flow. These genetic differentiation patterns reflect patterns of estuary use through life history. These results showed the power of GRAS-Di for fine-grained genetic analysis using field samples, including mangrove fishes.