Identification and characterisation of de novo germline structural variants in two commercial pig lines using trio-based whole genome sequencing

BACKGROUND

De novo mutations arising in the germline are a source of genetic variation and their discovery broadens our understanding of genetic disorders and evolutionary patterns. Although the number of de novo single nucleotide variants (dnSNVs) has been studied in a number of species, relatively little is known about the occurrence of de novo structural variants (dnSVs). In this study, we investigated 37 deeply sequenced pig trios from two commercial lines to identify dnSVs present in the offspring. The identified dnSVs were characterised by identifying their parent of origin, their functional annotations and characterizing sequence homology at the breakpoints.

RESULTS

We identified four swine germline dnSVs, all located in intronic regions of protein-coding genes. Our conservative, first estimate of the swine germline dnSV rate is 0.108 (95% CI 0.038-0.255) per generation (one dnSV per nine offspring), detected using short-read sequencing. Two detected dnSVs are clusters of mutations. Mutation cluster 1 contains a de novo duplication, a dnSNV and a de novo deletion. Mutation cluster 2 contains a de novo deletion and three de novo duplications, of which one is inverted. Mutation cluster 2 is 25 kb in size, whereas mutation cluster 1 (197 bp) and the other two individual dnSVs (64 and 573 bp) are smaller. Only mutation cluster 2 could be phased and is located on the paternal haplotype. Mutation cluster 2 originates from both micro-homology as well as non-homology mutation mechanisms, where mutation cluster 1 and the other two dnSVs are caused by mutation mechanisms lacking sequence homology. The 64 bp deletion and mutation cluster 1 were validated through PCR. Lastly, the 64 bp deletion and the 573 bp duplication were validated in sequenced offspring of probands with three generations of sequence data.

CONCLUSIONS

Our estimate of 0.108 dnSVs per generation in the swine germline is conservative, due to our small sample size and restricted possibilities of dnSV detection from short-read sequencing. The current study highlights the complexity of dnSVs and shows the potential of breeding programs for pigs and livestock species in general, to provide a suitable population structure for identification and characterisation of dnSVs.

Bayesian QTL analyses using pedigreed families of an outcrossing species, with application to fruit firmness in apple

Proof of concept of Bayesian integrated QTL analyses across pedigree-related families from breeding programs of an outbreeding species. Results include QTL confidence intervals, individuals' genotype probabilities and genomic breeding values. Bayesian QTL linkage mapping approaches offer the flexibility to study multiple full sib families with known pedigrees simultaneously. Such a joint analysis increases the probability of detecting these quantitative trait loci (QTL) and provide insight of the magnitude of QTL across different genetic backgrounds. Here, we present an improved Bayesian multi-QTL pedigree-based approach on an outcrossing species using progenies with different (complex) genetic relationships. Different modeling assumptions were studied in the QTL analyses, i.e., the a priori expected number of QTL varied and polygenic effects were considered. The inferences include number of QTL, additive QTL effect sizes and supporting credible intervals, posterior probabilities of QTL genotypes for all individuals in the dataset, and QTL-based as well as genome-wide breeding values. All these features have been implemented in the FlexQTL(™) software. We analyzed fruit firmness in a large apple dataset that comprised 1,347 individuals forming 27 full sib families and their known ancestral pedigrees, with genotypes for 87 SSR markers on 17 chromosomes. We report strong or positive evidence for 14 QTL for fruit firmness on eight chromosomes, validating our approach as several of these QTL were reported previously, though dispersed over a series of studies based on single mapping populations. Interpretation of linked QTL was possible via individuals' QTL genotypes. The correlation between the genomic breeding values and phenotypes was on average 90 %, but varied with the number of detected QTL in a family. The detailed posterior knowledge on QTL of potential parents is critical for the efficiency of marker-assisted breeding.

Multi-trait and multi-environment QTL analyses of yield and a set of physiological traits in pepper

A mixed model framework was defined for QTL analysis of multiple traits across multiple environments for a RIL population in pepper. Detection power for QTLs increased considerably and detailed study of QTL by environment interactions and pleiotropy was facilitated. For many agronomic crops, yield is measured simultaneously with other traits across multiple environments. The study of yield can benefit from joint analysis with other traits and relations between yield and other traits can be exploited to develop indirect selection strategies. We compare the performance of three multi-response QTL approaches based on mixed models: a multi-trait approach (MT), a multi-environment approach (ME), and a multi-trait multi-environment approach (MTME). The data come from a multi-environment experiment in pepper, for which 15 traits were measured in four environments. The approaches were compared in terms of number of QTLs detected for each trait, the explained variance, and the accuracy of prediction for the final QTL model. For the four environments together, the superior MTME approach delivered a total of 47 regions containing putative QTLs. Many of these QTLs were pleiotropic and showed quantitative QTL by environment interaction. MTME was superior to ME and MT in the number of QTLs, the explained variance and accuracy of predictions. The large number of model parameters in the MTME approach was challenging and we propose several guidelines to help obtain a stable final QTL model. The results confirmed the feasibility and strengths of novel mixed model QTL methodology to study the architecture of complex traits.

Multi-trait and multi-environment QTL analyses of yield and a set of physiological traits in pepper

A mixed model framework was defined for QTL analysis of multiple traits across multiple environments for a RIL population in pepper. Detection power for QTLs increased considerably and detailed study of QTL by environment interactions and pleiotropy was facilitated. For many agronomic crops, yield is measured simultaneously with other traits across multiple environments. The study of yield can benefit from joint analysis with other traits and relations between yield and other traits can be exploited to develop indirect selection strategies. We compare the performance of three multi-response QTL approaches based on mixed models: a multi-trait approach (MT), a multi-environment approach (ME), and a multi-trait multi-environment approach (MTME). The data come from a multi-environment experiment in pepper, for which 15 traits were measured in four environments. The approaches were compared in terms of number of QTLs detected for each trait, the explained variance, and the accuracy of prediction for the final QTL model. For the four environments together, the superior MTME approach delivered a total of 47 regions containing putative QTLs. Many of these QTLs were pleiotropic and showed quantitative QTL by environment interaction. MTME was superior to ME and MT in the number of QTLs, the explained variance and accuracy of predictions. The large number of model parameters in the MTME approach was challenging and we propose several guidelines to help obtain a stable final QTL model. The results confirmed the feasibility and strengths of novel mixed model QTL methodology to study the architecture of complex traits.

Bayesian reanalysis of a quantitative trait locus accounting for multiple environments by scaling in broilers1

A Bayesian method was developed to handle QTL analyses of multiple experimental data of outbred populations with heterogeneity of variance between sexes for all random effects. The method employed a scaled reduced animal model with random polygenic and QTL allelic effects. A parsimonious model specification was applied by choosing assumptions regarding the covariance structure to limit the number of parameters to estimate. Markov chain Monte Carlo algorithms were applied to obtain marginal posterior densities. Simulation demonstrated that joint analysis of multiple environments is more powerful than separate single trait analyses of each environment. Measurements on broiler BW obtained from 2 experiments concerning growth efficiency and carcass traits were used to illustrate the method. The population consisted of 10 full-sib families from a cross between 2 broiler lines. Microsatellite genotypes were determined on generations 1 and 2, and phenotypes were collected on groups of generation 3 animals. The model included a polygenic correlation, which had a posterior mean of 0.70 in the analyses. The reanalysis agreed on the presence of a QTL in marker bracket MCW0058-LEI0071 accounting for 34% of the genetic variation in males and 24% in females in the growth efficiency experiment. In the carcass experiment, this QTL accounted for 19% of the genetic variation in males and 6% in females.

Whole genome scan to detect chromosomal regions affecting multiple traits in dairy cattle

Chromosomal regions affecting multiple traits (multiple trait quantitative trait regions or MQR) in dairy cattle were detected using a method based on results from single trait analyses to detect quantitative trait loci (QTL). The covariance between contrasts for different traits in single trait regression analysis was computed. A chromosomal region was considered an MQR when the observed covariance between contrasts deviated from the expected covariance under the null hypothesis of no pleiotropy or close linkage. The expected covariance and the confidence interval for the expected covariance were determined by permutation of the data. Four categories of traits were analyzed: production (5 traits), udder conformation (6 traits), udder health (2 traits), and fertility (2 traits). The analysis of a granddaughter design involving 833 sons of 20 grandsires resulted in 59 MQR (alpha = 0.01, chromosomewise). Fifteen MQR were found on Bos taurus autosome (BTA) 14. Four or more MQR were found on BTA 6, 13, 19, 22, 23, and 25. Eight MQR involving udder conformation and udder health and 4 MQR involving production traits and udder health were found. Five MQR were identified for combinations of fertility and udder conformation traits, and another 5 MQR were identified for combinations of fertility and production traits. For 22 MQR, the difference between the correlation attributable to the MQR and the overall genetic correlation was >0.60. Although the false discovery rate was relatively high (0.52), it was considered important to present these results to assess potential consequences of using these MQR for marker-assisted selection.

Scaling to account for heterogeneous variances in a Bayesian analysis of broiler quantitative trait loci

A Bayesian method for QTL analysis that is capable of accounting for heterogeneity of variance between sexes, is introduced. The Bayesian method uses a parsimonious model that includes scaling parameters for polygenic and QTL allelic effects per sex. Furthermore, the method employs a reduced animal model to increase computational efficiency. Markov Chain Monte Carlo techniques were applied to obtain estimates of genetic parameters. In comparison with previous regression analyses, the Bayesian method 1) estimates dispersion parameters and polygenic effects, 2) uses individual observations instead of offspring averages, and 3) estimates fixed effect levels and covariates and heterogeneity of variance between sexes simultaneously with other parameters, taking uncertainties fully into account. Broiler data collected in a feed efficiency and a carcass experiment were used to illustrate QTL analysis based on the Bayesian method. The experiments were conducted in a population consisting of 10 full-sib families of a cross between two broiler lines. Microsatellite genotypes were determined on generation 1 and 2 animals and phenotypes were collected on third-generation offspring from mating members from different families. Chromosomal regions that seemed to contain a QTL in previous regression analyses and showed heterogeneity of variance were chosen. Traits analyzed in the feed efficiency experiment were BW at 48 d and growth, feed intake, and feed intake corrected for BW between 23 and 48 d. In the carcass experiment, carcass percentage was analyzed. The Bayesian method was successful in finding QTL in all regions previously detected.