The ASIP gene in the llama (Lama glama): Alternative transcripts, expression and relation with color phenotypes

•Seven ASIP transcripts with alternative 5′UTR were isolated.•A fusion transcript between ASIP and NCOA6 gene was found.•The novel non-coding exon 1A” was identified within the 5́UTR regions.•Multiple start sites and alternative splicing create diversity of ASIP transcripts.•ASIP expression patterns differed among color phenotypes. The Agouti gene (ASIP) is one of the most important genes for coat color determination in mammals. It has a complex structure with several promoters and alternative non-coding first exons that are transcribed into mRNAs with different 5'UTR. These mRNA isoforms regulate the temporal and spatial expression of the gene, producing diverse pigmentation patterns. Here, we studied ASIP transcriptional variants and their expression in the skin of llamas with different coat color phenotypes. We also described the ASIP locus, including promoter usage and the splicing events that originate each transcript variant. Using 5'RACE-PCR we isolated seven ASIP transcripts with alternative 5′UTR, where exons 1A, 1A’, 1C, 1D, and a novel non-coding exon 1A” were identified. Additionally, new alternative spliced forms were found. The diversity of ASIP 5′UTRs is originated by a complex pattern of alternative promoter usage, multiple transcription start sites and splicing events that include exon skipping and alternative 3′ splicing site selection. We found that ASIP was highly expressed in llamas with white and brown phenotypes while black animals presented very low expression. The main responsible for this difference was a fusion transcript between ASIP and NCOA6 genes, which was present in the skin of white and brown llamas but not in the black ones. The rest of ASIP transcripts presented very low expression in the skin, indicating that the main regulation point for ASIP gene expression is at the transcriptional level. Nevertheless, the characteristics of the 5′UTRs sequences suggest that alternative transcripts could be regulated differently at the protein synthesis level.