Two quantitative trait loci are associated with recapping of Varroa destructor-infested brood cells in Apis mellifera mellifera.

Recapping of Varroa destructor-infested brood cells is a trait that has recently attracted interest in honey bee breeding to select mite-resistant Apis mellifera colonies. To investigate the genetic architecture of this trait, we evaluated a sample of A. mellifera mellifera colonies (N = 155) from Switzerland and France and performed a genome-wide association study, using a pool of 500 workers per colony for next-generation sequencing. The results revealed that two QTL were significantly (P < 0.05) associated with recapping of V. destructor-infested brood cells. The best-associated QTL is located on chromosome 5 in a region previously found to be associated with grooming behaviour, a resistance trait against V. destructor, in A. mellifera and Apis cerana. The second best-associated QTL is located on chromosome 4 in an intron of the Dscam gene, which is involved in neuronal wiring. Previous research demonstrated that genes involved in neuronal wiring are associated with recapping and varroa sensitive hygiene. Therefore, our study confirms the role of a gene region on chromosome 5 in social immunity and simultaneously provides novel insights into genetic interactions between common mite resistance traits in honey bees.

Identification of quantitative trait loci associated with calmness and gentleness in honey bees using whole-genome sequences.

The identification of quantitative trait loci (QTL) through genome-wide association studies (GWAS) is a powerful method for unravelling the genetic background of selected traits and improving early-stage predictions. In honey bees (Apis mellifera), past genetic analyses have particularly focused on individual queens and workers. In this study, we used pooled whole-genome sequences to ascertain the genetic variation of the entire colony. In total, we sampled 216 Apis mellifera mellifera and 28 Apis mellifera carnica colonies. Different experts subjectively assessed the gentleness and calmness of the colonies using a standardised protocol. Conducting a GWAS for calmness on 211 purebred A. m. mellifera colonies, we identified three QTL, on chromosomes 8, 6, and 12. The two first QTL correspond to LOC409692 gene, coding for a disintegrin and metalloproteinase domain-containing protein 10, and to Abscam gene, coding for a Dscam family member Abscam protein, respectively. The last gene has been reported to be involved in the domestication of A. mellifera. The third QTL is located 13 kb upstream of LOC102655631, coding for a trehalose transporter. For gentleness, two QTL were identified on chromosomes 4 and 3. They are located within gene LOC413669, coding for a lap4 protein, and gene LOC413416, coding for a bicaudal C homolog 1-B protein, respectively. The identified positional candidate genes of both traits mainly affect the olfaction and nervous system of honey bees. Further research is needed to confirm the results and to better understand the genetic and phenotypic basis of calmness and gentleness.

Analysis of ROH patterns in the Noriker horse breed reveals signatures of selection for coat color and body size.

Overlapping runs of homozygosity (ROH islands) shared by the majority of a population are hypothesized to be the result of selection around a target locus. In this study we investigated the impact of selection for coat color within the Noriker horse on autozygosity and ROH patterns. We analyzed overlapping homozygous regions (ROH islands) for gene content in fragments shared by more than 50% of horses. Long-term assortative mating of chestnut horses and the small effective population size of leopard spotted and tobiano horses resulted in higher mean genome-wide ROH coverage (SROH ) within the range of 237.4-284.2 Mb, whereas for bay, black and roan horses, where rotation mating is commonly applied, lower autozygosity (SROH from 176.4-180.0 Mb) was determined. We identified seven common ROH islands considering all Noriker horses from our dataset. Specific islands were documented for chestnut, leopard spotted, roan and bay horses. The ROH islands contained, among others, genes associated with body size (ZFAT, LASP1 and LCORL/NCAPG), coat color (MC1R in chestnut and the factor PATN1 in leopard spotted horses) and morphogenesis (HOXB cluster in all color strains except leopard spotted horses). This study demonstrates that within a closed population sharing the same founders and ancestors, selection on a single phenotypic trait, in this case coat color, can result in genetic fragmentation affecting levels of autozygosity and distribution of ROH islands and enclosed gene content.

Novel insights into Sabino1 and splashed white coat color patterns in horses.

Within the framework of genome-wide analyses using the novel Axiom® genotyping array, we investigated the distribution of two previously described coat color patterns, namely sabino1 (SBI), associated with the KIT gene (KI16+1037A), and splashed white, associated with the PAX3 gene (ECA6:g.11429753C>T; PAX3C70Y ), including a total of 899 horses originating from eight different breeds (Achal Theke, Purebred Arabian, Partbred Arabian, Anglo-Arabian, Shagya Arabian, Haflinger, Lipizzan and Noriker). Based on the data we collected we were able to demonstrate that, besides Quarter horses, the PAX3C70Y allele is also present in Noriker (seven out of 189) and Lipizzan (three out of 329) horses. The SB1 allele was present in three breeds (Haflinger, 14 out of 98; Noriker, four out of 189; Lipizzan one out of 329). Furthermore, we examined the phenotypes of SB1- and PAX3C70Y -carrier horses for their characteristic white spotting patterns. None of the SB1/sb1-carrier horses met the criteria defining the Sabino1 pattern according to current applied protocols. From 10 heterozygous PAX3C70Y -carrier horses, two had nearly a splashed white phenotype. The results of this large-scale experiment on the genetic association of white spotting patterns in horses underline the influence of gene interactions and population differences on complex traits such as Sabino1 and splashed white.

Asian horses deepen the MSY phylogeny.

Humans have shaped the population history of the horse ever since domestication about 5500 years ago. Comparative analyses of the Y chromosome can illuminate the paternal origin of modern horse breeds. This may also reveal different breeding strategies that led to the formation of extant breeds. Recently, a horse Y-chromosomal phylogeny of modern horses based on 1.46 Mb of the male-specific Y (MSY) was generated. We extended this dataset with 52 samples from five European, two American and seven Asian breeds. As in the previous study, almost all modern European horses fall into a crown group, connected via a few autochthonous Northern European lineages to the outgroup, the Przewalski's Horse. In total, we now distinguish 42 MSY haplotypes determined by 158 variants within domestic horses. Asian horses show much higher diversity than previously found in European breeds. The Asian breeds also introduce a deep split to the phylogeny, preliminarily dated to 5527 ± 872 years. We conclude that the deep splitting Asian Y haplotypes are remnants of a far more diverse ancient horse population, whose haplotypes were lost in other lineages.

Genome-wide scans between two honeybee populations reveal putative signatures of human-mediated selection.

Human-mediated selection has left signatures in the genomes of many domesticated animals, including the European dark honeybee, Apis mellifera mellifera, which has been selected by apiculturists for centuries. Using whole-genome sequence information, we investigated selection signatures in spatially separated honeybee subpopulations (Switzerland, n = 39 and France, n = 17). Three different test statistics were calculated in windows of 2 kb (fixation index, cross-population extended haplotype homozygosity and cross-population composite likelihood ratio) and combined into a recently developed composite selection score. Applying a stringent false discovery rate of 0.01, we identified six significant selective sweeps distributed across five chromosomes covering eight genes. These genes are associated with multiple molecular and biological functions, including regulation of transcription, receptor binding and signal transduction. Of particular interest is a selection signature on chromosome 1, which corresponds to the WNT4 gene, the family of which is conserved across the animal kingdom with a variety of functions. In Drosophila melanogaster, WNT4 alleles have been associated with differential wing, cross vein and abdominal phenotypes. Defining phenotypic characteristics of different Apis mellifera ssp., which are typically used as selection criteria, include colour and wing venation pattern. This signal is therefore likely to be a good candidate for human mediated-selection arising from different applied breeding practices in the two managed populations.

A genome-wide association study for equine recurrent airway obstruction in European Warmblood horses reveals a suggestive new quantitative trait locus on chromosome 13.

Recurrent airway obstruction (RAO), also known as heaves, is an asthma-like respiratory disease. Its development is strongly influenced by environmental risk factors such as sensitization and exposure to moldy hay, straw bedding and stabling indoors. A hereditary component has been documented in previous studies; however, so far no causative genetic variant that influences the risk of developing RAO has been identified. In this study, we revised an existing dataset and selected 384 horses for genotyping on the Affymetrix high-density equine SNP array. We performed an allelic case-control genome-wide association study, which revealed a suggestively significant association on equine chromosome 13 at 32 843 309 bp. This SNP is located in the protein-coding gene TXNDC11, which is possibly involved in the folding process of the multiprotein complexes DUOX1 and DUOX2. In humans, these proteins are known to take part in regulating the production of H2 O2 in the respiratory tract epithelium as well as in MUC5AC mucin expression. Therefore, TXNDC11 may be considered a functional candidate gene, and further research is needed to explore its potential role in RAO-affected horses.

Association between population structure and allele frequencies of the glycogen synthase 1 mutation in the Austrian Noriker draft horse.

The aim of this study was to determine the allele frequency of the glycogen synthase 1 (GYS1) mutation associated with polysaccharide storage myopathy type 1 in the Austrian Noriker horse. Furthermore, we examined the influence of population substructures on the allele distribution. The study was based upon a comprehensive population sample (208 breeding stallions and 309 mares) and a complete cohort of unselected offspring from the year 2014 (1553 foals). The mean proportion of GYS1 carrier animals in the foal cohort was 33%, ranging from 15% to 50% according to population substructures based on coat colours. In 517 mature breeding horses the mutation carrier frequency reached 34%, ranging on a wider scale from 4% to 62% within genetic substructures. We could show that the occurrence of the mutated GYS1 allele is influenced by coat colour; genetic bottlenecks; and assortative, rotating and random mating strategies. Highest GYS1 carrier frequencies were observed in the chestnut sample comprising 50% in foals, 54% in mares and 62% in breeding stallions. The mean inbreeding of homozygous carrier animals reached 4.10%, whereas non-carrier horses were characterized by an inbreeding coefficient of 3.48%. Lowest GYS1 carrier frequencies were observed in the leopard spotted Noriker subpopulation. Here the mean carrier frequency reached 15% in foals, 17% in mares and 4% in stallions and inbreeding decreased from 3.28% in homozygous non-carrier horses to 2.70% in heterozygous horses and 0.94% in homozygous carriers. This study illustrates that lineage breeding and specified mating strategies result in genetic substructures, which affect the frequencies of the GYS1 gene mutation.

Selection signatures in Shetland ponies.

Shetland ponies were selected for numerous traits including small stature, strength, hardiness and longevity. Despite the different selection criteria, Shetland ponies are well known for their small stature. We performed a selection signature analysis including genome-wide SNPs of 75 Shetland ponies and 76 large-sized horses. Based upon this dataset, we identified a selection signature on equine chromosome (ECA) 1 between 103.8 Mb and 108.5 Mb. A total of 33 annotated genes are located within this interval including the IGF1R gene at 104.2 Mb and the ADAMTS17 gene at 105.4 Mb. These two genes are well known to have a major impact on body height in numerous species including humans. Homozygosity mapping in the Shetland ponies identified a region with increased homozygosity between 107.4 Mb and 108.5 Mb. None of the annotated genes in this region have so far been associated with height. Thus, we cannot exclude the possibility that the identified selection signature on ECA1 is associated with some trait other than height, for which Shetland ponies were selected.

Genome-wide association studies based on sequence-derived genotypes reveal new QTL associated with conformation and performance traits in the Franches-Montagnes horse breed.

To identify novel quantitative trait loci (QTL) within horses, we performed genome-wide association studies (GWAS) based on sequence-level genotypes for conformation and performance traits in the Franches-Montagnes (FM) horse breed. Sequence-level genotypes of FM horses were derived by re-sequencing 30 key founders and imputing 50K data of genotyped horses. In total, we included 1077 FM horses genotyped for ~4 million SNPs and their respective de-regressed breeding values of the traits in the analysis. Based on this dataset, we identified a total of 14 QTL associated with 18 conformation traits and one performance trait. Therefore, our results suggest that the application of sequence-derived genotypes increases the power to identify novel QTL which were not identified previously based on 50K SNP chip data.

Genetic diversity analyses reveal first insights into breed-specific selection signatures within Swiss goat breeds.

We used genotype data from the caprine 50k Illumina BeadChip for the assessment of genetic diversity within and between 10 local Swiss goat breeds. Three different cluster methods allowed the goat samples to be assigned to the respective breed groups, whilst the samples of Nera Verzasca and Tessin Grey goats could not be differentiated from each other. The results of the different genetic diversity measures show that Appenzell, Toggenburg, Valais and Booted goats should be prioritized in future conservation activities. Furthermore, we examined runs of homozygosity (ROH) and compared genomic inbreeding coefficients based on ROH (FROH ) with pedigree-based inbreeding coefficients (FPED ). The linear relationship between FROH and FPED was confirmed for goats by including samples from the three main breeds (Saanen, Chamois and Toggenburg goats). FROH appears to be a suitable measure for describing levels of inbreeding in goat breeds with missing pedigree information. Finally, we derived selection signatures between the breeds. We report a total of 384 putative selection signals. The 25 most significant windows contained genes known for traits such as: coat color variation (MITF, KIT, ASIP), growth (IGF2, IGF2R, HRAS, FGFR3) and milk composition (PITX2). Several other putative genes involved in the formation of populations, which might have been selected for adaptation to the alpine environment, are highlighted. The results provide a contemporary background for the management of genetic diversity in local Swiss goat breeds.

New genomic insights into the conformation of Lipizzan horses.

Conformation traits are important selection criteria in equine breeding, as they describe the exterior aspects of the horse (height, joint angles, shape). However, the genetic architecture of conformation is not well understood, as data of these traits mainly consist of subjective evaluation scores. Here, we performed genome-wide association studies on two-dimensional shape data of Lipizzan horses. Based on this data, we identified significant quantitative trait loci (QTL) associated with cresty neck on equine chromosome (ECA)16 within the MAGI1 gene, and with type, hereby differentiating heavy from light horses on ECA5 within the POU2F1 gene. Both genes were previously described to affect growth, muscling and fatty deposits in sheep, cattle and pigs. Furthermore, we pin-pointed another suggestive QTL on ECA21, near the PTGER4 gene, associated with human ankylosing spondylitis, for shape differences in the back and pelvis (roach back vs sway back). Further differences in the shape of the back and abdomen were suggestively associated with the RYR1 gene, involved in core muscle weakness in humans. Therefore, we demonstrated that horse shape space data enhance the genomic investigations of horse conformation.

Genome-wide association mapping of milk production traits in Braunvieh cattle.

A whole-genome association study of milk production traits: milk yield, protein yield, fat yield, protein percentage, and fat percentage, was performed on the population of Braunvieh cattle. Five hundred and fifty-four progeny-tested bulls and 36,219 autosomal single nucleotide polymorphism (SNP) markers on 29 Bos taurus autosomes (BTA) were included in the analysis. A principal component analysis was conducted to adjust for the effect of population stratification in the analyzed data set. For the principal component analysis, genome-wide relationships between individuals were calculated. Three different criteria (Horn's test, Kaiser's criterion, and Jolliffe's criterion) were tested to determine the number of significant principal components. Estimation of putative associations between SNP and milk production traits was carried out using a linear regression model in R software (R Foundation for Statistical Computing, Vienna, Austria). Significant principal components, adjusting for population stratification separately for each criterion and family relationships and genotypes at individual SNP were included as fixed effects in the model. The inflation factor λ and quantile-quantile plots were calculated to compare how the different criteria deal with stratification in our mapping population. Based on the analyses on all of the aforementioned criteria, we can conclude that Jolliffe's criterion deals the best with population stratification. Furthermore, significance thresholds for a given genome-wide false discovery rate of 5% were estimated and used for specific traits. Three of the analyzed traits showed genome-wide significant association with SNP. Two SNP had an effect on milk yield on BTA4, 2 SNP affected fat yield on BTA14 and BTA23, and 1 SNP was associated with fat percent on BTA1. Single nucleotide polymorphisms identified in this study as associated with milk production traits will further contribute to the mapping of corresponding quantitative trait loci and investigation of the genes responsible for polymorphisms in milk production traits in dairy cattle. Described comparison of different criteria for determination of significant principal components can provide important information for similar studies in stratified populations.