Comparative population genetics of congeneric limpets across a biogeographic transition zone reveals common patterns of genetic structure and demographic history.

The distribution of genetic diversity is often heterogeneous in space, and it usually correlates with environmental transitions or historical processes that affect demography. The coast of Chile encompasses two biogeographic provinces and spans a broad environmental gradient together with oceanographic processes linked to coastal topography that can affect species' genetic diversity. Here, we evaluated the genetic connectivity and historical demography of four Scurria limpets, S. scurra, S. variabilis, S. ceciliana and S. araucana, between ca. 19° S and 53° S in the Chilean coast using genome-wide SNPs markers. Genetic structure varied among species which was evidenced by species-specific breaks together with two shared breaks. One of the shared breaks was located at 22-25° S and was observed in S. araucana and S. variabilis, while the second break around 31-34° S was shared by three Scurria species. Interestingly, the identified genetic breaks are also shared with other low-disperser invertebrates. Demographic histories show bottlenecks in S. scurra and S. araucana populations and recent population expansion in all species. The shared genetic breaks can be linked to oceanographic features acting as soft barriers to dispersal and also to historical climate, evidencing the utility of comparing multiple and sympatric species to understand the influence of a particular seascape on genetic diversity.

An Out-of-Patagonia migration explains the worldwide diversity and distribution of Saccharomyces eubayanus lineages.

Population-level sampling and whole-genome sequences of different individuals allow one to identify signatures of hybridization, gene flow and potential molecular mechanisms of environmental responses. Here, we report the isolation of 160 Saccharomyces eubayanus strains, the cryotolerant ancestor of lager yeast, from ten sampling sites in Patagonia along 2,000 km of Nothofagus forests. Frequency of S. eubayanus isolates was higher towards southern and colder regions, demonstrating the cryotolerant nature of the species. We sequenced the genome of 82 strains and, together with 23 available genomes, performed a comprehensive phylogenetic analysis. Our results revealed the presence of five different lineages together with dozens of admixed strains. Various analytical methods reveal evidence of gene flow and historical admixture between lineages from Patagonia and Holarctic regions, suggesting the co-occurrence of these ancestral populations. Analysis of the genetic contribution to the admixed genomes revealed a Patagonian genetic origin of the admixed strains, even for those located in the North Hemisphere. Overall, the Patagonian lineages, particularly the southern populations, showed a greater global genetic diversity compared to Holarctic and Chinese lineages, in agreement with a higher abundance in Patagonia. Thus, our results are consistent with a likely colonization of the species from peripheral glacial refugia from South Patagonia. Furthermore, fermentative capacity and maltose consumption resulted negatively correlated with latitude, indicating better fermentative performance in northern populations. Our genome analysis, together with previous reports in the sister species S. uvarum suggests that a S. eubayanus ancestor was adapted to the harsh environmental conditions of Patagonia, a region that provides the ecological conditions for the diversification of these ancestral lineages.

Pre-domestication bottlenecks of the cultivated seaweed Gracilaria chilensis.

Gracilaria chilensis is the main cultivated seaweed in Chile. The low genetic diversity observed in the Chilean populations has been associated with the over-exploitation of natural beds and/or the founder effect that occurred during post-glacial colonization from New Zealand. How these processes have affected its evolutionary trajectory before farming and incipient domestication is poorly understood. In this study, we used 2232 single nucleotide polymorphisms (SNPs) to assess how the species' evolutionary history in New Zealand (its region of origin), the founder effect linked to transoceanic dispersion and colonization of South America, and the recent over-exploitation of natural populations have influenced the genetic architecture of G. chilensis in Chile. The contrasting patterns of genetic diversity and structure observed between the two main islands in New Zealand attest to the important effects of Quaternary glacial cycles on G. chilensis. Approximate Bayesian Computation (ABC) analyses indicated that Chatham Island and South America were colonized independently near the end of the Last Glacial Maximum and emphasized the importance of coastal and oceanic currents during that period. Furthermore, ABC analyses inferred the existence of a recent and strong genetic bottleneck in Chile, matching the period of over-exploitation of the natural beds during the 1970s, followed by rapid demographic expansion linked to active clonal propagation used in farming. Recurrent genetic bottlenecks strongly eroded the genetic diversity of G. chilensis prior to its cultivation, raising important challenges for the management of genetic resources in this incipiently domesticated species.

Genomic diversity and demographic history of the Dromiciops genus (Marsupialia: Microbiotheriidae).

Three orders represent the South American fauna of marsupials. Of these, Microbiotheria was until recently known as a monotypic genus with the only surviving species Dromiciops gliroides (monito del monte). The recent proposal of a new Dromiciops species (Dromiciops bozinovici), together with new information on the origin and diversification of living microbioterians has changed the prevailing paradigm around the evolutionary history of these emblematic marsupials. Here, we used a RADseq approach to test for evidence of admixture and past or current gene flow among both species of Dromiciops and evaluate the genetic structure within D. gliroides. We analyzed 127 samples of Dromiciops distributed across the known distribution range of both species. We also inferred the joint demographic history of these lineages, thus corroborating the status of D. bozinovici as a distinct species. Demographic history reconstruction indicated that D. bozinovici diverged from D. gliroides around 4my ago and has remained isolated and demographically stable ever since. In contrast, D. gliroides is subdivided into three subclades that experienced recent expansions and moderate gene flow among them (mostly from north to south). Furthermore, genetic distances among populations within D. gliroides were significantly correlated with geographic distances. These results suggest that some of the D. gliroides populations would have survived in glacial refuges, with posterior expansions after ice retreat. Our results have important implications for the systematics of the genus and have profound conservation consequences for the new species, especially considering the fragmentation level of the temperate rainforest.

DNA reconciles morphology and colouration in the drunk blenny genus Scartichthys (Teleostei: Blenniidae) and provides insights into their evolutionary history.

The blenniids of the genus Scartichthys are one of the most common fishes of Central and South American Pacific coastal reefs. This being said, Scartichthys spp. remain difficult to identify in the field, and identification is particularly challenging across the c. 6000 km where three of the four currently accepted species are known to occur in sympatry. A reason for this is that the main taxonomic characters from traditional taxonomy are indeed elusive. In addition, at the same time, species can display multiple colour patterns in the field, depending on their ontogenetic stage, habitat association and reproductive behaviour. Overall, molecular characterization is warranted to help address these issues. In this study, the authors have used a novel approach to revise the genus by combining colouration, morphological and molecular data of representative specimens of the four currently valid species and seven described colour patterns. From this, the authors show that only three of the four species should be considered as valid; Scartichthys gigas (Steindachner, 1876), Scartichthys variolatus (Valenciennes, 1836) and Scartichthys viridis (Valenciennes, 1836), whereas Scartichthys crapulatus (Williams, 1990) should be synonymized with S. viridis. In the same way, the analyses in this study show that one of the colour patterns attributed so far only to S. gigas is characteristic of the juvenile stages of S. viridis. The time-calibrated phylogeny of this study shows that this genus is relatively young and that the estimated time of divergence between S. gigas and S. viridis is c. 1.71 Ma. In comparison, the Desventuradas and Juan Fernandez Islands endemic S. variolatus diverged c. 1.95 Ma. The results of this study help to clarify the taxonomy of Scartichthys.

Stochastic nature of larval dispersal at sea.

The movement of individuals across landscapes remains a fundamental process in population and community ecology. All species have developed a capacity to disperse but this process remains elusive in organisms with complex life-cycles, and none more so than in the marine environment. Here, most organisms have developed a two-phased life-cycle, leaving the risky business of dispersing through the open ocean to their very small and intractable larval offspring. To this day, quantifying dispersal patterns in marine seascapes remains a significant challenge, and yet it is critical to the way we preserve marine ecosystems and the services they provide. In this issue of Molecular Ecology, Catalano et al. (2021) present one of the first longitudinal studies to demonstrate the stochastic nature of larval dispersal. Their work challenges some of our current ideas about marine population connectivity and provides new methodological insights to study its temporal dimension.

Genomic landscape of geographically structured colour polymorphism in a temperate marine fish.

The study of phenotypic variation patterns among populations is fundamental to elucidate the drivers of evolutionary processes. Empirical evidence that supports ongoing genetic divergence associated with phenotypic variation remains very limited for marine species where larval dispersal is a common homogenizing force. We present a genome-wide analysis of a marine fish, Labrus bergylta, comprising 144 samples distributed from Norway to Spain, a large geographical area that harbours a gradient of phenotypic differentiation. We analysed 39,602 biallelic single nucleotide polymorphisms and found a clear latitudinal gradient of genomic differentiation strongly correlated with the variation in phenotypic morph frequencies observed across the North Atlantic. We also detected a strong association between the latitude and the number of loci that appear to be under divergent selection, which increased with differences in coloration but not with overall genetic differentiation. Our results demonstrate that strong reproductive isolation is occurring between sympatric colour morphs of L. bergylta found at the southern areas and provide important new insights into the genomic changes shaping early stages of differentiation that might precede speciation with gene flow.

The biogeography of Dromiciops in southern South America: Middle Miocene transgressions, speciation and associations with Nothofagus.

The current distribution of the flora and fauna of southern South America is the result of drastic geological events that occurred during the last 20 million years, including marine transgressions, glaciations and active vulcanism. All these have been associated with fragmentation, isolation and subsequent expansion of the biota, south of 35°S, such as the temperate rainforest. This forest is mostly dominated by Nothofagus trees and is the habitat of the relict marsupial monito del monte, genus Dromiciops, sole survivor of the order Microbiotheria. Preliminary analyses using mtDNA proposed the existence of three main Dromiciops lineages, distributed latitudinally, whose divergence was initially attributed to recent Pleistocene glaciations. Using fossil-calibrated dating on nuclear and mitochondrial genes, here we reevaluate this hypothesis and report an older (Miocene) biogeographic history for the genus. We performed phylogenetic reconstructions using sequences from two mitochondrial DNA and four nuclear DNA genes in 159 specimens from 31 sites across Chile and Argentina. Our phylogenetic analysis resolved three main clades with discrete geographic distributions. The oldest and most differentiated clade corresponds to that of the northern distribution (35.2°S to 39.3°S), which should be considered a distinct species (D. bozinovici, sensu D'Elía et al. 2016). According to our estimations, this species shared a common ancestor with D. gliroides (southern clades) about ~13 million years ago. Divergence time estimates for the southern clades (39.6°S to 42.0°S) ranged from 9.57 to 6.5 Mya. A strong genetic structure was also detected within and between clades. Demographic analyses suggest population size stability for the northern clade (D. bozinovici), and recent demographic expansions for the central and southern clades. All together, our results suggest that the diversification of Dromiciops were initiated by the Middle Miocene transgression (MMT), the massive marine flooding that covered several lowlands of the western face of Los Andes between 37 and 48°S. The MMT resulted from an increase in global sea levels at the Miocene climatic optimum, which shaped the biogeographic origin of several species, including Nothofagus forests, the habitat of Dromiciops.

Strong habitat and weak genetic effects shape the lifetime reproductive success in a wild clownfish population.

The relative contributions of environmental, maternal and additive genetic factors to the Lifetime reproductive success (LRS) determine whether species can adapt to rapid environmental change. Yet to date, studies quantifying LRS across multiple generations in marine species in the wild are non-existent. Here we used 10-year pedigrees resolved for a wild orange clownfish population from Kimbe Island (PNG) and a quantitative genetic linear mixed model approach to quantify the additive genetic, maternal and environmental contributions to variation in LRS for the self-recruiting portion of the population. We found that the habitat of the breeder, including the anemone species and geographic location, made the greatest contribution to LRS. There were low to negligible contributions of genetic and maternal factors equating with low heritability and evolvability. Our findings imply that our population will be susceptible to short-term, small-scale changes in habitat structure and may have limited capacity to adapt to these changes.

Coping with Pleistocene climatic fluctuations: Demographic responses in remote endemic reef fishes.

Elucidating demographic history during the settlement of ecological communities is crucial for properly inferring the mechanisms that shape patterns of species diversity and their persistence through time. Here, we used genomic data and coalescent-based approaches to elucidate for the first time the demographic dynamics associated with the settlement by endemic reef fish fauna of one of the most remote peripheral islands of the Pacific Ocean, Rapa Nui (Easter Island). We compared the demographic history of nine endemic species in order to explore their demographic responses to Pleistocene climatic fluctuations. We found that species endemic to Rapa Nui share a common demographic history, as signatures of population expansions were retrieved for almost all of the species studied here, and synchronous demographic expansions initiated during the last glacial period were recovered for more than half of the studied species. These results suggest that eustatic fluctuations associated with Milankovitch cycles have played a central role in species demographic histories and in the final stage of the community assembly of many Rapa Nui reef fishes. Specifically, sea level lowstands resulted in the maximum reef habitat extension for Rapa Nui endemic species; we discuss the potential role of seamounts in allowing endemic species to cope with Pleistocene climatic fluctuations, and we highlight the importance of local historical processes over regional ones. Overall, our results shed light on the mechanisms by which endemism arises and is maintained in peripheral reef fish fauna.

First genealogy for a wild marine fish population reveals multigenerational philopatry.

Natal philopatry, the return of individuals to their natal area for reproduction, has advantages and disadvantages for animal populations. Natal philopatry may generate local genetic adaptation, but it may also increase the probability of inbreeding that can compromise persistence. Although natal philopatry is well documented in anadromous fishes, marine fish may also return to their birth site to spawn. How philopatry shapes wild fish populations is, however, unclear because it requires constructing multigenerational pedigrees that are currently lacking for marine fishes. Here we present the first multigenerational pedigree for a marine fish population by repeatedly genotyping all individuals in a population of the orange clownfish (Amphiprion percula) at Kimbe Island (Papua New Guinea) during a 10-y period. Based on 2927 individuals, our pedigree analysis revealed that longitudinal philopatry was recurrent over five generations. Progeny tended to settle close to their parents, with related individuals often sharing the same colony. However, successful inbreeding was rare, and genetic diversity remained high, suggesting occasional inbreeding does not impair local population persistence. Local reproductive success was dependent on the habitat larvae settled into, rather than the habitat they came from. Our study suggests that longitudinal philopatry can influence both population replenishment and local adaptation of marine fishes. Resolving multigenerational pedigrees during a relatively short period, as we present here, provides a framework for assessing the ability of marine populations to persist and adapt to accelerating climate change.

Widespread hybridization and bidirectional introgression in sympatric species of coral reef fish.

Coral reefs are highly diverse ecosystems, where numerous closely related species often coexist. How new species arise and are maintained in these high geneflow environments have been long-standing conundrums. Hybridization and patterns of introgression between sympatric species provide a unique insight into the mechanisms of speciation and the maintenance of species boundaries. In this study, we investigate the extent of hybridization between two closely related species of coral reef fish: the common coral trout (Plectropomus leopardus) and the bar-cheek coral trout (Plectropomus maculatus). Using a complementary set of 25 microsatellite loci, we distinguish pure genotype classes from first- and later-generation hybrids, identifying 124 interspecific hybrids from a collection of 2,991 coral trout sampled in inshore and mid-shelf reefs of the southern Great Barrier Reef. Hybrids were ubiquitous among reefs, fertile and spanned multiple generations suggesting both ecological and evolutionary processes are acting to maintain species barriers. We elaborate on these finding to investigate the extent of genomic introgression and admixture from 2,271 SNP loci recovered from a ddRAD library of pure and hybrid individuals. An analysis of genomic clines on recovered loci indicates that 261 SNP loci deviate from a model of neutral introgression, of which 132 indicate a pattern of introgression consistent with selection favouring both hybrid and parental genotypes. Our findings indicate genome-wide, bidirectional introgression between two sympatric species of coral reef fishes and provide further support to a growing body of evidence for the role of hybridization in the evolution of coral reef fishes.

The role of marine reserves in the replenishment of a locally impacted population of anemonefish on the Great Barrier Reef.

The development of parentage analysis to track the dispersal of juvenile offspring has given us unprecedented insight into the population dynamics of coral reef fishes. These tools now have the potential to inform fisheries management and species conservation, particularly for small fragmented populations under threat from exploitation and disturbance. In this study, we resolve patterns of larval dispersal for a population of the anemonefish Amphiprion melanopus in the Keppel Islands (southern Great Barrier Reef). Habitat loss and fishing appear to have impacted this population and a network of no-take marine reserves currently protects 75% of the potential breeders. Using parentage analysis, we estimate that 21% of recruitment in the island group was generated locally and that breeding adults living in reserves were responsible for 79% (31 of 39) of these of locally produced juveniles. Overall, the network of reserves was fully connected via larval dispersal; however, one reserve was identified as a critical source of larvae for the island group. The population in the Keppel Islands also appears to be well-connected to other source populations at least 60 km away, given that 79% (145 of 184) of the juveniles sampled remained unassigned in the parentage analysis. We estimated the effective size of the A. melanopus metapopulation to be 745 (582-993 95% CI) and recommend continued monitoring of its genetic status. Maintaining connectivity with populations beyond the Keppel Islands and recovery of local recruitment habitat, potentially through active restoration of host anemone populations, will be important for its long-term persistence.

Large-scale, multidirectional larval connectivity among coral reef fish populations in the Great Barrier Reef Marine Park.

Larval dispersal is the key process by which populations of most marine fishes and invertebrates are connected and replenished. Advances in larval tagging and genetics have enhanced our capacity to track larval dispersal, assess scales of population connectivity, and quantify larval exchange among no-take marine reserves and fished areas. Recent studies have found that reserves can be a significant source of recruits for populations up to 40 km away, but the scale and direction of larval connectivity across larger seascapes remain unknown. Here, we apply genetic parentage analysis to investigate larval dispersal patterns for two exploited coral reef groupers (Plectropomus maculatus and Plectropomus leopardus) within and among three clusters of reefs separated by 60-220 km within the Great Barrier Reef Marine Park, Australia. A total of 69 juvenile P. maculatus and 17 juvenile P. leopardus (representing 6% and 9% of the total juveniles sampled, respectively) were genetically assigned to parent individuals on reefs within the study area. We identified both short-distance larval dispersal within regions (200 m to 50 km) and long-distance, multidirectional dispersal of up to ~250 km among regions. Dispersal strength declined significantly with distance, with best-fit dispersal kernels estimating median dispersal distances of ~110 km for P. maculatus and ~190 km for P. leopardus. Larval exchange among reefs demonstrates that established reserves form a highly connected network and contribute larvae for the replenishment of fished reefs at multiple spatial scales. Our findings highlight the potential for long-distance dispersal in an important group of reef fishes, and provide further evidence that effectively protected reserves can yield recruitment and sustainability benefits for exploited fish populations.

Seascape and life-history traits do not predict self-recruitment in a coral reef fish.

The persistence and resilience of many coral reef species are dependent on rates of connectivity among sub-populations. However, despite increasing research efforts, the spatial scale of larval dispersal remains unpredictable for most marine metapopulations. Here, we assess patterns of larval dispersal in the angelfish Centropyge bicolor in Kimbe Bay, Papua New Guinea, using parentage and sibling reconstruction analyses based on 23 microsatellite DNA loci. We found that, contrary to previous findings in this system, self-recruitment (SR) was virtually absent at both the reef (0.4-0.5% at 0.15 km(2)) and the lagoon scale (0.6-0.8% at approx. 700 km(2)). While approximately 25% of the collected juveniles were identified as potential siblings, the majority of sibling pairs were sampled from separate reefs. Integrating our findings with earlier research from the same system suggests that geographical setting and life-history traits alone are not suitable predictors of SR and that high levels of localized recruitment are not universal in coral reef fishes.

Coral reef fish populations can persist without immigration.

Determining the conditions under which populations may persist requires accurate estimates of demographic parameters, including immigration, local reproductive success, and mortality rates. In marine populations, empirical estimates of these parameters are rare, due at least in part to the pelagic dispersal stage common to most marine organisms. Here, we evaluate population persistence and turnover for a population of orange clownfish, Amphiprion percula, at Kimbe Island in Papua New Guinea. All fish in the population were sampled and genotyped on five occasions at 2-year intervals spanning eight years. The genetic data enabled estimates of reproductive success retained in the same population (reproductive success to self-recruitment), reproductive success exported to other subpopulations (reproductive success to local connectivity), and immigration and mortality rates of sub-adults and adults. Approximately 50% of the recruits were assigned to parents from the Kimbe Island population and this was stable through the sampling period. Stability in the proportion of local and immigrant settlers is likely due to: low annual mortality rates and stable egg production rates, and the short larval stages and sensory capacities of reef fish larvae. Biannual mortality rates ranged from 0.09 to 0.55 and varied significantly spatially. We used these data to parametrize a model that estimated the probability of the Kimbe Island population persisting in the absence of immigration. The Kimbe Island population was found to persist without significant immigration. Model results suggest the island population persists because the largest of the subpopulations are maintained due to having low mortality and high self-recruitment rates. Our results enable managers to appropriately target and scale actions to maximize persistence likelihood as disturbance frequencies increase.

Seascape genetics along environmental gradients in the Arabian Peninsula: insights from ddRAD sequencing of anemonefishes.

Understanding the processes that shape patterns of genetic structure across space is a central aim of landscape genetics. However, it remains unclear how geographical features and environmental variables shape gene flow, particularly for marine species in large complex seascapes. Here, we evaluated the genomic composition of the two-band anemonefish Amphiprion bicinctus across its entire geographical range in the Red Sea and Gulf of Aden, as well as its close relative, Amphiprion omanensis endemic to the southern coast of Oman. Both the Red Sea and the Arabian Sea are complex and environmentally heterogeneous marine systems that provide an ideal scenario to address these questions. Our findings confirm the presence of two genetic clusters previously reported for A. bicinctus in the Red Sea. Genetic structure analyses suggest a complex seascape configuration, with evidence of both isolation by distance (IBD) and isolation by environment (IBE). In addition to IBD and IBE, genetic structure among sites was best explained when two barriers to gene flow were also accounted for. One of these coincides with a strong oligotrophic-eutrophic gradient at around 16-20˚N in the Red Sea. The other agrees with a historical bathymetric barrier at the straight of Bab al Mandab. Finally, these data support the presence of interspecific hybrids at an intermediate suture zone at Socotra and indicate complex patterns of genomic admixture in the Gulf of Aden with evidence of introgression between species. Our findings highlight the power of recent genomic approaches to resolve subtle patterns of gene flow in marine seascapes.

Linking local retention, self-recruitment, and persistence in marine metapopulations.

Three indices of larval retention have been used in the literature to assess the tendency for self-maintenance of local marine populations: local retention (LR), self-recruitment (SR), and relative local retention (RLR). Only one of these, LR, defined as the ratio of locally produced settlement to local egg production, has a clear relationship to self-persistence of individual sites. However, SR, the ratio of locally produced settlement to settlement of all origins at a site, is generally easier to measure experimentally. We use theoretical, simulation, and empirical approaches to bridge the gap between these different indices, and demonstrate that there is a proportional relationship between SR and LR, for metapopulations close to a stable state and with lifetime egg production (LEP) approximately uniform over space. Similarly, for systems where larval mortality rates are a relatively uniform function of release site, RLR (defined as the ratio of locally produced settlement to all settlement of local origin) and LR will also be proportional. Therefore, SR and RLR provide information on relative rates of LR for systems satisfying these conditions. Furthermore, the ratio between LR and SR can be used to evaluate global persistence of metapopulations, and therefore provides valuable information not necessarily available if only LR is considered.

Environmental gradients predict the genetic population structure of a coral reef fish in the Red Sea.

The relatively recent fields of terrestrial landscape and marine seascape genetics seek to identify the influence of biophysical habitat features on the spatial genetic structure of populations or individuals. Over the last few years, there has been accumulating evidence for the effect of environmental heterogeneity on patterns of gene flow and connectivity in marine systems. Here, we investigate the population genetic patterns of an anemonefish, Amphiprion bicinctus, along the Saudi Arabian coast of the Red Sea. We collected nearly one thousand samples from 19 locations, spanning approximately 1500 km, and genotyped them at 38 microsatellite loci. Patterns of gene flow appeared to follow a stepping-stone model along the northern and central Red Sea, which was disrupted by a distinct genetic break at a latitude of approximately 19°N. The Red Sea is characterized by pronounced environmental gradients along its axis, roughly separating the northern and central from the southern basin. Using mean chlorophyll-a concentrations as a proxy for this gradient, we ran tests of isolation by distance (IBD, R(2) = 0.52) and isolation by environment (IBE, R(2) = 0.64), as well as combined models using partial Mantel tests and multiple matrix regression with randomization (MMRR). We found that genetic structure across our sampling sites may be best explained by a combined model of IBD and IBE (Mantel: R(2) = 0.71, MMRR: R(2) = 0.86). Our results highlight the potential key role of environmental patchiness in shaping patterns of gene flow in species with pelagic larval dispersal. We support growing calls for the integration of biophysical habitat characteristics into future studies of population genetic structure.

Environmental DNA reveals temporal variation in mesophotic reefs of the Humboldt upwelling ecosystems of central Chile: Toward a baseline for biodiversity monitoring of unexplored marine habitats.

Temperate mesophotic reef ecosystems (TMREs) are among the least known marine habitats. Information on their diversity and ecology is geographically and temporally scarce, especially in highly productive large upwelling ecosystems. Lack of information remains an obstacle to understanding the importance of TMREs as habitats, biodiversity reservoirs and their connections with better-studied shallow reefs. Here, we use environmental DNA (eDNA) from water samples to characterize the community composition of TMREs on the central Chilean coast, generating the first baseline for monitoring the biodiversity of these habitats. We analyzed samples from two depths (30 and 60 m) over four seasons (spring, summer, autumn, and winter) and at two locations approximately 16 km apart. We used a panel of three metabarcodes, two that target all eukaryotes (18S rRNA and mitochondrial COI) and one specifically targeting fishes (16S rRNA). All panels combined encompassed eDNA assigned to 42 phyla, 90 classes, 237 orders, and 402 families. The highest family richness was found for the phyla Arthropoda, Bacillariophyta, and Chordata. Overall, family richness was similar between depths but decreased during summer, a pattern consistent at both locations. Our results indicate that the structure (composition) of the mesophotic communities varied predominantly with seasons. We analyzed further the better-resolved fish assemblage and compared eDNA with other visual methods at the same locations and depths. We recovered eDNA from 19 genera of fish, six of these have also been observed on towed underwater videos, while 13 were unique to eDNA. We discuss the potential drivers of seasonal differences in community composition and richness. Our results suggest that eDNA can provide valuable insights for monitoring TMRE communities but highlight the necessity of completing reference DNA databases available for this region.