Allergen Diversity and Abundance in Different Tissues of the Redclaw Crayfish (Cherax quadricarinatus)

Shellfish allergy affects ~2.5% of the global population and is a type I immune response resulting from exposure to crustacean and/or molluscan proteins. The Australian Redclaw crayfish (Cherax quadricarinatus) is a freshwater species endemic to and farmed in northern Australia and is becoming an aquaculture species of interest globally. Despite being consumed as food, allergenic proteins from redclaw have not been identified or characterised. In addition, as different body parts are often consumed, it is conceivable that redclaw tissues vary in allergenicity depending on tissue type and function. To better understand food-derived allergenicity, this study characterised allergenic proteins in various redclaw body tissues (the tail, claw, and cephalothorax) and how the stability of allergenic proteins was affected through cooking (raw vs. cooked tissues). The potential of redclaw allergens to cross-react and cause IgE-binding in patients allergic to other shellfish (i.e., shrimp) was also investigated. Raw and cooked extracts were prepared from each body part. SDS-PAGE followed by immunoblotting was performed to determine allergen-specific antibody reactivity to sarcoplasmic calcium-binding protein and hemocyanin, as well as to identify redclaw proteins binding to IgE antibodies from individual and pooled sera of shrimp-allergic patients. Liquid chromatography-mass spectrometry (LC/MS) was utilised to identify proteins and to determine the proportion within extracts. Known crustacean allergens were found in all tissues, with a variation in tissue distribution (e.g., higher levels of hemocyanin in the claw and cephalothorax than in the tail). The proportion of some allergens as a percentage of remaining heat-stable proteins increased in cooked tissues. Previously described heat-stable allergens (i.e., hemocyanin and sarcoplasmic calcium-binding protein) were found to be partially heat-labile. Immunoblotting indicated that shrimp-allergic patients cross-react to redclaw allergens. IgE-binding bands, analysed by LC/MS, identified up to 11 known shellfish allergens. The findings of this study provide fundamental knowledge into the diagnostic and therapeutic field of shellfish allergy.

Optimization of a non-activating medium for short-term chilled storage of barramundi (Lates calcarifer) testicular spermatozoa

Reliable short-term chilled sperm storage is a critical prerequisite to using advanced reproductive techniques for captive breeding of barramundi (Asian sea bass; Lates calcarifer). Marine Ringer's solution (MRS) is a common non-activating medium (NAM) and has previously been used to store sperm from wild-caught barramundi. However, MRS-stored spermatozoa from captive-bred barramundi were observed to lyse within 30 min incubation. Therefore, this study aimed to optimize the composition of NAM for short-term chilled storage by characterizing and mimicking the biochemical profile of seminal and blood plasma of captive-bred barramundi. To further understand the effect of each component, osmolality was first examined to determine its effect on sperm viability. Thereafter, the effects of NaHCO3, pH, and Na+ and K+ concentrations on sperm motility were investigated. Optimization of the NAM formula was achieved through iterative adaptions. The increase in NAM osmolality from 260 to 400 mOsm/kg led to a significant improvement in sperm viability. Moreover, using HEPES instead of NaHCO3 as buffering agent significantly enhanced sperm motility and velocity. As a result, sperm samples diluted with optimized NAM (185 mM NaCl, 5.1 mM KCl, 1.6 mM CaCl2·2H2O, 1.1 mM MgSO4·7H2O, 10.0 mM HEPES, 5.6 mM D+ glucose, 400 mOsm/kg, pH 7.4) and stored at 4 °C showed no significant loss in total motility for up to 48 h and retained progressive motility for up to 72 h. The optimized NAM developed in this study significantly extended the functional lifespan of spermatozoa during chilled storage, permitting the ongoing development of advanced reproductive technologies for barramundi.

Genetic parameters of color phenotypes of black tiger shrimp (Penaeus monodon)

Black tiger shrimp (Penaeus monodon) is the second most important aquaculture species of shrimp in the world. In addition to growth traits, uncooked and cooked body color of shrimp are traits of significance for profitability and consumer acceptance. This study investigated for the first time, the phenotypic and genetic variances and relationships for body weight and body color traits, obtained from image analyses of 838 shrimp, representing the progeny from 55 sires and 52 dams. The color of uncooked shrimp was subjectively scored on a scale from 1 to 4, with “1” being the lightest/pale color and “4” being the darkest color. For cooked shrimp color, shrimp were graded firstly by subjective scoring using a commercial grading score card, where the score ranged from 1 to 12 representing light to deep coloration which was subsequently found to not be sufficiently reliable with poor repeatability of measurement (r = 0.68–0.78) Therefore, all images of cooked color were regraded on a three-point scale from brightest and lightest colored cooked shrimp, to darkest and most color-intense, with a high repeatability (r = 0.80–0.92). Objective color of both cooked and uncooked color was obtained by measurement of RGB intensities (values range from 0 to 255) for each pixel from each shrimp. Using the “convertColor” function in “R”, the RGB values were converted to L*a*b* (CIE Lab) systems of color properties. This system of color space was established in 1976, by the International Commission of Illumination (CIE) where “L*” represents the measure of degree of lightness, values range from 0 to 100, where 0 = pure black and 100 = pure white. The value “a*” represents red to green coloration, where a positive value represents the color progression towards red and a negative value towards green. The value “b*” represents blue to yellow coloration, where a positive value refers to more yellowish and negative towards the blue coloration. In total, eight color-related traits were investigated. An ordinal mixed (threshold) model was adopted for manually (subjectively) scored color phenotypes, whereas all other traits were analyzed by linear mixed models using ASReml software to derive variance components and estimated breeding values (EBVs). Moderate to low heritability estimates (0.05–0.35) were obtained for body color traits. For subjectively scored cooked and uncooked color, EBV-based selection would result in substantial genetic improvement in these traits. The genetic correlations among cooked, uncooked and body weight traits were high and ranged from −0.88 to 0.81. These suggest for the first time that 1) cooked color can be improved indirectly by genetic selection based on color of uncooked/live shrimp, and 2) intensity of coloration is positively correlated with body weight traits and hence selection for body weight will also improve color traits in this population.

Genome-wide association study of VHSV-resistance trait in Paralichthys olivaceus

In flounder aquaculture, selective breeding plays a vital role in the development of disease-resistant traits and animals with high growth rates. Moreover, superior animals are required to achieve high profits. Unlike growth-related traits, disease-resistant experiments need to be conducted in a controlled environment, as the improper measurement of traits often leads to low genetic correlation and incorrect estimation of breeding values. In this study, viral hemorrhagic septicemia virus (VHSV) resistance was studied using a genome-wide association study (GWAS), and the genetic parameters were estimated. Genotyping was performed using a high-quality 70 K single nucleotide polymorphism (SNP) Affymetrix® Axiom® myDesign™ Genotyping Array of olive flounder. A heritability of ∼0.18 for resistance to VHSV was estimated using genomic information of the fish. According to the GWAS, significant SNPs were detected in chromosomes 21, 24, and contig AGQT02032065.1. Three SNPs showed significance at the genome-wide level (p < 1 × 10^-6), while others showed significance above the suggestive cutoff (p < 1 × 10^-4). The 3% phenotypic variation was explained by the highest significant SNP, named AX-419319631. Of the important genes for disease resistance, SNPs were associated with plcg1, epha4, clstn2, pik3cb, hes6, meis3, prx6, cep164, siae, and kirrel3b. Most of the genes associated with these SNPs have been previously reported with respect to viral entry, propagation, and immune mechanisms. Therefore, our study provides helpful information regarding VHSV resistance in olive flounder, which can be used for breeding applications.

Genome assembly of the Australian black tiger shrimp (Penaeus monodon) reveals a novel fragmented IHHNV EVE sequence

Shrimp are a valuable aquaculture species globally; however, disease remains a major hindrance to shrimp aquaculture sustainability and growth. Mechanisms mediated by endogenous viral elements have been proposed as a means by which shrimp that encounter a new virus start to accommodate rather than succumb to infection over time. However, evidence on the nature of such endogenous viral elements and how they mediate viral accommodation is limited. More extensive genomic data on Penaeid shrimp from different geographical locations should assist in exposing the diversity of endogenous viral elements. In this context, reported here is a PacBio Sequel-based draft genome assembly of an Australian black tiger shrimp (Penaeus monodon) inbred for 1 generation. The 1.89 Gbp draft genome is comprised of 31,922 scaffolds (N50: 496,398 bp) covering 85.9% of the projected genome size. The genome repeat content (61.8% with 30% representing simple sequence repeats) is almost the highest identified for any species. The functional annotation identified 35,517 gene models, of which 25,809 were protein-coding and 17,158 were annotated using interproscan. Scaffold scanning for specific endogenous viral elements identified an element comprised of a 9,045-bp stretch of repeated, inverted, and jumbled genome fragments of infectious hypodermal and hematopoietic necrosis virus bounded by a repeated 591/590 bp host sequence. As only near complete linear ∼4 kb infectious hypodermal and hematopoietic necrosis virus genomes have been found integrated in the genome of P. monodon previously, its discovery has implications regarding the validity of PCR tests designed to specifically detect such linear endogenous viral element types. The existence of joined inverted infectious hypodermal and hematopoietic necrosis virus genome fragments also provides a means by which hairpin double-stranded RNA could be expressed and processed by the shrimp RNA interference machinery.

Epigenetics underpins phenotypic plasticity of protandrous sex change in fish

Phenotypic plasticity is an important driver of species resilience. Often mediated by epigenetic changes, phenotypic plasticity enables individual genotypes to express variable phenotypes in response to environmental change. Barramundi (Lates calcarifer) are a protandrous (male‐first) sequential hermaphrodite that exhibits plasticity in length‐at‐sex change between geographic regions. This plasticity is likely to be mediated by changes in DNA methylation (DNAm), a well‐studied epigenetic modification. To investigate the relationships between length, sex, and DNAm in a sequential hermaphrodite, here, we compare DNAm in four conserved vertebrate sex‐determining genes in male and female barramundi of differing lengths from three geographic regions of northern Australia. Barramundi first mature as male and later sex change to female upon the attainment of a larger body size; however, a general pattern of increasing female‐specific DNAm markers with increasing length was not observed. Significant differences in DNAm between males and females of similar lengths suggest that female‐specific DNAm arises rapidly during sex change, rather than gradually with fish growth. The findings also reveal that region‐specific differences in length‐at‐sex change are accompanied by differences in DNAm and are consistent with variability in remotely sensed sea temperature and salinity. Together, these findings provide the first in situ evidence for epigenetically and environmentally mediated sex change in a protandrous hermaphrodite and offer significant insight into the molecular and ecological processes governing the marked and unique plasticity of sex in fish.

Population Structure, Genetic Connectivity, and Signatures of Local Adaptation of the Giant Black Tiger Shrimp (Penaeus monodon) throughout the Indo-Pacific Region

The giant black tiger shrimp (Penaeus monodon) is native to the Indo-Pacific and is the second most farmed penaeid shrimp species globally. Understanding genetic structure, connectivity, and local adaptation among Indo-Pacific black tiger shrimp populations is important for informing sustainable fisheries management and aquaculture breeding programs. Population genetic and outlier detection analyses were undertaken using 10,593 genome-wide single nucleotide polymorphisms (SNPs) from 16 geographically disparate Indo-Pacific P. monodon populations. Levels of genetic diversity were highest for Southeast Asian populations and were lowest for Western Indian Ocean (WIO) populations. Both neutral (n = 9,930) and outlier (n = 663) loci datasets revealed a pattern of strong genetic structure of P. monodon corresponding with broad geographical regions and clear genetic breaks among samples within regions. Neutral loci revealed seven genetic clusters and the separation of Fiji and WIO clusters from all other clusters, whereas outlier loci revealed six genetic clusters and high genetic differentiation among populations. The neutral loci dataset estimated five migration events that indicated migration to Southeast Asia from the WIO, with partial connectivity to populations in both oceans. We also identified 26 putatively adaptive SNPs that exhibited significant Pearson correlation (P < 0.05) between minor allele frequency and maximum or minimum sea surface temperature. Matched transcriptome contig annotations suggest putatively adaptive SNPs involvement in cellular and metabolic processes, pigmentation, immune response, and currently unknown functions. This study provides novel genome-level insights that have direct implications for P. monodon aquaculture and fishery management practices.

Scale Drop Disease Virus (SDDV) and Lates calcarifer Herpes Virus (LCHV) Coinfection Downregulate Immune-Relevant Pathways and Cause Splenic and Kidney Necrosis in Barramundi Under Commercial Farming Conditions

Marine farming of barramundi (Lates calcarifer) in Southeast Asia is currently severely affected by viral diseases. To better understand the biological implications and gene expression response of barramundi in commercial farming conditions during a disease outbreak, the presence of pathogens, comparative RNAseq, and histopathology targeting multiple organs of clinically “sick” and “healthy” juveniles were investigated. Coinfection of scale drop disease virus (SDDV) and L. calcarifer herpes virus (LCHV) were detected in all sampled fish, with higher SDDV viral loads in sick than in healthy fish. Histopathology showed that livers in sick fish often had moderate to severe abnormal fat accumulation (hepatic lipidosis), whereas the predominant pathology in the kidneys shows moderate to severe inflammation and glomerular necrosis. The spleen was the most severely affected organ, with sick fish presenting severe multifocal and coalescing necrosis. Principal component analysis (PC1 and PC2) explained 70.3% of the observed variance and strongly associated the above histopathological findings with SDDV loads and with the sick phenotypes, supporting a primary diagnosis of the fish being impacted by scale drop disease (SDD). Extracted RNA from kidney and spleen of the sick fish were also severely degraded likely due to severe inflammation and tissue necrosis, indicating failure of these organs in advanced stages of SDD. RNAseq of sick vs. healthy barramundi identified 2,810 and 556 differentially expressed genes (DEGs) in the liver and muscle, respectively. Eleven significantly enriched pathways (e.g., phagosome, cytokine-cytokine-receptor interaction, ECM-receptor interaction, neuroactive ligand-receptor interaction, calcium signaling, MAPK, CAMs, etc.) and gene families (e.g., tool-like receptor, TNF, lectin, complement, interleukin, chemokine, MHC, B and T cells, CD molecules, etc.) relevant to homeostasis and innate and adaptive immunity were mostly downregulated in sick fish. These DEGs and pathways, also previously identified in L. calcarifer as general immune responses to other pathogens and environmental stressors, suggest a failure of the clinically sick fish to cope and overcome the systemic inflammatory responses and tissue degeneration caused by SDD.

Novel Allergen Discovery through Comprehensive De Novo Transcriptomic Analyses of Five Shrimp Species

Shellfish allergy affects 2% of the world’s population and persists for life in most patients. The diagnosis of shellfish allergy, in particular shrimp, is challenging due to the similarity of allergenic proteins from other invertebrates. Despite the clinical importance of immunological cross-reactivity among shellfish species and between allergenic invertebrates such as dust mites, the underlying molecular basis is not well understood. Here we mine the complete transcriptome of five frequently consumed shrimp species to identify and compare allergens with all known allergen sources. The transcriptomes were assembled de novo, using Trinity, from raw RNA-Seq data of the whiteleg shrimp (Litopenaeus vannamei), black tiger shrimp (Penaeus monodon), banana shrimp (Fenneropenaeus merguiensis), king shrimp (Melicertus latisulcatus), and endeavour shrimp (Metapenaeus endeavouri). BLAST searching using the two major allergen databases, WHO/IUIS Allergen Nomenclature and AllergenOnline, successfully identified all seven known crustacean allergens. The analyses revealed up to 39 unreported allergens in the different shrimp species, including heat shock protein (HSP), alpha-tubulin, chymotrypsin, cyclophilin, beta-enolase, aldolase A, and glyceraldehyde-3-phosphate dehydrogenase (G3PD). Multiple sequence alignment (Clustal Omega) demonstrated high homology with allergens from other invertebrates including mites and cockroaches. This first transcriptomic analyses of allergens in a major food source provides a valuable resource for investigating shellfish allergens, comparing invertebrate allergens and future development of improved diagnostics for food allergy.

Comparing Genomic Signatures of Selection Between the Abbassa Strain and Eight Wild Populations of Nile Tilapia (Oreochromis niloticus) in Egypt

Domestication to captive rearing conditions, along with targeted selective breeding have genetic consequences that vary from those in wild environments. Nile tilapia (Oreochromis niloticus) is one of the most translocated and farmed aquaculture species globally, farmed throughout Asia, North and South America, and its African native range. In Egypt, a breeding program established the Abbassa Strain of Nile tilapia (AS) in 2002 based on local broodstock sourced from the Nile River. The AS has been intensively selected for growth and has gone through genetic bottlenecks which have likely shifted levels and composition of genetic diversity within the strain. Consequently, there are questions on the possible genetic impact AS escapees may have on endemic populations of Nile tilapia. However, to date there have been no genetic studies comparing genetic changes in the domesticated AS to local wild populations. This study used 9,827 genome-wide SNPs to investigate population genetic structure and signatures of selection in the AS (generations 9–11) and eight wild Nile tilapia populations from Egypt. SNP analyses identified two major genetic clusters (captive and wild populations), with wild populations showing evidence of isolation-by-distance among the Nile Delta and upstream riverine populations. Between genetic clusters, approximately 6.9% of SNPs were identified as outliers with outliers identified on all 22 O. niloticus chromosomes. A lack of localized outlier clustering on the genome suggests that no genes of major effect were presently detected. The AS has retained high levels of genetic diversity (H_{o_All} = 0.21 ± 0.01; H_{e_All} = 0.23 ± 0.01) when compared to wild populations (H_{o_All} = 0.18 ± 0.01; H_{e_All} = 0.17 ± 0.01) after 11 years of domestication and selective breeding. Additionally, 565 SNPs were unique within the AS line. While these private SNPs may be due to domestication signals or founder effects, it is suspected that introgression with blue tilapia (Oreochromis aureus) has occurred. This study highlights the importance of understanding the effects of domestication in addition to wild population structure to inform future management and dissemination decisions. Furthermore, by conducting a baseline genetic study of wild populations prior to the dissemination of a domestic line, the effects of aquaculture on these populations can be monitored over time.

Fine-scale population structure and evidence for local adaptation in Australian giant black tiger shrimp (Penaeus monodon) using SNP analysis

BACKGROUND

Restrictions to gene flow, genetic drift, and divergent selection associated with different environments are significant drivers of genetic differentiation. The black tiger shrimp (Penaeus monodon), is widely distributed throughout the Indian and Pacific Oceans including along the western, northern and eastern coastline of Australia, where it is an important aquaculture and fishery species. Understanding the genetic structure and the influence of environmental factors leading to adaptive differences among populations of this species is important for farm genetic improvement programs and sustainable fisheries management.

RESULTS

Based on 278 individuals obtained from seven geographically disparate Australian locations, 10,624 high-quality SNP loci were used to characterize genetic diversity, population structure, genetic connectivity, and adaptive divergence. Significant population structure and differentiation were revealed among wild populations (average F_ST = 0.001-0.107; p <  0.05). Eighty-nine putatively outlier SNPs were identified to be potentially associated with environmental variables by using both population differentiation (BayeScan and PCAdapt) and environmental association (redundancy analysis and latent factor mixed model) analysis methods. Clear population structure with similar spatial patterns were observed in both neutral and outlier markers with three genetically distinct groups identified (north Queensland, Northern Territory, and Western Australia). Redundancy, partial redundancy, and multiple regression on distance matrices analyses revealed that both geographical distance and environmental factors interact to generate the structure observed across Australian P. monodon populations.

CONCLUSION

This study provides new insights on genetic population structure of Australian P. monodon in the face of environmental changes, which can be used to advance sustainable fisheries management and aquaculture breeding programs.

Multi-species transcriptomics reveals evolutionary diversity in the mechanisms regulating shrimp tail muscle excitation-contraction coupling

The natural flight response in shrimp is powered by rapid contractions of the abdominal muscle fibres to propel themselves backwards away from perceived danger. This muscle contraction is dependent on repetitive depolarization of muscle plasma membrane, triggering tightly spaced cytoplasmic [Ca²⁺] transients and rapidly rising tetanic force responses. To achieve such high amplitude and high frequency of Ca²⁺ transients requires a high abundance of sarcoplasmic/endoplasmic reticulum Ca²⁺ ATPase (SERCA) to rapidly clear cytoplasmic Ca²⁺ between each transient and an efficient Ca²⁺ release system consisting of the Ryanodine Receptor (RyR), and voltage gated Ca²⁺ channels (Ca_Vs). With the aim to expand our knowledge of muscle gene function and identify orthologous genes regulating muscle excitation-contraction (EC) coupling, this study assembled nine Penaeid shrimp muscle transcriptomes. On average, the nine transcriptomes contained 27,000 contigs, with an annotation rate of 36% and a BUSCO completeness of 70%. Despite maintaining their function, the crustacean RyR and Ca_V proteins showed evidence of significant diversification from mammalian orthologs, while SERCA remained more conserved. Several key components of protein interaction were conserved, while others showed distinct crustacean specific evolutionary adaptations. Lastly, this study revealed approximately 1,000 orthologous genes involved in muscle specific processes present across all nine species.

Genomic Selection in Aquaculture: Application, Limitations and Opportunities With Special Reference to Marine Shrimp and Pearl Oysters

Within aquaculture industries, selection based on genomic information (genomic selection) has the profound potential to change genetic improvement programs and production systems. Genomic selection exploits the use of realized genomic relationships among individuals and information from genome-wide markers in close linkage disequilibrium with genes of biological and economic importance. We discuss the technical advances, practical requirements, and commercial applications that have made genomic selection feasible in a range of aquaculture industries, with a particular focus on molluscs (pearl oysters, Pinctada maxima) and marine shrimp (Litopenaeus vannamei and Penaeus monodon). The use of low-cost genome sequencing has enabled cost-effective genotyping on a large scale and is of particular value for species without a reference genome or access to commercial genotyping arrays. We highlight the pitfalls and offer the solutions to the genotyping by sequencing approach and the building of appropriate genetic resources to undertake genomic selection from first-hand experience. We describe the potential to capture large-scale commercial phenotypes based on image analysis and artificial intelligence through machine learning, as inputs for calculation of genomic breeding values. The application of genomic selection over traditional aquatic breeding programs offers significant advantages through being able to accurately predict complex polygenic traits including disease resistance; increasing rates of genetic gain; minimizing inbreeding; and negating potential limiting effects of genotype by environment interactions. Further practical advantages of genomic selection through the use of large-scale communal mating and rearing systems are highlighted, as well as presenting rate-limiting steps that impact on attaining maximum benefits from adopting genomic selection. Genomic selection is now at the tipping point where commercial applications can be readily adopted and offer significant short- and long-term solutions to sustainable and profitable aquaculture industries.

Genome-wide SNP analyses reveal high gene flow and signatures of local adaptation among the scalloped spiny lobster (Panulirus homarus) along the Omani coastline

BACKGROUND

The scalloped spiny lobster (Panulirus homarus) is a popular seafood commodity worldwide and an important export item from Oman. Annual catches in commercial fisheries are in serious decline, which has resulted in calls for the development of an integrated stock management approach. In Oman, the scalloped spiny lobster is currently treated as a single management unit (MU) or stock and there is an absence of information on the genetic population structure of the species that can inform management decisions, particularly at a fine-scale level. This work is the first to identify genome-wide single nucleotide polymorphisms (SNPs) for P. homarus using Diversity Arrays Technology sequencing (DArT-seq) and to elucidate any stock structure in the species.

RESULTS

After stringent filtering, 7988 high utility SNPs were discovered and used to assess the genetic diversity, connectivity and structure of P. homarus populations from Al Ashkharah, Masirah Island, Duqm, Ras Madrakah, Haitam, Ashuwaymiyah, Mirbat and Dhalkut landing sites. Pairwise FST estimates revealed low differentiation among populations (pairwise FST range = - 0.0008 - 0.0021). Analysis of genetic variation using putatively directional FST outliers (504 SNPs) revealed higher and significant pairwise differentiation (p < 0.01) for all locations, with Ashuwaymiyah being the most diverged population (Ashuwaymiyah pairwise FST range = 0.0288-0.0736). Analysis of population structure using Discriminant Analysis of Principal Components (DAPC) revealed a broad admixture among P. homarus, however, Ashuwaymiyah stock appeared to be potentially under local adaptive pressures. Fine scale analysis using Netview R provided further support for the general admixture of P. homarus.

CONCLUSIONS

Findings here suggested that stocks of P. homarus along the Omani coastline are admixed. Yet, fishery managers need to treat the lobster stock from Ashuwaymiyah with caution as it might be subject to local adaptive pressures. We emphasize further study with larger number of samples to confirm the genetic status of the Ashuwaymiyah stock. The approach utilised in this study has high transferability in conservation and management of other marine stocks with similar biological and ecological attributes.

Sex-specific dmrt1 and cyp19a1 methylation and alternative splicing in gonads of the protandrous hermaphrodite barramundi

Epigenetics is involved in sex differentiation of gonochoristic and hermaphroditic fish species, whereby two genes dmrt1 (pro-male) and cyp19a1 (pro-female) are known to play major roles. Barramundi, Lates calcarifer, is an important tropical aquaculture species that undergo natural and permanent male to female sex change, a process for which the exact underlying molecular mechanisms are still unknown. To elucidate whether DNA methylation is involved in sex control of barramundi, a next-generation bisulfite amplicon sequencing approach was used to target 146 CpG sites within proximal promoters and first exons of seven sex-related genes (dmrt1, cyp19a1, amh, foxl2, nr5a2, sox8 and sox9) of 24 testis and 18 ovaries of captive and wild adult barramundi. Moreover, comparative expression profiles of the key dmrt1 and cyp19a1 genes were further investigated using RT-qPCR and Sanger sequencing approaches, whereas expression levels of remaining targeted genes were based on available literature for the species. Results showed that cyp19a1 and amh were more methylated in males, whereas dmrt1 and nr5a2 were more methylated in females (P < 0.001), with no gender differences found for foxl2, sox8 or sox9 genes (P > 0.05). Sex-biased promoter DNA methylation was inversely related to gene expression only for dmrt1 and nr5a2, and directly related to amh expression, whereas no differences in cyp19a1 expression were found between testes and ovaries. Notably, unique sex-specific alternative splicing of dmrt1 and cyp19a1 were discovered, whereby males lacked the full-length aromatase coding cyp19a1 mRNA due to partial or total exon splicing, and females lacked the dmrt1 exon containing the DM-domain sequence. This study advances the current knowledge aiming to elucidate the genetic mechanisms within male and female gonads of this large protandrous hermaphrodite by providing the first evidence of epigenetics and alternative splicing simultaneously affecting key genes (cyp19a1 and dmrt1) central to sex differentiation pathways.

The State of "Omics" Research for Farmed Penaeids: Advances in Research and Impediments to Industry Utilization

Elucidating the underlying genetic drivers of production traits in agricultural and aquaculture species is critical to efforts to maximize farming efficiency. "Omics" based methods (i.e., transcriptomics, genomics, proteomics, and metabolomics) are increasingly being applied to gain unprecedented insight into the biology of many aquaculture species. While the culture of penaeid shrimp has increased markedly, the industry continues to be impeded in many regards by disease, reproductive dysfunction, and a poor understanding of production traits. Extensive effort has been, and continues to be, applied to develop critical genomic resources for many commercially important penaeids. However, the industry application of these genomic resources, and the translation of the knowledge derived from "omics" studies has not yet been completely realized. Integration between the multiple "omics" resources now available (i.e., genome assemblies, transcriptomes, linkage maps, optical maps, and proteomes) will prove critical to unlocking the full utility of these otherwise independently developed and isolated resources. Furthermore, emerging "omics" based techniques are now available to address longstanding issues with completing keystone genome assemblies (e.g., through long-read sequencing), and can provide cost-effective industrial scale genotyping tools (e.g., through low density SNP chips and genotype-by-sequencing) to undertake advanced selective breeding programs (i.e., genomic selection) and powerful genome-wide association studies. In particular, this review highlights the status, utility and suggested path forward for continued development, and improved use of "omics" resources in penaeid aquaculture.

Gene flow and genetic structure in Nile perch, Lates niloticus, from African freshwater rivers and lakes

BACKGROUND

Geological evolution of the African continent has been subject to complex processes including uplift, volcanism, desert formation and tectonic rifting. This complex geology has created substantial biogeographical barriers, and coupled with anthropogenic introductions of freshwater fishes, has influenced the genetic diversity, connectivity and sub-structuring of the teleost fauna. Nile perch, Lates niloticus, is an iconic fish in Africa and is of high commercial importance, both in the species' native range and where it has been translocated. However, the species is in decline and there is a need to understand its population genetic structure to facilitate sustainable management of the fishery and aquaculture development.

METHODOLOGY

Nile perch tissue samples were acquired from two West and four East (Lakes; Albert, Kyoga, Victoria and Turkana) African locations. Nineteen polymorphic microsatellite loci were used to study the genetic variation among populations across regions (West and East Africa), as well as between native and introduced environments within East Africa.

PRINCIPAL FINDINGS AND THEIR SIGNIFICANCE

Results revealed strong and significant genetic structuring among populations across the sampled distribution (divergence across regions, FCT = 0.26, P = 0.000). STRUCTURE analysis at a broad scale revealed K = 2 clusters, the West African individuals were assigned to one cluster, while all individuals from the East African region, regardless of whether native or introduced, were assigned to another cluster. The distinct genetic clusters identified in the current study between the West and East African Nile perch, appear to have been maintained by presence of biogeographic barriers and restricted gene flow between the two regions. Therefore, any translocations of Nile perch should be carefully considered across the regions of West and East Africa. Further analysis at a regional scale revealed further structuring of up to K = 3 genetic clusters in East African Nile perch. Significantly (P < 0.05) lower genetic diversity based on analysis of allelic richness (AR) was obtained for the two translocated populations of Lake Kyoga (AR = 3.61) and Lake Victoria (AR = 3.52), compared to Nile perch populations from their putative origins of Lakes Albert (AR = 4.12) and Turkana (AR = 4.43). The lower genetic diversity in the translocated populations may be an indication of previous bottlenecks and may also indicate a difficulty for these populations to persist and adapt to climatic changes and anthropogenic pressures that are currently present in the East African region.

Current status of parasitic ciliates Chilodonella spp. (Phyllopharyngea: Chilodonellidae) in freshwater fish aquaculture

Freshwater fish farming contributes to more than two-thirds of global aquaculture production. Parasitic ciliates are one of the largest causes of production loss in freshwater farmed fishes, with species from the genus Chilodonella being particularly problematic. While Chilodonella spp. include 'free-living' fauna, some species are involved in mortality events of fish, particularly in high-density aquaculture. Indeed, chilodonellosis causes major productivity losses in over 16 species of farmed freshwater fishes in more than 14 countries. Traditionally, Chilodonella species are identified based on morphological features; however, the genus comprises yet uncharacterized cryptic species, which indicates the necessity for molecular diagnostic methods. This review synthesizes current knowledge on the biology, ecology and geographic distribution of harmful Chilodonella spp. and examines pathological signs, diagnostic methods and treatments. Recent advances in molecular diagnostics and the ability to culture Chilodonella spp. in vitro will enable the development of preventative management practices and sustained freshwater fish aquaculture production.

Swept away: ocean currents and seascape features influence genetic structure across the 18,000 Km Indo-Pacific distribution of a marine invertebrate, the black-lip pearl oyster Pinctada margaritifera

BACKGROUND

Genetic structure in many widely-distributed broadcast spawning marine invertebrates remains poorly understood, posing substantial challenges for their fishery management, conservation and aquaculture. Under the Core-Periphery Hypothesis (CPH), genetic diversity is expected to be highest at the centre of a species' distribution, progressively decreasing with increased differentiation towards outer range limits, as populations become increasingly isolated, fragmented and locally adapted. The unique life history characteristics of many marine invertebrates such as high dispersal rates, stochastic survival and variable recruitment are also likely to influence how populations are organised. To examine the microevolutionary forces influencing population structure, connectivity and adaptive variation in a highly-dispersive bivalve, populations of the black-lip pearl oyster Pinctada margaritifera were examined across its ~18,000 km Indo-Pacific distribution.

RESULTS

Analyses utilising 9,624 genome-wide SNPs and 580 oysters, discovered differing patterns of significant and substantial broad-scale genetic structure between the Indian and Pacific Ocean basins. Indian Ocean populations were markedly divergent (F st = 0.2534-0.4177, p < 0.001), compared to Pacific Ocean oysters, where basin-wide gene flow was much higher (F st = 0.0007-0.1090, p < 0.001). Partitioning of genetic diversity (hierarchical AMOVA) attributed 18.1% of variance between ocean basins, whereas greater proportions were resolved within samples and populations (45.8% and 35.7% respectively). Visualisation of population structure at selectively neutral loci resolved three and five discrete genetic clusters for the Indian and Pacific Oceans respectively. Evaluation of genetic structure at adaptive loci for Pacific populations (89 SNPs under directional selection; F st = 0.1012-0.4371, FDR = 0.05), revealed five clusters identical to those detected at neutral SNPs, suggesting environmental heterogeneity within the Pacific. Patterns of structure and connectivity were supported by Mantel tests of isolation by distance (IBD) and independent hydrodynamic particle dispersal simulations.

CONCLUSIONS

It is evident that genetic structure and connectivity across the natural range of P. margaritifera is highly complex, and produced by the interaction of ocean currents, IBD and seascape features at a broad scale, together with habitat geomorphology and local adaptation at regional levels. Overall population organisation is far more elaborate than generalised CPH predictions, however valuable insights for regional fishery management, and a greater understanding of range-wide genetic structure in a highly-dispersive marine invertebrate have been gained.

Novel Polymorphic Microsatellite Markers for Panulirus ornatus and their Cross-species Primer Amplification in Panulirus homarus

Polymorphic microsatellite loci were isolated for Panulirus ornatus using 454 GS-FLX Titanium pyrosequencing. Fifteen markers containing perfect di-, tri-, tetra-, and penta-nucleotide motifs were consistently co-amplified in five multiplexes in a panel of 91 randomly selected samples. Observed number of alleles varied from 2 to 14 per locus. Observed and expected heterozygosity ranged from 0.090 to 0.79 and 0.08 to 0.87, respectively. Ten loci deviated from Hardy-Weinberg equilibrium after sequential Bonferroni correction. Genetic linkage disequilibrium analysis between all pairs of the loci showed significant departure from the null hypothesis between 11 loci. The microsatellite markers were also amplified successfully in related Panulirus homarus species with adequate level of polymorphism. The successful cross-species primer amplification of the 15 microsatellites indicates the potential of the developed markers to be transferred to other Panulirus species. The 15 novel microsatellite markers reported in this work add to the previously characterized markers by our group, exhibit adequate levels of polymorphism for wide range of future studies investigating population structure, genetic diversity, and evolutionary relationships among Panulirus species.

A Parallel Population Genomic and Hydrodynamic Approach to Fishery Management of Highly-Dispersive Marine Invertebrates: The Case of the Fijian Black-Lip Pearl Oyster Pinctada margaritifera

Fishery management and conservation of marine species increasingly relies on genetic data to delineate biologically relevant stock boundaries. Unfortunately for high gene flow species which may display low, but statistically significant population structure, there is no clear consensus on the level of differentiation required to resolve distinct stocks. The use of fine-scale neutral and adaptive variation, considered together with environmental data can offer additional insights to this problem. Genome-wide genetic data (4,123 SNPs), together with an independent hydrodynamic particle dispersal model were used to inform farm and fishery management in the Fijian black-lip pearl oyster Pinctada margaritifera, where comprehensive fishery management is lacking, and the sustainability of exploitation uncertain. Weak fine-scale patterns of population structure were detected, indicative of broad-scale panmixia among wild oysters, while a hatchery-sourced farmed population exhibited a higher degree of genetic divergence (Fst = 0.0850-0.102). This hatchery-produced population had also experienced a bottleneck (NeLD = 5.1; 95% C.I. = [5.1-5.3]); compared to infinite NeLD estimates for all wild oysters. Simulation of larval transport pathways confirmed the existence of broad-scale mixture by surface ocean currents, correlating well with fine-scale patterns of population structuring. Fst outlier tests failed to detect large numbers of loci supportive of selection, with 2-5 directional outlier SNPs identified (average Fst = 0.116). The lack of biologically significant population genetic structure, absence of evidence for local adaptation and larval dispersal simulation, all indicate the existence of a single genetic stock of P. margaritifera in the Fiji Islands. This approach using independent genomic and oceanographic tools has allowed fundamental insights into stock structure in this species, with transferability to other highly-dispersive marine taxa for their conservation and management.

Novel genomic microsatellite markers for genetic population and diversity studies of tropical scalloped spiny lobster (Panulirus homarus) and their potential application in related Panulirus species

Fourteen polymorphic microsatellites with perfect di-, tri-, and tetra-nucleotide repeats were identified for Panulirus homarus using Roche 454 whole-genome sequencing method. Microsatellites were efficiently co-amplified in four multiplexes and a singleplex, providing consistent and easily interpretable genotypes. The number of alleles per locus ranged from 2 to 11 with the observed and expected heterozygosity ranging between 0.000-0.532 and 0.031-0.836, respectively. A significant deviation from Hardy-Weinberg equilibrium was observed for majority of the loci, probably due to homozygote excess. Genetic linkage disequilibrium analysis between all the possible pairs of the loci showed significant departure from the null hypothesis in the loci pairs Pho-G11-Pho-G33 and Pho-G33-Pho-G57. High success in primer cross-species amplification of these microsatellite markers indicates their utility for genetic studies of different Panulirus species.

Chromosomal-Level Assembly of the Asian Seabass Genome Using Long Sequence Reads and Multi-layered Scaffolding

We report here the ~670 Mb genome assembly of the Asian seabass (Lates calcarifer), a tropical marine teleost. We used long-read sequencing augmented by transcriptomics, optical and genetic mapping along with shared synteny from closely related fish species to derive a chromosome-level assembly with a contig N50 size over 1 Mb and scaffold N50 size over 25 Mb that span ~90% of the genome. The population structure of L. calcarifer species complex was analyzed by re-sequencing 61 individuals representing various regions across the species' native range. SNP analyses identified high levels of genetic diversity and confirmed earlier indications of a population stratification comprising three clades with signs of admixture apparent in the South-East Asian population. The quality of the Asian seabass genome assembly far exceeds that of any other fish species, and will serve as a new standard for fish genomics.

Isolation, characterization, and multiplexing of novel microsatellite markers for the tropical scalloped spiny lobster (Panulirus homarus)

Of the various spiny lobster species in the tropical and subtropical Indo-West Pacific region, the tropical scalloped spiny lobster (Panulirus homarus) supports one of the most commercially valuable fishery resources in many coastal African and Asian countries. The last decade has witnessed a serious decline in the wild populations of this species. Knowledge of the genetic basis of spiny lobster population structure is a prerequisite to achieve a sustainable fisheries management for this species. Here, we describe 13 novel polymorphic microsatellite markers developed for P. homarus, using a cross-species primer design strategy based on P. ornatus Roche 454 shot-gun generated sequencing. Microsatellite polymorphisms were assessed in 96 unrelated P. homarus individuals of a natural population, with the number of alleles per locus varying from 2 to 14, the observed and expected heterozygosity from 0.00 to 0.78 and from 0.03 to 0.79, respectively, and with only four loci (Pho-G27, Pho-G32, Pho-G36, and Pho-G58) deviating from Hardy- Weinberg equilibrium. Genetic linkage disequilibrium analysis between all pairs of the loci showed significant departure from the null hypothesis between loci Pho-G22 - Pho-G30, and Pho-G30 - Pho-G35. The successful cross amplification of these microsatellites highlights the potential of the developed microsatellites for future population genetic research within the different Panulirus species.

Fishing for divergence in a sea of connectivity: The utility of ddRADseq genotyping in a marine invertebrate, the black-lip pearl oyster Pinctada margaritifera

Population genomic investigations on highly dispersive marine organisms typically require thousands of genome-wide SNP loci to resolve fine-scale population structure and detect signatures of selection. This information is important for species conservation efforts and stock management in both wild and captive populations, as well as genome mapping and genome wide association studies. Double digest Restriction site-Associated DNA Sequencing (ddRADseq) is a recent tool for delivering genome wide SNPs for non-model organisms. However, its application to marine invertebrate taxa has been limited, particularly given the complex and highly repetitive nature of many of these organisms' genomes. This study develops and evaluates an optimised ddRADseq technique together with associated analyses for generating genome-wide SNP data, and performs population genomic analyses to inform aquaculture and fishery management of a marine bivalve, the black-lip pearl oyster Pinctada margaritifera. A total of 5243 high-quality genome-wide SNP markers were detected, and used to assess population structure, genome diversity, detect Fst outliers and perform association testing in 156 individuals belonging to three wild and one hatchery produced populations from the Fiji Islands. Shallow but significant population structure was revealed among all wild populations (average pairwise Fst=0.046) when visualised with DAPC and an individual network analysis (NetView P), with clear evidence of a genetic bottleneck in the hatchery population (NeLD=6.1), compared to wild populations (NeLD>192.5). Fst outlier detection revealed 42-62 highly differentiated SNPs (p<0.02), while case-control association discovered up to 152 SNPs (p<0.001). Both analyses were able to successfully differentiate individuals between the orange and black tissue colour morphotypes characteristic of this species. BLAST searches revealed that five of these SNPs were associated with a melanin biosynthesis pathway, demonstrating their biological relevance. This study has produced highly informative SNP and population genomic data in P. margaritifera, and using the same approach promises to be of substantial value to a range of other non-model, broadcast-spawning or marine invertebrate taxa.

Domestication as a Novel Approach for Improving the Cultivation of Calanoid Copepods: A Case Study with Parvocalanus crassirostris

Calanoid copepods are an important food source for most fish larvae. Their role as a natural prey item means that it is important to develop culture technology for copepods to meet the requirements of larvae culture in aquaculture hatcheries. Copepods have been cultured successfully for some time; however, the implications of long-term cultivation or domestication on copepod life history traits have not yet been assessed. Therefore, the aim of this study was to determine if the domesticated and wild populations of Parvocalanus crassirostris are phenotypically or physiologically different. Wild-caught P. crassirostris were compared to a long-held domesticated strain (>2 years) for size of developmental stage, fecundity and lifespan of adult females, culture density, ingestion rates, faecal pellet production and fatty acid profiles. The domesticated strain was significantly different from the wild strain in size (eggs, nauplii, copepodites and adults were larger in the domesticated strain), egg production (112.3 ± 1.8 eggs female^-1 vs. 64.6 ± 3.3 eggs female^-1) and adult female lifespan (8.8 ± 0.1 days vs., 7.5 ± 0.1 days). At 1, 3 and 5 adults mL^-1, the domesticated strain performed significantly better than the wild strain in egg production (4189.8 ± 61.2, 11224.0 ± 71.7 and 21860.6 ± 103.6 eggs vs. 1319.5 ± 54.3, 2374.5 ± 80.9 and 4933.8 ± 269.5 eggs, respectively) and mean daily mortality rate (5.6% across all densities vs. 22.9 ± 1.6, 29.8 ± 1.2 and 31.3 ± 1.3%, respectively). The domesticated strain had significantly higher ingestion rates than the wild stain (888.4 ± 9.9 ng C l^-1 and 775.3 ± 11.2 ng C l^-1, respectively), while faecal pellet production was not significantly different between strains. Fatty acid profiles indicated higher levels (as % of total fatty acid) of long-chain polyunsaturated fatty acids in the domesticated strain (36.2±0.4%) than the wild strain (16.1±0.1%). Overall, this study found that the reproductive capacity and tolerance to the culture environment of the calanoid P. crassirostris have improved significantly due to domestication.

Oceanographic Currents and Local Ecological Knowledge Indicate, and Genetics Does Not Refute, a Contemporary Pattern of Larval Dispersal for The Ornate Spiny Lobster, Panulirus ornatus in the South-East Asian Archipelago

Here we utilize a combination of genetic data, oceanographic data, and local ecological knowledge to assess connectivity patterns of the ornate spiny lobster Panulirus ornatus (Fabricius, 1798) in the South-East Asian archipelago from Vietnam to Australia. Partial mitochondrial DNA control region and 10 polymorphic microsatellites did not detect genetic structure of 216 wild P. ornatus samples from Australia, Indonesia and Vietnam. Analyses show no evidence for genetic differentiation among populations (mtDNA control region sequences Φ_ST = -0.008; microsatellite loci F_ST = 0.003). A lack of evidence for regional or localized mtDNA haplotype clusters, or geographic clusters of microsatellite genotypes, reveals a pattern of high gene flow in P. ornatus throughout the South-East Asian Archipelago. This lack of genetic structure may be due to the oceanography-driven connectivity of the pelagic lobster larvae between spawning grounds in Papua New Guinea, the Philippines and, possibly, Indonesia. The connectivity cycle necessitates three generations. The lack of genetic structure of P. ornatus population in the South-East Asian archipelago has important implications for the sustainable management of this lobster in that the species within the region needs to be managed as one genetic stock.

Mitogenomic phylogeny of the Percichthyidae and Centrarchiformes (Percomorphaceae): comparison with recent nuclear gene-based studies and simultaneous analysis

Delineation of the fish family Percichthyidae (Percomorphaceae) has a long and convoluted history, with recent morphological-based studies restricting species members to South American and Australian freshwater and catadromous temperate perches. Four recent nuclear gene-based phylogenetic studies, however, found that the Percichthyidae was not monophyletic and was nested within a newly discovered inter-familial clade of Percomorphaceae, the Centrarchiformes, which comprises the Centrarchidae and 12 other families. Here, we reexamined the systematics of the Percichthyidae and Centrarchiformes based on new mitogenomic information. Our mitogenomic results are globally congruent with the recent nuclear gene-based studies although the overall amount of phylogenetic signal of the mitogenome is lower. They do not support the monophyly of the Percichthyidae, because the catadromous genus Percalates is not exclusively related to the freshwater percichthyids. The Percichthyidae (minus Percalates) and Percalates belong to a larger clade, equivalent to the Centrarchiformes, but their respective sister groups are unresolved. Because all recent analyses recover a monophyletic Centrarchiformes but with substantially different intra-relationships, we performed a simultaneous analysis for a character set combining the mitogenome and 19 nuclear genes previously published, for 22 centrarchiform taxa. This analysis furthermore indicates that the Centrarchiformes are divided into three lineages and the superfamily Cirrhitoidea is monophyletic as well as the temperate and freshwater centrarchiform perch-like fishes. It also clarifies some of the relationships within the freshwater Percichthyidae.

Differential IgE binding to isoallergens from Asian seabass (Lates calcarifer) in children and adults

Fish allergy is a common food allergy, with prevalence rates in the general population ranging between 0.2% and 2.3%. In both adults and children fish ranks in the top eight foods known to cause IgE mediated food allergy. Fish allergy is rarely outgrown and individuals with fish allergy may be allergic to some but not all species of fish. Whilst fish allergy occurs around the world, the characterization of allergenic components of individual species of fish has been largely confined to Northern hemisphere and European fish species. To date allergy to commonly consumed fish in the Asian-Pacific region including barramundi (Asian seabass; Lates calcarifer) have been less well investigated. The aim of this study was to identify and characterize allergenic proteins from barramundi in both fish allergic adult and pediatric patients. Serum from 17 fish allergic adults and children from Australia were characterized by immunoblotting and enzyme linked immunosorbent assays (ELISA) against raw and heated barramundi. Molecular analysis of identified allergens included genetic sequencing and generation of recombinant isoallergens. Two novel parvalbumin isoforms of the β-type were identified as the only allergens in barramundi and subsequently designated as Lat c 1.0101 and Lat c 1.0201 by the International Union of Immunological Societies. These two isoallergens do not differ in their ability to bind IgE antibodies, but are differentially expressed in barramundi tissue. This study characterized two novel heat stable parvalbumin allergens from barramundi, with differential IgE binding capacity between adults and pediatric patients.

Effects of growth phase and nitrogen starvation on expression of fatty acid desaturases and fatty acid composition of Isochrysis aff. galbana (TISO)

Very long-chain polyunsaturated fatty acids (VLC-PUFAs) are important dietary requirements for maintaining human health. Many marine microalgae are naturally high in ω-3 VLC-PUFAs, however, the molecular mechanisms underpinning fatty acid (FA) desaturation and elongation in algae are poorly understood. An advanced molecular understanding would facilitate improvements of this nascent industry. We aimed to investigate expression responses of four front-end fatty acid desaturase genes and downstream effects on FA profiles to nitrogen limitation and cultivation growth stage in Isochrysis aff. galbana (TISO). Cultures were grown in nitrogen-replete and -deplete medium; samples were harvested during logarithmic, late logarithmic and stationary growth phases to analyse FA content/composition and gene expression of ∆(6)-, ∆(8)-, ∆(5)- and ∆(4)-desaturases (d6FAD (putative), d8FAD, d5FAD and d4FAD, respectively). d6FAD (putative) exhibited no differential expression, while d8FAD, d5FAD and d4FAD were significantly upregulated during logarithmic growth of nutrient-replete cultures, coinciding with rapid cell division. In conclusion, it is demonstrated that expression of some FADs in I. aff. galbana varies with culture age and nitrogen status which has downstream consequences on FA desaturation levels. This has implications for the commercial production of VLC-PUFAs where a trade-off between total lipid yield and VLC-PUFAs has to be made.

Influence of drainage divides versus arid corridors on genetic structure and demography of a widespread freshwater turtle, Emydura macquarii krefftii, from Australia

The influence of Pleistocene climatic cycles on Southern Hemisphere biotas is not yet well understood. Australia's eastern coastal margin provides an ideal setting for examining the relative influence of landscape development, sea level fluctuation, and cyclic climatic aridity on the evolution of freshwater biodiversity. We examined the impact of climatic oscillations and physical biogeographic barriers on the evolutionary history of the wide-ranging Krefft's river turtle (Emydura macquarii krefftii), using range-wide sampling (649 individuals representing 18 locations across 11 drainages) and analysis of mitochondrial sequences (∼1.3-kb control region and ND4) and nuclear microsatellites (12 polymorphic loci). A range of phylogeographic (haplotype networks, molecular dating), demographic (neutrality tests, mismatch distributions), and population genetic analyses (pairwise F ST, analysis of molecular variance, Bayesian clustering analysis) were implemented to differentiate between competing demographic (local persistence vs. range expansion) and biogeographic (arid corridor vs. drainage divide) scenarios. Genetic data reveal population genetic structure in Krefft's river turtles primarily reflects isolation across drainage divides. Striking north-south regional divergence (2.2% ND4 p-distance; c. 4.73 Ma, 95% higher posterior density (HPD) 2.08-8.16 Ma) was consistent with long-term isolation across a major drainage divide, not an adjacent arid corridor. Ancient divergence among regional lineages implies persistence of northern Krefft's populations despite the recurrent phases of severe local aridity, but with very low contemporary genetic diversity. Stable demography and high levels of genetic diversity are inferred for southern populations, where aridity was less extreme. Range-wide genetic structure in Krefft's river turtles reflects contemporary and historical drainage architecture, although regional differences in the extent of Plio-Pleistocene climatic aridity may be reflected in current levels of genetic diversity.

Genome-wide SNP validation and mantle tissue transcriptome analysis in the silver-lipped pearl oyster, Pinctada maxima

Pearl oysters are not only farmed for their gemstone quality pearls worldwide, but they are also becoming important model organisms for investigating genetic mechanisms of biomineralisation. Despite their economic and scientific significance, limited genomic resources are available for this important group of bivalves, hampering investigations into identifying genes that regulate important pearl quality traits and unique biological characteristics (i.e. biomineralisation). The silver-lipped pearl oyster, Pinctada maxima, is one species where there is interest in understanding genes that regulate commercially important pearl traits, but presently, there is a dearth of genomic information. The objective of this study was to develop and validate a large number of type I genome-wide single nucleotide polymorphisms (SNPs) for P. maxima suitable for high-throughput genotyping. In addition, sequence annotations and Gene Ontology terms were assigned to a large mantle tissue 454 expressed sequence tag assembly (96,794 contigs) and information on known bivalve biomineralisation genes was incorporated into SNP discovery. The SNP discovery effort resulted in the de novo identification of 172,625 SNPs, of which 9,108 were identified as high value [minor allele frequency (MAF)≥ 0.15, read depth  ≥ 8]. Validation of 2,782 of these SNPs using Illumina iSelect Infinium genotyping technology returned some of the highest assay conversion (86.6 %) and validation (59.9 %; mean MAF 0.28) rates observed in aquaculture species to date. Genomic resources presented here will be pivotal to future research investigating the biological mechanisms behind biomineralisation and will form a strong foundation for genetic selective breeding programs in the P. maxima pearling industry.

A high-density SNP genetic linkage map for the silver-lipped pearl oyster, Pinctada maxima: a valuable resource for gene localisation and marker-assisted selection

Background

The silver-lipped pearl oyster, Pinctada maxima, is an important tropical aquaculture species extensively farmed for the highly sought "South Sea" pearls. Traditional breeding programs have been initiated for this species in order to select for improved pearl quality, but many economic traits under selection are complex, polygenic and confounded with environmental factors, limiting the accuracy of selection. The incorporation of a marker-assisted selection (MAS) breeding approach would greatly benefit pearl breeding programs by allowing the direct selection of genes responsible for pearl quality. However, before MAS can be incorporated, substantial genomic resources such as genetic linkage maps need to be generated. The construction of a high-density genetic linkage map for P. maxima is not only essential for unravelling the genomic architecture of complex pearl quality traits, but also provides indispensable information on the genome structure of pearl oysters.

Results

A total of 1,189 informative genome-wide single nucleotide polymorphisms (SNPs) were incorporated into linkage map construction. The final linkage map consisted of 887 SNPs in 14 linkage groups, spans a total genetic distance of 831.7 centimorgans (cM), and covers an estimated 96% of the P. maxima genome. Assessment of sex-specific recombination across all linkage groups revealed limited overall heterochiasmy between the sexes (i.e. 1.15:1 F/M map length ratio). However, there were pronounced localised differences throughout the linkage groups, whereby male recombination was suppressed near the centromeres compared to female recombination, but inflated towards telomeric regions. Mean values of LD for adjacent SNP pairs suggest that a higher density of markers will be required for powerful genome-wide association studies. Finally, numerous nacre biomineralization genes were localised providing novel positional information for these genes.

Conclusions

This high-density SNP genetic map is the first comprehensive linkage map for any pearl oyster species. It provides an essential genomic tool facilitating studies investigating the genomic architecture of complex trait variation and identifying quantitative trait loci for economically important traits useful in genetic selection programs within the P. maxima pearling industry. Furthermore, this map provides a foundation for further research aiming to improve our understanding of the dynamic process of biomineralization, and pearl oyster evolution and synteny.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-810) contains supplementary material, which is available to authorized users.

Temperature: a prolonged confounding factor on cholinesterase activity in the tropical reef fish Acanthochromis polyacanthus

Cholinesterase activity usually decreases in fish exposed to anticholinesterase compounds such as organophosphate and carbamate pesticides. Here we show that tropical reef fish Acanthochromis polyacanthus (or spiny damsel) also exhibits a decrease in ChE activity when exposed to elevated temperature from 28°C to 32°C or 34°C after 4 days. We further demonstrate that the decline persists even after 7 days of recovery at control temperature. This is the first report of a drop in ChE activity in fish as temperature increases. Our results strongly suggest the need for long-term monitoring of water temperature in the field prior to sampling A. polyacanthus for toxicology studies, as temperature is a prolonged and confounding factor for ChE activity in this species.

Next-generation transcriptome profiling reveals insights into genetic factors contributing to growth differences and temperature adaptation in Australian populations of barramundi (Lates calcarifer)

Identification of genetically-regulated adaptation in fish is a precursor to understanding how populations will respond to future climate induced stressors like temperature. Australian populations of barramundi (Lates calcarifer) show strong evidence of local adaptation to temperature. However, the phenotypic consequences of this adaptation are unknown and the genetic mechanisms underlying this response are poorly understood. In the current study, two populations of barramundi from temperature extremes of the species Australian distribution were communally reared at cool (22°C), control (28°C) and hot (36°C) water temperatures for 3.5months. Southern populations of barramundi originating from a cooler environment grew significantly faster at 22°C than northern populations of warm adapted barramundi. However, no difference in population growth was present at either 28°C or 36°C. The underlying transcriptome profile of barramundi was examined via Illumina mRNA deep sequencing to determine the major contributing gene categories giving rise to phenotypic differences in barramundi population growth. Gene ontology (GO) analysis revealed enrichment in categories relating to the regulation of peptidase activity as well as microtubule, cytoplasmic and cellular metabolic based processes. Further analysis of the GO category "microtubule based process" with associated genes from the "response to stress" category revealed an apparent re-organization of cytoskeletal elements in response to an induced cold stress in northern barramundi reared at 22°C, when compared with northern barramundi reared at 36°C. Between southern barramundi and northern barramundi reared at 36°C, an analysis of the "endopeptidase inhibitor activity" GO category in conjunction with stress genes indicated a suppression of the complement system in southern barramundi along with an increase in the cellular stress response. The results of the present study show that southern populations of barramundi exhibit underlying molecular adaptation to cooler water temperatures, but still retain a tolerance for warm water temperatures. Furthermore, GO profiling has revealed groups of genes that underlie population differences in temperature tolerance as a means to prioritize the analysis of differential gene expression in studies of local adaptation in the future.

The gene expression response of the catadromous perciform barramundi Lates calcarifer to an acute heat stress

The acute heat-shock response of the tropical estuarine fish species barramundi Lates calcarifer as indicated by the expression of genes within stress (hsp 90AA, hsp 90AB, hsp 70 and hsc 70), metabolic (cisy, cco II and ldh) and growth (igf1 and mstn 1) related pathways was examined following an increase in water temperature from 28 to 36° C over 30 min. Lates calcarifer were maintained at the acute stress temperature of 36° C for 1 h before being returned to 28° C and allowed to recover at this temperature for a further 2 weeks. Muscle tissue sampling over the experimental period allowed for the expression quantification of stress, metabolic and growth-related genes via quantitative real-time polymerase chain reaction (qrt-PCR) where a robust and reliable normalization approach identified both α-tub and Rpl8 as appropriate genes for the analysis of gene expression in response to an acute heat stress. hsp90AA and hsp70 of the inducible heat-shock response pathway showed a massive up-regulation of gene expression in response to heat stress, whilst the constitutive heat-shock genes hsp90AB and hsp70 showed no change over the course of the experiment and a small increase after 2 weeks of recovery, respectively. Of the three genes representing the metabolic pathway (cisy, cco II and ldh) only cco II changed significantly showing a decrease in gene expression, which may suggest a small suppression of aerobic metabolism. igf1 of the growth pathway showed no significant differences in response to an acute heat stress, whilst mstn1 increased at the beginning of the heat stress but returned to basal levels soon after. Overall, the results demonstrate that an acute heat stress in L. calcarifer caused a significant increase in the expression of genes from the stress response pathway and a possible decrease in aerobic metabolism with only relatively minor changes to the growth pathway highlighting the hardy nature of L. calcarifer and its resilience in coping with sudden temperature changes routinely encountered within its natural environment.

A robust flow-cytometric protocol for assessing growth rate of hatchery-reared barramundi Lates calcarifer larvae

In this study, a flow-cytometric cell cycle analysis method to assess instantaneous growth rate of whole larvae of the Australian barramundi Lates calcarifer was developed and validated. High-resolution DNA measurements of either fresh, frozen or RNAlater-preserved larvae (gap0-gap1, G(0) -G(1), coefficient of variation (c.v.) < 3, 4 and 5%, respectively) enabled the deconvolution of the DNA histogram and assignment of the proportion of nuclei into cell cycle compartments G(0) -G(1), S (DNA synthesis) and G(2) -M (Gap2-Mitosis). This technique can be also used for individual fish tissues such as brain, liver, fin and muscle. For the first time, the combined proportion of replicating nuclei (into S and G(2) -M phases) of whole fish larvae and absolute growth rate in length (mm day(-1)) has been correlated in commercial aquaculture conditions. Fast growing L. calcarifer larvae had an overall hyperplasia advantage as indicated by a greater proportion of cells in the S+G(2) -M phase compared with slow growing larvae, which might explain the increasing differences in size during culture. In a fasting trial, larvae ceased growth while maintaining the constant initial rates of cell division throughout a 6 day period. For a highly fed fast growing control group, cell division rates significantly increased after day 4. Flow-cytometric cell cycle analysis of whole fish larvae may provide fish biologists and aquaculturists with a better understanding of how cell division rates influence early growth in natural and artificial environments.