Low Genetic Connectivity of Strongly Inbred Ruppia brevipedunculata in Aquaculture Dominated Lagoons (Viet Nam)

Lagoonal environments exhibit high levels of instability depending on hydrological, climatic and ecological factors, thereby influencing the distribution and structure of submerged plant communities. Conditions typically fluctuate widely due to the interaction of freshwater from rivers with saltwater from the sea, as well as from aquaculture activities that together influence submerged hydrophyte community spatial and temporal variability depending on plant survival strategies. Ruppia species feature either underwater pollination mediated by an air bubble or by the release of pollen floating at the water surface, the former promoting self-pollination. Tropical Asian Ruppia brevipedunculata Yu and den Hartog was assumed to pollinate below the water surface and identified as a separate lineage among selfed Ruppia taxa. We used nine nuclear microsatellites to estimate inbreeding levels and connectivity of R. brevipedunculata within a large SE Asian lagoon complex. Ruppia brevipedunculata meadows were strongly inbred as could be derived from the many monomorphic or totally fixed loci for unique alleles in different parts of the lagoon, which appears consistent with selfing behavior. Those from aquaculture ponds were highly inbred (FIS = 0.620), though less than open lagoon sites that showed nearly total inbreeding (FIS = 0.942). Ruppia brevipedunculata from two major lagoon parts were highly differentiated with spatially structured gene pools and a strong barrier between parts of the lagoon over a 30 km distance. Migration-n analysis indicated unidirectional though limited gene flow and following potential hydrological connectivity. Overall, private alleles under homozygote conditions explained a stronger genetic differentiation of populations situated inside aquaculture ponds than of open lagoon populations. Kinship values were only relevant up to 5 km distance in the open lagoon. Within a confined area of aquaculture ponds featuring dense vegetation in stagnant water, there would be opportunity for mixed pollination, thereby explaining the higher diversity of unique multilocus genotypes of aquaculture pond habitats. Low connectivity prevents gene pools to homogenize however promoted sites with private alleles across the lagoon. Complex hydrodynamic systems and human-made habitats enclosed by physical structures impose barriers for propagule dispersal though may create refugia and contribute to conserving regional genetic diversity.

Lagoonal environments exhibit high levels of instability depending on hydrological, climatic and ecological factors, thereby influencing the distribution and structure of submerged plant communities. Conditions typically fluctuate widely due to the interaction of freshwater from rivers with saltwater from the sea, as well as from aquaculture activities that together influence submerged hydrophyte community spatial and temporal variability depending on plant survival strategies. Ruppia species feature either underwater pollination mediated by an air bubble or by the release of pollen floating at the water surface, the former promoting self-pollination. Tropical Asian Ruppia brevipedunculata Yu and den Hartog was assumed to pollinate below the water surface and identified as a separate lineage among selfed Ruppia taxa. We used nine nuclear microsatellites to estimate inbreeding levels and connectivity of R. brevipedunculata within a large SE Asian lagoon complex. Ruppia brevipedunculata meadows were strongly inbred as could be derived from the many monomorphic or totally fixed loci for unique alleles in different parts of the lagoon, which appears consistent with selfing behavior. Those from aquaculture ponds were highly inbred (F IS = 0.620), though less than open lagoon sites that showed nearly total inbreeding (F IS = 0.942). Ruppia brevipedunculata from two major lagoon parts were highly differentiated with spatially structured gene pools and a strong barrier between parts of the lagoon over a 30 km distance. Migration-n analysis indicated unidirectional though limited gene flow and following potential hydrological connectivity. Overall, private alleles under homozygote conditions explained a stronger genetic differentiation of populations situated inside aquaculture ponds than of open lagoon populations. Kinship values were only relevant up to 5 km distance in the open lagoon. Within a confined area of aquaculture ponds featuring dense vegetation in stagnant water, there would be opportunity for mixed pollination, thereby explaining the higher diversity of unique multilocus genotypes of aquaculture pond habitats. Low connectivity prevents gene pools to homogenize however promoted sites with private alleles across the lagoon. Complex hydrodynamic systems and human-made habitats enclosed by physical structures impose barriers for propagule dispersal though may create refugia and contribute to conserving regional genetic diversity.

INTRODUCTION
Physical alterations of the environment due to coastal development (Waycott et al., 2009;Grech et al., 2011;Unsworth et al., 2019), run-off from urban, agricultural and industrial areas (Waycott et al., 2009;Grech et al., 2011;Unsworth et al., 2019) together with ocean warming (Duarte et al., 2018;Kendrick et al., 2019) are the main stressors resulting in worldwide seagrass decline. Ruppia species appear to be more resistant toward environmental perturbations because of their high tolerance toward fluctuations in salinity, temperature and turbidity (Verhoeven, 1979;Mannino et al., 2015), their high colonization potential due the capability of self-pollination, large seed sets (Kendrick et al., 2012), fast rhizome growth rates (Marbà and Duarte, 1998), and the possibility of recovery after disturbance from seed banks (Strazisar et al., 2013). Consequently, Ruppia may replace more sensitive seagrass species and dominate in disturbed habitats (Cho et al., 2009). Nevertheless, with increasing coastal development, landform modifications and environmental perturbations, Ruppia populations may be particularly vulnerable in the future because of their proximity to human-induced land use and thus high exposure to direct and indirect effects of physical disturbance, sedimentation and eutrophication (Waycott et al., 2009;Quiros et al., 2017;Dierick et al., 2021).
The cosmopolitan genus Ruppia, found from northern to tropical regions, consists of submerged annual or perennial plants which occur in a broad array of coastal habitats ranging from brackish to hypersaline water (Short et al., 2007). Ruppia dominated seagrass meadows, which are often monospecific in saltwater lagoons of temperate regions or combined with other submerged macrophytes and seagrasses in brackish water lagoons in the tropics (Verhoeven, 1979), are important habitats for fish, invertebrates and birds and thus essential for the functioning of coastal ecosystems and related ecosystem services (Verhoeven, 1979;Triest et al., 2018a). Ruppia species are monoecious plants capable of both sexual reproduction by bisexual flowers, sometimes regarded as unisexual flowers which occur in pairs on a peduncle, as by asexual reproduction involving the formation of new shoots following rhizome growth (Verhoeven, 1979;den Hartog and Kuo, 2006). Pollination in Ruppia is characterized by hydrophily with either the release of pollen underwater in an air bubble or on the water surface (Cox and Knox, 1989;Kaul, 1993;Lacroix and Kemp, 1997;Yu and den Hartog, 2014). Flexibility in pollen release and capture location has been observed in R. maritima where the location depends on local environmental conditions and the time of pollen maturation (Lacroix and Kemp, 1997;Taylor et al., 2020). Ruppia fruits, which contain negatively buoyant seeds, sink almost immediately to the bottom when ripe and detached, and contribute to local seed banks which most likely have an important role in the yearly replenishment of Ruppia populations (Strazisar et al., 2013;Triest and Sierens, 2013;Mannino et al., 2015). The seed dispersal capacity of Ruppia species is therefore low with the most important dispersal vectors being currents, wave action and tidal flow, though birds can additionally play a role in propagule dispersal (Les et al., 2003;Ito et al., 2010;Triest and Sierens, 2013). In addition, the dispersal of seeds by rafting on detached floating reproductive shoots and fragments is considered to be an important dispersal mechanism especially for long-distance transport in Ruppia and may promote connectivity among distant populations (Orth et al., 1994;Källström et al., 2008;Triest et al., 2018b). Most population genetic studies focused on Ruppia species of Europe, more specifically in the Mediterranean, and recorded low levels of genetic diversity and strong inbreeding in R. maritima L. populations due to high levels of self-fertilization (Triest and Sierens, 2015). Consequently, selfing, which is promoted by bisexual flowers and the location of pollen release and capture, is considered to be the predominant reproduction strategy in those Ruppia species closely related to R. maritima (Triest et al., 2018a) though an incomplete protogyny that facilitates delayed selfing increases the potential for outcrossing before self-pollination can occur (Goodwillie and Weber, 2018;Taylor et al., 2020). Nevertheless, inter-and intraspecific variation in Ruppia has been shown in the contribution of self-fertilization, biparental inbreeding and outcrossing for population growth and sustainability (Triest and Sierens, 2015;Triest et al., 2018a;Beirinckx et al., 2020).
High phenotypic plasticity to cope with fluctuations in salinity, water level, temperature and irradiance (Verhoeven, 1979;Mannino et al., 2015), together with recurrent hybridizations (Beirinckx et al., 2020), historical misinterpretations and local, incomplete descriptions have complicated Ruppia genus taxonomy and phylogeny (Yu and den Hartog, 2014;den Hartog and Triest, 2020). As a result, Ruppia was thought to consist of only two cosmopolitan species (R. maritima and R. spiralis L. ex Dum.-the latter often cited under R. cirrhosa (Petagna) Grande) whose species range sizes were overestimated (Yu and den Hartog, 2014;den Hartog and Triest, 2020). Yu and den Hartog (2014) described two new (sub-)tropical species in Asia (China) using extensive morphological and molecular features (Yu et al., 2020) i.e., R. brevipedunculata Yu and den Hartog and R. sinensis Yu and den Hartog, which were previously considered as R. maritima and R. cirrhosa, respectively. The South East Asian R. brevipedunculata is found in shallow lagoons up to a depth of 1 m where it can grow as an annual together with other Ruppia, seagrass or aquatic plant species. It can be distinguished from related Ruppia species by its very short peduncle which restricts pollen release and capture under the water surface (Yu and den Hartog, 2014). Because of the sexual traits of R. brevipedunculata which promote self-pollination (i.e., bisexual flowers, short peduncle), selfing is expected to be an essential strategy for the maintenance of R. brevipedunculata meadows. We hypothesize that this strategy will result in reduced levels of genetic diversity and strong genetic differentiation among populations.
In this study, we determined the gene diversity and genetic structure within and among R. brevipedunculata meadows to assess the contribution of self-fertilization and inbreeding in a tropical aquaculture-dominated environment. We specifically tested for (1) evidence of selfing as mating system, and (2) connectivity levels and directionality between sites over various distances. Therefore, we developed nuclear microsatellites for the species and studied 12 sites in two distinct though hydrologically connected parts of a tropical lagoon complex (Hue, Central Viet Nam). In addition, we aim to test for a local spatial structure driven by enclosed aquaculture ponds. The consequence of a selfing plant mating system and the potential role of hydrodynamics in lagoons and aquaculture ponds is discussed in a conservation and coastal landscape context.

Study Sites and Plant Materials
Ruppia brevipedunculata plants were collected in 12 sites of the Hue lagoon complex, of which four sites in Tam Giang (TG) and eight sites of the Cau Hai (CH) lagoon (Table 1, Figure 1A). In each site, depending on the number of clumps we collected up to 24 individual shoots along a transect at distance intervals of more than 2 m, thereby putatively avoiding duplicates of a same individual. Individual clumps develop seasonally and were usually smaller than 1 m diameter, sometimes intermingled with other species, namely Najas malesiana Wild (often reported as Najas indica (Willd.) Cham., occasionally with Halophila beccarii Aschers. or Halophila ovalis (R.Br.) Hook.f. We sampled 230 Ruppia individual shoots that were dried on silica gel. The identity of this collection was previously checked on a single sample with cpDNA haplotypes (intron trnH-psbA) (Triest et al., 2018a), ITS sequences corresponding to the ITS-A variant (Triest and Sierens, 2014) and with nuclear microsatellites (Triest et al., 2018a). In the largest lagoon of South East Asia situated near Hue in Central Viet Nam, two main and distant parts of the Hue lagoon were sampled (Table 1, Figure 1A). The maximum geographical distance between sites and among both lagoons was 45 km whereas the distance between sites within a lagoon ranged from 0.94 to 12.6 km for TG and from 0.71 to 6.9 km for CH. All sites of the Hue lagoon were brackish to saltwater (Conductivity from 21.6 to 38.8 mS/cm; Salinity from 13.1 to 24.8 mg/l). The TG lagoon was most impacted by intensive aquaculture activities along shallowest borders (otherwise deep and less suitable habitats for submerged macrophytes) and comprised monospecific R. brevipedunculata vegetations inside the shrimp ponds with turbid water ( Table 1). The so-called high-tide aquaculture ponds were surrounded by an elevated earth wall and a narrow and temporary inlet/outlet. The CH lagoon is a more open system of very shallow (<2 m deep) clear water connected to the sea though influenced by freshwater river input and a brackish water connection (Thuy Tu) from major aquaculture areas of the TG lagoon. Total phosphorus ranged from (0.02-) 0.24 to 0.96 mg/l in TG and from 0.03 to 0.49 mg/l in CH ( Table 1). Using QGIS (v.3.10.3) (QGIS Development Team, 2021) we constructed a digital elevation model (DEM) map of the sampling area (Figures 2A,B). The topography of the surrounding area (<5 km) around both lagoons is characterized by lowlands never exceeding 50 m height except for Ba . ch Mã national park which borders the CH lagoon in the South and South-East. We mapped the water flow direction in both the high tide (Figure 2A) and low tide ( Figure 2B) during the dry season in the lagoon which was adapted from Tran (2005) who determined water flow directionality and velocity by the combination of direct measurements and modeling based on topographic maps. During the dry season at high tide, seawater enters the lagoon complex from the inlet in CH (Tu Hien) and in TG (Thuan An) with the latter generating a North-West current toward TG and a South-East current toward the narrow zone which connects TG and CH (Thuy Tu) (Figure 2A) (Nguyen Tien Lam, 2002;Tran, 2005). In contrast, the Perfume river and other small rivers discharges water in the Tam Giang-Cau Hai lagoon complex during low tide in the dry season and during the rainy season where water travels through the narrow connection from TG to CH at higher velocities (>0.200 m/s) and exits the lagoon complex at Tu Hien in TG and Thuan An in CH ( Figure 2B) (Nguyen Tien Lam, 2002;Tran, 2005).

DNA Extraction, Microsatellite Development, and Amplification
Genomic DNA was extracted from leaf tissue using the E.Z.N.A SP plant DNA Mini Kit (Omega bio-tek, Norcross, GA, USA) as in Triest and Sierens (2010). Microsatellite loci were developed and isolated from a microsatellite-enriched genomic library following an enrichment procedure with Dynabeads. We used RsaI and SspI as restriction enzymes and motifs AG, TG, GATA and GACA. Out of 96 sequences, 39 were identified as potential for designing primers with Primer3 (http://primer3.wi.mit.edu/). From these 39 primer pairs, the purity and quantity of the DNA were determined using a Nanodrop one Spectrophotometer (Thermo Fisher Scientific, Waltham, Massachusetts, USA). An Illumina paired-end library was constructed and sequenced using the Illumina HiSeq platform at Macrogen (Seoul, Republic of Korea). SSR_pipeline (Miller et al., 2013) was used to identify microsatellites and design primers. Out of 18.2 million 100 bp paired-end reads, 2.3 million pairs were successfully joined by the module joinseqs. The module SSR_search found 1547 dinucleotide SSRs with at least 14 repeats, 4,602 trinculeotide SSRs with at least 12 repeats and 3,776 tetranucleotide SSRs with at least 10 repeats. Batchprimer 3 (You et al., 2008) was used for designing primers and 48 primer-pairs were selected for synthesis  Table 1.
Frontiers in Conservation Science | www.frontiersin.org TABLE 1 | Characteristics of 12 sites in the Tam Giang (TG)-Cau Hai (CH) lagoon system (Viet Nam) with N: number of samples; decimal coordinates; conductivity (mS/cm); salinity (‰); total phosphorus (TP in mg/l); habitat features; D pond : distance to nearest aquaculture pond (m) and D ind : distance interval between individual samples (m).

Site
Locality  based on fragment length and number of repeats. We selected nine polymorphic loci presenting di-and trinucleotide repeats flanked by more than 10 nucleotides between primers region and motifs (Supplementary Table 1). Motifs were repeated 9-23 times within allele lengths of 103-189 bp. Using Multiplex Manager (Holleley and Geerts, 2009), one single multiplex reaction of nine amplifiable primer pairs with four different dye-labeled primers (6FAM/VIC/NED/PET) was designed.
Microsatellite amplification of all individuals was done using the abovementioned 9 loci. Multiplex polymerase chain reaction (PCR) conditions for each PCR reaction were: 6.25 µl master mix (Qiagen multiplex PCR kit master mix), 1.25 µl primer mix, 2.5µl H 2 O and 2.5µl of genomic DNA, making a total volume of 12.5 µl. PCR was performed in a thermal cycler  and Bio-Rad MyCycler) with the following conditions: an initial denaturation of 95 • C for 15 min followed by 35 cycles of 30 s denaturation at 95 • C, 90 s annealing at 57 • C and 80 s elongation at 72 • C followed by a final extension of 30 min at 60 • C. PCR products were run on an ABI3730XL sequencer (Macrogen, Seoul, Korea) and fragments were analyzed with GeneMarker (v.2.60) (Holland and Parson, 2011).

Data Treatment
Linkage Disequilibrium, Scoring Errors, and Null-Alleles A linkage test between all pairs of loci (1,000 permutations) was done across all populations to verify any genotypic disequilibrium at the 0.001 level using FSTAT (v.2.9.3) (Goudet, 2001) MICRO-CHECKER (Van Oosterhout et al., 2004) was used to verify data quality, potential null alleles, scoring errors or large allele dropout. This test, however, does not consider the effect of inbreeding situations and therefore INEST 2.2 (Chybicki and Burczyk, 2009) was applied as a Bayesian approach for estimating both null alleles and inbreeding simultaneously. Three parameters were used for the model comparison: n, null alleles; f, inbreeding; b, genotyping failure. With these three parameters, seven models are composed: f, n, b, nf, nb, bf, nfb (e.g., nf means the model includes null alleles and inbreeding but does not include genotyping failures). The models were run with 10,000 burn-ins and 50,000 cycles for each site. Summary values of the Deviance Information Criterion (DIC) were obtained for comparison of seven models to (1) verify whether several zero values scored could be explained from the presence of null alleles and (2) explain the large homogeneity of multilocus genotypes encountered. Post-processing was done for each sample on the full model "nbf " from the marginal posterior distribution of hyper-parameters, namely the empirical average inbreeding coefficient Avg(Fi) with 95% highest posterior density interval (HPD low-HPD high).

Shared Multilocus Genotypes
We calculated the number of multilocus genotypes (MLG's) which are shared among individuals using GenAlEx (v.6.5) (Peakall and Smouse, 2012). In addition, we determined the number and identity of shared MLG's among sites and visually represented the shared MLG's within sites by bar-charts using Poppr (Kamvar et al., 2014). To assess if shared MLG's are the result of an asexual reproduction event in the past or due to high levels of inbreeding, we calculated the probability of occurrence for each MLG (P gen_FIS ) and the probability for each identical pair within a population that they have arisen from distinct sexual reproduction events (p sex_FIS ) taking account of departures from H-W equilibrium using RClone (Bailleul et al., 2016). Moreover, we calculated and tested the significance of the aggregation index (A c ), the probability that closest sampled individuals have identical MLG's relative to randomly drawn pairs in a population.

Genetic Diversity
The probability of identity for a combination of 9 loci was estimated for the full data set of both lagoons using GenAlex. Deviations from Hardy-Weinberg equilibrium per locus and per population were tested using a Chi-Square test (GenAlex). At site and lagoon level, we estimated the mean number of alleles (N A ), effective number of alleles (Ae), observed heterozygosity (Ho), unbiased expected heterozygosity (uHe) corrected as 2N/(2N-1), number of monomorphic and fixed loci (i.e., loci with 1 or more alleles, but lacking any heterozygote conditions) using GenAlex. Allelic richness (Ar) for site and lagoon level, overall F ST and population inbreeding coefficient (F IS ) with p-values after randomization tests (1/1,000) were obtained from FSTAT. We tested for recent bottlenecks in each site under the two-phase model (TPM) with 95% single-step mutations and 5% multiple-step mutations (Wilcoxon's test 1-tailed) using BOTTLENECK 1.2.02 (Piry et al., 1999). Selfing rate (S) in TG and CH were estimated assuming a mixed mating model using SPAGeDi (v.1.5a) (Hardy and Vekemans, 2002). We explored the correlation between environmental features (conductivity, salinity, TP, distance to nearest aquaculture ponds) of 12 sites and their genetic variables (allelic richness, gene diversity, number of monomorphic or fixed loci, F IS ).

Genetic Differentiation and Population Structure
Genetic structure was assessed with a three-level analysis of molecular variance (AMOVA) (Michalakis and Excoffier, 1996) using hierarchical fixation indices for two regions, 12 sites and 230 individuals (GenAlex) giving F RT and the standardized F' RT between two lagoons. Pairwise F ST -values and standardized F' ST -values were calculated between all pairs of sites and F STvalues were tested for significant differentiation using 999 permutations. Isolation-by-distance was tested for pairwise F ST of sites and their geographical distances considering straight flight distances between sites in a Mantel test using 1,000 randomizations (Mantel, 1967) with GenAlex. Additionally, pairwise F ST values were tested for five distance classes (1,000 permutations) and regressed over full range (slope b tested with 1,000 permutations) using SPAGeDi. At the Hue lagoon scale, relative to the whole sample, we estimated Kinship values (F IJ according to (Loiselle et al., 1995) of individuals between sites and their spatial autocorrelation. Five distance classes (up to 2, 5, 20, 35, 45 km) were considered and represented an equal number of individual comparisons for each class, using 1,000 randomizations (SPAGeDi).
With the R-package Adegenet (Jombart, 2008; R Core Team, 2021) we performed a Discriminant Analysis of Principal Components (DAPC). We used the function "find.clusters" as a clustering algorithm to infer the genetic structure on both lagoon (with K = 2 clusters) as site level (with K = 4 clusters). In addition, we assigned individuals to 12 clusters and compared the membership of the newly assigned clusters to our 12 predefined sites using the function "assignplot." A Bayesian clustering method (Pritchard et al., 2000) was carried out using STRUCTURE version 2.3. We tested K in 20 independent runs from 1 to 15 (50,000 burn-in and 500,000 Markov chain Monte Carlo replicates in each run), without using sampling location as a prior. Runs were carried out assuming admixture and an independent model of allele frequencies. The optimal K was inferred with the K statistic (Evanno et al., 2005) and LnPK using Structure Harvester (Earl and vonHoldt, 2012) calculated with StructureSelector (Li and Liu, 2018). To determine zones of sharp genetic changes across the lagoon, the geographic coordinates and nine F ST data matrices of the individual loci were analyzed with the Monmonier's algorithm, implemented in BARRIER 2.2 software (Manni et al., 2004) and allowed for two barriers per locus. Three hypotheses to estimate gene flow were tested with Migrate-n (Beerli, 2006;Beerli and Palczewski, 2010) from the mutation-scaled population sizes (Theta) and immigration rates (M) assuming a constant population size through time. The Brownian model was tested locus by locus along with the product of all distributions of all loci. Uniform prior distribution settings (min, max, delta) were as follows for Theta = 0, 20, 2 and for M = 0, 20, 2. The number of recorded steps was 10 6 at a sampling frequency of 10 4 after an initial burnin. The effective number of immigrants per generation (N e m) was calculated as [Theta×M]/4. More precisely, we considered two pooled R. brevipedunculata populations of Tam Giang and Cau Hai lagoons to test for uni-and bidirectional migration models. The Brownian motion mutation model within each case was adopted for a subsample of 35 individuals of each lagoon part, following the abovementioned settings, computing two replicate chains (with different seed), and using the Bezier thermodynamic integration (Beerli and Palczewski, 2010) for calculation of the Bayes factors from marginal likelihoods giving model probabilities. To calculate the relative migration matrix between the TG and CH lagoon and between the 12 sites, an alternative method divMigrate-online (Sundqvist et al., 2016) was used for the parameter Nm (eqn. 12 from Alcala et al., 2014) considering 1,000 bootstraps (alpha = 0.5).

Microsatellite Interpretation, Data Quality, and Shared Multilocus Genotypes
A linkage test between all pairs of loci (1,000 permutations) was done across all populations to verify any genotypic disequilibrium (Supplementary Table 2) and gave a putative disequilibrium for RMV22 × RV32 (at p < 0.001) and for RV32 × RV40 (at p < 0.05). Both loci contain many homozygotes, though could not be considered as effectively linked. Data quality was verified with MICRO-CHECKER and gave neither scoring errors, nor evidence of large allele dropout. Potential presence of null alleles was suggested in eight out of nine loci. This can be better explained from a consistent very low heterozygosity due to strong inbreeding instead of null alleles because no single missing value that would refer to a homozygous null allele was recorded. Summary values of the Deviance Information Criterion (DIC) as obtained from INEST further allowed for comparison of seven models to explain the large homogeneity encountered or potential null alleles considering "f " the inbreeding; "n" the null alleles; and "b" genotyping failures and all combinations thereof. All models that included inbreeding (f, nf, fb, nbf) gave a comparable lowest DIC value within each of the 12 populations (Supplementary Table 3). For most sites, the "f " inbreeding model could not be calculated because of one or a few samples with zero score for alleles in a locus. Nonetheless, a comparison of those models with the lowest Deviance Information Criterion (DIC) selected nf, bf and nbf (in all 12 populations) as best fit, thereby outperforming the other models and evidencing an effect of inbreeding, hence putting inbreeding instead of null alleles as a most likely explanation for the lack of heterozygotes. The F IS of each site was corrected for putative null alleles and ranged from 0.473 to 0.999 (Table 2), though 10 out of 12 sites reached values exceeding 0.7.
Over the total dataset, 152 individuals shared identical MLG's with sites TG1, CH5, CH7, CH8, and CH12 containing an MLG being present in more than 10 individuals (Supplementary Figure 1). Identical genotypes were identified among sites in CH but none were found among sites in TG or between lagoons TG and CH (Supplementary Table 4). One dominant, fully homozygous multilocus genotype (MLG45) was identified within 71 individuals spanning all sites in CH and was composed of alleles with the highest frequencies at each locus when considering all individuals though only one individual within each site for shared MLG (Supplementary Table 4). A c ranged from 0.41 to 0.54 (mean = 0.47) and was not significant for all populations.
Overall Genetic Diversity at Lagoon Level and Its Two Main Parts (TG and CH) Nine polymorphic microsatellite loci gave interpretable patterns for a total of 49 alleles in the lagoon (Supplementary Table 1). We found a low probability of identity for 9 loci when considering both lagoons (PI = 0.000005 in TG; PI = 0.001 in CH). Regarding the total population, the allelic richness reached Ar = 5.4 though a low effective size (Ae = 2.3) and very low heterozygosity levels were observed (Ho = 0.047; He = 0.551) ( Table 2). A very high level of inbreeding with F IS = 0.991 might be subject to a Wahlund effect and therefore was reconsidered at smaller geographic scales.
Both major lagoon areas, Tam Giang (TG) and Cau Hai (CH), contained a nearly similar number of alleles (A = 38 and 37, respectively) and allelic richness (Ar = 4.2 and 3.9, respectively) though differed in observed gene diversity (Ho = 0.112 and 0.016 respectively; p = 0.012, one-sided test) ( Table 2). CH lagoon showed more cases of monomorph loci, of fixed loci under homozygous condition and higher levels of inbreeding than TG lagoon ( Table 2), with F IS = 0.953 (in CH) and F IS = 0.633 (in TG) at p = 0.014 (one-sided test). Selfing rate estimations were very high within both TG (S = 0.90, SE = 0.02; p < 0.001) and CH (S = 0.94, SE = 0.07; p < 0.001) lagoons. Significant recent bottlenecks could be detected for sites TG3, CH6, CH7, and CH11. Out of the 49 alleles, 23 were unique to either part of the Hue lagoon complex. In the TG lagoon 12 private alleles at q = 0.007-0.620 (mean q = 0.165) and in CH lagoon 11 private alleles at q = 0.003-0.884 (mean q = 0.253) could be found.

Genetic Diversity at Site Level (Transects)
At site (transect) level, R. brevipedunculata harbored 17-29 alleles out of the 49 for the whole lagoon and showed Ae ranging from 1.1 to 2.0 and Ar from 1.6 to 2.7 ( Table 2). R. brevipedunculata sites showed very low to low gene diversities with H O ranging from 0 to 0.196 and unbiased He from 0.110 to 0.476 (Table 2), thereby giving high inbreeding values with F IS ranging from 0.457 to 1 (i.e., total fixation). Within-site inbreeding values consistently were significantly high, especially in CH sites and can be explained from the occurrence of several polymorphic loci in which alleles were only present in homozygous conditions ( Table 2). A comparison between aquaculture pond habitats and open lagoon sites revealed lower levels of inbreeding (F IS = 0.620 and F IS = 0.942, respectively at p = 0.014) and consequently higher levels of observed heterozygosity (H O = 0.151 and H O = 0.018, respectively at p = 0.004) in aquaculture ponds. No correlation was found between environmental features (conductivity, salinity, TP, distance to nearest aquaculture ponds) of 12 sites and the genetic variables (allelic richness, gene diversity, number of monomorphic or fixed loci, F IS ).

Population Genetic Structure and Differentiation
A three-level AMOVA-F ST for 2 regions (TG and CH lagoon) and 12 sites resulted in strong differentiation at all hierarchical levels with 43% of the variance explained among lagoons and F RT = 0.472, F SR = 0.211, F ST = 0.548, F IS = 0.864, F IT = 0.938 (all at p < 0.001) ( Table 3). A maximum F RT of 0.553 gave a standardized F' RT = 0.772 between both lagoons ( A comparison test showed significant (p = 0.003) higher differentiation of the aquaculture sites within TG (F ST = 0.497) than for sites within the CH lagoon (F ST = 0.0035) due to many private alleles under homozygote conditions. Isolation-by-distance was obtained (Mantel test, p = 0.001), with high predictive value (y = 0.014x+0.034; R 2 = 0.80) (Supplementary Figure 2). Pairwise F ST testing showed significantly lower values than the average (F ST = 0.353) at shortest distance classes up to 3.8 km (F ST = 0.064; p = 0.002) and 5.4 km (F ST = 0.109; p = 0.004) though significantly higher values for large distance classes above 32.5 km (F ST = 0.348 p > 0.05) and up to 37 km (F ST = 0.543; p = 0.013) and 49 km (F ST = 0.677; p < 0.001). Regression with linear distance at locus level gave a significant positive slope for 8 out of 9 loci which strengthens the evidence of IBD (overall slope b = 0.015 with R 2 = 0.78; p < 0.001). Kinship values (F IJ ) of individuals among sites (thus omitting within-site comparisons, to avoid adverse effects from local inbreeding) were significantly higher than the average kinship value, relative to the whole sample, over shortest distances of 2 km (F IJ = 0.17; p = 0.002) and 5 km (F IJ = 0.14; p = 0.006) (Supplementary Figure 3).
DAPC results at lagoon level (K = 2) showed a strong genetic structure characterized by two genetically distinct clusters ( Figure 3A) whose membership was consistent with those sampled in both main parts in the lagoon (126 individuals (81.3%) of CH in cluster1, all 75 individuals of TG in cluster2), except for 29 individuals (18.70%) of CH which were assigned to cluster2 (Figure 3B). At site level (K = 4), a similar pattern was identified with cluster3, including only individuals of sites in CH, being genetically distinct from the other clusters containing mostly TG individuals (Figures 3C,D). These results verify the low level of connectivity between TG and CH although a low number of individuals in all CH populations (except CH6) were assigned to TG clusters potentially reflecting migrants from TG but not vice versa (Figures 3B,D). Genetically distinct clusters, nearly consistent with individuals in TG1 (cluster1), TG2 (cluster4), TG3 and TG4 (cluster2), show strong population structure within TG and consequently, sites inside aquaculture ponds appear to be genetically isolated (Figures 3C,D). With 12 new clusters, most individuals of TG were reassigned to their predefined sites (Figure 3E). In contrast, individuals of CH were assigned over CH sites without distinct patterns and with only a few individuals being fully reassigned to their predefined site. Low genetic differentiation of CH sites and the absence of a population structure across CH, together with a high proportion of shared and fully homozygous MLG's among CH sites indicate high levels of genetic connectivity among these open lagoon sites. The model-based clustering method as implemented in STRUCTURE showed K = 2 [LnP(K) = −2,762] because the K was highest and standard deviation was lowest. These two clusters corresponded largely to the TG and the CH lagoon, indicating their different gene pools ( Figure 1B). We further implicitly considered the inbred nature of the species. The LnP(K) values increased for K = 3 [LnP(K) = −2,353] and K = 4 [LnP(K) = −2,143] after which the larger K values converged to a plateau (Figures 1C,D). We also presented K = 3 and K = 4 because additional clusters distinguished isolated populations inside aquaculture ponds. The scenarios of K = 2, K = 3, and K = 4 allowed to indicate several individuals of another gene pool within the CH cluster, potentially reflecting migrants coming from TG which hereafter was tested in a migration analysis.
The strongest barriers for most loci and their highest pairwise F ST were indicated between TG and CH lagoon and between neighboring TG populations (TG1, TG2, and TG3) located inside aquaculture ponds (Figures 2A,B). The specific testing with migrate-n on directionality for R. brevipedunculata across both lagoon parts, indicated that a unidirectional stepping stone model from TG to CH (Bezier value = −1,345,084) appeared more likely compared to a unidirectional model from CH to TG (Bezier value = −1,385,917) and to the bidirectional model (Bezier value = −1,452,903) (Supplementary Table 6), although the estimated gene flow value from TG toward CH was low (N e m = 0.6). The migration between the TG and CH lagoon using divMigrate gave a significant direction from TG to CH with a relative value 1. The migration between all 12 sites of the area gave only a significant direction from CH10 to CH11 and with a relative value of 0.23 Nm.

High Levels of Selfing and Inbreeding in Open Lagoon Sites and Aquaculture Ponds
The very high levels of inbreeding in all sites with selfing rate estimations exceeding 0.90 are presumed to be the result of the self-pollination mating system of R. brevipedunculata. These findings are consistent with a broad-scale study on R. maritima populations in the Mediterranean where low genetic diversity, a high number of monomorphic and fixed loci and significant though lower levels of inbreeding were attributed to selffertilization (Triest and Sierens, 2015). R. cirrhosa populations in the Mediterranean were also significantly inbred but lacking fixed loci when compared to R. maritima which suggests that outcrossing, though with biparental inbreeding, is favored over self-pollination in R. cirrhosa (Triest et al., 2018b). In contrast, populations of Halophila beccarii studied in the same lagoon complex (Cau Hai) did not show evidence of inbreeding or selfing and were maintained by local seed recruitment and clonal extension (Phan et al., 2017). The strong ability of selffertilization in Ruppia is strongly determined by the location of pollen release, with hypohydrophily hypothesized to be less efficient for outcrossing than pollen transfer along the water surface (Taylor et al., 2020). Due to the very short peduncle, R. brevipedunculata may be restricted to pollen release and pollen capture under the water surface which most likely explains selffertilization as the most prevalent pollination strategy for R. brevipedunculata in our study. The ability of selfing in Ruppia may increase seed-set (Les, 1988;Taylor et al., 2020) and provides reproductive assurance when mates are scarce or vectors for pollen flow are absent or inadequate (Wright et al., 2013;Razanajatovo et al., 2016). Autogamy could therefore promote the establishment and sustainability of R. brevipedunculata populations in enclosed and isolated habitats.
Somewhat less elevated inbreeding values within high tide ponds in aquaculture sites when compared to open lagoon populations indicate that pollen flow among distinct but potentially related individuals appears to be promoted within confined aquaculture areas. Dense vegetation in stagnant water may explain the elevated levels of pollen flow between individuals in aquaculture ponds whereas pollen in open lagoon sites could be lost when not immediately used upon release for self-fertilization (Van Vierssen et al., 1984;Reusch, 2003). Geitonogamy could be an additional and important form of self-fertilization alongside autogamy (Reusch, 2001;Taylor et al., 2020), especially when pollen are released between entangled branches of a same plant. Due to the physical enclosure and dense patches within high-tide ponds, both local outcrossing, evidenced by the lower inbreeding values and the higher number of multilocus genotypes in these particular man-made habitats, and geitonogamy can be hypothesized. If outcrossing is possible, cases of geitonogamy may also occur.
Thus far, we have interpreted the obtained levels and patterns of genetic diversity and structure of R. brevipedunculata as the result of an overarching effect of selfing and inbreeding. Nevertheless, this species showed a high proportion of shared multilocus genotypes (MLGs) within sites as well as across open lagoon sites, a phenomenon usually indicative of asexual reproduction and vegetative dispersal (Arnaud-Haond et al., 2007;McMahon et al., 2017;Dierick et al., 2021). Due to high P gen_FIS and P sex _ FIS values, non-significant aggregation of identical MLGs and the presence of several repeated genotypes that are completely homozygous, although they are widely shared among individuals, we assume that multilocus genotypes similarities in our study are the result from inbreeding events. This assumption is further reinforced by the knowledge about the mating system of R. brevipedunculata, the description of R. brevipedunculata being a species that grows annually from seeds (Yu and den Hartog, 2014) and our sampling strategy that avoided clumps of shoots. The possibility of inadvertent re-sampling or the existence of larger individuals cannot be completely ruled out, but after careful consideration, we maintain the view that repetitive MLGs are the result of selfing/inbreeding events based on the observations and interpretations discussed below.
Nearly all MLGs were found in a fully homozygous state for all 9 loci, despite an intrinsic allelic diversity at site and lagoon level. Only twice, within the same site, did a repeated genotype (three shoots) contain a single heterozygous locus. Exceptional to observe was that one MLG occurred 71 times in all eight sampled meadows of the open lagoon, although it showed complete homozygous conditions for the most common allele from each locus (when considering the allele frequencies of a subset omitting the repeated MLGs, this to avoid bias when including so many repeats). On the contrary, in case of clonal growth and vegetative propagation, one should expect a much higher number of heterozygous clones (i.e., heterozygous in at least one locus but rather any combination of up to 9 loci) than from fully homozygous clones. Thus, shared MLGs were extremely skewed toward a homozygous state, despite the presence of numerous unique MLGs harboring one to multiple polymorphic loci. It seems highly unlikely that the shared MLGs were so disproportionate by chance or selection. Such homogeneity of the common alleles as we observed in R. brevipedunculata is therefore best explained by selfing behavior and inbreeding over several generations rather than clonal growth and asexual spread. Thus, in addition to a lack of heterozygosity, R. brevipedunculata meadows were also characterized by completely fixed and monomorphic loci, a combination of characteristics as might be expected from an inbreeding behavior of different strengths and historical origins. This situation is very similar to the inbred R. maritima from Europe (Triest and Sierens, 2015) but completely different from the outcrossed R. spiralis (Triest et al., 2018b). An alternative explanation, such as bias in the sampling design performed (e.g., spatially clustered patches) is not expected to hold here, as individual shoots were taken at a distance interval of 2 m which is sufficient to avoid repeated sampling of the same individual. Meadows of Ruppia were spatially monitored over a gradient of distance classes. The 12 sites were at least 0.7 km but usually several km or tens of km apart and in fact represented replica sites within two lagoons, all sites of which exhibited nearly similar genetic diversity levels, similar patterns, and consistently elevated levels of inbreeding levels, giving some widespread occurrence of the homozygous state of shared MLGs. Occasionally, population genetic studies may suffer from an unfortunate choice of uninformative microsatellite markers yielding large amounts of homozygous genotypes, e.g., a microsatellite locus with a very common allele at near maximum frequency and one or a few alleles at very low frequencies. However, we used nine microsatellite loci each containing 2-8 alleles in a lagoon and were recovered under different heterozygous conditions and in all locations. Allele lengths from 9 to 23 motif repeats were appropriately sized for consistent PCR amplification. Scoring failures or null alleles was not a problem in this study as shown by the INEST analysis. Crossamplification of 24 primers for microsatellite loci of other Ruppia species on R. brevipedunculata was possible but gave only a single amplification product (Triest et al., 2018a), hence an indication of a reorganized nuclear genome and monomorphism of R. brevipedunculata.
Therefore, in view of the above discourse, we hold the view that R. brevipedunculata in the studied lagoon system was subject to selfing and inbreeding events explaining homozygosity and homogeneity of MLGs. Alternative explanations such as clonal growth and dispersal, severe bias in sampling design and uninformative microsatellite loci are considered less likely or even unlikely to explain the patterns obtained. We further hypothesize that subsequent self-pollination could lead to a range of inbred lines, some even spanning an entire genome, and this should be explored in future studies using nuclear single nucleotide polymorphic markers covering a wide range of the genome.

Low Connectivity Among Tam Giang and Cau Hai and Between Aquaculture Ponds
Despite being physically connected by water, the main lagoon parts of the Tam Giang-Cau Hai complex are highly differentiated which is confirmed by an overall high fixation index among these regions, high values of pairwise genetic differentiation among sites between TG and CH, nearly half of the total alleles being private for one of both lagoon parts and a strong spatial structure based on the DAPC (K = 2) and STRUCTURE results. In both TG and CH the mean number of alleles is higher than 4 with the allele number for each locus exceeding 2. Consequently, we hypothesize that the presence of viable, but highly inbred and differentiated populations within the lagoon and aquaculture ponds are most likely the result of various founders who colonized empty patches and established populations predominantly by self-fertilization as a result of a small number of independent historical introductions in CH and TG or multiple historical recurrent introductions in the total Tam Giang-Cau Hai lagoon complex followed by subsequent inbreeding which led to a different subset of alleles in both parts. The subsequent yearly replenishment and recruitment may then originate from local seed banks (Strazisar et al., 2013). A pattern of isolation-by-distance was detected with lowest F ST -values and highest relative migration estimates among closely sampled sites in CH, thereby largely explaining the lower levels of genetic differentiation over small distances. Four sites at close vicinity within CH (CH8, CH10, CH11, and CH12) were not differentiated from each other, being genetically identical or nearly identical inbred exhibiting homozygous loci. Consequently, highly inbred seeds are dispersed and mixed among open lagoon sites in CH leading to the absence of population structure and an overall higher relatedness between individuals up to 5 km compared to the total Tam Giang-Cau Hai system. In contrast, sites within aquaculture ponds in TG were strongly structured and highly differentiated from each other except for closely sampled sites TG3 and TG4. Hence, aquaculture ponds such as these high-tide systems surrounded by an elevated earth wall, act as strong barriers for propagule dispersal of R. brevipedunculata.
The low genetic connectivity between both main lagoon parts over a distance exceeding 30 km's for R. brevipedunculata is most likely explained by the low seed dispersal capacity of the Ruppia genus which fully ripened fruits almost directly sink to the bottom when detached (Triest and Sierens, 2013), however, seedlings may become detached and get distributed throughout a water column (Koch et al., 2010). Besides the relatively large distance between TG and CH, outlets to the open sea and water discharge from the Perfume river in TG during high tide and the rainy season may impose physical barriers for the dispersal of seedlings or ripening fruit-bearing shoots and most likely contributes to the low observed genetic connectivity between both parts. Nevertheless, migration was proposed from TG to CH but the estimated number of migrants per generation appeared low (N e m = 0.6). Such one-way migration of dispersal propagules aligns with the pattern of unidirectional water flow (from TG to CH) found in the narrow connection during both low and high tide which makes it nearly impossible for propagules to disperse from CH to TG. An estimated dispersal distance of 20 km for R. cirrhosa in the Mediterranean Sea was attributed to detached seedlings, rooted shoot fragments containing ripening fruits that connected sites over these larger distances but only within hydrologically connected lagoon areas (Triest et al., 2018b). Likewise, detached fragments may promote gene flow for R. brevipedunculata among distant lagoon parts though no direct evidence is available to support this. Next to currents and tidal waves, birds are reported to be important or occasional dispersal vectors for Ruppia (Ito et al., 2010;Martínez-Garrido et al., 2017). A high number of bird species (S = 73), with a large proportion being migratory (nearly 50%), have been recorded in the estuaries in the North-West of TG (mainly Anas sp. and Egretta sp.) (Tuan 2012). We suggest that the dispersal of seeds by birds cannot fully explain the pattern of unidirectional gene flow between TG and CH. Still, we do not rule out birds as potential dispersal vectors, though their prevalence and effect on population structure is probably low and difficult to detect compared to the hydrodynamics within the lagoon (Triest and Sierens, 2013). In addition, a large proportion of people that live in and around the lagoon depends on fishing activities which may promote connectivity within and between lagoon parts, though more information is needed on the adherence and transport capacity of Ruppia fragments by boats, nets and other fishing gear.

Ecological Significance and Conservation
From a phylogeographical perspective, it is important to consider coastal landform changes in recent historical times as they might shed light on the dynamics that populations have experienced. For a long time, Tu Hien (CH) was the only direct connection to the open sea in the lagoon complex and water was mainly transported unidirectionally from the estuaries in TG to the outlet in CH (Nguyen Tien Lam, 2002). Since the early fifteenth century, a new inlet in TG (Thuan An) was created which changed the hydrodynamics in the lagoon due to the presence of an additional exit for water discharge during low tide and floods and the formation of two reverse currents during high tide, one toward the estuaries in TG and another toward CH (Nguyen Tien Lam, 2002). Our study indicates that this new man-made connection along with the creation of high-tide aquaculture ponds affect the population structure and genetic diversity within and among lagoon parts by decreasing genetic connectivity, reducing population viability and adaption potential toward environmental perturbations (Takebayashi and Morrell, 2001;Wright et al., 2013). However, a higher number of unique multilocus genotypes found in Tam Giang compared to the open lagoon system in Cau Hai may indicate that man-made habitats enclosed by physical structures can provide refuge for rare alleles and genotypes and thus enhance genotypic diversity relative to highly inbred populations. Nevertheless, we assume that the positive effects of preserving rare alleles at isolated sites can be neglected compared to the negative effects of reduced gene flow.
The wide variety in survival strategies (e.g., selfing, outcrossing, clonal persistent, etc.) and their effect on population dynamics and diversification have been overlooked for the genus Ruppia due to largely focusing on two wrongly assumed cosmopolitan Ruppia species in the past (Mannino et al., 2015;Triest et al., 2018a). Correct taxon identification and comprehensive molecular studies at different spatial scales focusing on a variety of Ruppia lineages are therefore needed to understand how Ruppia populations are maintained and persist, which is essential information for the design and application of effective conservation measures to mitigate the effects of habitat deterioration and climate change. Our study highlights the importance of self-fertilization and local recruitment for population maintenance and sustainability for R. brevipedunculata, as recorded in other Ruppia species, though the potential for connectivity beyond the genetic neighborhood should be explored in future studies (Triest and Sierens, 2015;Triest et al., 2018a). Recurrent selfing and mating among relatives may result in inbreeding depression and the loss of populations in the long term due to low effective population sizes and the accumulation of deleterious mutations (Takebayashi and Morrell, 2001;Wright et al., 2013). Therefore, information on the genetic and demographic consequences of long-term inbreeding for R. brevipedunculata is essential to assess if self-fertilization is an advantageous strategy both in the short and long term.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

FUNDING
This research was supported by the Vrije Universiteit Brussel (BAS42), the VLIR-International University Cooperation program 2013-2022 with Hue University and the Erasmus Mundus Lotus Action II Program (NQH research fellowship at Vrije Universiteit Brussel.