ABSTRACT
Lymphoma is the most common haematological malignancy affecting dogs and has a high incidence in the Bullmastiff breed. The aim of this study was to identify risk loci predisposing this breed to the disease. The average age of lymphoma diagnosis in 55 cases was less than 6 years, similar to the median age of 64 cases from our clinical and pathology databases. When fine-scale population structure was explored using NETVIEW, cases were distributed throughout an extended pedigree. When genotyped cases (n = 49) and dogs from the control group (n = 281) were compared in a genome-wide association analysis of lymphoma risk, the most prominent associated regions were detected on CFA13 and CFA33. The top SNPs in a 5.4 Mb region on CFA13 were significant at a chromosome-wide level, and the region was fine-mapped to ~1.2 Mb (CFA13: 25.2-26.4 Mb; CanFam3.1) with four potential functional candidates, including the MYC proto-oncogene bHLH transcription factor (MYC) and a region syntenic with the human and mouse lncRNA Pvt1 oncogene (PVT1). A 380 Kb associated region at CFA33: 7.7-8.1 Mb contained the coding sequence for SUMO specific peptidase7 (SENP7) and NFK inhibitor zeta (NFKBIZ) genes. These genes have annotations related to cancer, amongst others, and both have functional links to MYC regulation. Genomic signatures identified in lymphoma cases suggest that increased risk contributed by the regions identified by GWAS may complement a complex predisposing genetic background.
ABSTRACT
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 genotypingby-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.
ABSTRACT
The black tiger shrimp, Penaeus monodon, is the second most economically important aquaculture shrimp species in the world, and in Australia it is one of the most commonly farmed shrimp species. Despite its economic significance, very few studies have reported the genetic evaluation of economically important morphological size and shape traits of shrimp grown in commercial grow-out environments. In this study we obtained genetic parameter estimates and evaluated genotype-by-environment interaction (GxE) for nine body morphological traits of shrimp derived from images. The data set contained image and body weight (BW) records of 5,308 shrimp, from 64 sires and 54 dams, reared in eight grow-out ponds for an average of 133 days. From the images, landmark based morphological distances were computed from which novel morphological traits associated with size and shape were derived for genetic evaluation. These traits included body weight (BW), body length (BL), body size (BS), head size (HS), Abdominal size (AS), abdominal percentage (AP), tail tip (TT), front by back ratio (FBR), condition factor (CF) and condition factor length (CFL). We also evaluated G×E interaction effects of these traits for shrimp reared in different ponds. The heritability estimates for growth related morphological and body weight traits were moderately high (BW: h 2 = 0.32 ± 0.05; BL: h 2 = 0.36 ± 0.06; BS: h 2 = 0.32 ± 0.05; HS: h 2 = 0.31 ± 0.05; AS: h 2 = 0.32 ± 0.05; and TT: h 2 = 0.28 ± 0.05) and low for abdominal percentage and body shape traits (AP: h 2 = 0.09 ± 0.02; FBR: h 2 = 0.08 ± 0.02; CF: h 2 = 0.06 ± 0.02; and CFL: h 2 = 0.003 ± 0.004). G × E interaction were negligible for all traits for shrimp reared in different ponds, suggesting re-ranking is not prevalent for this population. Genetic correlations among growth related morphological traits were high ranging from 0.36 to 0.99, suggesting these traits can be simultaneously improved through indirect genetic selection. However, negative genetic correlations were observed for FBR & CF shape traits with major growth traits (ranged -0.08 to -0.95), suggesting genetic selection for rapid growth will likely result in thick/fatty shrimp over generations. Our study showed image-based landmark data can be successfully employed for genetic evaluation of complex morphological traits of shrimp and is potentially amenable to machine-learning derived parameters in semi-automated high volume phenotyping systems needed under commercial conditions.
ABSTRACT
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.
ABSTRACT
OBJECTIVE: This study was aimed to estimate the genetic parameters, including genetic and phenotypic correlations, of milk yield, lactation curve traits and milk composition of Thai dairy cattle from three government research farms. METHODS: The data of 25,789 test-day milk yield and milk composition records of 1,468 cattle from lactation 1 to 3 of Holstein Friesian (HF) and crossbred HF dairy cattle calved between 1990 and 2015 from three government research farms in Thailand were analysed. 305-day milk yield was estimated by the Wood model and a test interval method. The Wood model was used for estimating cumulative 305-day milk yield, peak milk yield, days to peak milk yield and persistency. Genetic parameters were estimated using linear mixed models with herd, breed group, year and season of calving as fixed effects, and animals linked to a pedigree as random effects, together with a residual error. Univariate models were used to estimate variance components, heritability, estimated breeding values (EBVs) and repeatability of each trait, while pairwise bivariate models were used to estimate covariance components and correlations between traits in the same lactation and in the same trait across lactations. RESULTS: The heritability of 305-day milk yield, peak milk yield and protein percentage have moderate to high estimates ranging from 0.19 to 0.45 while days to peak milk yield, persistency and fat percentage have low heritability ranging from 0.08 to 0.14 in lactation 1 cows. Further, heritability of most traits considered was higher in lactation 1 compared with lactations 2 and 3. For cows in lactation 1, high genetic correlations were found between 305-day milk yield and peak milk yield (0.86±0.07) and days to peak milk yield and persistency (0.99±0.02) while estimates of genetic correlations between the remaining traits were imprecise due to the high standard errors. The genetic correlations within the traits across lactation were high. There was no consistent trend of EBVs for most traits in the first lactation over the study period. CONCLUSION: Both the Wood model and test interval method can be used for milk yield estimates in these herds. However, the Wood model has advantages over the test interval method as it can be fitted using fewer test-day records and the estimated model parameters can be used to derive estimates of other lactation curve parameters. Milk yield, peak milk yield and protein percentage can be improved by a selection and mating program while days to peak milk yield, persistency and fat percentage can be improved by including into a selection index.
ABSTRACT
The Arabian oryx was the first species to be rescued from extinction in the wild by the concerted efforts of captive programmes in zoos and private collections around the world. Reintroduction efforts have used two main sources: the 'World Herd', established at the Phoenix Zoo, and private collections in Saudi Arabia. The breeding programme at the Al-Wusta Wildlife Reserve (WWR) in Oman has played a central role in the rescue of the oryx. Individuals from the 'World Herd' and the United Arab Emirates have been the main source for the WWR programme. However, no breeding strategies accounting for genetic diversity have been implemented. To address this, we investigated the diversity of the WWR population and historical samples using mitochondrial DNA (mtDNA) and single nucleotide polymorphisms (SNPs). We found individuals at WWR contain 58% of the total mtDNA diversity observed globally. Inference of ancestry and spatial patterns of SNP variation shows the presence of three ancestral sources and three different groups of individuals. Similar levels of diversity and low inbreeding were observed between groups. We identified individuals and groups that could most effectively contribute to maximizing genetic diversity. Our results will be valuable to guide breeding and reintroduction programmes at WWR.
ABSTRACT
COVID-19 has affected all countries. Its containment represents a unique challenge for India due to a large population (> 1.38 billion) across a wide range of population densities. Assessment of the COVID-19 disease burden is required to put the disease impact into context and support future pandemic policy development. Here, we present the national-level burden of COVID-19 in India in 2020 that accounts for differences across urban and rural regions and across age groups. Input data were collected from official records or published literature. The proportion of excess COVID-19 deaths was estimated using the Institute for Health Metrics and Evaluation, Washington data. Disability-adjusted life years (DALY) due to COVID-19 were estimated in the Indian population in 2020, comprised of years of life lost (YLL) and years lived with disability (YLD). YLL was estimated by multiplying the number of deaths due to COVID-19 by the residual standard life expectancy at the age of death due to the disease. YLD was calculated as a product of the number of incident cases of COVID-19, disease duration and disability weight. Scenario analyses were conducted to account for excess deaths not recorded in the official data and for reported COVID-19 deaths. The direct impact of COVID-19 in 2020 in India was responsible for 14,100,422 (95% uncertainty interval [UI] 14,030,129-14,213,231) DALYs, consisting of 99.2% (95% UI 98.47-99.64%) YLLs and 0.80% (95% UI 0.36-1.53) YLDs. DALYs were higher in urban (56%; 95% UI 56-57%) than rural areas (44%; 95% UI 43.4-43.6) and in men (64%) than women (36%). In absolute terms, the highest DALYs occurred in the 51-60-year-old age group (28%) but the highest DALYs per 100,000 persons were estimated for the 71-80 years old age group (5481; 95% UI 5464-5500 years). There were 4,815,908 (95% UI 4,760,908-4,924,307) DALYs after considering reported COVID-19 deaths only. The DALY estimations have direct and immediate implications not only for public policy in India, but also internationally given that India represents one sixth of the world's population.
Subject(s)
COVID-19/prevention & control , Disability-Adjusted Life Years , Public Health/statistics & numerical data , Quality-Adjusted Life Years , Adolescent , Adult , Aged , Aged, 80 and over , COVID-19/epidemiology , COVID-19/virology , Child , Female , Humans , India/epidemiology , Male , Middle Aged , Pandemics/prevention & control , Public Health/methods , Rural Population/statistics & numerical data , SARS-CoV-2/physiology , Urban Population/statistics & numerical data , Young AdultABSTRACT
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.
Subject(s)
Densovirinae , Penaeidae , Animals , Australia , Densovirinae/genetics , Genome, Viral , Penaeidae/genetics , Polymerase Chain ReactionABSTRACT
[This corrects the article DOI: 10.1371/journal.pone.0238697.].
ABSTRACT
Bovine brucellosis is a neglected zoonotic disease prevalent in several developing countries including India. It has been successfully controlled in many developed countries by using vaccination in conjunction with extensive surveillance and test-and-cull approaches, but some of these approaches do not suit Indian culture and norms. This study was conducted to investigate the feasibility and social acceptability of various bovine brucellosis control strategies in India. Focus group discussions and key-informant interviews were conducted with veterinarians, para-veterinarians, veterinary academics, farmers and other stakeholders. Vaccination with the Brucella strain 19 vaccine was considered feasible, but the participants were concerned about the risk of self-inoculation, the inability to vaccinate pregnant and male animals, the difficulty to differentiate vaccinated from diseased animals and the challenges of maintaining the vaccine cold chain in India. As expected, the test-and-cull approach was not considered feasible as cattle are considered sacred by Hindus and their slaughter is banned in most states. Although the test-and-segregation approach appears reasonable in theory, it would have low acceptability, if implemented without providing any compensation to farmers. Negligible biosecurity was implemented by farmers: almost no biosecurity procedures were performed for visitors entering a farm, and testing of animals was rarely undertaken before introducing them to a farm. However, the participants considered that improving biosecurity would be more acceptable and feasible than both the test-and-cull and the test-and-segregation approaches. Similarly, inadequate personal protection was used by veterinary personnel for handling parturition, retention of placenta and abortion cases; this was considered as another area of possible improvement. Farmers and veterinarians expressed serious concerns about stray cattle as many of them could potentially be infected with brucellosis, and thus could spread the infection between farms. This study recommends using vaccination and biosecurity along with some ancillary strategies to control brucellosis in India. Information from the study could be used to develop an evidence-based disease control program for the disease in the country.
Subject(s)
Brucellosis, Bovine , Communicable Disease Control/methods , Abortion, Veterinary , Animals , Brucellosis, Bovine/prevention & control , Cattle , Farmers , Feasibility Studies , Female , Humans , India , Male , Pregnancy , Vaccination/veterinaryABSTRACT
Niemann-Pick type C disease is a lysosomal storage disease affecting primarily the nervous system that results in premature death. Here we present the first report and investigation of Niemann-Pick type C disease in Australian Angus/Angus-cross calves. After a preliminary diagnosis of Niemann-Pick type C, samples from two affected calves and two obligate carriers were analysed using single nucleotide polymorphism genotyping and homozygosity mapping, and NPC1 was considered as a positional candidate gene. A likely causal missense variant on chromosome 24 in the NPC1 gene (NM_174758.2:c.2969C>G) was identified by Sanger sequencing of cDNA. SIFT analysis, protein alignment and protein modelling predicted the variant to be deleterious to protein function. Segregation of the variant with disease was confirmed in two additional affected calves and two obligate carrier dams. Genotyping of 403 animals from the original herd identified an estimated allele frequency of 3.5%. The Niemann-Pick type C phenotype was additionally confirmed via biochemical analysis of Lysotracker Green, cholesterol, sphingosine and glycosphingolipids in fibroblast cell cultures originating from two affected calves. The identification of a novel missense variant for Niemann-Pick type C disease in Angus/Angus-cross cattle will enable improved breeding and management of this disease in at-risk populations. The results from this study offer a unique opportunity to further the knowledge of human Niemann-Pick type C disease through the potential availability of a bovine model of disease.
Subject(s)
Niemann-Pick Disease, Type C/genetics , Niemann-Pick Disease, Type C/pathology , Amino Acid Sequence , Animals , Base Sequence , Cattle , Cells, Cultured , Cholera Toxin/metabolism , Cholesterol/metabolism , DNA, Complementary/genetics , Disease Models, Animal , Fibroblasts/pathology , G(M1) Ganglioside/metabolism , Homozygote , Mutation/genetics , Niemann-Pick C1 Protein/chemistry , Niemann-Pick C1 Protein/genetics , Phenotype , Polymorphism, Single Nucleotide/genetics , Polysaccharides/metabolism , RNA, Messenger/genetics , RNA, Messenger/metabolismABSTRACT
BACKGROUND: 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). RESULTS: Based on a single reduced-representation library, 31,262 polymorphic SNPs were identified across 650 individuals obtained from Australian wild stocks and commercial aquaculture populations. After filtering to remove SNPs with low read depth, low MAF, low call rate, deviation from HWE, and non-Mendelian inheritance, 7542 high-quality SNPs were retained. From these, 4236 high-quality genome-wide loci were selected for baits-probe development and 4194 SNPs were included within a finalized target-capture genotype-by-sequence assay (DArTcap). This assay was designed for routine and cost effective commercial application in large scale breeding programs, and demonstrates higher confidence in genotype calls through increased call rate (from 80.2 ± 14.7 to 93.0% ± 3.5%), increased read depth (from 20.4 ± 15.6 to 80.0 ± 88.7), as well as a 3-fold reduction in cost over traditional genotype-by-sequencing approaches. CONCLUSION: Importantly, this assay equips the P. monodon industry with the ability to simultaneously assign parentage of communally reared animals, undertake genomic relationship analysis, manage mate pairings between cryptic family lines, as well as undertake advance studies of genome and trait architecture. Critically this assay can be cost effectively applied as P. monodon breeding programs transition to undertaking genomic selection.
Subject(s)
Penaeidae , Animals , Australia , Genome , Genomics , Genotype , Penaeidae/genetics , Polymorphism, Single NucleotideABSTRACT
Camel is an important domestic animal that is well adapted to extremely harsh environments. Due to its multi-purpose role, the camel is gaining importance, particularly as a meat producer. Body weights at different ages and growth curves have been extensively studied in several livestock species. However, the information on the camel is very limited. We investigate patterns of variation in birth weight and weaning weight of two prominent breeds of camels in Pakistan, namely Marecha and Lassi. A further objective is not only to develop growth models which will account for systematic differences due to breed and sex but also to allow for individual variation in growth, by the application of appropriate statistical models. Birth weight, weaning, and monthly weights for a period of up to 48 months were recorded for each camel. Growth records of 137 camels consisting of 108 Marecha and 29 Lassi were obtained. Differences were detected between the two breeds, Marecha and Lassi, as well as differences between sexes, but only after 2 years of age. Beyond that, female Marecha camels had slower growth than males, whereas Lassi (mostly females) had faster growth than female Marechas. In addition, there is a sizable variation between individuals within the same breed-sex combination. The results of the present study will be of use to explore the potential of the camel as a meat producer in Pakistan.
Subject(s)
Camelus/growth & development , Aging , Animals , Camelus/genetics , Female , Male , Models, Biological , PakistanABSTRACT
Unraveling and exploiting mechanisms of disease resistance in cereal crops is currently limited by their large repeat-rich genomes and the lack of genetic recombination or cultivar (cv)-specific sequence information. We cloned the first leaf rust resistance gene Rph1 (Rph1 a) from cultivated barley (Hordeum vulgare) using "MutChromSeq," a recently developed molecular genomics tool for the rapid cloning of genes in plants. Marker-trait association in the CI 9214/Stirling doubled haploid population mapped Rph1 to the short arm of chromosome 2H in a physical region of 1.3 megabases relative to the barley cv Morex reference assembly. A sodium azide mutant population in cv Sudan was generated and 10 mutants were confirmed by progeny-testing. Flow-sorted 2H chromosomes from Sudan (wild type) and six of the mutants were sequenced and compared to identify candidate genes for the Rph1 locus. MutChromSeq identified a single gene candidate encoding a coiled-coil nucleotide binding site Leucine-rich repeat (NLR) receptor protein that was altered in three different mutants. Further Sanger sequencing confirmed all three mutations and identified an additional two independent mutations within the same candidate gene. Phylogenetic analysis determined that Rph1 clustered separately from all previously cloned NLRs from the Triticeae and displayed highest sequence similarity (89%) with a homolog of the Arabidopsis (Arabidopsis thaliana) disease resistance protein 1 protein in Triticum urartu In this study we determined the molecular basis for Rph1-mediated resistance in cultivated barley enabling varietal improvement through diagnostic marker design, gene editing, and gene stacking technologies.
Subject(s)
Hordeum/physiology , Host-Pathogen Interactions , NLR Proteins/physiology , Chromosome Mapping , Genes, Plant , Plant Proteins/physiology , Sequence Analysis, DNAABSTRACT
The black tiger shrimp (Penaeus monodon) remains the second most widely cultured shrimp species globally; however, issues with disease and domestication have seen production levels stagnate over the past two decades. To help identify innovative solutions needed to resolve bottlenecks hampering the culture of this species, it is important to generate genetic and genomic resources. Towards this aim, we have produced the most complete publicly available P. monodon transcriptome database to date based on nine adult tissues and eight early life-history stages (BUSCO - Complete: 98.2% [Duplicated: 51.3%], Fragmented: 0.8%, Missing: 1.0%). The assembly resulted in 236,388 contigs, which were then further segregated into 99,203 adult tissue specific and 58,678 early life-history stage specific clusters. While annotation rates were low (approximately 30%), as is typical for a non-model organisms, annotated transcript clusters were successfully mapped to several hundred functional KEGG pathways. Transcripts were clustered into groups within tissues and early life-history stages, providing initial evidence for their roles in specific tissue functions, or developmental transitions. We expect the transcriptome to provide an essential resource to investigate the molecular basis of commercially relevant-significant traits in P. monodon and other shrimp species.
Subject(s)
Gene Expression Regulation, Developmental , Genome/genetics , Penaeidae/genetics , Transcriptome/genetics , Animals , Aquaculture , Gene Expression Profiling , Multigene Family/genetics , Quantitative Trait Loci/genetics , RNA, Long Noncoding/geneticsABSTRACT
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.
ABSTRACT
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.
Subject(s)
Genetic Linkage , Genomics , Linkage Disequilibrium , Penaeidae/classification , Penaeidae/genetics , Quantitative Trait Loci , Animals , Breeding , Chromosome Mapping , Computational Biology/methods , Female , Gene Ontology , Genetic Heterogeneity , Genomics/methods , High-Throughput Nucleotide Sequencing , Male , Molecular Sequence Annotation , Polymorphism, Single Nucleotide , Recombination, Genetic , Selection, Genetic , Sex FactorsABSTRACT
Evaluating the genetic contribution of individuals to population structure is essential to select informative individuals for genome sequencing, genotype imputation and to ascertain complex population structures. Existing methods for the selection of informative individuals for genomic imputation solely focus on the identification of key ancestors, which can lead to a loss of phasing accuracy of the reference population. Currently many methods are independently applied to investigate complex population structures. Based on the Eigenvalue Decomposition (EVD) of a genomic relationship matrix we describe a novel approach to evaluate the genetic contribution of individuals to population structure. We combined the identification of key contributors with model-based clustering and population network visualization into an integrated three-step approach, which allows identification of high-resolution population structures and substructures around such key contributors. The approach was applied and validated in four disparate datasets including a simulated population (5,100 individuals and 10,000 SNPs), a highly structured experimental sheep population (1,421 individuals and 44,693 SNPs) and two large complex pedigree populations namely horse (1,077 individuals and 38,124 SNPs) and cattle (2,457 individuals and 45,765 SNPs). In the simulated and experimental sheep dataset, our method, which is unsupervised, successfully identified all known key contributors. Applying our three-step approach to the horse and cattle populations, we observed high-resolution population substructures including the absence of obvious important key contributors. Furthermore, we show that compared to commonly applied strategies to select informative individuals for genotype imputation including the computation of marginal gene contributions (Pedig) and the optimization of genetic relatedness (Rel), the selection of key contributors provided the highest phasing accuracies within the selected reference populations. The presented approach opens new perspectives in the characterization and informed management of populations in general, and in areas such as conservation genetics and selective animal breeding in particular, where assessing the genetic contribution of influential and admixed individuals is crucial for research and management applications. As such, this method provides a valuable complement to common applied tools to visualize complex population structures and to select individuals for re-sequencing.
Subject(s)
Genetics, Population , Models, Genetic , Algorithms , Animals , Cattle , Computer Simulation , Horses , Reproducibility of Results , Sheep , WorkflowABSTRACT
Since domestication, significant genetic improvement has been achieved for many traits of commercial importance in cattle, including adaptation, appearance and production. In response to such intense selection pressures, the bovine genome has undergone changes at the underlying regions of functional genetic variants, which are termed "selection signatures". This article reviews 64 recent (2009-2015) investigations testing genomic diversity for departure from neutrality in worldwide cattle populations. In particular, we constructed a meta-assembly of 16,158 selection signatures for individual breeds and their archetype groups (European, African, Zebu and composite) from 56 genome-wide scans representing 70,743 animals of 90 pure and crossbred cattle breeds. Meta-selection-scores (MSS) were computed by combining published results at every given locus, within a sliding window span. MSS were adjusted for common samples across studies and were weighted for significance thresholds across and within studies. Published selection signatures show extensive coverage across the bovine genome, however, the meta-assembly provides a consensus profile of 263 genomic regions of which 141 were unique (113 were breed-specific) and 122 were shared across cattle archetypes. The most prominent peaks of MSS represent regions under selection across multiple populations and harboured genes of known major effects (coat color, polledness and muscle hypertrophy) and genes known to influence polygenic traits (stature, adaptation, feed efficiency, immunity, behaviour, reproduction, beef and dairy production). As the first meta-assembly of selection signatures, it offers novel insights about the hotspots of selective sweeps in the bovine genome, and this method could equally be applied to other species.