Little is known about co-transcriptional or post-transcriptional regulatory mechanisms linking noncoding variation to variation in organismal traits. To begin addressing this gap, we used 3’ Seq to study the impact of genetic variation on alternative polyadenylation (APA) in the nuclear and total mRNA fractions of 52 HapMap Yoruba human lymphoblastoid cell lines. We mapped 602 APA quantitative trait loci (apaQTLs) at 10% FDR, of which 152 were nuclear specific. Effect sizes at intronic apaQTLs are negatively correlated with eQTL effect sizes. These observations suggest genetic variants can decrease mRNA expression levels by increasing usage of intronic PAS. We also identified 24 apaQTLs associated with protein levels, but not mRNA expression. Finally, we found that 19% of apaQTLs can be associated with disease. Thus, our work demonstrates that APA links genetic variation to variation in gene expression, protein expression, and disease risk, and reveals uncharted modes of genetic regulation.

While comparative functional genomic studies have shown that inter-species differences in gene expression can be explained by corresponding inter-species differences in genetic and epigenetic regulatory mechanisms, co-transcriptional mechanisms, such as alternative polyadenylation (APA), have received little attention. We characterized APA in lymphoblastoid cell lines from six humans and six chimpanzees by identifying and estimating usage for 44,432 polyadenylation sites (PAS) in 9,518 genes. Although APA is largely conserved, 1,705 genes showed significantly different PAS usage (FDR 0.05) between species. Genes with divergent APA also tend to be differentially expressed, are enriched among genes showing differences in protein translation, and can explain a subset of observed inter-species protein expression differences that do not differ at the transcript level. Finally, we found that genes with a dominant PAS, which is used more often than other PAS, are particularly enriched for differentially expressed genes.

With the exception of mRNA splicing, little is known about co-transcriptional or post-transcriptional regulatory mechanisms that link noncoding variation to variation in organismal traits. To begin addressing this gap, we used 3' Seq to characterize alternative polyadenylation (APA) in the nuclear and total RNA fractions of 52 HapMap Yoruba lymphoblastoid cell lines, which we have studied extensively in the past. We identified thousands of polyadenylation sites that are differentially detected in nuclear mRNA and whole cell mRNA, and found that APA is an important mediator of genetic effects on gene regulation and complex traits. Specifically, we mapped 602 apaQTLs at 10% FDR, of which 152 were found only in the nuclear fraction. Nuclear-specific apaQTLs are highly enriched in introns and are also often associated with changes in steady-state expression levels, suggesting a widespread mechanism whereby genetic variants decrease mRNA expression levels by increasing usage of intronic PAS. We identified 24 apaQTLs associated with protein expression levels, but not mRNA expression, and found that eQTLs that are not associated with chromatin QTLs are enriched in apaQTLs. These findings support multiple independent pathways through which genetic effects on APA can impact gene regulation. Finally, we found that 19% of apaQTLs were also previously associated with disease. Thus, our work demonstrates that APA links genetic variation to variation in gene expression levels, protein expression levels, and disease risk, and reveals uncharted modes of genetic regulation.

No abstract

No abstract