Whole-genome sequencing of two captive black soldier fly populations: Implications for commercial production.

Black soldier fly (BSF; Hermetia illucens) is a promising insect species for food and feed production as its larvae can convert different organic waste to high-value protein. Selective breeding is one way to optimize production, but the potential of breeding is only starting to be explored and not yet utilized for BSF. To assist in monitoring a captive population and implementing a breeding program, genomics tools are imperative. We conducted whole genome sequencing of two captive populations separated by geographical distance - Denmark (DK) and Texas, USA (TX). Various population genetics analyses revealed a moderate genetic differentiation between two populations. Moreover, we observed higher inbreeding in the DK population, and the detection of a subpopulation within DK population aligned well with the recent foundation of the DK population from two captive populations. Additionally, we generated gene ontology annotation and variants annotation for wider potential applications. Our findings establish a robust marker set for research in population genetics, facilitating the monitoring of inbreeding and laying the groundwork for practical breeding programs for BSF.

Estimation of genetic parameters for the implementation of selective breeding in commercial insect production.

BACKGROUND: There is a burgeoning interest in using insects as a sustainable source of food and feed, particularly by capitalising on various waste materials and by-products that are typically considered of low value. Enhancing the commercial production of insects can be achieved through two main approaches: optimising environmental conditions and implementing selective breeding strategies. In order to successfully target desirable traits through selective breeding, having a thorough understanding of the genetic parameters pertaining to those traits is essential. In this study, a full-sib half-sib mating design was used to estimate variance components and heritabilities for larval size and survival at day seven of development, development time and survival from egg to adult, and to estimate correlations between these traits, within an outbred population of house flies (Musca domestica), using high-throughput phenotyping for data collection. RESULTS: The results revealed low to intermediate heritabilities and positive genetic correlations between all traits except development time and survival to day seven of development and from egg to adulthood. Surprisingly, larval size at day seven exhibited a comparatively low heritability (0.10) in contrast to development time (0.25), a trait that is believed to have a stronger association with overall fitness. A decline in family numbers resulting from low mating success and high overall mortality reduced the amount of available data which resulted in large standard errors for the estimated parameters. Environmental factors made a substantial contribution to the phenotypic variation, which was overall high for all traits. CONCLUSIONS: There is potential for genetic improvement in all studied traits and estimates of genetic correlations indicate a partly shared genetic architecture among the traits. All estimates have large standard errors. Implementing high-throughput phenotyping is imperative for the estimation of genetic parameters in fast developing insects, and facilitates age synchronisation, which is vital in a breeding population. In spite of endeavours to minimise non-genetic sources of variation, all traits demonstrated substantial influences from environmental components. This emphasises the necessity of thorough attention to the experimental design before breeding is initiated in insect populations.

Impact of kinship matrices on genetic gain and inbreeding with optimum contribution selection in a genomic dairy cattle breeding program.

BACKGROUND: Genomic selection has increased genetic gain in dairy cattle, but in some cases it has resulted in higher inbreeding rates. Therefore, there is need for research on efficient management of inbreeding in genomically-selected dairy cattle populations, especially for local breeds with a small population size. Optimum contribution selection (OCS) minimizes the increase in average kinship while it maximizes genetic gain. However, there is no consensus on how to construct the kinship matrix used for OCS and whether it should be based on pedigree or genomic information. VanRaden's method 1 (VR1) is a genomic relationship matrix in which centered genotype scores are scaled with the sum of 2p(1-p) where p is the reference allele frequency at each locus, and VanRaden's method 2 (VR2) scales each locus with 2p(1-p), thereby giving greater weight to loci with a low minor allele frequency. We compared the effects of nine kinship matrices on genetic gain, kinship, inbreeding, genetic diversity, and minor allele frequency when applying OCS in a simulated small dairy cattle population. We used VR1 and VR2, each using base animals, all genotyped animals, and the current generation of animals to compute reference allele frequencies. We also set the reference allele frequencies to 0.5 for VR1 and the pedigree-based relationship matrix. We constrained OCS to select a fixed number of sires per generation for all scenarios. Efficiency of the different matrices were compared by calculating the rate of genetic gain for a given rate of increase in average kinship. RESULTS: We found that: (i) genomic relationships were more efficient than pedigree-based relationships at managing inbreeding, (ii) reference allele frequencies computed from base animals were more efficient compared to reference allele frequencies computed from recent animals, and (iii) VR1 was slightly more efficient than VR2, but the difference was not statistically significant. CONCLUSIONS: Using genomic relationships for OCS realizes more genetic gain for a given amount of kinship and inbreeding than using pedigree relationships when the number of sires is fixed. For a small genomic dairy cattle breeding program, we recommend that the implementation of OCS uses VR1 with reference allele frequencies estimated either from base animals or old genotyped animals.

Genome-wide association study identifies functional genomic variants associated with young stock survival in Nordic Red Dairy Cattle.

Identifying quantitative trait loci (QTL) associated with calf survival is essential for both reducing economic loss in cattle industry and understanding the genetic basis of the trait. To identify mutations and genes underlying young stock survival (YSS), we performed GWAS using de-regressed estimated breeding values of a YSS index and its component traits defined by sex and age in 3,077 Nordic Red Dairy Cattle (RDC) bulls and 2 stillbirth traits (first lactation and later lactations) in 5,141 RDC bulls. Two associated QTL regions on Bos taurus autosome (BTA) 4 and 6 were identified for the YSS index. The results of 4 YSS component traits indicate that same QTL regions were associated with bull and heifer calf mortality, but the effects were different over the growing period and suggested an additional QTL on BTA23. The GWAS on stillbirth identified 3 additional QTL regions on BTA5, 14, and 24 compared with YSS and its component traits. The conditional test of BTA6 showed at least 2 closely located QTL segregating for YSS component traits and stillbirth. We found 2 independent QTL for stillbirth on BTA23. The post-GWAS revealed LCORL, PPM1K, SSP1, MED28, and LAP3 are putative causal genes on BTA6, and a frame shift variant within LCORL, BTA6:37401770 (rs384548488) could be the putative causal variant. On BTA4, the GRB10 gene is the putative causal gene and BTA4:5296018 is the putative causal variant. In addition, NDUFA9 and FGF23 on BTA5, LYN on BTA14, and KCNK5 on BTA23 are putative causal genes for QTL for stillbirth. The gene analysis also proposed several candidate genes. Our findings shed new light on the candidate genes affecting calf survival, and the knowledge could be utilized to reduce calf mortality and thereby enhance welfare of dairy cattle.

Large-scale association study on daily weight gain in pigs reveals overlap of genetic factors for growth in humans.

BACKGROUND: Imputation from genotyping array to whole-genome sequence variants using resequencing of representative reference populations enhances our ability to map genetic factors affecting complex phenotypes in livestock species. The accumulation of knowledge about gene function in human and laboratory animals can provide substantial advantage for genomic research in livestock species. RESULTS: In this study, 201,388 pigs from three commercial Danish breeds genotyped with low to medium (8.5k to 70k) SNP arrays were imputed to whole genome sequence variants using a two-step approach. Both imputation steps achieved high accuracies, and in total this yielded 26,447,434 markers on 18 autosomes. The average estimated imputation accuracy of markers with minor allele frequency ≥ 0.05 was 0.94. To overcome the memory consumption of running genome-wide association study (GWAS) for each breed, we performed within-breed subpopulation GWAS then within-breed meta-analysis for average daily weight gain (ADG), followed by a multi-breed meta-analysis of GWAS summary statistics. We identified 15 quantitative trait loci (QTL). Our post-GWAS analysis strategy to prioritize of candidate genes including information like gene ontology, mammalian phenotype database, differential expression gene analysis of high and low feed efficiency pig and human GWAS catalog for height, obesity, and body mass index, we proposed MRAP2, LEPROT, PMAIP1, ENSSSCG00000036234, BMP2, ELFN1, LIG4 and FAM155A as the candidate genes with biological support for ADG in pigs. CONCLUSION: Our post-GWAS analysis strategy helped to identify candidate genes not just by distance to the lead SNP but also by multiple sources of biological evidence. Besides, the identified QTL overlap with genes which are known for their association with human growth-related traits. The GWAS with this large data set showed the power to map the genetic factors associated with ADG in pigs and have added to our understanding of the genetics of growth across mammalian species.

Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions.

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY and DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes, Y1A, Y1B, Y2A and Y2B, with a marked geographical partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and seven wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during migrations into northern Europe, eastern and southern Asia, and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographical range.

Application of mixed linear models for the estimation of functional effects on bovine stature based on SNP summary statistics from a whole-genome association study.

Genome-wide association studies (GWAS) help identify polymorphic sites or genes linked to phenotypic variance, but a few identified genes and/or single nucleotide polymorphisms (SNPs) are unlikely to explain a large part of the phenotypic variability of complex traits. In this study, the focus was moved from single loci to functional units, expressed by the metabolic pathways as defined in the Kyoto Encyclopaedia of Genes and Genomes (KEGG) database. Consequently, the aim of this study was to estimate KEGG effects on stature in three Nordic dairy cattle breeds using SNP effects from GWAS as the dependent variable. The SNPs were annotated to genes, then the genes to KEGG pathways. The effects of KEGG pathways were estimated separately for each breed using a mixed linear model incorporating the similarity between pathways expressed by common genes. The KEGG pathway D-amino acid metabolism (map00473) was estimated to be significant for stature in two of the analysed breeds and revealed a borderline significance in the third breed. Thus, we demonstrate that the approach to statistical modelling of higher order functional effects on complex traits is useful, and provides evidence of the importance of D-amino acids for growth in cattle.

Genome-wide association study with imputed whole-genome sequence variants including large deletions for female fertility in 3 Nordic dairy cattle breeds.

Fertility is an economically important trait in livestock. Poor fertility in dairy cattle can be due to loss-of-function variants affecting any essential gene that causes early embryonic mortality in homozygotes. To identify fertility-associated quantitative trait loci, we performed single-marker association analyses for 8 fertility traits in Holstein, Jersey, and Nordic Red Dairy cattle using imputed whole-genome sequence variants including SNPs, indels, and large deletion. We then performed stepwise selection of independent markers from GWAS loci using conditional and joint association analyses. From single-marker analyses for fertility traits, we reported genome-wide significant associations of 30,384 SNPs, 178 indels, and 3 deletions in Holstein; 23,481 SNPs, 189 indels, and 13 deletions in Nordic Red; and 17 SNPs in Jersey cattle. Conditional and joint association analyses identified 37 and 23 independent associations in Holstein and Nordic Red Dairy cattle, respectively. Fertility-associated GWAS loci were enriched for developmental and cellular processes (Gene Ontology enrichment, false discovery rate < 0.05). For these quantitative trait loci regions (top marker and 500 kb of surrounding regions), we proposed several candidate genes with functional annotations corresponding to embryonic lethality and various fertility-related phenotypes in mouse and cattle. The inclusion of these top markers in future releases of the custom SNP chip used for genomic evaluations will enable their validation in independent populations and improve the accuracy of genomic predictions.

Host genetics and the rumen microbiome jointly associate with methane emissions in dairy cows.

Cattle and other ruminants produce large quantities of methane (~110 million metric tonnes per annum), which is a potent greenhouse gas affecting global climate change. Methane (CH4) is a natural by-product of gastro-enteric microbial fermentation of feedstuffs in the rumen and contributes to 6% of total CH4 emissions from anthropogenic-related sources. The extent to which the host genome and rumen microbiome influence CH4 emission is not yet well known. This study confirms individual variation in CH4 production was influenced by individual host (cow) genotype, as well as the host's rumen microbiome composition. Abundance of a small proportion of bacteria and archaea taxa were influenced to a limited extent by the host's genotype and certain taxa were associated with CH4 emissions. However, the cumulative effect of all bacteria and archaea on CH4 production was 13%, the host genetics (heritability) was 21% and the two are largely independent. This study demonstrates variation in CH4 emission is likely not modulated through cow genetic effects on the rumen microbiome. Therefore, the rumen microbiome and cow genome could be targeted independently, by breeding low methane-emitting cows and in parallel, by investigating possible strategies that target changes in the rumen microbiome to reduce CH4 emissions in the cattle industry.

A selection signatures study among Middle Eastern and European sheep breeds.

Selection, both natural and artificial, leaves patterns on the genome during domestication of animals and leads to changes in allele frequencies among populations. Detecting genomic regions influenced by selection in livestock may assist in understanding the processes involved in genome evolution and discovering genomic regions related to traits of economic and ecological interests. In the current study, genetic diversity analyses were conducted on 34,206 quality-filtered SNP positions from 450 individuals in 15 sheep breeds, including six indigenous breeds from the Middle East, namely Iranian Balouchi, Afshari, Moghani, Qezel, Karakas and Norduz, and nine breeds from Europe, namely East Friesian Sheep, Ile de France, Mourerous, Romane, Swiss Mirror, Spaelsau, Suffolk, Comisana and Engadine Red Sheep. The SNP genotype data generated by the Illumina OvineSNP50 Genotyping BeadChip array were used in this analysis. We applied two complementary statistical analyses, FST (fixation index) and xp-EHH (cross-population extended haplotype homozygosity), to detect selection signatures in Middle Eastern and European sheep populations. FST and xp-EHH detected 629 and 256 genes indicating signatures of selection, respectively. Genomic regions identified using FST and xp-EHH contained the CIDEA, HHATL, MGST1, FADS1, RTL1 and DGKG genes, which were reported earlier to influence a number of economic traits. Both FST and xp-EHH approaches identified 60 shared genes as the signatures of selection, including four candidate genes (NT5E, ADA2, C8A and C8B) that were enriched for two significant Gene Ontology (GO) terms associated with the adenosine metabolic procedure. Knowledge about the candidate genomic regions under selective pressure in sheep breeds may facilitate identification of the underlying genes and enhance our understanding on these genes role in local adaptation.

Novel approach to incorporate information about recessive lethal genes increases the accuracy of genomic prediction for mortality traits.

The genetic underpinnings of calf mortality can be partly polygenic and partly due to deleterious effects of recessive lethal alleles. Prediction of the genetic merits of selection candidates should thus take into account both genetic components contributing to calf mortality. However, simultaneously modeling polygenic risk and recessive lethal allele effects in genomic prediction is challenging due to effects that behave differently. In this study, we present a novel approach where mortality risk probabilities from polygenic and lethal allele components are predicted separately to compute the total risk probability of an individual for its future offspring as a basis for selection. We present methods for transforming genomic estimated breeding values of polygenic effect into risk probabilities using normal density and cumulative distribution functions and show computations of risk probability from recessive lethal alleles given sire genotypes and population recessive allele frequencies. Simulated data were used to test the novel approach as implemented in probit, logit, and linear models. In the simulation study, the accuracy of predicted risk probabilities was computed as the correlation between predicted mortality probabilities and observed calf mortality for validation sires. The results indicate that our novel approach can greatly increase the accuracy of selection for mortality traits compared with the accuracy of predictions obtained without distinguishing polygenic and lethal gene effects.

Distinguishing pleiotropy from linked QTL between milk production traits and mastitis resistance in Nordic Holstein cattle.

BACKGROUND: Production and health traits are central in cattle breeding. Advances in next-generation sequencing technologies and genotype imputation have increased the resolution of gene mapping based on genome-wide association studies (GWAS). Thus, numerous candidate genes that affect milk yield, milk composition, and mastitis resistance in dairy cattle are reported in the literature. Effect-bearing variants often affect multiple traits. Because the detection of overlapping quantitative trait loci (QTL) regions from single-trait GWAS is too inaccurate and subjective, multi-trait analysis is a better approach to detect pleiotropic effects of variants in candidate genes. However, large sample sizes are required to achieve sufficient power. Multi-trait meta-analysis is one approach to deal with this problem. Thus, we performed two multi-trait meta-analyses, one for three milk production traits (milk yield, protein yield and fat yield), and one for milk yield and mastitis resistance. RESULTS: For highly correlated traits, the power to detect pleiotropy was increased by multi-trait meta-analysis compared with the subjective assessment of overlapping of single-trait QTL confidence intervals. Pleiotropic effects of lead single nucleotide polymorphisms (SNPs) that were detected from the multi-trait meta-analysis were confirmed by bivariate association analysis. The previously reported pleiotropic effects of variants within the DGAT1 and MGST1 genes on three milk production traits, and pleiotropic effects of variants in GHR on milk yield and fat yield were confirmed. Furthermore, our results suggested that variants in KCTD16, KCNK18 and ENSBTAG00000023629 had pleiotropic effects on milk production traits. For milk yield and mastitis resistance, we identified possible pleiotropic effects of variants in two genes, GC and DGAT1. CONCLUSIONS: Multi-trait meta-analysis improves our ability to detect pleiotropic interactions between milk production traits and identifies variants with pleiotropic effects on milk production traits and mastitis resistance. In particular, this should contribute to better understand the biological mechanisms that underlie the unfavorable genetic correlation between milk yield and mastitis.

Weighted single-step genomic best linear unbiased prediction integrating variants selected from sequencing data by association and bioinformatics analyses.

BACKGROUND: Sequencing data enable the detection of causal loci or single nucleotide polymorphisms (SNPs) highly linked to causal loci to improve genomic prediction. However, until now, studies on integrating such SNPs using a single-step genomic best linear unbiased prediction (ssGBLUP) model are scarce. We investigated the integration of sequencing SNPs selected by association (1262 SNPs) and bioinformatics (2359 SNPs) analyses into the currently used 54K-SNP chip, using three ssGBLUP models which make different assumptions on the distribution of SNP effects: a basic ssGBLUP model, a so-called featured ssGBLUP (ssFGBLUP) model that considered selected sequencing SNPs as a feature genetic component, and a weighted ssGBLUP (ssWGBLUP) model in which the genomic relationship matrix was weighted by the SNP variances estimated from a Bayesian whole-genome regression model, with every 1, 30, or 100 adjacent SNPs within a chromosome region sharing the same variance. We used data on milk production and female fertility in Danish Jersey. In total, 15,823 genotyped and 528,981‬ non-genotyped females born between 1990 and 2013 were used as reference population and 7415 genotyped females and 33,040 non-genotyped females born between 2014 and 2016 were used as validation population. RESULTS: With basic ssGBLUP, integrating SNPs selected from sequencing data improved prediction reliabilities for milk and protein yields, but resulted in limited or no improvement for fat yield and female fertility. Model performances depended on the SNP set used. When using ssWGBLUP with the 54K SNPs, reliabilities for milk and protein yields improved by 0.028 for genotyped animals and by 0.006 for non-genotyped animals compared with ssGBLUP. However, with the SNP set that included SNPs selected from sequencing data, no statistically significant difference in prediction reliability was observed between the three ssGBLUP models. CONCLUSIONS: In summary, when using 54K SNPs, a ssWGBLUP model with a common weight on the SNPs in a given region is a feasible approach for single-trait genetic evaluation. Integrating relevant SNPs selected from sequencing data into the standard SNP chip can improve the reliability of genomic prediction. Based on such SNP data, a basic ssGBLUP model was suggested since no significant improvement was observed from using alternative models such as ssWGBLUP and ssFGBLUP.

Novel haplotypes responsible for prenatal death in Nordic Red and Danish Jersey cattle.

Haplotypes that are common in a population but not observed as homotypes in living animals may harbor lethal alleles that compromise embryo survival. In this study, we searched for homozygous-deficient haplotypes in the genomes of 19,309 Nordic Red Dairy (RDC) and 4,291 Danish Jersey (JER) cattle genotyped using the Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA). For statistically significant deficient haplotypes, we evaluated the effect on nonreturn rate in at-risk matings (mating between carrier bull and daughter of carrier sire) versus not-at-risk matings (mating between noncarrier bull and daughter of noncarrier sire). Next, we analyzed whole-genome sequence variants from the 1000 Bull Genomes Project to identify putative causal variants underlying these haplotypes. In RDC, we identified 3 homozygous-deficient regions (HDR) that overlapped with known recessive lethal mutations: a 662-kb deletion on chromosome 12 in RDC [Online Mendelian Inheritance in Animals (OMIA) 001901-9913), a missense mutation in TUBD1, g.11063520T>C, in Braunvieh cattle (OMIA 001939-9913), and a 525-kb deletion on chromosome 23 in RDC (OMIA 001991-9913)]. In addition, we identified 15 novel HDR and their tag haplotypes for the underlying causative variants. The tag haplotype located between 39.2 and 40.3 Mbp on chromosome 18 had a negative effect on nonreturn rate in at-risk mating, confirming embryonic lethality. In Danish Jersey, we identified 12 novel HDR and their tag haplotypes for underlying causative variants. For 3 of these 12 tag haplotypes, insemination records of at-risk mating showed a negative effect on nonreturn rate, confirming the association with early embryonic mortality. Cattle that had both genotype and whole-genome sequence data were analyzed to detect the causative variants underlying each tag haplotype. However, none of the functional variants or deletions showed concordance with carrier status of the novel tag haplotypes. Carrier status of these detected haplotypes can be used to select bulls to reduce the frequencies of lethal alleles in the population and to avoid at-risk matings.

Predictive ability of host genetics and rumen microbiome for subclinical ketosis.

Subclinical metabolic disorders such as ketosis cause substantial economic losses for dairy farmers in addition to the serious welfare issues they pose for dairy cows. Major hurdles in genetic improvement against metabolic disorders such as ketosis include difficulties in large-scale phenotype recording and low heritability of traits. Milk concentrations of ketone bodies, such as acetone and ß-hydroxybutyric acid (BHB), might be useful indicators to select cows for low susceptibility to ketosis. However, heritability estimates reported for milk BHB and acetone in several dairy cattle breeds were low. The rumen microbial community has been reported to play a significant role in host energy homeostasis and metabolic and physiologic adaptations. The current study aims at investigating the effects of cows' genome and rumen microbial composition on concentrations of acetone and BHB in milk, and identifying specific rumen microbial taxa associated with variation in milk acetone and BHB concentrations. We determined the concentrations of acetone and BHB in milk using nuclear magnetic resonance spectroscopy on morning milk samples collected from 277 Danish Holstein cows. Imputed high-density genotype data were available for these cows. Using genomic and microbial prediction models with a 10-fold resampling strategy, we found that rumen microbial composition explains a larger proportion of the variation in milk concentrations of acetone and BHB than do host genetics. Moreover, we identified associations between milk acetone and BHB with some specific bacterial and archaeal operational taxonomic units previously reported to have low to moderate heritability, presenting an opportunity for genetic improvement. However, higher covariation between specific microbial taxa and milk acetone and BHB concentrations might not necessarily indicate a causal relationship; therefore further validation is needed before considering implementation in selection programs.

Prioritizing candidate genes for fertility in dairy cows using gene-based analysis, functional annotation and differential gene expression.

BACKGROUND: An unfavorable genetic correlation between milk production and fertility makes simultaneous improvement of milk production and fertility difficult in cattle breeding. Rapid genetic improvement in milk production traits in dairy cattle has been accompanied by decline in cow fertility. The genetic basis of this correlation remains poorly understood. Expanded reference populations and large sets of sequenced animals make genome-wide association studies (GWAS) with imputed markers possible for large populations and thereby studying genetic architecture of complex traits. RESULTS: In this study, we associated 15,551,021 SNPs with female fertility index in 5038 Nordic Holstein cattle. We have identified seven quantitative trait loci (QTL) on six chromosomes in cattle. Along with nearest genes to GWAS hits, we used gene-based analysis and spread of linkage disequilibrium (LD) information to generate a list of potential candidate genes affecting fertility in cattle. Subsequently, we used prior knowledge on gene related to fertility from Gene Ontology terms, Kyoto Encyclopedia of Genes and Genomes pathway analysis, mammalian phenotype database, and public available RNA-seq data to refine the list of candidate genes for fertility. We used variant annotations to investigate candidate mutations within the prioritized candidate genes. Using multiple source of information, we proposed candidate genes with biological relevance underlying each of these seven QTL. On chromosome 1, we have identified ten candidate genes for two QTL. For the rest of chromosomes, we proposed one candidate gene for each QTL. In the candidate genes list, differentially expressed genes from different studies support FRAS1, ITGB5, ADCY5, and SEMA5B as candidate genes for cow fertility. CONCLUSION: The GWAS result not only confirmed previously mapped QTL, but also made new findings. Our findings contributes towards dissecting the genetics for female fertility in cattle. Moreover, this study shows the usefulness of adding independent information to pick candidate genes during post-GWAS analysis.

Whole genome sequence and de novo assembly revealed genomic architecture of Indian Mithun (Bos frontalis).

BACKGROUND: Mithun (Bos frontalis), also called gayal, is an endangered bovine species, under the tribe bovini with 2n = 58 XX chromosome complements and reared under the tropical rain forests region of India, China, Myanmar, Bhutan and Bangladesh. However, the origin of this species is still disputed and information on its genomic architecture is scanty so far. We trust that availability of its whole genome sequence data and assembly will greatly solve this problem and help to generate many information including phylogenetic status of mithun. Recently, the first genome assembly of gayal, mithun of Chinese origin, was published. However, an improved reference genome assembly would still benefit in understanding genetic variation in mithun populations reared under diverse geographical locations and for building a superior consensus assembly. We, therefore, performed deep sequencing of the genome of an adult female mithun from India, assembled and annotated its genome and performed extensive bioinformatic analyses to produce a superior de novo genome assembly of mithun. RESULTS: We generated ≈300 Gigabyte (Gb) raw reads from whole-genome deep sequencing platforms and assembled the sequence data using a hybrid assembly strategy to create a high quality de novo assembly of mithun with 96% recovered as per BUSCO analysis. The final genome assembly has a total length of 3.0 Gb, contains 5,015 scaffolds with an N50 value of 1 Mb. Repeat sequences constitute around 43.66% of the assembly. The genomic alignments between mithun to cattle showed that their genomes, as expected, are highly conserved. Gene annotation identified 28,044 protein-coding genes presented in mithun genome. The gene orthologous groups of mithun showed a high degree of similarity in comparison with other species, while fewer mithun specific coding sequences were found compared to those in cattle. CONCLUSION: Here we presented the first de novo draft genome assembly of Indian mithun having better coverage, less fragmented, better annotated, and constitutes a reasonably complete assembly compared to the previously published gayal genome. This comprehensive assembly unravelled the genomic architecture of mithun to a great extent and will provide a reference genome assembly to research community to elucidate the evolutionary history of mithun across its distinct geographical locations.

A splice donor variant in CCDC189 is associated with asthenospermia in Nordic Red dairy cattle.

BACKGROUND: Cattle populations are highly amenable to the genetic mapping of male reproductive traits because longitudinal data on ejaculate quality and dense microarray-derived genotypes are available for thousands of artificial insemination bulls. Two young Nordic Red bulls delivered sperm with low progressive motility (i.e., asthenospermia) during a semen collection period of more than four months. The bulls were related through a common ancestor on both their paternal and maternal ancestry. Thus, a recessive mode of inheritance of asthenospermia was suspected. RESULTS: Both bulls were genotyped at 54,001 SNPs using the Illumina BovineSNP50 Bead chip. A scan for autozygosity revealed that they were identical by descent for a 2.98 Mb segment located on bovine chromosome 25. This haplotype was not found in the homozygous state in 8557 fertile bulls although five homozygous haplotype carriers were expected (P = 0.018). Whole genome-sequencing uncovered that both asthenospermic bulls were homozygous for a mutation that disrupts a canonical 5' splice donor site of CCDC189 encoding the coiled-coil domain containing protein 189. Transcription analysis showed that the derived allele activates a cryptic splice site resulting in a frameshift and premature termination of translation. The mutated CCDC189 protein is truncated by more than 40%, thus lacking the flagellar C1a complex subunit C1a-32 that is supposed to modulate the physiological movement of the sperm flagella. The mutant allele occurs at a frequency of 2.5% in Nordic Red cattle. CONCLUSIONS: Our study in cattle uncovered that CCDC189 is required for physiological movement of sperm flagella thus enabling active progression of spermatozoa and fertilization. A direct gene test may be implemented to monitor the asthenospermia-associated allele and prevent the birth of homozygous bulls that are infertile. Our results have been integrated in the Online Mendelian Inheritance in Animals (OMIA) database ( https://omia.org/OMIA002167/9913/ ).

Dissecting closely linked association signals in combination with the mammalian phenotype database can identify candidate genes in dairy cattle.

BACKGROUND: Genome-wide association studies (GWAS) have been successfully implemented in cattle research and breeding. However, moving from the associations to identify the causal variants and reveal underlying mechanisms have proven complicated. In dairy cattle populations, we face a challenge due to long-range linkage disequilibrium (LD) arising from close familial relationships in the studied individuals. Long range LD makes it difficult to distinguish if one or multiple quantitative trait loci (QTL) are segregating in a genomic region showing association with a phenotype. We had two objectives in this study: 1) to distinguish between multiple QTL segregating in a genomic region, and 2) use of external information to prioritize candidate genes for a QTL along with the candidate variants. RESULTS: We observed fixing the lead SNP as a covariate can help to distinguish additional close association signal(s). Thereafter, using the mammalian phenotype database, we successfully found candidate genes, in concordance with previous studies, demonstrating the power of this strategy. Secondly, we used variant annotation information to search for causative variants in our candidate genes. The variant information successfully identified known causal mutations and showed the potential to pinpoint the causative mutation(s) which are located in coding regions. CONCLUSIONS: Our approach can distinguish multiple QTL segregating on the same chromosome in a single analysis without manual input. Moreover, utilizing information from the mammalian phenotype database and variant effect predictor as post-GWAS analysis could benefit in candidate genes and causative mutations finding in cattle. Our study not only identified additional candidate genes for milk traits, but also can serve as a routine method for GWAS in dairy cattle.

Weighting sequence variants based on their annotation increases the power of genome-wide association studies in dairy cattle.

BACKGROUND: Genome-wide association studies (GWAS) are widely used to identify regions of the genome that harbor genetic determinants of quantitative traits. However, the multiple-testing burden from scanning tens of millions of whole-genome sequence variants reduces the power to identify associated variants, especially if sample size is limited. In addition, factors such as inaccuracy of imputation, complex linkage disequilibrium structures, and multiple closely-located causal variants may result in an identified causative mutation not being the most significant single nucleotide polymorphism in a particular genomic region. Therefore, the use of information from different sources, particularly variant annotations, was proposed to enhance the fine-mapping of causal variants. Here, we tested whether applying significance thresholds based on variant annotation categories increases the power of GWAS compared with a flat Bonferroni multiple-testing correction. RESULTS: Whole-genome sequence variants in dairy cattle were categorized according to type and predicted impact. Then, GWAS between markers and 17 quantitative traits were analyzed for enrichment for association of each annotation category. By using annotation categories that were determined with the variants effect predictor software and datasets indicating regions of open chromatin, "low impact" variants were found to be highly enriched. Moreover, when the variants annotated as "modifier" and not located at open chromatin regions were further classified into different types of potential regulatory elements, the high impact variants, moderate impact variants, variants located in the 3' and 5' untranslated regions, and variants located in potential non-coding RNA regions exhibited relatively more enrichment. In contrast, a similar study on human GWAS data reported that enrichment of association signals was highest with high impact variants. We observed an increase in power when these variant category-based significance thresholds were applied for GWAS results on stature in Nordic Holstein cattle, as more candidate genes from previous large GWAS meta-analysis for cattle stature were confirmed. CONCLUSIONS: Use of variant category-based genome-wide significance thresholds can marginally increase the power to detect the candidate genes in cattle. With the continued improvements in annotation of the bovine genome, we anticipate that the growing usefulness of variant category-based significance thresholds will be demonstrated.