Genome-wide association studies detects candidate genes for wool traits by re-sequencing in Chinese fine-wool sheep

Correlation of 20 Single-Nucleotide Polymorphisms with Weight and Wool Traits in Alpine Merino Sheep

SNPs associated with important traits of fine-wool sheep that were previously obtained through genome-wide association analysis screening were verified and analyzed. A total of 20 SNPs related to birth weight, bundle strength, cleaning rate, and fiber diameter were screened using whole-genome resequencing, and the SNPshot assay was used to detect and analyze polymorphisms. This study found that, among the 20 SNPs associated with important traits in Alpine Merino sheep, 8 were monomorphic and 12 were polymorphic, of which 6 showed moderate polymorphisms and 6 showed low polymorphisms. The heterozygosity of the 12 polymorphic loci ranged from 0.10 to 0.49, the effective number of alleles ranged from 1.11 to 1.98, and the polymorphic information content ranged from 0.09 to 0.37. The chi-square test showed that only RHPN2:g.42678119T>G and RALYL:g.90030866A>G were in Hardy-Weinberg disequilibrium (p < 0.05); the other loci were in equilibrium (p > 0.05). These SNPs associated with important traits in Alpine Merino sheep provide a theoretical basis for genomic selection and molecular design breeding.

Association between the cashmere production performance, milk production performance, and body size traits and polymorphism of COL6A5 and LOC102181374 genes in Liaoning cashmere goats

The purpose of this study was to analyze the relationship between COL6A5 (collagen type VI alpha 5 chain) and LOC102181374 (alcohol dehydrogenase 1) genes and the production performance of Liaoning cashmere goats by single nucleotide polymorphism (SNP). We have searched for SNP loci of COL6A5 and LOC102181374 genes through sequence alignment and PCR experiments, and have used SPSS and SHEsis software to analyze production data. We obtained five SNP loci in total, including three SNP loci (G50985A, G51140T, G51175A) in COL6A5 gene and two SNP loci (A10067G, T10108C) in LOC102181374 gene. The genotypes G50985A (AG), G51140T (GT), G51175A (AA), A10067G (AA), and T10108C (CC) of these loci have certain advantages in improving the production performance of Liaoning cashmere goats. The haplotype combinations that can improve production performance in COL6A5 gene were H1H5:AGGGAG, H4H4:GGGGAA, and H4H4:GGGGAA. H3H3:GGCC and H2H4:AGTT were the dominant combinations in LOC102181374 gene. At G51175A and A10067G loci, we found that H1H2:AAAG and H1H3:AGAA have dominant effects. These results may provide some support for the molecular breeding of production traits in Liaoning cashmere goats.

Association of SLIT3 and ZNF280B Gene Polymorphisms with Wool Fiber Diameter

The SLIT3 gene encodes a secreted protein, and the ZNF280B gene is a member of the transcription factor family. Both genes have multiple biological functions. This study was conducted to investigate the association between SLIT3 and ZNF280B gene polymorphisms and wool fiber diameter and to determine potential molecular marker sites for breeding sheep with fine wool. We used Kompetitive Allele-Specific PCR to type the single nucleotide polymorphism (SNP) loci in the SLIT3 and ZNF280B genes within 1081 Alpine Merino sheep and associated these SNPs with wool fiber diameter. The results revealed one SNP in SLIT3 and ZNF280B, which were each related to sheep fiber diameter. The wool fiber diameters of sheep with the CC genotype in SLIT3 g.478807C>G and AA genotype in ZNF280B g.677G>A were the smallest and differed significantly from the diameters of other genotypes (p < 0.05). These results suggest potential molecular marker sites for fine-wool sheep breeding.

Genome-Wide Association Study of Fiber Diameter in Alpacas

The aim of this study was the identification of candidate genomic regions associated with fiber diameter in alpacas. DNA samples were collected from 1011 female Huacaya alpacas from two geographical Andean regions in Peru (Pasco and Puno), and three alpaca farms within each region. The samples were genotyped using an Affymetrix Custom Alpaca genotyping array containing 76,508 SNPs. After the quality controls, 960 samples and 51,742 SNPs were retained. Three association study methodologies were performed. The GWAS based on a linear model allowed us to identify 11 and 35 SNPs (-log10(p-values) > 4) using information on all alpacas and alpacas with extreme values of fiber diameter, respectively. The haplotype and marker analysis method allowed us to identify nine haplotypes with standardized haplotype heritability higher than six standard deviations. The selection signatures based on cross-population extended haplotype homozygosity (XP-EHH) allowed us to identify 180 SNPs with XP-EHH values greater than |3|. Four candidate regions with adjacent SNPs identified via two association methods of analysis are located on VPA6, VPA9, VPA29 and one chromosomally unassigned scaffold. This study represents the first analysis of alpaca whole genome association with fiber diameter, using a recently assembled alpaca SNP microarray.

Whole-Genome Resequencing Reveals Selection Signal Related to Sheep Wool Fineness

Wool fineness affects the quality of wool, and some studies have identified about forty candidate genes that affect sheep wool fineness, but these genes often reveal only a certain proportion of the variation in wool thickness. We further explore additional genes associated with the fineness of sheep wool. Whole-genome resequencing of eight sheep breeds was performed to reveal selection signals associated with wool fineness, including four coarse wool and four fine/semi-fine wool sheep breeds. Multiple methods to reveal selection signals (Fst and θπ Ratio and XP-EHH) were applied for sheep wool fineness traits. In total, 269 and 319 genes were annotated in the fine wool (F vs. C) group and the coarse wool (C vs. F) group, such as LGR4, PIK3CA, and SEMA3C and NFIB, OPHN1, and THADA. In F vs. C, 269 genes were enriched in 15 significant GO Terms (p < 0.05) and 38 significant KEGG Pathways (p < 0.05), such as protein localization to plasma membrane (GO: 0072659) and Inositol phosphate metabolism (oas 00562). In C vs. F, 319 genes were enriched in 21 GO Terms (p < 0.05) and 16 KEGG Pathways (p < 0.05), such as negative regulation of focal adhesion assembly (GO: 0051895) and Axon guidance (oas 04360). Our study has uncovered genomic information pertaining to significant traits in sheep and has identified valuable candidate genes. This will pave the way for subsequent investigations into related traits.

Hungarian indigenous Tsigai, a promising breed for excellent heat tolerance and immunity

The adverse effects of climate change on sheep farming have become more noticeable in recent decades. Extensive efforts have been made to untangle the complex relationship between heat tolerance, animal health, and productivity, also to identify a resilient and economically suitable breed for selection that can be resilient to future climate change conditions. Using quantitative real-time polymerase chain reaction (qRT-PCR), we observed the seasonal variations in the expression of several important genes related to heat stress and immunity (HSP70, IL10, TLR2, TLR4, and TLR8) in three of the most widely kept sheep breeds in Hungary: The indigenous Tsigai, Hungarian Merino, and White Dorper. We found that the seasonal stressor affected the relative gene expression of all genes in this study. Notably, The Hungarian indigenous Tsigai was the most robust breed adapted to the Hungarian continental (hot summer, cold winter) environment, with excellent thermotolerance and immunity. Furthermore, despite suffering from heat stress in the summer, Hungarian Merino maintained their robust immune system well throughout the year.

Investigation of SNP markers for the melatonin production trait in the Hu sheep with bulked segregant analysis

BACKGROUND

As an important reproductive hormone, melatonin plays an important role in regulating the reproductive activities of sheep and other mammals. Hu sheep is a breed favoring for meat, with prolific traits. In order to explore the relationship between melatonin and reproductive function of Hu sheep, 7,694,759 SNPs were screened out through the whole genome sequencing analysis from high and low melatonin production Hu sheep.

RESULTS

A total of 68,673 SNPs, involving in 1126 genes, were identified by ED association analysis. Correlation analysis of SNPs of AANAT/ASMT gene and MTNR1A/MTNR1B gene were carried out. The melatonin level of CG genotype 7,981,372 of AANAT, GA genotype 7,981,866 of ASMT and GG genotype 17,355,171 of MTNR1A were higher than the average melatonin level of 1.64 ng/mL. High melatonin Hu sheep appear to have better multiple reproductive performance.

CONCLUSIONS

By using different methods, three SNPs which are associated with high melatonin production trait have been identified in Hu sheep. These 3 SNPs are located in melatonin synthetase AANAT/ASMT and receptor MTNR1A, respectively. Considering the positive association between melatonin production and reproductive performance in ruminants, these three SNPs can be served as the potential molecular markers for breading Hu sheep with the desirable reproductive traits.

Transcriptome analysis reveals genes associated with wool fineness in merinos

Hair/wool usually plays an important role in maintaining animal physiological activities, and the economic value of wool cannot be ignored. At present, people set higher demands on wool fineness. Hence, improving wool fineness is the concern of fine wool sheep breeding. Using RNA-Seq to screen the potential candidate genes that associate with wool fineness can provide theoretical references for fine-wool sheep breeding, and also provide us new ideas for further understand the molecular regulation mechanism of hair growth. In this study, we compared the expression pattern difference of genome-wide genes between the skin transcriptomes of Subo and Chinese Merinos. The results showed that, 16 candidate differentially expressed genes (DEGs) (Included: CACNA1S, GP5, LOC101102392, HSF5, SLITRK2, LOC101104661, CREB3L4, COL1A1, PTPRR, SFRP4, LOC443220, COL6A6, COL6A5, LAMA1, LOC114115342 and LOC101116863 genes) that may associate with wool fineness were screened, and these genes were located in signaling pathways that regulate hair follicle development, cycle or hair growth. It is worth noting that, among the 16 DEGs, COL1A1 gene has the highest expression level in Merino skins, and the fold change of LOC101116863 gene is the highest, and the structures of these two genes are both highly conserved in different species. In conclusion, we speculate that these two genes may play a key role in regulating wool fineness and respectively have similar and conserved functions in different species.

Genome-Wide Scan of Wool Production Traits in Akkaraman Sheep

The objective of this study was to uncover the genetic background of wool quality, a production trait, by estimating genomic heritability and implementing GWAS in Akkaraman sheep. The wool characteristics measured included fibre diameter (FD) and staple length (SL) at the age of 8 months and yearling fibre diameter (YFD), yearling staple length (YSL) and yearling greasy fleece weight (YGFW) at 18 months of age. Animals were genotyped using the Axiom 50 K Ovine Genotyping Array. Maximum likelihood estimations of a linear mixed model (LMM) were used to estimate genomic heritability, where GWAS was conducted following a score test of each trait. Genomic heritability estimates for the traits ranged between 0.22 and 0.63, indicating that phenotypes have a moderate range of heritability. One genome- and six chromosome-wide significant SNPs were associated with the wool traits in Akkaraman lambs. Accordingly, TRIM2, MND1, TLR2, RNF175, CEP290, TMTC3, RERE, SLC45A1, SOX2, MORN1, SKI, FAAP20, PRKCZ, GABRD, CFAP74, CALML6 and TMEM52 genes as well as nine uncharacterized regions (LOC101118971, LOC105609137, LOC105603067, LOC101122892, LOC106991694, LOC106991467, LOC106991455, LOC105616534 and LOC105609719) were defined as plausible candidates. The findings of this study shed light on the genetics of wool quality and yield for the Akkaraman breed and suggests targets for breeders during systematic breeding programmes.

Hair Follicle Transcriptome Analysis Reveals Differentially Expressed Genes That Regulate Wool Fiber Diameter in Angora Rabbits

Wool fiber diameter (WFD) is an important index of wool traits and the main determinant of wool quality and value. However, the genetic determinants of fiber diameter have not yet been fully elucidated. Here, coarse and fine wool of Wan strain Angora rabbits and their hair follicle traits were characterized. The results indicated significant differences in the diameters of wool fibers and their hair follicles. The RNA sequencing (RNA-Seq) technique was used to identify differences in gene expression in hair follicles between coarse and fine wool. In total, 2574 differentially expressed genes (DEGs) were found between the two hair follicle groups. Transcription factors, keratin-associated protein (KAP) and keratin (KRT) families, and ECM-related genes may control the structure of fine fibers in rabbits. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that skin development, epidermal cell and keratinocyte differentiation, epithelium development, and Notch and ribosome signaling pathways were significantly enriched, respectively. GSEA further filtered six important pathways and related core genes. PPI analysis also mined functional DEGs associated with hair structure, including LEF1, FZD3, SMAD3, ITGB6, and BMP4. Our findings provide valuable information for researching the molecular mechanisms regulating wool fiber and could facilitate enhanced selection of super-fine wool rabbits through gene-assisted selection in the future.

Genetic association of wool quality characteristics in United States Rambouillet sheep

Introduction: Fine wool production is an important source of revenue, accounting for up to 13% of total revenue in extensively managed wool sheep production systems of the United States. The Rambouillet are a predominant breed that excels in wool quality characteristics. Understanding the genetic basis of wool quality characteristics would aid in the development of genomic breeding strategies to facilitate genetic improvement. Methods: Wool characteristics and DNA were collected for rams enrolled in the North Dakota State University and University of Wyoming annual central performance ram tests over a three-year period (2019-2021, N = 313). The relationships of wool quality characteristics including grease fleece weight adjusted 365 days (wt. 365 adj.), clean fleece wt. 365 adj., staple length 365 adj., average fiber diameter, face wool cover, amount of skin wrinkles and belly wool were evaluated through genome-wide association studies (GWAS), Pearson correlation and ANOVA. Results: The GWAS identified four genome-wide significant genetic markers (p-value <1.19e-06) and five chromosome-wide significant markers (p-value <1.13e-05) on chromosomes 1, 2, 4, 15, and 19. Significant markers were associated with genes notable for relevant wool biological functions, including the gene ABCC8 which codes for SUR1, an ATP-sensitive potassium channel known to affect hair growth and 60S ribosomal protein L17-like, previously found to be expressed during follicle formation. The strongest Pearson correlation coefficients were identified between clean fleece wt. 365 adj. and grease fleece wt. 365 adj. (r = 0.83) and between clean fleece wt. 365 adj. and staple length 365 adj. (r = 0.53). Additionally, clean fleece wt. 365 adj. was correlated with final body weight (r = 0.35) and scrotal circumference (r = 0.16). Staple length 365 adj. (p-value = 5e-04), average fiber diameter (p-value = .0053) and clean fleece wt. 365 adj. (p-value = .014) were significantly associated with belly wool score. Discussion: The results of this study provide important insight into the relationships between wool quality characteristics and report specific markers that Rambouillet sheep producers may use to help inform selection and breeding decisions for improved wool quality.

Genomic Regions Associated with Wool, Growth and Reproduction Traits in Uruguayan Merino Sheep

The aim of this study was to identify genomic regions and genes associated with the fiber diameter (FD), clean fleece weight (CFW), live weight (LW), body condition score (BCS), pregnancy rate (PR) and lambing potential (LP) of Uruguayan Merino sheep. Phenotypic records of approximately 2000 mixed-age ewes were obtained from a Merino nucleus flock. Genome-wide association studies were performed utilizing single-step Bayesian analysis. For wool traits, a total of 35 genomic windows surpassed the significance threshold (PVE ≥ 0.25%). The proportion of the total additive genetic variance explained by those windows was 4.85 and 9.06% for FD and CFW, respectively. There were 42 windows significantly associated with LWM, which collectively explained 43.2% of the additive genetic variance. For BCS, 22 relevant windows accounted for more than 40% of the additive genetic variance, whereas for the reproduction traits, 53 genomic windows (24 and 29 for PR and LP, respectively) reached the suggestive threshold of 0.25% of the PVE. Within the top 10 windows for each trait, we identified several genes showing potential associations with the wool (e.g., IGF-1, TGFB2R, PRKCA), live weight (e.g., CAST, LAP3, MED28, HERC6), body condition score (e.g., CDH10, TMC2, SIRPA, CPXM1) or reproduction traits (e.g., ADCY1, LEPR, GHR, LPAR2) of the mixed-age ewes.

Multi-locus genome-wide association studies reveal genomic regions and putative candidate genes associated with leaf spot diseases in African groundnut (Arachis hypogaea L.) germplasm

Early leaf spot (ELS) and late leaf spot (LLS) diseases are the two most destructive groundnut diseases in Ghana resulting in ≤ 70% yield losses which is controlled largely by chemical method. To develop leaf spot resistant varieties, the present study was undertaken to identify single nucleotide polymorphism (SNP) markers and putative candidate genes underlying both ELS and LLS. In this study, six multi-locus models of genome-wide association study were conducted with the best linear unbiased predictor obtained from 294 African groundnut germplasm screened for ELS and LLS as well as image-based indices of leaf spot diseases severity in 2020 and 2021 and 8,772 high-quality SNPs from a 48 K SNP array Axiom platform. Ninety-seven SNPs associated with ELS, LLS and five image-based indices across the chromosomes in the 2 two sub-genomes. From these, twenty-nine unique SNPs were detected by at least two models for one or more traits across 16 chromosomes with explained phenotypic variation ranging from 0.01 - 62.76%, with exception of chromosome (Chr) 08 (Chr08), Chr10, Chr11, and Chr19. Seventeen potential candidate genes were predicted at ± 300 kbp of the stable/prominent SNP positions (12 and 5, down- and upstream, respectively). The results from this study provide a basis for understanding the genetic architecture of ELS and LLS diseases in African groundnut germplasm, and the associated SNPs and predicted candidate genes would be valuable for breeding leaf spot diseases resistant varieties upon further validation.

Association of SNPs within PTPN3 gene with wool production and growth traits in a dual-purpose sheep population

Protein tyrosine phosphatase non-receptor type 3 (PTPN3), a member of the membrane-associated non-receptor protein tyrosine phosphatase (PTP) family, plays significant roles in the cytoplasm and affects the development and growth of skin and hair. A recent study identified the PTPN3 as the potential gene related to sheep wool quality. To detect single-nucleotide polymorphisms (SNPs) of PTPN3 and elucidate its association with wool production and growth traits in fine wool sheep a total of 644 healthy SG (South African mutton merino♂ × Gansu alpine fine-wool sheep♀, SG) and SSG (South African mutton merino♂ × SG♀, SSG) hybrid sheep were selected. Pooled-DNA sequencing and SNPscan methods were used to scan and genotype SNPs within PTPN3. Association analyses between SNPs and wool production and growth traits were implemented. Consequently, the results revealed that PTPN3 has six SNPs (two missense mutations, one synonymous mutation, and three intron mutations), of which four loci (SNP2, SNP3, SNP4, and SNP5) were significantly positively correlated with growth and wool traits (p < 0.05). SNP4 was significantly (p < 0.05) linked with thigh wool length, and SNP6 was significantly (p < 0.05) associated with abdomen wool length. Moreover, one strongly linked SNP block was identified to be correlated with wool production and growth traits (body weight and body size). The significant SNPs founded by this study could serve as useful genetic markers for breeding fine-wool sheep.^®.