Genome‐wide association study for semen traits of the bulls in Chinese Holstein

A genome‐wide association study (GWAS) was performed to identify markers and candidate genes for five semen traits in the Holstein bull population in China. The analyzed dataset consisted of records from 692 bulls from eight bull stations; each bull was genotyped using the Illumina BovineSNP50 BeadChip. Association tests between each trait and the 41 188 informative high‐quality SNPs were achieved with gapit software. In total, 19 suggestive significant SNPs, partly located within the reported QTL regions or within or close to the reported candidate genes, associated with five semen traits were detected. By combining our GWAS results with the biological functions of these genes, eight novel promising candidate genes, including ETNK1, PDE3A, PDGFRB, CSF1R, WT1, DSCAML1, SOD1 and RUNX2, were identified that potentially relate to semen traits. Our findings may provide a basis for further research on the genetic mechanism of semen traits and marker‐assisted selection of such traits in Holstein bulls.     

Breeding soundness and libido examination of Belgian Blue and Holstein Friesian artificial insemination bulls in Belgium and The Netherlands

Data on breeding soundness and libido evaluations in Belgian Blue (BB) bulls are scarce. The present study compared results of breeding soundness and libido evaluations of young BB bulls to young Holstein Friesian (HF) bulls prior to acceptance into an AI program. Breed differences for breeding soundness exist between BB and HF bulls, as 93.7% of the young BB bulls failed the breeding soundness examination (BSE) compared to 59.3% of the HF bulls (P=0.0005). Within the BB breed, differences were present between bulls of different ages, and bull selection for better fertility with increasing age apparently influenced the results. The number of reasons for which bulls failed the test differed between the age groups in the BB breed, whereas a tendency for more failure reasons in the BB breed was noticed in the breed comparison. The most important reasons for failure were sperm morphology and scrotal circumference (SC), but far more BB bulls failed for these traits compared to the HF breed (82.8% versus 56.0% and 43.8% versus 17.6% in the BB and the HF breed for sperm morphology (P=0.0005) and SC (P<0.0001), respectively). The high proportion of BB bulls with a substandard SC and poor sperm morphology might suggest an increased prevalence of testicular hypoplasia or degeneration within this breed. Concerning libido, the reaction time did not differ either between breeds or between age groups within the BB breed, whereas mounting enthusiasm, although not different between the two breeds, did decline with increasing age, probably due to the greater mating experience of the older bulls. All in all, libido did not seem to be different between the breeds.     

Y chromosome haplotype diversity of domestic sheep (Ovis aries) in northern Eurasia

Variation in two SNPs and one microsatellite on the Y chromosome was analyzed in a total of 663 rams representing 59 breeds from a large geographic range in northern Eurasia. SNPA‐oY1 showed the highest allele frequency (91.55%) across the breeds, whereas SNPG‐oY1 was present in only 56 samples. Combined genotypes established seven haplotypes (H4, H5, H6, H7, H8, H12 and H19). H6 dominated in northern Eurasia, and H8 showed the second‐highest frequency. H4, which had been earlier reported to be absent in European breeds, was detected in one European breed (Swiniarka), whereas H7, which had been previously identified to be unique to European breeds, was present in two Chinese breeds (Ninglang Black and Large‐tailed Han), one Buryatian (Transbaikal Finewool) and two Russian breeds (North Caucasus Mutton‐Wool and Kuibyshev). H12, which had been detected only in Turkish breeds, was also found in Chinese breeds in this work. An overall low level of haplotype diversity (median h = 0.1288) was observed across the breeds with relatively higher median values in breeds from the regions neighboring the Near Eastern domestication center of sheep. H6 is the dominant haplotype in northwestern and eastern China, in which the haplotype distribution could be explained by the historical translocations of the H4 and H8 Y chromosomes to China via the Mongol invasions followed by expansions to northwestern and eastern China. Our findings extend previous results of sheep Y chromosomal genetic variability and indicate probably recent paternal gene flows between sheep breeds from distinct major geographic regions.     

Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds

A dataset consisting of 787 animals with high‐density SNP chip genotypes (346 774 SNPs) and 939 animals with medium‐density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between F_PED and F_ROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean F_ROH estimates from medium‐density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make‐up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds.     

Resequencing and signatures of selection scan in two Siberian native sheep breeds point to candidate genetic variants for adaptation and economically important traits

Russian sheep breeds represent an important economic asset by providing meat and wool, whilst being adapted to extreme climates. By resequencing two Russian breeds from Siberia: Tuva (n = 20) and Baikal (n = 20); and comparing them with a European (UK) sheep outgroup (n = 14), 41 million variants were called, and signatures of selection were identified. High‐frequency missense mutations on top of selection peaks were found in genes related to immunity (LOC101109746) in the Baikal breed and wool traits (IDUA), cell differentiation (GLIS1) and fat deposition (AADACL3) in the Tuva breed. In addition, genes found under selection owing to haplotype frequency changes were related to wool traits (DSC2), parasite resistance (CLCA1), insulin receptor pathway (SOCS6) and DNA repair (DDB2) in the Baikal breed, and vision (GPR179) in the Tuva breed. Our results present candidate genes and SNPs for future selection programmes, which are necessary to maintain and increase socioeconomic gain from Siberian breeds.     

Beef production potential of Norwegian Red and Holstein-Friesian bulls slaughtered at two ages

There is a paucity of data on the beef production potential of Norwegian Red (NOR) compared with 'modern' Holstein-Friesian (HF) cattle. The present study used a total of 64 bulls in a 2 × 2 factorial design study encompassing two breeds (HF and NOR) and two slaughter ages (485; E, and 610; L, days). The mean initial age and live weight of the HF bulls were 179 (s.d. 47.1) days and 203 (s.d. 64.0) kg, while the corresponding data for the NOR bulls were 176 (s.d. 39.7) days and 185 (s.d. 63.6) kg, respectively. Bulls were offered a 50 : 50 mixture (dry matter (DM) basis) of grass silage and concentrates. No breed × slaughter group interactions were recorded for any parameters evaluated (P > 0.05). HF bulls had higher (P < 0.001) DM intake and poorer (P < 0.01) efficiency of conversion of food to carcass gain than NOR bulls. HF bulls tended (P = 0.07) to have a higher rate of live-weight gain and were heavier (P < 0.001) at slaughter than NOR bulls, though both carcass weight and rate of carcass gain did not differ between the breeds (P > 0.05). NOR bulls had higher (P < 0.001) dressing proportion and carcass conformation score than HF bulls, while breed of bull had no influence (P > 0.05) on carcass fat classification, depth of subcutaneous fat, marbling score or on the weight of fat in the internal depots. Daily food intakes did not differ (P > 0.05) across the two slaughter age groups, though efficiency of conversion of food to carcass gain was poorer (P < 0.05) in the L compared with E bulls. Rate of live-weight gain was lower (P < 0.01) for L bulls, although rate of carcass gain did not differ (P > 0.05) between the E and L bulls. Increasing age at slaughter increased (P < 0.01 or greater) dressing proportion, carcass fat class, depth of subcutaneous fat, marbling score and internal fat depots, but had no effect (P > 0.05) on the carcass conformation score. Instrumental measures of meat quality indicated that meat from NOR bulls was tougher (P < 0.01) than meat from HF bulls, while delaying slaughter increased (P < 0.001) a* and C*ab, and decreased (P < 0.01) h0, indicating improved redness. It is concluded that NOR bulls have higher food efficiency and produce more highly conformed carcasses than HF bulls, but HF bulls produce more tender meat.     

Searching new signals for production traits through gene‐based association analysis in three Italian cattle breeds

Genome‐wide association studies (GWAS) have been widely applied to disentangle the genetic basis of complex traits. In cattle breeds, classical GWAS approaches with medium‐density marker panels are far from conclusive, especially for complex traits. This is due to the intrinsic limitations of GWAS and the assumptions that are made to step from the association signals to the functional variations. Here, we applied a gene‐based strategy to prioritize genotype–phenotype associations found for milk production and quality traits with classical approaches in three Italian dairy cattle breeds with different sample sizes (Italian Brown n = 745; Italian Holstein n = 2058; Italian Simmental n = 477). Although classical regression on single markers revealed only a single genome‐wide significant genotype–phenotype association, for Italian Holstein, the gene‐based approach identified specific genes in each breed that are associated with milk physiology and mammary gland development. As no standard method has yet been established to step from variation to functional units (i.e., genes), the strategy proposed here may contribute to revealing new genes that play significant roles in complex traits, such as those investigated here, amplifying low association signals using a gene‐centric approach.     

Genetic diversity analysis of French goat populations reveals selective sweeps involved in their differentiation

After domestication 11 000 years ago in Asia Minor, the goat followed human migration to Europe and Asia. It was then introduced in Africa and is now raised all over the world. In this study, we exploited a dataset composed of 54 000 SNPs (Illumina goat DNA chip) to analyze the genetic diversity of 223 individuals belonging to eight French breeds (Alpine, Angora, Corse, Fossés, Poitevine, Provençale, Pyrénées and Saanen). Analyses carried out included individual‐based approaches (principal component analysis and population structure) and population‐based approaches (phylogenetic tree constructions). The results of the genetic diversity analyses revealed that French breeds are clearly differentiated, in particular, the Angora breed that originates from south west Asia. The Provençale breed shows a very original genetic pattern that could be the result of ancient admixture. Then, selection signatures were detected by identifying regions of outlying genetic differentiation between populations. Five genomic regions were detected under selection on chromosomes 5, 6, 11, 13 and 20, revealing mainly soft selective sweeps and a few hard selective sweeps and highlighting candidate genes that had been selected for during the evolutionary history of these breeds. Among them, two coat coloration genes (ADAMTS20 and ASIP) and one gene related to milk composition (CSN1S1) were involved.     

Haplotype analysis of key genes governing grain yield and quality traits across 3K RG panel reveals scope for the development of tailor‐made rice with enhanced genetic gains

Though several genes governing various major traits have been reported in rice, their superior haplotype combinations for developing ideal variety remains elusive. In this study, haplotype analysis of 120 previously functionally characterized genes, influencing grain yield (87 genes) and grain quality (33 genes) revealed significant variations in the 3K rice genome (RG) panel. For selected genes, meta‐expression analysis using already available datasets along with co‐expression network provided insights at systems level. Also, we conducted candidate gene based association study for the 120 genes and identified 21 strongly associated genes governing 10‐grain yield and quality traits. We report superior haplotypes upon phenotyping the subset of 3K RG panel, SD1‐H8 with haplotype frequency (HF) of 30.13% in 3K RG panel, MOC1‐H9 (HF: 23.08%), IPA1‐H14 (HF: 6.64%), DEP3‐H2 (HF: 5.59%), DEP1‐H2 (HF: 37.53%), SP1‐H3 (HF: 5.05%), LAX1‐H5 (HF: 1.56%), LP‐H13 (3.64%), OSH1‐H4 (5.52%), PHD1‐H14 (HF: 15.21%), AGO7‐H15 (HF: 3.33%), ROC5‐H2 (31.42%), RSR1‐H8 (HF: 4.20%) and OsNAS3‐H2 (HF: 1.00%). For heading date, Ghd7‐H8 (HF: 3.08%), TOB1‐H10 (HF: 4.60%) flowered early, Ghd7‐H14 (HF: 42.60%), TRX1‐H9 (HF: 27.97%), OsVIL3‐H14 (HF: 1.72%) for medium duration flowering, while Ghd7‐H6 (HF: 1.65%), SNB‐H9 (HF: 9.35%) were late flowering. GS5‐H4 (HF: 65.84%) attributed slender, GS5‐H5 (HF: 29.00%), GW2‐H2 (HF: 4.13%) were medium slender and GS5‐H9 (HF: 2.15%) for bold grains. Furthermore, haplotype analysis explained possible genetic basis for superiority of selected mega‐varieties. Overall, this study suggests the possibility for developing next‐generation tailor‐made rice with superior haplotype combinations of target genes suiting future food and nutritional demands via haplotype‐based breeding.     

Transcriptomic investigation of meat tenderness in two Italian cattle breeds

Our objectives for this study were to understand the biological basis of meat tenderness and to provide an overview of the gene expression profiles related to meat quality as a tool for selection. Through deep mRNA sequencing, we analyzed gene expression in muscle tissues of two Italian cattle breeds: Maremmana and Chianina. We uncovered several differentially expressed genes that encode for proteins belonging to a family of tripartite motif proteins, which are involved in growth, cell differentiation and apoptosis, such as TRIM45, or play an essential role in regulating skeletal muscle differentiation and the regeneration of adult skeletal muscle, such as TRIM32. Other differentially expressed genes (SCN2B, SLC9A7 and KCNK3) emphasize the involvement of potassium–sodium pumps in tender meat. By mapping splice junctions in RNA‐Seq reads, we found significant differences in gene isoform expression levels. The PRKAG3 gene, which is involved in the regulation of energy metabolism, showed four isoforms that were differentially expressed. This distinct pattern of PRKAG3 gene expression could indicate impaired glycogen storage in skeletal muscle, and consequently, this gene very likely has a role in the tenderization process. Furthermore, with this deep RNA‐sequencing, we captured a high number of expressed SNPs, for example, we found 1462 homozygous SNPs showing the alternative allele with a 100% frequency when comparing tender and tough meat. SNPs were then classified into categories by their position and also by their effect on gene coding (174 non‐synonymous polymorphisms) based on the available UMD_3.1 annotations.     

Taking advantage from phenotype variability in a local animal genetic resource: identification of genomic regions associated with the hairless phenotype in Casertana pigs

Casertana is an endangered autochthonous pig breed (raised in south‐central Italy) that is considered to be the descendant of the influential Neapolitan pig population that was used to improve British breeds in the 19th century. Casertana pigs are characterized by a typical, almost complete, hairless phenotype, even though a few Casertana pigs are normal haired. In this work, using Illumina PorcineSNP60 BeadChip data, we carried out a genome‐wide association study and an F_ST analysis with this breed by comparing animals showing the classical hairless phenotype (n = 81) versus pigs classified as haired (n = 15). Combining the results obtained with the two approaches, we identified two significant regions: one on porcine chromosome (SSC) 7 and one on SSC15. The SSC7 region contains the forkhead box N3 (FOXN3) gene, the most plausible candidate gene of this region, considering that mutations in another gene of the same family (forkhead box N1; Foxn1 or FOXN1) are responsible for the nude locus in rodents and alopecia in humans. Another potential candidate gene, rho guanine nucleotide exchange factor 10 (ARHGEF10), is located in the SSC15 region. FOXN3 and ARHGEF10 have been detected as differentially expressed in androgenetic and senescent alopecia respectively. This study on an autochthonous pig breed contributes to shed some light on novel genes potentially involved in hair development and growth and demonstrates that local animal breeds can be valuable genetic resources for disclosing genetic factors affecting unique traits, taking advantage of phenotype variability segregating in small populations.     

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 (F_ROH) with pedigree‐based inbreeding coefficients (F_PED). The linear relationship between F_ROH and F_PED was confirmed for goats by including samples from the three main breeds (Saanen, Chamois and Toggenburg goats). F_ROH 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.     

Signatures of selection in five Italian cattle breeds detected by a 54K SNP panel

In this study we used a medium density panel of SNP markers to perform population genetic analysis in five Italian cattle breeds. The BovineSNP50 BeadChip was used to genotype a total of 2,935 bulls of Piedmontese, Marchigiana, Italian Holstein, Italian Brown and Italian Pezzata Rossa breeds. To determine a genome-wide pattern of positive selection we mapped the F _st values against genome location. The highest F _st peaks were obtained on BTA6 and BTA13 where some candidate genes are located. We identified selection signatures peculiar of each breed which suggest selection for genes involved in milk or meat traits. The genetic structure was investigated by using a multidimensional scaling of the genetic distance matrix and a Bayesian approach implemented in the STRUCTURE software. The genotyping data showed a clear partitioning of the cattle genetic diversity into distinct breeds if a number of clusters equal to the number of populations were given. Assuming a lower number of clusters beef breeds group together. Both methods showed all five breeds separated in well defined clusters and the Bayesian approach assigned individuals to the breed of origin. The work is of interest not only because it enriches the knowledge on the process of evolution but also because the results generated could have implications for selective breeding programs.     

Genome‐wide analysis reveals artificial selection on coat colour and reproductive traits in Chinese domestic pigs

Pigs from Asia and Europe were independently domesticated from c. 9000 years ago. During this period, strong artificial selection has led to dramatic phenotypic changes in domestic pigs. However, the genetic basis underlying these morphological and behavioural adaptations is relatively unknown, particularly for indigenous Chinese pigs. Here, we performed a genome‐wide analysis to screen 196 regions with selective sweep signals in Tongcheng pigs, which are a typical indigenous Chinese breed. Genes located in these regions have been found to be involved in lipid metabolism, melanocyte differentiation, neural development and other biological processes, which coincide with the evolutionary phenotypic changes in this breed. A synonymous substitution, c.669T>C, in ESR1, which colocalizes with a major quantitative trait locus for litter size, shows extreme differences in allele frequency between Tongcheng pigs and wild boars. Notably, the variant C allele in this locus exhibits high allele frequency in most Chinese populations, suggesting a consequence of positive selection. Five genes (PRM1, PRM2, TNP2, GPR149 and JMJD1C) related to reproductive traits were found to have high haplotype similarity in Chinese breeds. Two selected genes, MITF and EDNRB, are implied to shape the two‐end black colour trait in Tongcheng pig. Subsequent SNP microarray studies of five Chinese white‐spotted breeds displayed a concordant signature at both loci, suggesting that these two genes are responsible for colour variations in Chinese breeds. Utilizing massively parallel sequencing, we characterized the candidate sites that adapt to artificial and environmental selections during the Chinese pig domestication. This study provides fundamental proof for further research on the evolutionary adaptation of Chinese pigs.