Unraveling the genotype by environment interaction in a thermosensitive fish with a polygenic sex determination system.

In most animals, sex determination occurs at conception, when sex chromosomes are segregated following Mendelian laws. However, in multiple reptiles and fishes, this genetic sex can be overridden by external factors after fertilization or birth. In some species, the genetic sex may also be governed by multiple genes, further limiting our understanding of sex determination in such species. We used the European sea bass (Dicentrarchus labrax) as a model and combined genomic (using a single nucleotide polymorphism chip) and transcriptomic (RNA-Sequencing) approaches to thoroughly depict this polygenic sex determination system and its interaction with temperature. We estimated genetic sex tendency (eGST), defined as the estimated genetic liability to become a given sex under a liability threshold model for sex determination, which accurately predicts the future phenotypic sex. We found evidence that energetic pathways, concerning the regulation of lipids and glucose, are involved in sex determination and could explain why females tend to exhibit higher energy levels and improved growth compared to males. Besides, early exposure to high-temperature up-regulated sox3, followed by sox9a in individuals with intermediate eGST, but not in individuals showing highly female-biased eGST, providing the most parsimonious explanation for temperature-induced masculinization. This gonadal state was maintained likely by DNA methylation and the up-regulation of several genes involved in histone modifications, including jmjd1c Overall, we describe a sex determination system resulting from continuous genetic and environmental influences in an animal. Our results provide significant progress in our understanding of the mechanisms underlying temperature-induced masculinization in fish.

In a marine teleost, the significance of oxygen supply for acute thermal tolerance depends upon the context and the endpoint used.

Eight juvenile European seabass were exposed to two thermal ramping protocols with different levels of aerobic activity and tolerance endpoint: the critical thermal maximum for swimming (CTSmax) while exercising aerobically until fatigue and the critical thermal maximum (CTmax) under static conditions until loss of equilibrium (LOE). In the CTSmax protocol, warming caused a profound increase in the rate of oxygen uptake (MO2), culminating in a gait transition from steady aerobic towards unsteady anaerobic swimming, then fatigue at 30.3±0.4°C (mean±s.e.m.). Gait transition and fatigue presumably indicate an oxygen limitation, which reflects the inability to meet the combined demands of swimming plus warming. The CTmax protocol also elicited an increase in MO2, culminating in LOE at 34.0±0.4°C, which is significantly warmer than fatigue at CTSmax. The maximum MO2 achieved in the CTmax protocol was, however, less than 30% of that achieved in the CTSmax protocol. Therefore, the static CTmax did not exploit full cardiorespiratory capacity for oxygen supply, indicating that LOE was not caused by systemic oxygen limitation. Consequently, systemic oxygen supply can be significant for tolerance of acute warming in seabass but this depends upon the physiological context and the endpoint used.

Whole-genome sequencing identifies interferon-induced protein IFI6/IFI27-like as a strong candidate gene for VNN resistance in European sea bass.

BACKGROUND: Viral nervous necrosis (VNN) is a major disease that affects European sea bass, and understanding the biological mechanisms that underlie VNN resistance is important for the welfare of farmed fish and sustainability of production systems. The aim of this study was to identify genomic regions and genes that are associated with VNN resistance in sea bass. RESULTS: We generated a dataset of 838,451 single nucleotide polymorphisms (SNPs) identified from whole-genome sequencing (WGS) in the parental generation of two commercial populations (A: 2371 individuals and B: 3428 individuals) of European sea bass with phenotypic records for binary survival in a VNN challenge. For each population, three cohorts were submitted to a red-spotted grouper nervous necrosis virus (RGNNV) challenge by immersion and genotyped on a 57K SNP chip. After imputation of WGS SNPs from their parents, quantitative trait loci (QTL) were mapped using a Bayesian sparse linear mixed model (BSLMM). We found several QTL regions that were specific to one of the populations on different linkage groups (LG), and one 127-kb QTL region on LG12 that was shared by both populations and included the genes ZDHHC14, which encodes a palmitoyltransferase, and IFI6/IFI27-like, which encodes an interferon-alpha induced protein. The most significant SNP in this QTL region was only 1.9 kb downstream of the coding sequence of the IFI6/IFI27-like gene. An unrelated population of four large families was used to validate the effect of the QTL. Survival rates of susceptible genotypes were 40.6% and 45.4% in populations A and B, respectively, while that of the resistant genotype was 66.2% in population B and 78% in population A. CONCLUSIONS: We have identified a genomic region that carries a major QTL for resistance to VNN and includes the ZDHHC14 and IFI6/IFI27-like genes. The potential involvement of the interferon pathway, a well-known anti-viral defense mechanism in several organisms (chicken, human, or fish), in survival to VNN infection is of particular interest. Our results can lead to major improvements for sea bass breeding programs through marker-assisted genomic selection to obtain more resistant fish.

Natural cortisol production is not linked to the sexual fate of European sea bass.

In this study, we aimed to investigate the relationship between cortisol and the determination of sexual fate in the commercially important European sea bass (Dicentrarchus labrax). To test our hypothesis, we designed two temperature-based experiments (19 ℃, 21 ℃ and 23 ℃, experiment 1; 16 ℃ and 21 ℃, experiment 2) to assess the effects of these thermal treatments on European sea bass sex determination and differentiation. In the fish from the first experiment, we evaluated whether blood cortisol levels and expression of stress key regulatory genes were different between differentiating (149 to 183 dph) males and females. In the second experiment, we assessed whether cortisol accumulated in scales over time during the labile period for sex determination as well as the neuroanatomical localisation of brain cells expressing brain aromatase (cyp19a1b) and corticotropin-releasing factor (crf) differed between males and females undergoing molecular sex differentiation (117 to 124 dph). None of the gathered results allowed to detect differences between males and females regarding cortisol production and regulatory mechanisms. Altogether, our data provide strong physiological, molecular and histochemical evidence, indicating that in vivo cortisol regulation has no major effects on the sex of European sea bass.

Population-specific variations of the genetic architecture of sex determination in wild European sea bass Dicentrarchus labrax L.

Polygenic sex determination (PSD) may show variations in terms of genetic and environmental components between populations of fish species exposed/adapted to different environments. The European sea bass (Dicentrarchus labrax) is an interesting model, combining both a PSD system and a genetic subdivision into an Atlantic and a Mediterranean lineage, with genetic substructures within the Mediterranean Sea. Here, we produced experimental progeny crosses (N = 927) from broodstock sampled in four wild populations (North Atlantic, NAT; Western Mediterranean, WEM; North-Eastern Mediterranean, NEM; South-Eastern Mediterranean, SEM). We found less females than males in the progeny, both in the global dataset (32.5%) and within each paternal group (from 25.1% for NEM to 39.0% for WEM), with significant variation among populations, dams, and sires. Sex, body weight (BW), and body length (BL) showed moderate heritability (0.52 ± 0.17, 0.46 ± 0.17, 0.34 ± 0.15, respectively) and sex was genetically correlated with BW and BL, with rAsex/BW = 0.69 ± 0.12 and rA sex/BL = 0.66 ± 0.13. A weighted GWAS performed both on the global dataset and within each paternal group revealed a different genetic architecture of sex determination between Atlantic and Mediterranean populations (9 QTLs found in NAT, 7 in WEM, 5 in NEM, and 4 in SEM, with a cumulated variance explained of 27.04%, 21.87%, 15.89%, and 12.10%, respectively) and a more similar genetic architecture among geographically close populations compared to geographically distant populations, consistent with the hypothesis of a population-specific evolution of polygenic sex determination systems in different environments.

Genome-wide association and genomic prediction of resistance to viral nervous necrosis in European sea bass (Dicentrarchus labrax) using RAD sequencing.

BACKGROUND: European sea bass (Dicentrarchus labrax) is one of the most important species for European aquaculture. Viral nervous necrosis (VNN), commonly caused by the redspotted grouper nervous necrosis virus (RGNNV), can result in high levels of morbidity and mortality, mainly during the larval and juvenile stages of cultured sea bass. In the absence of efficient therapeutic treatments, selective breeding for host resistance offers a promising strategy to control this disease. Our study aimed at investigating genetic resistance to VNN and genomic-based approaches to improve disease resistance by selective breeding. A population of 1538 sea bass juveniles from a factorial cross between 48 sires and 17 dams was challenged with RGNNV with mortalities and survivors being recorded and sampled for genotyping by the RAD sequencing approach. RESULTS: We used genome-wide genotype data from 9195 single nucleotide polymorphisms (SNPs) for downstream analysis. Estimates of heritability of survival on the underlying scale for the pedigree and genomic relationship matrices were 0.27 (HPD interval 95%: 0.14-0.40) and 0.43 (0.29-0.57), respectively. Classical genome-wide association analysis detected genome-wide significant quantitative trait loci (QTL) for resistance to VNN on chromosomes (unassigned scaffolds in the case of 'chromosome' 25) 3, 20 and 25 (P < 1e06). Weighted genomic best linear unbiased predictor provided additional support for the QTL on chromosome 3 and suggested that it explained 4% of the additive genetic variation. Genomic prediction approaches were tested to investigate the potential of using genome-wide SNP data to estimate breeding values for resistance to VNN and showed that genomic prediction resulted in a 13% increase in successful classification of resistant and susceptible animals compared to pedigree-based methods, with Bayes A and Bayes B giving the highest predictive ability. CONCLUSIONS: Genome-wide significant QTL were identified but each with relatively small effects on the trait. Tests of genomic prediction suggested that incorporating genome-wide SNP data is likely to result in higher accuracy of estimated breeding values for resistance to VNN. RAD sequencing is an effective method for generating such genome-wide SNPs, and our findings highlight the potential of genomic selection to breed farmed European sea bass with improved resistance to VNN.

How could fully scaled carps appear in natural waters in Madagascar?

The capacity of organisms to rapidly evolve in response to environmental changes is a key feature of evolution, and studying mutation compensation is a way to evaluate whether alternative routes of evolution are possible or not. Common carps (Cyprinus carpio) carrying a homozygous loss-of-function mutation for the scale cover gene fgfr1a1, causing the 'mirror' reduced scale cover, were introduced in Madagascar a century ago. Here we show that carps in Malagasy natural waters are now predominantly covered with scales, though they still all carry the homozygous mutation. We also reveal that the number of scales in mutated carps is under strong polygenic genetic control, with a heritability of 0.49. As a whole, our results suggest that carps submitted to natural selection could evolve a wild-type-like scale cover in less than 40 generations from standing polygenic genetic variation, confirming similar findings mainly retrieved from model organisms.

APIS: an updated parentage assignment software managing triploids induced from diploid parents.

In aquaculture, sterile triploids are commonly used for production as sterility gives them potential gains in growth, yields and quality. However, they cannot be reproduced, and DNA parentage assignment to their diploid or tetraploid parents is required to estimate breeding values for triploid phenotypes. No publicly available software has the ability to assign triploids to their parents. Here, we updated the R package APIS to support triploids induced from diploid parents. First, we created new exclusion and likelihood tables that account for the double allelic contribution of the dam and the recombination that can occur during female meiosis. As the effective recombination rate of each marker with the centromere is usually unknown, we set it at 0.5 and found that this value maximises the assignment rate even for markers with high or low recombination rates. The number of markers needed for a high true assignment rate did not strongly depend on the proportion of missing parental genotypes. The assignment power was however affected by the quality of the markers (minor allele frequency, call rate). Altogether, 96 to 192 SNPs were required to have a high parentage assignment rate in a real rainbow trout dataset of 1232 triploid progenies from 288 parents. The likelihood approach was more efficient than exclusion when the power of the marker set was limiting. When more markers were used, exclusion was more advantageous, with sensitivity reaching unity, very low False Discovery Rate (<0.01) and excellent specificity (0.96-0.99). Thus, APIS provides an efficient solution to assign triploids to their diploid parents.

Functional Additives in a Selected European Sea Bass (Dicentrarchus labrax) Genotype: Effects on the Stress Response and Gill Antioxidant Response to Hydrogen Peroxide (H2O2) Treatment.

Functional ingredients have profiled as suitable candidates for reinforcing the fish antioxidant response and stress tolerance. In addition, selective breeding strategies have also demonstrated a correlation between fish growth performance and susceptibility to stressful culture conditions as a key component in species domestication processes. The aim of the present study is to evaluate the ability of a selected high-growth genotype of 300 days post-hatch European sea bass (Dicentrarchus labrax) juveniles to use different functional additives as endogenous antioxidant capacity and stress resistance boosters when supplemented in low fish meal (FM) and fish oil (FO) diets. Three isoenergetic and isonitrogenous diets (10% FM/6% FO) were supplemented with 200 ppm of a blend of garlic and Labiatae plant oils (PHYTO0.02), 1000 ppm of a mixture of citrus flavonoids and Asteraceae and Labiatae plant essential oils (PHYTO0.1) or 5000 ppm of galactomannan-oligosaccharides (GMOS0.5). A reference diet was void of supplementation. The fish were fed the experimental diets for 72 days and subjected to a H2O2 exposure oxidative stress challenge. The fish stress response was evaluated through measuring the circulating plasma cortisol levels and the fish gill antioxidant response by the relative gene expression analysis of nfΚß2, il-1b, hif-1a, nd5, cyb, cox, sod, cat, gpx, tnf-1α and caspase 9. After the oxidative stress challenge, the genotype origin determined the capacity of the recovery of basal cortisol levels after an acute stress response, presenting GS fish with a better pattern of recovery. All functional diets induced a significant upregulation of cat gill gene expression levels compared to fish fed the control diet, regardless of the fish genotype. Altogether, suggesting an increased capacity of the growth selected European sea bass genotype to cope with the potential negative side-effects associated to an H2O2 bath exposure.

Genomic Selection in Aquaculture Species.

To date, genomic prediction has been conducted in about 20 aquaculture species, with a preference for intra-family genomic selection (GS). For every trait under GS, the increase in accuracy obtained by genomic estimated breeding values instead of classical pedigree-based estimation of breeding values is very important in aquaculture species ranging from 15% to 89% for growth traits, and from 0% to 567% for disease resistance. Although the implementation of GS in aquaculture is of little additional investment in breeding programs already implementing sib testing on pedigree, the deployment of GS remains sparse, but could be boosted by adaptation of cost-effective imputation from low-density panels. Moreover, GS could help to anticipate the effect of climate change by improving sustainability-related traits such as production yield (e.g., carcass or fillet yields), feed efficiency or disease resistance, and by improving resistance to environmental variation (tolerance to temperature or salinity variation). This chapter synthesized the literature in applications of GS in finfish, crustaceans and molluscs aquaculture in the present and future breeding programs.

Adaptation to drought in two wild tomato species: the evolution of the Asr gene family.

Wild tomato species are a valuable system in which to study local adaptation to drought: they grow in diverse environments ranging from mesic to extremely arid conditions. Here, we investigate the evolution of members of the Asr (ABA/water stress/ripening induced) gene family, which have been reported to be involved in the water stress response. We analysed molecular variation in the Asr gene family in populations of two closely related species, Solanum chilense and Solanum peruvianum. We concluded that Asr1 has evolved under strong purifying selection. In contrast to previous reports, we did not detect evidence for positive selection at Asr2. However, Asr4 shows patterns consistent with local adaptation in an S. chilense population that lives in an extremely dry environment. We also discovered a new member of the gene family, Asr5. Our results show that the Asr genes constitute a dynamic gene family and provide an excellent example of tandemly arrayed genes that are of importance in adaptation. Taking the potential distribution of the species into account, it appears that S. peruvianum can cope with a great variety of environmental conditions without undergoing local adaptation, whereas S. chilense undergoes local adaptation more frequently.

Molecular data reveal isolation by distance and past population expansion for the shea tree (Vitellaria paradoxa C.F. Gaertn) in West Africa.

While the genetic structure of many tree species in temperate, American and Asian regions is largely explained by climatic oscillations and subsequent habitat contractions and expansions, little is known about Africa. We investigated the genetic diversity and structure of shea tree (Vitellaria paradoxa,) in Western Africa, an economically important tree species in the Sudano-Sahelian zone. Eleven nuclear microsatellites (nuc) were used to genotype 673 trees selected in 38 populations. They revealed moderate to high within-population diversity: allelic richness ranged from R(nuc) = 3.99 to 5.63. This diversity was evenly distributed across West Africa. Populations were weakly differentiated (F(STnuc) = 0.085; P < 0.0001) and a pattern of isolation by distance was noted. No phylogeographic signal could be detected across the studied sample. Additionally, two chloroplast microsatellite loci, leading to 11 chlorotypes, were used to analyse a sub-set of 370 individuals. Some variation in chloroplast allelic richness among populations could be detected (R(cp) = 0.00 to 4.36), but these differences were not significant. No trend with latitude and longitude were observed. Differentiation was marked (G(STcp) = 0.553; P < 0.0001), but without a significant phylogeographical signal. Population expansion was detected considering the total population using approximate Bayesian computation (nuclear microsatellites) and mismatch distribution (chloroplast microsatellites) methods. This expansion signal and the isolation by distance pattern could be linked to the past climatic conditions in West Africa during the Pleistocene and Holocene which should have been favourable to shea tree development. In addition, human activities through agroforestry and domestication (started 10,000 bp) have probably enhanced gene flow and population expansion.

Case report of osteomyelitis of the mandible in osteopetrosis and management considerations.

INTRODUCTION AND IMPORTANCE: Osteopetrosis is a poorly known and probably underdiagnosed pathology. It is caused by various genetic abnormalities resulting in osteoclast dysfunction. Functional and aesthetic consequences have a major impact on the patient's quality of life. Ten percent of osteopetrosis cases develop osteomyelitis that usually involves the mandible. Management of this complication remains complex and often unsatisfactory. CASE PRESENTATION: We report a case of a 62-year-old woman with osteopetrosis, complicated by mandibular osteomyelitis with intra-oral bone exposure and submental fistulas. Management was performed with antibiotic therapy and surgical necrotic resection. This cured the fistulas but the bone exposure persisted. DISCUSSION: This case report highlights the difficulty of achieving complete healing of osteomyelitis in osteopetrosis. Antibiotic therapy, surgical management, or even hyperbaric oxygen therapy are required, but must be adapted to the case. A free flap procedure is undesirable but, when it is necessary, a bone marrow transplant could be considered to restore osteoclast function. CONCLUSION: The management of mandibular osteomyelitis in patients with osteopetrosis must adapt to the situation and severity. To avoid most cases of osteomyelitic complications in patients suffering from osteopetrosis, we propose that a preventive strategy of better dental care should be considered.

Parental selection for growth and early-life low stocking density increase the female-to-male ratio in European sea bass.

In European sea bass (Dicentrarchus labrax), as in many other fish species, temperature is known to influence the sex of individuals, with more males produced at relatively high temperatures. It is however unclear to what extent growth or stress are involved in such a process, since temperature is known to influence both growth rate and cortisol production. Here, we designed an experiment aiming at reducing stress and affecting early growth rate. We exposed larvae and juveniles originating from both captive and wild parents to three different treatments: low stocking density, food supplemented with tryptophan and a control. Low stocking density and tryptophan treatment respectively increased and decreased early growth rate. Each treatment influenced the stress response depending on the developmental stage, although no clear pattern regarding the whole-body cortisol concentration was found. During sex differentiation, fish in the low-density treatment exhibited lower expression of gr1, gr2, mr, and crf in the hypothalamus when compared to the control group. Fish fed tryptophan displayed lower crf in the hypothalamus and higher level of serotonin in the telencephalon compared to controls. Overall, fish kept at low density produced significantly more females than both control and fish fed tryptophan. Parents that have been selected for growth for three generations also produced significantly more females than parents of wild origin. Our findings did not allow to detect a clear effect of stress at the group level and rather point out a key role of early sexually dimorphic growth rate in sex determination.

Sex dimorphism in European sea bass (Dicentrarchus labrax L.): New insights into sex-related growth patterns during very early life stages.

The European sea bass (Dicentrarchus labrax) exhibits female-biased sexual size dimorphism (SSD) early in development. New tagging techniques provide the opportunity to monitor individual sex-related growth during the post-larval and juvenile stages. We produced an experimental population through artificial fertilization and followed a rearing-temperature protocol (~16°C from hatching to 112 days post-hatching, dph; ~20°C from 117 to 358 dph) targeting a roughly balanced sex ratio. The fish were tagged with microchips between 61 and 96 dph in five tagging trials of 50 fish each; individual standard length (SL) was recorded through repeated biometric measurements performed between 83 to 110 dph via image analyses. Body weight (BW) was modelled using the traits measured on the digital pictures (i.e. area, perimeter and volume). At 117 dph, the fish were tagged with microtags and regularly measured for SL and BW until 335 dph. The experiment ended at 358 dph with the sexing of the fish. The sex-ratio at the end of the experiment was significantly in favor of the females (65.6% vs. 34.4%). The females were significantly longer and heavier than the males from 103 dph (~30 mm SL, ~0.44 g BW) to 165 dph, but the modeling of the growth curves suggests that differences in size already existed at 83 dph. A significant difference in the daily growth coefficient (DGC) was observed only between 96 and 103 dph, suggesting a physiological or biological change occurring during this period. The female-biased SSD pattern in European sea bass is thus strongly influenced by very early growth differences between sexes, as already shown in previous studies, and in any case long before gonadal sex differentiation has been started, and thus probably before sex has been determined. This leads to the hypothesis that early growth might be a cause rather than a consequence of sex differentiation in sea bass.

Optimization of Genomic Selection to Improve Disease Resistance in Two Marine Fishes, the European Sea Bass (Dicentrarchus labrax) and the Gilthead Sea Bream (Sparus aurata).

Disease outbreaks are a major threat to the aquaculture industry, and can be controlled by selective breeding. With the development of high-throughput genotyping technologies, genomic selection may become accessible even in minor species. Training population size and marker density are among the main drivers of the prediction accuracy, which both have a high impact on the cost of genomic selection. In this study, we assessed the impact of training population size as well as marker density on the prediction accuracy of disease resistance traits in European sea bass (Dicentrarchus labrax) and gilthead sea bream (Sparus aurata). We performed a challenge to nervous necrosis virus (NNV) in two sea bass cohorts, a challenge to Vibrio harveyi in one sea bass cohort and a challenge to Photobacterium damselae subsp. piscicida in one sea bream cohort. Challenged individuals were genotyped on 57K-60K SNP chips. Markers were sampled to design virtual SNP chips of 1K, 3K, 6K, and 10K markers. Similarly, challenged individuals were randomly sampled to vary training population size from 50 to 800 individuals. The accuracy of genomic-based (GBLUP model) and pedigree-based estimated breeding values (EBV) (PBLUP model) was computed for each training population size using Monte-Carlo cross-validation. Genomic-based breeding values were also computed using the virtual chips to study the effect of marker density. For resistance to Viral Nervous Necrosis (VNN), as one major QTL was detected, the opportunity of marker-assisted selection was investigated by adding a QTL effect in both genomic and pedigree prediction models. As training population size increased, accuracy increased to reach values in range of 0.51-0.65 for full density chips. The accuracy could still increase with more individuals in the training population as the accuracy plateau was not reached. When using only the 6K density chip, accuracy reached at least 90% of that obtained with the full density chip. Adding the QTL effect increased the accuracy of the PBLUP model to values higher than the GBLUP model without the QTL effect. This work sets a framework for the practical implementation of genomic selection to improve the resistance to major diseases in European sea bass and gilthead sea bream.

Corrigendum: Optimization of Genomic Selection to Improve Disease Resistance in Two Marine Fishes, The European Sea Bass (Dicentrarchus labrax) and the Gilthead Sea Bream (Sparus aurata).

[This corrects the article DOI: 10.3389/fgene.2021.665920.].

An investigation of links between metabolic rate and feed efficiency in European sea bass Dicentrarchus labrax.

Feed efficiency (FE) is the amount of body weight gain for a given feed intake. Improving FE through selective breeding is key for sustainable finfish aquaculture but its evaluation at individual level is technically challenging. We therefore investigated whether individual routine metabolic rate (RMR) was a predictor of individual FE in the European sea bass Dicentrarchus labrax, a major species in European mariculture. The European sea bass has three genetically distinct populations across its geographical range, namely Atlantic (AT), West Mediterranean (WM), and East Mediterranean (EM). We compared FE and RMR of fish from these three populations at 18 or 24 °C. We held 200 fish (62 AT, 66 WM, and 72 EM) in individual aquaria and fed them from ad libitum down to fasting. FI was assessed for an ad libitum feeding rate and for a fixed restricted ration (1% of metabolic body weight·day-1, with metabolic body weight = body weight0.8). After being refed 12 wk in a common tank, individual RMR was measured over 36 h by intermittent flow respirometry. There was a significant effect of temperature whereby fish at 18 °C had greater mean FE (P < 0.05) and lower RMR (P < 0.001). There was also a significant effect of population, where AT fish had lower FE (P < 0.05) and greater RMR (P < 0.001) than WM and EM, at both temperatures. Despite these differences in temperature and population means, individual FE and RMR were not significantly correlated (P > 0.05). Therefore, although the results provide evidence of an association between metabolic rate and FE, RMR was not a predictor of individual FE, for reasons that require further investigation.

The contribution of ancient admixture to reproductive isolation between European sea bass lineages.

Understanding how new species arise through the progressive establishment of reproductive isolation (RI) barriers between diverging populations is a major goal in Evolutionary Biology. An important result of speciation genomics studies is that genomic regions involved in RI frequently harbor anciently diverged haplotypes that predate the reconstructed history of species divergence. The possible origins of these old alleles remain much debated, as they relate to contrasting mechanisms of speciation that are not yet fully understood. In the European sea bass (Dicentrarchus labrax), the genomic regions involved in RI between Atlantic and Mediterranean lineages are enriched for anciently diverged alleles of unknown origin. Here, we used haplotype-resolved whole-genome sequences to test whether divergent haplotypes could have originated from a closely related species, the spotted sea bass (Dicentrarchus punctatus). We found that an ancient admixture event between D. labrax and D. punctatus is responsible for the presence of shared derived alleles that segregate at low frequencies in both lineages of D. labrax. An exception to this was found within regions involved in RI between the two D. labrax lineages. In those regions, archaic tracts originating from D. punctatus locally reached high frequencies or even fixation in Atlantic genomes but were almost absent in the Mediterranean. We showed that the ancient admixture event most likely occurred between D. punctatus and the D. labrax Atlantic lineage, while Atlantic and Mediterranean D. labrax lineages were experiencing allopatric isolation. Our results suggest that local adaptive introgression and/or the resolution of genomic conflicts provoked by ancient admixture have probably contributed to the establishment of RI between the two D. labrax lineages.

APIS: An auto-adaptive parentage inference software that tolerates missing parents.

In the context of parentage assignment using genomic markers, key issues are genotyping errors and an absence of parent genotypes because of sampling, traceability or genotyping problems. Most likelihood-based parentage assignment software programs require a priori estimates of genotyping errors and the proportion of missing parents to set up meaningful assignment decision rules. We present here the R package APIS, which can assign offspring to their parents without any prior information other than the offspring and parental genotypes, and a user-defined, acceptable error rate among assigned offspring. Assignment decision rules use the distributions of average Mendelian transmission probabilities, which enable estimates of the proportion of offspring with missing parental genotypes. APIS has been compared to other software (CERVUS, VITASSIGN), on a real European seabass (Dicentrarchus labrax) single nucleotide polymorphism data set. The type I error rate (false positives) was lower with APIS than with other software, especially when parental genotypes were missing, but the true positive rate was also lower, except when the theoretical exclusion power reached 0.99999. In general, APIS provided assignments that satisfied the user-set acceptable error rate of 1% or 5%, even when tested on simulated data with high genotyping error rates (1% or 3%) and up to 50% missing sires. Because it uses the observed distribution of Mendelian transmission probabilities, APIS is best suited to assigning parentage when numerous offspring (>200) are genotyped. We have demonstrated that APIS is an easy-to-use and reliable software for parentage assignment, even when up to 50% of sires are missing.