Understanding gene regulation and function requires a genome-wide method capable of capturing both gene expression levels and isoform diversity at the single-cell level. Short-read RNAseq is limited in its ability to resolve complex isoforms because it fails to sequence full-length cDNA copies of RNA molecules. Here, we investigate whether RNAseq using the long-read single-molecule Oxford Nanopore MinION sequencer is able to identify and quantify complex isoforms without sacrificing accurate gene expression quantification. After benchmarking our approach, we analyse individual murine B1a cells using a custom multiplexing strategy. We identify thousands of unannotated transcription start and end sites, as well as hundreds of alternative splicing events in these B1a cells. We also identify hundreds of genes expressed across B1a cells that display multiple complex isoforms, including several B cell-specific surface receptors. Our results show that we can identify and quantify complex isoforms at the single cell level.

Emerging influenza viruses are a serious threat to human health because of their pandemic potential. A promising target for the development of novel anti-influenza therapeutics is the PA protein, whose endonuclease activity is essential for viral replication. Translation of viral mRNAs by the host ribosome requires mRNA capping for recognition and binding, and the necessary mRNA caps are cleaved or “snatched” from host pre-mRNAs by the PA endonuclease. The structure-based development of inhibitors that target PA endonuclease is now possible with the recent crystal structure of the PA catalytic domain. In this study, we sought to understand the molecular mechanism of inhibition by several compounds that are known or predicted to block endonuclease-dependent polymerase activity. Using an in vitro endonuclease activity assay, we show that these compounds block the enzymatic activity of the isolated PA endonuclease domain. Using X-ray crystallography, we show how these inhibitors coordinate the two-metal endonuclease active site and engage the active site residues. Two structures also reveal an induced-fit mode of inhibitor binding. The structures allow a molecular understanding of the structure-activity relationship of several known influenza inhibitors and the mechanism of drug resistance by a PA mutation. Taken together, our data reveal new strategies for structure-based design and optimization of PA endonuclease inhibitors.

Chromatin structure plays a pivotal role in facilitating proper control of gene expression. Transcription factor (TF) binding of cis-elements is often associated with accessible chromatin regions. Therefore, the ability to identify these accessible regions throughout plant genomes will advance understanding of the relationship between TF binding, chromatin status and the regulation of gene expression. Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) is a recently developed technique used to map open chromatin zones in animal genomes. However, in plants, the existence of cell walls, subcellular organelles and the lack of stable cell lines have prevented routine application of this technique. Here, we describe an assay combining ATAC-seq with fluorescence-activated nuclei sorting (FANS) to identify and map open chromatin and TF-binding sites in plant genomes. FANS-ATAC-seq compares favorably with published DNaseI sequencing (DNase-seq) results and it requires less than 50 000 nuclei for accurate identification of accessible genomic regions.Summary: Application of ATAC-seq to sorted nuclei identifies accessible regions genome-wide.

Highly pathogenic avian influenza viruses of the H5N1 subtype continue to threaten agriculture and human health. Here, we use biochemistry and x-ray crystallography to reveal how amino-acid variations in the hemagglutinin (HA) protein contribute to the pathogenicity of H5N1 influenza virus in chickens. HA proteins from highly pathogenic (HP) A/chicken/Hong Kong/YU562/2001 and moderately pathogenic (MP) A/goose/Hong Kong/437-10/1999 isolates of H5N1 were found to be expressed and cleaved in similar amounts, and both proteins had similar receptor-binding properties. However, amino-acid variations at positions 104 and 115 in the vestigial esterase sub-domain of the HA1 receptor-binding domain (RBD) were found to modulate the pH of HA activation such that the HP and MP HA proteins are activated for membrane fusion at pH 5.7 and 5.3, respectively. In general, an increase in H5N1 pathogenicity in chickens was found to correlate with an increase in the pH of HA activation for mutant and chimeric HA proteins in the observed range of pH 5.2 to 6.0. We determined a crystal structure of the MP HA protein at 2.50 Å resolution and two structures of HP HA at 2.95 and 3.10 Å resolution. Residues 104 and 115 that modulate the acid stability of the HA protein are situated at the N- and C-termini of the 110-helix in the vestigial esterase sub-domain, which interacts with the B loop of the HA2 stalk domain. Interactions between the 110-helix and the stalk domain appear to be important in regulating HA protein acid stability, which in turn modulates influenza virus replication and pathogenesis. Overall, an optimal activation pH of the HA protein is found to be necessary for high pathogenicity by H5N1 influenza virus in avian species.

Little is known about the three-dimensional organization of rubella virus, which causes a relatively mild measles-like disease in children but leads to serious congenital health problems when contracted in utero. Although rubella virus belongs to the same family as the mosquito-borne alphaviruses, in many respects it is more similar to other aerosol-transmitted human viruses such as the agents of measles and mumps. Although the use of the triple MMR (measles, mumps and rubella) live vaccine has limited its incidence in western countries, congenital rubella syndrome remains an important health problem in the developing world. Here we report the 1.8 Å resolution crystal structure of envelope glycoprotein E1, the main antigen and sole target of neutralizing antibodies against rubella virus. E1 is the main player during entry into target cells owing to its receptor-binding and membrane-fusion functions. The structure reveals the epitope and the neutralization mechanism of an important category of protecting antibodies against rubella infection. It also shows that rubella virus E1 is a class II fusion protein, which had hitherto only been structurally characterized for the arthropod-borne alphaviruses and flaviviruses. In addition, rubella virus E1 has an extensive membrane-fusion surface that includes a metal site, reminiscent of the T-cell immunoglobulin and mucin family of cellular proteins that bind phosphatidylserine lipids at the plasma membrane of cells undergoing apoptosis. Such features have not been seen in any fusion protein crystallized so far. Structural comparisons show that the class II fusion proteins from alphaviruses and flaviviruses, despite belonging to different virus families, are closer to each other than they are to rubella virus E1. This suggests that the constraints on arboviruses imposed by alternating cycles between vertebrates and arthropods resulted in more conservative evolution. By contrast, in the absence of this constraint, the strictly human rubella virus seems to have drifted considerably into a unique niche as sole member of the Rubivirus genus.