Quantitative trait locus analysis of hybrid pedigrees: variance-components model, inbreeding parameter, and power

<p>Abstract</p> <p>Background</p> <p>For the last years reliable mapping of quantitative trait loci (QTLs) has become feasible through linkage analysis based on the variance-components method. There are now many approaches to the QTL analysis of various types of crosses within one population (breed) as well as crosses between divergent populations (breeds). However, to analyse a complex pedigree with dominance and inbreeding, when the pedigree's founders have an inter-population (hybrid) origin, it is necessary to develop a high-powered method taking into account these features of the pedigree.</p> <p>Results</p> <p>We offer a universal approach to QTL analysis of complex pedigrees descended from crosses between outbred parental lines with different QTL allele frequencies. This approach improves the established variance-components method due to the consideration of the genetic effect conditioned by inter-population origin and inbreeding of individuals. To estimate model parameters, namely additive and dominant effects, and the allelic frequencies of the QTL analysed, and also to define the QTL positions on a chromosome with respect to genotyped markers, we used the maximum-likelihood method. To detect linkage between the QTL and the markers we propose statistics with a non-central χ²-distribution that provides the possibility to deduce analytical expressions for the power of the method and therefore, to estimate the pedigree's size required for 80% power. The method works for arbitrarily structured pedigrees with dominance and inbreeding.</p> <p>Conclusion</p> <p>Our method uses the phenotypic values and the marker information for each individual of the pedigree under observation as initial data and can be valuable for fine mapping purposes. The power of the method is increased if the QTL effects conditioned by inter-population origin and inbreeding are enhanced. Several improvements can be developed to take into account fixed factors affecting trait formation, such as age and sex.</p

Exploring deep learning for complex trait genomic prediction in polyploid outcrossing species

Genomic prediction (GP) is the procedure whereby the genetic merits of untested candidates are predicted using genome wide marker information. Although numerous examples of GP exist in plants and animals, applications to polyploid organisms are still scarce, partly due to limited genome resources and the complexity of this system. Deep learning (DL) techniques comprise a heterogeneous collection of machine learning algorithms that have excelled at many prediction tasks. A potential advantage of DL for GP over standard linear model methods is that DL can potentially take into account all genetic interactions, including dominance and epistasis, which are expected to be of special relevance in most polyploids. In this study, we evaluated the predictive accuracy of linear and DL techniques in two important small fruits or berries: strawberry and blueberry. The two datasets contained a total of 1,358 allopolyploid strawberry (2n=8x=112) and 1,802 autopolyploid blueberry (2n=4x=48) individuals, genotyped for 9,908 and 73,045 single nucleotide polymorphism (SNP) markers, respectively, and phenotyped for five agronomic traits each. DL depends on numerous parameters that influence performance and optimizing hyperparameter values can be a critical step. Here we show that interactions between hyperparameter combinations should be expected and that the number of convolutional filters and regularization in the first layers can have an important effect on model performance. In terms of genomic prediction, we did not find an advantage of DL over linear model methods, except when the epistasis component was important. Linear Bayesian models were better than convolutional neural networks for the full additive architecture, whereas the opposite was observed under strong epistasis. However, by using a parameterization capable of taking into account these non-linear effects, Bayesian linear models can match or exceed the predictive accuracy of DL. A semiautomatic implementation of the DL pipeline is available at https://github.com/lauzingaretti/deepGP/

Genetical genomics: spotlight on QTL hotspots

No abstract available

Association between plasma metabolites and gene expression profiles in five porcine endocrine tissues

Background: Endocrine tissues play a fundamental role in maintaining homeostasis of plasma metabolites such as non-esterified fatty acids and glucose, the levels of which reflect the energy balance or the health status of animals. However, the relationship between the transcriptome of endocrine tissues and plasma metabolites has been poorly studied. Methods: We determined the blood levels of 12 plasma metabolites in 27 pigs belonging to five breeds, each breed consisting of both females and males. The transcriptome of five endocrine tissues i.e. hypothalamus, adenohypophysis, thyroid gland, gonads and backfat tissues from 16 out of the 27 pigs was also determined. Sex and breed effects on the 12 plasma metabolites were investigated and associations between genes expressed in the five endocrine tissues and the 12 plasma metabolites measured were analyzed. A probeset was defined as a quantitative trait transcript (QTT) when its association with a particular metabolic trait achieved a nominal P value < 0.01. Results: A larger than expected number of QTT was found for non-esterified fatty acids and alanine aminotransferase in at least two tissues. The associations were highly tissue-specific. The QTT within the tissues were divided into co-expression network modules enriched for genes in Kyoto Encyclopedia of Genes and Genomes or gene ontology categories that are related to the physiological functions of the corresponding tissues. We also explored a multi-tissue co-expression network using QTT for non-esterified fatty acids from the five tissues and found that a module, enriched in hypothalamus QTT, was positioned at the centre of the entire multi-tissue network. Conclusions: These results emphasize the relationships between endocrine tissues and plasma metabolites in terms of gene expression. Highly tissue-specific association patterns suggest that candidate genes or gene pathways should be investigated in the context of specific tissues

Sheep Genome Functional Annotation Reveals Proximal Regulatory Elements Contributed to The Evolution of Modern Breeds

Domestication fundamentally reshaped animal morphology, physiology and behaviour, offering the opportunity to investigate the molecular processes driving evolutionary change. Here we assess sheep domestication and artificial selection by comparing genome sequence from 43 modern breeds (Ovis aries) and their Asian mouflon ancestor (O. orientalis) to identify selection sweeps. Next, we provide a comparative functional annotation of the sheep genome, validated using experimental ChIP-Seq of sheep tissue. Using these annotations, we evaluate the impact of selection and domestication on regulatory sequences and find that sweeps are significantly enriched for protein coding genes, proximal regulatory elements of genes and genome features associated with active transcription. Finally, we find individual sites displaying strong allele frequency divergence are enriched for the same regulatory features. Our data demonstrate that remodelling of gene expression is likely to have been one of the evolutionary forces that drove phenotypic diversification of this common livestock species

DNA Status on Thawed Semen from Fighting Bull: A Comparison Between the SCD and the SCSA Tests

P. 424-431The assessment of sperm chromatin status is compulsory in a complete spermiogram. Here we applied the sperm chromatin structure assay (SCSA) and the sperm chromatin dispersion (SCD) test to assess the chromatin status of three fighting bulls. Cryopreserved semen (two straws/bull) were analysed by duplicate after thawing and after 6 h at 37°C with and without oxidative stress (1 mm FE2+). Results (SCD: percentage of spermatozoa with halo; SCSA: SD‐DFI, %DFI and HDS) were analysed for differences between bulls and treatments, sensitivity and specificity (receiver operating characteristic curves) and repeatability (repeatability coefficients as 2SD of duplicate differences).%DFI for the three bulls was below 2% at 0 h, indicating no risk for fertility according to previous reports. It increased slightly for two of the bulls after FE2+ treatment (%DFI < 5%) and more pronouncedly for the other bull (C, %DFI∼10%), which merits further investigation. SCD rendered higher percentage of halos for bull C, but could not discriminate between samples with and without oxidizing treatment (AUC: 0.52). SCSA (%DFI) showed a high discriminating ability between treatments (AUC: 0.96). The repeatability coefficient was also higher for SCD (5.9) than for %DFI (1.8), indicating lower repeatability for SCD. Overall, %DFI might be the most useful parameter for assessing sperm chromatin on fighting bull. SCD might yield different information than SCSA, hence further research is warranted.S

Fine mapping and replication of QTL in outbred chicken advanced intercross lines

Background: Linkage mapping is used to identify genomic regions affecting the expression of complex traits. However, when experimental crosses such as F2 populations or backcrosses are used to map regions containing a Quantitative Trait Locus (QTL), the size of the regions identified remains quite large, i.e. 10 or more Mb. Thus, other experimental strategies are needed to refine the QTL locations. Advanced Intercross Lines (AIL) are produced by repeated intercrossing of F2 animals and successive generations, which decrease linkage disequilibrium in a controlled manner. Although this approach is seen as promising, both to replicate QTL analyses and fine-map QTL, only a few AIL datasets, all originating from inbred founders, have been reported in the literature. Methods: We have produced a nine-generation AIL pedigree (n = 1529) from two outbred chicken lines divergently selected for body weight at eight weeks of age. All animals were weighed at eight weeks of age and genotyped for SNP located in nine genomic regions where significant or suggestive QTL had previously been detected in the F2 population. In parallel, we have developed a novel strategy to analyse the data that uses both genotype and pedigree information of all AIL individuals to replicate the detection of and fine-map QTL affecting juvenile body weight. Results: Five of the nine QTL detected with the original F2 population were confirmed and fine-mapped with the AIL, while for the remaining four, only suggestive evidence of their existence was obtained. All original QTL were confirmed as a single locus, except for one, which split into two linked QTL. Conclusions: Our results indicate that many of the QTL, which are genome-wide significant or suggestive in the analyses of large intercross populations, are true effects that can be replicated and fine-mapped using AIL. Key factors for success are the use of large populations and powerful statistical tools. Moreover, we believe that the statistical methods we have developed to efficiently study outbred AIL populations will increase the number of organisms for which in-depth complex traits can be analyzed

A gene frequency model for QTL mapping using Bayesian inference

<p>Abstract</p> <p>Background</p> <p>Information for mapping of quantitative trait loci (QTL) comes from two sources: linkage disequilibrium (non-random association of allele states) and cosegregation (non-random association of allele origin). Information from LD can be captured by modeling conditional means and variances at the QTL given marker information. Similarly, information from cosegregation can be captured by modeling conditional covariances. Here, we consider a Bayesian model based on gene frequency (BGF) where both conditional means and variances are modeled as a function of the conditional gene frequencies at the QTL. The parameters in this model include these gene frequencies, additive effect of the QTL, its location, and the residual variance. Bayesian methodology was used to estimate these parameters. The priors used were: logit-normal for gene frequencies, normal for the additive effect, uniform for location, and inverse chi-square for the residual variance. Computer simulation was used to compare the power to detect and accuracy to map QTL by this method with those from least squares analysis using a regression model (LSR).</p> <p>Results</p> <p>To simplify the analysis, data from unrelated individuals in a purebred population were simulated, where only LD information contributes to map the QTL. LD was simulated in a chromosomal segment of 1 cM with one QTL by random mating in a population of size 500 for 1000 generations and in a population of size 100 for 50 generations. The comparison was studied under a range of conditions, which included SNP density of 0.1, 0.05 or 0.02 cM, sample size of 500 or 1000, and phenotypic variance explained by QTL of 2 or 5%. Both 1 and 2-SNP models were considered. Power to detect the QTL for the BGF, ranged from 0.4 to 0.99, and close or equal to the power of the regression using least squares (LSR). Precision to map QTL position of BGF, quantified by the mean absolute error, ranged from 0.11 to 0.21 cM for BGF, and was better than the precision of LSR, which ranged from 0.12 to 0.25 cM.</p> <p>Conclusions</p> <p>In conclusion given a high SNP density, the gene frequency model can be used to map QTL with considerable accuracy even within a 1 cM region.</p