Genome-wide association studies for additive and dominance effects for body composition traits in commercial crossbred Piétrain pigs.

Fat depth (FD) and muscle depth (MD) are economically important traits and used to estimate carcass lean content (LMP), which is one of the main breeding objectives in pig breeding programmes. We assessed the genetic architectures of body composition traits for additive and dominance effects in commercial crossbred Piétrain pigs using both 50 K array and sequence genotypes. We first performed a genome-wide association study (GWAS) using single-marker association analysis with a false discovery rate of 0.1. Then, we estimated the additive and dominance effects of the most significant variant in the quantitative trait loci (QTL) regions. It was investigated whether the use of whole-genome sequence (WGS) will improve the QTL detection (both additive and dominance) with a higher power compared with lower density SNP arrays. Our results showed that more QTL regions were detected by WGS compared with 50 K array (n = 54 vs. n = 17). Of the novel associated regions associated with FD and LMP and detected by WGS, the most pronounced peak was on SSC13, situated at ~116-118, 121-127 and 129-134 Mbp. Additionally, we found that only additive effects contributed to the genetic architecture of the analysed traits and no significant dominance effects were found for the tested SNPs at QTL regions, regardless of panel density. The associated SNPs are located in or near several relevant candidate genes. Of these genes, GABRR2, GALR1, RNGTT, CDH20 and MC4R have been previously reported as being associated with fat deposition traits. However, the genes on SSC1 (ZNF292, ORC3, CNR1, SRSF12, MDN1, TSHZ1, RELCH and RNF152) and SSC18 (TTC26 and KIAA1549) have not been reported previously to our best knowledge. Our current findings provide insights into the genomic regions influencing composition traits in Piétrain pigs.

Demography and mating system shape the genome-wide impact of purifying selection in Arabis alpina.

Plant mating systems have profound effects on levels and structuring of genetic variation and can affect the impact of natural selection. Although theory predicts that intermediate outcrossing rates may allow plants to prevent accumulation of deleterious alleles, few studies have empirically tested this prediction using genomic data. Here, we study the effect of mating system on purifying selection by conducting population-genomic analyses on whole-genome resequencing data from 38 European individuals of the arctic-alpine crucifer Arabis alpina We find that outcrossing and mixed-mating populations maintain genetic diversity at similar levels, whereas highly self-fertilizing Scandinavian A. alpina show a strong reduction in genetic diversity, most likely as a result of a postglacial colonization bottleneck. We further find evidence for accumulation of genetic load in highly self-fertilizing populations, whereas the genome-wide impact of purifying selection does not differ greatly between mixed-mating and outcrossing populations. Our results demonstrate that intermediate levels of outcrossing may allow efficient selection against harmful alleles, whereas demographic effects can be important for relaxed purifying selection in highly selfing populations. Thus, mating system and demography shape the impact of purifying selection on genomic variation in A. alpina These results are important for an improved understanding of the evolutionary consequences of mating system variation and the maintenance of mixed-mating strategies.

A survey of functional genomic variation in domesticated chickens.

BACKGROUND: Deleterious genetic variation can increase in frequency as a result of mutations, genetic drift, and genetic hitchhiking. Although individual effects are often small, the cumulative effect of deleterious genetic variation can impact population fitness substantially. In this study, we examined the genome of commercial purebred chicken lines for deleterious and functional variations, combining genotype and whole-genome sequence data. RESULTS: We analysed over 22,000 animals that were genotyped on a 60 K SNP chip from four purebred lines (two white egg and two brown egg layer lines) and two crossbred lines. We identified 79 haplotypes that showed a significant deficit in homozygous carriers. This deficit was assumed to stem from haplotypes that potentially harbour lethal recessive variations. To identify potentially deleterious mutations, a catalogue of over 10 million variants was derived from 250 whole-genome sequenced animals from three purebred white-egg layer lines. Out of 4219 putative deleterious variants, 152 mutations were identified that likely induce embryonic lethality in the homozygous state. Inferred deleterious variation showed evidence of purifying selection and deleterious alleles were generally overrepresented in regions of low recombination. Finally, we found evidence that mutations, which were inferred to be evolutionally intolerant, likely have positive effects in commercial chicken populations. CONCLUSIONS: We present a comprehensive genomic perspective on deleterious and functional genetic variation in egg layer breeding lines, which are under intensive selection and characterized by a small effective population size. We show that deleterious variation is subject to purifying selection and that there is a positive relationship between recombination rate and purging efficiency. In addition, multiple putative functional coding variants were discovered in selective sweep regions, which are likely under positive selection. Together, this study provides a unique molecular perspective on functional and deleterious variation in commercial egg-laying chickens, which can enhance current genomic breeding practices to lower the frequency of undesirable variants in the population.

Deciphering the genetic control of fruit texture in apple by multiple family-based analysis and genome-wide association.

Fruit texture is a complex feature composed of mechanical and acoustic properties relying on the modifications occurring in the cell wall throughout fruit development and ripening. Apple is characterized by a large variation in fruit texture behavior that directly impacts both the consumer's appreciation and post-harvest performance. To decipher the genetic control of fruit texture comprehensively, two complementing quantitative trait locus (QTL) mapping approaches were employed. The first was represented by a pedigree-based analysis (PBA) carried out on six full-sib pedigreed families, while the second was a genome-wide association study (GWAS) performed on a collection of 233 apple accessions. Both plant materials were genotyped with a 20K single nucleotide polymorphism (SNP) array and phenotyped with a sophisticated high-resolution texture analyzer. The overall QTL results indicated the fundamental role of chromosome 10 in controlling the mechanical properties, while chromosomes 2 and 14 were more associated with the acoustic response. The latter QTL, moreover, showed a consistent relationship between the QTL-estimated genotypes and the acoustic performance assessed among seedlings. The in silico annotation of these intervals revealed interesting candidate genes potentially involved in fruit texture regulation, as suggested by the gene expression profile. The joint integration of these approaches sheds light on the specific control of fruit texture, enabling important genetic information to assist in the selection of valuable fruit quality apple varieties.

Detecting QTLs and putative candidate genes involved in budbreak and flowering time in an apple multiparental population.

In temperate trees, growth resumption in spring time results from chilling and heat requirements, and is an adaptive trait under global warming. Here, the genetic determinism of budbreak and flowering time was deciphered using five related full-sib apple families. Both traits were observed over 3 years and two sites and expressed in calendar and degree-days. Best linear unbiased predictors of genotypic effect or interaction with climatic year were extracted from mixed linear models and used for quantitative trait locus (QTL) mapping, performed with an integrated genetic map containing 6849 single nucleotide polymorphisms (SNPs), grouped into haplotypes, and with a Bayesian pedigree-based analysis. Four major regions, on linkage group (LG) 7, LG10, LG12, and LG9, the latter being the most stable across families, sites, and years, explained 5.6-21.3% of trait variance. Co-localizations for traits in calendar days or growing degree hours (GDH) suggested common genetic determinism for chilling and heating requirements. Homologs of two major flowering genes, AGL24 and FT, were predicted close to LG9 and LG12 QTLs, respectively, whereas Dormancy Associated MADs-box (DAM) genes were near additional QTLs on LG8 and LG15. This suggests that chilling perception mechanisms could be common among perennial and annual plants. Progenitors with favorable alleles depending on trait and LG were identified and could benefit new breeding strategies for apple adaptation to temperature increase.

FaRXf1: a locus conferring resistance to angular leaf spot caused by Xanthomonas fragariae in octoploid strawberry.

KEY MESSAGE: Angular leaf spot is a devastating bacterial disease of strawberry. Resistance from two wild accessions is highly heritable and controlled by a major locus on linkage group 6D. Angular leaf spot caused by Xanthomonas fragariae is the only major bacterial disease of cultivated strawberry (Fragaria ×ananassa). While this disease may cause reductions of up to 8 % of marketable yield in Florida winter annual production, no resistant cultivars have been commercialized. Wild accessions US4808 and US4809 were previously identified as resistant to the four genetic clades of X. fragariae, and introgression of the trait into commercial quality perennial-type germplasm was initiated. Previous reports indicated high heritability for the trait but proposed both single-locus and multi-locus inheritance models. The objective of this study was to determine the mode of inheritance of resistance, to identify causal loci, and to begin introgression of resistance into Florida-adapted germplasm. Resistance was observed in two years of field trials with inoculated plants that assayed four full-sib families descended from US4808 to US4809. Resistance segregated 1:1 in all families indicating control by a dominant allele at a single locus. Using a selective genotyping approach with the IStraw90 Axiom(®) SNP array and pedigree-based QTL detection, a single major-effect QTL was identified in two full-sib families, one descended from each resistant accession. High-resolution melt curve analysis validated the presence of the QTL in separate populations. The QTL was delimited to the 33.1-33.6 Mbp (F. vesca vesca v1.1 reference) and 34.8-35.3 Mbp (F. vesca bracteata v2.0 reference) regions of linkage group 6D for both resistance sources and was designated FaRXf1. Characterization of this locus will facilitate marker-assisted selection toward the development of resistant cultivars.

Genomic prediction using imputed whole-genome sequence data in Holstein Friesian cattle.

BACKGROUND: In contrast to currently used single nucleotide polymorphism (SNP) panels, the use of whole-genome sequence data is expected to enable the direct estimation of the effects of causal mutations on a given trait. This could lead to higher reliabilities of genomic predictions compared to those based on SNP genotypes. Also, at each generation of selection, recombination events between a SNP and a mutation can cause decay in reliability of genomic predictions based on markers rather than on the causal variants. Our objective was to investigate the use of imputed whole-genome sequence genotypes versus high-density SNP genotypes on (the persistency of) the reliability of genomic predictions using real cattle data. METHODS: Highly accurate phenotypes based on daughter performance and Illumina BovineHD Beadchip genotypes were available for 5503 Holstein Friesian bulls. The BovineHD genotypes (631,428 SNPs) of each bull were used to impute whole-genome sequence genotypes (12,590,056 SNPs) using the Beagle software. Imputation was done using a multi-breed reference panel of 429 sequenced individuals. Genomic estimated breeding values for three traits were predicted using a Bayesian stochastic search variable selection (BSSVS) model and a genome-enabled best linear unbiased prediction model (GBLUP). Reliabilities of predictions were based on 2087 validation bulls, while the other 3416 bulls were used for training. RESULTS: Prediction reliabilities ranged from 0.37 to 0.52. BSSVS performed better than GBLUP in all cases. Reliabilities of genomic predictions were slightly lower with imputed sequence data than with BovineHD chip data. Also, the reliabilities tended to be lower for both sequence data and BovineHD chip data when relationships between training animals were low. No increase in persistency of prediction reliability using imputed sequence data was observed. CONCLUSIONS: Compared to BovineHD genotype data, using imputed sequence data for genomic prediction produced no advantage. To investigate the putative advantage of genomic prediction using (imputed) sequence data, a training set with a larger number of individuals that are distantly related to each other and genomic prediction models that incorporate biological information on the SNPs or that apply stricter SNP pre-selection should be considered.

Genome-wide distribution of genetic diversity and linkage disequilibrium in a mass-selected population of maritime pine.

BACKGROUND: The accessibility of high-throughput genotyping technologies has contributed greatly to the development of genomic resources in non-model organisms. High-density genotyping arrays have only recently been developed for some economically important species such as conifers. The potential for using genomic technologies in association mapping and breeding depends largely on the genome wide patterns of diversity and linkage disequilibrium in current breeding populations. This study aims to deepen our knowledge regarding these issues in maritime pine, the first species used for reforestation in south western Europe. RESULTS: Using a new map merging algorithm, we first established a 1,712 cM composite linkage map (comprising 1,838 SNP markers in 12 linkage groups) by bringing together three already available genetic maps. Using rigorous statistical testing based on kernel density estimation and resampling we identified cold and hot spots of recombination. In parallel, 186 unrelated trees of a mass-selected population were genotyped using a 12k-SNP array. A total of 2,600 informative SNPs allowed to describe historical recombination, genetic diversity and genetic structure of this recently domesticated breeding pool that forms the basis of much of the current and future breeding of this species. We observe very low levels of population genetic structure and find no evidence that artificial selection has caused a reduction in genetic diversity. By combining these two pieces of information, we provided the map position of 1,671 SNPs corresponding to 1,192 different loci. This made it possible to analyze the spatial pattern of genetic diversity (He) and long distance linkage disequilibrium (LD) along the chromosomes. We found no particular pattern in the empirical variogram of He across the 12 linkage groups and, as expected for an outcrossing species with large effective population size, we observed an almost complete lack of long distance LD. CONCLUSIONS: These results are a stepping stone for the development of strategies for studies in population genomics, association mapping and genomic prediction in this economical and ecologically important forest tree species.

Longitudinal individual predictions from irregular repeated measurements data.

Intensive longitudinal data can be used to explore important associations and patterns between various types of inputs and outcomes. Nonlinear relations and irregular measurement occasions can pose problems to develop an accurate model for these kinds of data. This paper focuses on the development, fitting and evaluation of a prediction model with irregular intensive longitudinal data. A three-step process for developing a prediction tool for (daily) monitoring and prediction is outlined and illustrated for intensive weight measurements in piglets. Step 1 addresses a nonlinear relation in the data by developing and applying a normalizing transformation. Step 2 addresses the intermittent nature of the time points by aligning the measurement times to a common time grid with a broken-stick model. Step 3 addresses the prediction problem by selecting and evaluating inputs and covariates in the model to obtain accurate predictions. The final model predicts future outcomes accurately, while allowing for nonlinearities between input and output and for different measurement histories of individuals. The methodology described can be used to develop a tool to deal with intensive irregular longitudinal data that uses the available information in an optimal way. The resulting tool demonstrated to perform well for piglet weight prediction and can be adapted to many different applications.

Mutations on a conserved distal enhancer in the porcine C-reactive protein gene impair its expression in liver.

C-reactive protein (CRP) is an evolutionary highly conserved protein. Like humans, CRP acts as a major acute phase protein in pigs. While CRP regulatory mechanisms have been extensively studied in humans, little is known about the molecular mechanisms that control pig CRP gene expression. The main goal of the present work was to study the regulatory mechanisms and identify functional genetic variants regulating CRP gene expression and CRP blood levels in pigs. The characterization of the porcine CRP proximal promoter region revealed a high level of conservation with both cow and human promoters, sharing binding sites for transcription factors required for CRP expression. Through genome-wide association studies and fine mapping, the most associated variants with both mRNA and protein CRP levels were localized in a genomic region 39.3 kb upstream of CRP. Further study of the region revealed a highly conserved putative enhancer that contains binding sites for several transcriptional regulators such as STAT3, NF-kB or C/EBP-ß. Luciferase reporter assays showed the necessity of this enhancer-promoter interaction for the acute phase induction of CRP expression in liver, where differences in the enhancer sequences significantly modified CRP activity. The associated polymorphisms disrupted the putative binding sites for HNF4α and FOXA2 transcription factors. The high correlation between HNF4α and CRP expression levels suggest the participation of HNF4α in the regulatory mechanism of porcine CRP expression through the modification of its binding site in liver. Our findings determine, for the first time, the relevance of a distal regulatory element essential for the acute phase induction of porcine CRP in liver and identify functional polymorphisms that can be included in pig breeding programs to improve immunocompetence.

QTL linkage analysis of connected populations using ancestral marker and pedigree information.

The common assumption in quantitative trait locus (QTL) linkage mapping studies that parents of multiple connected populations are unrelated is unrealistic for many plant breeding programs. We remove this assumption and propose a Bayesian approach that clusters the alleles of the parents of the current mapping populations from locus-specific identity by descent (IBD) matrices that capture ancestral marker and pedigree information. Moreover, we demonstrate how the parental IBD data can be incorporated into a QTL linkage analysis framework by using two approaches: a Threshold IBD model (TIBD) and a Latent Ancestral Allele Model (LAAM). The TIBD and LAAM models are empirically tested via numerical simulation based on the structure of a commercial maize breeding program. The simulations included a pilot dataset with closely linked QTL on a single linkage group and 100 replicated datasets with five linkage groups harboring four unlinked QTL. The simulation results show that including parental IBD data (similarly for TIBD and LAAM) significantly improves the power and particularly accuracy of QTL mapping, e.g., position, effect size and individuals' genotype probability without significantly increasing computational demand.

Pedimap: software for the visualization of genetic and phenotypic data in pedigrees.

Pedimap is a user-friendly software tool for visualizing phenotypic and genotypic data for related individuals linked in pedigrees. Genetic data can include marker scores, Identity-by-Descent probabilities, and marker linkage map positions, allowing the visualization of haplotypes through lineages. The pedigrees can accommodate all types of inheritance, including selfing, cloning, and repeated backcrossing, and all ploidy levels are supported. Visual association of the genetic data with phenotypic data simplifies the exploration of large data sets, thereby improving breeding decision making. Data are imported from text files; in addition data exchange with other software packages (FlexQTL(TM) and GenomeStudio(TM)) is possible. Instructions for use and an executable version compatible with the Windows platform are available for free from http://www.plantbreeding.wur.nl/UK/software_pedimap.html.

Prediction performance of linear models and gradient boosting machine on complex phenotypes in outbred mice.

We compared the performance of linear (GBLUP, BayesB, and elastic net) methods to a nonparametric tree-based ensemble (gradient boosting machine) method for genomic prediction of complex traits in mice. The dataset used contained genotypes for 50,112 SNP markers and phenotypes for 835 animals from 6 generations. Traits analyzed were bone mineral density, body weight at 10, 15, and 20 weeks, fat percentage, circulating cholesterol, glucose, insulin, triglycerides, and urine creatinine. The youngest generation was used as a validation subset, and predictions were based on all older generations. Model performance was evaluated by comparing predictions for animals in the validation subset against their adjusted phenotypes. Linear models outperformed gradient boosting machine for 7 out of 10 traits. For bone mineral density, cholesterol, and glucose, the gradient boosting machine model showed better prediction accuracy and lower relative root mean squared error than the linear models. Interestingly, for these 3 traits, there is evidence of a relevant portion of phenotypic variance being explained by epistatic effects. Using a subset of top markers selected from a gradient boosting machine model helped for some of the traits to improve the accuracy of prediction when these were fitted into linear and gradient boosting machine models. Our results indicate that gradient boosting machine is more strongly affected by data size and decreased connectedness between reference and validation sets than the linear models. Although the linear models outperformed gradient boosting machine for the polygenic traits, our results suggest that gradient boosting machine is a competitive method to predict complex traits with assumed epistatic effects.

Adding gene transcripts into genomic prediction improves accuracy and reveals sampling time dependence.

Recent developments allowed generating multiple high-quality 'omics' data that could increase the predictive performance of genomic prediction for phenotypes and genetic merit in animals and plants. Here, we have assessed the performance of parametric and nonparametric models that leverage transcriptomics in genomic prediction for 13 complex traits recorded in 478 animals from an outbred mouse population. Parametric models were implemented using the best linear unbiased prediction, while nonparametric models were implemented using the gradient boosting machine algorithm. We also propose a new model named GTCBLUP that aims to remove between-omics-layer covariance from predictors, whereas its counterpart GTBLUP does not do that. While gradient boosting machine models captured more phenotypic variation, their predictive performance did not exceed the best linear unbiased prediction models for most traits. Models leveraging gene transcripts captured higher proportions of the phenotypic variance for almost all traits when these were measured closer to the moment of measuring gene transcripts in the liver. In most cases, the combination of layers was not able to outperform the best single-omics models to predict phenotypes. Using only gene transcripts, the gradient boosting machine model was able to outperform best linear unbiased prediction for most traits except body weight, but the same pattern was not observed when using both single nucleotide polymorphism genotypes and gene transcripts. Although the GTCBLUP model was not able to produce the most accurate phenotypic predictions, it showed the highest accuracies for breeding values for 9 out of 13 traits. We recommend using the GTBLUP model for prediction of phenotypes and using the GTCBLUP for prediction of breeding values.

Imputation to whole-genome sequence and its use in genome-wide association studies for pork colour traits in crossbred and purebred pigs.

Imputed whole-genome sequence (WGS) has been proposed to improve genome-wide association studies (GWAS), since all causative mutations responsible for phenotypic variation are expected to be present in the data. This approach was applied on a large number of purebred (PB) and crossbred (CB) pigs for 18 pork color traits to evaluate the impact of using imputed WGS relative to medium-density marker panels. The traits included Minolta A*, B*, and L* for fat (FCOL), quadriceps femoris muscle (QFCOL), thawed loin muscle (TMCOL), fresh ham gluteus medius (GMCOL), ham iliopsoas muscle (ICOL), and longissimus dorsi muscle on the fresh loin (FMCOL). Sequence variants were imputed from a medium-density marker panel (61K for CBs and 50K for PBs) in all genotyped pigs using BeagleV5.0. We obtained high imputation accuracy (average of 0.97 for PBs and 0.91 for CBs). GWAS were conducted for three datasets: 954 CBs and 891 PBs, and the combined CBs and PBs. For most traits, no significant associations were detected, regardless of panel density or population type. However, quantitative trait loci (QTL) regions were only found for a few traits including TMCOL Minolta A* and GMCOL Minolta B* (CBs), FMCOL Minolta B*, FMCOL Minolta L*, and ICOL Minolta B* (PBs) and FMCOL Minolta A*, FMCOL Minolta B*, GMCOL Minolta B*, and ICOL Minolta B* (Combined dataset). More QTL regions were identified with WGS (n = 58) relative to medium-density marker panels (n = 22). Most of the QTL were linked to previously reported QTLs or candidate genes that have been previously reported to be associated with meat quality, pH and pork color; e.g., VIL1, PRKAG3, TTLL4, and SLC11A1, USP37. CTDSP1 gene on SSC15 has not been previously associated with meat color traits in pigs. The findings suggest any added value of WGS was only for detecting novel QTL regions when the sample size is sufficiently large as with the Combined dataset in this study. The percentage of phenotypic variance explained by the most significant SNPs also increased with WGS compared with medium-density panels. The results provide additional insights into identification of a number of candidate regions and genes for pork color traits in different pig populations.

A new haplotype-resolved turkey genome to enable turkey genetics and genomics research.

BACKGROUND: The domesticated turkey (Meleagris gallopavo) is a species of significant agricultural importance and is the second largest contributor, behind broiler chickens, to world poultry meat production. The previous genome is of draft quality and partly based on the chicken (Gallus gallus) genome. A high-quality reference genome of M. gallopavo is essential for turkey genomics and genetics research and the breeding industry. RESULTS: By adopting the trio-binning approach, we were able to assemble a high-quality chromosome-level F1 assembly and 2 parental haplotype assemblies, leveraging long-read technologies and genome-wide chromatin interaction data (Hi-C). From a total of 40 chromosomes (2n = 80), we captured 35 chromosomes in a single scaffold, showing much improved genome completeness and continuity compared to the old assembly build. The 3 assemblies are of higher quality than the previous draft quality assembly and comparable to the chicken assemblies (GRCg7) shown by the largest contig N50 (26.6 Mb) and comparable BUSCO gene set completeness scores (96-97%). Comparative analyses confirm a previously identified large inversion of around 19 Mbp on the Z chromosome not found in other Galliformes. Structural variation between the parent haplotypes was identified, which poses potential new target genes for breeding. CONCLUSIONS: We contribute a new high-quality turkey genome at the chromosome level, benefiting turkey genetics and other avian genomics research as well as the turkey breeding industry.

Identification of transcriptional regulatory variants in pig duodenum, liver, and muscle tissues.

BACKGROUND: In humans and livestock species, genome-wide association studies (GWAS) have been applied to study the association between variants distributed across the genome and a phenotype of interest. To discover genetic polymorphisms affecting the duodenum, liver, and muscle transcriptomes of 300 pigs from 3 different breeds (Duroc, Landrace, and Large White), we performed expression GWAS between 25,315,878 polymorphisms and the expression of 13,891 genes in duodenum, 12,748 genes in liver, and 11,617 genes in muscle. RESULTS: More than 9.68 × 1011 association tests were performed, yielding 14,096,080 significantly associated variants, which were grouped in 26,414 expression quantitative trait locus (eQTL) regions. Over 56% of the variants were within 1 Mb of their associated gene. In addition to the 100-kb region upstream of the transcription start site, we identified the importance of the 100-kb region downstream of the 3'UTR for gene regulation, as most of the cis-regulatory variants were located within these 2 regions. We also observed 39,874 hotspot regulatory polymorphisms associated with the expression of 10 or more genes that could modify the protein structure or the expression of a regulator gene. In addition, 2 motifs (5'-GATCCNGYGTTGCYG-3' and a poly(A) sequence) were enriched across the 3 tissues within the neighboring sequences of the most significant single-nucleotide polymorphisms in each cis-eQTL region. CONCLUSIONS: The 14 million significant associations obtained in this study are publicly available and have enabled the identification of expression-associated cis-, trans-, and hotspot regulatory variants within and across tissues, thus shedding light on the molecular mechanisms of regulatory variations that shape end-trait phenotypes.

High-quality, genome-wide SNP genotypic data for pedigreed germplasm of the diploid outbreeding species apple, peach, and sweet cherry through a common workflow.

High-quality genotypic data is a requirement for many genetic analyses. For any crop, errors in genotype calls, phasing of markers, linkage maps, pedigree records, and unnoticed variation in ploidy levels can lead to spurious marker-locus-trait associations and incorrect origin assignment of alleles to individuals. High-throughput genotyping requires automated scoring, as manual inspection of thousands of scored loci is too time-consuming. However, automated SNP scoring can result in errors that should be corrected to ensure recorded genotypic data are accurate and thereby ensure confidence in downstream genetic analyses. To enable quick identification of errors in a large genotypic data set, we have developed a comprehensive workflow. This multiple-step workflow is based on inheritance principles and on removal of markers and individuals that do not follow these principles, as demonstrated here for apple, peach, and sweet cherry. Genotypic data was obtained on pedigreed germplasm using 6-9K SNP arrays for each crop and a subset of well-performing SNPs was created using ASSIsT. Use of correct (and corrected) pedigree records readily identified violations of simple inheritance principles in the genotypic data, streamlined with FlexQTL software. Retained SNPs were grouped into haploblocks to increase the information content of single alleles and reduce computational power needed in downstream genetic analyses. Haploblock borders were defined by recombination locations detected in ancestral generations of cultivars and selections. Another round of inheritance-checking was conducted, for haploblock alleles (i.e., haplotypes). High-quality genotypic data sets were created using this workflow for pedigreed collections representing the U.S. breeding germplasm of apple, peach, and sweet cherry evaluated within the RosBREED project. These data sets contain 3855, 4005, and 1617 SNPs spread over 932, 103, and 196 haploblocks in apple, peach, and sweet cherry, respectively. The highly curated phased SNP and haplotype data sets, as well as the raw iScan data, of germplasm in the apple, peach, and sweet cherry Crop Reference Sets is available through the Genome Database for Rosaceae.

Predicting Flowering Behavior and Exploring Its Genetic Determinism in an Apple Multi-family Population Based on Statistical Indices and Simplified Phenotyping.

Irregular flowering over years is commonly observed in fruit trees. The early prediction of tree behavior is highly desirable in breeding programmes. This study aims at performing such predictions, combining simplified phenotyping and statistics methods. Sequences of vegetative vs. floral annual shoots (AS) were observed along axes in trees belonging to five apple related full-sib families. Sequences were analyzed using Markovian and linear mixed models including year and site effects. Indices of flowering irregularity, periodicity and synchronicity were estimated, at tree and axis scales. They were used to predict tree behavior and detect QTL with a Bayesian pedigree-based analysis, using an integrated genetic map containing 6,849 SNPs. The combination of a Biennial Bearing Index (BBI) with an autoregressive coefficient (γ g ) efficiently predicted and classified the genotype behaviors, despite few misclassifications. Four QTLs common to BBIs and γ g and one for synchronicity were highlighted and revealed the complex genetic architecture of the traits. Irregularity resulted from high AS synchronism, whereas regularity resulted from either asynchronous locally alternating or continual regular AS flowering. A relevant and time-saving method, based on a posteriori sampling of axes and statistical indices is proposed, which is efficient to evaluate the tree breeding values for flowering regularity and could be transferred to other species.

Genome-Wide Association Mapping of Flowering and Ripening Periods in Apple.

Deciphering the genetic control of flowering and ripening periods in apple is essential for breeding cultivars adapted to their growing environments. We implemented a large Genome-Wide Association Study (GWAS) at the European level using an association panel of 1,168 different apple genotypes distributed over six locations and phenotyped for these phenological traits. The panel was genotyped at a high-density of SNPs using the Axiom®Apple 480 K SNP array. We ran GWAS with a multi-locus mixed model (MLMM), which handles the putatively confounding effect of significant SNPs elsewhere on the genome. Genomic regions were further investigated to reveal candidate genes responsible for the phenotypic variation. At the whole population level, GWAS retained two SNPs as cofactors on chromosome 9 for flowering period, and six for ripening period (four on chromosome 3, one on chromosome 10 and one on chromosome 16) which, together accounted for 8.9 and 17.2% of the phenotypic variance, respectively. For both traits, SNPs in weak linkage disequilibrium were detected nearby, thus suggesting the existence of allelic heterogeneity. The geographic origins and relationships of apple cultivars accounted for large parts of the phenotypic variation. Variation in genotypic frequency of the SNPs associated with the two traits was connected to the geographic origin of the genotypes (grouped as North+East, West and South Europe), and indicated differential selection in different growing environments. Genes encoding transcription factors containing either NAC or MADS domains were identified as major candidates within the small confidence intervals computed for the associated genomic regions. A strong microsynteny between apple and peach was revealed in all the four confidence interval regions. This study shows how association genetics can unravel the genetic control of important horticultural traits in apple, as well as reduce the confidence intervals of the associated regions identified by linkage mapping approaches. Our findings can be used for the improvement of apple through marker-assisted breeding strategies that take advantage of the accumulating additive effects of the identified SNPs.