ABSTRACT
Knowledge of genetic structure at the finest level is essential for the conservation of genetic resources. Despite no visible barriers limiting gene flow, significant genetic structure has been shown in marine species. The common cockle (Cerastoderma edule) is a bivalve of great commercial and ecological value inhabiting the Northeast Atlantic Ocean. Previous population genomics studies demonstrated significant structure both across the Northeast Atlantic, but also within small geographic areas, highlighting the need to investigate fine-scale structuring. Here, we analysed two geographic areas that could represent opposite models of structure for the species: (1) the SW British Isles region, highly fragmented due to biogeographic barriers, and (2) Galicia (NW Spain), a putative homogeneous region. A total of 9250 SNPs genotyped by 2b-RAD on 599 individuals from 22 natural beds were used for the analysis. The entire SNP dataset mostly confirmed previous observations related to genetic diversity and differentiation; however, neutral and divergent SNP outlier datasets enabled disentangling physical barriers from abiotic environmental factors structuring both regions. While Galicia showed a homogeneous structure, the SW British Isles region was split into four reliable genetic regions related to oceanographic features and abiotic factors, such as sea surface salinity and temperature. The information gathered supports specific management policies of cockle resources in SW British and Galician regions also considering their particular socio-economic characteristics; further, these new data will be added to those recently reported in the Northeast Atlantic to define sustainable management actions across the whole distribution range of the species.
Subject(s)
Cardiidae , Humans , Animals , Atlantic Ocean , Spain , Genotype , Genetic StructuresABSTRACT
Multiple parasites can infect a single host, creating a dynamic environment where each parasite must compete over host resources. Such interactions can cause greater harm to the host than single infections and can also have negative consequences for the parasites themselves. In their first intermediate hosts, trematodes multiply asexually and can eventually reach up to 20% of the host's biomass. In most species, it is unclear whether this biomass results from a single infection or co-infection by 2 or more infective stages (miracidia), the latter being more likely a priori in areas where prevalence of infection is high. Using as model system the trematode Bucephalus minimus and its first intermediate host cockles, we examined the genetic diversity of the cytochrome c oxidase subunit I region in B. minimus from 3 distinct geographical areas and performed a phylogeographic study of B. minimus populations along the Northeast Atlantic coast. Within localities, the high genetic variability found across trematodes infecting different individual cockles, compared to the absence of variability within the same host, suggests that infections could be generally originating from a single miracidium. On a large spatial scale, we uncovered significant population structure of B. minimus, specifically between the north and south of Bay of Biscay. Although other explanations are possible, we suggest this pattern may be driven by the population structure of the final host.
Subject(s)
Cardiidae , Parasites , Trematoda , Animals , Host-Parasite Interactions , Trematoda/genetics , Phylogeography , Cardiidae/parasitologyABSTRACT
Raphidascarid nematodes have been the focus of several studies, mainly due to the zoonotic potential of some species, even though the cases are underreported. Due to the difficulty in identifying their larvae, the use of diagnostic techniques involving morphological and molecular analyses has grown in the last 20 years. The present study had as objective the morphological and molecular characterization of the L3 larval types of Hysterothylacium collected in Pomatomus saltatrix and Pagrus pagrus from the Brazilian coast, close to the municipality of Santos, State of São Paulo. Twenty specimens of P. saltatrix were necropsied and Hysterothylacium type V (n = 257) and Hysterothylacium type X (n = 5) larvae were found. Five specimens of P. pagrus were necropsied and all were parasitized by Hysterothylacium type V larvae. The analyses showed a genetic proximity relationship between Hysterothylacium types V with other Hysterothylacium V and with H. deardorffoverstreetorum, although this is a species inquirenda. Haplotypes for Hysterothylacium type X found in the present study formed a monophyletic group with other Hysterothylacium X, H. amoyense, and H. zhoushanense. Through this study, new hosts and localities were registered for Hysterothylacium type V and Hysterothylacium type X.
Subject(s)
Ascaridoidea , Fish Diseases , Perciformes , Animals , Brazil , Ascaridoidea/anatomy & histology , Ascaridoidea/genetics , Larva/genetics , FishesABSTRACT
BACKGROUND: The irruption of Next-generation sequencing (NGS) and restriction site-associated DNA sequencing (RAD-seq) in the last decade has led to the identification of thousands of molecular markers and their genotyping for refined genomic screening. This approach has been especially useful for non-model organisms with limited genomic resources. Many building-loci pipelines have been developed to obtain robust single nucleotide polymorphism (SNPs) genotyping datasets using a de novo RAD-seq approach, i.e. without reference genomes. Here, the performances of two building-loci pipelines, STACKS 2 and Meyer's 2b-RAD v2.1 pipeline, were compared using a diverse set of aquatic species representing different genomic and/or population structure scenarios. Two bivalve species (Manila clam and common edible cockle) and three fish species (brown trout, silver catfish and small-spotted catshark) were studied. Four SNP panels were evaluated in each species to test both different building-loci pipelines and criteria for SNP selection. Furthermore, for Manila clam and brown trout, a reference genome approach was used as control. RESULTS: Despite different outcomes were observed between pipelines and species with the diverse SNP calling and filtering steps tested, no remarkable differences were found on genetic diversity and differentiation within species with the SNP panels obtained with a de novo approach. The main differences were found in brown trout between the de novo and reference genome approaches. Genotyped vs missing data mismatches were the main genotyping difference detected between the two building-loci pipelines or between the de novo and reference genome comparisons. CONCLUSIONS: Tested building-loci pipelines for selection of SNP panels seem to have low influence on population genetics inference across the diverse case-study scenarios here studied. However, preliminary trials with different bioinformatic pipelines are suggested to evaluate their influence on population parameters according with the specific goals of each study.
Subject(s)
Metagenomics , Polymorphism, Single Nucleotide , Animals , Benchmarking , Genome , Sequence Analysis, DNAABSTRACT
Among several marine pelagic species of the Brazilian coast, Scomber colias Gmelin, 1789 (Perciformes: Scombridae) stands out for having great economic importance, since it is widely used as a food resource and presents moderate vulnerability. Twenty specimens of S. colias were purchased from October 2015 to October 2016 from the coast of Santa Catarina, Brazil. In the present study, we recorded Hysterothylacium fabri (Rudolphi, 1819) (Deardorff and Overstreet, Proc Biol Soc Wash 93(4):1035-1079 1980) from the S. colias intestine using an integrative taxonomy approach, where morphological data are used in combination with partial sequences of the ITS gene, to validate the taxonomic status of the species and establish their relationships with other members of the genus. This species is being recorded for the first time in the South Atlantic and S. colias. The specimens of H. fabri collected in this study parasitizing S. colias presented morphology similar to the other specimens already registered parasitizing other hosts. The distance matrix generated showed that the partial sequences obtained in this study were more similar to sequences of Hysterothylacium sp. collected in China. In phylogenetic analysis, the two detected haplotypes of this study were grouped with H. fabri haplotypes deposited in GenBank in a monophyletic subclade.
Subject(s)
Ascaridoidea/isolation & purification , Perciformes/parasitology , Animals , Ascaridoidea/classification , Ascaridoidea/ultrastructure , Atlantic Ocean , Perciformes/classification , PhylogenyABSTRACT
Several methods have been described to quantify the first ray mobility. They all have certain disadvantages (great size, sophistication, or lack of validation). The objective of this work was to study the validity and reliability of a new instrument for the measurement of first ray mobility. Anterior-posterior radiographs were obtained from 25 normal feet and 24 hallux valgus feet, with the first ray in a neutral position, maximally dorsiflexed and maximally plantarflexed. The first ray mobility was radiographicaly measured in both groups, and was also manually examined with the new device. A cluster analysis determined whether normal and hallux valgus feet were correctly classified, and a graphic analysis of Bland-Altman was performed to compare the radiographic and manual measurement techniques. Based on the radiographs, the first ray mobility only showed significant differences in dorsiflexion between both groups (P = 0.015). First ray dorsiflexion, plantarflexion and total range of motion measured with the new device were different between both groups (P = 0.040, P = 0.011 and P = 0.006, respectively). The silhouette measure of the cohesion and separation coefficients from the cluster analysis was greater than 0.50 for the dorsiflexion, plantarflexion and total range of motion obtained from the radiographs and from the new device. The Bland-Altman graph suggested that 96% of the data presented agreement between both measurement methods. These results suggested that the new instrument was valid and reliable.
ABSTRACT
BACKGROUND: The superficial lymphatic component of vascular malformations poses a significant treatment challenge. It is responsible for the majority of symptoms presented, and to date, there is no consensus regarding treatment. OBJECTIVE: To evaluate the effectiveness of topical rapamycin in treating superficial lymphatic malformations (LM). METHODS: A case series study was performed of patients with superficial LM, treated with topical rapamycin. The clinical characteristics of patients and the concentration and application mode of the drug were recorded. The changes in the signs and symptoms observed and associated adverse effects were noted and analyzed. RESULTS: The study population consisted of 11 patients of an average age of 10.5 years. All were treated with topical rapamycin: 6 patients with a 1% concentration, 1 with a 0.8% concentration, and 4 with a 0.4% concentration. Changes in the clinical appearance of the lesions were observed in all patients. The associated symptoms, present in 9 of 11 patients, improved in every case. The mean follow-up time was 16.1 months. LIMITATIONS: This study is retrospective, with a small sample size and considerable heterogeneity of lesions and treatment approaches. CONCLUSION: Treatment with topical rapamycin modifies the clinical appearance and alleviates symptoms of superficial LM.
Subject(s)
Immunosuppressive Agents/therapeutic use , Lymphatic Abnormalities/diagnosis , Lymphatic Abnormalities/drug therapy , Sirolimus/therapeutic use , Administration, Topical , Adolescent , Child , Child, Preschool , Cohort Studies , Female , Humans , Male , Prognosis , Retrospective Studies , Severity of Illness Index , Treatment OutcomeABSTRACT
Biological invasions rank among the main global threats for biodiversity. The Eastern mosquitofish (Gambusia holbrooki) is considered one of the 100 world worst invasive species due to its high adaptation capability to new environments. Using the restriction-site-associated DNA tags (RADtags), introduced European locations were compared against native US mosquitofish populations to analyse genomic changes that occurred during invasive process of European locations. After filtering, 7724 RADtags containing only one SNP were retained for population studies. Comparative genomics indicated that 186 of these RADtags matched sequences in the transcriptome of Xyphophorus maculatus, the most closely related genome available. Genomic analyses showed that invasive populations show high reductions in diversity. Further, analyses of population structuring based on these data are concordant with previous analyses based on microsatellites. It is concluded that during the invasion process genetic drift was the main evolutionary force affecting patterns of diversity and population structure. While recognizing that positive selection could be masked by the strong drift during founder events, adaptive processes were evidenced in a reduced number of RADtags (<2%), with only one of these in a putative coding region. Surprisingly, balancing selection was detected in several coding RADtags, suggesting that the preservation of polymorphism in specific genes could be more important than the average population diversity for the population maintenance at any location, particularly for the survival of introduced populations.
Subject(s)
Cyprinodontiformes/genetics , Evolution, Molecular , Genetic Drift , Genetics, Population , Introduced Species , Animals , Genetic Variation , Genomics , Italy , Microsatellite Repeats , Polymorphism, Single Nucleotide , Selection, Genetic , Sequence Analysis, DNA , Spain , Transcriptome , United StatesABSTRACT
While the Philippines aspires to be one of the top tourist destinations in Southeast Asia, self-inflicted wounds like the failure of the government to comply with international aviation safety standards may derail the country from achieving its goals. This article estimates the short- and long-term impact of the US FAA downgrade of the Philippine civil aviation system in 2008 and the EU ban of Philippine carriers in 2010 on tourist expenditures, arrivals, and length of stay using monthly time series data. The econometric model, consisting of three equations due to the endogeneity of the tourist arrivals and length of stay variables in the tourist expenditures equation, is estimated simultaneously using the generalized method of moments. The results indicate that the US FAA downgrade and the EU ban impact monthly tourist receipts negatively in the short term while the downgrade also impacts tourist expenditures in the long term. Moreover, the ban impacts length of stay negatively in the short and long term while the downgrade impacts length of stay negatively only in the long term. The substantial decline in tourism receipts from 2008 to 2010 despite an increasing trend in tourist arrivals is due to the shorter stay of tourists, indicating that high-spending tourists have not returned following the downgrade and ban.
ABSTRACT
BACKGROUND: The microbial landscape within termite guts varies across termite families. The gut microbiota of lower termites (LT) is dominated by cellulolytic flagellates that sequester wood particles in their digestive vacuoles, whereas in the flagellate-free higher termites (HT), cellulolytic activity has been attributed to fiber-associated bacteria. However, little is known about the role of individual lineages in fiber digestion, particularly in LT. RESULTS: We investigated the lignocellulolytic potential of 2223 metagenome-assembled genomes (MAGs) recovered from the gut metagenomes of 51 termite species. In the flagellate-dependent LT, cellulolytic enzymes are restricted to MAGs of Bacteroidota (Dysgonomonadaceae, Tannerellaceae, Bacteroidaceae, Azobacteroidaceae) and Spirochaetota (Breznakiellaceae) and reflect a specialization on cellodextrins, whereas their hemicellulolytic arsenal features activities on xylans and diverse heteropolymers. By contrast, the MAGs derived from flagellate-free HT possess a comprehensive arsenal of exo- and endoglucanases that resembles that of termite gut flagellates, underlining that Fibrobacterota and Spirochaetota occupy the cellulolytic niche that became vacant after the loss of the flagellates. Furthermore, we detected directly or indirectly oxygen-dependent enzymes that oxidize cellulose or modify lignin in MAGs of Pseudomonadota (Burkholderiales, Pseudomonadales) and Actinomycetota (Actinomycetales, Mycobacteriales), representing lineages located at the hindgut wall. CONCLUSIONS: The results of this study refine our concept of symbiotic digestion of lignocellulose in termite guts, emphasizing the differential roles of specific bacterial lineages in both flagellate-dependent and flagellate-independent breakdown of cellulose and hemicelluloses, as well as a so far unappreciated role of oxygen in the depolymerization of plant fiber and lignin in the microoxic periphery during gut passage in HT. Video Abstract.
Subject(s)
Gastrointestinal Microbiome , Isoptera , Lignin , Animals , Isoptera/microbiology , Lignin/metabolism , Phylogeny , Bacteria/genetics , Bacteria/classification , Bacteria/isolation & purification , Bacteria/enzymology , Metagenome , Polysaccharides/metabolism , Genome, Bacterial , Genomics , Cellulose/metabolismABSTRACT
The turbot (Scophthalmus maximus) is a commercially valuable flatfish and one of the most promising aquaculture species in Europe. Two transcriptome 454-pyrosequencing runs were used in order to detect Single Nucleotide Polymorphisms (SNPs) in genes related to immune response and gonad differentiation. A total of 866 true SNPs were detected in 140 different contigs representing 262,093 bp as a whole. Only one true SNP was analyzed in each contig. One hundred and thirteen SNPs out of the 140 analyzed were feasible (genotyped), while Ш were polymorphic in a wild population. Transition/transversion ratio (1.354) was similar to that observed in other fish studies. Unbiased gene diversity (He) estimates ranged from 0.060 to 0.510 (mean = 0.351), minimum allele frequency (MAF) from 0.030 to 0.500 (mean = 0.259) and all loci were in Hardy-Weinberg equilibrium after Bonferroni correction. A large number of SNPs (49) were located in the coding region, 33 representing synonymous and 16 non-synonymous changes. Most SNP-containing genes were related to immune response and gonad differentiation processes, and could be candidates for functional changes leading to phenotypic changes. These markers will be useful for population screening to look for adaptive variation in wild and domestic turbot.
ABSTRACT
The European brown trout, Salmo trutta, is a cold-adapted fish reported as a Least Concern species in the IUCN Red List. This species colonized new territories from southern refuges during the last glacial melting, but during the 20th century suffered from anthropic impacts on its habitats. The long-time survival of the species relies on the genetic diversity within and among populations. Brown trout is among the genetically most diverse vertebrate species; however, native populations in Mediterranean rivers have dramatically suffered of introgressive hybridization from extensive releases of evolutionary distant non-native Atlantic stocks. In addition, in Mediterranean rivers climate change will result in unsuitable conditions for the species during the 21st century. Using brown trout populations at the headstreams of a Pyrenean river as a model, this paper revised how hatchery releases have affected the native gene pools and how environmental and climatic variables controlled the amount of local introgression at intra-basin level. Introgressive hybridization was detected in all studied sites. Ten times larger divergence was observed among populations at tributaries than among populations along the main stem. A highly impacted population distributed in a long transect in the main stem suggested that hatchery fish move towards the main stem wherever released. From already highly impacted populations and despite the cessation of hatchery releases, warmer temperatures and lower precipitation expected from climate change will extend the introgressive hybridization along the basin, contributing to the extinction of the native gene pools. Based on available morphological distinction of native, hatchery and hybrid brown trout, we advocate the involvement of regional social groups (e.g. riverside dwellers, anglers, conservationists, hikers) in citizen science programs to detect the spread of non-native phenotypes along the rivers. These are cheap and fast methods to collaborate with fishery managers in the preservation and recovery of the regional native populations.
Subject(s)
Gene Pool , Rivers , Animals , Trout/genetics , Ecosystem , Hong KongABSTRACT
The common cockle is a valuable bivalve species inhabiting the Atlantic European coasts. The parasite Marteilia cochillia has devastated cockle beds in the southern Galician (NW Spain) rias since 2012. Previous data suggested that cockles from Ría de Arousa acquired some resilience to this parasite through natural selection after consecutive annual marteiliosis outbreaks and candidate markers associated with marteiliosis resilience were identified using population genomics and transcriptomics approaches. Here, a common garden experiment was performed using a naïve stock (from Ría de Muros-Noia) and an affected stock (from Ría de Arousa) to test this hypothesis. Breeders from both stocks were used to produce seed cohorts at hatchery, which were pre-grown in a raft (outdoor nursery stage) and deployed in two shellfish beds affected by marteiliosis in Ría de Arousa (growing-out stage). In both beds, the naïve stock showed high marteiliosis prevalence and was fully depleted in a short period, while the affected stock barely showed evidence of marteiliosis. A set of 45 SNPs putatively associated with marteiliosis resilience were fitted for MassARRAY genotyping to check their role in the differential resilience detected between both stocks. Though no significant differentiation was found between the naïve and the affected stocks with neutral markers, 28 SNPs showed significant divergence between them, suggesting that these SNPs were involved in directional selection during eight generations (to the most) of marteiliosis pressure (long-term selection). Furthermore, signals of selection were also detected in the naïve stock along the marteiliosis outbreak in the growing-out stage (short-term selection) and six SNPs, all shared with the long-term evaluation, showed consistent signals of differentiation according to the infection severity. Some of these SNPs were located within immune genes pertaining to families such as proteasome, ubiquitin, tumor necrosis factor, and glutathione S-transferase. These resilience-associated markers will be useful to recover cockle production in Galicia.
ABSTRACT
BACKGROUND: The turbot (Scophthalmus maximus) is a relevant species in European aquaculture. The small turbot genome provides a source for genomics strategies to use in order to understand the genetic basis of productive traits, particularly those related to sex, growth and pathogen resistance. Genetic maps represent essential genomic screening tools allowing to localize quantitative trait loci (QTL) and to identify candidate genes through comparative mapping. This information is the backbone to develop marker-assisted selection (MAS) programs in aquaculture. Expressed sequenced tag (EST) resources have largely increased in turbot, thus supplying numerous type I markers suitable for extending the previous linkage map, which was mostly based on anonymous loci. The aim of this study was to construct a higher-resolution turbot genetic map using EST-linked markers, which will turn out to be useful for comparative mapping studies. RESULTS: A consensus gene-enriched genetic map of the turbot was constructed using 463 SNP and microsatellite markers in nine reference families. This map contains 438 markers, 180 EST-linked, clustered at 24 linkage groups. Linkage and comparative genomics evidences suggested additional linkage group fusions toward the consolidation of turbot map according to karyotype information. The linkage map showed a total length of 1402.7 cM with low average intermarker distance (3.7 cM; ~2 Mb). A global 1.6:1 female-to-male recombination frequency (RF) ratio was observed, although largely variable among linkage groups and chromosome regions. Comparative sequence analysis revealed large macrosyntenic patterns against model teleost genomes, significant hits decreasing from stickleback (54%) to zebrafish (20%). Comparative mapping supported particular chromosome rearrangements within Acanthopterygii and aided to assign unallocated markers to specific turbot linkage groups. CONCLUSIONS: The new gene-enriched high-resolution turbot map represents a useful genomic tool for QTL identification, positional cloning strategies, and future genome assembling. This map showed large synteny conservation against model teleost genomes. Comparative genomics and data mining from landmarks will provide straightforward access to candidate genes, which will be the basis for genetic breeding programs and evolutionary studies in this species.
Subject(s)
Chromosome Mapping/methods , Expressed Sequence Tags , Fishes/genetics , Flatfishes/genetics , Animals , Female , Male , Microsatellite Repeats , Quantitative Trait Loci , SyntenyABSTRACT
In this study, we quantified the three key biological processes, growth, recruitment, and dispersal pattern, which are necessary for a better understanding of the population dynamics of the blue and red shrimp Aristeus antennatus. This marine exploited crustacean shows sex-related distribution along the water column, being females predominate in the middle slope. The present study attempts to fill the existing gap in the females' genetic demography, as scarce knowledge is available despite being the most abundant sex in catches. We analyzed morphometric data and genotyped 12 microsatellite loci in 665 A. antennatus females collected in two consecutive seasons, winter and summer 2016, at the main Mediterranean fishing ground as a model. Almost every female in summer was inseminated. Five modal groups were observed in both seasons, from 0+ to 4+ in winter and from 1+ to 5+ in summer. Commercial-sized sorting based on fishermen's experience resulted in a moderate-to-high assertive method concerning cohort determination. Genetic data pointed out females' horizontal movement between neighboring fishing grounds, explaining the low genetic divergence detected among western Mediterranean grounds. Our results could represent critical information for the future implementation of management measures to ensure long-time conservation of the A. antennatus populations.
Subject(s)
Decapoda , Penaeidae , Animals , Decapoda/genetics , Female , Genotype , Humans , Male , Microsatellite Repeats/genetics , Penaeidae/genetics , Population DynamicsABSTRACT
Brown trout (Salmo trutta L.) populations have been restocked during recent decades to satisfy angling demand and counterbalance the decline of wild populations. Millions of fertile brown trout individuals were released into Mediterranean and Atlantic rivers from hatcheries with homogeneous central European stocks. Consequently, many native gene pools have become endangered by introgressive hybridization with those hatchery stocks. Different genetic tools have been used to identify and evaluate the degree of introgression starting from pure native and restocking reference populations (e.g., LDH-C* locus, microsatellites). However, due to the high genetic structuring of brown trout, the definition of the "native pool" is hard to achieve. Additionally, although the LDH-C* locus is useful for determining the introgression degree at the population level, its consistency at individual level is far from being accurate, especially after several generations were since releases. Accordingly, the development of a more powerful and cost-effective tool is essential for an appropriate monitoring to recover brown-trout-native gene pools. Here, we used the 2b restriction site-associated DNA sequencing (2b-RADseq) and Stacks 2 with a reference genome to identify single-nucleotide polymorphisms (SNPs) diagnostic for hatchery-native fish discrimination in the Atlantic and Mediterranean drainages of the Iberian Peninsula. A final set of 20 SNPs was validated in a MassARRAY® System genotyping by contrasting data with the whole SNP dataset using samples with different degree of introgression from those previously recorded. Heterogeneous introgression impact was confirmed among and within river basins, and was the highest in the Mediterranean Slope. The SNP tool reported here should be assessed in a broader sample scenario in Southern Europe considering its potential for monitoring recovery plans.
Subject(s)
Polymorphism, Single Nucleotide , Rivers , Animals , Gene Pool , Humans , Microsatellite Repeats , Trout/geneticsABSTRACT
European flat oyster (Ostrea edulis) is an ecologically and economically important marine bivalve, that has been severely affected by the intracellular parasite Bonamia ostreae. In this study, a flat oyster SNP array (~14,000 SNPs) was used to validate previously reported outlier loci for divergent selection associated with B. ostreae exposure in the Northeast Atlantic Area. A total of 134 wild and hatchery individuals from the North Sea, collected in naïve (NV) and long-term affected (LTA) areas, were analysed. Genetic diversity and differentiation were related to the sampling origin (wild vs. hatchery) when using neutral markers, and to bonamiosis status (NV vs. LTA) when using outlier loci for divergent selection. Two genetic clusters appeared intermingled in all sampling locations when using outlier loci, and their frequency was associated with their bonamiosis status. When both clusters were compared, outlier data sets showed high genetic divergence (F ST > 0.25) unlike neutral loci (F ST not ≠ 0). Moreover, the cluster associated with LTA samples showed much higher genetic diversity and significant heterozygote excess with outlier loci, but not with neutral data. Most outliers mapped on chromosome 8 (OE-C8) of the flat oyster genome, supporting a main genomic region underlying resilience to bonamiosis. Furthermore, differentially expressed genes previously reported between NV and LTA strains showed higher mapping density on OE-C8. A range of relevant immune functions were specifically enriched among genes annotated on OE-C8, providing hypotheses for resilience mechanisms to an intracellular parasite. The results suggest that marker-assisted selection could be applied to breed resilient strains of O. edulis to bonamiosis, if lower parasite load and/or higher viability of the LTA genetic cluster following B. ostreae infection is demonstrated.
ABSTRACT
Knowledge on correlations between environmental factors and genome divergence between populations of marine species is crucial for sustainable management of fisheries and wild populations. The edible cockle (Cerastoderma edule) is a marine bivalve distributed along the Northeast Atlantic coast of Europe and is an important resource from both commercial and ecological perspectives. We performed a population genomics screening using 2b-RAD genotyping on 9309 SNPs localized in the cockle's genome on a sample of 536 specimens pertaining to 14 beds in the Northeast Atlantic Ocean to analyse the genetic structure with regard to environmental variables. Larval dispersal modelling considering species behaviour and interannual/interseasonal variation in ocean conditions was carried out as an essential background to which compare genetic information. Cockle populations in the Northeast Atlantic displayed low but significant geographical differentiation between populations (F ST = 0.0240; p < 0.001), albeit not across generations. We identified 742 and 36 outlier SNPs related to divergent and balancing selection in all the geographical scenarios inspected, and sea temperature and salinity were the main environmental correlates suggested. Highly significant linkage disequilibrium was detected at specific genomic regions against the very low values observed across the whole genome. Two main genetic groups were identified, northwards and southwards of French Brittany. Larval dispersal modelling suggested a barrier for larval dispersal linked to the Ushant front that could explain these two genetic clusters. Further genetic subdivision was observed using outlier loci and considering larval advection. The northern group was divided into the Irish/Celtic Seas and the English Channel/North Sea, while the southern group was divided into three subgroups. This information represents the baseline for the management of cockles, designing conservation strategies, founding broodstock for depleted beds and producing suitable seed for aquaculture production.
ABSTRACT
Tambaqui or cachama (Colossoma macropomum) is one of the most important neotropical freshwater fish used for aquaculture in South America, and its production is concentrated at low latitudes (close to the Equator, 0°), where the water temperature is warm. Therefore, understanding how selection shapes genetic variations and structure in farmed populations is of paramount importance in evolutionary biology. High-throughput sequencing to generate genome-wide data for fish species allows for elucidating the genomic basis of adaptation to local or farmed conditions and uncovering genes that control the phenotypes of interest. The present study aimed to detect genomic selection signatures and analyze the genetic variability in farmed populations of tambaqui in South America using single-nucleotide polymorphism (SNP) markers obtained with double-digest restriction site-associated DNA sequencing. Initially, 199 samples of tambaqui farmed populations from different locations (located in Brazil, Colombia, and Peru), a wild population (Amazon River, Brazil), and the base population of a breeding program (Aquaculture Center, CAUNESP, Jaboticabal, SP, Brazil) were genotyped. Observed and expected heterozygosity was 0.231-0.350 and 0.288-0.360, respectively. Significant genetic differentiation was observed using global FST analyses of SNP loci (FST = 0.064, p < 0.050). Farmed populations from Colombia and Peru that differentiated from the Brazilian populations formed distinct groups. Several regions, particularly those harboring the genes of significance to aquaculture, were identified to be under positive selection, suggesting local adaptation to stress under different farming conditions and management practices. Studies aimed at improving the knowledge of genomics of tambaqui farmed populations are essential for aquaculture to gain deeper insights into the evolutionary history of these fish and provide resources for the establishment of breeding programs.