Development of Poly(A)-ClickSeq as a tool enabling simultaneous genome-wide poly(A)-site identification and differential expression analysis

Regulation of the Drosophila transcriptome by Pumilio and the CCR4-NOT deadenylase complex

The sequence-specific RNA-binding protein Pumilio (Pum) controls Drosophila development; however, the network of mRNAs that it regulates remains incompletely characterized. In this study, we use knockdown and knockout approaches coupled with RNA-seq to measure the impact of Pum on the transcriptome of Drosophila cells in culture. We also use an improved RNA coimmunoprecipitation method to identify Pum-bound mRNAs in Drosophila embryos. Integration of these data sets with the locations of Pum-binding motifs across the transcriptome reveals novel direct Pum target genes involved in neural, muscle, wing, and germ cell development and in cellular proliferation. These genes include components of Wnt, TGF-β, MAPK/ERK, and Notch signaling pathways, DNA replication, and lipid metabolism. We identify the mRNAs regulated by the CCR4-NOT deadenylase complex, a key factor in Pum-mediated repression, and observe concordant regulation of Pum:CCR4-NOT target mRNAs. Computational modeling reveals that Pum binding, binding site number, clustering, and sequence context are important determinants of regulation. In contrast, we show that the responses of direct mRNA targets to Pum-mediated repression are not influenced by the content of optimal synonymous codons. Moreover, contrary to a prevailing model, we do not detect a role for CCR4-NOT in the degradation of mRNAs with low codon optimality. Together, the results of this work provide new insights into the Pum regulatory network and mechanisms and the parameters that influence the efficacy of Pum-mediated regulation.

Regulation of the Drosophila transcriptome by Pumilio and CCR4-NOT deadenylase

The sequence-specific RNA-binding protein Pumilio controls development of Drosophila; however, the network of mRNAs that it regulates remains incompletely characterized. In this study, we utilize knockdown and knockout approaches coupled with RNA-Seq to measure the impact of Pumilio on the transcriptome of Drosophila cells. We also used an improved RNA co-immunoprecipitation method to identify Pumilio bound mRNAs in Drosophila embryos. Integration of these datasets with the content of Pumilio binding motifs across the transcriptome revealed novel direct Pumilio target genes involved in neural, muscle, wing, and germ cell development, and cellular proliferation. These genes include components of Wnt, TGF-beta, MAPK/ERK, and Notch signaling pathways, DNA replication, and lipid metabolism. Additionally, we identified the mRNAs regulated by the CCR4-NOT deadenylase complex, a key factor in Pumilio-mediated repression, and observed concordant regulation of Pumilio:CCR4-NOT target mRNAs. Computational modeling revealed that Pumilio binding, binding site number, density, and sequence context are important determinants of regulation. Moreover, the content of optimal synonymous codons in target mRNAs exhibits a striking functional relationship to Pumilio and CCR4-NOT regulation, indicating that the inherent translation efficiency and stability of the mRNA modulates their response to these trans-acting regulatory factors. Together, the results of this work provide new insights into the Pumilio regulatory network and mechanisms, and the parameters that influence the efficacy of Pumilio-mediated regulation.

Sex chromosome mechanisms in cardiac development and disease

Many human diseases, including cardiovascular disease, show differences between men and women in pathology and treatment outcomes. In the case of cardiac disease, sex differences are exemplified by differences in the frequency of specific types of congenital and adult-onset heart disease. Clinical studies have suggested that gonadal hormones are a factor in sex bias. However, recent research has shown that gene and protein networks under non-hormonal control also account for cardiac sex differences. In this review, we describe the sex chromosome pathways that lead to sex differences in the development and function of the heart and highlight how these findings affect future care and treatment of cardiac disease.

Alternative polyadenylation alters protein dosage by switching between intronic and 3'UTR sites

Alternative polyadenylation (APA) creates distinct transcripts from the same gene by cleaving the pre-mRNA at poly(A) sites that can lie within the 3' untranslated region (3'UTR), introns, or exons. Most studies focus on APA within the 3'UTR; however, here, we show that CPSF6 insufficiency alters protein levels and causes a developmental syndrome by deregulating APA throughout the transcript. In neonatal humans and zebrafish larvae, CPSF6 insufficiency shifts poly(A) site usage between the 3'UTR and internal sites in a pathway-specific manner. Genes associated with neuronal function undergo mostly intronic APA, reducing their expression, while genes associated with heart and skeletal function mostly undergo 3'UTR APA and are up-regulated. This suggests that, under healthy conditions, cells toggle between internal and 3'UTR APA to modulate protein expression.

Macrophage differentiation is marked by increased abundance of the mRNA 3' end processing machinery, altered poly(A) site usage, and sensitivity to the level of CstF64

Regulation of mRNA polyadenylation is important for response to external signals and differentiation in several cell types, and results in mRNA isoforms that vary in the amount of coding sequence or 3' UTR regulatory elements. However, its role in differentiation of monocytes to macrophages has not been investigated. Macrophages are key effectors of the innate immune system that help control infection and promote tissue-repair. However, overactivity of macrophages contributes to pathogenesis of many diseases. In this study, we show that macrophage differentiation is characterized by shortening and lengthening of mRNAs in relevant cellular pathways. The cleavage/polyadenylation (C/P) proteins increase during differentiation, suggesting a possible mechanism for the observed changes in poly(A) site usage. This was surprising since higher C/P protein levels correlate with higher proliferation rates in other systems, but monocytes stop dividing after induction of differentiation. Depletion of CstF64, a C/P protein and known regulator of polyadenylation efficiency, delayed macrophage marker expression, cell cycle exit, attachment, and acquisition of structural complexity, and impeded shortening of mRNAs with functions relevant to macrophage biology. Conversely, CstF64 overexpression increased use of promoter-proximal poly(A) sites and caused the appearance of differentiated phenotypes in the absence of induction. Our findings indicate that regulation of polyadenylation plays an important role in macrophage differentiation.

Click Chemistry Enables Rapid Amplification of Full-Length Reverse Transcripts for Long-Read Third Generation Sequencing

Here we describe the development of a novel click chemistry-based method for the generation and amplification of full-length cDNA libraries from total RNA, while avoiding the need for problematic template-switching (TS) reactions. Compared with prior efforts, our method involves neither random priming nor stochastic cDNA termination, thus enabling amplification of transcripts that were previously inaccessible via related click chemistry-based RNA sequencing techniques. A key modification involving the use of PCR primers containing two overhanging 3'-nucleotides substantially improved the read-through compatibility of the 1,4-disubstituted 1,2,3-triazole-containing cDNA, where such modifications typically hinder amplification. This allowed us to more than double the possible insert size compared with the state-of-the art click chemistry-based technique, PAC-seq. Furthermore, our method performed on par with a commercially available PCR-cDNA RNA sequencing kit, as determined by Oxford Nanopore sequencing. Given the known advantages of PAC-seq, namely, suppression of PCR artifacts, we anticipate that our contribution could enable diverse applications including improved analyses of mRNA splicing variants and fusion transcripts.

Zika Virus Infection Alters Gene Expression and Poly-Adenylation Patterns in Placental Cells

Flaviviruses are small RNA viruses that are mainly transmitted via arthropod vectors and are found in tropic and sub-tropical regions. Most infections are asymptomatic (90-95%), but symptoms can be as severe as hemorrhagic fever and encephalitis. One recently emerged flavivirus is Zika virus (ZIKV), which was originally isolated from rhesus monkeys in Uganda roughly 70 years ago but has recently spread east, reaching S. America in 2015-2016. This outbreak was associated with the development of Guillain-Barré syndrome in adults and microcephaly in infants born to expectant mothers infected early in pregnancy. ZIKV must traverse the placenta to impact the development of the fetus, but the mechanisms responsible are unknown. While flaviviruses are known to disrupt splicing patterns in host cells, little is known about how flaviviruses such as ZIKV impact the alternative polyadenylation (APA) of host transcripts. This is important as APA is well-established as a mechanism in the regulation of mRNA metabolism and translation. Thus, we sought to characterize transcriptomic changes including APA in human placental (JEG3) cells in response to ZIKV infection using Poly(A)-ClickSeq (PAC-Seq). We used our differential Poly(A)-cluster (DPAC) analysis pipeline to characterize changes in differential gene expression, alternative poly-adenylation (APA) and the use of alternative terminal exons. We identified 98 upregulated genes and 28 downregulated genes. Pathway enrichment analysis indicated that many RNA processing and immune pathways were upregulated in ZIKV-infected JEG3 cells. We also updated DPAC to provide additional metrics of APA including the percentage-distal usage index (PDUI), which revealed that APA was extensive and the 3' UTRs of 229 genes were lengthened while 269 were shortened. We further found that there were 214 upregulated and 59 downregulated poly(A)-clusters (PACs). We extracted the nucleotide sequences surrounding these PACs and found that the canonical signals for poly-adenylation (binding site for poly-A binding protein (PABP) upstream and a GU-rich region down-stream of the PAC) were only enriched in the downregulated PACs. These results indicate that ZIKV infection makes JEG3 cells more permissive to non-canonical poly-adenylation signals.

QTQTN motif upstream of the furin-cleavage site plays a key role in SARS-CoV-2 infection and pathogenesis

The furin cleavage site (FCS), an unusual feature in the SARS-CoV-2 spike protein, has been spotlighted as a factor key to facilitating infection and pathogenesis by increasing spike processing. Similarly, the QTQTN motif directly upstream of the FCS is also an unusual feature for group 2B coronaviruses (CoVs). The QTQTN deletion has consistently been observed in in vitro cultured virus stocks and some clinical isolates. To determine whether the QTQTN motif is critical to SARS-CoV-2 replication and pathogenesis, we generated a mutant deleting the QTQTN motif (ΔQTQTN). Here, we report that the QTQTN deletion attenuates viral replication in respiratory cells in vitro and attenuates disease in vivo. The deletion results in a shortened, more rigid peptide loop that contains the FCS and is less accessible to host proteases, such as TMPRSS2. Thus, the deletion reduced the efficiency of spike processing and attenuates SARS-CoV-2 infection. Importantly, the QTQTN motif also contains residues that are glycosylated, and disruption of its glycosylation also attenuates virus replication in a TMPRSS2-dependent manner. Together, our results reveal that three aspects of the S1/S2 cleavage site-the FCS, loop length, and glycosylation-are required for efficient SARS-CoV-2 replication and pathogenesis.

Generation and characterization of human Fetal membrane and Decidual cell lines for reproductive biology experiments†

Human fetal membrane and maternal decidua parietalis form one of the major feto-maternal interfaces during pregnancy. Studies on this feto-maternal interface is limited as several investigators have limited access to the placenta, and experience difficulties to isolate and maintain primary cells. Many cell lines that are currently available do not have the characteristics or properties of their primary cells of origin. Therefore, we created, characterized the immortalized cells from primary isolates from fetal membrane-derived amnion epithelial cells, amnion and chorion mesenchymal cells, chorion trophoblast cells and maternal decidua parietalis cells. Primary cells were isolated from a healthy full-term, not in labor placenta. Primary cells were immortalized using either a HPV16E6E7 retroviral or a SV40T lentiviral system. The immortalized cells were characterized for the morphology, cell type-specific markers, and cell signalling pathway activation. Genomic stability of these cells was tested using RNA seq, karyotyping, and short tandem repeats DNA analysis. Immortalized cells show their characteristic morphology, and express respective epithelial, mesenchymal and decidual markers similar to that of primary cells. Gene expression of immortalized and primary cells were highly correlated (R = 0.798 to R = 0.974). Short tandem repeats DNA analysis showed in the late passage number (>P30) of cell lines matched 84-100% to the early passage number (<P10) of the cell lines revealing there were no genetic drift over the passages. Karyotyping also revealed no chromosomal anomalies. Creation of these cell lines can standardize experimental approaches, eliminate subject to subject variabilities, and benefit the reproductive biological studies on pregnancies by using these cells.

Tau Modulates mRNA Transcription, Alternative Polyadenylation Profiles of hnRNPs, Chromatin Remodeling and Spliceosome Complexes

Tau protein is a known contributor in several neurodegenerative diseases, including Alzheimer’s disease (AD) and frontotemporal dementia (FTD). It is well-established that tau forms pathological aggregates and fibrils in these diseases. Tau has been observed within the nuclei of neurons, but there is a gap in understanding regarding the mechanism by which tau modulates transcription. We are interested in the P301L mutation of tau, which has been associated with FTD and increased tau aggregation. Our study utilized tau-inducible HEK (iHEK) cells to reveal that WT and P301L tau distinctively alter the transcription and alternative polyadenylation (APA) profiles of numerous nuclear precursors mRNAs, which then translate to form proteins involved in chromatin remodeling and splicing. We isolated total mRNA before and after over-expressing tau and then performed Poly(A)-ClickSeq (PAC-Seq) to characterize mRNA expression and APA profiles. We characterized changes in Gene Ontology (GO) pathways using EnrichR and Gene Set Enrichment Analysis (GSEA). We observed that P301L tau up-regulates genes associated with reactive oxygen species responsiveness as well as genes involved in dendrite, microtubule, and nuclear body/speckle formation. The number of genes regulated by WT tau is greater than the mutant form, which indicates that the P301L mutation causes loss-of-function at the transcriptional level. WT tau up-regulates genes contributing to cytoskeleton-dependent intracellular transport, microglial activation, microtubule and nuclear chromatin organization, formation of nuclear bodies and speckles. Interestingly, both WT and P301L tau commonly down-regulate genes responsible for ubiquitin-proteosome system. In addition, WT tau significantly down-regulates several genes implicated in chromatin remodeling and nucleosome organization. Although there are limitations inherent to the model systems used, this study will improve understanding regarding the nuclear impact of tau at the transcriptional and post-transcriptional level. This study also illustrates the potential impact of P301L tau on the human brain genome during early phases of pathogenesis.

An atlas of alternative polyadenylation quantitative trait loci contributing to complex trait and disease heritability

Genome-wide association studies have identified thousands of noncoding variants associated with human traits and diseases. However, the functional interpretation of these variants is a major challenge. Here, we constructed a multi-tissue atlas of human 3'UTR alternative polyadenylation (APA) quantitative trait loci (3'aQTLs), containing approximately 0.4 million common genetic variants associated with the APA of target genes, identified in 46 tissues isolated from 467 individuals (Genotype-Tissue Expression Project). Mechanistically, 3'aQTLs can alter poly(A) motifs, RNA secondary structure and RNA-binding protein-binding sites, leading to thousands of APA changes. Our CRISPR-based experiments indicate that such 3'aQTLs can alter APA regulation. Furthermore, we demonstrate that mapping 3'aQTLs can identify APA regulators, such as La-related protein 4. Finally, 3'aQTLs are colocalized with approximately 16.1% of trait-associated variants and are largely distinct from other QTLs, such as expression QTLs. Together, our findings show that 3'aQTLs contribute substantially to the molecular mechanisms underlying human complex traits and diseases.

A computational pipeline to infer alternative poly-adenylation from 3' sequencing data

An increasing number of investigations have established alternative polyadenylation (APA) as a key mechanism of gene regulation through altering the length of 3' untranslated region (UTR) and generating distinct mRNA termini. Further, appreciation for the significance of APA in disease contexts propelled the development of several 3' sequencing techniques. While these RNA sequencing technologies have advanced APA analysis, the intrinsic limitation of 3' read coverage and lack of appropriate computational tools constrain precise mapping and quantification of polyadenylation sites. Notably, Poly(A)-ClickSeq (PAC-seq) overcomes limiting factors such as poly(A) enrichment and 3' linker ligation steps using click-chemistry. Here we provide an updated PolyA-miner protocol, a computational approach to analyze PAC-seq or other 3'-Seq datasets. As a key practical constraint, we also provide a detailed account on the impact of sequencing depth on the number of detected polyadenylation sites and APA changes. This protocol is also updated to handle unique molecular identifiers used to address PCR duplication potentially observed in PAC-seq.

SAM homeostasis is regulated by CFIm-mediated splicing of MAT2A

S-adenosylmethionine (SAM) is the methyl donor for nearly all cellular methylation events. Cells regulate intracellular SAM levels through intron detention of MAT2A, the only SAM synthetase expressed in most cells. The N6-adenosine methyltransferase METTL16 promotes splicing of the MAT2A detained intron by an unknown mechanism. Using an unbiased CRISPR knock-out screen, we identified CFIm25 (NUDT21) as a regulator of MAT2A intron detention and intracellular SAM levels. CFIm25 is a component of the cleavage factor Im (CFIm) complex that regulates poly(A) site selection, but we show it promotes MAT2A splicing independent of poly(A) site selection. CFIm25-mediated MAT2A splicing induction requires the RS domains of its binding partners, CFIm68 and CFIm59 as well as binding sites in the detained intron and 3´ UTR. These studies uncover mechanisms that regulate MAT2A intron detention and reveal a previously undescribed role for CFIm in splicing and SAM metabolism.

Application and design considerations for 3'-end sequencing using click-chemistry

Over the past 15 years, investigations into alternative polyadenylation (APA) and its function in cellular physiology and pathology have greatly expanded due to the emergent appreciation of its key role in driving transcriptomic diversity. This growth has necessitated the development of new technologies capable of monitoring cleavage and polyadenylation events genome-wide. Advancements in approaches include both the creation of computational tools to re-analyze RNA-seq to identify APA events as well as targeted sequencing approaches customized to focus on the 3'-end of mRNA. Here we describe a streamlined protocol for polyA-Click-seq (PAC-seq), which utilizes click-chemistry to create mRNA 3'-ends sequencing libraries. Importantly, we offer additional considerations not present in our previous protocols including the use of spike-ins, unique molecular identifier primers, and guidance for appropriate depth of PAC-seq. In conjunction with the companion chapter on PolyA-miner (Yalamanchili et al., 2021) to computationally analyze PAC-seq data, we provide a complete experimental pipeline to analyze mRNA 3'-end usage in eukaryotic cells.

Human Pumilio proteins directly bind the CCR4-NOT deadenylase complex to regulate the transcriptome

Pumilio paralogs, PUM1 and PUM2, are sequence-specific RNA-binding proteins that are essential for vertebrate development and neurological functions. PUM1&2 negatively regulate gene expression by accelerating degradation of specific mRNAs. Here, we determined the repression mechanism and impact of human PUM1&2 on the transcriptome. We identified subunits of the CCR4-NOT (CNOT) deadenylase complex required for stable interaction with PUM1&2 and to elicit CNOT-dependent repression. Isoform-level RNA sequencing revealed broad coregulation of target mRNAs through the PUM-CNOT repression mechanism. Functional dissection of the domains of PUM1&2 identified a conserved amino-terminal region that confers the predominant repressive activity via direct interaction with CNOT. In addition, we show that the mRNA decapping enzyme, DCP2, has an important role in repression by PUM1&2 amino-terminal regions. Our results support a molecular model of repression by human PUM1&2 via direct recruitment of CNOT deadenylation machinery in a decapping-dependent mRNA decay pathway.

PolyA-miner: accurate assessment of differential alternative poly-adenylation from 3'Seq data using vector projections and non-negative matrix factorization

Almost 70% of human genes undergo alternative polyadenylation (APA) and generate mRNA transcripts with varying lengths, typically of the 3′ untranslated regions (UTR). APA plays an important role in development and cellular differentiation, and its dysregulation can cause neuropsychiatric diseases and increase cancer severity. Increasing awareness of APA’s role in human health and disease has propelled the development of several 3′ sequencing (3′Seq) techniques that allow for precise identification of APA sites. However, despite the recent data explosion, there are no robust computational tools that are precisely designed to analyze 3′Seq data. Analytical approaches that have been used to analyze these data predominantly use proximal to distal usage. With about 50% of human genes having more than two APA isoforms, current methods fail to capture the entirety of APA changes and do not account for non-proximal to non-distal changes. Addressing these key challenges, this study demonstrates PolyA-miner, an algorithm to accurately detect and assess differential alternative polyadenylation specifically from 3′Seq data. Genes are abstracted as APA matrices, and differential APA usage is inferred using iterative consensus non-negative matrix factorization (NMF) based clustering. PolyA-miner accounts for all non-proximal to non-distal APA switches using vector projections and reflects precise gene-level 3′UTR changes. It can also effectively identify novel APA sites that are otherwise undetected when using reference-based approaches. Evaluation on multiple datasets—first-generation MicroArray Quality Control (MAQC) brain and Universal Human Reference (UHR) PolyA-seq data, recent glioblastoma cell line NUDT21 knockdown Poly(A)-ClickSeq (PAC-seq) data, and our own mouse hippocampal and human stem cell-derived neuron PAC-seq data—strongly supports the value and protocol-independent applicability of PolyA-miner. Strikingly, in the glioblastoma cell line data, PolyA-miner identified more than twice the number of genes with APA changes than initially reported. With the emerging importance of APA in human development and disease, PolyA-miner can significantly improve data analysis and help decode the underlying APA dynamics.

Partial loss of CFIm25 causes learning deficits and aberrant neuronal alternative polyadenylation

We previously showed that NUDT21-spanning copy-number variations (CNVs) are associated with intellectual disability (Gennarino et al., 2015). However, the patients' CNVs also included other genes. To determine if reduced NUDT21 function alone can cause disease, we generated Nudt21+/- mice to mimic NUDT21-deletion patients. We found that although these mice have 50% reduced Nudt21 mRNA, they only have 30% less of its cognate protein, CFIm25. Despite this partial protein-level compensation, the Nudt21+/- mice have learning deficits, cortical hyperexcitability, and misregulated alternative polyadenylation (APA) in their hippocampi. Further, to determine the mediators driving neural dysfunction in humans, we partially inhibited NUDT21 in human stem cell-derived neurons to reduce CFIm25 by 30%. This induced APA and protein level misregulation in hundreds of genes, a number of which cause intellectual disability when mutated. Altogether, these results show that disruption of NUDT21-regulated APA events in the brain can cause intellectual disability.

DPAC: A Tool for Differential Poly(A)-Cluster Usage from Poly(A)-Targeted RNAseq Data

Poly(A)-tail targeted RNAseq approaches, such as 3′READS, PAS-Seq and Poly(A)-ClickSeq, are becoming popular alternatives to random-primed RNAseq to focus sequencing reads just to the 3′ ends of polyadenylated RNAs to identify poly(A)-sites and characterize changes in their usage. Additionally, we and others have demonstrated that these approaches perform similarly to other RNAseq strategies for differential gene expression analysis, while saving on the volume of sequencing data required and providing a simpler library synthesis strategy. Here, we present DPAC (Differential Poly(A)-Clustering); a streamlined pipeline for the preprocessing of poly(A)-tail targeted RNAseq data, mapping of poly(A)-sites, poly(A)-site clustering and annotation, and determination of differential poly(A)-cluster usage using DESeq2. Changes in poly(A)-cluster usage is simultaneously used to report differential gene expression, differential terminal exon usage and alternative polyadenylation (APA).