Alternative splicing of the sheep MITF gene: novel transcripts detectable in skin

Genetic mapping identifies SNP mutations in MITF-M promoter associated with melanin formation in Putian black duck

The coloration of plumage in poultry species has substantial economic significance. Putian black ducks encompass 2 distinct strains characterized by black and white plumage variations resulting from selective breeding. This study aimed to identify the molecular mechanisms responsible for plumage coloration in these 2 distinct strains. A comprehensive genome-wide association study was conducted using DNA data sourced from a F2 segregating population, consisting of 71 individuals with black plumage and 39 individuals with white plumage, derived from these distinct 2 strains. This analysis revealed that 894 nucleotide polymorphisms and identified 58 candidate genes. Subsequent Gene Ontology and Kyoto Encyclopedia of Genes and Genomes coenrichment analyses identified MITF as a key candidate gene implicated in melanin biosynthesis. Furthermore, extensive screening of significant polymorphic loci within MITF was carried out via mass spectrometry in 3 distinct populations: 100 individuals with black plumage and 100 individuals with white plumage from the F0 generation; and 50 with black plumage form the F1 generation). Eighteen candidate polymorphic loci were identified demonstrating significant associations with variations in black and white plumage. Notably, 8 of these loci were located within the 2,000 bp region upstream of MITF-M. To validate the critical regulatory role of MITF-M in black and white plumage formation, a dual-fluorescence reporter system was constructed, and dual-fluorescence activity was assessed. The results revealed that the fluorescence activity at wild-type sites (corresponding to black plumage) was significantly higher than that at the mutant-type sites (associated with white plumage) (P < 0.01). To corroborate the pivotal role of MITF-M in black and white plumage formation, qPCR was employed to evaluate the expression levels of various MITF variants in black and white feather bulbs. This analysis demonstrated that only MITF-M exhibited specific expression in black feather bulbs. These results elucidate the central role of polymorphic mutations within the MITF promoter region in the regulation of black and white plumage coloration in Putian black ducks. This study extends our understanding of mechanisms governing duck plumage coloration and provides valuable molecular markers for future research in duck production and breeding based on plumage coloration.

Whole genome sequencing revealed genetic diversity, population structure, and selective signature of Panou Tibetan sheep

BACKGROUND

The detection of selective traits in different populations can not only reveal current mechanisms of artificial selection for breeding, but also provide new insights into phenotypic variation in new varieties and the search for genes associated with important traits. Panou sheep is a cultivated breed of Tibetan sheep in China with stable genetic performance, consistent appearance and fast growth and development after decades of artificial selection and cultivation. Due to long-term adaptation to the high altitude, cold and hypoxic environment in the plateau area, they may have formed a unique gene pool that is different from other Tibetan sheep breeds. To explore the genetic resources of Panou sheep, we used next-generation sequencing technology for the first time to investigate the genome-wide population structure, genetic diversity, and candidate signatures of positive selection in Panou sheep.

RESULTS

Comparative genomic analysis with the closely related species Oula sheep (a native breed of Tibetan sheep in China) was used to screen the population selection signal of Panou sheep. Principal component analysis and neighbor joining tree showed that Panou sheep and Oula sheep had differences in population differentiation. Furthermore, analyses of population structure, they came from the same ancestor, and when K = 2, the two populations could be distinguished. Panou sheep exhibit genetic diversity comparable to Oula sheep, as shown by observed heterozygosity, expected heterozygosity and runs of homozygosity. Genome-wide scanning using the Fst and π ratio methods revealed a list of potentially selected related genes in Panou sheep compared to Oula sheep, including histone deacetylase 9 (HDAC9), protein tyrosine kinase 2 (PTK2), microphthalmia-related transcription factor (MITF), vesicular amine transporter 1 (VAT1), trichohyalin-like 1 (TCHHL1), amine oxidase, copper containing 3 (AOC3), interferon-inducible protein 35 (IFI35).

CONCLUSIONS

The results suggest that traits related to growth and development and plateau adaptation may be selection targets for the domestication and breeding improvement of Tibetan sheep. This study provides the fundamental footprints for Panou sheep breeding and management.

Embryonic transcriptome unravels mechanisms and pathways underlying embryonic development with respect to muscle growth, egg production, and plumage formation in native and broiler chickens

Background: Muscle development, egg production, and plumage colors are different between native and broiler chickens. The study was designed to investigate why improved Aseel (PD4) is colorful, stronger, and grew slowly compared with the control broiler (CB).

Methods: A microarray was conducted using the 7th-day embryo (7EB) and 18th-day thigh muscle (18TM) of improved Aseel and broiler, respectively. Also, we have selected 24 Gallus gallus candidate reference genes from NCBI, and total RNA was isolated from the broiler, improved Aseel embryo tissues, and their expression profiles were studied by real-time quantitative PCR (qPCR). Furthermore, microarray data were validated with qPCR using improved Aseel and broiler embryo tissues.

Results: In the differential transcripts screening, all the transcripts obtained by microarray of slow and fast growth groups were screened by fold change ≥ 1 and false discovery rate (FDR) ≤ 0.05. In total, 8,069 transcripts were differentially expressed between the 7EB and 18TM of PD4 compared to the CB. A further analysis showed that a high number of transcripts are differentially regulated in the 7EB of PD4 (6,896) and fewer transcripts are differentially regulated (1,173) in the 18TM of PD4 compared to the CB. On the 7th^- and 18th-day PD4 embryos, 3,890, 3,006, 745, and 428 transcripts were up- and downregulated, respectively. The commonly up- and downregulated transcripts are 91 and 44 between the 7th- and 18th-day of embryos. In addition, the best housekeeping gene was identified. Furthermore, we validated the differentially expressed genes (DEGs) related to muscle growth, myostatin signaling and development, and fatty acid metabolism genes in PD4 and CB embryo tissues by qPCR, and the results correlated with microarray expression data.

Conclusion: Our study identified DEGs that regulate the myostatin signaling and differentiation pathway; glycolysis and gluconeogenesis; fatty acid metabolism; Jak-STAT, mTOR, and TGF-β signaling pathways; tryptophan metabolism; and PI3K-Akt signaling pathways in PD4. The results revealed that the gene expression architecture is present in the improved Aseel exhibiting embryo growth that will help improve muscle development, differentiation, egg production, protein synthesis, and plumage formation in PD4 native chickens. Our findings may be used as a model for improving the growth in Aseel as well as optimizing the growth in the broiler.

Genomic Uniqueness of Local Sheep Breeds From Morocco

Sheep farming is a major source of meat in Morocco and plays a key role in the country's agriculture. This study aims at characterizing the whole-genome diversity and demographic history of the main Moroccan sheep breeds, as well as to identify selection signatures within and between breeds. Whole genome data from 87 individuals representing the five predominant local breeds were used to estimate their level of neutral genetic diversity and to infer the variation of their effective population size over time. In addition, we used two methods to detect selection signatures: either for detecting selective sweeps within each breed separately or by detecting differentially selected regions by contrasting different breeds. We identified hundreds of genomic regions putatively under selection, which related to several biological terms involved in local adaptation or the expression of zootechnical performances such as Growth, UV protection, Cell maturation or Feeding behavior. The results of this study revealed selection signatures in genes that have an important role in traits of interest and increased our understanding of how genetic diversity is distributed in these local breeds. Thus, Moroccan local sheep breeds exhibit both a high genetic diversity and a large set of adaptive variations, and therefore, represent a valuable genetic resource for the conservation of sheep in the context of climate change.

Exon-1 skipping and intron-1 retaining by alternative splicing of the c-KIT gene encodes a novel splice variant in the skin of Merino sheep (Ovis aries)

c-KIT, a type III receptor protein tyrosine kinase, plays an essential role in melanocyte development, migration, and survival. Mutations within the c-KIT gene are previously shown to cause the white coat color phenotypes in pigs, mice, goats, and humans. However, up so far, the splicing isoform(s), genomic architecture of c-KIT have not been characterized well in merino sheep. Reverse transcriptase (RT)-PCR analysis with molecular prediction identified two basic splice variants: Transcript Variant-1, 2 for 12 bp insertion coding sequences (CDS) corresponding to the four amino acids 'GNSK', respectively. Using 5' RACE, here we report for the first time a novel c-KIT 'Transcript Variant-3' from the skin of merino sheep by comparative genome analyses at exon(1)-intron(1)-exon(2) boundaries. In contrast, a single product of 795 bp was characterized by 3' RACE. We also demonstrated that the c-KIT gene expression at the transcript level is not mediated via an intron-9 splicing event. Overall, beyond what was observed in other mammals, our data provide novel insights into the molecular structure of the c-KIT gene in sheep.

Transcriptome analysis identifies candidate genes associated with skin color variation in Triplophysa siluroides

In vertebrates, skin pigmentation is the most diverse phenotypic trait, and it is produced by a complex biological process that is often genetically controlled. Recently, two different colors (the typical brown and orange varieties) of Triplophysa siluroides, a species restricted to Yellow River drainage in China, were discovered. In the present study, the skin, brain and liver transcriptomes of T. siluroides of both colors were sequenced to search for genes related to skin pigmentation. Transcriptome sequencing generated 1,484,197,774 clean reads, resulting in a total of 222.6 Gb of sequence. The reads were assembled into 470,788 unigenes with a mean length of 1550 bp and an N50 size of 2944 bp. Functional annotation of the unigene dataset showed that 214,507, 304,161, 112,886, 179,074, 180,064, 184,837 and 82,081 unigenes were significantly matched to entries in the Nr protein, Nt, KO, Swiss-Prot, Pfam, GO and KOG databases, respectively. A differential expression analysis revealed that 2774, 3552 and 1529 unigenes were upregulated and 2720, 2663 and 1103 unigenes were downregulated in the skin, brain and liver of orange-skinned T. siluroides, respectively. Several genes that play key roles in pigmentation, i.e., Agouti, Slc45a2, Cbs, Mift and Slc7a11, showed significantly differential expression between brown and orange fish. In addition, we detected 158,863 simple sequence repeats (SSRs) in the T. siluroides transcriptome, and a total of 201,338 single-nucleotide polymorphisms (SNPs) were discovered in the different transcriptomes. The present results will facilitate further study of the molecular mechanisms of skin pigmentation and marker-assisted breeding of fish with valuable skin colors.

TYR Gene in Llamas: Polymorphisms and Expression Study in Different Color Phenotypes

Tyrosinase, encoded by TYR gene, is an enzyme that plays a major role in mammalian pigmentation. It catalyzes the oxidation of L-dihydroxy-phenylalanine (DOPA) to DOPA quinone, a precursor of both types of melanin: eumelanin and pheomelanin. TYR is commonly known as the albino locus since mutations in this gene result in albinism in several species. However, many other TYR mutations have been found to cause diluted phenotypes, like the Himalayan or chinchilla phenotypes in mice. The llama (Lama glama) presents a wide variety of coat colors ranging from non-diluted phenotypes (eumelanic and pheomelanic), through different degrees of dilution, to white. To investigate the possible contribution of TYR gene to coat color variation in llamas, we sequenced TYR exons and their flanking regions and genotyped animals with diluted, non-diluted, and white coat, including three blue-eyed white individuals. Moreover, we analyzed mRNA expression levels in skin biopsies by qPCR. TYR coding region presented nine SNPs, of which three were non-synonymous, c.428A > G, c.859G > T, and c.1490G > T. We also identified seven polymorphisms in non-coding regions, including two microsatellites, an homopolymeric repeat, and five SNPs: one in the promoter region (c.1-26C > T), two in the 3'-UTR, and two flanking the exons. Although no complete association was found between coat color and SNPs, c.1-26C > T was partially associated to diluted phenotypes. Additionally, the frequency of the G allele from c.428A > G was significantly higher in white compared to non-diluted. Results from qPCR showed that expression levels of TYR in white llamas were significantly lower (p < 0.05) than those in diluted and non-diluted phenotypes. Screening for variation in regulatory regions of TYR did not reveal polymorphisms that explain such differences. However, data from this study showed that TYR expression levels play a role in llama pigmentation.

Characterization and expression analysis of KIT and MITF-M genes in llamas and their relation to white coat color

The llama (Lama glama) is a fiber-producing species that presents a wide range of coat colors, among which white is one of the most important for the textile industry. However, there is little information about the molecular mechanisms that control the white phenotype in this species. In domestic mammals, a white coat is usually produced by mutations in the KIT proto-oncogene receptor tyrosine kinase (KIT) and microphthalmia-associated transcription factor (MITF) genes. In this work we have sequenced and described the coding regions of KIT and MITF-M, the melanocyte-specific isoform, and the two transcriptional variants MITF-M(-) and MITF-M(+). Moreover, we studied the expression of these genes in the skin of white and colored llamas. Although no variants were revealed to be associated with white coat color, significant differences between phenotypes were observed in the expression levels of KIT and MITF-M. Interestingly, white llamas expressed less MITF-M(+) than did colored ones, which is consistent with a consequent reduction in the synthesis of melanin. Even though our results indicate that downregulation of KIT and MITF-M expression is involved in white phenotype production in llamas, the causative gene of white coat color remains unknown.

A wide variety of Mitf transcript variants are expressed in the Xenopus laevis periodic albino mutant

The periodic albino mutant of Xenopus laevis has been used to study the development of pigment cells because both the retinal pigment epithelium (RPE) and melanophores are affected. In this mutant, "white pigment cells" containing both melanophore-specific and iridophore-specific pigment organelles appear. The present experiments were designed to investigate the structural organization and expression of microphthalmia-associated transcription factor (Mitf) in the mutant since Mitf is known to regulate the development of melanocytes and RPE. The exon structures of X. laevis mitf genes (mitf.L and mitf.S) were defined using newly obtained Mitf transcripts and X. laevis genomic data. Compared to mouse mitf, exons 3 and 6a were absent in X. laevis mitf. The four exons between exons 4 and 6b in X. laevis mitf were named 5α, 5β, 5γ, and 5δ. Exons 5α and 5δ were specific to X. laevis mitf, whereas the continuous exons 5β/γ were identical to exon 5 of mouse mitf. A wide variety of A-Mitf and M-Mitf transcript variants lacking one or more exons were found in X. laevis; however, different types of Mitf transcripts were expressed in the mutant. In addition, white pigment cells and melanophores expressed both the mitf and dopachrome tautomerase (dct) genes.

Identification of novel alternative splicing transcript and expression analysis of bovine TMEM95 gene

Transmembrane protein 95 (TMEM95) is closely related to male reproductive performance in cattle, but does not affect semen quality. Alternative splicing plays an important role in regulating biological function as well as in generating proteomic and functional diversity in metazoan organisms. Thus, the aim of this study was to clone and identify transcripts of the TMEM95 gene in cattle using RT-PCR, characterize them via bioinformatics analysis, and detect their expression patterns using qRT-PCR. Two transcripts of TMEM95 were identified in cattle, including TMEM95-SV1 and TMEM95-SV2. Bioinformatics predicted that TMEM95-SV1 has a leucine-rich repeat C-terminal domain and a Pfam: IZUMO. These regions are closely related to protein interactions and the acrosome reaction, respectively. Interestingly, the two transcripts were exclusively expressed in the testes and brain in male fetus cattle, and TMEM95-SV1 was expressed in the brain at significantly higher levels than in the testis (P<0.05, 4.06-fold) and TMEM95-SV2 in the brain (P<0.05, 4.95-fold). These findings enrich the understanding of the TMEM95 gene function and benefit for enhancing male reproduction in cattle industry.

Skin transcriptome profiles associated with skin color in chickens

Nutritional and medicinal benefits have been attributed to the consumption of tissues from the black-boned chickens in oriental countries. Lueyang black-boned chicken is one of the native chicken breeds. However, some birds may instead have white or lighter skin, which directly causes economic losses every year. Previous studies of pigmentation have focused on a number of genes that may play important roles in coat color regulation. Illumina2000 sequencing technology was used to catalog the global gene expression profiles in the skin of the Lueyang chicken with white versus black skin. A total of 18,608 unigenes were assembled from the reads obtained from the skin of the white and black chickens. A total of 649 known genes were differentially expressed in the black versus white chickens, with 314 genes that were up regulated and 335 genes that were down-regulated, and a total of 162 novel genes were differentially expressed in the black versus white chickens, consisting of 73 genes that were up-regulated (including 4 highly expressed genes that were expressed exclusively in the skin of the black chickens) and 89 genes that were down-regulated. There were also a total of 8 known coat-color genes expressed in previous studies (ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R). In this study, 4 of which showed greater expression in the black chickens, and several were up-regulated, such as KIT, ASIP, TYR and OCA2. To our surprise, KITLG, MITF and MC1R showed no significant difference in expression between the black- and white-skinned chickens, and the expression of TYRP1 was not detected in either skin color. The expression of ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R was validated by real-time quantitative polymerase chain reaction (qPCR), and the results of the qPCR were consistent with the RNA-seq. This study provides several candidate genes that may be associated with the development of black versus white skin. More importantly, the fact that the MC1R gene showed no significant difference in expression between the black and white chickens is of particular interest for future studies that aim to elucidate its functional role in the regulation of skin color.