Functional Segregation of Overlapping Genes in HIV.

Overlapping genes pose an evolutionary dilemma as one DNA sequence evolves under the selection pressures of multiple proteins. Here, we perform systematic statistical and mutational analyses of the overlapping HIV-1 genes tat and rev and engineer exhaustive libraries of non-overlapped viruses to perform deep mutational scanning of each gene independently. We find a "segregated" organization in which overlapped sites encode functional residues of one gene or the other, but never both. Furthermore, this organization eliminates unfit genotypes, providing a fitness advantage to the population. Our comprehensive analysis reveals the extraordinary manner in which HIV minimizes the constraint of overlapping genes and repurposes that constraint to its own advantage. Thus, overlaps are not just consequences of evolutionary constraints, but rather can provide population fitness advantages.

Development of a novel, highly sensitive assay for quantification of minimal residual B cells in autoimmune disease and comparison to traditional methods across B-cell-depleting agents.

Targeted B-cell depletion is a useful therapy for many diseases, including autoimmune disorders and certain cancers. We developed a sensitive blood B-cell depletion assay, MRB 1.1, compared its performance with the T-cell/B-cell/NK-cell (TBNK) assay, and assessed B-cell depletion with different therapies. The empirically defined lower limit of quantification (LLOQ) for CD19+ cells in the TBNK assay was 10 cells/µL, and 0.441 cells/µL for the MRB 1.1 assay. The TBNK LLOQ was used to compare differences between B-cell depletion in similar lupus nephritis patient populations who received rituximab (LUNAR), ocrelizumab (BELONG), or obinutuzumab (NOBILITY). After 4 weeks, 10% of patients treated with rituximab retained detectable B cells vs 1.8% with ocrelizumab and 1.7% for obinutuzumab; at 24 weeks 93% of patients who received obinutuzumab remained below LLOQ vs 63% for rituximab. More-sensitive measurements of B cells may reveal differences in potency among anti-CD20 agents, which may associate with clinical outcomes.

Effect of Ocrelizumab on B- and T-Cell Receptor Repertoire Diversity in Patients With Relapsing Multiple Sclerosis From the Randomized Phase III OPERA Trial.

BACKGROUND AND OBJECTIVES: The B cell-depleting anti-CD20 antibody ocrelizumab (OCR) effectively reduces MS disease activity and slows disability progression. Given the role of B cells as antigen-presenting cells, the primary goal of this study was to evaluate the effect of OCR on the T-cell receptor repertoire diversity. METHODS: To examine whether OCR substantially alters the molecular diversity of the T-cell receptor repertoire, deep immune repertoire sequencing (RepSeq) of CD4+ and CD8+ T-cell receptor ß-chain variable regions was performed on longitudinal blood samples. The IgM and IgG heavy chain variable region repertoire was also analyzed to characterize the residual B-cell repertoire under OCR treatment. RESULTS: Peripheral blood samples for RepSeq were obtained from 8 patients with relapsing MS enrolled in the OPERA I trial over a period of up to 39 months. Four patients each were treated with OCR or interferon ß1-a during the double-blind period of OPERA I. All patients received OCR during the open-label extension. The diversity of the CD4+/CD8+ T-cell repertoires remained unaffected in OCR-treated patients. The expected OCR-associated B-cell depletion was mirrored by reduced B-cell receptor diversity in peripheral blood and a shift in immunoglobulin gene usage. Despite deep B-cell depletion, longitudinal persistence of clonally related B-cells was observed. DISCUSSION: Our data illustrate that the diversity of CD4+/CD8+ T-cell receptor repertoires remained unaltered in OCR-treated patients with relapsing MS. Persistence of a highly diverse T-cell repertoire suggests that aspects of adaptive immunity remain intact despite extended anti-CD20 therapy. TRIAL REGISTRATION INFORMATION: This is a substudy (BE29353) of the OPERA I (WA21092; NCT01247324) trial. Date of registration, November 23, 2010; first patient enrollment, August 31, 2011.

Profiling immunoglobulin repertoires across multiple human tissues using RNA sequencing.

Profiling immunoglobulin (Ig) receptor repertoires with specialized assays can be cost-ineffective and time-consuming. Here we report ImReP, a computational method for rapid and accurate profiling of the Ig repertoire, including the complementary-determining region 3 (CDR3), using regular RNA sequencing data such as those from 8,555 samples across 53 tissues types from 544 individuals in the Genotype-Tissue Expression (GTEx v6) project. Using ImReP and GTEx v6 data, we generate a collection of 3.6 million Ig sequences, termed the atlas of immunoglobulin repertoires (TAIR), across a broad range of tissue types that often do not have reported Ig repertoires information. Moreover, the flow of Ig clonotypes and inter-tissue repertoire similarities across immune-related tissues are also evaluated. In summary, TAIR is one of the largest collections of CDR3 sequences and tissue types, and should serve as an important resource for studying immunological diseases.

Author Correction: Profiling immunoglobulin repertoires across multiple human tissues using RNA sequencing.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Ebola virus, but not Marburg virus, replicates efficiently and without required adaptation in snake cells.

Ebola virus (EBOV) disease is a viral hemorrhagic fever with a high case-fatality rate in humans. This disease is caused by four members of the filoviral genus Ebolavirus, including EBOV. The natural hosts reservoirs of ebolaviruses remain to be identified. Glycoprotein 2 of reptarenaviruses, known to infect only boa constrictors and pythons, is similar in sequence and structure to ebolaviral glycoprotein 2, suggesting that EBOV may be able to infect reptilian cells. Therefore, we serially passaged EBOV and a distantly related filovirus, Marburg virus (MARV), in boa constrictor JK cells and characterized viral infection/replication and mutational frequency by confocal imaging and sequencing. We observed that EBOV efficiently infected and replicated in JK cells, but MARV did not. In contrast to most cell lines, EBOV-infected JK cells did not result in an obvious cytopathic effect. Surprisingly, genomic characterization of serial-passaged EBOV in JK cells revealed that genomic adaptation was not required for infection. Deep sequencing coverage (>10,000×) demonstrated the existence of only a single nonsynonymous variant (EBOV glycoprotein precursor pre-GP T544I) of unknown significance within the viral population that exhibited a shift in frequency of at least 10 per cent over six serial passages. In summary, we present the first reptilian cell line that replicates a filovirus at high titers, and for the first time demonstrate a filovirus genus-specific restriction to MARV in a cell line. Our data suggest the possibility that there may be differences between the natural host spectra of ebolaviruses and marburgviruses.

ROP: dumpster diving in RNA-sequencing to find the source of 1 trillion reads across diverse adult human tissues.

High-throughput RNA-sequencing (RNA-seq) technologies provide an unprecedented opportunity to explore the individual transcriptome. Unmapped reads are a large and often overlooked output of standard RNA-seq analyses. Here, we present Read Origin Protocol (ROP), a tool for discovering the source of all reads originating from complex RNA molecules. We apply ROP to samples across 2630 individuals from 54 diverse human tissues. Our approach can account for 99.9% of 1 trillion reads of various read length. Additionally, we use ROP to investigate the functional mechanisms underlying connections between the immune system, microbiome, and disease. ROP is freely available at https://github.com/smangul1/rop/wiki .

Prominent features of the amino acid mutation landscape in cancer.

Cancer can be viewed as a set of different diseases with distinctions based on tissue origin, driver mutations, and genetic signatures. Accordingly, each of these distinctions have been used to classify cancer subtypes and to reveal common features. Here, we present a different analysis of cancer based on amino acid mutation signatures. Non-negative Matrix Factorization and principal component analysis of 29 cancers revealed six amino acid mutation signatures, including four signatures that were dominated by either arginine to histidine (Arg>His) or glutamate to lysine (Glu>Lys) mutations. Sample-level analyses reveal that while some cancers are heterogeneous, others are largely dominated by one type of mutation. Using a non-overlapping set of samples from the COSMIC somatic mutation database, we validate five of six mutation signatures, including signatures with prominent arginine to histidine (Arg>His) or glutamate to lysine (Glu>Lys) mutations. This suggests that our classification of cancers based on amino acid mutation patterns may provide avenues of inquiry pertaining to specific protein mutations that may generate novel insights into cancer biology.

Cancer-associated arginine-to-histidine mutations confer a gain in pH sensing to mutant proteins.

The intracellular pH (pHi) of most cancers is constitutively higher than that of normal cells and enhances proliferation and cell survival. We found that increased pHi enabled the tumorigenic behaviors caused by somatic arginine-to-histidine mutations, which are frequent in cancer and confer pH sensing not seen with wild-type proteins. Experimentally raising the pHi increased the activity of R776H mutant epidermal growth factor receptor (EGFR-R776H), thereby increasing proliferation and causing transformation in fibroblasts. An Arg-to-Gly mutation did not confer these effects. Molecular dynamics simulations of EGFR suggested that decreased protonation of His776 at high pH causes conformational changes in the αC helix that may stabilize the active form of the kinase. An Arg-to-His, but not Arg-to-Lys, mutation in the transcription factor p53 (p53-R273H) decreased its transcriptional activity and attenuated the DNA damage response in fibroblasts and breast cancer cells with high pHi. Lowering pHi attenuated the tumorigenic effects of both EGFR-R776H and p53-R273H. Our data suggest that some somatic mutations may confer a fitness advantage to the higher pHi of cancer cells.

Statistical inference of a convergent antibody repertoire response to influenza vaccine.

BACKGROUND: Vaccines dramatically affect an individual's adaptive immune system and thus provide an excellent means to study human immunity. Upon vaccination, the B cells that express antibodies (Abs) that happen to bind the vaccine are stimulated to proliferate and undergo mutagenesis at their Ab locus. This process may alter the composition of B cell lineages within an individual, which are known collectively as the antibody repertoire (AbR). Antibodies are also highly expressed in whole blood, potentially enabling RNA sequencing (RNA-seq) technologies to query this diversity. Less is known about the diversity of AbR responses across individuals to a given vaccine and if individuals tend to yield a similar response to the same antigenic stimulus. METHODS: Here we implement a bioinformatic pipeline that extracts the AbR information from a time-series RNA-seq dataset of five patients who were administered a seasonal trivalent influenza vaccine (TIV). We harness the detailed time-series nature of this dataset and use methods based in functional data analysis (FDA) to identify the Abs that respond to the vaccine. We then design and implement rigorous statistical tests in order to ask whether or not these patients exhibit a convergent AbR response to the same TIV. RESULTS: We find that high-resolution time-series data can be used to help identify the Abs that respond to an antigenic stimulus and that this response can exhibit a convergent nature across patients inoculated with the same vaccine. However, correlations in AbR diversity among individuals prior to inoculation can confound inference of a convergent signal unless it is taken into account. CONCLUSIONS: We developed a framework to identify the elements of an AbR that respond to an antigen. This information could be used to understand the diversity of different immune responses in different individuals, as well as to gauge the effectiveness of the immune response to a given stimulus within an individual. We also present a framework for testing a convergent hypothesis between AbRs; a hypothesis that is more difficult to test than previously appreciated. Our discovery of a convergent signal suggests that similar epitopes do select for antibodies with similar sequence characteristics.

BMI1 represses Ink4a/Arf and Hox genes to regulate stem cells in the rodent incisor.

The polycomb group gene Bmi1 is required for maintenance of adult stem cells in many organs. Inactivation of Bmi1 leads to impaired stem cell self-renewal due to deregulated gene expression. One critical target of BMI1 is Ink4a/Arf, which encodes the cell-cycle inhibitors p16(Ink4a) and p19(Arf). However, deletion of Ink4a/Arf only partially rescues Bmi1-null phenotypes, indicating that other important targets of BMI1 exist. Here, using the continuously growing mouse incisor as a model system, we report that Bmi1 is expressed by incisor stem cells and that deletion of Bmi1 resulted in fewer stem cells, perturbed gene expression and defective enamel production. Transcriptional profiling revealed that Hox expression is normally repressed by BMI1 in the adult, and functional assays demonstrated that BMI1-mediated repression of Hox genes preserves the undifferentiated state of stem cells. As Hox gene upregulation has also been reported in other systems when Bmi1 is inactivated, our findings point to a general mechanism whereby BMI1-mediated repression of Hox genes is required for the maintenance of adult stem cells and for prevention of inappropriate differentiation.