A comprehensive coding and microRNA transcriptome of vertebral bone in postlarvae and juveniles of Senegalese sole (Solea senegalensis).

Understanding vertebral bone development is essential to prevent skeletal malformations in farmed fish related to genetic and environmental factors. This is an important issue in Solea senegalensis, with special impact of spinal anomalies in postlarval and juvenile stages. Vertebral bone transcriptomics in farmed fish mainly comes from coding genes, and barely on miRNA expression. Here, we used RNA-seq of spinal samples to obtain the first comprehensive coding and miRNA transcriptomic repertoire for postlarval and juvenile vertebral bone, covering different vertebral phenotypes and egg-incubation temperatures related to skeleton health in S. senegalensis. Coding genes, miRNA and pathways regulating bone development and growth were identified. Differential transcriptomic profiles and suggestive mRNA-miRNA interactions were found between postlarvae and juveniles. Bone-related genes and functions were associated with the extracellular matrix, development and regulatory processes, calcium binding, retinol and lipid metabolism or response to stimulus, including those revealed by the miRNA targets related to signaling, cellular and metabolic processes, growth, cell proliferation and biological adhesion. Pathway enrichment associated with fish skeleton were identified when comparing postlarvae and juveniles: growth and bone development functions in postlarvae, while actin cytoskeleton, focal adhesion and proteasome related to bone remodeling in juveniles. The transcriptome data disclosed candidate coding and miRNA gene markers related to bone cell processes, references for functional studies of the anosteocytic bone of S. senegalensis. This study establishes a broad transcriptomic foundation to study healthy and anomalous spines under early thermal conditions across life-stages in S. senegalensis, and for comparative analysis of skeleton homeostasis and pathology in fish and vertebrates.

A common garden experiment supports a genetic component underlying the increased resilience of common cockle (Cerastoderma edule) to the parasite Marteilia cochillia.

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.

Heterogeneous microgeographic genetic structure of the common cockle (Cerastoderma edule) in the Northeast Atlantic Ocean: biogeographic barriers and environmental factors.

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.

Non-synonymous variation and protein structure of candidate genes associated with selection in farm and wild populations of turbot (Scophthalmus maximus).

Non-synonymous variation (NSV) of protein coding genes represents raw material for selection to improve adaptation to the diverse environmental scenarios in wild and livestock populations. Many aquatic species face variations in temperature, salinity and biological factors throughout their distribution range that is reflected by the presence of allelic clines or local adaptation. The turbot (Scophthalmus maximus) is a flatfish of great commercial value with a flourishing aquaculture which has promoted the development of genomic resources. In this study, we developed the first atlas of NSVs in the turbot genome by resequencing 10 individuals from Northeast Atlantic Ocean. More than 50,000 NSVs where detected in the ~ 21,500 coding genes of the turbot genome, and we selected 18 NSVs to be genotyped using a single Mass ARRAY multiplex on 13 wild populations and three turbot farms. We detected signals of divergent selection on several genes related to growth, circadian rhythms, osmoregulation and oxygen binding in the different scenarios evaluated. Furthermore, we explored the impact of NSVs identified on the 3D structure and functional relationship of the correspondent proteins. In summary, our study provides a strategy to identify NSVs in species with consistently annotated and assembled genomes to ascertain their role in adaptation.

A chromosome-level genome assembly enables the identification of the follicule stimulating hormone receptor as the master sex-determining gene in the flatfish Solea senegalensis.

Sex determination (SD) shows huge variation among fish and a high evolutionary rate, as illustrated by the Pleuronectiformes (flatfishes). This order is characterized by its adaptation to demersal life, compact genomes and diversity of SD mechanisms. Here, we assembled the Solea senegalensis genome, a flatfish of great commercial value, into 82 contigs (614 Mb) combining long- and short-read sequencing, which were next scaffolded using a highly dense genetic map (28,838 markers, 21 linkage groups), representing 98.9% of the assembly. Further, we established the correspondence between the assembly and the 21 chromosomes by using BAC-FISH. Whole genome resequencing of six males and six females enabled the identification of 41 single nucleotide polymorphism variants in the follicle stimulating hormone receptor (fshr) consistent with an XX/XY SD system. The observed sex association was validated in a broader independent sample, providing a novel molecular sexing tool. The fshr gene displayed differential expression between male and female gonads from 86 days post-fertilization, when the gonad is still an undifferentiated primordium, concomitant with the activation of amh and cyp19a1a, testis and ovary marker genes, respectively, in males and females. The Y-linked fshr allele, which included 24 nonsynonymous variants and showed a highly divergent 3D protein structure, was overexpressed in males compared to the X-linked allele at all stages of gonadal differentiation. We hypothesize a mechanism hampering the action of the follicle stimulating hormone driving the undifferentiated gonad toward testis.

Sex separation unveils the functional plasticity of the vomeronasal organ in rabbits.

Chemosensory cues are vital for social and sexual behaviours and are primarily detected and processed by the vomeronasal system (VNS), whose plastic capacity has been investigated in mice. However, studying chemosensory plasticity outside of laboratory conditions may give a more realistic picture of how the VNS adapts to a changing environment. Rabbits are a well-described model of chemocommunication since the discovery of the rabbit mammary pheromone and their vomeronasal organ (VNO) transcriptome was recently characterised, a first step to further study plasticity-mediated transcriptional changes. In this study, we assessed the plastic capacity of the rabbit male and female VNO under sex-separation vs. sex-combined scenarios, including adults and juveniles, to determine whether the rabbit VNO is plastic and, if so, whether such plasticity is already established at early stages of life. First, we characterised the number of differentially expressed genes (DEGs) between the VNO of rabbit male and female under sex-separation and compared it to sex-combined individuals, both in adults and juveniles, finding that differences between male and female were larger in a sex-separated scenario. Secondly, we analysed the number of DEGs between sex-separated and sex-combined scenarios, both in males and females. In adults, both sexes showed a high number of DEGs while in juveniles only females showed differences. Additionally, the vomeronasal receptor genes were strikingly downregulated in sex-separated adult females, whereas in juveniles upregulation was shown for the same condition, suggesting a role of VRs in puberty onset. Finally, we described the environment-modulated plastic capacity of genes involved in reproduction, immunity and VNO functional activity, including G-protein coupled receptors. Our results show that sex-separation induces sex- and stage-specific gene expression differences in the VNO of male and female rabbit, both in adults and juveniles. These results bring out for the first time the plastic capacity of the rabbit VNO, supporting its functional adaptation to specifically respond to a continuous changing environment. Finally, species-specific differences and individual variability should always be considered in VNO studies and overall chemocommunication research.

The first high-density genetic map of common cockle (Cerastoderma edule) reveals a major QTL controlling shell color variation.

Shell color shows broad variation within mollusc species and despite information on the genetic pathways involved in shell construction and color has recently increased, more studies are needed to understand its genetic architecture. The common cockle (Cerastoderma edule) is a valuable species from ecological and commercial perspectives which shows important variation in shell color across Northeast Atlantic. In this study, we constructed a high-density genetic map, as a tool for screening common cockle genome, which was applied to ascertain the genetic basis of color variation in the species. The consensus genetic map comprised 19 linkage groups (LGs) in accordance with the cockle karyotype (2n = 38) and spanned 1073 cM, including 730 markers per LG and an inter-marker distance of 0.13 cM. Five full-sib families showing segregation for several color-associated traits were used for a genome-wide association study and a major QTL on chromosome 13 associated to different color-traits was detected. Mining on this genomic region revealed several candidate genes related to shell construction and color. A genomic region previously reported associated with divergent selection in cockle distribution overlapped with this QTL suggesting its putative role on adaptation.

A single genomic region involving a putative chromosome rearrangement in flat oyster (Ostrea edulis) is associated with differential host resilience to the parasite Bonamia ostreae.

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.

Integration of Maps Enables a Cytogenomics Analysis of the Complete Karyotype in Solea senegalensis.

The Pleuronectiformes order, which includes several commercially-important species, has undergone extensive chromosome evolution. One of these species is Solea senegalensis, a flatfish with 2n = 42 chromosomes. In this study, a cytogenomics approach and integration with previous maps was applied to characterize the karyotype of the species. Synteny analysis of S. senegalensis was carried out using two flatfish as a reference: Cynoglossus semilaevis and Scophthalmus maximus. Most S. senegalensis chromosomes (or chromosome arms for metacentrics and submetacentrics) showed a one-to-one macrosyntenic pattern with the other two species. In addition, we studied how repetitive sequences could have played a role in the evolution of S. senegalensis bi-armed (3, and 5-9) and acrocentric (11, 12 and 16) chromosomes, which showed the highest rearrangements compared with the reference species. A higher abundance of TEs (Transposable Elements) and other repeated elements was observed adjacent to telomeric regions on chromosomes 3, 7, 9 and 16. However, on chromosome 11, a greater abundance of DNA transposons was detected in interstitial BACs. This chromosome is syntenic with several chromosomes of the other two flatfish species, suggesting rearrangements during its evolution. A similar situation was also found on chromosome 16 (for microsatellites and low complexity sequences), but not for TEs (retroelements and DNA transposons). These differences in the distribution and abundance of repetitive elements in chromosomes that have undergone remodeling processes during the course of evolution also suggest a possible role for simple repeat sequences in rearranged regions.

Full-Length Transcriptome Sequences Provide Insight Into Hermaphroditism of Freshwater Pearl Mussel Hyriopsis schlegelii.

The freshwater mussel Hyriopsis schlegelii is a cultured bivalve in China, and the quality of the pearls produced is affected by the type of gonads. However, because of the lack of a published genome and the complexity of sex determination, research on sex reversal and development of this species is limited. In this study, Illumina RNA-seq and PacBio Isoform Sequencing (Iso-Seq) were combined to analyze the gonads of H. schlegelii. A total of 201,481 high-quality transcripts were generated. The study identified 7,922 differentially expressed genes in three comparison group (females versus males, hermaphrodites versus females, and hermaphrodites versus males). Twenty-four genes were identified as potential sex-related genes, including sox9 and wnt4 involved in sex determination, and vtg, cyp17a1 and 17ß-hsd2 involved in gonadal development. We also speculated a possible pathways for the formation of hermaphroditism in H. schlegelii. The data provide a clear view of the transcriptome for H. schlegelii gonads and will be valuable in elucidating the mechanisms of gonad development.

Genomic survey of edible cockle (Cerastoderma edule) in the Northeast Atlantic: A baseline for sustainable management of its wild resources.

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.

Assessment of Biostimulation Methods Based on Chemical Communication in Female Doe Reproduction.

Biostimulation is an animal management practice that helps improve reproductive parameters by modulating animal sensory systems. Chemical signals, mostly known as pheromones, have a great potential in this regard. This study was conducted to determine the influence of short-term female rabbit exposure to different conditions, mainly pheromone-mediated, on reproductive parameters of inseminated does. Groups of 60 females/each were exposed to (1) female urine, (2) male urine, (3) seminal plasma and (4) female-female (F-F) separated, just before artificial insemination, and compared to a 'golden method' female-female interaction. The following reproductive parameters were analyzed for each group: receptivity (vulvar color), fertility (kindling rate), prolificacy and number of born alive and dead kits/litter. Our results showed that the biostimulation methods employed in this experiment did not significantly improve any of the analyzed parameters. However, female doe exposure to urine, especially to male urine, showed no significant higher fertility values (95.4%) when compared to the rest of the experimental conditions (on average 92.4%). Female-female interaction before artificial insemination, which is a common practice in rabbit farms, showed similar results as not establishing social interaction (F-F separated), which suggests that F-F interaction could be replaced by F-F separated, therefore avoiding unnecessary animal management and time cost. On the other hand, fertility ranges were lower for animals with a pale vulvar color whereas no differences were noticed among the other three colors which measure receptivity (pink, red, purple), thus suggesting that these three colors could be grouped together. Future studies should aim at determining potential chemical cues/pheromones released through bodily secretions that influence reproduction in rabbits, therefore contributing to animal welfare and to a natural image of animal production.

Genomic Hatchery Introgression in Brown Trout (Salmo trutta L.): Development of a Diagnostic SNP Panel for Monitoring the Impacted Mediterranean Rivers.

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.

Caracterización de la lesión renal aguda en un hospital de la República Dominicana

Resumen: Objetivo: Caracterizar la lesión renal aguda (LRA) para describir el patrón etiológico, manifestaciones clínicas y de laboratorio. Material y métodos: Diseño observacional, retrospectivo tipo serie de casos. Se analizaron los registros de 88 pacientes, atendidos en la Unidad de Nefrología de un Hospital en la República Dominicana, sin importar el estadio de la enfermedad. Resultados: Encontramos que la edad media es 64 años con una TFG 17 mL/min/1.73 m2. El principal motivo de ingreso como de interconsulta es la "elevación de azoados" con 7.51% (n= 19 / 253) y 29.90% (n=58 / 194), respectivamente. El diagnóstico final más común fue la anemia con un 17.09% (n=50 / 396), la hipertensión arterial con un 34.16% (n=55 / 151) fue la comorbilidad predominante. La mayoría de pacientes no tuvo ningún antecedente medicamentoso con un 31.17% (n=48 / 154), la intervención más común fue el uso de fluidos en un 40.83% (n=49 / 120). El hallazgo clínico más descrito fueron los campos pulmonares hipoventilados 42,72% (n=44 / 103). Conclusiones: Recomendamos poner en funcionamiento acciones que ayuden con un oportuno diagnóstico y tratamiento temprano de la LRA en el hospital.

Abstract: Objective: To characterize acute kidney injury (AKI) and describe its etiological pattern, clinical and laboratory manifestations. Material and methods: Retrospective, observational design, case series type. We analyzed the records of 88 patients treated in the nephrology unit of a hospital in the Dominican Republic, regardless of the stage of the disease. Results: We found that the average age is 64 years with a GFR 17mL/min/1.73 m2. The main reason for admission and consultation is the "elevation of BUN and creatinine levels" with 7.51% (n= 19 / 253) and 29.90% (n=58 / 194), respectively. The most common final diagnosis was anemia 17.09% (n=50 / 396), hypertension with 34.16% (n=55 / 151) was the predominant comorbidity. The majority of patients had no medication history with 31.17% (n=48 / 154), the most common intervention performed was fluid use with 40.83% (n=49 / 120). The most frequently described clinical finding was hypoventilated lungs with 42, 72% (n=44 / 103). Conclusions: We recommend the implementation of actions to improve the timely diagnosis and management of AKI in the hospital.

Resumo: Objetivo: Caracterizar a lesão renal aguda (LRA) prá descrever o padrão etiológico, manifestações clinicas e do laboratório. Métodos: Modelo observacional, retrospectivo tipo série de casos. Foram analisados os registros de 88 pacientes, antendidos na Unidade de Nefrología dum Hospital na República Dominicana, sem importar o estágio da doença. Resultados: Encontramos que a edade média é 64 anos com uma TFG de 17 mL/min/1.73 m2. A razão principal da admissão e interconsulta é a "elevação de azoados" com 7.51% (n=19/253) e 29.90% (n=58/194), respectivamente. O diagnóstico final mais comum foi a anemia com um 17.09% (n=50/396), a hipertensão arterial com um 34.16% (n=55/151) foi a comorbidade predominante. A maioria dos pacientes não tinha histórico medicinal com um 31.17% (n=48/154), a intervenção mais comum foi o uso de fluidos num 40.83% (n=49/120). O achado clínico mais comum foi campos pulmonares hipoventilados 47.72% (n=44/103). Conclusões: Recomendamos colocar em operação ações que ajudam com um diagnóstico oportuno e um tratamento precoce da LRA no hospital.

Estimates of recent and historical effective population size in turbot, seabream, seabass and carp selective breeding programmes.

BACKGROUND: The high fecundity of fish species allows intense selection to be practised and therefore leads to fast genetic gains. Based on this, numerous selective breeding programmes have been started in Europe in the last decades, but in general, little is known about how the base populations of breeders have been built. Such knowledge is important because base populations can be created from very few individuals, which can lead to small effective population sizes and associated reductions in genetic variability. In this study, we used genomic information that was recently made available for turbot (Scophthalmus maximus), gilthead seabream (Sparus aurata), European seabass (Dicentrarchus labrax) and common carp (Cyprinus carpio) to obtain accurate estimates of the effective size for commercial populations. METHODS: Restriction-site associated DNA sequencing data were used to estimate current and historical effective population sizes. We used a novel method that considers the linkage disequilibrium spectrum for the whole range of genetic distances between all pairs of single nucleotide polymorphisms (SNPs), and thus accounts for potential fluctuations in population size over time. RESULTS: Our results show that the current effective population size for these populations is small (equal to or less than 50 fish), potentially putting the sustainability of the breeding programmes at risk. We have also detected important drops in effective population size about five to nine generations ago, most likely as a result of domestication and the start of selective breeding programmes for these species in Europe. CONCLUSIONS: Our findings highlight the need to broaden the genetic composition of the base populations from which selection programmes start, and suggest that measures designed to increase effective population size within all farmed populations analysed here should be implemented in order to manage genetic variability and ensure the sustainability of the breeding programmes.

Detection of Genomic Regions with Pleiotropic Effects for Growth and Carcass Quality Traits in the Rubia Gallega Cattle Breed.

The breeding scheme in the Rubia Gallega cattle population is based upon traits measured in farms and slaughterhouses. In recent years, genomic evaluation has been implemented by using a ssGBLUP (single-step Genomic Best Linear Unbiased Prediction). This procedure can reparameterized to perform ssGWAS (single-step Genome Wide Association Studies) by backsolving the SNP (single nucleotide polymorphisms) effects. Therefore, the objective of this study was to identify genomic regions associated with the genetic variability in growth and carcass quality traits. We implemented a ssGBLUP by using a database that included records for Birth Weight (BW-327,350 records-), Weaning Weight (WW-83,818-), Cold Carcass Weight (CCW-91,621-), Fatness (FAT-91,475-) and Conformation (CON-91,609-). The pedigree included 464,373 individuals, 2449 of which were genotyped. After a process of filtering, we ended up using 43,211 SNP markers. We used the GBLUP and SNPBLUP model equivalences to obtain the effects of the SNPs and then calculated the percentage of variance explained by the regions of the genome between 1 Mb. We identified 7 regions of the genome for CCW; 8 regions for BW, WW, FAT and 9 regions for CON, which explained the percentage of variance above 0.5%. Furthermore, a number of the genome regions had pleiotropic effects, located at: BTA1 (131-132 Mb), BTA2 (1-11 Mb), BTA3 (32-33 Mb), BTA6 (36-38 Mb), BTA16 (24-26 Mb), and BTA 21 (56-57 Mb). These regions contain, amongst others, the following candidate genes: NCK1, MSTN, KCNA3, LCORL, NCAPG, and RIN3.

Blood Transcriptomics of Turbot Scophthalmus maximus: A Tool for Health Monitoring and Disease Studies.

Blood transcriptomics is emerging as a relevant tool to monitor the status of the immune system and assist in diagnosis, prognosis, treatment and pathogenesis studies of diseases. In fish pathology, the potential of transcriptome profiling of blood is still poorly explored. Here, RNA sequencing was applied to analyze the blood transcriptional profile of turbot (Scophthalmus maximus), the most important farmed flatfish. The study was conducted in healthy specimens and specimens parasitized by the myxozoan Enteromyxum scophthalmi, which causes one of the most devastating diseases in turbot aquaculture. The blood of healthy turbot showed a transcriptomic profile mainly related to erythrocyte gas transportation function, but also to antigen processing and presentation. In moderately infected turbot, the blood reflected a broad inhibition of the immune response. Particularly, down-regulation of the B cell receptor signaling pathway was shared with heavily parasitized fish, which showed larger transcriptomic changes, including the activation of the inflammatory response. Turbot response to enteromyxosis proved to be delayed, dysregulated and ineffective in stopping the infection. The study evinces that blood transcriptomics can contribute to a better understanding of the teleost immune system and serve as a reliable tool to investigate the physiopathological status of fish.

A genome-wide association study, supported by a new chromosome-level genome assembly, suggests sox2 as a main driver of the undifferentiatiated ZZ/ZW sex determination of turbot (Scophthalmus maximus).

BACKGROUND: Understanding sex determination (SD) across taxa is a major challenge for evolutionary biology. The new genomic tools are paving the way to identify genomic features underlying SD in fish, a group frequently showing limited sex chromosome differentiation and high SD evolutionary turnover. Turbot (Scophthalmus maximus) is a commercially important flatfish with an undifferentiated ZW/ZZ SD system and remarkable sexual dimorphism. Here we describe a new long-read turbot genome assembly used to disentangle the genetic architecture of turbot SD by combining genomics and classical genetics approaches. RESULTS: The new turbot genome assembly consists of 145 contigs (N50 = 22.9 Mb), 27 of them representing >95% of its estimated genome size. A genome wide association study (GWAS) identified a ~ 6.8 Mb region on chromosome 12 associated with sex in 69.4% of the 36 families analyzed. The highest associated markers flanked sox2, the only gene in the region showing differential expression between sexes before gonad differentiation. A single SNP showed consistent differences between Z and W chromosomes. The analysis of a broad sample of families suggested the presence of additional genetic and/or environmental factors on turbot SD. CONCLUSIONS: The new chromosome-level turbot genome assembly, one of the most contiguous fish assemblies to date, facilitated the identification of sox2 as a consistent candidate gene putatively driving SD in this species. This chromosome SD system barely showed any signs of differentiation, and other factors beyond the main QTL seem to control SD in a certain proportion of families.

Low impact of different SNP panels from two building-loci pipelines on RAD-Seq population genomic metrics: case study on five diverse aquatic species.

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.

First description outside Europe of the emergent pathogen Vibrio europaeus in shellfish aquaculture.

Vibrio europaeus is an emergent pathogen affecting the most important bivalve species reared in Spanish and French hatcheries. Using a genomic approach, we identified V. europaeus outside Europe for the first time from massive larval mortalities of scallop (Argopecten purpuratus) in Chile and from seawater near a shellfish hatchery in the US West Coast. Results show the worldwide spreading and potential impact of V. europaeus for aquaculture; these four countries are among the 10 major producers of mollusks. Pathogenicity of V. europaeus was demonstrated for the first time towards scallop, the second most important species for Chilean mariculture.