Polypyrimdine tract binding protein (PTB) is a regulator of alternative splicing, mRNA 3' end formation, mRNA stability and localization, and IRES-mediated translation. Transient overexpression of PTB can influence alternative splicing, sometimes resulting in nonphysiological splicing patterns. Here, we show that alternative skipping of PTB exon 11 leads to an mRNA that is removed by NMD and that this pathway consumes at least 20% of the PTB mRNA in HeLa cells. We also show that exon 11 skipping is itself promoted by PTB in a negative feedback loop. This autoregulation may serve both to prevent disruptively high levels of PTB expression and to restore nuclear levels when PTB is mobilized to the cytoplasm. Our findings suggest that alternative splicing can act not only to generate protein isoform diversity but also to quantitatively control gene expression and complement recent bioinformatic analyses, indicating a high prevalence of human alternative splicing leading to NMD.

Transcription factors are key mediators of human complex disease processes. Identifying the target genes of transcription factors will increase our understanding of the biological network leading to disease risk. The prediction of transcription factor binding sites (TFBSs) is one method to identify these target genes; however, current prediction methods need improvement. We chose the transcription factor upstream stimulatory factor 1 ( USF1 ) to evaluate the performance of our novel TFBS prediction method because of its known genetic association with coronary artery disease (CAD) and the recent availability of USF1 chromatin immunoprecipitation microarray (ChIP-chip) results. The specific goals of our study were to develop a novel and accurate genome-scale method for predicting USF1 binding sites and associated target genes to aid in the study of CAD. Previously published USF1 ChIP-chip data for 1 per cent of the genome were used to develop and evaluate several kernel logistic regression prediction models. A combination of genomic features (phylogenetic conservation, regulatory potential, presence of a CpG island and DNaseI hypersensitivity), as well as position weight matrix (PWM) scores, were used as variables for these models. Our most accurate predictor achieved an area under the receiver operator characteristic curve of 0.827 during cross-validation experiments, significantly outperforming standard PWM-based prediction methods. When applied to the whole human genome, we predicted 24,010 USF1 binding sites within 5 kilobases upstream of the transcription start site of 9,721 genes. These predictions included 16 of 20 genes with strong evidence of USF1 regulation. Finally, in the spirit of genomic convergence, we integrated independent experimental CAD data with these USF1 binding site prediction results to develop a prioritised set of candidate genes for future CAD studies. We have shown that our novel prediction method, which employs genomic features related to the presence of regulatory elements, enables more accurate and efficient prediction of USF1 binding sites. This method can be extended to other transcription factors identified in human disease studies to help further our understanding of the biology of complex disease.

Fuchs corneal dystrophy (FCD) is a degenerative genetic disorder of the corneal endothelium that represents one of the most common causes of corneal transplantation in the United States. Despite its high prevalence (4% over the age of 40), the underlying genetic basis of FCD is largely unknown. Here we report missense mutations in TCF8, a transcription factor whose haploinsufficiency causes posterior polymorphous corneal dystrophy (PPCD), in a cohort of late-onset FCD patients. In contrast to PPCD-causing mutations, all of which are null, FCD-associated mutations encode rare missense changes suggested to cause loss of function by an in vivo complementation assay. Importantly, segregation of a recurring p.Q840P mutation in a large, multigenerational FCD pedigree showed this allele to be sufficient but not necessary for pathogenesis. Execution of a genome-wide scan conditioned for the presence of the 840P allele identified an additional late-onset FCD locus on chromosome 9p, whereas haplotype analysis indicated that the presence of the TCF8 allele and the disease haplotype on 9p leads to a severe FCD manifestation with poor prognosis. Our data suggest that PPCD and FCD are allelic variants of the same disease continuum and that genetic interaction between genes that cause corneal dystrophies can modulate the expressivity of the phenotype.

The RNA modification N6-methyladenosine (m6A) post-transcriptionally regulates RNA function. The cellular machinery that controls m6A includes methyltransferases and demethylases that add or remove this modification, as well as m6A-binding YTHDF proteins that promote the translation or degradation of m6A-modified mRNA. We demonstrate that m6A modulates infection by hepatitis C virus (HCV). Depletion of m6A methyltransferases or an m6A demethylase, respectively, increases or decreases infectious HCV particle production. During HCV infection, YTHDF proteins relocalize to lipid droplets, sites of viral assembly, and their depletion increases infectious viral particles. We further mapped m6A sites across the HCV genome and determined that inactivating m6A in one viral genomic region increases viral titer without affecting RNA replication. Additional mapping of m6A on the RNA genomes of other Flaviviridae, including dengue, Zika, yellow fever, and West Nile virus, identifies conserved regions modified by m6A. Altogether, this work identifies m6A as a conserved regulatory mark across Flaviviridae genomes.

Given the importance of cardiovascular disease (CVD) to public health and the demonstrated heritability of both disease status and its related risk factors, identifying the genetic variation underlying these susceptibilities is a critical step in understanding the pathogenesis of CVD and informing prevention and treatment strategies. Although one can look for genetic variation underlying susceptibility to CVD per se, it can be difficult to define the disease phenotype for such a qualitative analysis and CVD itself represents a convergence of diverse etiologic pathways. Alternatively, one can study the genetics of intermediate traits that are known risk factors for CVD, which can be measured quantitatively. Using the latter strategy, we have measured 21 cardiovascular-related biomarkers in an extended multigenerational pedigree, the CARRIAGE family (Carolinas Region Interaction of Aging, Genes, and Environment). These biomarkers belong to inflammatory and immune, connective tissue, lipid, and hemostasis pathways. Of these, 18 met our quality control standards. Using the pedigree and biomarker data, we have estimated the broad sense heritability (H2) of each biomarker (ranging from 0.09-0.56). A genome-wide panel of 6,015 SNPs was used subsequently to map these biomarkers as quantitative traits. Four showed noteworthy evidence for linkage in multipoint analysis (LOD score ≥ 2.6): paraoxonase (chromosome 8p11, 21), the chemokine RANTES (22q13.33), matrix metalloproteinase 3 (MMP3, 17p13.3), and granulocyte colony stimulating factor (GCSF, 8q22.1). Identifying the causal variation underlying each linkage score will help to unravel the genetic architecture of these quantitative traits and, by extension, the genetic architecture of cardiovascular risk.

Dengue virus is the most prevalent cause of arthropod-borne infection worldwide. Due to the limited coding capacity of the viral genome and the complexity of the viral life cycle, host cell proteins play essential roles throughout the course of viral infection. Host RNA-binding proteins mediate various aspects of virus replication through their physical interactions with viral RNA. Here we describe a technique designed to identify such interactions in the context of infected cells using UV cross-linking followed by antisense-mediated affinity purification and mass spectrometry. Using this approach, we identified interactions, several of them novel, between host proteins and dengue viral RNA in infected Huh7 cells. Most of these interactions were subsequently validated using RNA immunoprecipitation. Using small interfering RNA (siRNA)-mediated gene silencing, we showed that more than half of these host proteins are likely involved in regulating virus replication, demonstrating the utility of this method in identifying biologically relevant interactions that may not be identified using traditional in vitro approaches.

Survival bias may unduly impact genetic association with complex diseases; gene-specific survival effects may further complicate such investigations. Coronary artery disease (CAD) is a complex phenotype for which little is understood about gene-specific survival effects; yet, such information can offer insight into refining genetic associations, improving replications, and can provide candidate genes for both mortality risk and improved survivorship in CAD. Building on our previous work, the purpose of this current study was to: evaluate LSAMP SNP-specific hazards for all-cause mortality post-catheterization in a larger cohort of our CAD cases; and, perform additional replication in an independent dataset. We examined two LSAMP SNPs-rs1462845 and rs6788787-using CAD case-only Cox proportional hazards regression for additive genetic effects, censored on time-to-all-cause mortality or last follow-up among Caucasian subjects from the Catheterization Genetics Study (CATHGEN; n = 2,224) and the Intermountain Heart Collaborative Study (IMHC; n = 3,008). Only after controlling for age, sex, body mass index, histories of smoking, type 2 diabetes, hyperlipidemia and hypertension (HR = 1.11, 95%CI = 1.01-1.22, p = 0.032), rs1462845 conferred significantly increased hazards of all-cause mortality among CAD cases. Even after controlling for multiple covariates, but in only the primary cohort, rs6788787 conferred significantly improved survival (HR = 0.80, 95% CI = 0.69-0.92, p = 0.002). Post-hoc analyses further stratifying by sex and disease severity revealed replicated effects for rs1462845: even after adjusting for aforementioned covariates and coronary interventional procedures, males with severe burden of CAD had significantly amplified hazards of death with the minor variant of rs1462845 in both cohorts (HR = 1.29, 95% CI = 1.08-1.55, p = 0.00456; replication HR = 1.25, 95% CI = 1.05-1.49, p = 0.013). Kaplan-Meier curves revealed unique cohort-specific genotype effects on survival. Additional analyses demonstrated that the homozygous risk genotype ('A/A') fully explained the increased hazard in both cohorts. None of the post-hoc analyses in control subjects were significant for any model. This suggests that genetic effects of rs1462845 on survival are unique to CAD presence. This represents formal, replicated evidence of genetic contribution of rs1462845 to increased risk for all-cause mortality; the contribution is unique to CAD case status and specific to males with severe burden of CAD.

Coronary artery disease (CAD), and one of its intermediate risk factors, dyslipidemia, possess a demonstrable genetic component, although the genetic architecture is incompletely defined. We previously reported a linkage peak on chromosome 5q31-33 for early-onset CAD where the strength of evidence for linkage was increased in families with higher mean low density lipoprotein-cholesterol (LDL-C). Therefore, we sought to fine-map the peak using association mapping of LDL-C as an intermediate disease-related trait to further define the etiology of this linkage peak. The study populations consisted of 1908 individuals from the CATHGEN biorepository of patients undergoing cardiac catheterization; 254 families (N = 827 individuals) from the GENECARD familial study of early-onset CAD; and 162 aorta samples harvested from deceased donors. Linkage disequilibrium-tagged SNPs were selected with an average of one SNP per 20 kb for 126.6-160.2 MB (region of highest linkage) and less dense spacing (one SNP per 50 kb) for the flanking regions (117.7-126.6 and 160.2-167.5 MB) and genotyped on all samples using a custom Illumina array. Association analysis of each SNP with LDL-C was performed using multivariable linear regression (CATHGEN) and the quantitative trait transmission disequilibrium test (QTDT; GENECARD). SNPs associated with the intermediate quantitative trait, LDL-C, were then assessed for association with CAD (i.e., a qualitative phenotype) using linkage and association in the presence of linkage (APL; GENECARD) and logistic regression (CATHGEN and aortas).

Amplification of the q21-q23 region on chromosome 1 is frequently found in sarcomas and a variety of other solid tumours. Previous analyses of sarcomas have indicated the presence of at least two separate amplicons within this region, one located in 1q21 and one located near the apolipoprotein A-II (APOA2) gene in 1q23. In this study we have mapped and characterized the amplicon in 1q23 in more detail.

Multiple sclerosis (MS) is an autoimmune disorder where T cells attack neurons in the central nervous system (CNS) leading to demyelination and neurological deficits. A driver of increased MS risk is the soluble form of the interleukin-7 receptor alpha chain gene (sIL7R) produced by alternative splicing of IL7R exon 6. Here, we identified the RNA helicase DDX39B as a potent activator of this exon and consequently a repressor of sIL7R, and we found strong genetic association of DDX39B with MS risk. Indeed, we showed that a genetic variant in the 5' UTR of DDX39B reduces translation of DDX39B mRNAs and increases MS risk. Importantly, this DDX39B variant showed strong genetic and functional epistasis with allelic variants in IL7R exon 6. This study establishes the occurrence of biological epistasis in humans and provides mechanistic insight into the regulation of IL7R exon 6 splicing and its impact on MS risk.

Hundreds of cellular host factors are required to support dengue virus infection, but their identity and roles are incompletely characterized. Here, we identify human host dependency factors required for efficient dengue virus-2 (DENV2) infection of human cells. We focused on two, TTC35 and TMEM111, which we previously demonstrated to be required for yellow fever virus (YFV) infection and others subsequently showed were also required by other flaviviruses. These proteins are components of the human endoplasmic reticulum membrane protein complex (EMC), which has roles in ER-associated protein biogenesis and lipid metabolism. We report that DENV, YFV and Zika virus (ZIKV) infections were strikingly inhibited, while West Nile virus infection was unchanged, in cells that lack EMC subunit 4. Furthermore, targeted depletion of EMC subunits in live mosquitoes significantly reduced DENV2 propagation in vivo. Using a novel uncoating assay, which measures interactions between host RNA-binding proteins and incoming viral RNA, we show that EMC is required at or prior to virus uncoating. Importantly, we uncovered a second and important role for the EMC. The complex is required for viral protein accumulation in a cell line harboring a ZIKV replicon, indicating that EMC participates in the complex process of viral protein biogenesis.

Background DNA methylation is implicated in many chronic diseases and may contribute to mortality. Therefore, we conducted an epigenome-wide association study (EWAS) for all-cause mortality with whole-transcriptome data in a cardiovascular cohort (CATHGEN [Catheterization Genetics]). Methods and Results Cases were participants with mortality≥7 days postcatheterization whereas controls were alive with≥2 years of follow-up. The Illumina Human Methylation 450K and EPIC arrays (Illumina, San Diego, CA) were used for the discovery and validation sets, respectively. A linear model approach with empirical Bayes estimators adjusted for confounders was used to assess difference in methylation (Δβ). In the discovery set (55 cases, 49 controls), 25 629 (6.5%) probes were differently methylated (P<0.05). In the validation set (108 cases, 108 controls), 3 probes were differentially methylated with a false discovery rate-adjusted P<0.10: cg08215811 (SLC4A9; log2 fold change=-0.14); cg17845532 (MATK; fold change=-0.26); and cg17944110 (castor zinc finger 1 [CASZ1]; FC=0.26; P<0.0001; false discovery rate-adjusted P=0.046-0.080). Meta-analysis identified 6 probes (false discovery rate-adjusted P<0.05): the 3 above, cg20428720 (intergenic), cg17647904 (NCOR2), and cg23198793 (CAPN3). Messenger RNA expression of 2 MATK isoforms was lower in cases (fold change=-0.24 [P=0.007] and fold change=-0.61 [P=0.009]). The CASZ1, NCOR2, and CAPN3 transcripts did not show differential expression (P>0.05); the SLC4A9 transcript did not pass quality control. The cg17944110 probe is located within a potential regulatory element; expression of predicted targets (using GeneHancer) of the regulatory element, UBIAD1 (P=0.01) and CLSTN1 (P=0.03), were lower in cases. Conclusions We identified 6 novel methylation sites associated with all-cause mortality. Methylation in CASZ1 may serve as a regulatory element associated with mortality in cardiovascular patients. Larger studies are necessary to confirm these observations.

Zika virus (ZIKV) is an emerging flavivirus transmitted to humans by Aedes mosquitos. ZIKV can be transmitted from mother to fetus during pregnancy and can cause microcephaly and other birth defects. Effective vaccines for Zika are yet to be approved. Detection of the ZIKV is based on serological testing that often shows cross-reactivity with the Dengue virus (DENV) and other flaviviruses. We aimed to assemble a highly specific anti-Zika antibody panel to be utilized in the development of a highly specific and cost-effective ZIKV rapid quantification assay for viral load monitoring at point-of-care settings. To this end, we tested the affinity and specificity of twenty one commercially available monoclonal and polyclonal antibodies against ZIKV and DENV envelope proteins utilizing nine ZIKV and twelve DENV strains. We finalized and tested a panel of five antibodies for the specific detection and differentiation of ZIKV and DENV infected samples.

The ribosomal stalk proteins, RPLP1 and RPLP2 (RPLP1/2), which form the ancient ribosomal stalk, were discovered decades ago but their functions remain mysterious. We had previously shown that RPLP1/2 are exquisitely required for replication of dengue virus (DENV) and other mosquito-borne flaviviruses. Here, we show that RPLP1/2 function to relieve ribosome pausing within the DENV envelope coding sequence, leading to enhanced protein stability. We evaluated viral and cellular translation in RPLP1/2-depleted cells using ribosome profiling and found that ribosomes pause in the sequence coding for the N-terminus of the envelope protein, immediately downstream of sequences encoding two adjacent transmembrane domains (TMDs). We also find that RPLP1/2 depletion impacts a ribosome density for a small subset of cellular mRNAs. Importantly, the polarity of ribosomes on mRNAs encoding multiple TMDs was disproportionately affected by RPLP1/2 knockdown, implying a role for RPLP1/2 in multi-pass transmembrane protein biogenesis. These analyses of viral and host RNAs converge to implicate RPLP1/2 as functionally important for ribosomes to elongate through ORFs encoding multiple TMDs. We suggest that the effect of RPLP1/2 at TMD associated pauses is mediated by improving the efficiency of co-translational folding and subsequent protein stability.

The balance of myeloid populations and lymphoid populations must be well controlled. Here we found that osteopontin (OPN) skewed this balance during pathogenic conditions such as infection and autoimmunity. Notably, two isoforms of OPN exerted distinct effects in shifting this balance through cell-type-specific regulation of apoptosis. Intracellular OPN (iOPN) diminished the population size of myeloid progenitor cells and myeloid cells, and secreted OPN (sOPN) increase the population size of lymphoid cells. The total effect of OPN on skewing the leukocyte population balance was observed as host sensitivity to early systemic infection with Candida albicans and T cell-mediated colitis. Our study suggests previously unknown detrimental roles for two OPN isoforms in causing the imbalance of leukocyte populations.

Globally re-emerging dengue viruses are transmitted from human-to-human by Aedes mosquitoes. While viral determinants of human pathogenicity have been defined, there is a lack of knowledge of how dengue viruses influence mosquito transmission. Identification of viral determinants of transmission can help identify isolates with high epidemiological potential. Additionally, mechanistic understanding of transmission will lead to better understanding of how dengue viruses harness evolution to cycle between the two hosts. Here, we identified viral determinants of transmission and characterized mechanisms that enhance production of infectious saliva by inhibiting immunity specifically in salivary glands. Combining oral infection of Aedes aegypti mosquitoes and reverse genetics, we identified two 3' UTR substitutions in epidemic isolates that increased subgenomic flaviviral RNA (sfRNA) quantity, infectious particles in salivary glands and infection rate of saliva, which represents a measure of transmission. We also demonstrated that various 3'UTR modifications similarly affect sfRNA quantity in both whole mosquitoes and human cells, suggesting a shared determinism of sfRNA quantity. Furthermore, higher relative quantity of sfRNA in salivary glands compared to midgut and carcass pointed to sfRNA function in salivary glands. We showed that the Toll innate immune response was preferentially inhibited in salivary glands by viruses with the 3'UTR substitutions associated to high epidemiological fitness and high sfRNA quantity, pointing to a mechanism for higher saliva infection rate. By determining that sfRNA is an immune suppressor in a tissue relevant to mosquito transmission, we propose that 3'UTR/sfRNA sequence evolution shapes dengue epidemiology not only by influencing human pathogenicity but also by increasing mosquito transmission, thereby revealing a viral determinant of epidemiological fitness that is shared between the two hosts.

Lymphoblastoid cell lines (LCLs) are generated by transforming primary B cells with Epstein-Barr virus (EBV) and are used extensively as model systems in viral oncology, immunology, and human genetics research. In this study, we characterized single-cell transcriptomic profiles of five LCLs and present a simple discrete-time simulation to explore the influence of stochasticity on LCL clonal evolution. Single-cell RNA sequencing (scRNA-seq) revealed substantial phenotypic heterogeneity within and across LCLs with respect to immunoglobulin isotype; virus-modulated host pathways involved in survival, activation, and differentiation; viral replication state; and oxidative stress. This heterogeneity is likely attributable to intrinsic variance in primary B cells and host-pathogen dynamics. Stochastic simulations demonstrate that initial primary cell heterogeneity, random sampling, time in culture, and even mild differences in phenotype-specific fitness can contribute substantially to dynamic diversity in populations of nominally clonal cells.

Ebola virus (EBOV) genome and mRNAs contain long, structured regions that could hijack host RNA-binding proteins to facilitate infection. We performed RNA affinity chromatography coupled with mass spectrometry to identify host proteins that bind to EBOV RNAs and identified four high-confidence proviral host factors, including Staufen1 (STAU1), which specifically binds both 3' and 5' extracistronic regions of the EBOV genome. We confirmed that EBOV infection rate and production of infectious particles were significantly reduced in STAU1-depleted cells. STAU1 was recruited to sites of EBOV RNA synthesis upon infection and enhanced viral RNA synthesis. Furthermore, STAU1 interacts with EBOV nucleoprotein (NP), virion protein 30 (VP30), and VP35; the latter two bridge the viral polymerase to the NP-coated genome, forming the viral ribonucleoprotein (RNP) complex. Our data indicate that STAU1 plays a critical role in EBOV replication by coordinating interactions between the viral genome and RNA synthesis machinery.IMPORTANCE Ebola virus (EBOV) is a negative-strand RNA virus with significant public health importance. Currently, no therapeutics are available for Ebola, which imposes an urgent need for a better understanding of EBOV biology. Here we dissected the virus-host interplay between EBOV and host RNA-binding proteins. We identified novel EBOV host factors, including Staufen1, which interacts with multiple viral factors and is required for efficient viral RNA synthesis.

Paramyxoviruses and pneumoviruses have similar life cycles and share the respiratory tract as a point of entry. In comparative genome-scale siRNA screens with wild-type-derived measles, mumps, and respiratory syncytial viruses in A549 cells, a human lung adenocarcinoma cell line, we identified vesicular transport, RNA processing pathways, and translation as the top pathways required by all three viruses. As the top hit in the translation pathway, ABCE1, a member of the ATP-binding cassette transporters, was chosen for further study. We found that ABCE1 supports replication of all three viruses, confirming its importance for viruses of both families. More detailed characterization revealed that ABCE1 is specifically required for efficient viral but not general cellular protein synthesis, indicating that paramyxoviral and pneumoviral mRNAs exploit specific translation mechanisms. In addition to providing a novel overview of cellular proteins and pathways that impact these important pathogens, this study highlights the role of ABCE1 as a host factor required for efficient paramyxovirus and pneumovirus translation.IMPORTANCE The Paramyxoviridae and Pneumoviridae families include important human and animal pathogens. To identify common host factors, we performed genome-scale siRNA screens with wild-type-derived measles, mumps, and respiratory syncytial viruses in the same cell line. A comparative bioinformatics analysis yielded different members of the coatomer complex I, translation factors ABCE1 and eIF3A, and several RNA binding proteins as cellular proteins with proviral activity for all three viruses. A more detailed characterization of ABCE1 revealed its essential role for viral protein synthesis. Taken together, these data sets provide new insight into the interactions between paramyxoviruses and pneumoviruses and host cell proteins and constitute a starting point for the development of broadly effective antivirals.

During pandemics, individuals exhibit differences in risk and clinical outcomes. Here, we developed single-cell high-throughput human in vitro susceptibility testing (scHi-HOST), a method for rapidly identifying genetic variants that confer resistance and susceptibility. We applied this method to influenza A virus (IAV), the cause of four pandemics since the start of the 20th century. scHi-HOST leverages single-cell RNA sequencing (scRNA-seq) to simultaneously assign genetic identity to cells in mixed infections of cell lines of European, African, and Asian origin, reveal associated genetic variants for viral burden, and identify expression quantitative trait loci. Integration of scHi-HOST with human challenge and experimental validation demonstrated that a missense variant in endoplasmic reticulum aminopeptidase 1 (ERAP1; rs27895) increased IAV burden in cells and human volunteers. rs27895 exhibits population differentiation, likely contributing to greater permissivity of cells from African populations to IAV. scHi-HOST is a broadly applicable method and resource for decoding infectious-disease genetics.