Hot alkaline lysis gDNA extraction from formalin-fixed archival tissues.

Formalin fixation of natural history specimens and histopathological material has historically been viewed as an impediment to successful genomic analysis. However, the development of extraction methods specifically tailored to contend with heavily crosslinked archival tissues, re-contextualises millions of previously overlooked specimens as viable molecular assets. Here, we present an easy-to-follow protocol for screening archival wet specimens for molecular viability and subsequent genomic DNA extraction suitable for sequencing. The protocol begins with non-destructive assessment of specimen degradation and preservation media conditions to allow both museum curators and researchers to select specimens most likely to yield an acceptable proportion (20-60%) of mappable endogenous DNA during short-read DNA sequencing. The extraction protocol uses hot alkaline lysis in buffer (0.1M NaOH, 1% SDS, pH 13) to simultaneously lyse and de-crosslink the tissue. To maximise DNA recovery, phenol:chloroform extraction is coupled with a small-fragment optimised SPRI bead clean up. Applied to well-preserved archival tissues, the protocol can yield 1-2 µg DNA per 50 mg of tissue with mean fragment sizes typically ranging from 50-150 bp, which is suitable to recover genomic DNA sufficient to reconstruct complete mitochondrial genomes and achieve up to 25X nuclear genome coverage. We provide guidance for read mapping to a reference genome and discuss the limitations of relying on small fragments for SNP genotyping and de novo genome assembly. This protocol opens the door to broader-scale genetic and phylogenetic analysis of historical specimens, contributing to a deeper understanding of evolutionary trends and adaptation in response to changing environments.

Low levels of sibship encourage use of larvae in western Atlantic bluefin tuna abundance estimation by close-kin mark-recapture.

Globally, tunas are among the most valuable fish stocks, but are also inherently difficult to monitor and assess. Samples of larvae of Western Atlantic bluefin tuna Thunnus thynnus (Linnaeus, 1758) from standardized annual surveys in the northern Gulf of Mexico provide a potential source of "offspring" for close-kin mark-recapture (CKMR) estimates of abundance. However, the spatial patchiness and highly skewed numbers of larvae per tow suggest sampled larvae may come from a small number of parents, compromising the precision of CKMR. We used high throughput genomic profiling to study sibship within and among larval tows from the 2016 standardized Gulf-wide survey compared to targeted sampling carried out in 2017. Full- and half-siblings were found within both years, with 12% of 156 samples in 2016 and 56% of 317 samples in 2017 having at least one sibling. There were also two pairs of cross cohort half-siblings. Targeted sampling increased the number of larvae collected per sampling event but resulted in a higher proportion of siblings. The combined effective sample size across both years was about 75% of the nominal size, indicating that Gulf of Mexico larval collections could be a suitable source of juveniles for CKMR in Western Atlantic bluefin tuna.

Lack of genetic differentiation in yellowfin tuna has conservation implications in the Eastern Pacific Ocean.

Yellowfin tuna, Thunnus albacares, is an important global fishery and of particular importance in the Eastern Pacific Ocean (EPO). According to the 2019 Inter-American Tropical Tuna Commission (IATTC) assessment, yellowfin tuna within the EPO is a single stock, and is being managed as one stock. However, previous studies indicate site fidelity, or limited home ranges, of yellowfin tuna which suggests the potential for multiple yellowfin tuna stocks within the EPO, which was supported by a population genetic study using microsatellites. If numerous stocks are present, management at the wrong spatial scales could cause the loss of minor yellowfin tuna populations in the EPO. In this study we used double digestion RADseq to assess the genetic structure of yellowfin tuna in the EPO. A total of 164 yellowfin tuna from Cabo San Lucas, México, and the Galápagos Islands and Santa Elena, Ecuador, were analysed using 18,011 single nucleotide polymorphisms. Limited genetic differentiation (FST = 0.00058-0.00328) observed among the sampling locations (México, Ecuador, Peru, and within Ecuador) is consistent with presence of a single yellowfin tuna population within the EPO. Our findings are consistent with the IATTC assessment and provide further evidence of the need for transboundary cooperation for the successful management of this important fishery throughout the EPO.

Accounting for kin sampling reveals genetic connectivity in Tasmanian and New Zealand school sharks, Galeorhinus galeus.

Fishing represents a major problem for conservation of chondrichthyans, with a quarter of all species being overexploited. School sharks, Galeorhinus galeus, are targeted by commercial fisheries in Australia and New Zealand. The Australian stock has been depleted to below 20% of its virgin biomass, and the species is recorded as Conservation Dependent within Australia. Individuals are known to move between both countries, but it is disputed whether the stocks are reproductively linked. Accurate and unbiased determination of stock and population connectivity is crucial to inform effective management. In this study, we assess the genetic composition and population connectivity between Australian and New Zealand school sharks using genome-wide SNPs, while accounting for non-random kin sampling. Between 2009 and 2013, 88 neonate and juvenile individuals from Tasmanian and New Zealand nurseries were collected and genotyped. Neutral loci were analyzed to detect fine-scale signals of reproductive connectivity. Seven full-sibling groups were identified and removed for unbiased analysis. Based on 6,587 neutral SNPs, pairwise genetic differentiation from Tasmanian and New Zealand neonates was non-significant (F ST = 0.0003, CI95 = [-0.0002, 0.0009], p = 0.1163; D est = 0.0006 ± 0.0002). This pattern was supported by clustering results. In conclusion, we show a significant effect of non-random sampling of kin and identify fine-scale reproductive connectivity between Australian and New Zealand school sharks. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://doi.org/10.5061/dryad.pd8612j.

Robust estimates of a high Ne/N ratio in a top marine predator, southern bluefin tuna.

Genetic studies of several marine species with high fecundity have produced "tiny" estimates (≤10-3) of the ratio of effective population size (Ne) to adult census size (N), suggesting that even very large populations might be at genetic risk. A recent study using close-kin mark-recapture methods estimated adult abundance at N ≈ 2 × 106 for southern bluefin tuna (SBT), a highly fecund top predator that supports a lucrative (~$1 billion/year) fishery. We used the same genetic and life history data (almost 13,000 fish collected over 5 years) to generate genetic and demographic estimates of Ne per generation and Nb (effective number of breeders) per year and the Ne/N ratio. Demographic estimates, which accounted for age-specific vital rates, skip breeding, variation in fecundity at age, and persistent individual differences in reproductive success, suggest that Ne/N is >0.1 and perhaps about 0.5. The genetic estimates supported this conclusion. Simulations using true Ne = 5 × 105 (Ne/N = 0.25) produced results statistically consistent with the empirical genetic estimates, whereas simulations using Ne = 2 × 104 (Ne/N = 0.01) did not. Our results show that robust estimates of Ne and Ne/N can be obtained for large populations, provided sufficiently large numbers of individuals and genetic markers are used and temporal replication (here, 5 years of adult and juvenile samples) is sufficient to provide a distribution of estimates. The high estimated Ne/N ratio in SBT is encouraging and suggests that the species will not be compromised by a lack of genetic diversity in responding to environmental change and harvest.

Reliably discriminating stock structure with genetic markers: Mixture models with robust and fast computation.

Delineating naturally occurring and self-sustaining subpopulations (stocks) of a species is an important task, especially for species harvested from the wild. Despite its central importance to natural resource management, analytical methods used to delineate stocks are often, and increasingly, borrowed from superficially similar analytical tasks in human genetics even though models specifically for stock identification have been previously developed. Unfortunately, the analytical tasks in resource management and human genetics are not identical-questions about humans are typically aimed at inferring ancestry (often referred to as "admixture") rather than breeding stocks. In this article, we argue, and show through simulation experiments and an analysis of yellowfin tuna data, that ancestral analysis methods are not always appropriate for stock delineation. In this work, we advocate a variant of a previously introduced and simpler model that identifies stocks directly. We also highlight that the computational aspects of the analysis, irrespective of the model, are difficult. We introduce some alternative computational methods and quantitatively compare these methods to each other and to established methods. We also present a method for quantifying uncertainty in model parameters and in assignment probabilities. In doing so, we demonstrate that point estimates can be misleading. One of the computational strategies presented here, based on an expectation-maximization algorithm with judiciously chosen starting values, is robust and has a modest computational cost.

Inferring contemporary and historical genetic connectivity from juveniles.

Measuring population connectivity is a critical task in conservation biology. While genetic markers can provide reliable long-term historical estimates of population connectivity, scientists are still limited in their ability to determine contemporary patterns of gene flow, the most practical time frame for management. Here, we tackled this issue by developing a new approach that only requires juvenile sampling at a single time period. To demonstrate the usefulness of our method, we used the Speartooth shark (Glyphis glyphis), a critically endangered species of river shark found only in tropical northern Australia and southern Papua New Guinea. Contemporary adult and juvenile shark movements, estimated with the spatial distribution of kin pairs across and within three river systems, was contrasted with historical long-term connectivity patterns, estimated from mitogenomes and genome-wide SNP data. We found strong support for river fidelity in juveniles with the within-cohort relationship analysis. Male breeding movements were highlighted with the cross-cohort relationship analysis, and female reproductive philopatry to the river systems was revealed by the mitogenomic analysis. We show that accounting for juvenile river fidelity and female philopatry is important in population structure analysis and that targeted sampling in nurseries and juvenile aggregations should be included in the genomic toolbox of threatened species management.

Absolute abundance of southern bluefin tuna estimated by close-kin mark-recapture.

Southern bluefin tuna is a highly valuable, severely depleted species, whose abundance and productivity have been difficult to assess with conventional fishery data. Here we use large-scale genotyping to look for parent-offspring pairs among 14,000 tissue samples of juvenile and adult tuna collected from the fisheries, finding 45 pairs in total. Using a modified mark-recapture framework where 'recaptures' are kin rather than individuals, we can estimate adult abundance and other demographic parameters such as survival, without needing to use contentious fishery catch or effort data. Our abundance estimates are substantially higher and more precise than previously thought, indicating a somewhat less-depleted and more productive stock. More broadly, this technique of 'close-kin mark-recapture' has widespread utility in fisheries and wildlife conservation. It estimates a key parameter for management-the absolute abundance of adults-while avoiding the expense of independent surveys or tag-release programmes, and the interpretational problems of fishery catch rates.

Suppression and restoration of primordial germ cell marker gene expression in channel catfish, Ictalurus punctatus, using knockdown constructs regulated by copper transport protein gene promoters: Potential for reversible transgenic sterilization.

Complementary DNA overexpression and short hairpin RNA interference approaches were evaluated for decreasing expression of primordial germ cell (PGC) marker genes and thereby sterilizing channel catfish, Ictalurus punctatus, by delivering knockdown constructs driven by a constitutive promoter from yeast and a copper transport protein gene into fish embryos by electroporation. Two PGC marker genes, nanos and dead end, were the target knockdown genes, and their expressions, along with that of an off-target gene, vasa, were evaluated temporally using real-time polymerase chain reaction. Copper sulfate was evaluated as a repressor compound. Some of the constructs knocked down PGC marker gene expression, and some of the constructs were partially repressed by application of 0.1-ppm copper sulfate. When the rate of sexual maturity was compared for three-year-old broodfish that had been exposed to the sterilizing constructs during embryologic development and controls that had not been exposed, several treatments had reduced sexual maturity for the exposed fish. Of two promoter systems evaluated, the one which had been designed to be less sensitive to copper generally was more effective at achieving sterilization and more responsive to repression. Knockdown constructs based on 3' nanos short hairpin RNA interference appeared to result in the best repression and restoration of normal sexual maturity. We conclude that these copper-based systems exhibited good potential for repressible transgenic sterilization. Optimization of this system could allow environmentally safe application of transgenic technology and might be applicable to other applications for aquatic organisms.

Whole mitogenome of the Endangered dwarf sawfish Pristis clavata (Rajiformes: Pristidae).

In this study, we describe the first complete mitochondrial sequence for the Endangered dwarf sawfish Pristis clavata. The base composition of the 16,804 bp long mitogenome is 31.9% A, 26.5% C, 13.3% G and 28.3% T and the gene arrangement and transcriptional direction are the same as those found in most vertebrates. All protein-coding genes start with ATG except the COI gene, which starts with GTG. Stop codons include incomplete T, AGG and TAA; however, TAG is not found in the mitogenome of this euryhaline elasmobranch species.

Complete mitogenomic sequence of the Critically Endangered Northern River Shark Glyphis garricki (Carcharhiniformes: Carcharhinidae).

In this study we describe the first complete mitochondrial sequence for the Critically Endangered Northern River shark Glyphis garricki. The complete mitochondrial sequence is 16,702 bp in length, contains 37 genes and one control region with the typical gene order and transcriptional direction of vertebrate mitogenomes. The overall base composition is 31.5% A, 26.3% C, 12.9% G and 29.3% T. The length of 22 tRNA genes ranged from 68 (tRNA-Ser2 and tRNA-Cys) to 75 (tRNA-Leu1) bp. The control region of G. garricki was 1067 bp in length with high A+T (67.9%) and poor G (12.6%) content. The mitogenomic characters (base composition, codon usage and gene length) of G. garricki were very similar to Glyphis glyphis.

Mitogenomics of the Speartooth Shark challenges ten years of control region sequencing.

BACKGROUND: Mitochondrial DNA markers have long been used to identify population boundaries and are now a standard tool in conservation biology. In elasmobranchs, evolutionary rates of mitochondrial genes are low and variation between distinct populations can be hard to detect with commonly used control region sequencing or other single gene approaches. In this study we sequenced the whole mitogenome of 93 Critically Endangered Speartooth Shark Glyphis glyphis from the last three river drainages they inhabit in northern Australia. RESULTS: Genetic diversity was extremely low (π =0.00019) but sufficient to demonstrate the existence of barriers to gene flow among river drainages (AMOVA Φ ST =0.28283, P <0.00001). Surprisingly, the comparison with single gene sub-datasets revealed that ND5 and 12S were the only ones carrying enough information to detect similar levels of genetic structure. The control region exhibited only one mutation, which was not sufficient to detect any structure among river drainages. CONCLUSIONS: This study strongly supports the use of single river drainages as discrete management units for the conservation of G. glyphis. Furthermore when genetic diversity is low, as is often the case in elasmobranchs, our results demonstrate a clear advantage of using the whole mitogenome to inform population structure compared to single gene approaches. More specifically, this study questions the extensive use of the control region as the preferential marker for elasmobranch population genetic studies and whole mitogenome sequencing will probably uncover a large amount of cryptic population structure in future studies.

Complete mitochondrial genome of the Critically Endangered speartooth shark Glyphis glyphis (Carcharhiniformes: Carcharhinidae).

In this study we present the first complete mitogenome for the speartooth shark Glyphis glyphis, a rare euryhaline elasmobranch from northern Australia and Papua New Guinea. The mitogenome is 16,702 bp in length and the overall base composition is 31.5% A; 26.0% C; 13.0% G and 29.5% T. It includes 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, 13 protein-coding genes and a putative 1066 bp long control region. The COI gene is initiated by GTG codon whereas the remaining protein-coding genes started with the ATG codon. This study will help elucidate the taxonomy of this poorly known group of sharks.

A PL10 vasa-like gene in the kuruma shrimp, Marsupenaeus japonicus, expressed during development and in adult gonad.

A PL10 vasa-like gene was isolated from the Kuruma shrimp Marsupenaeus japonicus and therefore called Mjpl10. It is differentially expressed during embryonic, larval, and postlarval development, and in female and male gonads. Using absolute real-time reverse transcriptase-polymerase chain reaction (RT-PCR), we demonstrate that Mjpl10 transcripts are present in the two-cell embryo, suggesting it is maternally expressed, and continually at low levels throughout embryogenesis. Mjpl10 expression increases significantly in the first 25 h after hatching (nauplii IV) and then decreases in a linear fashion by 316-fold over the next 52-day period. Its continued expression throughout embryonic and larval development is compatible with a conserved role in early germ cell specification. Transcript levels of Mjpl10 are also detected in the ovary and testes of mature adults.