Integrated multi-omics analysis of zinc-finger proteins uncovers roles in RNA regulation.

RNA interactome studies have revealed that hundreds of zinc-finger proteins (ZFPs) are candidate RNA-binding proteins (RBPs), yet their RNA substrates and functional significance remain largely uncharacterized. Here, we present a systematic multi-omics analysis of the DNA- and RNA-binding targets and regulatory roles of more than 100 ZFPs representing 37 zinc-finger families. We show that multiple ZFPs are previously unknown regulators of RNA splicing, alternative polyadenylation, stability, or translation. The examined ZFPs show widespread sequence-specific RNA binding and preferentially bind proximal to transcription start sites. Additionally, several ZFPs associate with their targets at both the DNA and RNA levels. We highlight ZNF277, a C2H2 ZFP that binds thousands of RNA targets and acts as a multi-functional RBP. We also show that ZNF473 is a DNA/RNA-associated protein that regulates the expression and splicing of cell cycle genes. Our results reveal diverse roles for ZFPs in transcriptional and post-transcriptional gene regulation.

Codon optimality is a major determinant of mRNA stability.

mRNA degradation represents a critical regulated step in gene expression. Although the major pathways in turnover have been identified, accounting for disparate half-lives has been elusive. We show that codon optimality is one feature that contributes greatly to mRNA stability. Genome-wide RNA decay analysis revealed that stable mRNAs are enriched in codons designated optimal, whereas unstable mRNAs contain predominately non-optimal codons. Substitution of optimal codons with synonymous, non-optimal codons results in dramatic mRNA destabilization, whereas the converse substitution significantly increases stability. Further, we demonstrate that codon optimality impacts ribosome translocation, connecting the processes of translation elongation and decay through codon optimality. Finally, we show that optimal codon content accounts for the similar stabilities observed in mRNAs encoding proteins with coordinated physiological function. This work demonstrates that codon optimization exists as a mechanism to finely tune levels of mRNAs and, ultimately, proteins.

Characterization and implementation of the MarathonRT template-switching reaction to expand the capabilities of RNA-Seq.

End-to-end RNA sequencing methods that capture 5'-sequence content without cumbersome library manipulations are of great interest, particularly for analysis of long RNAs. While template-switching methods have been developed for RNA sequencing by distributive short-read RTs, such as the MMLV RT enzymes used in SMART-Seq methods, they have not been adapted to leverage the power of ultraprocessive RTs, such as those that derive from group II self-splicing introns. To facilitate this transition, we dissected the individual processes that guide the enzymatic specificity and efficiency of the multi-step template switching reaction carried out by RT enzymes, in this case, by a well-characterized enzyme known as MarathonRT. Remarkably, this is the first study of its kind, for any RT. First, we characterized and optimized the enzymatic nontemplated addition (NTA) reaction that occurs when the RT enzyme extends past the RNA 5'-terminus, and we determined the nucleotide specificity of the NTA reaction. We then evaluated the binding specificity of specialized template-switching oligonucleotides, optimizing their sequences and chemical properties to guide efficient template switching reaction. Having dissected and optimized these individual steps, we then unified them into a procedure for performing RNA sequencing with MarathonRT enzymes, using a well-characterized RNA reference set. The resulting reads span a six-log range in transcript concentration and accurately represent the input RNA identities in both length and composition. We also performed RNA-seq starting from total human RNA and poly(A)-enriched RNA, with short and long-read sequencing demonstrating that MarathonRT enhances the discovery of unseen RNA molecules by conventional RT. Altogether, by employing mechanistic enzymology on RT enzymes and using them to modify RNA-seq technologies, we have generated a new pipeline for rapid, accurate sequencing of complex RNA libraries containing mixtures of long RNA transcripts.

A large-scale binding and functional map of human RNA-binding proteins.

Many proteins regulate the expression of genes by binding to specific regions encoded in the genome1. Here we introduce a new data set of RNA elements in the human genome that are recognized by RNA-binding proteins (RBPs), generated as part of the Encyclopedia of DNA Elements (ENCODE) project phase III. This class of regulatory elements functions only when transcribed into RNA, as they serve as the binding sites for RBPs that control post-transcriptional processes such as splicing, cleavage and polyadenylation, and the editing, localization, stability and translation of mRNAs. We describe the mapping and characterization of RNA elements recognized by a large collection of human RBPs in K562 and HepG2 cells. Integrative analyses using five assays identify RBP binding sites on RNA and chromatin in vivo, the in vitro binding preferences of RBPs, the function of RBP binding sites and the subcellular localization of RBPs, producing 1,223 replicated data sets for 356 RBPs. We describe the spectrum of RBP binding throughout the transcriptome and the connections between these interactions and various aspects of RNA biology, including RNA stability, splicing regulation and RNA localization. These data expand the catalogue of functional elements encoded in the human genome by the addition of a large set of elements that function at the RNA level by interacting with RBPs.

Late-life shift in caloric intake affects fly metabolism and longevity.

The prevalence of obesity is increasing in older adults and contributes to age-related decline. Caloric restriction (CR) alleviates obesity phenotypes and delays the onset of age-related changes. However, how late in life organisms benefit from switching from a high-(H) to a low-calorie (L) diet is unclear. We transferred male flies from a H to a L (HL) diet or vice versa (LH) at different times during life. Both shifts immediately change fly rate of aging even when applied late in life. HL shift rapidly reduces fly mortality rate to briefly lower rate than in flies on a constant L diet, and extends lifespan. Transcriptomic analysis uncovers that flies aged on H diet have acquired increased stress response, which may have temporal advantage over flies aged on L diet and leads to rapid decrease in mortality rate after HL switch. Conversely, a LH shift increases mortality rate, which is temporarily higher than in flies aged on a H diet, and shortens lifespan. Unexpectedly, more abundant transcriptomic changes accompanied LH shift, including increase in ribosome biogenesis, stress response and growth. These changes reflect protection from sudden release of ROS, energy storage, and use of energy to growth, which all likely contribute to higher mortality rate. As the beneficial effects of CR on physiology and lifespan are conserved across many organisms, our study provides framework to study underlying mechanisms of CR interventions that counteract the detrimental effects of H diets and reduce rate of aging even when initiated later in life.

RBP Image Database: A resource for the systematic characterization of the subcellular distribution properties of human RNA binding proteins.

RNA binding proteins (RBPs) are central regulators of gene expression implicated in all facets of RNA metabolism. As such, they play key roles in cellular physiology and disease etiology. Since different steps of post-transcriptional gene expression tend to occur in specific regions of the cell, including nuclear or cytoplasmic locations, defining the subcellular distribution properties of RBPs is an important step in assessing their potential functions. Here, we present the RBP Image Database, a resource that details the subcellular localization features of 301 RBPs in the human HepG2 and HeLa cell lines, based on the results of systematic immuno-fluorescence studies conducted using a highly validated collection of RBP antibodies and a panel of 12 markers for specific organelles and subcellular structures. The unique features of the RBP Image Database include: (i) hosting of comprehensive representative images for each RBP-marker pair, with ∼250,000 microscopy images; (ii) a manually curated controlled vocabulary of annotation terms detailing the localization features of each factor; and (iii) a user-friendly interface allowing the rapid querying of the data by target or annotation. The RBP Image Database is freely available at https://rnabiology.ircm.qc.ca/RBPImage/.

Type III-A CRISPR systems as a versatile gene knockdown technology.

CRISPR-Cas systems are functionally diverse prokaryotic antiviral defense systems, which encompass six distinct types (I-VI) that each encode different effector Cas nucleases with distinct nucleic acid cleavage specificities. By harnessing the unique attributes of the various CRISPR-Cas systems, a range of innovative CRISPR-based DNA and RNA targeting tools and technologies have been developed. Here, we exploit the ability of type III-A CRISPR-Cas systems to carry out RNA-guided and sequence-specific target RNA cleavage for establishment of research tools for post-transcriptional control of gene expression. Type III-A systems from three bacterial species (L. lactis, S. epidermidis, and S. thermophilus) were each expressed on a single plasmid in E. coli, and the efficiency and specificity of gene knockdown was assessed by northern blot and transcriptomic analysis. We show that engineered type III-A modules can be programmed using tailored CRISPR RNAs to efficiently knock down gene expression of both coding and noncoding RNAs in vivo. Moreover, simultaneous degradation of multiple cellular mRNA transcripts can be directed by utilizing a CRISPR array expressing corresponding gene-targeting crRNAs. Our results demonstrate the utility of distinct type III-A modules to serve as specific and effective gene knockdown platforms in heterologous cells. This transcriptome engineering technology has the potential to be further refined and exploited for key applications including gene discovery and gene pathway analyses in additional prokaryotic and perhaps eukaryotic cells and organisms.

RNA-guided RNA cleavage by a CRISPR RNA-Cas protein complex.

Compelling evidence indicates that the CRISPR-Cas system protects prokaryotes from viruses and other potential genome invaders. This adaptive prokaryotic immune system arises from the clustered regularly interspaced short palindromic repeats (CRISPRs) found in prokaryotic genomes, which harbor short invader-derived sequences, and the CRISPR-associated (Cas) protein-coding genes. Here, we have identified a CRISPR-Cas effector complex that is comprised of small invader-targeting RNAs from the CRISPR loci (termed prokaryotic silencing (psi)RNAs) and the RAMP module (or Cmr) Cas proteins. The psiRNA-Cmr protein complexes cleave complementary target RNAs at a fixed distance from the 3' end of the integral psiRNAs. In Pyrococcus furiosus, psiRNAs occur in two size forms that share a common 5' sequence tag but have distinct 3' ends that direct cleavage of a given target RNA at two distinct sites. Our results indicate that prokaryotes possess a unique RNA silencing system that functions by homology-dependent cleavage of invader RNAs.

Transcriptional regulation by NR5A2 links differentiation and inflammation in the pancreas.

Chronic inflammation increases the risk of developing one of several types of cancer. Inflammatory responses are currently thought to be controlled by mechanisms that rely on transcriptional networks that are distinct from those involved in cell differentiation. The orphan nuclear receptor NR5A2 participates in a wide variety of processes, including cholesterol and glucose metabolism in the liver, resolution of endoplasmic reticulum stress, intestinal glucocorticoid production, pancreatic development and acinar differentiation. In genome-wide association studies, single nucleotide polymorphisms in the vicinity of NR5A2 have previously been associated with the risk of pancreatic adenocarcinoma. In mice, Nr5a2 heterozygosity sensitizes the pancreas to damage, impairs regeneration and cooperates with mutant Kras in tumour progression. Here, using a global transcriptomic analysis, we describe an epithelial-cell-autonomous basal pre-inflammatory state in the pancreas of Nr5a2+/- mice that is reminiscent of the early stages of pancreatitis-induced inflammation and is conserved in histologically normal human pancreases with reduced expression of NR5A2 mRNA. In Nr5a2+/-mice, NR5A2 undergoes a marked transcriptional switch, relocating from differentiation-specific to inflammatory genes and thereby promoting gene transcription that is dependent on the AP-1 transcription factor. Pancreatic deletion of Jun rescues the pre-inflammatory phenotype, as well as binding of NR5A2 to inflammatory gene promoters and the defective regenerative response to damage. These findings support the notion that, in the pancreas, the transcriptional networks involved in differentiation-specific functions also suppress inflammatory programmes. Under conditions of genetic or environmental constraint, these networks can be subverted to foster inflammation.

Distal Alternative Last Exons Localize mRNAs to Neural Projections.

Spatial restriction of mRNA to distinct subcellular locations enables local regulation and synthesis of proteins. However, the organizing principles of mRNA localization remain poorly understood. Here we analyzed subcellular transcriptomes of neural projections and soma of primary mouse cortical neurons and two neuronal cell lines and found that alternative last exons (ALEs) often confer isoform-specific localization. Surprisingly, gene-distal ALE isoforms were four times more often localized to neurites than gene-proximal isoforms. Localized isoforms were induced during neuronal differentiation and enriched for motifs associated with muscleblind-like (Mbnl) family RNA-binding proteins. Depletion of Mbnl1 and/or Mbnl2 reduced localization of hundreds of transcripts, implicating Mbnls in localization of mRNAs to neurites. We provide evidence supporting a model in which the linkage between genomic position of ALEs and subcellular localization enables coordinated induction of localization-competent mRNA isoforms through a post-transcriptional regulatory program that is induced during differentiation and reversed in cellular reprogramming and cancer.

Resources for the Comprehensive Discovery of Functional RNA Elements.

Transcriptome-wide maps of RNA binding protein (RBP)-RNA interactions by immunoprecipitation (IP)-based methods such as RNA IP (RIP) and crosslinking and IP (CLIP) are key starting points for evaluating the molecular roles of the thousands of human RBPs. A significant bottleneck to the application of these methods in diverse cell lines, tissues, and developmental stages is the availability of validated IP-quality antibodies. Using IP followed by immunoblot assays, we have developed a validated repository of 438 commercially available antibodies that interrogate 365 unique RBPs. In parallel, 362 short-hairpin RNA (shRNA) constructs against 276 unique RBPs were also used to confirm specificity of these antibodies. These antibodies can characterize subcellular RBP localization. With the burgeoning interest in the roles of RBPs in cancer, neurobiology, and development, these resources are invaluable to the broad scientific community. Detailed information about these resources is publicly available at the ENCODE portal (https://www.encodeproject.org/).

Direct tracking of reverse-transcriptase speed and template sensitivity: implications for sequencing and analysis of long RNA molecules.

Although reverse-transcriptase (RT) enzymes are critical reagents for research and biotechnology, their mechanical properties are not well understood. In particular, we know little about their relative speed and response to structural obstacles in the template. Commercial retroviral RTs stop at many positions along mixed sequence templates, resulting in truncated cDNA products that complicate downstream analysis. By contrast, group II intron-encoded RTs appear to copy long RNAs with high processivity and minimal stops. However, their speed, consistency and pausing behavior have not been explored. Here, we analyze RT velocity as the enzyme moves through heterogeneous sequences and structures that are embedded within a long noncoding RNA transcript. We observe that heterogeneities in the template are highly disruptive to primer extension by retroviral RTs. However, sequence composition and template structure have negligible effects on behavior of group II intron RTs, such as MarathonRT (MRT). Indeed, MRT copies long RNAs in a single pass, and displays synchronized primer extension at a constant speed of 25 nt/sec. In addition, it passes through stable RNA structural motifs without perturbation of velocity. Taken together, the results demonstrate that consistent, robust translocative behavior is a hallmark of group II intron-encoded RTs, some of which operate at high velocity.

Investigation of CRISPR-Independent Phage Resistance Mechanisms Reveals a Role for FtsH in Phage Adsorption to Streptococcus thermophilus.

Prokaryotes are under constant pressure from phage infection and thus have evolved multiple means of defense or evasion. While CRISPR-Cas constitutes a robust immune system and appears to be the predominant means of survival for Streptococcus thermophilus when facing lytic phage infection, other forms of phage resistance coexist in this species. Here, we show that S. thermophilus strains with deleted CRISPR-Cas loci can still give rise to phage-resistant clones following lytic phage challenge. Notably, non-CRISPR phage-resistant survivors had multiple mutations which would truncate or recode a membrane-anchored host protease, FtsH. Phage adsorption was dramatically reduced in FtsH mutants, implicating this protein in phage attachment. Phages were isolated which could bypass FtsH-based resistance through mutations predicted to alter tape measure protein translation. Together, these results identify key components in phage propagation that are subject to mutation in the molecular arms race between phage and host cell. IMPORTANCE Streptococcus thermophilus is an important organism for production of cultured dairy foods, but it is susceptible to lytic phages which can lead to failed products. Consequently, mechanisms for phage resistance are an active area of research. One such mechanism is CRISPR-Cas, and S. thermophilus is a model organism for the study of this form of adaptive immunity. Here, we expand on known mechanisms with our finding that spontaneous mutations in ftsH, a gene encoding a membrane-anchored protease, protected against phage infection by disrupting phage adsorption. In turn, mutations in phage tail protein genes allowed phages to overcome ftsH-based resistance. Our results identified components in phage propagation that are subject to mutation in the molecular arms race between phage and host.

Loss of the seipin gene perturbs eggshell formation in Caenorhabditiselegans.

Seipin, an evolutionary conserved protein, plays pivotal roles during lipid droplet (LD) biogenesis and is associated with various human diseases with unclear mechanisms. Here, we analyzed Caenorhabditis elegans mutants deleted of the sole SEIPIN gene, seip-1 Homozygous seip-1 mutants displayed penetrant embryonic lethality, which is caused by the disruption of the lipid-rich permeability barrier, the innermost layer of the C. elegans embryonic eggshell. In C. elegans oocytes and embryos, SEIP-1 is associated with LDs and is crucial for controlling LD size and lipid homeostasis. The seip-1 deletion mutants reduced the ratio of polyunsaturated fatty acids (PUFAs) in their embryonic fatty acid pool. Interestingly, dietary supplementation of selected n-6 PUFAs rescued the embryonic lethality and defective permeability barrier. Accordingly, we propose that SEIP-1 may maternally regulate LD biogenesis and lipid homeostasis to orchestrate the formation of the permeability barrier for eggshell synthesis during embryogenesis. A lipodystrophy allele of seip-1 resulted in embryonic lethality as well and could be rescued by PUFA supplementation. These experiments support a great potential for using C. elegans to model SEIPIN-associated human diseases.

Risk factors for endometrial cancer in Black women.

PURPOSE: The incidence of endometrial cancer (EC) has been increasing faster among Black women than among other racial/ethnic groups in the United States. Although the mortality rate is nearly twice as high among Black than White women, there is a paucity of literature on risk factors for EC among Black women, particularly regarding menopausal hormone use and severe obesity. METHODS: We pooled questionnaire data on 811 EC cases and 3,124 controls from eight studies with data on self-identified Black women (4 case-control and 4 cohort studies). We analyzed cohort studies as nested case-control studies with up to 4 controls selected per case. We used logistic regression to estimate multivariable-adjusted odds ratios (ORs) and 95% confidence intervals (CIs). RESULTS: We observed a positive association between BMI and EC incidence (Ptrend < 0.0001) The OR comparing BMI ≥ 40 vs. < 25 kg/m2 was 3.92 (95% CI 2.91, 5.27). Abdominal obesity among those with BMI < 30 kg/m2 was not appreciably associated with EC risk (OR 1.21, 95% CI 0.74, 1.99). Associations of reproductive history with EC were similar to those observed in studies of White women. Long-term use of estrogen-only menopausal hormones was associated with an increased risk of EC (≥ 5 years vs. never use: OR 2.08, 95% CI: 1.06, 4.06). CONCLUSIONS: Our results suggest that the associations of established risk factors with EC are similar between Black and White women. Other explanations, such as differences in the prevalence of known risk factors or previously unidentified risk factors likely underlie the recent increases in EC incidence among Black women.

A region-based gene association study combined with a leave-one-out sensitivity analysis identifies SMG1 as a pancreatic cancer susceptibility gene.

Pancreatic adenocarcinoma (PC) is a lethal malignancy that is familial or associated with genetic syndromes in 10% of cases. Gene-based surveillance strategies for at-risk individuals may improve clinical outcomes. However, familial PC (FPC) is plagued by genetic heterogeneity and the genetic basis for the majority of FPC remains elusive, hampering the development of gene-based surveillance programs. The study was powered to identify genes with a cumulative pathogenic variant prevalence of at least 3%, which includes the most prevalent PC susceptibility gene, BRCA2. Since the majority of known PC susceptibility genes are involved in DNA repair, we focused on genes implicated in these pathways. We performed a region-based association study using the Mixed-Effects Score Test, followed by leave-one-out characterization of PC-associated gene regions and variants to identify the genes and variants driving risk associations. We evaluated 398 cases from two case series and 987 controls without a personal history of cancer. The first case series consisted of 109 patients with either FPC (n = 101) or PC at ≤50 years of age (n = 8). The second case series was composed of 289 unselected PC cases. We validated this discovery strategy by identifying known pathogenic BRCA2 variants, and also identified SMG1, encoding a serine/threonine protein kinase, to be significantly associated with PC following correction for multiple testing (p = 3.22x10-7). The SMG1 association was validated in a second independent series of 532 FPC cases and 753 controls (p<0.0062, OR = 1.88, 95%CI 1.17-3.03). We showed segregation of the c.4249A>G SMG1 variant in 3 affected relatives in a FPC kindred, and we found c.103G>A to be a recurrent SMG1 variant associating with PC in both the discovery and validation series. These results suggest that SMG1 is a novel PC susceptibility gene, and we identified specific SMG1 gene variants associated with PC risk.

The p.Ser64Leu and p.Pro104Leu missense variants of PALB2 identified in familial pancreatic cancer patients compromise the DNA damage response.

PALB2 has been identified as a breast and pancreatic cancer susceptibility gene. Utilizing a targeted sequencing approach, we discovered two novel germline missense PALB2 variants c.191C>T and c.311C>T, encoding p.Ser64Leu and p.Pro104Leu, respectively, in individuals in a pancreatic cancer registry. No missense PALB2 variants from familial pancreatic cancer patients, and few PALB2 variants overall, have been functionally characterized. Given the known role of PALB2, we tested the impact of p.Ser64Leu and p.Pro104Leu variants on DNA damage responses. Neither p.Ser64Leu nor p.Pro104Leu have clear effects on interactions with BRCA1 and KEAP1, which are mediated by adjacent motifs in PALB2. However, both variants are associated with defective recruitment of PALB2, and the RAD51 recombinase downstream, to DNA damage foci. Furthermore, p.Ser64Leu and p.Pro104Leu both largely compromise DNA double-strand break-initiated homologous recombination, and confer increased cellular sensitivity to ionizing radiation (IR) and the poly (ADP-ribose) polymerase (PARP) inhibitor Olaparib. Taken together, our results represent the first demonstration of functionally deleterious PALB2 missense variants associated with familial pancreatic cancer and of deleterious variants in the N-terminus outside of the coiled-coil domain. Furthermore, our results suggest the possibility of personalized treatments, using IR or PARP inhibitor, of pancreatic and other cancers that carry a deleterious PALB2 variant.

Pregnancy outcomes and risk of endometrial cancer: A pooled analysis of individual participant data in the Epidemiology of Endometrial Cancer Consortium.

A full-term pregnancy is associated with reduced endometrial cancer risk; however, whether the effect of additional pregnancies is independent of age at last pregnancy is unknown. The associations between other pregnancy-related factors and endometrial cancer risk are less clear. We pooled individual participant data from 11 cohort and 19 case-control studies participating in the Epidemiology of Endometrial Cancer Consortium (E2C2) including 16 986 women with endometrial cancer and 39 538 control women. We used one- and two-stage meta-analytic approaches to estimate pooled odds ratios (ORs) for the association between exposures and endometrial cancer risk. Ever having a full-term pregnancy was associated with a 41% reduction in risk of endometrial cancer compared to never having a full-term pregnancy (OR = 0.59, 95% confidence interval [CI] 0.56-0.63). The risk reduction appeared the greatest for the first full-term pregnancy (OR = 0.78, 95% CI 0.72-0.84), with a further ~15% reduction per pregnancy up to eight pregnancies (OR = 0.20, 95% CI 0.14-0.28) that was independent of age at last full-term pregnancy. Incomplete pregnancy was also associated with decreased endometrial cancer risk (7%-9% reduction per pregnancy). Twin births appeared to have the same effect as singleton pregnancies. Our pooled analysis shows that, while the magnitude of the risk reduction is greater for a full-term pregnancy than an incomplete pregnancy, each additional pregnancy is associated with further reduction in endometrial cancer risk, independent of age at last full-term pregnancy. These results suggest that the very high progesterone level in the last trimester of pregnancy is not the sole explanation for the protective effect of pregnancy.

A pooled genome-wide association study identifies pancreatic cancer susceptibility loci on chromosome 19p12 and 19p13.3 in the full-Jewish population.

Jews are estimated to be at increased risk of pancreatic cancer compared to non-Jews, but their observed 50-80% excess risk is not explained by known non-genetic or genetic risk factors. We conducted a GWAS in a case-control sample of American Jews, largely Ashkenazi, including 406 pancreatic cancer patients and 2332 controls, identified in the dbGaP, PanScan I/II, PanC4 and GERA data sets. We then examined resulting SNPs with P < 10-7 in an expanded sample set, of 539 full- or part-Jewish pancreatic cancer patients and 4117 full- or part-Jewish controls from the same data sets. Jewish ancestries were genetically determined using seeded FastPCA. Among the full Jews, a novel genome-wide significant association was detected on chromosome 19p12 (rs66562280, per-allele OR = 1.55, 95% CI = 1.33-1.81, P = 10-7.6). A suggestive relatively independent association was detected on chromosome 19p13.3 (rs2656937, OR = 1.53, 95% CI = 1.31-1.78, P = 10-7.0). Similar associations were seen for these SNPs among the full and part Jews combined. This is the first GWAS conducted for pancreatic cancer in the increased-risk Jewish population. The SNPs rs66562280 and rs2656937 are located in introns of ZNF100-like and ARRDC5, respectively, and are known to alter regulatory motifs of genes that play integral roles in pancreatic carcinogenesis.