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.

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 mud crab, Scylla serrata, is currently being evaluated for its feasibility for mass aquaculture in Australia. As part of the evaluation process, pathogens that can affect this species need to be identified. This research note describes a possible new parasite that infects the eggs of S. serrata. The parasite was discovered in two separate cases (three months apart) in a broodstock research program and appears to cause 100% egg mortality. The parasite consists of a cluster of cells with rhizoids that appear to function as an anchorage and a feeding organ. The individual cells range from 3 to 6 microm with a single nucleus. The parasite could not be classified to a phylum by morphology alone. However, BLAST analysis of the DNA sequence from a PCR amplicon generated using universal 18S ribosomal RNA gene primers indicated similarity to pathogenic thraustochytrids, Dermocystidium sp. and Rhinosporidium seeberi. All except R. seeberi are protistan parasites of fish and crustaceans. A neighbour-joining phylogenetic tree confirmed the homology with the thraustochytrids; however, further molecular techniques need to be exploited for absolute classification of this new parasite.

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.

Evolutionary relationships in the black-lip pearl oyster Pinctada margaritifera which is highly valued for pearl production remain poorly understood. This species possesses an 18,000 km Indo-Pacific natural distribution, and its current description includes six subspecies defined exclusively on morphological characters. To evaluate its taxonomic identity using molecular data, 14 populations in both the Indian and Pacific Oceans (n = 69), and the congeneric taxa P. maxima and P. mazatlanica (n = 29 and n = 10, respectively) were sampled. Phylogenomic reconstruction was carried out using both 8,308 genome-wide SNPs and 10,000 dominant loci (DArTseq PAVs). Reconstructions using neighbour-joining (Nei's 1972 distances), maximum likelihood and Bayesian approaches all indicate that the taxonomy of P. margaritifera is quite complex, with distinct evolutionary significant units (ESUs) identified within Tanzanian and Iranian populations. Contrastingly, phylogenies generated for Pacific Ocean oysters resolved a large monophyletic clade, suggesting little support for two current morphological subspecies classifications. Furthermore, P. mazatlanica formed a basal clade closest to French Polynesian P. margaritifera, suggesting it may be conspecific. Collectively, these findings provide evidence that P. margaritifera comprises a species complex, perhaps as a result of population fragmentation and increased divergence at range limits.

Asian seabass (or commonly known as barramundi), Lates calcarifer, is a bony euryhaline teleost from the Family Latidae, inhabiting nearshore, estuarine, and marine connected freshwaters throughout the tropical Indo-West Pacific region. The species is catadromous, whereby adults spawn in salinities between 28 and 34 ppt at the mouth of estuaries, with resultant juveniles usually moving into brackish and freshwater systems to mature, before returning to the sea to spawn again as adults. The species lives in both marine and freshwater habitats and can move quickly between the two; thus, the species' ability to tolerate changes in salinity makes it a good candidate for studying the salinity acclimation response in teleosts. In this study, the transcriptome of two major osmoregulatory organs (gills and kidneys) of young juvenile Asian seabass reared in freshwater and seawater were compared. The euryhaline nature of Asian seabass was found to be highly pliable and the moldability of the trait was further confirmed by histological analyses of gills and kidneys. Differences in major expression pathways were observed, with differentially expressed genes including those related to osmoregulation, tissue/organ morphogenesis, and cell volume regulation as central to the osmo-adaptive response. Additionally, genes coding for mucins were upregulated specifically under saline conditions, whereas several genes important for growth and development, as well as circadian entrainment were specifically enriched in fish reared in freshwater. Routing of the circadian rhythm mediated by salinity changes could be the initial step in salinity acclimation and possibly migration in euryhaline fish species such as the Asian seabass.

Pseudomonas aeruginosa ST308 clone has been reported to carry carbapenemase genes such as blaIMP and blaVIM but has been rarely associated with blaNDM-1. A total of 199 P. aeruginosa ST308 clinical and environmental isolates obtained between April 2019 and November 2020 from a tertiary-care hospital in Singapore were characterized using whole-genome sequencing. In addition, 71 blaNDM-1-positive ST308 whole-genome sequences from two other local tertiary-care hospitals in Singapore and 83 global blaNDM-1-negative ST308 whole-genome sequences in public databases were included to assess phylogenetic relationships and perform genome analyses. Phylogenetic analysis and divergent time estimation revealed that blaNDM-1-positive P. aeruginosa ST308 was introduced into Singapore in 2005 (95 % highest posterior density: 2001 to 2008). Core genome, resistome, and analyses of all local blaNDM-1-positive ST308 isolates showed chromosomal integration of multiple antibiotic resistance genes (ARGs) [aac(3)-Id, aac(6')-Il, aadA6, aadA11, dfrB5, msr(E), floR, sul2, and qnrVC1], which was absent in global blaNDM-1-negative ST308 sequences. Most ARGs and virulence genes were conserved across isolates originating from the three different local hospitals. Close genetic relatedness of the blaNDM-1-positive ST308 clinical and environmental isolates suggests cocirculation between the hospital environment and human hosts with the hospital environment as a potential reservoir. Core genome single nucleotide polymorphism analyses revealed possible clonal transmission of blaNDM-1-positive ST308 isolates between the three hospitals over 7 years. Bloodstream isolates accounted for six of 95 (6.3%) clinical isolates. This study reports the introduction of a pathogenic blaNDM-1-positive P. aeruginosa ST308 more than a decade ago in Singapore and warrants surveillance for wider dissemination. IMPORTANCE P. aeruginosa is a Gram-negative opportunistic pathogen ubiquitously found in the environment and a major cause of nosocomial infections. While the P. aeruginosa ST308 clone has been known to bear blaIMP and blaVIM among global isolates, reports of blaNDM-1-positive P. aeruginosa ST308 are rare. The local blaNDM-1-positive P. aeruginosa ST308 isolates detected in this study appear to be unique to this region, with evidence of chromosomal acquisition of multiple ARGs compared to global blaNDM-1-negative P. aeruginosa ST308 isolates. Surveillance in Singapore and beyond for dissemination is essential to determine whether existing measures are sufficient to control the spread of this ST308 clone.

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.

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.

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.

The Pacific whiteleg shrimp, Litopenaeus vannamei, is the most farmed aquaculture species worldwide with global production exceeding 3 million tonnes annually. Litopenaeus vannamei has been the focus of many selective breeding programs aiming to improve growth and disease resistance. However, these have been based primarily on phenotypic measurements and omit potential gains by integrating genetic selection into existing breeding programs. Such integration of genetic information has been hindered by the limited available genomic resources, background genetic parameters and knowledge on the genetic architecture of commercial traits for L. vannamei. This study describes the development of a comprehensive set of genomic gene-based resources including the identification and validation of 234,452 putative single nucleotide polymorphisms in-silico, of which 8,967 high value SNPs were incorporated into a commercially available Illumina Infinium ShrimpLD-24 v1.0 genotyping array. A framework genetic linkage map was constructed and combined with locus ordering by disequilibrium methodology to generate an integrated genetic map containing 4,817 SNPs, which spanned a total of 4552.5 cM and covered an estimated 98.12% of the genome. These gene-based genomic resources will not only be valuable for identifying regions underlying important L. vannamei traits, but also as a foundational resource in comparative and genome assembly activities.

Receptor-like cytoplasmic kinases (RLCKs) in plants belong to the super family of receptor-like kinases (RLKs). These proteins show homology to RLKs in kinase domain but lack the transmembrane domain. Some of the functionally characterized RLCKs from plants have been shown to play roles in development and stress responses. Previously, 149 and 187 RLCK encoding genes were identified from Arabidopsis and rice, respectively. By using HMM-based domain structure and phylogenetic relationships, we have identified 379 OsRLCKs from rice. OsRLCKs are distributed on all 12 chromosomes of rice and some members are located on duplicated chromosomal segments. Several OsRLCKs probably also undergo alternative splicing, some having evidence only in the form of gene models. To understand their possible functions, expression patterns during landmark stages of vegetative and reproductive development as well as abiotic and biotic stress using microarray and MPSS-based data were analyzed. Real-time PCR-based expression profiling for a selected few genes confirmed the outcome of microarray analysis. Differential expression patterns observed for majority of OsRLCKs during development and stress suggest their involvement in diverse functions in rice. Majority of the stress-responsive OsRLCKs were also found to be localized within mapped regions of abiotic stress QTLs. Outcome of this study would help in selecting organ/development stage specific OsRLCK genes/targets for functional validation studies.

Myostatin (MSTN) is a member of the transforming growth factor-beta superfamily that negatively regulates growth of skeletal muscle tissue. The gene encoding for the MSTN peptide is a consolidate candidate for the enhancement of productivity in terrestrial livestock. This gene potentially represents an important target for growth improvement of cultured finfish.

The understanding of algal phylogeny is being impeded by an unknown number of events of horizontal gene transfer (HGT), and primary and secondary/tertiary endosymbiosis. Through these events, previously heterotrophic eukaryotes developed photosynthesis and acquired new biochemical pathways. Acetyl-CoA carboxylase (ACCase) is a key enzyme in the fatty acid synthesis and elongation pathways in algae, where ACCase exists in two locations (cytosol and plastid) and in two forms (homomeric and heteromeric). All algae contain nucleus-encoded homomeric ACCase in the cytosol, independent of the origin of the plastid. Nucleus-encoded homomeric ACCase is also found in plastids of algae that arose from a secondary/tertiary endosymbiotic event. In contrast, plastids of algae that arose from a primary endosymbiotic event contain heteromeric ACCase, which consists of three nucleus-encoded and one plastid-encoded subunits. These properties of ACCase provide the potential to inform on the phylogenetic relationships of hosts and their plastids, allowing different hypothesis of endosymbiotic events to be tested. Alveolata (Dinoflagellata and Apicomplexa) and Chromista (Stramenopiles, Haptophyta and Cryptophyta) have traditionally been grouped together as Chromalveolata, forming the red lineage. However, recent genetic evidence groups the Stramenopiles, Alveolata and green plastid containing Rhizaria as SAR, excluding Haptophyta and Cryptophyta. Sequences coding for plastid and cytosol targeted homomeric ACCases were isolated from Isochrysis aff. galbana (TISO), Chromera velia and Nannochloropsis oculata, representing three taxonomic groups for which sequences were lacking. Phylogenetic analyses show that cytosolic ACCase strongly supports the SAR grouping. Conversely, plastidial ACCase groups the SAR with the Haptophyta, Cryptophyta and Prasinophyceae (Chlorophyta). These two ACCase based, phylogenetic relationships suggest that the plastidial homomeric ACCase was acquired by the Haptophyta, Cryptophyta and SAR, before the photosynthetic Rhizaria acquired their green plastid. Additionally, plastidial ACCase was derived by HGT from an ancestor or relative of the Prasinophyceae and not by duplication of cytosolic ACCase.

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.

The Asian seabass is an important marine food fish that has been cultured for several decades in Asia Pacific. However, the lack of a high quality reference genome has hampered efforts to improve its selective breeding. A 3D BAC pool set generated in this study was screened using 22 SSR markers located on linkage group 2 which contains a growth-related QTL region. Seventy-two clones corresponding to 22 FPC contigs were sequenced by Illumina MiSeq technology. We co-assembled the MiSeq-derived scaffolds from each FPC contig with error-corrected PacBio reads, resulting in 187 sequences covering 9.7 Mb. Eleven genes annotated within this region were found to be potentially associated with growth and their tissue-specific expression was investigated. Correlation analysis demonstrated that SNPs in ctsb, skp1 and ppp2ca can be potentially used as markers for selecting fast-growing fingerlings. Conserved syntenies between seabass LG2 and five other teleosts were identified. This study i) provided a 10 Mb targeted genome assembly; ii) demonstrated NGS of BAC pools as a potential approach for mining candidates underlying QTLs of this species; iii) detected eleven genes potentially responsible for growth in the QTL region; and iv) identified useful SNP markers for selective breeding programs of Asian seabass.

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.

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.

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 [Ca2+] transients and rapidly rising tetanic force responses. To achieve such high amplitude and high frequency of Ca2+ transients requires a high abundance of sarcoplasmic/endoplasmic reticulum Ca2+ ATPase (SERCA) to rapidly clear cytoplasmic Ca2+ between each transient and an efficient Ca2+ release system consisting of the Ryanodine Receptor (RyR), and voltage gated Ca2+ channels (CaVs). 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 CaV 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.

The development of genome-wide genotyping resources has provided terrestrial livestock and crop industries with the unique ability to accurately assess genomic relationships between individuals, uncover the genetic architecture of commercial traits, as well as identify superior individuals for selection based on their specific genetic profile. Utilising recent advancements in de-novo genome-wide genotyping technologies, it is now possible to provide aquaculture industries with these same important genotyping resources, even in the absence of existing genome assemblies. Here, we present the development of a genome-wide SNP assay for the Black Tiger shrimp (Penaeus monodon) through utilisation of a reduced-representation whole-genome genotyping approach (DArTseq).