Estimating heritability using family-pooled phenotypic and genotypic data: a simulation study applied to aquaculture

Estimating heritability based on individual phenotypic and genotypic measurements can be expensive and labour-intensive in commercial aquaculture breeding. Here, the feasibility of estimating heritability using within-family means of phenotypes and allelic frequencies was investigated. Different numbers of full-sib families and family sizes across ten generations with phenotypic and genotypic information on 10 K SNPs were analysed in ten replicates. Three scenarios, representing differing numbers of pools per family (one, two and five) were considered. The results showed that using one pool per family did not reliably estimate the heritability of family means. Using simulation parameters appropriate for aquaculture, at least 200 families of 60 progeny per family divided equally in two pools per family was required to estimate the heritability of family means effectively. Although application of five pools generated more within- and between- family relationships, it reduced the number of individuals per pool and increased within-family residual variation, hence, decreased the heritability of family means. Moreover, increasing the size of pools resulted in increasing the heritability of family means towards one. In addition, heritability of family mean estimates were higher than family heritabilities obtained from Falconer’s formula due to lower intraclass correlation estimate compared to the coefficient of relationship

Quantitative trait loci mapping in dairy cattle: review and meta-analysis

From an extensive review of public domain information on dairy cattle quantitative trait loci (QTL), we have prepared a draft online QTL map for dairy production traits. Most publications (45 out of 55 reviewed) reported QTL for the major milk production traits (milk, fat and protein yield, and fat and protein concentration (%)) and somatic cell score. Relatively few QTL studies have been reported for more complex traits such as mastitis, fertility and health. The collated QTL map shows some chromosomal regions with a high density of QTL, as well as a substantial number of QTL at single chromosomal locations. To extract the most information from these published records, a meta-analysis was conducted to obtain consensus on QTL location and allelic substitution effect of these QTL. This required modification and development of statistical methodologies. The meta-analysis indicated a number of consensus regions, the most striking being two distinct regions affecting milk yield on chromosome 6 at 49 cM and 87 cM explaining 4.2 and 3.6 percent of the genetic variance of milk yield, respectively. The first of these regions (near marker BM143) affects five separate milk production traits (protein yield, protein percent, fat yield, fat percent, as well as milk yield)

Prawn Morphometrics and Weight Estimation from Images using Deep Learning for Landmark Localization

Accurate morphometric analyses and weight estimation are useful in aquaculture for optimizing feeding, predicting harvest yields, identifying desirable traits for selective breeding, grading processes, and monitoring the health status of production animals. However, the collection of phenotypic data through traditional manual approaches at industrial scales and in real-time is time-consuming, labour-intensive, and prone to errors. Digital imaging of individuals and subsequent training of prediction models using Deep Learning (DL) has the potential to rapidly and accurately acquire phenotypic data from aquaculture species. In this study, we applied a novel DL approach to automate morphometric analysis and weight estimation using the black tiger prawn (Penaeus monodon) as a model crustacean. The DL approach comprises two main components: a feature extraction module that efficiently combines low-level and high-level features using the Kronecker product operation; followed by a landmark localization module that then uses these features to predict the coordinates of key morphological points (landmarks) on the prawn body. Once these landmarks were extracted, weight was estimated using a weight regression module based on the extracted landmarks using a fully connected network. For morphometric analyses, we utilized the detected landmarks to derive five important prawn traits. Principal Component Analysis (PCA) was also used to identify landmark-derived distances, which were found to be highly correlated with shape features such as body length, and width. We evaluated our approach on a large dataset of 8164 images of the Black tiger prawn (Penaeus monodon) collected from Australian farms. Our experimental results demonstrate that the novel DL approach outperforms existing DL methods in terms of accuracy, robustness, and efficiency

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

Linkage disequilibrium on chromosome 6 in Australian Holstein-Friesian cattle

We analysed linkage disequilibrium (LD) in Australian Holstein-Friesian cattle by genotyping a sample of 45 bulls for 15 closely-spaced microsatellites on two regions of BTA6 reported to carry important QTL for dairy traits. The order and distance of markers were based on the USDA-MARC linkage map. Frequencies of haplotypes were estimated using the E-M approach and a more computationally-intensive Bayesian approach as implemented in PHASE. LD was then estimated using the Hedrick multiallelic extension of Lewontin normalised coefficient D'. Estimates of D' from the two approaches were in close agreement (r = 0.91). The mean estimates of D' for marker pairs with an inter-marker distance of less than 5 cM (n = 13) are 0.57 and 0.51, and for distances more than 20 cM (n = 44) are 0.29 and 0.17, estimated from the E-M and Bayesian approaches, respectively. The Malecot model was fitted for the exponential decline of LD with map distance between markers. The swept radii (the distance at which LD has declined to 1/e (~37%) of its initial value) are 11.6 and 13.7 cM for the above two methods, respectively. The Malecot model was also fitted using map distance in Mb from the bovine integrated map (bovine location database, bLDB) in addition to cM from the MARC map. Overall, the results indicate a high level of LD on chromosome 6 in Australian dairy cattle

Genome wide association study for growth in Pakistani dromedary camels using genotyping-by-sequencing

Objective Growth performance and growth-related traits have a crucial role in livestock due to their influence on productivity. This genome-wide association study (GWAS) in Pakistani dromedary camels was conducted to identify single nucleotide polymorphisms (SNPs) associated with growth at specific camel ages, and for selected SNPs, to investigate in detail how their effects change with increasing camel age. This is the first GWAS conducted on dromedary camels in this region. Methods Two Pakistani breeds, Marecha and Lassi, were selected for this study. A genotyping-by-sequencing method was used, and a total of 65,644 SNPs were identified. For GWAS, weight records data with several body weight traits, namely, birthweight, weaning weight, and weights of camels at 1, 2, 4, and 6 years of age were analysed by using model-based growth curve analysis. Age-specific weight data were analysed with a linear mixed model that included fixed effects of SNP genotype as well as sex. Results Based on the q-value method for false discovery control, for Marecha camels, five SNPs at q<0.01 and 96 at q<0.05 were significantly associated with the weight traits considered, while three (q<0.01) and seven (q<0.05) SNP associations were identified for Lassi camels. Several candidate genes harbouring these SNP were discovered. Conclusion These results will help to better understand the genetic architecture of growth including how these genes are expressed at different phases of their life. This will serve to lay the foundations for applied breeding programs of camels by allowing the genetic selection of superior animals

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

Considerations for maintaining family diversity in commercially mass-spawned penaeid shrimp: a case study on Penaeus monodon

Skewed family distributions are common in aquaculture species that are highly fecund, communally (mass) spawned, and/or communally reared. The magnitude of skews pose challenges for maintaining family-specific genetic diversity, as increased resources are required to detect individuals from underrepresented families, or reliably determine relative survival as a measure of family performance. There is limited understanding of family skews or changes in family proportion of communally reared shrimp under commercial rearing conditions and particularly how this may affect genotyping strategies to recover family performance data in breeding programs. In this study, three separate batches of shrimp, Penaeus monodon, were communally spawned and reared, and then sampled as larvae when ponds were stocked at 30 days of culture (DOC) and as juveniles from commercial ponds during harvest at 150 DOC. A total of 199 broodstock contributed to the 5,734 progeny that were genotyped with a custom multiplex single nucleotide polymorphism (SNP) panel, and family assignments were cross-referenced using two parentage assignment methods, CERVUS and COLONY. A total of 121 families were detected, with some families contributing up to 11% of progeny at 30 DOC and up to 18% of progeny at harvest. Significant changes were detected for 20% of families from 30 to 150 DOC, with up to a 9% change in relative contribution. Family skew data was applied in several models to determine the optimal sample size to detect families, along with the ability to detect changes in relative family contribution over time. Results showed that an order of magnitude increase in sampling was required to capture the lowest represented 25% of families, as well as significantly improve the accuracy to determine changes in family proportion from 30 to 150 DOC. Practical measures may be implemented at the hatchery to reduce family skews; a cost-effective measure may be to address the initial magnitude differences in viable progeny produced among families, by pooling equal quantities of hatched larvae from each family. This study demonstrates the relationships between skews in families under commercial conditions, the ability to accurately detect families, and the balance of sampling effort and genotyping cost in highly fecund species such as shrimp

A comparative integrated gene-based linkage and locus ordering by linkage disequilibrium map for the Pacific white shrimp, Litopenaeus vannamei

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 incorporatedm 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

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