Genomic resources for the Yellowfin tuna Thunnus albacares.

BACKGROUND: The Yellowfin tuna (Thunnus albacares) is a large tuna exploited by major fisheries in tropical and subtropical waters of all oceans except the Mediterranean Sea. Genomic studies of population structure, adaptive variation or of the genetic basis of phenotypic traits are needed to inform fisheries management but are currently limited by the lack of a reference genome for this species. Here we report a draft genome assembly and a linkage map for use in genomic studies of T. albacares. METHODS AND RESULTS: Illumina and PacBio SMRT sequencing were used in combination to generate a hybrid assembly that comprises 743,073,847 base pairs contained in 2,661 scaffolds. The assembly has a N50 of 351,587 and complete and partial BUSCO scores of 86.47% and 3.63%, respectively. Double-digest restriction associated DNA (ddRAD) was used to genotype the 2 parents and 164 of their F1 offspring resulting from a controlled breeding cross, retaining 19,469 biallelic single nucleotide polymorphism (SNP) loci. The SNP loci were used to construct a linkage map that features 24 linkage groups that represent the 24 chromosomes of yellowfin tuna. The male and female maps span 1,243.8 cM and 1,222.9 cM, respectively. The map was used to anchor the assembly in 24 super-scaffolds that contain 79% of the yellowfin tuna genome. Gene prediction identified 46,992 putative genes 20,203 of which could be annotated via gene ontology. CONCLUSIONS: The draft reference will be valuable to interpret studies of genome wide variation in T. albacares and other Scombroid species.

Nursery origin of yellowfin tuna in the western Atlantic Ocean: significance of Caribbean Sea and trans-Atlantic migrants.

Natural geochemical markers in the otolith of yellowfin tuna (Thunnus albacares) were used to establish nursery-specific signatures for investigating the origin of fish captured in the western Atlantic Ocean (WAO). Two classes of chemical markers (trace elements, stable isotopes) were used to first establish nursery-specific signatures of age-0 yellowfin tuna from four primary production zones in the Atlantic Ocean: Gulf of Mexico, Caribbean Sea, Cape Verde, and Gulf of Guinea. Next, mixture and individual assignment methods were applied to predict the origin of sub-adult and adult yellowfin tuna from two regions in the WAO (Gulf of Mexico, Mid Atlantic Bight) by relating otolith core signatures (corresponding to age-0 period) to baseline signatures of age-0 fish from each nursery. Significant numbers of migrants from Caribbean Sea and eastern Atlantic Ocean (EAO) production zones (Gulf of Guinea, Cape Verde) were detected in the WAO, suggesting that fisheries in this region were subsidized by outside spawning/nursery areas. Contributions from local production (Gulf of Mexico) were also evident in samples from both WAO fisheries, but highly variable from year to year. High levels of mixing by yellowfin tuna from the different production zones and pronounced interannual trends in nursery-specific contribution rates in the WAO emphasize the complex and dynamic nature of this species' stock structure and population connectivity. Given that geographic shifts in distribution across national or political boundaries leads to governance and management challenges, this study highlights the need for temporally resolved estimates of nursery origin to refine assessment models and promote the sustainable harvest of this species.

Development and Evaluation of High-Density SNP Arrays for the Eastern Oyster Crassostrea virginica.

The eastern oyster Crassostrea virginica is a major aquaculture species for the USA. The sustainable development of eastern oyster aquaculture depends upon the continued improvement of cultured stocks through advanced breeding technologies. The Eastern Oyster Breeding Consortium (EOBC) was formed to advance the genetics and breeding of the eastern oyster. To facilitate efficient genotyping needed for genomic studies and selection, the consortium developed two single-nucleotide polymorphism (SNP) arrays for the eastern oyster: one screening array with 566K SNPs and one breeders' array with 66K SNPs. The 566K screening array was developed based on whole-genome resequencing data from 292 oysters from Atlantic and Gulf of Mexico populations; it contains 566,262 SNPs including 47K from protein-coding genes with a marker conversion rate of 48.34%. The 66K array was developed using best-performing SNPs from the screening array, which contained 65,893 oyster SNPs including 22,984 genic markers with a calling rate of 99.34%, a concordance rate of 99.81%, and a much-improved marker conversion rate of 92.04%. Null alleles attributable to large indels were found in 13.1% of the SNPs, suggesting that copy number variation is pervasive. Both arrays provided easy identification and separation of selected stocks from wild progenitor populations. The arrays contain 31 mitochondrial SNPs that allowed unambiguous identification of Gulf mitochondrial genotypes in some Atlantic populations. The arrays also contain 756 probes from 13 oyster and human pathogens for possible detection. Our results show that marker conversion rate is low in high polymorphism species and that the two-step process of array development can greatly improve array performance. The two arrays will advance genomic research and accelerate genetic improvement of the eastern oyster by delineating genetic architecture of production traits and enabling genomic selection. The arrays also may be used to monitor pedigree and inbreeding, identify selected stocks and their introgression into wild populations, and assess the success of oyster restoration.

Population structure of blackfin tuna (Thunnus atlanticus) in the western Atlantic Ocean inferred from microsatellite loci.

The blackfin tuna, Thunnus atlanticus, is a small tropical tuna exploited by recreational and commercial fisheries in various parts of its range. Information on stock structure is needed to develop management plans for this species but is currently lacking. In this work, 470 blackfin tuna from nine geographic populations were assayed at 13 homologous microsatellite markers to provide a first assessment of stock structure across the species range. The overall divergence among locality samples was very low (overall FST = 0.0004) indicating high connectivity of blackfin tuna across their range. No clear grouping of localities in differentiated units was inferred but structuring followed a weak isolation by distance pattern (r = 0.16, P = 0.032). Pairwise exact tests and spatial analysis of molecular variance suggested divergence of the sample collected offshore Baía Formosa (Brazil) possibly reflecting reproductive isolation of Brazilian populations from those in the Caribbean region and further north. Further study of the status of Brazilian populations and the transition between this region and the Caribbean is warranted. Cryptic subdivision within the Northern Hemisphere part of the range is possible and should be evaluated using increased marker density and a more comprehensive geographic coverage.

Sperm Repository for a Breeding Program of the Eastern Oyster Crassostrea virginica: Sample Collection, Processing, Cryopreservation, and Data Management Plan.

The Eastern oyster Crassostrea virginica (Family Ostreidae) is one of the most important fishery and aquaculture species in the U.S. and is a keystone species for coastal reefs. A breeding program was initiated in 2019 to support the fast-growing aquaculture industry culturing this species in the Gulf of Mexico. Oysters from 17 wild populations in embayment along the U.S. Gulf of Mexico coast from southwest Florida to the Matagorda Bay, Texas were used as broodstock for the program to maximize genetic diversity in the base population. A sperm repository of the broodstock was established to support the breeding project. The goal of this study was to demonstrate the sperm sample collection, processing, cryopreservation, and the data management plan involved in the establishment of a sperm germplasm repository of base populations. The supporting objectives were to: (1) develop a data management plan for the sperm repository; (2) streamline the procedure for sample collection, processing, and cryopreservation; (3) incorporate sperm quality analysis into the procedure, and (4) archive the cryopreserved samples as a repository for future use in the breeding program. This sperm repository included a total of 102 male oysters from the 17 collection sites (six oysters per site). A data management plan was developed with six categories, including sample collection, phenotype, fresh sperm, genotype, cryopreservation, and post-thaw sperm, as guide for data collection. Sperm collection was accomplished by strip spawn, and fresh sperm production, motility, and fertility were recorded for quality analysis. Cryopreserved sperm samples were sorted, labelled, archived, and stored in liquid nitrogen for future use. Post-thaw motility (1-30%) and plasm membrane integrity (15.34-70.36%) were recorded as post-thaw quality parameters. Overall, this study demonstrated a streamlined procedure of oyster sperm collection, processing, and cryopreservation for establishing a sperm repository that can serve as a template for construction of oyster germplasm repositories for breeding programs.

The rise and fall of the ancient northern pike master sex-determining gene.

The understanding of the evolution of variable sex determination mechanisms across taxa requires comparative studies among closely related species. Following the fate of a known master sex-determining gene, we traced the evolution of sex determination in an entire teleost order (Esociformes). We discovered that the northern pike (Esox lucius) master sex-determining gene originated from a 65 to 90 million-year-old gene duplication event and that it remained sex linked on undifferentiated sex chromosomes for at least 56 million years in multiple species. We identified several independent species- or population-specific sex determination transitions, including a recent loss of a Y chromosome. These findings highlight the diversity of evolutionary fates of master sex-determining genes and the importance of population demographic history in sex determination studies. We hypothesize that occasional sex reversals and genetic bottlenecks provide a non-adaptive explanation for sex determination transitions.

Community composition and antibiotic resistance of bacteria in bottlenose dolphins Tursiops truncatus - Potential impact of 2010 BP Oil Spill.

Aquatic contamination, oil spills in particular, could lead to the accumulation of antibiotic resistance by promoting selection for and/or transfer of resistance genes. However, there have been few studies on antibiotic resistance in marine mammals in relation to environmental disturbances, specifically oil contaminations. Here we initiated a study on antibiotic resistance bacteria in bottlenose dolphins Tursiops truncatus in relation to oil contamination following the 2010 BP Oil Spill in the northern Gulf of Mexico. Bacterial communities and antibiotic resistance prevalence one year after the 2010 BP Oil Spill were compared between Barataria Bay (BB) and Sarasota Bay (SB) by applying the rarefaction curve method, and (generalized) linear mixed models. The results showed that the most common bacteria included Vibrio, Shewanella, Bacillus and Pseudomonas. The prevalence of antibiotic resistance was high in the bacterial isolates at both bays. Though bacterial diversity did not differ significantly among water or dolphin samples, and antibiotic resistance did not differ significantly among water samples between the two bays, antibiotic resistance and multi-drug resistance in dolphin samples was significantly higher in the BB than in the SB, mainly attributed to the resistance to E, CF, FEP and SXT. We also found sulfamethoxazole-trimethoprim-resistant Stenotrophomonas maltophilia the first time in the natural aquatic environment. The higher antibiotic resistance in the dolphins in BB is likely attributed to 2010 BP Oil Spill as we expected SB, a more urbanized bay area, would have had higher antibiotic resistance based on the previous studies. The antibiotic resistance data gathered in this research will fill in the important data gaps and contributes to the broader spatial-scale emerging studies on antibiotic resistance in aquatic environments.

Development and characterization of genomic resources for a non-model marine teleost, the red snapper (Lutjanus campechanus, Lutjanidae): Construction of a high-density linkage map, anchoring of genome contigs and comparative genomic analysis.

The red snapper Lutjanus campechanus is an exploited reef fish of major economic importance in the Gulf of Mexico region. Studies of genome wide genetic variation are needed to understand the structure of wild populations and develop breeding programs for aquaculture but interpretation of these genome scans is limited by the absence of reference genome. In this work, the first draft of a reference genome was developed and characterized for the red snapper. P-454 and Illumina sequencing were conducted to produce paired-end reads that were assembled into reference contigs and scaffolds. The current assembly spans over 770 Mb, representing an estimated 69% of the red snapper genome in 67,254 scaffolds (N50 = 16,803 bp). The genome contigs were applied to map double digest Restriction-Site Associated DNA Tags and characterize Single Nucleotide Polymorphisms (SNPs) in five outbred full-sib families. The identified SNPs and 97 microsatellite loci were used to generate a high-density linkage map that includes 7,420 markers distributed across 24 linkage groups and spans 1,346.64 cM with an average inter-marker distance of 0.18 cM. Sex-specific maps revealed a 1.10:1 female to male map length ratio. A total of 4,422 genome contigs (10.5% of the assembly) were anchored to the map and used in a comparative genomic analysis of the red snapper and two model teleosts. Red snapper showed a high degree of chromosome level syntenic conservation with both medaka and spotted green puffer and a near one to one correspondence between the 24 red snapper linkage groups and corresponding medaka chromosomes was observed. This work established the first draft of a reference genome for a lutjanid fish. The obtained genomic resources will serve as a framework for the interpretation of genome scans during studies of wild populations and captive breeding programs.

Temperature effects and genotype-temperature interactions on sex determination in the European sea bass (Dicentrarchus labrax L.).

The effect of temperature on sex-ratios in 27 families of sea bass reared in the same tank from the fertilization stage onward was investigated. An excess of males (68%) was found in the groups that were reared at high temperature (mean +/- standard deviation: 20+/-1 degrees C) until they reached the mean size of 8.1 cm (Standard Length, 149 days post-fertilization [p.f.]). Masculinization was higher (89% of males) in the groups maintained at low temperature (13 degrees C), from fertilization to a mean length of 6.5 cm (346 days p.f.). Shifts from high to low temperature at 8.1cm and from low to high temperature at 6.5 cm had no consequence on the sex-ratio. The percentage of males showing intratesticular oocytes was higher at low temperature (63%) than at high temperature (36%), suggesting that these males may be sensitive fish that have been masculinized by environmental factors. Fish sampled in the groups reared at high (2,200 fish) and low (500 fish) temperature were genotyped on three microsatellite loci. This allowed them to be assigned to the breeders used in the crossing design, thus permitting an analysis of parental influence on sex-ratios. In groups reared at high temperature, both parents had a significant additive effect on the percentage of females, and the interaction between sire and dam was not significant. Genotype temperature interactions were also detected and their existence suggests the interesting possibility of selecting nonsensitive genotypes in breeding programs.