GABPα and CREB1 Binding to Double Nucleotide Polymorphisms of Their Consensus Motifs and Cooperative Binding to the Composite ETS ⇔ CRE Motif (ACCGGAAGTGACGTCA).

Previously, cooperative binding of the bZIP domain of CREB1 and the ETS domain of GABPα was observed for the composite DNA ETS ⇔ CRE motif (A 0 C 1 C 2 G 3 G 4 A 5 A 6 G 7 T 8 G 9 A 10 C 11 G 12 T 13 C 14 A 15 ). Single nucleotide polymorphisms (SNPs) at the beginning and end of the ETS motif (ACCGGAAGT) increased cooperative binding. Here, we use an Agilent microarray of 60-mers containing all double nucleotide polymorphisms (DNPs) of the ETS ⇔ CRE motif to explore GABPα and CREB1 binding to their individual motifs and their cooperative binding. For GABPα, all DNPs were bound as if each SNP acted independently. In contrast, CREB1 binding to some DNPs was stronger or weaker than expected, depending on the locations of each SNP. CREB1 binding to DNPs where both SNPs were in the same half site, T 8 G 9 A 10 or T 13 C 14 A 15 , was greater than expected, indicating that an additional SNP cannot destroy binding as much as expected, suggesting that an individual SNP is enough to abolish sequence-specific DNA binding of a single bZIP monomer. If a DNP contains SNPs in each half site, binding is weaker than expected. Similar results were observed for additional ETS and bZIP family members. Cooperative binding between GABPα and CREB1 to the ETS ⇔ CRE motif was weaker than expected except for DNPs containing A 7 and SNPs at the beginning of the ETS motif.

Metagenomic Discovery of 83 New Human Papillomavirus Types in Patients with Immunodeficiency.

Several immunodeficiencies are associated with high susceptibility to persistent and progressive human papillomavirus (HPV) infection leading to a wide range of cutaneous and mucosal lesions. However, the HPV types most commonly associated with such clinical manifestations in these patients have not been systematically defined. Here, we used virion enrichment, rolling circle amplification, and deep sequencing to identify circular DNA viruses present in skin swabs and/or wart biopsy samples from 48 patients with rare genetic immunodeficiencies, including patients with warts, hypogammaglobulinemia, infections, myelokathexis (WHIM) syndrome, or epidermodysplasia verruciformis (EV). Their profiles were compared with the profiles of swabs from 14 healthy adults and warts from 6 immunologically normal children. Individual patients were typically infected with multiple HPV types; up to 26 different types were isolated from a single patient (multiple anatomical sites, one time point). Among these, we identified the complete genomes of 83 previously unknown HPV types and 35 incomplete genomes representing possible additional new types. HPV types in the genus Gammapapillomavirus were common in WHIM patients, whereas EV patients mainly shed HPVs from the genus Betapapillomavirus. Preliminary evidence based on three WHIM patients treated with plerixafor, a leukocyte mobilizing agent, suggest that longer-term therapy may correlate with decreased HPV diversity and increased predominance of HPV types associated with childhood skin warts.IMPORTANCE Although some members of the viral family Papillomaviridae cause benign skin warts (papillomas), many human papillomavirus (HPV) infections are not associated with visible symptoms. For example, most healthy adults chronically shed Gammapapillomavirus (Gamma) virions from apparently healthy skin surfaces. To further explore the diversity of papillomaviruses, we performed viromic surveys on immunodeficient individuals suffering from florid skin warts. Our results nearly double the number of known Gamma HPV types and suggest that WHIM syndrome patients are uniquely susceptible to Gamma HPV-associated skin warts. Preliminary results suggest that treatment with the drug plerixafor may promote resolution of the unusual Gamma HPV skin warts observed in WHIM patients.

Characterization of BK Polyomaviruses from Kidney Transplant Recipients Suggests a Role for APOBEC3 in Driving In-Host Virus Evolution.

BK polyomavirus (BKV) frequently causes nephropathy (BKVN) in kidney transplant recipients (KTRs). BKV has also been implicated in the etiology of bladder and kidney cancers. We characterized BKV variants from two KTRs who developed BKVN followed by renal carcinoma. Both patients showed a swarm of BKV sequence variants encoding non-silent mutations in surface loops of the viral major capsid protein. The temporal appearance and disappearance of these mutations highlights the intra-patient evolution of BKV. Some of the observed mutations conferred resistance to antibody-mediated neutralization. The mutations also modified the spectrum of receptor glycans engaged by BKV during host cell entry. Intriguingly, all observed mutations were consistent with DNA damage caused by antiviral APOBEC3 cytosine deaminases. Moreover, APOBEC3 expression was evident upon immunohistochemical analysis of renal biopsies from KTRs. These results provide a snapshot of in-host BKV evolution and suggest that APOBEC3 may drive BKV mutagenesis in vivo.

An Efficient Single-Cell RNA-Seq Approach to Identify Neoantigen-Specific T Cell Receptors.

The adoptive transfer of neoantigen-reactive tumor-infiltrating lymphocytes (TILs) can result in tumor regression in patients with metastatic cancer. To improve the efficacy of adoptive T cell therapy targeting these tumor-specific mutations, we have proposed a new therapeutic strategy, which involves the genetic modification of autologous T cells with neoantigen-specific T cell receptors (TCRs) and the transfer of these modified T cells back to cancer patients. However, the current techniques to isolate neoantigen-specific TCRs are labor intensive, time consuming, and technically challenging, not suitable for clinical applications. To facilitate this process, a new approach was developed, which included the co-culture of TILs with tandem minigene (TMG)-transfected or peptide-pulsed autologous antigen-presenting cells (APCs) and the single-cell RNA sequencing (RNA-seq) analysis of T cells to identify paired TCR sequences associated with cells expressing high levels of interferon-γ (IFN-γ) and interleukin-2 (IL-2). Following this new approach, multiple TCRs were identified, synthesized, cloned into a retroviral vector, and then transduced into donor T cells. These transduced T cells were shown to specifically recognize the neoantigens presented by autologous APCs. In conclusion, this approach provides an efficient procedure to isolate neoantigen-specific TCRs for clinical applications, as well as for basic and translational research.

Methylated Cytosines Mutate to Transcription Factor Binding Sites that Drive Tetrapod Evolution.

In mammals, the cytosine in CG dinucleotides is typically methylated producing 5-methylcytosine (5mC), a chemically less stable form of cytosine that can spontaneously deaminate to thymidine resulting in a T•G mismatched base pair. Unlike other eukaryotes that efficiently repair this mismatched base pair back to C•G, in mammals, 5mCG deamination is mutagenic, sometimes producing TG dinucleotides, explaining the depletion of CG dinucleotides in mammalian genomes. It was suggested that new TG dinucleotides generate genetic diversity that may be critical for evolutionary change. We tested this conjecture by examining the DNA sequence properties of regulatory sequences identified by DNase I hypersensitive sites (DHSs) in human and mouse genomes. We hypothesized that the new TG dinucleotides generate transcription factor binding sites (TFBS) that become tissue-specific DHSs (TS-DHSs). We find that 8-mers containing the CG dinucleotide are enriched in DHSs in both species. However, 8-mers containing a TG and no CG dinucleotide are preferentially enriched in TS-DHSs when compared with 8-mers with neither a TG nor a CG dinucleotide. The most enriched 8-mer with a TG and no CG dinucleotide in tissue-specific regulatory regions in both genomes is the AP-1 motif ( TG: A(C)/GT CA: N), and we find evidence that TG dinucleotides in the AP-1 motif arose from CG dinucleotides. Additional TS-DHS-enriched TFBS containing the TG/CA dinucleotide are the E-Box motif (G CA: GC TG: C), the NF-1 motif (GG CATG: CC), and the GR (glucocorticoid receptor) motif (G-A CATG: T-C). Our results support the suggestion that cytosine methylation is mutagenic in tetrapods producing TG dinucleotides that create TFBS that drive evolution.

Hamburger polyomaviruses.

Epidemiological studies have suggested that consumption of beef may correlate with an increased risk of colorectal cancer. One hypothesis to explain this proposed link might be the presence of a carcinogenic infectious agent capable of withstanding cooking. Polyomaviruses are a ubiquitous family of thermostable non-enveloped DNA viruses that are known to be carcinogenic. Using virion enrichment, rolling circle amplification (RCA) and next-generation sequencing, we searched for polyomaviruses in meat samples purchased from several supermarkets. Ground beef samples were found to contain three polyomavirus species. One species, bovine polyomavirus 1 (BoPyV1), was originally discovered as a contaminant in laboratory FCS. A previously unknown species, BoPyV2, occupies the same clade as human Merkel cell polyomavirus and raccoon polyomavirus, both of which are carcinogenic in their native hosts. A third species, BoPyV3, is related to human polyomaviruses 6 and 7. Examples of additional DNA virus families, including herpesviruses, adenoviruses, circoviruses and gyroviruses were also detected either in ground beef samples or in comparison samples of ground pork and ground chicken. The results suggest that the virion enrichment/RCA approach is suitable for random detection of essentially any DNA virus with a detergent-stable capsid. It will be important for future studies to address the possibility that animal viruses commonly found in food might be associated with disease.

Genome Sequence of a Fish-Associated Polyomavirus, Black Sea Bass (Centropristis striata) Polyomavirus 1.

All known polyomaviruses are associated with mammals or birds. Using virion enrichment, random-primed rolling circle amplification, and deep sequencing, we identified a polyomavirus associated with black sea bass (Centropristis striata). The virus has a variety of novel genetic features, suggesting a long evolutionary separation from polyomaviruses of terrestrial animals.

A divergent variant of the eleventh human polyomavirus species, saint louis polyomavirus.

Saint Louis polyomavirus (STLPyV) was recently discovered in human feces. Using random-primed rolling circle amplification combined with deep sequencing, we have found a divergent variant of STLPyV in a sanitized human skin wart specimen. The result strongly suggests that STLPyV directly infects humans and is not simply a dietary contaminant.

Nucleosome-free region dominates histone acetylation in targeting SWR1 to promoters for H2A.Z replacement.

The histone variant H2A.Z is a genome-wide signature of nucleosomes proximal to eukaryotic regulatory DNA. Whereas the multisubunit chromatin remodeler SWR1 is known to catalyze ATP-dependent deposition of H2A.Z, the mechanism of SWR1 recruitment to S. cerevisiae promoters has been unclear. A sensitive assay for competitive binding of dinucleosome substrates revealed that SWR1 preferentially binds long nucleosome-free DNA and the adjoining nucleosome core particle, allowing discrimination of gene promoters over gene bodies. Analysis of mutants indicates that the conserved Swc2/YL1 subunit and the adenosine triphosphatase domain of Swr1 are mainly responsible for binding to substrate. SWR1 binding is enhanced on nucleosomes acetylated by the NuA4 histone acetyltransferase, but recognition of nucleosome-free and nucleosomal DNA is dominant over interaction with acetylated histones. Such hierarchical cooperation between DNA and histone signals expands the dynamic range of genetic switches, unifying classical gene regulation by DNA-binding factors with ATP-dependent nucleosome remodeling and posttranslational histone modifications.

Contribution of nucleosome binding preferences and co-occurring DNA sequences to transcription factor binding.

BACKGROUND: Chromatin plays a critical role in regulating transcription factors (TFs) binding to their canonical transcription factor binding sites (TFBS). Recent studies in vertebrates show that many TFs preferentially bind to genomic regions that are well bound by nucleosomes in vitro. Co-occurring secondary motifs sometimes correlated with functional TFBS. RESULTS: We used a logistic regression to evaluate how well the propensity for nucleosome binding and co-occurrence of a secondary motif identify which canonical motifs are bound in vivo. We used ChIP-seq data for three transcription factors binding to their canonical motifs: c-Jun binding the AP-1 motif (TGA(C)/(G)TCA), GR (glucocorticoid receptor) binding the GR motif (G-ACA---(T)/(C)GT-C), and Hoxa2 (homeobox a2) binding the Pbx (Pre-B-cell leukemia homeobox) motif (TGATTGAT). For all canonical TFBS in the mouse genome, we calculated intrinsic nucleosome occupancy scores (INOS) for its surrounding 150-bps DNA and examined the relationship with in vivo TF binding. In mouse mammary 3134 cells, c-Jun and GR proteins preferentially bound regions calculated to be well-bound by nucleosomes in vitro with the canonical AP-1 and GR motifs themselves contributing to the high INOS. Functional GR motifs are enriched for AP-1 motifs if they are within a nucleosome-sized 150-bps region. GR and Hoxa2 also bind motifs with low INOS, perhaps indicating a different mechanism of action. CONCLUSION: Our analysis quantified the contribution of INOS and co-occurring sequence to the identification of functional canonical motifs in the genome. This analysis revealed an inherent competition between some TFs and nucleosomes for binding canonical TFBS. GR and c-Jun cooperate if they are within 150-bps. Binding of Hoxa2 and a fraction of GR to motifs with low INOS values suggesting they are not in competition with nucleosomes and may function using different mechanisms.

Overlapping ETS and CRE Motifs ((G/C)CGGAAGTGACGTCA) preferentially bound by GABPα and CREB proteins.

Previously, we identified 8-bps long DNA sequences (8-mers) that localize in human proximal promoters and grouped them into known transcription factor binding sites (TFBS). We now examine split 8-mers consisting of two 4-mers separated by 1-bp to 30-bps (X(4)-N(1-30)-X(4)) to identify pairs of TFBS that localize in proximal promoters at a precise distance. These include two overlapping TFBS: the ETS⇔ETS motif ((C/G)CCGGAAGCGGAA) and the ETS⇔CRE motif ((C/G)CGGAAGTGACGTCAC). The nucleotides in bold are part of both TFBS. Molecular modeling shows that the ETS⇔CRE motif can be bound simultaneously by both the ETS and the B-ZIP domains without protein-protein clashes. The electrophoretic mobility shift assay (EMSA) shows that the ETS protein GABPα and the B-ZIP protein CREB preferentially bind to the ETS⇔CRE motif only when the two TFBS overlap precisely. In contrast, the ETS domain of ETV5 and CREB interfere with each other for binding the ETS⇔CRE. The 11-mer (CGGAAGTGACG), the conserved part of the ETS⇔CRE motif, occurs 226 times in the human genome and 83% are in known regulatory regions. In vivo GABPα and CREB ChIP-seq peaks identified the ETS⇔CRE as the most enriched motif occurring in promoters of genes involved in mRNA processing, cellular catabolic processes, and stress response, suggesting that a specific class of genes is regulated by this composite motif.

Stepwise histone replacement by SWR1 requires dual activation with histone H2A.Z and canonical nucleosome.

Histone variant H2A.Z-containing nucleosomes are incorporated at most eukaryotic promoters. This incorporation is mediated by the conserved SWR1 complex, which replaces histone H2A in canonical nucleosomes with H2A.Z in an ATP-dependent manner. Here, we show that promoter-proximal nucleosomes are highly heterogeneous for H2A.Z in Saccharomyces cerevisiae, with substantial representation of nucleosomes containing one, two, or zero H2A.Z molecules. SWR1-catalyzed H2A.Z replacement in vitro occurs in a stepwise and unidirectional fashion, one H2A.Z-H2B dimer at a time, producing heterotypic nucleosomes as intermediates and homotypic H2A.Z nucleosomes as end products. The ATPase activity of SWR1 is specifically stimulated by H2A-containing nucleosomes without ensuing histone H2A eviction. Remarkably, further addition of free H2A.Z-H2B dimer leads to hyperstimulation of ATPase activity, eviction of nucleosomal H2A-H2B, and deposition of H2A.Z-H2B. These results suggest that the combination of H2A-containing nucleosome and free H2A.Z-H2B dimer acting as both effector and substrate for SWR1 governs the specificity and outcome of the replacement reaction.

The conserved tyrosine residues 401 and 1044 in ATP sites of human P-glycoprotein are critical for ATP binding and hydrolysis: evidence for a conserved subdomain, the A-loop in the ATP-binding cassette.

Each nucleotide-binding domain (NBD) of mammalian P-glycoproteins (Pgps) and human ATP-binding cassette (ABC) B subfamily members contains a tyrosine residue approximately 25 residues upstream of the Walker A domain. To assess the role of the conserved Y401 and Y1044 residues of human Pgp, we substituted these residues with F, W, C, or A either singly or together. The mutant proteins were expressed in a Vaccinia virus-based transient expression system as well as in baculovirus-infected HighFive insect cells. The Y401F, Y401W, Y1044F, Y1044W, or Y401F/Y1004F mutants transported fluorescent substrates similar to the wild-type protein. On the other hand, Y401L and Y401C exhibited partial (30-50%) function, and transport was completely abolished in Y401A, Y1044A, and Y401A/Y1044A mutant Pgps. Similarly, in Y401A, Y1044A, and Y401A/Y1044A mutants, TNP-ATP binding, vanadate-induced trapping of nucleotide, and ATP hydrolysis were completely abolished. Thus, an aromatic residue upstream of the Walker A motif in ABC transporters is critical for binding of ATP. Additionally, the crystal structures of several NBDs in the nucleotide-bound form, data mining, and alignment of 18,514 ABC domains with the consensus conserved sequence in a database of all nonredundant proteins indicate that an aromatic residue is highly conserved in approximately 85% of ABC proteins. Although the role of this aromatic residue has previously been studied in a few ABC proteins, we provide evidence for a near-universal structural and functional role for this residue and recognize its presence as a conserved subdomain approximately 25 amino acids upstream of the Walker A motif that is critical for ATP binding. We named this subdomain the "A-loop" (aromatic residue interacting with the adenine ring of ATP).

The mouse S100A15 ortholog parallels genomic organization, structure, gene expression, and protein-processing pattern of the human S100A7/A15 subfamily during epidermal maturation.

The calcium-binding proteins of the human S100A7/A15 (hS100A7/A15) subfamily are differentially expressed in normal and pathological epidermis. The hS100A7 (psoriasin) and S100A15 reside in a chromosomal cluster of highly similar paralogs. To exploit the power of mouse models for determining functions of gene products, the corresponding S100A7/A15 ortholog was cloned and examined in murine skin. The single mouse S100A15 (mS100A15) gene encodes a protein of 104 amino acids with a predicted molecular weight of 12,870 Da and two EF-hand calcium binding sites. Using gene-specific primers and specific antibodies, expression of mS100A15 in both skin and isolated keratinocytes is confined to differentiating granular and cornified epidermal cells. Immunoblotting of epidermal extracts revealed a series of high molecular weight bands that are also recognized by an antibody for transglutaminase-mediated protein crosslinks. mS100A15 expression is upregulated in cultured keratinocytes induced to differentiate by calcium or phorbol esters. Maximal induction occurs concordantly with expression of late differentiation markers. Induction is enhanced in keratinocytes overexpressing protein kinase Calpha and is dependent on activator protein-1 transcription factors. The regulation, expression pattern and crosslinking of mS100A15 are consistent with the characteristics of the human orthologs, providing a valid surrogate model to study changes in these proteins associated with cutaneous pathologies.

Identification of tandem duplicate regulatory small RNAs in Pseudomonas aeruginosa involved in iron homeostasis.

In many bacteria, iron homeostasis is controlled primarily by the ferric uptake regulator (Fur), a transcriptional repressor. However, some genes, including those involved in iron storage, are positively regulated by Fur. A Fur-repressed regulatory small RNA (sRNA), RyhB, has been identified in Escherichia coli, and it has been demonstrated that negative regulation of genes by this sRNA is responsible for the positive regulation of some genes by Fur. No RyhB sequence homologs were found in Pseudomonas aeruginosa, despite the identification of genes positively regulated by its Fur homolog. A bioinformatics approach identified two tandem sRNAs in P. aeruginosa that were candidates for functional homologs of RyhB. These sRNAs (PrrF1 and PrrF2) are >95% identical to each other, and a functional Fur box precedes each. Their expression is induced under iron limitation. Deletion of both sRNAs is required to affect the iron-dependent regulation of an array of genes, including those involved in resistance to oxidative stress, iron storage, and intermediary metabolism. As in E. coli, induction of the PrrF sRNAs leads to the rapid loss of mRNAs for sodB (superoxide dismutase), sdh (succinate dehydrogenase), and a gene encoding a bacterioferritin. Thus, the PrrF sRNAs are the functional homologs of RyhB sRNA. At least one gene, bfrB, is positively regulated by Fur and Fe(2+), even in the absence of the PrrF sRNAs. This work suggests that the role of sRNAs in bacterial iron homeostasis may be broad, and approaches similar to those described here may identify these sRNAs in other organisms.

Reactivity of human sera in a sensitive, high-throughput pseudovirus-based papillomavirus neutralization assay for HPV16 and HPV18.

Sensitive high-throughput neutralization assays, based upon pseudoviruses carrying a secreted alkaline phosphatase (SEAP) reporter gene, were developed and validated for human papillomavirus (HPV)16, HPV18, and bovine papillomavirus 1 (BPV1). SEAP pseudoviruses were produced by transient transfection of codon-modified papillomavirus structural genes into an SV40 T antigen expressing line derived from 293 cells, yielding sufficient pseudovirus from one flask for thousands of titrations. In a 96-well plate format, in this initial characterization, the assay was reproducible and appears to be as sensitive as, but more specific than, a standard papillomavirus-like particle (VLP)-based enzyme-linked immunosorbent assay (ELISA). The neutralization assay detected type-specific HPV16 or HPV18 neutralizing antibodies (titers of 160-10240) in sera of the majority of a group of women infected with the corresponding HPV type, but not in virgin women. Sera from HPV16 VLP vaccinees had high anti-HPV16 neutralizing titers (mean: 45000; range: 5120-163840), but no anti-HPV18 neutralizing activity. The SEAP pseudovirus-based neutralization assay should be a practical method for quantifying potentially protective antibody responses in HPV natural history and prophylactic vaccine studies.