Characterization of the porcine KIT ligand gene: expression analysis, genomic structure, polymorphism detection and association with coat colour traits

Selection signatures of wool color in Gangba sheep revealed by genome-wide SNP discovery

BACKGROUND

Gangba sheep as a famous breed of Tibetan sheep, its wool color is mainly white and black. Gangba wool is economically important as a high-quality raw material for Tibetan blankets and Tibetan serge. However, relatively few studies have been conducted on the wool color of Tibetan sheep.

RESULTS

To fill this research gap, this study conducted an in-depth analysis of two populations of Gangba sheep (black and white wool color) using whole genome resequencing to identify genetic variation associated with wool color. Utilizing PCA, Genetic Admixture, and N-J Tree analyses, the present study revealed a consistent genetic relationship and structure between black and white wool colored Gangba sheep populations, which is consistent with their breed history. Analysis of selection signatures using multiple methods (FST, π ratio, Tajima's D), 370 candidate genes were screened in the black wool group (GBB vs GBW); among them, MC1R, MLPH, SPIRE2, RAB17, SMARCA4, IRF4, CAV1, USP7, TP53, MYO6, MITF, MC2R, TET2, NF1, JAK1, GABRR1 genes are mainly associated with melanin synthesis, melanin delivery, and distribution. The enrichment results of the candidate genes identified 35 GO entries and 19 KEGG pathways associated with the formation of the black phenotype. 311 candidate genes were screened in the white wool group (GBW vs GBB); among them, REST, POU2F1, ADCY10, CCNB1, EP300, BRD4, GLI3, and SDHA genes were mainly associated with interfering with the differentiation of neural crest cells into melanocytes, affecting the proliferation of melanocytes, and inhibiting melanin synthesis. 31 GO entries and 22 KEGG pathways were associated with the formation of the white phenotype.

CONCLUSIONS

This study provides important information for understanding the genetic mechanism of wool color in Gangba, and provides genetic knowledge for improving and optimizing the wool color of Tibetan sheep. Genetic improvement and selective breeding to produce wool of specific colors can meet the demand for a diversity of wool products in the Tibetan wool textile market.

Therapeutic modulation of KIT ligand in melanocytic disorders with implications for mast cell diseases

KIT ligand and its associated receptor KIT serve as a master regulatory system for both melanocytes and mast cells controlling survival, migration, proliferation and activation. Blockade of this pathway results in cell depletion, while overactivation leads to mastocytosis or melanoma. Expression defects are associated with pigmentary and mast cell disorders. KIT ligand regulation is complex but efficient targeting of this system would be of significant benefit to those suffering from melanocytic or mast cell disorders. Herein, we review the known associations of this pathway with cutaneous diseases and the regulators of this system both in skin and in the more well-studied germ cell system. Exogenous agents modulating this pathway will also be presented. Ultimately, we will review potential therapeutic opportunities to help our patients with melanocytic and mast cell disease processes potentially including vitiligo, hair greying, melasma, urticaria, mastocytosis and melanoma.

De novo mutation in KITLG gene causes a variant of Familial Progressive Hyper- and Hypo-pigmentation (FPHH)

Familial Progressive Hyper‐ and Hypopigmentation is a pigmentary disorder characterized by a mix of hypo‐ and hyperpigmented lesions, café‐au‐lait spots and hypopigmented ash‐leaf macules. The disorder was previously linked to KITLG and various mutations have been reported to segregate in different families. Furthermore, association between KITLG mutations and malignancies was also suggested. Exome and SANGER sequencing were performed for identification of KITLG mutations. Functional in silico analyses were additionally performed to assess the findings. We identified a de novo mutation in exon 4 of KITLG gene causing NM_000899.4:c.[329A>T] (chr12:88912508A>T) leading to NP_000890.1:p.(Asp110Val) substitution in the 3rd alpha helix. It was predicted as pathogenic, located in a conserved region and causing an increase in hydrophobicity in the KITLG protein. Our findings clearly confirm an additional hot spot of KITLG mutations in the 3rd alpha helix, which very likely increases the risk of malignancies. To our knowledge the present study provides the strongest evidence of association of the KITLG mutation with both Familial Progressive Hyper‐ and Hypopigmentation and malignancy due to its’ location on somatic cancer mutation locus. Additionally we also address difficulties with classification of the unique phenotype and propose a subtype within broader diagnosis.

Genomic Characteristics and Selection Signatures in Indigenous Chongming White Goat (Capra hircus)

The Chongming white goat (CM) is an indigenous goat breed exhibits unique traits that are adapted to the local environment and artificial selection. By performing whole-genome re-sequencing, we generated 14-20× coverage sequences from 10 domestic goat breeds to explore the genomic characteristics and selection signatures of the CM breed. We identified a total of 23,508,551 single-nucleotide polymorphisms (SNPs) and 2,830,800 insertion-deletion mutations (indels) after read mapping and variant calling. We further specifically identified 1.2% SNPs (271,713) and 0.9% indels (24,843) unique to the CM breed in comparison with the other nine goat breeds. Missense (SIFT < 0.05), frameshift, splice-site, start-loss, stop-loss, and stop-gain variants were identified in 183 protein-coding genes of the CM breed. Of the 183, 36 genes, including AP4E1, FSHR, COL11A2, and DYSF, are involved in phenotype ontology terms related to the nervous system, short stature, and skeletal muscle morphology. Moreover, based on genome-wide F ST and pooled heterozygosity (Hp) calculation, we further identified selection signature genes between the CM and the other nine goat breeds. These genes are significantly associated with the nervous system (C2CD3, DNAJB13, UCP2, ZMYND11, CEP126, SCAPER, and TSHR), growth (UCP2, UCP3, TSHR, FGFR1, ERLIN2, and ZNF703), and coat color (KITLG, ASIP, AHCY, RALY, and MC1R). Our results suggest that the CM breed may be differentiated from other goat breeds in terms of nervous system owing to natural or artificial selection. The whole-genome analysis provides an improved understanding of genetic diversity and trait exploration for this indigenous goat breed.

Genome-Wide Identification of Discriminative Genetic Variations in Beef and Dairy Cattle via an Information-Theoretic Approach

Analyzing the associations between genotypic changes and phenotypic traits on a genome-wide scale can contribute to understanding the functional roles of distinct genetic variations during breed development. We performed a whole-genome analysis of Angus and Jersey cattle breeds using conditional mutual information, which is an information-theoretic method estimating the conditional independency among multiple factor variables. The proposed conditional mutual information-based approach allows breed-discriminative genetic variations to be explicitly identified from tens of millions of SNP (single nucleotide polymorphism) positions on a genome-wide scale while minimizing the usage of prior knowledge. Using this data-driven approach, we identified biologically relevant functional genes, including breed-specific variants for cattle traits such as beef and dairy production. The identified lipid-related genes were shown to be significantly associated with lipid and triglyceride metabolism, fat cell differentiation, and muscle development. In addition, we confirmed that milk-related genes are involved in mammary gland development, lactation, and mastitis-associated processes. Our results provide the distinct properties of Angus and Jersey cattle at a genome-wide level. Moreover, this study offers important insights into discovering unrevealed genetic variants for breed-specific traits and the identification of genetic signatures of diverse cattle breeds with respect to target breed-specific properties.

Discovery of selection-driven genetic differences of Duroc, Landrace, and Yorkshire pig breeds by EigenGWAS and Fst analyses

Pigs are one of the earliest domesticated animals and multiple breeds have been developed to meet the various demands of consumers. EigenGWAS is a novel strategy to identify candidate genes that underlying population genetic differences and to infer candidate regions under selection as well. In this study, EigenGWAS and F_st analyses were performed using the public re-sequencing data of three typical commercial pig breeds, Duroc, Landrace and Yorkshire. The intersection of genome-wide significant SNPs detected by EigenGWAS and top-ranked 1% SNPs of F_st results were treated as signals under selection. Using the data of all three breeds, 3062 signals under selection were detected and the nearby genomic regions within 300 kb upstream and downstream covered 6.54% of whole genome. Pairs of breeds were analysed along with the pathway analysis. The gene function enrichment results indicated that many candidate genes located in the genomic regions of the signals under selection were associated with biological processes related to growth, metabolism, reproduction, sensory perception, etc. Among the candidate genes, the FSHB, AHR, PTHLH, KDR and FST genes were reported to be associated with reproductive performance; the KIT, KITLG, MITF, MC1R and EDNRB genes were previously identified to affect coat colour; the RETREG1, TXNIP, BMP5, PPARD and RBP4 genes were reported to be associated with lipid metabolism and growth traits. The identified genetic differences across the three commercial breeds will advance understanding of the artificial selection history of pigs and the signals under selection will suggest potential uses in pig genomic breeding programmes.

Comparative Transcriptome Analysis of Mink (Neovison vison) Skin Reveals the Key Genes Involved in the Melanogenesis of Black and White Coat Colour

Farmed mink (Neovison vison) is one of the most important fur-bearing species worldwide, and coat colour is a crucial qualitative characteristic that contributes to the economic value of the fur. To identify additional genes that may play important roles in coat colour regulation, Illumina/Solexa high-throughput sequencing technology was used to catalogue the global gene expression profiles in mink skin with two different coat colours (black and white). RNA-seq analysis indicated that a total of 12,557 genes were differentially expressed in black versus white minks, with 3,530 genes up-regulated and 9,027 genes down-regulated in black minks. Significant differences were not observed in the expression of MC1R and TYR between the two different coat colours, and the expression of ASIP was not detected in the mink skin of either coat colour. The expression levels of KITLG, LEF1, DCT, TYRP1, PMEL, Myo5a, Rab27a and SLC7A11 were validated by qRT-PCR, and the results were consistent with RNA-seq analysis. This study provides several candidate genes that may be associated with the development of two coat colours in mink skin. These results will expand our understanding of the complex molecular mechanisms underlying skin physiology and melanogenesis in mink and will provide a foundation for future studies.

Whole-genome sequencing of eight goat populations for the detection of selection signatures underlying production and adaptive traits

The goat (Capra hircus) is one of the first farm animals that have undergone domestication and extensive natural and artificial selection by adapting to various environments, which in turn has resulted in its high level of phenotypic diversity. Here, we generated medium-coverage (9–13×) sequences from eight domesticated goat breeds, representing morphologically or geographically specific populations, to identify genomic regions representing selection signatures. We discovered ~10 million single nucleotide polymorphisms (SNPs) for each breed. By combining two approaches, ZH_p and di values, we identified 22 genomic regions that may have contributed to the phenotypes in coat color patterns, body size, cashmere traits, as well as high altitude adaptation in goat populations. Candidate genes underlying strong selection signatures including coloration (ASIP, KITLG, HTT, GNA11, and OSTM1), body size (TBX15, DGCR8, CDC25A, and RDH16), cashmere traits (LHX2, FGF9, and WNT2), and hypoxia adaptation (CDK2, SOCS2, NOXA1, and ENPEP) were identified. We also identified candidate functional SNPs within selected genes that may be important for each trait. Our results demonstrated the potential of using sequence data in identifying genomic regions that are responsible for agriculturally significant phenotypes in goats, which in turn can be used in the selection of goat breeds for environmental adaptation and domestication.

Signatures of diversifying selection in European pig breeds

Following domestication, livestock breeds have experienced intense selection pressures for the development of desirable traits. This has resulted in a large diversity of breeds that display variation in many phenotypic traits, such as coat colour, muscle composition, early maturity, growth rate, body size, reproduction, and behaviour. To better understand the relationship between genomic composition and phenotypic diversity arising from breed development, the genomes of 13 traditional and commercial European pig breeds were scanned for signatures of diversifying selection using the Porcine60K SNP chip, applying a between-population (differentiation) approach. Signatures of diversifying selection between breeds were found in genomic regions associated with traits related to breed standard criteria, such as coat colour and ear morphology. Amino acid differences in the EDNRB gene appear to be associated with one of these signatures, and variation in the KITLG gene may be associated with another. Other selection signals were found in genomic regions including QTLs and genes associated with production traits such as reproduction, growth, and fat deposition. Some selection signatures were associated with regions showing evidence of introgression from Asian breeds. When the European breeds were compared with wild boar, genomic regions with high levels of differentiation harboured genes related to bone formation, growth, and fat deposition.

The domestic pig, an important source of protein worldwide, was domesticated from the ancestral wild boar in multiple locations throughout the world. In Europe, local types were developed following domestication, but phenotypically distinct breeds only arose in the eighteenth century with the advent of systematic breeding. Recently developed molecular tools for pigs (as well as other livestock species) now allow a genetic characterisation of breed histories, including identification of regions of the genome that have been under selection in the establishment of breeds. We have applied these tools to identify genomic regions associated with breed development in a set of commercial and traditional pig breeds. We found strong evidence of genetic differentiation between breeds near genes associated with traits that are used to define breed standards, such as ear morphology and coat colour, as well as in regions of the genome that are associated with pork production traits. It is well documented that crosses with Asian pigs have been used to modify European breeds. We have found evidence of genetic influence from Asian pigs in European breeds, again in regions of the genome associated with breed standard characteristics, including ear shape and coat colour, as well as production traits.

Skipping of exons by premature termination of transcription and alternative splicing within intron-5 of the sheep SCF gene: a novel splice variant

Stem cell factor (SCF) is a growth factor, essential for haemopoiesis, mast cell development and melanogenesis. In the hematopoietic microenvironment (HM), SCF is produced either as a membrane-bound (-) or soluble (+) forms. Skin expression of SCF stimulates melanocyte migration, proliferation, differentiation, and survival. We report for the first time, a novel mRNA splice variant of SCF from the skin of white merino sheep via cloning and sequencing. Reverse transcriptase (RT)-PCR and molecular prediction revealed two different cDNA products of SCF. Full-length cDNA libraries were enriched by the method of rapid amplification of cDNA ends (RACE-PCR). Nucleotide sequencing and molecular prediction revealed that the primary 1519 base pair (bp) cDNA encodes a precursor protein of 274 amino acids (aa), commonly known as 'soluble' isoform. In contrast, the shorter (835 and/or 725 bp) cDNA was found to be a 'novel' mRNA splice variant. It contains an open reading frame (ORF) corresponding to a truncated protein of 181 aa (vs 245 aa) with an unique C-terminus lacking the primary proteolytic segment (28 aa) right after the D(175)G site which is necessary to produce 'soluble' form of SCF. This alternative splice (AS) variant was explained by the complete nucleotide sequencing of splice junction covering exon 5-intron (5)-exon 6 (948 bp) with a premature termination codon (PTC) whereby exons 6 to 9/10 are skipped (Cassette Exon, CE 6-9/10). We also demonstrated that the Northern blot analysis at transcript level is mediated via an intron-5 splicing event. Our data refine the structure of SCF gene; clarify the presence (+) and/or absence (-) of primary proteolytic-cleavage site specific SCF splice variants. This work provides a basis for understanding the functional role and regulation of SCF in hair follicle melanogenesis in sheep beyond what was known in mice, humans and other mammals.