Pronounced genetic differentiation and recent secondary contact in the mangrove tree Lumnitzera racemosa revealed by population genomic analyses

Glacial vicariance and oceanic circulation shape population structure of the coastal legume Derris trifoliata in the Indo-West Pacific

PREMISE

The phylogeography of coastal plant species is shaped by contemporary and historical biogeographic processes. In this study, we aim to decipher the phylogeography of Derris trifoliata, a woody legume of relatively recent origin and wide distribution, in coastal areas in the Indo-West Pacific (IWP) region.

METHODS

Genetic diversity and population structure were assessed by analyzing six nuclear and three chloroplast DNA sequences from 30 populations across the species' range. Phylogeography was inferred by estimating gene flow, divergence time, historical population size changes, and historical habitat suitability using paleoclimatic niche modeling.

RESULTS

High genetic diversity was observed at the species level. The populations of three oceanic regions included in this study (i.e., Indian Ocean, South China Sea, and Pacific Ocean) formed distinct clades and likely diverged during the late Pleistocene. Potential barriers to gene flow were identified, including the Sunda and Sahul shelves, geographic distance, and current patterns of oceanic circulation. Analysis of changes in population size supported the bottleneck model, which was strengthened by estimates of habitat suitability across paleoclimatic conditions.

CONCLUSIONS

The once widespread distribution of D. trifoliata was fragmented by changes in climatic suitability and biogeographic barriers that arose following sea-level changes during the Pleistocene. In addition, contemporary patterns of oceanic circulation and geographic distance between populations appear to maintain genetic differentiation across its distribution in the IWP.

Establishing community-wide DNA barcode references for conserving mangrove forests in China

BACKGROUND

Mangrove ecosystems have been the focus of global attention for their crucial role in sheltering coastal communities and retarding global climate change by sequestering 'blue carbon'. China is relatively rich in mangrove diversity, with one-third of the ca. 70 true mangrove species and a number of mangrove associate species occurring naturally along the country's coasts. Mangrove ecosystems, however, are widely threatened by intensifying human disturbances and rising sea levels. DNA barcoding technology may help protect mangrove ecosystems by providing rapid species identification.

RESULTS

To investigate this potential, 898 plant specimens were collected from 33 major mangrove sites in China. Based on the morphologic diagnosis, the specimens were assigned to 72 species, including all 28 true mangrove species and all 12 mangrove associate species recorded in China. Three chloroplast DNA markers rbcL, trnH-psbA, matK, and one nuclear marker ITS2 were chosen to investigate the utility of using barcoding to identify these species. According to the criteria of barcoding gaps in genetic distance, sequence similarity, and phylogenetic monophyly, we propose that a single marker, ITS2, is sufficient to barcode the species of mangroves and their associates in China. Furthermore, rbcL or trnH-psbA can also be used to gather supplement confirming data. In using these barcodes, we revealed a very low level of genetic variation among geographic locations in the mangrove species, which is an alert to their vulnerability to climate and anthropogenic disturbances.

CONCLUSION

We suggest using ITS2 to barcode mangrove species and terrestrial coastal plants in South China. The DNA barcode sequences we obtained would be valuable in monitoring biodiversity and the restoration of ecosystems, which are essential for mangrove conservation.

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.

Liao L, Papenbrock J. Analysis of rDNA reveals a high genetic diversity of Halophila major in the Wallacea region

The genus Halophila shows the highest species diversity within the seagrass genera. Southeast Asian countries where several boundary lines exist were considered as the origin of seagrasses. We hypothesize that the boundary lines, such as Wallace's and Lydekker's Lines, may act as marine geographic barriers to the population structure of Halophila major. Seagrass samples were collected at three islands in Vietnamese waters and analyzed by the molecular maker ITS. These sequences were compared with published ITS sequences from seagrasses collected in the whole region of interest. In this study, we reveal the haplotype and nucleotide diversity, linking population genetics, phylogeography, phylogenetics and estimation of relative divergence times of H. major and other members of the Halophila genus. The morphological characters show variation. The results of the ITS marker analysis reveal smaller groups of H. major from Myanmar, Shoalwater Bay (Australia) and Okinawa (Japan) with high supporting values. The remaining groups including Sri Lanka, Viet Nam, the Philippines, Thailand, Malaysia, Indonesia, Two Peoples Bay (Australia) and Tokushima (Japan) showed low supporting values. The Wallacea region shows the highest haplotype and also nucleotide diversity. Non-significant differences were found among regions, but significant differences were presented among populations. The relative divergence times between some members of section Halophila were estimated 2.15-6.64 Mya.

Contrasting Phylogeographic Patterns in Lumnitzera Mangroves Across the Indo-West Pacific

Mangroves are ecologically important forest communities in tropical and subtropical coasts, the effective management of which requires understanding of their phylogeographic patterns. However, these patterns often vary among different species, even among ecologically similar taxa or congeneric species. Here, we investigated the levels and patterns of genetic variation within Lumnitzera consisting of two species (L. racemosa and L. littorea) with nearly sympatric ranges across the Indo-West Pacific (IWP) region by sequencing three chloroplast DNA regions (for both species) and genotyping 11 nuclear microsatellite loci (for L. littorea). Consistent with findings in studies on other mangrove species, we found that both L. racemosa and L. littorea showed relatively high genetic variation among populations but low genetic variation within populations. Haplotype network and genetic clustering analyses indicated two well-differentiated clades in both L. racemosa and L. littorea. The relationship between geographic and genetic distances and divergence time estimates of the haplotypes indicated that limited dispersal ability of the propagules, emergence of land barriers during ancient sea-level changes, and contemporary oceanic circulation pattern in the IWP influenced the current population structure of the two species. However, the position of genetic break was found to vary between the two species: in L. racemosa, strong divergence was observed between populations from the Indian Ocean and the Pacific Ocean possibly due to land barrier effect of the Malay Peninsula; in L. littorea, the phylogeographic pattern was created by a more eastward genetic break along the biogeographic barrier identified as the Huxley's line. Overall, our findings strongly supported previous hypothesis of mangrove species divergence and revealed that the two Lumnitzera species have different phylogeographic patterns despite their close genetic relationship and similar current geographic distribution. The findings also provided references for the management of Lumnitzera mangroves, especially for the threatened L. littorea.

Low genetic diversity indicating the threatened status of Rhizophora apiculata (Rhizophoraceae) in Malaysia: declined evolution meets habitat destruction

Worldwide, many mangrove species are experiencing significant population declines, including Rhizophora apiculata, which is one of the most widespread and economically important species in tropical Asia. In Malaysia, there has been an alarming decline in R. apiculata populations driven primarily by anthropogenic activities. However, the lack of genetic and demographic information on this species has hampered local efforts to conserve it. To address these gaps, we generated novel genetic information for R. apiculata, based on 1,120 samples collected from 39 natural populations in Peninsular Malaysia. We investigated its genetic diversity and genetic structure with 19 transcriptome and three nuclear microsatellite markers. Our analyses revealed a low genetic diversity (mean H_e: 0.352) with significant genetic differentiation (F_ST: 0.315) among populations of R. apiculata. Approximately two-third of the populations showed significant excess of homozygotes, indicating persistent inbreeding which might be due to the decrease in population size or fragmentation. From the cluster analyses, the populations investigated were divided into two distinct clusters, comprising the west and east coasts of Peninsular Malaysia. The western cluster was further divided into two sub-clusters with one of the sub-clusters showing strong admixture pattern that harbours high levels of genetic diversity, thus deserving high priority for conservation.

Chloroplast Genome Sequences and Comparative Analyses of Combretaceae Mangroves with Related Species

In the Combretaceae family, only two species of Lumnitzera and one species of Laguncularia belong to mangroves. Among them, Lumnitzera littorea (Jack) Voigt. is an endangered mangrove plant in China for the limited occurrence and seed abortion. In contrast, Lumnitzera racemosa Willd. is known as the most widespread mangrove plant in China. Laguncularia racemosa C. F. Gaertn., an exotic mangrove in China, has the fast growth and high adaptation ability. To better understand the phylogenetic positions of these mangroves in Combretaceae and in Myrtales and to provide information for studies on evolutionary adaptation for intertidal habitat, the complete chloroplast (cp) genomes of Lu. racemosa and La. racemosa were sequenced. Furthermore, we present here the results from the assembly and annotation of the two cp genomes, which were further subjected to the comparative analysis with Lu. littorea cp genomes we published before and other eleven closely related species within Myrtales. The chloroplast genomes of the three Combretaceae mangrove species: Lu. littorea, Lu. racemosa, and La. racemosa are 159,687 bp, 159,473 bp, and 158,311 bp in size. All three cp genomes host 130 genes including 85 protein-coding genes, 37 tRNAs, and 4 rRNAs. A comparative analysis of those three genomes revealed the high similarity of genes in coding-regions and conserved gene order in the IR and LSC/SSC regions. The differences between Lumnitzera and Laguncularia cp genomes are the locations of rps19 and rpl2 genes in the IR/SC boundary regions. Investigating the effects of selection events on shared protein-coding genes showed a relaxed selection had acted on the ycf2, ycf1, and matK genes of Combretaceae mangroves compared to the nonmangrove species Eucalyptus aromaphloia. The phylogenetic analysis based on the whole chloroplast genome sequence with one outgroup species strongly supported three Combretaceae mangroves together with other two Combretaceae species formed a cluster in Combretaceae. This study is the first report on the comparative analysis of three Combretaceae mangrove chloroplast genomes, which will provide the significant information for understanding photosynthesis and evolution in Combretaceae mangrove plants.

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.

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.

Comparative transcriptional profiling of tildipirosin-resistant and sensitive Haemophilus parasuis

Numerous studies have been conducted to examine the molecular mechanism of Haemophilus parasuis resistance to antibiotic, but rarely to tildipirosin. In the current study, transcriptional profiling was applied to analyse the variation in gene expression of JS0135 and tildipirosin-resistant JS32. The growth curves showed that JS32 had a higher growth rate but fewer bacteria than JS0135. The cell membranes of JS32 and a resistant clinical isolate (HB32) were observed to be smoother than those of JS0135. From the comparative gene expression profile 349 up- and 113 downregulated genes were observed, covering 37 GO and 63 KEGG pathways which are involved in biological processes (11), cellular components (17), molecular function (9), cellular processes (1), environmental information processing (4), genetic information processing (9) and metabolism (49) affected in JS32. In addition, the relative overexpression of genes of the metabolism pathway (HAPS_RS09315, HAPS_RS09320), ribosomes (HAPS_RS07815) and ABC transporters (HAPS_RS10945) was detected, particularly the metabolism pathway, and verified with RT-qPCR. Collectively, the gene expression profile in connection with tildipirosin resistance factors revealed unique and highly resistant determinants of H. parasuis to macrolides that warrant further attention due to the significant threat of bacterial resistance.

Comparative Analysis of the Pattern of Population Genetic Diversity in Three Indo-West Pacific Rhizophora Mangrove Species

Rhizophora species are the most widely distributed mangrove trees in the Indo-West Pacific (IWP) region. Comparative studies of these species with shared life history traits can help identify evolutionary factors that have played most important roles in determining genetic diversity within and between populations in ocean-current dispersed mangrove tree species. We sampled 935 individuals from 54 natural populations for genotyping with 13 microsatellite markers to investigate the level of genetic variation, population structure, and gene flow on a broad geographic scale in Rhizophora apiculata, Rhizophora mucronata, and Rhizophora stylosa across the IWP region. In contrast to the pattern expected of long-lived woody plants with predominant wind-pollination, water-dispersed seeds and wide geographic range, genetic variation within populations was generally low in all the three species, especially in those peripheral populations from geographic range limits. Although the large water-buoyant propagules of Rhizophora have capacity for long distance dispersal, such events might be rare in reality, as reflected by the low level of gene flow and high genetic differentiation between most of population pairs within each species. Phylogeographic separation of Australian and Pacific island populations from SE Asian lineages previously revealed with DNA sequence data was still detectable in R. apiculata based on genetic distances, but this pattern of disjunction was not always evident in R. mucronata and R. stylosa, suggesting that fast-evolving molecular markers could be more suitable for detecting contemporary genetic structure but not deep evolutionary divergence caused by historical vicariance. Given that mangrove species generally have small effective population sizes, we conclude that genetic drift coupled with limited gene flow have played a dominant role in producing the current pattern of population genetic diversity in the IWP Rhizophora species, overshadowing the effects of their life history traits. Recent population fragmentation and disturbances arising from human activities could further endanger genetic diversity in mangrove trees.

Ancient Geographical Barriers Drive Differentiation among Sonneratia caseolaris Populations and Recent Divergence from S. lanceolata

Glacial vicariance is thought to influence population dynamics and speciation of many marine organisms. Mangroves, a plant group inhabiting intertidal zones, were also profoundly influenced by Pleistocene glaciations. In this study, we investigated phylogeographic patterns of a widespread mangrove species Sonneratia caseolaris and a narrowly distributed, closely related species S. lanceolata to infer their divergence histories and related it to historical geological events. We sequenced two chloroplast fragments and five nuclear genes for one population of S. lanceolata and 12 populations of S. caseolaris across the Indo-West Pacific (IWP) region to evaluate genetic differentiation and divergence time among them. Phylogenetic analysis based on sequences of nuclear ribosomal internal transcribed spacer and a nuclear gene rpl9 for all Sonneratia species indicate that S. lanceolata individuals are nested within S. caseolaris. We found strong genetic structure among geographic regions (South China Sea, the Indian Ocean, and eastern Australia) inhabited by S. caseolaris. We estimated that divergence between the Indo-Malesia and Australasia populations occurred 4.035 million years ago (MYA), prior to the onset of Pleistocene. BARRIERS analysis suggested that complex geographic features in the IWP region had largely shaped the phylogeographic patterns of S. caseolaris. Furthermore, haplotype analyses provided convincing evidence for secondary contact of the South China Sea and the Indian Ocean lineages at the Indo-Pacific boundary. Demographic history inference under isolation and migration (IM) model detected substantial gene flow from the Sri Lanka populations to the populations in the Java Island. Moreover, multi-locus sequence analysis indicated that S. lanceolata was most closely related to the Indian Ocean populations of S. caseolaris and the divergence time between them was 2.057 MYA, coinciding with the onset of the Pleistocene glaciation. Our results suggest that geographic isolation driven by the Pleistocene ice age resulted in the recent origin of S. lanceolata.