Search

Early development stages in marine bivalve are critical periods where larvae transition from pelagic free‐life to sessile mature individuals. The successive metamorphosis requires the expression of key genes, the functions of which might be under high selective pressure, hence understanding larval development represents key knowledge for both fundamental and applied research. Phenotypic larvae development is well known, but the underlying molecular mechanisms such as associated gene expression dynamic and molecular cross‐talks remains poorly described for several nonmodel species, such as P. margaritifera. We designed a whole transcriptome RNA‐sequencing analysis to describe such gene expression dynamics following four larval developmental stages: d‐shape, Veliger, Umbo and Eye‐spot. Larval gene expression and annotated functions drastically diverge. Metabolic function (gene expression related to lipid, amino acid and carbohydrate use) is highly upregulated in the first development stages, with increasing demand from d‐shape to umbo. Morphogenesis and larval transition are partly ordered by Thyroid hormones and Wnt signaling. While larvae shells show some similar characteristic to adult shells, the cause of initialization of biomineralization differ from the one found in adults. The present study provides a global overview of Pinctada margaritifera larval stages transitioning through gene expression dynamics, molecular mechanisms and ontogeny of biomineralization, immune system, and sensory perception processes.

Tahitian pearls, artificially cultivated from the black-lipped pearl oyster Pinctada margaritifera, are renowned for their unique color and large size, making the pearl industry vital for the French Polynesian economy. Understanding the mechanisms of pearl formation is essential for enabling quality and sustainable production. In this paper, we explore the process of pearl formation by studying pearl rotation. Here we show, using a deep convolutional neural network, a direct link between the rotation of the pearl during its formation in the oyster and its final shape. We propose a new method for non-invasive pearl monitoring and a model for predicting the final shape of the pearl from rotation data with 81.9% accuracy. These novel resources provide a fresh perspective to study and enhance our comprehension of the overall mechanism of pearl formation, with potential long-term applications for improving pearl production and quality control in the industry.

Identifying the molecular mechanisms involved in rapid adaptation to novel environments and determining their predictability, are central questions in evolutionary biology and pressing issues due to rapid global changes. Complementary to genetic responses to selection, faster epigenetic variations such as modifications of DNA methylation may play a substantial role in rapid adaptation. In the context of rampant urbanization, joint examinations of genomic and epigenomic mechanisms are still lacking. Here, we investigated genomic (SNP) and epigenomic (CpG methylation) responses to urban life in a passerine bird, the Great tit (Parus major). To test whether urban evolution is predictable (i.e parallel) or involves mostly non-parallel molecular processes among cities, we analysed both SNP and CpG methylation variations across three distinct pairs of city and forest Great tit populations in Europe. Our analyses reveal a polygenic response to urban life, with both many genes putatively under weak divergent selection and multiple differentially methylated regions (DMRs) between forest and city great tits. DMRs mainly overlapped transcription start sites and promotor regions, suggesting their importance in modulating gene expression. Both genomic and epigenomic outliers were found in genomic regions enriched for genes with biological functions related to the nervous system, immunity, or behavioural, hormonal and stress responses. Interestingly, comparisons across the three pairs of city-forest populations suggested little parallelism in both genetic and epigenetic responses. Our results confirm, at both the genetic and epigenetic levels, hypotheses of polygenic and largely non-parallel mechanisms of rapid adaptation in novel environments such as urbanized areas.

Human activities and resource exploitation led to a massive decline of wild salmonid populations, consequently numerous conservation programs have been developed to supplement wild populations. However, many studies documented reduced fitness of hatchery-born relative to wild fish. Here, by using both RNA sequencing and Whole Genome Bisulfite Sequencing (WGBS), we show that of hatchery and wild born adult Coho Salmon (Oncorhynchus kisutch) originating from two previously studied river systems, early-life hatchery rearing environment induced significant and parallel gene expression differentiation is maintained until Coho come back to their natal river for reproduction. A total of 3,643 genes differentially expressed and 859 co-expressed genes were down-regulated in parallel in hatchery born fish from both rivers relative to their wild congeners. Among those genes, 26 displayed a significant relationship between gene expression and the median gene body methylation and 669 single CpG displayed a significant correlation between methylation level and the associated gene expression. The link between methylation and gene expression was weak suggesting that DNA methylation is not the only player in mediating hatchery-related expression differences. Yet, significant gene expression differentiation was observed despite 18 month spent in a common environment (i.e. the sea). Finally, the differentiation is observed in parallel in two different river system, highlighting the fact that early life environment may account for at least some of the reduced fitness of the hatchery salmon in the wild. These results illustrate the relevance and importance of considering both epigenome and transcriptome to evaluate the costs and benefits of large-scale supplementation programs.

Studies on the impact of ocean acidification on marine organisms involve exposing organisms to future acidification scenarios which has limited relevance for coastal calcifiers living in a mosaic of habitats. Identification of tipping points beyond which detrimental effects are observed is a widely generalizable proxy of acidification susceptibility at the populational level. This approach is limited to a handful of studies that focus on only a few macro-physiological traits, thus overlooking the whole organism response. Here we develop a framework to analyze the broad macro-physiological and molecular responses over a wide pH range in juvenile oyster. We identify low tipping points for physiological traits at pH 7.3-6.9 that coincide with a major reshuffling in membrane lipids and transcriptome. In contrast, a drop in pH affects shell parameters above tipping points, likely impacting animal fitness. These findings were made possible by the development of an innovative methodology to synthesize and identify the main patterns of variations in large -omic datasets, fitting them to pH and identifying molecular tipping-points. We propose the broad application of our framework to the assessment of effects of global change on other organisms.

Ocean warming challenges marine organisms' resilience, especially for species experiencing temperatures close to their upper thermal limits. A potential increase in thermal tolerance might significantly reduce the risk of population decline, which is intrinsically linked to variability in local habitat temperatures. Our goal was to assess the plastic and genetic potential of response to elevated temperatures in a tropical bivalve model, Pinctada margaritifera. We benefit from two ecotypes for which local environmental conditions are characterized by either large diurnal variations in the tide pools (Marquesas archipelago) or low mean temperature with stable to moderate seasonal variations (Gambier archipelago). We explored the physiological basis of individual responses to elevated temperature, genetic divergence as well as plasticity and acclimation by combining lipidomic and transcriptomic approaches. We show that P. margaritifera has certain capacities to adjust to long-term elevated temperatures that was thus far largely underestimated. Genetic variation across populations overlaps with gene expression and involves the mitochondrial respiration machinery, a central physiological process that contributes to species thermal sensitivity and their distribution ranges. Our results present evidence for acclimation potential in P. margaritifera and urge for longer term studies to assess populations resilience in the face of climate change.

Background Tenacibaculum maritimum is a fish pathogen known for causing serious damage to a broad range of wild and farmed marine fish populations worldwide. The recently sequenced genome of T. maritimum strain NCIMB 2154T provided unprecedented information on the possible molecular mechanisms involved in the virulence of this species. However, little is known about the dynamic of infection in vivo, and information is lacking on both the intrinsic host response (gene expression) and its associated microbiota. Here, we applied complementary omic approaches, including dual RNAseq and 16S rRNA gene metabarcoding sequencing using Nanopore and short-read Illumina technologies to unravel the host–pathogen interplay in an experimental infection system using the tropical fish Platax orbicularis as model. Results We showed that the infection of the host is characterised by an enhancement of functions associated with antibiotic and glucans catabolism functions but a reduction of sulfate assimilation process in T. maritimum. The fish host concurrently displays a large panel of immune effectors, notably involving innate response and triggering acute inflammatory response. In addition, our results suggest that fish activate an adaptive immune response visible through the stimulation of T-helper cells, Th17, with congruent reduction of Th2 and T-regulatory cells. Fish were, however, largely sensitive to infection, and less than 25% survived after 96 hpi. These surviving fish showed no evidence of stress (cortisol levels) or significant difference in microbiome diversity compared with controls at the same sampling time. The presence of T. maritimum in resistant fish skin and the total absence of any skin lesions suggest that these fish did not escape contact with the pathogen, but rather that some mechanisms prevented pathogens entry. In resistant individuals, we detected up-regulation of specific immune-related genes differentiating resistant individuals from controls at 96 hpi

Human activities induce direct or indirect selection pressure on natural population and may ultimately affect population’s integrity. While numerous conservation programs aimed to minimize human‐induced genomic variation, human‐induced environmental variation may generate epigenomic variation potentially affecting fitness through phenotypic modifications. Major questions remain pertaining to how much epigenomic variation arises from environmental heterogeneity, whether this variation can persist throughout life, and whether it can be transmitted across generations. We performed whole genome bisulfite sequencing (WGBS) on the sperm of genetically indistinguishable hatchery and wild born migrating adults of Coho salmon (Oncorhynchus kisutch) from two geographically distant rivers at different epigenome scales. Our results showed that coupling WGBS with fine scale analyses (local and chromosomal) allows the detection of parallel early‐life hatchery‐induced epimarks that differentiate wild from hatchery‐reared salmon. Four chromosomes and 183 differentially methylated regions (DMRs) displayed a significant signal of methylation differentiation between hatchery and wild born Coho salmon. Moreover, those early‐life epimarks persisted in germ‐line cells despite about 1.5 year spent in the ocean following release from hatchery, opening the possibility for transgenerational inheritance. Our results strengthen the hypothesis that epigenomic modifications environmentally‐induced during early‐life development persist in germ cells of adults until reproduction, which could potentially impact their fitness.

Genotype-by-environment (GxE) interactions are non-parallel reaction norms among individuals with different genotypes in response to different environmental conditions. GxE interactions are an extension of phenotypic plasticity and consequently studying such interactions improves our ability to predict effects of different environments on phenotype as well as the fitness of genetically distinct organisms and their capacity to interact with ecosystems. Growth hormone transgenic coho salmon grow much faster than non-transgenics when raised in tank environments, but show little difference in growth when reared in nature-like streams. We used this model system to evaluate potential mechanisms underlying this growth rate GxE interaction, performing RNA-seq to measure gene transcription and whole-genome bisulfite sequencing to measure gene methylation in liver tissue. Gene ontology (GO) term analysis revealed stress as an important biological process potentially influencing growth rate GxE interactions. While few genes with transcription differences also had methylation differences, in promoter or gene regions, many genes were differentially methylated between tank and stream environments. A GO term analysis of differentially methylated genes between tank and stream environments revealed increased methylation in the stream environment of more than 95% of the differentially methylated genes, many with biological processes unrelated to liver function. The lower nutritional condition of the stream environment may cause increased negative regulation of genes less vital for liver tissue function than when fish are reared in tanks with unlimited food availability. These data show a large effect of rearing environment both on gene expression and methylation, but it is less clear that the detected epigenetic marks are responsible for the observed altered growth and physiological responses.

Identifying the molecular mechanisms involved in rapid adaptation to novel environments and determining their predictability are central questions in Evolutionary Biology and pressing issues due to rapid global changes. Complementary to genetic responses to selection, faster epigenetic variations such as modifications of DNA methylation may play a substantial role in rapid adaptation. In the context of rampant urbanization, joint examinations of genomic and epigenomic mechanisms are still lacking. Here, we investigated genomic (SNP) and epigenomic (CpG methylation) responses to urban life in a passerine bird, the Great tit (Parus major). To test whether urban evolution is predictable (i.eparallel) or involves mostly non-parallel molecular processes among cities, we analysed three distinct pairs of city and forest Great tit populations across Europe. Results reveal a polygenic response to urban life, with both many genes putatively under weak divergent selection and multiple differentially methylated regions (DMRs) between forest and city great tits. DMRs mainly overlapped transcription start sites and promotor regions, suggesting their importance in the modulation gene expression. Both genomic and epigenomic outliers were found in genomic regions enriched for genes with biological functions related to nervous system, immunity, behaviour, hormonal and stress responses. Interestingly, comparisons across the three pairs of city-forest populations suggested little parallelism in both genetic and epigenetic responses. Our results confirm, at both the genetic and epigenetic levels, hypotheses of polygenic and largely non-parallel mechanisms of rapid adaptation in new environments such as urbanized areas.

Background Albino mutations are commonly observed in the animal kingdom, including in bivalves. In the black-lipped pearl oyster Pinctada margaritifera, albino specimens are characterized by total or partial absence of colouration resulting in typical white shell phenotype expression. The relationship of shell colour with resulting cultured pearl colour is of great economic interest in P. margaritifera, on which a pearl industry is based. Hence, the albino phenotype provides a useful way to examine the molecular mechanisms underlying pigmentation. Results Whole transcriptome RNA-sequencing analysis comparing albino and black wild-type phenotypes at three stages over the culture cycle of P. margaritifera revealed a total of 1606, 798 and 187 differentially expressed genes in whole juvenile, adult mantle and pearl sac tissue, respectively. These genes were found to be involved in five main molecular pathways, tightly linked to known pigmentation pathways: melanogenesis, calcium signalling pathway, Notch signalling pathway, pigment transport and biomineralization. Additionally, significant phenotype-associated SNPs were selected (N = 159), including two located in the Pif biomineralization gene, which codes for nacre formation. Interestingly, significantly different transcript splicing was detected between juvenile (N = 1366) and adult mantle tissue (N = 313) in, e.g., the tyrosinase Tyr-1 gene, which showed more complex regulation in mantle, and the Notch1 encoding gene, which was upregulated in albino juveniles. Conclusion This multiple RNA-seq approach provided new knowledge about genes associated with the P. margaritifera albino phenotype, highlighting: 1) new molecular pathways, such as the Notch signalling pathway in pigmentation, 2) associated SNP markers with biomineraliszation gene of interest like Pif for marker-assisted selection and prevention of inbreeding, and 3) alternative gene splicing for melanin biosynthesis implicating tyrosinase.

Understanding the genomic basis of adaptative intraspecific phenotypic variation is a central goal in conservation genetics and evolutionary biology. Lake trout (Salvelinus namaycush) are an excellent species for addressing the genetic basis for adaptive variation because they express a striking degree of ecophenotypic variation across their range; however, necessary genomic resources are lacking. Here we utilize recently-developed analytical methods and sequencing technologies to (1) construct a high-density linkage and centromere map for lake trout, (2) identify loci underlying variation in traits that differentiate lake trout ecophenotypes and populations, (3) determine the location of the lake trout sex determination locus, and (4) identify chromosomal homologies between lake trout and other salmonids of varying divergence. The resulting linkage map contains 15,740 single nucleotide polymorphisms (SNPs) mapped to 42 linkage groups, likely representing the 42 lake trout chromosomes. Female and male linkage group lengths ranged from 43.07 to 134.64 centimorgans, and 1.97 to 92.87 centimorgans, respectively. We improved the map by determining coordinates for 41 of 42 centromeres, resulting in a map with 8 metacentric chromosomes and 34 acrocentric or telocentric chromosomes. We use the map to localize the sex determination locus and multiple quantitative trait loci (QTL) associated with intraspecific phenotypic divergence including traits related to growth and body condition, patterns of skin pigmentation, and two composite geomorphometric variables quantifying body shape. Two QTL for the presence of vermiculations and spots mapped with high certainty to an arm of linkage group Sna3, growth related traits mapped to two QTL on linkage groups Sna1 and Sna12, and putative body shape QTL were detected on six separate linkage groups. The sex determination locus was mapped to Sna4 with high confidence. Synteny analysis revealed that lake trout and congener Arctic char (S

Seawater temperature rise in French Polynesia has repeatedly resulted in the bleaching of corals and giant clams. Because giant clams possess distinctive ectosymbiotic features, they represent a unique and powerful model for comparing molecular pathways involved in 1) maintenance of symbiosis and 2) acquisition of thermo‐tolerance among coral reef organisms. Herein, we explored the physiological and transcriptomic responses of the clam hosts and their photosynthetically active symbionts over a 65‐day experiment in which clams were exposed to either normal or environmentally relevant elevated seawater temperatures. Additionally, we used metabarcoding data coupled with in situ sampling/survey data to explore the relative importance of holobiont adaptation (i.e., a symbiont community shift) versus acclimation (i.e., physiological changes at the molecular level) in the clams’ responses to environmental change. We finally compared transcriptomic data to publicly available genomic datasets for Symbiodiniaceae dinoflagellates (both cultured and in hospite with the coral Pocillopora damicornis) to better tease apart the responses of both hosts and specific symbiont genotypes in this mutualistic association. Gene module preservation analysis revealed that the function of the symbionts’ photosystem II was impaired at high temperature, and this response was also found across all holobionts and Symbiodiniaceae lineages examined. Similarly, epigenetic modulation appeared to be a key response mechanism for symbionts in hospite with giant clams exposed to high temperatures, and such modulation was able to distinguish thermo‐tolerant from thermo‐sensitive Cladocopium goreaui ecotypes; epigenetic processes may, then, represent a promising research avenue for those interested in coral reef conservation in this era of changing global climate.

Understanding the genetic and evolutionary impacts of stocking on wild fish populations has long been of interest as negative consequences such as reduced fitness and loss of genetic diversity are commonly reported outcomes. In an attempt to sustain a fishery, managers implemented nearly five decades of extensive stocking of over a million Muskellunge (Esox masquinongy), a native species in the Lower St. Lawrence River (Québec, Canada). We investigated the effect of this stocking on population genetic structure and allelic diversity in the St. Lawrence River in addition to tributaries and several stocked inland lakes. Using Genotype‐By‐Sequencing (GBS), we genotyped 643 individuals representing 22 locations and combined this information with forward simulations to investigate the genetic consequences of long‐term stocking. Individuals native to the St. Lawrence watershed were genetically differentiated from stocking sources and tributaries and inland lakes were naturally differentiated from the main river. Empirical data and simulations within the St. Lawrence River revealed weak stocking effects on admixture patterns. Our data suggest that the genetic structure associated with stocked fish was diluted into its relatively large effective population size. This interpretation is also consistent with a hypothesis that selection against introgression was in operation and relatively efficient within the large St. Lawrence River system. In contrast smaller populations from adjacent tributaries and lakes displayed greater stocking‐related admixture that resulted in comparatively higher heterozygosity than the St. Lawrence. Finally, individuals from inland lakes that were established by stocking maintained a close affinity with their source populations. This study illustrated a benefit of combining extensive genomic data with forward simulations for improved inference regarding population‐level genetic effects of long‐term stocking, and its relevance for fishery management

Background Cultured pearls are unique gems produced by living organisms, mainly molluscs of the Pinctada genus, through the biomineralization properties of pearl sac tissue. Improvement of P. margaritifera pearl quality is one of the biggest challenges that Polynesian research has faced to date. To achieve this goal, a better understanding of the complex mechanisms related to nacre and pearl formation is essential and can now be approached through the use of massive parallel sequencing technologies. The aim of this study was to use RNA-seq to compare whole transcriptome expression of pearl sacs that had producing pearls with high and low quality. For this purpose, a comprehensive reference transcriptome of P. margaritifera was built based on multi-tissue sampling (mantle, gonad, whole animal), including different living stages (juvenile, adults) and phenotypes (colour morphotypes, sex). Results Strikingly, few genes were found to be up-regulated for high quality pearls (n = 16) compared to the up-regulated genes in low quality pearls (n = 246). Biomineralization genes up-regulated in low quality pearls were specific to prismatic and prism-nacre layers. Alternative splicing was further identified in several key biomineralization genes based on a recent P. margaritifera draft genome. Conclusion This study lifts the veil on the multi-level regulation of biomineralization genes 38 associated with pearl quality determination.

The genomic shock hypothesis stipulates that the stress associated with divergent genome admixture can cause transposable element (TE) derepression, which could act as a postzygotic isolation mechanism. TEs affect gene structure, expression patterns, and chromosome organization and may have deleterious consequences when released. For these reasons, they are silenced by heterochromatin formation, which includes DNA methylation. Here, we show that a significant proportion of TEs are differentially methylated between the “dwarf” (limnetic) and the “normal” (benthic) whitefish, two nascent species that diverged some 15,000 generations ago within the Coregonus clupeaformis species complex. Moreover, TEs are overrepresented among loci that were demethylated in hybrids, indicative of their transcriptional derepression. These results are consistent with earlier studies in this system that revealed TE transcriptional derepression causes abnormal embryonic development and death of hybrids. Hence, this supports a role of DNA methylation reprogramming and TE derepression in postzygotic isolation of nascent animal species.

Identifying genes and pathways involved in domestication is critical to understand how species change in response to human-induced selection pressures, such as increased temperatures. Given the profound influence of temperature on fish metabolism and organismal performance, a comparison of how temperature affects wild and domestic strains of snapper is an important question to address. We experimentally manipulated temperature conditions for F1-hatchery and wild Australasian snapper (Chrysophrys auratus) for 18 days to mimic seasonal extremes and measured differences in growth, white muscle RNA transcription and hematological parameters. Over 2.2 Gb paired-end reads were assembled de novo for a total set of 33,017 transcripts (N50 = 2,804). We found pronounced growth and gene expression differences between wild and domesticated individuals related to global developmental and immune pathways. Temperature-modulated growth responses were linked to major pathways affecting metabolism, cell regulation and signaling. This study is the first step toward gaining an understanding of the changes occurring in the early stages of domestication, and the mechanisms underlying thermal adaptation and associated growth in poikilothermic vertebrates. Our study further provides the first transcriptome resources for studying biological questions in this non-model fish species.

In fisheries management, intensive stocking programs are commonly used to enhance population abundance and maintain stock productivity. However, such practices are increasingly raising concerns as multiple studies documented adverse genetic and evolutionary impacts of stocking on wild populations. Improvement of stocking management relies on a better understanding of the dynamic of introgressive hybridization between wild and domestic population and on assessment of the genetic state of wild populations after stocking cessation. In Québec, Canada, over five million captive-reared Brook Charr (Salvelinus fontinalis) are stocked every year to support recreational fishing activities. Here, we investigated how variation in stocking history and environmental variables, including water temperature, pH, and dissolved oxygen, may influence the impact of stocking practices on the genetic integrity of wild Brook Charr populations. We collected DNA samples (n = 862, average of 30 individuals per lake) from 29 lakes that underwent different stocking intensity through time and also collected environmental parameters for each sampled lake. An average of 4,580 high-quality filtered SNPs was obtained for each population using genotyping by sequencing (GBS), which were then used to quantify the mean domestic membership of each sampled population. An exhaustive process of model selection was conducted to obtain a best-fitted model that explained 56% of the variance observed in mean domestic genetic membership. The number of years since the mean year of stocking was the best explanatory variable to predict variation in mean domestic genetic membership whereas environmental characteristics had little influence on observed patterns of admixture. Our model predictions also revealed that each sampled wild population could potentially return to a wild genetic state (absence of domestic genetic background) after stocking cessation. Overall, our study provides new insights on factors determi

Migration is a ubiquitous life history trait with profound evolutionary and ecological consequences. Recent developments in telemetry and genomics, when combined, can bring significant insights on the migratory ecology of non-model organisms in the wild. Here, we used this integrative approach to document dispersal, gene flow and potential for local adaptation in anadromous Arctic Char from six rivers in the Canadian Arctic. Acoustic telemetry data from 124 tracked individuals indicated asymmetric dispersal, with a large proportion of fish (72%) tagged in three different rivers migrating up the same short river in the fall. Population genomics data from 6,136 SNP markers revealed weak, albeit significant, population differentiation (average pairwise FST = 0.011) and asymmetric dispersal was also revealed by population assignments. Approximate Bayesian Computation simulations suggested the presence of asymmetric gene flow, although in the opposite direction to that observed from the telemetry data, suggesting that dispersal does not necessarily lead to gene flow. These observations suggested that Arctic Char home to their natal river to spawn, but may overwinter in rivers with the shortest migratory route to minimize the costs of migration in non-breeding years. Genome scans and genetic-environment associations identified 90 outlier markers putatively under selection, 23 of which were in or near a gene. Of these, at least four were involved in muscle and cardiac function, consistent with the hypothesis that migratory harshness could drive local adaptation. Our study illustrates the power of integrating genomics and telemetry to study migrations in non-model organisms in logistically challenging environments such as the Arctic.

Using massively parallel sequencing data from two species with different life history traits, American lobster (Homarus americanus) and Arctic Char (Salvelinus alpinus), we highlight how an unbalanced sex ratio in the samples and a few sex-linked markers may lead to false interpretations of population structure and thus to potentially erroneous management recommendations. Here, multivariate analyses revealed two genetic clusters separating samples by sex instead of by expected spatial variation: inshore and offshore locations in lobster, or east and west locations in Arctic Char. To further investigate this, we created several subsamples artificially varying the sex ratio in the inshore/offshore and east/west groups and then demonstrated that significant genetic differentiation could be observed despite panmixia in lobster, and that FST values were overestimated in Arctic Char. This pattern was due to 12 and 94 sex-linked markers driving differentiation for lobster and Arctic Char, respectively. Removing sex-linked markers led to nonsignificant genetic structure in lobster and a more accurate estimation of FST in Arctic Char. The locations of these markers and putative identities of genes containing or nearby the markers were determined using available transcriptomic and genomic data, and this provided new information related to sex determination in both species. Given that only 9.6% of all marine/diadromous population genomic studies to date have reported sex information, we urge researchers to collect and consider individual sex information. Sex information is therefore relevant for avoiding unexpected biases due to sex-linked markers as well as for improving our knowledge of sex determination systems in nonmodel species.

The Cuban gar (Atractosteus tristoechus) is an endemic lepisosteid living in Cuba. Among gars, this species is one of the most threatened and has the smallest natural distribution range. Lepisosteids are air-breathing fishes belonging to the Holostean, a basal non-teleost clade of actinopterygians. Recent studies have indicated that these fishes could be a ‘bridge between tetrapods and teleost biomedical models’. Herein, we sequenced and assembled the first complete mitochondrial genome of A. tristoechus. The total length of the mt genome is 16,290 bp, containing the typical 13 protein-coding genes, two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a 537 bp length control region.