Polymorphism of alternative splicing of major histocompatibility complex transcripts in wild tiger salamanders

Alternative splicing (AS) of mRNA transcripts is increasingly recognized as a source of transcriptome diversity. To date, most AS studies have focused either on comparisons across taxa or on intragenomic comparisons across gene families. We generated a novel data set that represents one of the first population genetic comparisons of AS across individuals. In ambystomatid salamanders, AS of the major histocompatibility complex (MHC) class IIbeta gene (Amti-DAB) produces two transcripts, one full-length and one truncated. The full-length transcript is functional, but the truncated transcript is missing the critical beta1 domain that forms half of the peptide binding region in the intact MHC class II molecule. We captured wild salamander larvae (Ambystoma tigrinum tigrinum) and genotyped them at Amti-DAB via DNA sequencing. From these same larvae, we extracted RNA from gill and spleen and evaluated the relative expression level of Amti-DAB in each tissue. Across individuals, 21% of the transcripts were truncated (alternatively spliced), and the absolute level of alternative transcript expression was higher in gill. The high level of nucleotide variation among seven Amti-DAB alleles provides the ability to detect substitutions (or linked DNA polymorphisms) that might have influenced AS. The data reveal no correlation between AS and haplotype, allele, or zygosity. However, indirect evidence (comparative expression patterns across 3 million years of evolution) suggests that the truncated Amti-DAB transcript may be functional and maintained by natural selection.