Uveal melanoma (UM) is the most common primary cancer of the eye, resulting not only in vision loss, but also in metastatic death. This study attempts to identify changes in the patterns of gene expression that lead to malignant transformation and proliferation of normal uveal melanocytes (UVM) using the Suppressive Subtractive Hybridization (SSH) technique.

Adenosine-to-inosine (A-to-I) editing modifies RNA transcripts from their genomic blueprint. A prerequisite for this process is a double-stranded RNA (dsRNA) structure. Such dsRNAs are formed as part of the microRNA (miRNA) maturation process, and it is therefore expected that miRNAs are affected by A-to-I editing. Editing of miRNAs has the potential to add another layer of complexity to gene regulation pathways, especially if editing occurs within the miRNA-mRNA recognition site. Thus, it is of interest to study the extent of this phenomenon. Current reports in the literature disagree on its extent; while some reports claim that it may be widespread, others deem the reported events as rare. Utilizing a next-generation sequencing (NGS) approach supplemented by an extensive bioinformatic analysis, we were able to systematically identify A-to-I editing events in mature miRNAs derived from human brain tissues. Our algorithm successfully identified many of the known editing sites in mature miRNAs and revealed 17 novel human sites, 12 of which are in the recognition sites of the miRNAs. We confirmed most of the editing events using in vitro ADAR overexpression assays. The editing efficiency of most sites identified is very low. Similar results are obtained for publicly available data sets of mouse brain-regions tissues. Thus, we find that A-to-I editing does alter several miRNAs, but it is not widespread.

The B cell response to Salmonella typhimurium (STm) occurs massively at extrafollicular sites, without notable germinal centers (GCs). Little is known in terms of its specificity. To expand the knowledge of antigen targets, we screened plasmablast (PB)-derived monoclonal antibodies (mAbs) for Salmonella specificity, using ELISA, flow cytometry, and antigen microarray. Only a small fraction (0.5%-2%) of the response appeared to be Salmonella-specific. Yet, infection of mice with limited B cell receptor (BCR) repertoires impaired the response, suggesting that BCR specificity was important. We showed, using laser microdissection, that somatic hypermutation (SHM) occurred efficiently at extrafollicular sites leading to affinity maturation that in turn led to detectable STm Ag-binding. These results suggest a revised vision of how clonal selection and affinity maturation operate in response to Salmonella. Clonal selection initially is promiscuous, activating cells with virtually undetectable affinity, yet SHM and selection occur during the extrafollicular response yielding higher affinity, detectable antibodies.

The B cell response to Ehrlichia muris is dominated by plasmablasts (PBs), with few-if any-germinal centers (GCs), yet it generates protective immunoglobulin M (IgM) memory B cells (MBCs) that express the transcription factor T-bet and harbor V-region mutations. Because Ehrlichia prominently infects the liver, we investigated the nature of liver B cell response and that of the spleen. B cells within infected livers proliferated and underwent somatic hypermutation (SHM). Vh-region sequencing revealed trafficking of clones between the spleen and liver and often subsequent local clonal expansion and intraparenchymal localization of T-bet+ MBCs. T-bet+ MBCs expressed MBC subset markers CD80 and PD-L2. Many T-bet+ MBCs lacked CD11b or CD11c expression but had marginal zone (MZ) B cell phenotypes and colonized the splenic MZ, revealing T-bet+ MBC plasticity. Hence, liver and spleen are generative sites of B cell responses, and they include V-region mutation and result in liver MBC localization.

Snakebite envenomings represent a neglected public health issue in many parts of the rural tropical world. Animal-derived antivenoms have existed for more than a hundred years and are effective in neutralizing snake venom toxins when timely administered. However, the low immunogenicity of many small but potent snake venom toxins represents a challenge for obtaining a balanced immune response against the medically relevant components of the venom. Here, we employ high-throughput sequencing of the immunoglobulin (Ig) transcriptome of mice immunized with a three-finger toxin and a phospholipase A2 from the venom of the Central American coral snake, Micrurus nigrocinctus. Although exploratory in nature, our indicate results showed that only low frequencies of mRNA encoding IgG isotypes, the most relevant isotype for therapeutic purposes, were present in splenocytes of five mice immunized with 6 doses of the two types of toxins over 90 days. Furthermore, analysis of Ig heavy chain transcripts showed that no particular combination of variable (V) and joining (J) gene segments had been selected in the immunization process, as would be expected after a strong humoral immune response to a single antigen. Combined with the titration of toxin-specific antibodies in the sera of immunized mice, these data support the low immunogenicity of three-finger toxins and phospholipases A2 found in M. nigrocinctusvenoms, and highlight the need for future studies analyzing the complexity of antibody responses to toxins at the molecular level.

Inference of germline polymorphisms in immunoglobulin genes from B cell receptor repertoires is complicated by somatic hypermutations, sequencing/PCR errors, and by varying length of reference alleles. The light chain inference is particularly challenging owing to large gene duplications and absence of D genes. We analyzed the light chain cDNA sequences from naïve B cell receptor repertoires from 100 individuals. We optimized light chain allele inference by tweaking parameters of the TIgGER functions, extending the germline reference sequences, and establishing mismatch frequency patterns at polymorphic positions to filter out false-positive candidates. We identified 48 previously unreported variants of light chain variable genes. We selected 14 variants for validation and successfully validated 11 by Sanger sequencing. Clustering of light chain 5'UTR, L-PART1, and L-PART2 revealed partial intron retention in 11 kappa and 9 lambda V alleles. Our results provide insight into germline variation in human light chain immunoglobulin loci.

GATA binding protein 3 (GATA3) is a regulator of mammary luminal cell differentiation, and an estrogen receptor (ER) associated marker in breast cancer. Tumor suppressor functions of GATA3 have been demonstrated primarily in basal-like breast cancers. Here, we focused on its function in luminal breast cancer, where GATA3 is frequently mutated, and its levels are significantly elevated.

Long noncoding RNAs (lncRNAs) are major regulators of many cellular processes including cell cycle progression and tumorigenesis. In this study, we identify a novel lncRNA, MA-linc1, and reveal its effects on cell cycle progression and cancer growth. Inhibition of MA-linc1 expression alters cell cycle distribution, leading to a decrease in the number of G1 cells and a concomitant increase in all other stages of the cell cycle, and in particular G2/M, suggesting its involvement in the regulation of M phase. Accordingly, knock down of MA-linc1 inhibits M phase exit upon release from a mitotic block. We further demonstrate that MA-linc1 predominantly functions in cis to repress expression of its neighboring gene, Purα, which is often deleted in human cancers and whose ectopic expression inhibits cell cycle progression. Knock down of Purα partially rescues the MA-linc1 dependent inhibition of M phase exit. In agreement with its suggested role in M phase, inhibition of MA-linc1 enhances apoptotic cell death induced by the antimitotic drug, Paclitaxel and this enhancement of apoptosis is rescued by Purα knockdown. Furthermore, high levels of MA-linc1 are associated with reduced survival in human breast and lung cancer patients.Taken together, our data identify MA-linc1 as a novel lncRNA regulator of cell cycle and demonstrate its potential role in cancer progression and treatment.

Broadly neutralizing HIV antibodies (bnAbs) are typically highly somatically mutated, raising doubts as to whether they can be elicited by vaccination. We used 454 sequencing and designed a novel phylogenetic method to model lineage evolution of the bnAbs PGT121-134 and found a positive correlation between the level of somatic hypermutation (SHM) and the development of neutralization breadth and potency. Strikingly, putative intermediates were characterized that show approximately half the mutation level of PGT121-134 but were still capable of neutralizing roughly 40-80% of PGT121-134 sensitive viruses in a 74-virus panel at median titers between 15- and 3-fold higher than PGT121-134. Such antibodies with lower levels of SHM may be more amenable to elicitation through vaccination while still providing noteworthy coverage. Binding characterization indicated a preference of inferred intermediates for native Env binding over monomeric gp120, suggesting that the PGT121-134 lineage may have been selected for binding to native Env at some point during maturation. Analysis of glycan-dependent neutralization for inferred intermediates identified additional adjacent glycans that comprise the epitope and suggests changes in glycan dependency or recognition over the course of affinity maturation for this lineage. Finally, patterns of neutralization of inferred bnAb intermediates suggest hypotheses as to how SHM may lead to potent and broad HIV neutralization and provide important clues for immunogen design.

Many microbial communities are characterized by high genetic diversity. 16S ribosomal RNA sequencing can determine community members, and metagenomics can determine the functional diversity, but resolving the functional role of individual cells in high throughput remains an unsolved challenge. Here, we describe epicPCR (Emulsion, Paired Isolation and Concatenation PCR), a new technique that links functional genes and phylogenetic markers in uncultured single cells, providing a throughput of hundreds of thousands of cells with costs comparable to one genomic library preparation. We demonstrate the utility of our technique in a natural environment by profiling a sulfate-reducing community in a freshwater lake, revealing both known sulfate reducers and discovering new putative sulfate reducers. Our method is adaptable to any conserved genetic trait and translates genetic associations from diverse microbial samples into a sequencing library that answers targeted ecological questions. Potential applications include identifying functional community members, tracing horizontal gene transfer networks and mapping ecological interactions between microbial cells.

The mechanistic target of rapamycin, (mTOR) kinase plays a pivotal role in controlling critical cellular growth and survival pathways, and its aberrant induction is implicated in cancer pathogenesis. Therefore, suppression of active mTOR signaling has been of great interest to researchers; several mTOR inhibitors have been discovered to date. Ethanol (EtOH), similar to pharmacologic mTOR inhibitors, has been shown to suppress the mTOR signaling pathway, though in a non-catalytic manner. Despite population studies showing that the consumption of EtOH has a protective effect against hematological malignancies, the mechanisms behind EtOH's modulation of mTOR activity in cells and its downstream consequences are largely unknown. Here we evaluated the effects of EtOH on the mTOR pathway, in comparison to the active-site mTOR inhibitor INK128, and compared translatome analysis of their downstream effects in diffuse large B-cell lymphoma (DLBCL).

The adaptive immune receptor repertoire (AIRR) contains information on an individuals' immune past, present and potential in the form of the evolving sequences that encode the B cell receptor (BCR) repertoire. AIRR sequencing (AIRR-seq) studies rely on databases of known BCR germline variable (V), diversity (D), and joining (J) genes to detect somatic mutations in AIRR-seq data via comparison to the best-aligning database alleles. However, it has been shown that these databases are far from complete, leading to systematic misidentification of mutated positions in subsets of sample sequences. We previously presented TIgGER, a computational method to identify subject-specific V gene genotypes, including the presence of novel V gene alleles, directly from AIRR-seq data. However, the original algorithm was unable to detect alleles that differed by more than 5 single nucleotide polymorphisms (SNPs) from a database allele. Here we present and apply an improved version of the TIgGER algorithm which can detect alleles that differ by any number of SNPs from the nearest database allele, and can construct subject-specific genotypes with minimal prior information. TIgGER predictions are validated both computationally (using a leave-one-out strategy) and experimentally (using genomic sequencing), resulting in the addition of three new immunoglobulin heavy chain V (IGHV) gene alleles to the IMGT repertoire. Finally, we develop a Bayesian strategy to provide a confidence estimate associated with genotype calls. All together, these methods allow for much higher accuracy in germline allele assignment, an essential step in AIRR-seq studies.

Analysis of antibody repertoires by high-throughput sequencing is of major importance in understanding adaptive immune responses. Our knowledge of variations in the genomic loci encoding immunoglobulin genes is incomplete, resulting in conflicting VDJ gene assignments and biased genotype and haplotype inference. Haplotypes can be inferred using IGHJ6 heterozygosity, observed in one third of the people. Here, we propose a robust novel method for determining VDJ haplotypes by adapting a Bayesian framework. Our method extends haplotype inference to IGHD- and IGHV-based analysis, enabling inference of deletions and copy number variations in the entire population. To test this method, we generated a multi-individual data set of naive B-cell repertoires, and found allele usage bias, as well as a mosaic, tiled pattern of deleted IGHD and IGHV genes. The inferred haplotypes may have clinical implications for genetic disease predispositions. Our findings expand the knowledge that can be extracted from antibody repertoire sequencing data.

Hepatitis C virus (HCV) is a major public health concern, with over 70 million people infected worldwide, who are at risk for developing life-threatening liver disease. No vaccine is available, and immunity against the virus is not well-understood. Following the acute stage, HCV usually causes chronic infections. However, ~30% of infected individuals spontaneously clear the virus. Therefore, using HCV as a model for comparing immune responses between spontaneous clearer (SC) and chronically infected (CI) individuals may empower the identification of mechanisms governing viral infection outcomes. Here, we provide the first in-depth analysis of adaptive immune receptor repertoires in individuals with current or past HCV infection. We demonstrate that SC individuals, in contrast to CI patients, develop clusters of antibodies with distinct properties. These antibodies' characteristics were used in a machine learning framework to accurately predict infection outcome. Using combinatorial antibody phage display library technology, we identified HCV-specific antibody sequences. By integrating these data with the repertoire analysis, we constructed two antibodies characterized by high neutralization breadth, which are associated with clearance. This study provides insight into the nature of effective immune response against HCV and demonstrates an innovative approach for constructing antibodies correlating with successful infection clearance. It may have clinical implications for prognosis of the future status of infection, and the design of effective immunotherapies and a vaccine for HCV.

Although structural studies of individual T cell receptors (TCRs) have revealed important roles for both the α and β chain in directing MHC and antigen recognition, repertoire-level immunogenomic analyses have historically examined the β chain alone. To determine the amount of useful information about TCR repertoire function encoded within αβ pairings, we analyzed paired TCR sequences from nearly 100,000 unique CD4+ and CD8+ T cells captured using two different high-throughput, single-cell sequencing approaches. Our results demonstrate little overlap in the healthy CD4+ and CD8+ repertoires, with shared TCR sequences possessing significantly shorter CDR3 sequences corresponding to higher generation probabilities. We further utilized tools from information theory and machine learning to show that while α and β chains are only weakly associated with lineage, αβ pairings appear to synergistically drive TCR-MHC interactions. Vαβ gene pairings were found to be the TCR feature most informative of T cell lineage, supporting the existence of germline-encoded paired αβ TCR-MHC interaction motifs. Finally, annotating our TCR pairs using a database of sequences with known antigen specificities, we demonstrate that approximately a third of the T cells possess α and β chains that each recognize different known antigens, suggesting that αβ pairing is critical for the accurate inference of repertoire functionality. Together, these findings provide biological insight into the functional implications of αβ pairing and highlight the utility of single-cell sequencing in immunogenomics.

Germline variations in immunoglobulin genes influence the repertoire of B cell receptors and antibodies, and such polymorphisms may impact disease susceptibility. However, the knowledge of the genomic variation of the immunoglobulin loci is scarce. Here, we report 25 potential novel germline IGHV alleles as inferred from rearranged naïve B cell cDNA repertoires of 98 individuals. Thirteen novel alleles were selected for validation, out of which ten were successfully confirmed by targeted amplification and Sanger sequencing of non-B cell DNA. Moreover, we detected a high degree of variability upstream of the V-REGION in the 5'UTR, L-PART1 and L-PART2 sequences, and found that identical V-REGION alleles can differ in upstream sequences. Thus, we have identified a large genetic variation not only in the V-REGION but also in the upstream sequences of IGHV genes. Our findings provide a new perspective for annotating immunoglobulin repertoire sequencing data.

The partial success of tumor immunotherapy induced by checkpoint blockade, which is not antigen-specific, suggests that the immune system of some patients contain antigen receptors able to specifically identify tumor cells. Here we focused on T-cell receptor (TCR) repertoires associated with spontaneous breast cancer. We studied the alpha and beta chain CDR3 domains of TCR repertoires of CD4 T cells using deep sequencing of cell populations in mice and applied the results to published TCR sequence data obtained from human patients. We screened peripheral blood T cells obtained monthly from individual mice spontaneously developing breast tumors by 5 months. We then looked at identical TCR sequences in published human studies; we used TCGA data from tumors and healthy tissues of 1,256 breast cancer resections and from 4 focused studies including sequences from tumors, lymph nodes, blood and healthy tissues, and from single cell dataset of 3 breast cancer subjects. We now report that mice spontaneously developing breast cancer manifest shared, Public CDR3 regions in both their alpha and beta and that a significant number of women with early breast cancer manifest identical CDR3 sequences. These findings suggest that the development of breast cancer is associated, across species, with biomarker, exclusive TCR repertoires.

Hepatocellular carcinoma (HCC) is a model of a diverse spectrum of cancers because it is induced by well-known etiologies, mainly hepatitis C virus (HCV) and hepatitis B virus. Here, we aimed to identify HCV-specific mutational signatures and explored the link between the HCV-related regional variation in mutations rates and HCV-induced alterations in genome-wide chromatin organization.