An important goal of marine population genetics is to understand how spatial connectivity patterns are influenced by historical and evolutionary factors. In this study, we evaluate the demographic history and population structure of Littoraria flava, a highly dispersive marine gastropod in the Brazilian intertidal zone. To test the hypotheses that the species has (1) historically high levels of gene flow on a macrogeographical spatial scale and (2) a distribution in rocky shores that consists of subpopulations, we collected specimens along the Brazilian coastline and combined different sets of genetic markers (mitochondrial DNA, ITS-2 and single nucleotide polymorphisms) with niche-based modelling to predict its palaeodistribution. Low genetic structure was observed, as well as high gene flow over long distances. The demographic analyses suggest that L. flava has had periods of population bottlenecks followed by expansion. According to both palaeodistribution and coalescent simulations, these expansion events occurred during the Pleistocene interglacial cycles (21 kya) and the associated climatic changes were the probable drivers of the distribution of the species. This is the first phylogeographical study of a marine gastropod on the South American coast based on genomic markers associated with niche modelling.

Understanding how selection shapes population differentiation and local adaptation in marine species remains one of the greatest challenges in the field of evolutionary biology. The selection of genes in response to environment-specific factors and microenvironmental variation often results in chaotic genetic patchiness, which is commonly observed in rocky shore organisms. To identify these genes, the expression profile of the marine gastropod Littoraria flava collected from four Southeast Brazilian locations in ten rocky shore sites was analyzed. In this first L. flava transcriptome, 250,641 unigenes were generated, and 24% returned hits after functional annotation. Independent paired comparisons between 1) transects, 2) sites within transects, and 3) sites from different transects were performed for differential expression, detecting 8,622 unique differentially expressed genes. Araçá (AR) and São João (SJ) transect comparisons showed the most divergent gene products. For local adaptation, fitness-related differentially expressed genes were chosen for selection tests. Nine and 24 genes under adaptative and purifying selection, respectively, were most related to biomineralization in AR and chaperones in SJ. The biomineralization-genes perlucin and gigasin-6 were positively selected exclusively in the site toward the open ocean in AR, with sequence variants leading to pronounced protein structure changes. Despite an intense gene flow among L. flava populations due to its planktonic larva, gene expression patterns within transects may be the result of selective pressures. Our findings represent the first step in understanding how microenvironmental genetic variation is maintained in rocky shore populations and the mechanisms underlying local adaptation in marine species.

Background: An important goal of population genetics studies in marine ecosystems is understanding how connectivity patterns, both spatial and temporal, are influenced by historical and evolutionary factors. When it comes to dispersion and connectivity in marine ecosystems, the role of historical and evolutionary factors over population dynamics andstructure still remains enigmatic. We evaluated the demographic history and population structure of Littoraria flava, a highly dispersive and widely distributed marine gastropod on the Brazilian intertidal zone, to predict the effects of such factors on intrapopulation divergence. To test the hypotheses that (1) the species has historically high levels of geneflow on a macrogeographic spatial scale and (2) the species distribution in rocky shores consists of subpopulations due to high degrees of environmental heterogeneity, we collected specimens along the Brazilian coastline and combined different sets of genetic markers (mitochondrial DNA, nuclear internal transcribed spacer 2, and single nucleotide polymorphisms) with niche-based modeling to predict species paleodistribution.Results: Low genetic structure was observed along the coastline, and all clustering and migration analyses supported the high gene flow over long distances hypothesis (> 3,000 km). Three genetic clusters were identified by the assignment test, each mostly composed of individuals from the three sampled regions. No fine-scale variation was observed for any location. The neutrality tests and the haplotype networks suggest that L. flava had experienced population bottleneck followed by population expansion. Both paleodistribution and coalescent simulations highlight that expansion events occurred in the Southeastern coastline during the Pleistocene interglacial cycles (21 kya).Conclusions: This is the first study on the South American coast that highlights the demographic history on a marine gastropod based on genomic markers associated with niche modelling. We found that climatic changes since the interglacial periods are potentially relevant drivers for the species distribution in the past. Our findings could enhance the understanding of the population dynamics under an evolutionary view for widely distributed non-model marine organisms.

Aim Seascape genomic studies aim to understand how environmental variables shape species diversity through genotype–environment associations. Identifying these effects on lecithotrophic larval species that live in intertidal zones is particularly challenging because they are subject to environmental heterogeneity and anthropogenic events. Here, we evaluate how biotic and abiotic features in the Southwest Atlantic littoral zone can affect a high dispersal species' present and historical demography. Location Brazilian coast, covering more than 3000 km. Taxon Perinereis ponteni. Methods We investigated population genetic diversity, connectivity and past dynamics using 23,300 single‐nucleotide polymorphisms (SNPs) generated using genotyping by sequencing. We tested whether environmental abiotic variables could explain the variance found in genotype frequencies using isolation‐by‐environment (IBE) and landscape association approaches. These data, combined with palaeodistribution simulations and oceanic circulation modelling, were used to infer species demographic history and connectivity patterns. Results Along with high levels of connectivity detected, we found a genetic boundary in the southeastern region of Brazil around Cabo Frio (Rio de Janeiro), and a cline trend for some loci. The palaeodistribution simulations reveal a spatial refuge in the southeast during the Last Glacial Maximum (21 kya), with the expansion of the northern region. We identified 1421 SNPs with frequencies associated with eight environmental variables, most of which were related to temperature—the main environmental factor determining IBE. Main Conclusions Perinereis ponteni, a polychaete with high gene flow capability responds to biogeographic barriers, highlighting the importance of biotic and abiotic factors in shaping population connectivity. Furthermore, the effect of temperature indicates that future climate change and ocean warming can hugely impact this species.

Evolutionary history leads to genome changes over time, especially for species that have experienced intense selective pressures over a short period. Here, we investigated the genomic evolution of Bos species by searching for potential selection signatures, focusing on Nelore, an economically relevant cattle breed in Brazil. We assessed the genomic processes determining the molecular evolution across Nelore and thirteen other related taxa by evaluating (i) amino acid sequence conservation, (ii) the dN/dS ratio, and (iii) gene families’ turnover rate (λ). Low conserved regions potentially associated with fatty acid metabolism seem to reflect differences in meat fat content in taxa with different evolutionary histories. All Bos species presented genes under positive selection, especially B. indicus and Nelore, which include transport protein cobalamin, glycolipid metabolism, and hormone signaling. These findings could be explained by constant selective pressures to obtain higher immune resistance and efficient metabolism. The gene contraction rate across the Nelore + B. indicus branch was almost nine times higher than that in other lineages (λ = 0.01043 vs. 0.00121), indicating gene losses during the domestication process. Amino acid biosynthesis, reproductive and innate immune system-related pathways were associated with genes recognized within the most frequent rapidly evolving gene families and in genes under positive selection, supporting the substantial relevance of such traits from a domestication perspective. Our data provide new insights into how the genome may respond to intense artificial selection in distinct taxa, and reinforces the presence of selective pressures on traits potentially relevant for future animal breeding investments.