The complete sequence of the mitochondrial genome of the crustacean Penaeus monodon: are malacostracan crustaceans more closely related to insects than to branchiopods?

The complete sequence of the mitochondrial genome of the giant tiger prawn, Penaeus monodon (Arthropoda, Crustacea, Malacostraca), is presented. The gene content and gene order are identical to those observed in Drosophila yakuba. The overall AT composition is lower than that observed in the known insect mitochondrial genomes, but higher than that observed in the other two crustaceans for which complete mitochondrial sequence is available. Analysis of the effect of nucleotide bias on codon composition across the Arthropoda reveals a trend with the crustaceans represented showing the lowest proportion of AT-rich codons in mitochondrial protein genes. Phylogenetic analysis among arthropods using concatenated protein-coding sequences provides further support for the possibility that Crustacea are paraphyletic. Furthermore, in contrast to data from the nuclear gene EF1alpha, the first complete sequence of a malacostracan mitochondrial genome supports the possibility that Malacostraca are more closely related to Insecta than to Branchiopoda.

Mitochondrial DNA variation in Indo-Pacific populations of the giant tiger prawn, Penaeus monodon

Surveys of mitochondrial DNA (mtDNA) variation in the giant tiger prawn, Penaeus monodon, using restriction fragment length polymorphisms have provided the first clear evidence that the Indo-West Pacific region is a site of accumulation of genetic diversity rather than a site of origin of genetic diversity. No haplotyes were found in common between a group of five southeast African populations and a group of five Australian (including Western Australia) and three southeast Asian populations. The dominant haplotype was different in the Australian and southeast Asian population groups. Genetic diversity (pi) was greatest in Indonesia (pi averaged 0.05), less in the Philippines and Australia (pi averaged 0.01), and markedly less in the southeast African and the West Australian populations (pi averaged 0.003). The high diversity of the southeast Asian populations resulted from the occurrence in those populations of a set of haplotypes found only in southeast Asia but derived from the southeast African haplotypes. These genetic variants therefore evolved in the Indian Ocean and later migrated into the Indo-West Pacific region. Low genetic variation in the geographically marginal populations in southeast Africa and Western Australia is considered to be the result of bottlenecks, but mismatch distributions suggest that large population sizes have been maintained in Indonesian populations for long periods.

Genomics in marine monitoring: new opportunities for assessing marine health status

This viewpoint paper explores the potential of genomics technology to provide accurate, rapid, and cost efficient observations of the marine environment. The use of such approaches in next generation marine monitoring programs will help achieve the goals of marine legislation implemented world-wide. Genomic methods can yield faster results from monitoring, easier and more reliable taxonomic identification, as well as quicker and better assessment of the environmental status of marine waters. A summary of genomic methods that are ready or show high potential for integration into existing monitoring programs is provided (e.g. qPCR, SNP based methods, DNA barcoding, microarrays, metagenetics, metagenomics, transcriptomics). These approaches are mapped to existing indicators and descriptors and a series of case studies is presented to assess the cost and added value of these molecular techniques in comparison with traditional monitoring systems. Finally, guidelines and recommendations are suggested for how such methods can enter marine monitoring programs in a standardized manner.

Phospholipase A2 in Cnidaria

Phospholipase A2 (PLA2) is an enzyme present in snake and other venoms and body fluids. We measured PLA2 catalytic activity in tissue homogenates of 22 species representing the classes Anthozoa, Hydrozoa, Scyphozoa and Cubozoa of the phylum Cnidaria. High PLA2 levels were found in the hydrozoan fire coral Millepora sp. (median 735 U/g protein) and the stony coral Pocillopora damicornis (693 U/g) that cause skin irritation upon contact. High levels of PLA2 activity were also found in the acontia of the sea anemone Adamsia carciniopados (293 U/g). Acontia are long threads containing nematocysts and are used in defense and aggression by the animal. Tentacles of scyphozoan and cubozoan species had high PLA2 activity levels: those of the multitentacled box jellyfish Chironex fleckeri contained 184 U/g PLA2 activity. The functions of cnidarian PLA2 may include roles in the capture and digestion of prey and defense of the animal. The current observations support the idea that cnidarian PLA2 may participate in the sting site irritation and systemic envenomation syndrome resulting from contact with cnidarians.

Gene flow and population history in high dispersal marine invertebrates: mitochondrial DNA analysis of Holothuria nobilis (Echinodermata: Holothuroidea) populations from the Indo-Pacific

The sea cucumber, Holothuria nobilis, has a long-lived planktotrophic larvae, and previous allozyme surveys have suggested that high dispersal is realized. In contrast, recent ecological studies indicate that dispersal is low. To reconcile these data, and to investigate the evolution of this Indo-Pacific species, we screened geographical variation in 559 bp of a mitochondrial gene (COI) in 360 samples from the Australasian region and La Réunion. Sequences from La Réunion differed by > 7% from others and may constitute another species. Haplotype diversity in other samples was high (0.942, SD = 0.007), but haplotypes were closely related (mean nucleotide diversity: 0.0075, SD = 0.0041). AMOVA, pairwise FST values and exact tests did not detect significant population structure. Nested clade analysis showed that one of two main clades was over-represented in west Australia, whereas the other was more common in the northern Great Barrier Reef. Isolation-by-distance was identified as the main determinant of population structure at several clade levels. Contiguous range expansion was inferred for evolutionary older clade levels and this may correspond to a late Pleistocene (88 000-193 000 years ago) population expansion inferred from haplotype mismatch distributions. Thus, the population genetic structures detected are likely to be formed prior to the last ice age, with some indications for high dispersal on shorter time scales.

A threat to coral reefs multiplied? Four species of crown-of-thorns starfish

In the face of ever-increasing threats to coral reef ecosystems, it is essential to understand the impact of natural predators in order to devise appropriate management strategies. Destructive population explosions of the crown-of-thorns starfish Acanthaster planci have devastated coral reefs throughout the Indo-Pacific for decades. But despite extensive research, the causes of outbreaks are still unclear. An important consideration in this research is that A. planci has been regarded as a single taxonomic entity. Using molecular data from its entire distribution, we find that A. planci is in fact a species complex. This discovery has important consequences for future coral reef research, and might prove critical for successful reef conservation management.

The Effects of Sperm Concentration, Sperm:Egg Ratio, and Gamete Age on Fertilization Success in Crown-of-Thorns Starfish (Acanthaster planci) in the Laboratory

Laboratory experiments varying gamete concentrations and gamete age demonstrated significant reductions in fertilization success of the starfish Acanthaster planci (L.) with decreasing sperm concentration and increasing age of both eggs and sperm. The effect of aging in sperm was faster than that of eggs, and the speed of sperm aging increased with increasing dilution of sperm. Fertilization success was high over a wide range of sperm: egg ratios but declined rapidly at ratios less than 50, particularly at low sperm concentrations. A. planci gametes aged more slowly, and the loss of fertilizing capacity of sperm with dilution (the respiratory dilution effect) was far less, than in sea urchins. These characteristics provide a mechanism for enhanced fertilization success at given sperm concentrations and at greater distances and times from the point of gamete release, and may explain the higher fertilization rates achieved over longer distances in the wild by A. planci relative to sea urchins. Gametes would remain competent for longer periods at more dilute concentrations and so better achieve long-distance fertilization. Gametes obtained at the end of the breeding season were qualitatively different from those obtained early in the breeding season and showed reduced fertilization success for a given combination of variables, and different fertilization dynamics.

The distribution, abundance, and the effects of fire on mound building termites (Trinervitermes and Cubitermes spp., Isoptera: Termitidae) in northern guinea savanna West Africa

Termite mound densities in typical guinea savanna, Detarium, and grassland (boval) habitats in northern guinea savanna were determined by random quadratting of 2-3 sites in each habitat (100, 10x10 m quadrats per habitat). Dominant species in guinea savanna were T. geminatus (46 mounds ha^-1) and T. oeconomus (21 mounds ha^-1), in Detarium T. geminatus (59 mounds ha^-1) and C. curtatus (45 mounds ha^-1) and in boval C. curtatus (72 mounds ha^-1) and T. geminatus (22 mounds ha^-1). Only C. curtatus densities and total densities differed significantly between sites within habitats, but all species differed significantly in abundance between habitats. The composition of each community was related to general environment but no particular environmental variable was shown to be a major determinant of termite distribution. Evidence for the limitation of termite populations was obtained from indirect evidence of competition between colonies in Detarium, and by experimental manipulation of fire regimes in the typical guinea savanna habitat. Harvester termites increased four-five fold over two years in fire-protected plots as a result of increased food supplies. Total termite densities in the fire-protected community equilibrated to the new population density (100 mounds ha^-1) after only two-three years.

MAJOR GENETIC DIFFERENCES BETWEEN CROWN-OF-THORNS STARFISH (ACANTHASTER PLANCI) POPULATIONS IN THE INDIAN AND PACIFIC OCEANS

Spatial variation in allelic frequencies at nine allozyme loci were assayed in 20 populations of the crown-of-thorns starfish, Acanthaster planci, collected throughout the Pacific and Indian Oceans. These data were analyzed together with published data, for the same loci, from an additional 19 populations, giving a total sample size of approximately 1800 individuals. There was a marked discontinuity between the Indian and Pacific Ocean populations, but those off Western Australia and from the Southeast Asian region had a strong Pacific affinity. The genetic groups were congruent with the distributions of two color morph groups: gray-green to red-brown forms in the Pacific and a blue to pale red form in the Indian Ocean. These patterns of genetic structure are similar to those described for the starfish Linckia laevigata, which has similar life-history characteristics. Vicariant events may have influenced some populations within the Pacific, but the allozyme data cannot resolve the effects of these events clearly. Patterns of variation within regions were consistent with isolation by distance, but, at larger scales, were obscured by regional vicariance and some outliers, particularly by apparently high levels of gene flow between Japan and the Great Barrier Reef, Australia. Apparent gene flow between population pairs was not closely related to present-day ocean currents. The results demonstrate a strong influence of allopatric separation on genetic divergence at large geographic scales, but also show evidence of slow rates of change in gene frequencies consistent with the large population sizes of this species. Low levels of divergence between groups demonstrate the genetic structure is recent (Pleistocene) and are likely responses to changes in climate and sea level.

Quantifying the genetic parameters of feed efficiency in juvenile Nile tilapia Oreochromis niloticus

BACKGROUND

Improving feed efficiency in fish is crucial at the economic, social and environmental levels with respect to developing a more sustainable aquaculture. The important contribution of genetic improvement to achieve this goal has been hampered by the lack of accurate basic information on the genetic parameters of feed efficiency in fish. We used video assessment of feed intake on individual fish reared in groups to estimate the genetic parameters of six growth traits, feed intake, feed conversion ratio (FCR) and residual feed intake in 40 pedigreed families of the GIFT strain of Nile tilapia, Oreochromis niloticus. Feed intake and growth were measured on juvenile fish (22.4 g mean body weight) during 13 consecutive meals, representing 7 days of measurements. We used these data to estimate the FCR response to different selection criteria to assess the potential of genetics as a means of increasing FCR in tilapia.

RESULTS

Our results demonstrate genetic control for FCR in tilapia, with a heritability estimate of 0.32 ± 0.11. Response to selection estimates showed FCR could be efficiently improved by selective breeding. Due to low genetic correlations, selection for growth traits would not improve FCR. However, weight loss at fasting has a high genetic correlation with FCR (0.80 ± 0.25) and a moderate heritability (0.23), and could be an easy to measure and efficient criterion to improve FCR by selective breeding in tilapia.

CONCLUSION

At this age, FCR is genetically determined in Nile tilapia. A selective breeding program could be possible and could help enabling the development of a more sustainable aquaculture production.

Small-Scale Dispersion of Eggs and Sperm of the Crown-of-Thorns Starfish (Acanthaster planci) in a Shallow Coral Reef Habitat

The dispersal of eggs and sperm of crown-of-thorns starfish, Acanthaster planci (L.), was measured in the field using an array of collectors up to 10 m downstream of a spawning starfish. Hydrodynamic measurements, gamete dispersal numerical models, and the gamete cloud dispersal measurements for the first time quantified the relationship between hydrodynamic conditions and the dispersion of eggs and sperm in the field. In general, gamete concentrations fell rapidly and logarithmically with distance from the spawning starfish; egg concentrations at 3 m were 1% of those near the starfish. Simplified dispersal models showed a good correspondence with these field data, and confirmed the observation that eggs rose higher in the water column and spread more laterally at low current speeds over the short spatial scales being considered. Fertilization rates, predicted from laboratory measurements of fertilization success and the gamete concentrations measured in the field, were estimated to be 90-100% within 1 m and 70-100% at 10 m. These results are explained by high success rates of fertilization (fertilizing capacity) at the measured dilutions, and were similar to fertilization rates previously measured by others for crown-of-thorns starfish in the field. Although the eggs were observed to spread upwards into the water column due to turbulence, laboratory measurements of sinking rates showed eggs to be very slightly negatively buoyant (median fall velocity of 0.072 mm·s-1), whereas sperm were neutrally buoyant. A significant fraction of eggs also entered the seabed near the starfish; the proportion decreased with increasing current strength. This process may provide a mechanism for enhanced fertilization of these gametes and/or a mechanism for self-recruitment to a given reef population.

Development of Diversity Arrays Technology markers as a tool for rapid genomic assessment in Nile tilapia, Oreochromis niloticus

The development of genomic markers is described for Nile tilapia, Oreochromis niloticus, using the Diversity Arrays Technology (DArT) genotype-by-sequencing platform. A total of 13 215 single nucleotide polymorphism (SNP) markers and 12 490 silicoDArT (dominant) markers were identified from broodstock of two selective breeding programs [Genetically Improved Farmed Tilapia (GIFT) strain from Malaysia and the Abbassa strain from Egypt]. Over 10 000 SNPs were polymorphic in either strain, and 2985 and 3087 showed strain-specific polymorphisms for the GIFT and Abbassa strains respectively. We demonstrate the potential utility of these markers for rapid genomic screening and use in breeding programs.

Sex determination in the GIFT strain of tilapia is controlled by a locus in linkage group 23

Background

Tilapias (Family Cichlidae) are the second most important group of aquaculture species in the world. They have been the subject of much research on sex determination due to problems caused by early maturation in culture and their complex sex-determining systems. Different sex-determining loci (linkage group 1, 20 and 23) have been detected in various tilapia stocks. The ‘genetically improved farmed tilapia’ (GIFT) stock, founded from multiple Nile tilapia (Oreochromis niloticus) populations, with some likely to have been introgressed with O. mossambicus, is a key resource for tilapia aquaculture. The sex-determining mechanism in the GIFT stock was unknown, but potentially complicated due to its multiple origins.

Results

A bulk segregant analysis (BSA) version of double-digest restriction-site associated DNA sequencing (BSA-ddRADseq) was developed and used to detect and position sex-linked single nucleotide polymorphism (SNP) markers in 19 families from the GIFT strain breeding nucleus and two Stirling families as controls (a single XY locus had been previously mapped to LG1 in the latter). About 1500 SNPs per family were detected across the genome. Phenotypic sex in Stirling families showed strong association with LG1, whereas only SNPs located in LG23 showed clear association with sex in the majority of the GIFT families. No other genomic regions linked to sex determination were apparent. This region was validated using a series of LG23-specific DNA markers (five SNPs with highest association to sex from this study, the LG23 sex-associated microsatellite UNH898 and ARO172, and the recently isolated amhy marker for individual fish (n = 284).

Conclusions

Perhaps surprisingly given its multiple origins, sex determination in the GIFT strain breeding nucleus was associated only with a locus in LG23. BSA-ddRADseq allowed cost-effective analysis of multiple families, strengthening this conclusion. This technique has potential to be applied to other complex traits. The sex-linked SNP markers identified will be useful for potential marker-assisted selection (MAS) to control sex-ratio in GIFT tilapia to suppress unwanted reproduction during growout.

Breeding for robustness: investigating the genotype-by-environment interaction and micro-environmental sensitivity of Genetically Improved Farmed Tilapia (Oreochromis niloticus)

Robustness has become a highly desirable breeding goal in the globalized agricultural market. Both genotype‐by‐environment interaction (G × E) and micro‐environmental sensitivity are important robustness components of aquaculture production, in which breeding stock is often disseminated to different environments. The objectives of this study were (i) to quantify the degree of G × E by assessing the growth performance of Genetically Improved Farmed Tilapia (GIFT) across three countries (Malaysia, India and China) and (ii) to quantify the genetic heterogeneity of environmental variance for body weight at harvest (BW) in GIFT as a measure of micro‐environmental sensitivity. Selection for BW was carried out for 13 generations in Malaysia. Subsets of 60 full‐sib families from Malaysia were sent to China and India after five and nine generations respectively. First, a multi‐trait animal model was used to analyse the BW in different countries as different traits. The results indicate a strong G × E. Second, a genetically structured environmental variance model, implemented using Bayesian inference, was used to analyse micro‐environmental sensitivity of BW in each country. The analysis revealed the presence of genetic heterogeneity of both BW and its environmental variance in all environments. The presence of genetic variation in residual variance of BW implies that the residual variance can be modified by selection. Incorporating both G × E and micro‐environmental sensitivity information may help in selecting robust genotypes with high performance across environments and resilience to environmental fluctuations.

Genomic Selection and Genome-wide Association Study for Feed-Efficiency Traits in a Farmed Nile Tilapia (Oreochromis niloticus) Population

Nile tilapia is a key aquaculture species with one of the highest production volumes globally. Genetic improvement of feed efficiency via selective breeding is an important goal, and genomic selection may expedite this process. The aims of this study were to 1) dissect the genetic architecture of feed-efficiency traits in a Nile tilapia breeding population, 2) map the genomic regions associated with these traits and identify candidate genes, 3) evaluate the accuracy of breeding value prediction using genomic data, and 4) assess the impact of the genetic marker density on genomic prediction accuracies. Using an experimental video recording trial, feed conversion ratio (FCR), body weight gain (BWG), residual feed intake (RFI) and feed intake (FI) traits were recorded in 40 full-sibling families from the GIFT (Genetically Improved Farmed Tilapia) Nile tilapia breeding population. Fish were genotyped with a ThermoFisher Axiom 65 K Nile tilapia SNP array. Significant heritabilities, ranging from 0.12 to 0.22, were estimated for all the assessed traits using the genomic relationship matrix. A negative but favourable genetic correlation was found between BWG and the feed-efficiency related traits; -0.60 and -0.63 for FCR and RFI, respectively. While the genome-wide association analyses suggested a polygenic genetic architecture for all the measured traits, there were significant QTL identified for BWG and FI on chromosomes seven and five respectively. Candidate genes previously found to be associated with feed-efficiency traits were located in these QTL regions, including ntrk3a, ghrh and eif4e3. The accuracy of breeding value prediction using the genomic data was up to 34% higher than using pedigree records. A SNP density of approximately 5,000 SNPs was sufficient to achieve similar prediction accuracy as the full genotype data set. Our results highlight the potential of genomic selection to improve feed efficiency traits in Nile tilapia breeding programmes.

MHC-mediated spatial distribution in brown trout (Salmo trutta) fry

Major histocompatibility complex (MHC) class I-linked microsatellite data and parental assignment data for a group of wild brown trout (Salmo trutta L.) provide evidence of closer spatial aggregation among fry sharing greater numbers of MHC class I alleles under natural conditions. This result confirms predictions from laboratory experiments demonstrating a hierarchical preference for association of fry sharing MHC alleles. Full-siblings emerge from the same nest (redd), and a passive kin association pattern arising from limited dispersal from the nest (redd effect) would predict that all such pairs would have a similar distribution. However, this study demonstrates a strong, significant trend for reduced distance between pairs of full-sibling fry sharing more MHC class I alleles reflecting their closer aggregation (no alleles shared, 311.5 ± (s.e.)21.03 m; one allele shared, 222.2 ± 14.49 m; two alleles shared, 124.9 ± 23.88 m; P<0.0001). A significant trend for closer aggregation among fry sharing more MHC class I alleles was also observed in fry pairs, which were known to have different mothers and were otherwise unrelated (ML-r = 0) (no alleles: 457.6 ± 3.58 m; one allele (422.4 ± 3.86 m); two alleles (381.7 ± 10.72 m); P<0.0001). These pairs are expected to have emerged from different redds and a passive association would then be unlikely. These data suggest that sharing MHC class I alleles has a role in maintaining kin association among full-siblings after emergence. This study demonstrates a pattern consistent with MHC-mediated kin association in the wild for the first time.

Single Nucleotide Polymorphisms (SNPs) Reveal Sibship Among Founders of a Bangladeshi Rohu (Labeo rohita) Breeding Population

Rohu (Labeo rohita) is a significant freshwater aquaculture species with approximately 1.8 Mt produced annually. Fin clips obtained from the founders of a newly established Bangladesh-based breeding population (∼140 fish from each of the Halda, Jamuna, and Padma rivers) were used to identify 9157 SNPs and 14 411 silicoDArT markers using the Diversity Arrays Technology (DArT) genotyping-by-sequencing platform known as DArTseq. After quality control, 1985 SNPs were retained and used to examine population structure within and among river systems. Examination of genomic relationships revealed evidence of full- and half-sibling relationships among founders. Accordingly, sibship and dummy parents were assigned within each river population using a maximum likelihood approach with COLONY software. Founders that had no dummy parents in common were then identified for population genetic analyses. Only 40 unique dummy parents and 17 founders with no common dummy parents were identified from the Halda river, compared with 206 (96) from the Jamuna and 184 (83) from the Padma. Overall pairwise FST estimates among rivers were low (<0.005) and the optimum number of clusters using unsupervised K-means clustering was one, indicating little genetic divergence among the river populations in our SNPs. These results suggest that observed sibship among founders should be accounted for in future pedigree-based analyses and it cannot be assumed that fertilized spawn collections are representative samples of river populations.

Sibship assignment to the founders of a Bangladeshi Catla catla breeding population

Catla catla (Hamilton) fertilised spawn was collected from the Halda, Jamuna and Padma rivers in Bangladesh from which approximately 900 individuals were retained as 'candidate founders' of a breeding population. These fish were fin-clipped and genotyped using the DArTseq platform to obtain, 3048 single nucleotide polymorphisms (SNPs) and 4726 silicoDArT markers. Using SNP data, individuals that shared no putative parents were identified using the program COLONY, i.e. 140, 47 and 23 from the Halda, Jamuna and Padma rivers, respectively. Allele frequencies from these individuals were considered as representative of those of the river populations, and genomic relationship matrices were generated. Then, half-sibling and full-sibling relationships between individuals were assigned manually based on the genomic relationship matrices. Many putative half-sibling and full-sibling relationships were found between individuals from the Halda and Jamuna rivers, which suggests that catla sampled from rivers as spawn are not necessarily representative of river populations. This has implications for the interpretation of past population genetics studies, the sampling strategies to be adopted in future studies and the management of broodstock sourced as river spawn in commercial hatcheries. Using data from individuals that shared no putative parents, overall multi-locus pairwise estimates of Wright's fixation index (FST) were low (≤ 0.013) and the optimum number of clusters using unsupervised K-means clustering was equal to 1, which indicates little genetic divergence among the SNPs included in our study within and among river populations.