Investigating the role of genetic variation in vgll3 and six6 in the domestication of gilthead seabream (Sparus aurata Linnaeus) and European seabass (Dicentrarchus labrax Linnaeus)

Gene function conservation is crucial in molecular ecology, especially for key traits like growth and maturation in teleost fish. The vgll3 and six6 genes are known to influence age-at-maturity in Atlantic salmon, but their impact on other fish species is poorly understood. Here, we investigated the association of vgll3 and six6 in the domestication of gilthead seabream and European seabass, both undergoing selective breeding for growth-related traits in the Mediterranean. We analysed two different sets of samples using two different genotyping approaches. The first dataset comprised farmed and wild populations from Greece, genotyped for SNPs within the two genes ('gene-level genotyping'). The second dataset examined 300-600 k SNPs located in the chromosomes of the two genes, derived from a meta-analysis of a Pool-Seq experiment involving farmed and wild populations distributed widely across the Mediterranean ('chromosome-level genotyping'). The gene-level analysis revealed a statistically significant allele frequency differences between farmed and wild populations on both genes in each species. This finding was partially supported by the chromosome-level analysis, identifying highly differentiated regions may be involved in the domestication process at varying distances from the candidate genes. Noteworthy genomic features were found, such as a CpG island in gilthead seabream and novel candidate genes in European seabass, warranting further investigation. These findings support a putative role of vgll3 and six6 in the maturation and growth of gilthead seabream and European seabass, emphasizing the need for further research on their conserved function.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture

BACKGROUND

The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.

RESULTS

We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.

CONCLUSIONS

We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Genetic Identification and Traceability of Insect Meals

Insects have been proposed as a rich alternative source of protein for the partial or total replacement of fishmeal in aquaculture. For maximum safety and effectiveness of insect meals, control of the quality composition of these products is considered mandatory. The aim of this study was the genetic analysis of the composition of commercially available insect meals at the species level. Commercially available Hermetia illucens, Tenebrio molitor and Musca domestica individuals, as well as nine insect meals produced from these species, were analyzed. The genetic identification of insects at the species level was based on a COI fragment, and analysis of the insect meals' composition was performed with the processes of cloning and colony PCR. Genetic analysis indicated that the commercially available larvae morphologically identified as Musca domestica belonged to the species Muscina stabulans. In the commercially available insect meals, no other animal species was identified beyond the expected one. However, in the insect meal produced for research purposes, fungal growth was detected. The used methodology, herein, allows for the qualitative genetic identification of insect meals and could be included in the methods of traceability of products containing insects and other animal species.

Advances in understanding the mitogenic, metabolic, and cell death signaling in teleost development: the case of greater amberjack (Seriola dumerili, Risso 1810)

Cell growth and differentiation signals of insulin-like growth factor-1 (IGF-1), a key regulator in embryonic and postnatal development, are mediated through the IGF-1 receptor (IGF-1R), which activates several downstream pathways. The present study aims to address crucial organogenesis and development pathways including Akt, MAPKs, heat shock response, apoptotic and autophagic machinery, and energy metabolism in relation to IGF-1R activation during five developmental stages of reared Seriola dumerili: 1 day prior to hatching fertilized eggs (D-1), hatching day (D0), 3 days post-hatching larvae (D3), 33 (D33) and 46 (D46) days post-hatching juveniles. During both the fertilized eggs stage and larval-to-juvenile transition, IGF-1R/Akt pathway activation may mediate the hypertrophic signaling, while p44/42 MAPK phosphorylation was apparent at S. dumerili post-hatching processes and juvenile organs completion. On the contrary, apoptosis was induced during embryogenesis and autophagy at hatching day indicating a potential involvement in morphogenetic rearrangements and yolk-sac reserves depletion. Larvae morphogenesis was accompanied by a metabolic turnover with increased substantial energy consumption. The findings of the present study demonstrate the developmental stages-specific shift in critical signaling pathways during the ontogeny of reared S. dumerili.

DNA barcoding identification of Greek freshwater fishes

Biodiversity is a key factor for the functioning and efficiency of an ecosystem. Greece, though covering a relatively small surface area, hosts a great deal of species diversity. This is especially true for freshwater fishes. In recent years, the traditional methods of species identification have been supplemented by the use of molecular markers. The present study therefore aims to extensively produce DNA barcodes for Greek freshwater fish species and investigate thoroughly if the presently accepted species classification is in agreement with molecular data. A 624-bases long fragment of the COI gene was sequenced, from 406 freshwater fish specimens belonging to 24 genera and originating from 18 lake and river sites. These sequences were used along with 596 sequences from the same genera, recovered from BOLD, for the construction of phylogenetic trees and the estimation of genetic distances between individuals. In total, 1002 sequences belonging to 72 species were analyzed. The method was found to be effective for 55 of 72 studied species. 17 closely related species with low interspecific genetic distances were observed, for which further study is proposed. It should also be noted that, in four cases, cryptic diversity was observed, where groups originally identified as one species exhibited genetic distance great enough to be separated into discrete species. Region specific haplotypes were also detected within populations of 14 species, giving the possibility to identify even the geographic origin of a species. Our findings are discussed in the light of the rich history of the Balkan peninsula and provide a significant steppingstone for the further study of Greek and European freshwater fish biodiversity.

FoodOmicsGR_RI. A Consortium for Comprehensive Molecular Characterisation of Food Products

The national infrastructure FoodOmicsGR_RI coordinates research efforts from eight Greek Universities and Research Centers in a network aiming to support research and development (R&D) in the agri-food sector. The goals of FoodOmicsGR_RI are the comprehensive in-depth characterization of foods using cutting-edge omics technologies and the support of dietary/nutrition studies. The network combines strong omics expertise with expert field/application scientists (food/nutrition sciences, plant protection/plant growth, animal husbandry, apiculture and 10 other fields). Human resources involve more than 60 staff scientists and more than 30 recruits. State-of-the-art technologies and instrumentation is available for the comprehensive mapping of the food composition and available genetic resources, the assessment of the distinct value of foods, and the effect of nutritional intervention on the metabolic profile of biological samples of consumers and animal models. The consortium has the know-how and expertise that covers the breadth of the Greek agri-food sector. Metabolomics teams have developed and implemented a variety of methods for profiling and quantitative analysis. The implementation plan includes the following research axes: development of a detailed database of Greek food constituents; exploitation of "omics" technologies to assess domestic agricultural biodiversity aiding authenticity-traceability control/certification of geographical/genetic origin; highlighting unique characteristics of Greek products with an emphasis on quality, sustainability and food safety; assessment of diet's effect on health and well-being; creating added value from agri-food waste. FoodOmicsGR_RI develops new tools to evaluate the nutritional value of Greek foods, study the role of traditional foods and Greek functional foods in the prevention of chronic diseases and support health claims of Greek traditional products. FoodOmicsGR_RI provides access to state-of-the-art facilities, unique, well-characterised sample sets, obtained from precision/experimental farming/breeding (milk, honey, meat, olive oil and so forth) along with more than 20 complementary scientific disciplines. FoodOmicsGR_RI is open for collaboration with national and international stakeholders.

Genome-wide analysis clarifies the population genetic structure of wild gilthead sea bream (Sparus aurata)

Gilthead sea bream is an important target for both recreational and commercial fishing in Europe, where it is also one of the most important cultured fish. Its distribution ranges from the Mediterranean to the African and European coasts of the North-East Atlantic. Until now, the population genetic structure of this species in the wild has largely been studied using microsatellite DNA markers, with minimal genetic differentiation being detected. In this geographically widespread study, 958 wild gilthead sea bream from 23 locations within the Mediterranean Sea and Atlantic Ocean were genotyped at 1159 genome-wide SNP markers by RAD sequencing. Outlier analyses identified 18 loci potentially under selection. Neutral marker analyses identified weak subdivision into three genetic clusters: Atlantic, West, and East Mediterranean. The latter group could be further subdivided into an Ionian/Adriatic and an Aegean group using the outlier markers alone. Seascape analysis suggested that this differentiation was mainly due to difference in salinity, this being also supported by preliminary genomic functional analysis. These results are of fundamental importance for the development of proper management of this species in the wild and are a first step toward the study of the potential genetic impact of the sea bream aquaculture industry.

Genetic structure and divergence of tench Tinca tinca European populations

The tench Tinca tinca is a freshwater species with human-mediated translocations, aquaculture interest and limited information on its genetic structure. mtDNA sequencing analysis of control region and two genes in 50 individuals from five European populations identified two phylogroups, with greater variability than that reported until now, and a hybridization zone in the Danube River region. Restriction analyses of additional samples reveal the complicated genetic structure characteristics of tench's wild and translocated populations, supporting future breeding practices.

The genetic population structure and temporal genetic stability of gilthead sea bream Sparus aurata populations in the Aegean and Ionian Seas, using microsatellite DNA markers

We examined 662 gilthead sea bream Sparus aurata from wild samples of the species in the Aegean and Ionian Seas, using 20 EST-linked microsatellite markers, in three multiplex panels, as well as seven anonymous loci. Most of the markers were revealed to be highly polymorphic. We found low genetic differentiation between the sampling stations/areas with total F_ST 0.002 (P < 0.05). Based on comparison of five temporal samples, our results indicate genetic data consistency over time for all tested samples, pointing to stable populations, despite reported repeated escape events. Our results confirm the genetic population structure previously observed in these specific areas, using by far more markers than in previous studies in both coding and non-coding DNA loci. The limited genetic structure and the temporal genetic stability indicate neither major genetic differentiation of local populations by geographic isolation nor influence from anthropogenic factors. These results provide a baseline for future reference in any management programme of both wild and farmed population of S. aurata as well as of other aquaculture species with a potential introgression among farmed and wild populations.

Performance and precision of double digestion RAD (ddRAD) genotyping in large multiplexed datasets of marine fish species

The development of Genotyping-By-Sequencing (GBS) technologies enables cost-effective analysis of large numbers of Single Nucleotide Polymorphisms (SNPs), especially in "non-model" species. Nevertheless, as such technologies enter a mature phase, biases and errors inherent to GBS are becoming evident. Here, we evaluated the performance of double digest Restriction enzyme Associated DNA (ddRAD) sequencing in SNP genotyping studies including high number of samples. Datasets of sequence data were generated from three marine teleost species (>5500 samples, >2.5 × 10¹² bases in total), using a standardized protocol. A common bioinformatics pipeline based on STACKS was established, with and without the use of a reference genome. We performed analyses throughout the production and analysis of ddRAD data in order to explore (i) the loss of information due to heterogeneous raw read number across samples; (ii) the discrepancy between expected and observed tag length and coverage; (iii) the performances of reference based vs. de novo approaches; (iv) the sources of potential genotyping errors of the library preparation/bioinformatics protocol, by comparing technical replicates. Our results showed use of a reference genome and a posteriori genotype correction improved genotyping precision. Individual read coverage was a key variable for reproducibility; variance in sequencing depth between loci in the same individual was also identified as an important factor and found to correlate to tag length. A comparison of downstream analysis carried out with ddRAD vs single SNP allele specific assay genotypes provided information about the levels of genotyping imprecision that can have a significant impact on allele frequency estimations and population assignment. The results and insights presented here will help to select and improve approaches to the analysis of large datasets based on RAD-like methodologies.

Greece: a Balkan subrefuge for a remnant red deer (cervus elaphus) population

A number of phylogeographic studies have revealed the existence of multiple ice age refugia within the Balkan Peninsula, marking it as a biodiversity hotspot. Greece has been reported to harbor genetically differentiated lineages from the rest of Balkans for a number of mammal species. We therefore searched for distinct red deer lineages in Greece, by analyzing 78 samples originating from its last population in Parnitha Mountain (Central Greece). Additionally, we tested the impact of human-induced practices on this population. The presence of 2 discrete mtDNA lineages was inferred: 1) an abundant one not previously sampled in the Balkans and 2) a more restricted one shared with other Balkan populations, possibly the result of successful translocations of Eastern European individuals. Microsatellite-based analyses of 14 loci strongly support the existence of 2 subpopulations with relative frequencies similar to mitochondrial analyses. This study stresses the biogeographic importance of Central Greece as a separate Last Glacial Maximum period refugium within the Balkans. It also delineates the possible effects that recent translocations of red deer populations had on the genetic structuring within Parnitha. We suggest that the Greek red deer population of Parnitha is genetically distinct, and restocking programs should take this genetic evidence into consideration.