Transcriptomics approaches to identify genes and markers for production traits in giant freshwater prawn Macrobrachium rosenbergii

This study used next generation sequencing technologies to investigate growth in a cultured crustacean. The objective was to identify and characterise specific gene loci that contribute important phenotypic variation to growth as well as to develop a large set of SNP markers in candidate genes for assessing correlations between specific mutations and individual growth performance. The genomic dataset generated provides a fundamental resource for application in future crustacean stock improvement programs. Ultimately, the data can be applied to development of culture lines with improved growth performance.

Genetic improvement of giant freshwater prawn in Vietnam

The giant freshwater prawn (Macrobrachium rosenbergii) or GFP is one of the most important freshwater crustacean species in the inland aquaculture sector of many tropical and subtropical countries. Since the 1990’s, there has been rapid global expansion of freshwater prawn farming, especially in Asian countries, with an average annual rate of increase of 48% between 1999 and 2001 (New, 2005). In Vietnam, GFP is cultured in a variety of culture systems, typically in integrated or rotational rice-prawn culture (Phuong et al., 2006) and has become one of the most common farmed aquatic species in the country, due to its ability to grow rapidly and to attract high market price and high demand. Despite potential for expanded production, sustainability of freshwater prawn farming in the region is currently threatened by low production efficiency and vulnerability of farmed stocks to disease. Commercial large scale and small scale GFP farms in Vietnam have experienced relatively low stock productivity, large size and weight variation, a low proportion of edible meat (large head to body ratio), scarcity of good quality seed stock. The current situation highlights the need for a systematic stock improvement program for GFP in Vietnam aimed at improving economically important traits in this species. This study reports on the breeding program for fast growth employing combined (between and within) family selection in giant freshwater prawn in Vietnam. The base population was synthesized using a complete diallel cross including 9 crosses from two local stocks (DN and MK strains) and a third exotic stock (Malaysian strain - MY). In the next three selection generations, matings were conducted between genetically unrelated brood stock to produce full-sib and (paternal) half-sib families. All families were produced and reared separately until juveniles in each family were tagged as a batch using visible implant elastomer (VIE) at a body size of approximately 2 g. After tags were verified, 60 to 120 juveniles chosen randomly from each family were released into two common earthen ponds of 3,500 m2 pond for a grow-out period of 16 to 18 weeks. Selection applied at harvest on body weight was a combined (between and within) family selection approach. 81, 89, 96 and 114 families were produced for the Selection line in the F0, F1, F2 and F3 generations, respectively. In addition to the Selection line, 17 to 42 families were produced for the Control group in each generation. Results reported here are based on a data set consisting of 18,387 body and 1,730 carcass records, as well as full pedigree information collected over four generations. Variance and covariance components were estimated by restricted maximum likelihood fitting a multi-trait animal model. Experiments assessed performance of VIE tags in juvenile GFP of different size classes and individuals tagged with different numbers of tags showed that juvenile GFP at 2 g were of suitable size for VIE tags with no negative effects evident on growth or survival. Tag retention rates were above 97.8% and tag readability rates were 100% with a correct assignment rate of 95% through to mature animal size of up to 170 g. Across generations, estimates of heritability for body traits (body weight, body length, cephalothorax length, abdominal length, cephalothorax width and abdominal width) and carcass weight traits (abdominal weight, skeleton-off weight and telson-off weight) were moderate and ranged from 0.14 to 0.19 and 0.17 to 0.21, respectively. Body trait heritabilities estimated for females were significantly higher than for males whereas carcass weight trait heritabilities estimated for females and males were not significantly different (P > 0.05). Maternal and common environmental effects for body traits accounted for 4 to 5% of the total variance and were greater in females (7 to 10%) than in males (4 to 5%). Genetic correlations among body traits were generally high in both sexes. Genetic correlations between body and carcass weight traits were also high in the mixed sexes. Average selection response (% per generation) for body weight (transformed to square root) estimated as the difference between the Selection and the Control group was 7.4% calculated from least squares means (LSMs), 7.0% from estimated breeding values (EBVs) and 4.4% calculated from EBVs between two consecutive generations. Favourable correlated selection responses (estimated from LSMs) were detected for other body traits (12.1%, 14.5%, 10.4%, 15.5% and 13.3% for body length, cephalothorax length, abdominal length, cephalothorax width and abdominal width, respectively) over three selection generations. Data in the second selection generation showed positive correlated responses for carcass weight traits (8.8%, 8.6% and 8.8% for abdominal weight, skeleton-off weight and telson-off weight, respectively). Data in the third selection generation showed that heritability for body traits were moderate and ranged from 0.06 to 0.11 and 0.11 to 0.22 at weeks 10 and 18, respectively. Body trait heritabilities estimated at week 10 were not significantly lower than at week 18. Genetic correlations between body traits within age and genetic correlations for body traits between ages were generally high. Overall our results suggest that growth rate responds well to the application of family selection and carcass weight traits can also be improved in parallel, using this approach. Moreover, selection for high growth rate in GFP can be undertaken successfully before full market size has been reached. The outcome of this study was production of an improved culture strain of GFP for the Vietnamese culture industry that will be trialed in real farm production environments to confirm the genetic gains identified in the experimental stock improvement program.

Genetic diversity of cultured and wild populations of the giant freshwater prawn Macrobrachium rosenbergii (de Man, 1879) based on microsatellite analysis

Freshwater prawn (Macrobrachium rosenbergii) culture in the Western Hemisphere is primarily, if not entirely, derived from 36 individual prawns originally introduced to Hawaii from Malaysia in 1965 and 1966. Little information is available regarding genetic variation within and among cultured prawn stocks worldwide. The goal of the current study was to characterize genetic diversity in various prawn populations with emphasis on those cultured in North America. Five microsatellite loci were screened to estimate genetic diversity in two wild (Myanmar and India-wild) and seven cultured (Hawaii-1, Hawaii-2, India-cultured, Israel, Kentucky, Mississippi and Texas) populations. Average allelic richness ranged from 3.96 (Israel) to 20.45 (Myanmar). Average expected heterozygosity ranged from 0.580 (Israel) to 0.935 (Myanmar). Many of the cultured populations exhibited reduced genetic diversity when compared with the Myanmar and the India-cultured populations. Significant deficiency in heterozygotes was detected in the India-cultured, Mississippi and Kentucky populations (overall Fis estimated of 0.053, 0.067 and 0.108 respectively) reflecting moderate levels of inbreeding. Overall estimate of fixation index (Fst = 0.1569) revealed moderately high levels of differentiation among the populations. Outcome of this study provide a baseline assessment of genetic diversity in some available strains that will be useful for the development of breeding programmes.

Quantitative genetic parameter estimates for body and carcass traits in a cultured stock of giant freshwater prawn (Macrobrachium rosenbergii) selected for harvest weight in Vietnam

We estimated the heritability and correlations between body and carcass weight traits in a cultured stock of giant freshwater prawn (GFP) (Macrobrachium rosenbergii) selected for harvest body weight in Vietnam. The data set consisted of 18,387 body and 1,730 carcass records, as well as full pedigree information collected over four generations. Variance and covariance components were estimated by restricted maximum likelihood fitting a multi-trait animal model. Across generations, estimates of heritability for body and carcass weight traits were moderate and ranged from 0.14 to 0.19 and 0.17 to 0.21, respectively. Body trait heritabilities estimated for females were significantly higher than for males whereas carcass weight trait heritabilities estimated for females and males were not significantly different (P>. 0.05). Maternal effects for body traits accounted for 4 to 5% of the total variance and were greater in females than in males. Genetic correlations among body traits were generally high in the mixed sexes. Genetic correlations between body and carcass weight traits were also high. Although some issues remain regarding the best statistical model to be fitted to GFP data, our results suggest that selection for high harvest body weight based on breeding values estimated by fitting an animal model to the data can significantly improve mean body and carcass weight in GFP.

Genetic response to combined family selection for improved mean harvest weight in giant freshwater prawn (Macrobrachium rosenbergii) in Vietnam

We estimated genetic changes in body and carcass weight traits in a giant freshwater prawn (GFP) (Macrobrachium rosenbergii) population selected for increased body weight at harvest in Vietnam. The data set consisted of 18,387 individual body and 1730 carcass weight records, as well as full pedigree information collected over four generations. Average selection response (per generation) in body weight at harvest (transformed to square root) estimated as the difference between the Selection line and the Control group was 7.4% calculated from least squares mean (LSMs), 7.0% from estimated breeding values (EBVs) and 4.4% calculated from EBVs between two consecutive generations. Favorable correlated selection responses (estimated from LSMs) were found for other body traits including: total length, cephalothorax length, abdominal length, cephalothorax width, and abdominal width (12.1%, 14.5%, 10.4%, 15.5% and 13.3% over three selection generations, respectively). Data in the second generation of selection showed positive correlated responses for carcass weight traits including: abdominal weight, exoskeleton-off weight, and telson-off weight of 8.8%, 8.6% and 8.8%, respectively. We conclude that body weight at harvest responded well to the application of combined (between and within) family selection and correlated responses in carcass weight traits were favorable.

Impacts of climatic factors on evolution of molecular diversity and the natural distribution of wild stocks of the giant freshwater prawn (Macrobrachium rosenbergii)

Natural distributions of most freshwater taxa are restricted geographically, a pattern that reflects dispersal limitation. Macrobrachium rosenbergii is unusual because it occurs naturally in rivers from near Pakistan in the west, across India and Bangladesh to the Malay Peninsula, and across the Sunda Shelf and Indonesian archipelago to western Java. Individuals cannot tolerate full marine conditions, so dispersal between river drainage basins must occur at limited geographical scales when ecological or climatic factors are favorable. We examined molecular diversity in wild populations of M. rosenbergii across its complete natural range to document patterns of diversity and to relate them to factors that have driven evolution of diversity in this species. We found 3 clades in the mitochondrial deoxyribonucleic acid (mtDNA) data set that corresponded geographically with eastern, central, and western sets of haplotypes that last shared a common ancestor 1 × 106 y ago. The eastern clade was closest to the common ancestor of all 3 clades and to the common ancestor with its congener, Macrobrachium spinipes, distributed east of Huxley's Line. Macrobrachium rosenbergii could have evolved in the western Indonesian archipelago and spread westward during the early to mid-Pleistocene to India and Sri Lanka. Additional groups identified in the nuclear DNA data set in the central and western clades probably indicate secondary contact via dispersal between regions and modern introductions that have mixed nuclear and mtDNA genes. Pleistocene sea-level fluctuations can explain dispersal across the Indonesian archipelago and parts of mainland southeastern Asia via changing river drainage connections in shallow seas on wide continental shelves. At the western end of the modern distribution where continental shelves are smaller, intermittent freshwater plumes from large rivers probably permitted larval dispersal across inshore areas of lowered salinity.

Quantitative genetic parameters for body traits at different ages in a cultured stock of giant freshwater prawn (Macrobrachium rosenbergii) selected for fast growth

The aim of the current study was to estimate heritabilities and correlations for body traits at different ages (Weeks 10 and 18 after stocking) in a giant freshwater prawn (Macrobrachium rosenbergii) population selected for fast growth rate in Vietnam. The dataset consisted of 4650 body records (2432 and 2218 records collected at Weeks 10 and 18, respectively) in the full pedigree comprising a total of 18 387 records. Variance and covariance components were estimated using restricted maximum likelihood fitting a multi-trait animal model. Estimates of heritability for body traits (bodyweight, body length, cephalothorax length, abdominal length, cephalothorax width and abdominal width) were moderate and ranged from 0.06 to 0.11 and from 0.11 to 0.22 at Weeks 10 and 18, respectively. Body-trait heritabilities estimated at Week 10 were not significantly lower than at Week 18. Genetic correlations between body traits within age and genetic correlations for body traits between ages were generally high. Our results suggested that selection for high growth rate in GFP can be undertaken successfully before full market size has been reached.

Impacts of vessel capacity reduction programmes on efficiency in fisheries : the case of Australia’s multispecies northern prawn fishery

Capacity reduction programmes, in the form of buybacks or decommissioning, have had relatively widespread application in fisheries in the US, Europe and Australia. A common criticism of such programmes is that they remove the least efficient vessels first, resulting in an increase in average efficiency of the remaining fleet, which tends to increase the effective fishing power of the remaining fleet. In this paper, the effects of a buyback programme on average technical efficiency in Australia’s Northern Prawn Fishery are examined using a multi-output production function approach with an explicit inefficiency model. As expected, the results indicate that average efficiency of the remaining vessels was generally greater than that of the removed vessels. Further, there was some evidence of an increase in average scale efficiency in the fleet as the remaining vessels were closer, on average, to the optimal scale. Key factors affecting technical efficiency included company structure and the number of vessels fishing. In regard to fleet size, our model suggests positive externalities associated with more boats fishing at any point in time (due to information sharing and reduced search costs), but also negative externalities due to crowding, with the latter effect dominating the former. Hence, the buyback resulted in a net increase in the individual efficiency of the remaining vessels due to reduced crowding, as well as raising average efficiency through removal of less efficient vessels.

Phylogeography of two freshwater prawn species from far-northern Queensland

The distribution of a freshwater species is often dependent on its ability to disperse within the riverine system. Species with high dispersal abilities tend to be widespread, whereas those with restricted dispersal tend to be geographically restricted and are usually given higher conservation priority. Population structure was compared between a widespread freshwater prawn species, Macrobrachium australiense, and a narrow-range endemic freshwater prawn, Macrobrachium koombooloomba. The distribution of M. australiense and M. koombooloomba did not overlap, although suggested historical river-boundary rearrangements indicate that there has been the potential for dispersal into neighbouring catchments. A fragment of the mtDNA CO1 gene was analysed and a Mantel test revealed a significant isolation by distance effect for both species. Significant overall FST values confirmed that both species exhibited low levels of dispersal, a prediction for populations inhabiting a fragmented upland environment. The level of structure in M. australiense is surprising for a widely distributed species. Not all M. australiense populations conformed to the stream-hierarchy model, with results being best explained by historical river realignment or cross-catchment dispersal. The fact that both species show limited dispersal highlights the importance of conservation in highland areas for both endemic and widely spread species.

Identification of a putative cellulase gene in the giant freshwater prawn, Macrobrachium rosenbergii (De Man, 1879)

Nutrition plays an important role in the development of all organisms and in particular that of farmed aquatic species where costs associated with feed can often exceed 60% of total production costs. Crustacean species in addition, have the added metabolic requirement for regular moulting to allow normal growth and this requires large amounts of energy in the form of sugars (glucose). The current study explored the capacity of the giant freshwater prawn to produce endogenous cellulose-degrading enzymes capable of extracting nutrients (simple sugars) from plant sources in formulated feeds used in the prawn aquaculture industry. We identified a putative cellulase cDNA fragment in the target organism of 1576 base pairs in length of non-microbial origin that after protein modelling exhibited a TM-score of 0.916 with a described cellulase reported from another crustacean species. The functional role of cellulase enzymes is to hydrolyse cellulose to glucose and the fragment identified in GFP was highly expressed in the hepatopancreas, the site of primary food digestion and absorption in crustaceans. Hepatopancreatic tissue from Macrobrachium rosenbergii also showed active digestion of cellulose to glucose following an endoglucanase assay. Cellulase gene(s) are present in the genomes of many invertebrate taxa and play an active role in the conversion of cellulose to available energy. Identification and characterization of endogenous cellulase gene(s) in giant freshwater prawn can assist development of the culture industry because the findings confirm that potentially greater levels of low-cost plant-material could be included in artificial formulated diets in the future without necessarily compromising individual growth performance. Ultimately, this development may contribute to more efficient, cost-effective production systems for freshwater prawn culture stocks that meet the animal's basic nutritional requirements and that also support good individual growth rates.

Genetic evidence for natural hybridization between Red Snow Crab (Chionoecetes japonicus) and Snow Crab (Chionoecetes opilio) in Korea

Both red snow crab (Chionoecetes japonicus Rathbun, 1932) and snow crab (Chionoecetes opilio Fabricius, 1788) are commercially important species in Korea. The geographical ranges of the two species overlap in the East Sea, where both species are fished commercially. Morphological identification of the two species and putative hybrids can be difficult because of their overlapping morphological characteristics. The presence of putative hybrids can affect the total allowable catch (TAC) of C. japonicus and C. opilio, and causes problems managing C. japonicus and C. opilio wild resources. To date, however, no natural hybridization has been reported between C. japonicus and C. opilio, despite their overlapping distributions along the coast of the East Sea. In this study, the internal transcribed spacer (ITS) region of major ribosomal RNA genes from the nuclear genome and the cytochrome oxidase I (CO I) gene from the mitochondrial genome were sequenced to determine whether natural hybridization occurs between the two species. Our results revealed that all putative hybrids identified using morphological traits had two distinct types of ITS sequences corresponding to those of both parental species. Mitochondrial CO I gene sequencing showed that all putative hybrids had sequences identical to C. japonicus. A genotyping assay based on single nucleotide polymorphisms in the ITS1 region and the CO I gene produced the most efficient and accurate identification of all hybrid individuals. Molecular data clearly demonstrate that natural hybridization does occur between C. japonicus and C. opilio, but only with C. japonicus as the maternal parent.

A candidate gene association study for growth performance in an improved Giant Freshwater Prawn (Macrobrachium rosenbergii) culture line

A candidate gene approach using type I single nucleotide polymorphism (SNP) markers can provide an effective method for detecting genes and gene regions that underlie phenotypic variation in adaptively significant traits. In the absence of available genomic data resources, transcriptomes were recently generated in Macrobrachium rosenbergii to identify candidate genes and markers potentially associated with growth. The characterisation of 47 candidate loci by ABI re-sequencing of four cultured and eight wild samples revealed 342 putative SNPs. Among these, 28 SNPs were selected in 23 growth-related candidate genes to genotype in 200 animals selected for improved growth performance in an experimental GFP culture line in Vietnam. The associations between SNP markers and individual growth performance were then examined. For additive and dominant effects, a total of three exonic SNPs in glycogen phosphorylase (additive), heat shock protein 90 (additive and dominant) and peroxidasin (additive), and a total of six intronic SNPs in ankyrin repeats-like protein (additive and dominant), rolling pebbles (dominant), transforming growth factor-β induced precursor (dominant), and UTP-glucose-1-phosphate uridylyltransferase 2 (dominant) genes showed significant associations with the estimated breeding values in the experimental animals (P =0.001−0.031). Individually, they explained 2.6−4.8 % of the genetic variance (R2=0.026−0.048). This is the first large set of SNP markers reported for M. rosenbergii and will be useful for confirmation of associations in other samples or culture lines as well as having applications in marker-assisted selection in future breeding programs.

Modeling meat yield based on measurements of body traits in genetically improved giant freshwater prawn (GFP) Macrobrachium rosenbergii

A single-generation dataset consisting of 1,730 records from a selection program for high growth rate in giant freshwater prawn (GFP, Macrobrachium rosenbergii) was used to derive prediction equations for meat weight and meat yield. Models were based on body traits [body weight, total length and abdominal width (AW)] and carcass measurements (tail weight and exoskeleton-off weight). Lengths and width were adjusted for the systematic effects of selection line, male morphotypes and female reproductive status, and for the covariables of age at slaughter within sex and body weight. Body and meat weights adjusted for the same effects (except body weight) were used to calculate meat yield (expressed as percentage of tail weight/body weight and exoskeleton-off weight/body weight). The edible meat weight and yield in this GFP population ranged from 12 to 15 g and 37 to 45 %, respectively. The simple (Pearson) correlation coefficients between body traits (body weight, total length and AW) and meat weight were moderate to very high and positive (0.75–0.94), but the correlations between body traits and meat yield were negative (−0.47 to −0.74). There were strong linear positive relationships between measurements of body traits and meat weight, whereas relationships of body traits with meat yield were moderate and negative. Step-wise multiple regression analysis showed that the best model to predict meat weight included all body traits, with a coefficient of determination (R 2) of 0.99 and a correlation between observed and predicted values of meat weight of 0.99. The corresponding figures for meat yield were 0.91 and 0.95, respectively. Body weight or length was the best predictor of meat weight, explaining 91–94 % of observed variance when it was fitted alone in the model. By contrast, tail width explained a lower proportion (69–82 %) of total variance in the single trait models. It is concluded that in practical breeding programs, improvement of meat weight can be easily made through indirect selection for body trait combinations. The improvement of meat yield, albeit being more difficult, is possible by genetic means, with 91 % of the variation in the trait explained by the body and carcass traits examined in this study.

Deriving optimal fishing effort for managing Australia's Moreton Bay multispecies trawl fishery with aggregated effort data

Deriving an estimate of optimal fishing effort or even an approximate estimate is very valuable for managing fisheries with multiple target species. The most challenging task associated with this is allocating effort to individual species when only the total effort is recorded. Spatial information on the distribution of each species within a fishery can be used to justify the allocations, but often such information is not available. To determine the long-term overall effort required to achieve maximum sustainable yield (MSY) and maximum economic yield (MEY), we consider three methods for allocating effort: (i) optimal allocation, which optimally allocates effort among target species; (ii) fixed proportions, which chooses proportions based on past catch data; and (iii) economic allocation, which splits effort based on the expected catch value of each species. Determining the overall fishing effort required to achieve these management objectives is a maximizing problem subject to constraints due to economic and social considerations. We illustrated the approaches using a case study of the Moreton Bay Prawn Trawl Fishery in Queensland (Australia). The results were consistent across the three methods. Importantly, our analysis demonstrated the optimal total effort was very sensitive to daily fishing costs-the effort ranged from 9500-11 500 to 6000-7000, 4000 and 2500 boat-days, using daily cost estimates of $0, $500, $750, and $950, respectively. The zero daily cost corresponds to the MSY, while a daily cost of $750 most closely represents the actual present fishing cost. Given the recent debate on which costs should be factored into the analyses for deriving MEY, our findings highlight the importance of including an appropriate cost function for practical management advice. The approaches developed here could be applied to other multispecies fisheries where only aggregated fishing effort data are recorded, as the literature on this type of modelling is sparse.

A retrospective evaluation of sustainable yields for Australia's Northern Prawn Fishery

The fundamental aim in fisheries management is to determine an optimal fishing effort for sustainably harvesting from a replenishable resource. The current management objective of Australia's Northern Prawn Fishery is to maximize the long-term net economic return following Australian government policy, resulting in an average recent catch of tiger prawn species of about 1,250 tons only. However, the maximum sustainable catch stated from different studies is around 3,000-4,700 tons. We also evaluated the net profit assuming that there was no buyback scheme in 2005 and the fishing fleet was kept at 89 vessels since 2005 and concluded that 40% more catch on average (2006-2009) and an additional total profit of A$ 17 million ( excluding crew cost) could have been gained in addition to the many millions of dollars of savings in the buyback scheme. These findings have great implications for future management in Australia and elsewhere because there is a grave concern of overfishing worldwide.