Genome-wide association study for wool production traits in a Chinese Merino sheep population

Runs of homozygosity reveal genome-wide autozygosity in Italian sheep breeds

The availability of dense single nucleotide polymorphism (SNP) assays allows for the determination of autozygous segments based on runs of consecutive homozygous genotypes (ROH). The aim of the present study was to investigate the occurrence and distribution of ROH in 21 Italian sheep breeds using medium-density SNP genotypes in order to characterize autozygosity and identify genomic regions that frequently appeared in ROH within individuals, namely ROH islands. After filtering, the final number of animals and SNPs retained for analyses were 502 and 46 277 respectively. A total of 12 302 ROH were identified. The mean number of ROH per breed ranged from 10.58 (Comisana) to 44.54 (Valle del Belice). The average length of ROH across breeds was 4.55 Mb and ranged from 3.85 Mb (Biellese) to 5.51 Mb (Leccese). Valle del Belice showed the highest value of inbreeding on the basis of ROH (F_ROH  = 0.099), whereas Comisana showed the lowest (F_ROH  = 0.016), and high standard deviation values revealed high variability in autozygosity levels within each breed. Differences also existed in the length of ROH. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were <10 Mb in length, with a few long ROH >25 Mb. The levels of ROH that we estimated here reflect the inbreeding history of the investigated sheep breeds. These results also highlight that ancient and recent inbreeding have had an impact on the genome of the Italian sheep breeds and suggest that several animals have experienced recent autozygosity events. Comisana and Bergamasca appeared as the less consanguineous breeds, whereas Barbaresca, Leccese and Valle del Belice showed ROH patterns typically produced by recent inbreeding. Moreover, within the genomic regions most commonly associated with ROH, several candidate genes were detected.

High-resolution analysis of selection sweeps identified between fine-wool Merino and coarse-wool Churra sheep breeds

BACKGROUND

With the aim of identifying selection signals in three Merino sheep lines that are highly specialized for fine wool production (Australian Industry Merino, Australian Merino and Australian Poll Merino) and considering that these lines have been subjected to selection not only for wool traits but also for growth and carcass traits and parasite resistance, we contrasted the OvineSNP50 BeadChip (50 K-chip) pooled genotypes of these Merino lines with the genotypes of a coarse-wool breed, phylogenetically related breed, Spanish Churra dairy sheep. Genome re-sequencing datasets of the two breeds were analyzed to further explore the genetic variation of the regions initially identified as putative selection signals.

RESULTS

Based on the 50 K-chip genotypes, we used the overlapping selection signals (SS) identified by four selection sweep mapping analyses (that detect genetic differentiation, reduced heterozygosity and patterns of haplotype diversity) to define 18 convergence candidate regions (CCR), five associated with positive selection in Australian Merino and the remainder indicating positive selection in Churra. Subsequent analysis of whole-genome sequences from 15 Churra and 13 Merino samples identified 142,400 genetic variants (139,745 bi-allelic SNPs and 2655 indels) within the 18 defined CCR. Annotation of 1291 variants that were significantly associated with breed identity between Churra and Merino samples identified 257 intragenic variants that caused 296 functional annotation variants, 275 of which were located across 31 coding genes. Among these, four synonymous and four missense variants (NPR2_His847Arg, NCAPG_Ser585Phe, LCORL_Asp1214Glu and LCORL_Ile1441Leu) were included.

CONCLUSIONS

Here, we report the mapping and genetic variation of 18 selection signatures that were identified between Australian Merino and Spanish Churra sheep breeds, which were validated by an additional contrast between Spanish Merino and Churra genotypes. Analysis of whole-genome sequencing datasets allowed us to identify divergent variants that may be viewed as candidates involved in the phenotypic differences for wool, growth and meat production/quality traits between the breeds analyzed. The four missense variants located in the NPR2, NCAPG and LCORL genes may be related to selection sweep regions previously identified and various QTL reported in sheep in relation to growth traits and carcass composition.

A comparison of transcriptomic patterns measured in the skin of Chinese fine and coarse wool sheep breeds

We characterised wool traits, and skin gene expression profiles of fine wool Super Merino (SM) and coarse wool Small Tail Han (STH) sheep. SM sheep had a significantly higher total density of wool follicles, heavier fleeces, finer fibre diameter, and increased crimp frequency, staple length and wool grease (lanolin) production. We found 435 genes were expressed at significantly different levels in the skin of the two breeds (127 genes more highly in SM and 308 genes more highly in STH sheep). Classification of the genes more highly expressed in SM sheep revealed numerous lipid metabolic genes as well as genes encoding keratins, keratin-associated proteins, and wool follicle stem cell markers. In contrast, mammalian epidermal development complex genes and other genes associated with skin cornification and muscle function were more highly expressed in STH sheep. Genes identified in this study may be further evaluated for inclusion in breeding programs, or as targets for therapeutic or genetic interventions, aimed at altering wool quality or yield. Expression of the lipid metabolic genes in the skin of sheep may be used as a novel trait with the potential to alter the content or properties of lanolin or the fleece.

Multiple-trait QTL mapping and genomic prediction for wool traits in sheep

BACKGROUND

The application of genomic selection to sheep breeding could lead to substantial increases in profitability of wool production due to the availability of accurate breeding values from single nucleotide polymorphism (SNP) data. Several key traits determine the value of wool and influence a sheep's susceptibility to fleece rot and fly strike. Our aim was to predict genomic estimated breeding values (GEBV) and to compare three methods of combining information across traits to map polymorphisms that affect these traits.

METHODS

GEBV for 5726 Merino and Merino crossbred sheep were calculated using BayesR and genomic best linear unbiased prediction (GBLUP) with real and imputed 510,174 SNPs for 22 traits (at yearling and adult ages) including wool production and quality, and breech conformation traits that are associated with susceptibility to fly strike. Accuracies of these GEBV were assessed using fivefold cross-validation. We also devised and compared three approximate multi-trait analyses to map pleiotropic quantitative trait loci (QTL): a multi-trait genome-wide association study and two multi-trait methods that use the output from BayesR analyses. One BayesR method used local GEBV for each trait, while the other used the posterior probabilities that a SNP had an effect on each trait.

RESULTS

BayesR and GBLUP resulted in similar average GEBV accuracies across traits (~0.22). BayesR accuracies were highest for wool yield and fibre diameter (>0.40) and lowest for skin quality and dag score (<0.10). Generally, accuracy was higher for traits with larger reference populations and higher heritability. In total, the three multi-trait analyses identified 206 putative QTL, of which 20 were common to the three analyses. The two BayesR multi-trait approaches mapped QTL in a more defined manner than the multi-trait GWAS. We identified genes with known effects on hair growth (i.e. FGF5, STAT3, KRT86, and ALX4) near SNPs with pleiotropic effects on wool traits.

CONCLUSIONS

The mean accuracy of genomic prediction across wool traits was around 0.22. The three multi-trait analyses identified 206 putative QTL across the ovine genome. Detailed phenotypic information helped to identify likely candidate genes.

Detection of QTL for greasy fleece weight in sheep using a 50 K single nucleotide polymorphism chip

Genome-wide association studies (GWAS) have introduced an influential tool in the search for quantitative trait loci (QTL) influencing economically important traits in sheep. To identify QTL associated with greasy fleece weight, a GWAS with 50 K single nucleotide polymorphisms (SNPs) was performed in a Baluchi sheep population. Association with greasy fleece weights was tested using the software Plink. The results of our GWAS provided three novel SNP markers and candidate genes associated with greasy fleece weight. A total of three chromosome-wide significant associations were detected for SNP on chromosomes 17 and 20 affecting greasy fleece weight across the four shearing. One of the significant SNP markers was located within ovine known genes namely FAM101A. Further investigation of these identified regions in validation studies will facilitate the identification of strong candidate genes for wool production in sheep.

Estimates of linkage disequilibrium and effective population sizes in Chinese Merino (Xinjiang type) sheep by genome-wide SNPs

Knowledge of linkage disequilibrium (LD) is important for effective genome-wide association studies and accurate genomic prediction. Chinese Merino (Xinjiang type) is well-known fine wool sheep breed. However, the extent of LD across the genome remains unexplored. In this study, we calculated autosomal LD based on genome-wide SNPs of 635 Chinese Merino (Xinjiang type) sheep by Illumina Ovine SNP50 BeadChip. A moderate level of LD (r² ≥ 0.25) across the whole genome was observed at short distances of 0–10 kb. Further, the ancestral effective population size (N_e) was analyzed by extent of LD and found that N_e increased with the increase of generations and declined rapidly within the most recent 50 generations, which is consistent with the history of Chinese Merino sheep breeding, initiated in 1971. We also noted that even when the effective population size was estimated across different single chromosomes, N_e only ranged from 140.36 to 183.33 at five generations in the past, exhibiting a rapid decrease compared with that at ten generations in the past. These results indicated that the genetic diversity in Chinese Merino sheep recently decreased and proper protective measures should be taken to maintain the diversity. Our datasets provided essential genetic information to track molecular variations which potentially contribute to phenotypic variation in Chinese Merino sheep.

Polymorphisms of FST gene and their association with wool quality traits in Chinese Merino sheep

Follistatin (FST) is involved in hair follicle morphogenesis. However, its effects on hair traits are not clear. This study was designed to investigate the effects of FST gene single nucleotide polymorphisms (SNP) on wool quality traits in Chinese Merino sheep (Junken Type). We performed gene expression analysis, SNP detection, and association analysis of FST gene with sheep wool quality traits. The real-time RT-PCR analysis showed that FST gene was differentially expressed in adult skin between Chinese Merino sheep (Junken Type) and Suffolk sheep. Immunostaining showed that FST was localized in inner root sheath (IRS) and matrix of hair follicle (HF) in both SF and Suffolk sheep. Sequencing analysis identified a total of seven SNPs (termed SNPs 1-7) in the FST gene in Chinese Merino sheep (Junken Type). Association analysis showed that SNP2 (Chr 16. 25,633,662 G>A) was significantly associated with average wool fiber diameter, wool fineness SD, and wool crimp (P < 0.05). SNP4 (Chr 16. 25,633,569 C>T) was significantly associated with wool fineness SD and CV of fiber diameter (P < 0.05). Similarly, the haplotypes derived from these seven identified SNPs were also significantly associated with average wool fiber diameter, wool fineness SD, CV of fiber diameter, and wool crimp (P < 0.05). Our results suggest that FST influences wool quality traits and its SNPs 2 and 4 might be useful markers for marker-assisted selection and sheep breeding.

Unveiling genomic regions that underlie differences between Afec-Assaf sheep and its parental Awassi breed

BACKGROUND

Sheep production in Israel has improved by crossing the fat-tailed local Awassi breed with the East Friesian and later, with the Booroola Merino breed, which led to the formation of the highly prolific Afec-Assaf strain. This strain differs from its parental Awassi breed in morphological traits such as tail and horn size, coat pigmentation and wool characteristics, as well as in production, reproductive and health traits. To identify major genes associated with the formation of the Afec-Assaf strain, we genotyped 41 Awassi and 141 Afec-Assaf sheep using the Illumina Ovine SNP50 BeadChip array, and analyzed the results with PLINK and EMMAX software. The detected variable genomic regions that differed between Awassi and Afec-Assaf sheep (variable genomic regions; VGR) were compared to selection signatures that were reported in 48 published genome-wide association studies in sheep. Because the Afec-Assaf strain, but not the Awassi breed, carries the Booroola mutation, association analysis of BMPR1B used as the test gene was performed to evaluate the ability of this study to identify a VGR that includes such a major gene.

RESULTS

Of the 20 detected VGR, 12 were novel to this study. A ~7-Mb VGR was identified on Ovies aries chromosome OAR6 where the Booroola mutation is located. Similar to other studies, the most significant VGR was detected on OAR10, in a region that contains candidate genes affecting horn type (RXFP2), climate adaptation (ALOX5AP), fiber diameter (KATNAl1), coat pigmentation (FRY) and genes associated with fat distribution. The VGR on OAR2 included BNC2, which is also involved in controlling coat pigmentation in sheep. Six other VGR contained genes that were shown to be involved in coat pigmentation by analyzing their mammalian orthologues. Genes associated with fat distribution in humans, including GRB14 and COBLL1, were located in additional VGR. Sequencing DNA from Awassi and Afec-Assaf individuals revealed non-synonymous mutations in some of these candidate genes.

CONCLUSIONS

Our results highlight VGR that differentiate the Awassi breed from the Afec-Assaf strain, some of which may include genes that confer an advantage to Afec-Assaf and Assaf over Awassi sheep with respect to intensive sheep production under Mediterranean conditions.

Mutations in TSPEAR, Encoding a Regulator of Notch Signaling, Affect Tooth and Hair Follicle Morphogenesis

Despite recent advances in our understanding of the pathogenesis of ectodermal dysplasias (EDs), the molecular basis of many of these disorders remains unknown. In the present study, we aimed at elucidating the genetic basis of a new form of ED featuring facial dysmorphism, scalp hypotrichosis and hypodontia. Using whole exome sequencing, we identified 2 frameshift and 2 missense mutations in TSPEAR segregating with the disease phenotype in 3 families. TSPEAR encodes the thrombospondin-type laminin G domain and EAR repeats (TSPEAR) protein, whose function is poorly understood. TSPEAR knock-down resulted in altered expression of genes known to be regulated by NOTCH and to be involved in murine hair and tooth development. Pathway analysis confirmed that down-regulation of TSPEAR in keratinocytes is likely to affect Notch signaling. Accordingly, using a luciferase-based reporter assay, we showed that TSPEAR knock-down is associated with decreased Notch signaling. In addition, NOTCH1 protein expression was reduced in patient scalp skin. Moreover, TSPEAR silencing in mouse hair follicle organ cultures was found to induce apoptosis in follicular epithelial cells, resulting in decreased hair bulb diameter. Collectively, these observations indicate that TSPEAR plays a critical, previously unrecognized role in human tooth and hair follicle morphogenesis through regulation of the Notch signaling pathway.

Ectodermal dysplasias refer to a large group of inherited disorders characterized by developmental defects in tissues of ectodermal origin. The study of these conditions has been instrumental in the discovery of biological pathways involved in the regulation of epithelial tissue morphogenesis. In this report, through the delineation of the molecular basis of a novel form of autosomal recessive ectodermal dysplasia, we identified a new key player in ectodermal development. We detected a number of mutations in TSPEAR co-segregating with abnormal hair and tooth development in three families. TSPEAR encodes the thrombospondin-type laminin G domain and EAR repeats (TSPEAR) protein, whose function is poorly understood. TSPEAR was found to be strongly expressed in murine hair and tooth. Using a reporter assay, we showed that it regulates Notch activity. Accordingly, NOTCH1 expression was altered in patient skin, and NOTCH1, as well as many of its known targets, was down-regulated in TSPEAR deficient keratinocytes. Moreover, Tspear silencing in mouse hair follicle organ cultures was found to induce apoptosis in follicular epithelial cells, resulting in decreased hair bulb diameter. Collectively, these observations indicate that TSPEAR plays a critical, previously unrecognized role in human tooth and hair follicle morphogenesis through regulation of the Notch pathway. As such, these new data are likely to lead to further investigations aimed at characterizing the role of Notch signaling pathway in other forms of ectodermal dysplasias as well as acquired hair and tooth pathologies.

Association of methionine synthase gene polymorphisms with wool production and quality traits in Chinese Merino population

Methionine synthase (MTR) plays a crucial role in maintaining homeostasis of intracellular methionine, folate, and homocysteine, and its activity correlates with DNA methylation in many mammalian tissues. Our previous genomewide association study identified that 1 SNP located in the gene was associated with several wool production and quality traits in Chinese Merino. To confirm the potential involvement of the gene in sheep wool production and quality traits, we performed sheep tissue expression profiling, SNP detection, and association analysis with sheep wool production and quality traits. The semiquantitative reverse transcription PCR analysis showed that the gene was differentially expressed in skin from Merino and Kazak sheep. The sequencing analysis identified a total of 13 SNP in the gene from Chinese Merino sheep. Comparison of the allele frequencies revealed that these 13 identified SNP were significantly different among the 6 tested Chinese Merino strains ( < 0.001). Linkage disequilibrium analysis showed that SNP 3 to 11 were strongly linked in a single haplotype block in the tested population. Association analysis showed that SNP 2 to 11 were significantly associated with the average wool fiber diameter and the fineness SD and that SNP 4 to 11 were significantly associated with the CV of fiber diameter trait ( < 0.05). Single nucleotide polymorphism 2 and SNP 5 to 12 were weakly associated with wool crimp. Similarly, the haplotypes derived from these 13 identified SNP were also significantly associated with the average wool fiber diameter, fineness SD, and the CV of fiber diameter ( < 0.05). Our results suggest that is a candidate gene for sheep wool production and quality traits, and the identified SNP might be used in sheep breeding.

Wool Keratin-Associated Protein Genes in Sheep-A Review

The importance of sheep's wool in making textiles has inspired extensive research into its structure and the underlying genetics since the 1960s. Wool keratin-associated proteins (KAPs) are a key structural component of the wool fibre. The characterisation of the genes encoding these proteins has progressed rapidly with advances in the nucleotide and protein sequencing. This review describes our knowledge of ovine KAPs, their categorisation into families, polymorphism in the proteins and genes, the clustering and chromosomal location of the genes, some characteristics of gene expression and some potential effects of the KAPs on wool traits. The extent and nature of genetic variation in wool KAP genes and its association with fibre characteristics, provides an opportunity for the development of gene-markers for selective breeding of sheep to produce better wool with properties highly matched to specific end-uses.

Genome-Wide Specific Selection in Three Domestic Sheep Breeds

BACKGROUND

Commercial sheep raised for mutton grow faster than traditional Chinese sheep breeds. Here, we aimed to evaluate genetic selection among three different types of sheep breed: two well-known commercial mutton breeds and one indigenous Chinese breed.

RESULTS

We first combined locus-specific branch lengths and di statistical methods to detect candidate regions targeted by selection in the three different populations. The results showed that the genetic distances reached at least medium divergence for each pairwise combination. We found these two methods were highly correlated, and identified many growth-related candidate genes undergoing artificial selection. For production traits, APOBR and FTO are associated with body mass index. For meat traits, ALDOA, STK32B and FAM190A are related to marbling. For reproduction traits, CCNB2 and SLC8A3 affect oocyte development. We also found two well-known genes, GHR (which affects meat production and quality) and EDAR (associated with hair thickness) were associated with German mutton merino sheep. Furthermore, four genes (POL, RPL7, MSL1 and SHISA9) were associated with pre-weaning gain in our previous genome-wide association study.

CONCLUSIONS

Our results indicated that combine locus-specific branch lengths and di statistical approaches can reduce the searching ranges for specific selection. And we got many credible candidate genes which not only confirm the results of previous reports, but also provide a suite of novel candidate genes in defined breeds to guide hybridization breeding.