Technological advances have enabled transcriptome characterization of cell types at the single-cell level providing new biological insights. New methods that enable simple yet high-throughput single-cell expression profiling are highly desirable.

The human epidermal growth factor receptor (HER) family of tyrosine kinases is deregulated in multiple cancers either through amplification, overexpression, or mutation. ERBB3/HER3, the only member with an impaired kinase domain, although amplified or overexpressed in some cancers, has not been reported to carry oncogenic mutations. Here, we report the identification of ERBB3 somatic mutations in ~11% of colon and gastric cancers. We found that the ERBB3 mutants transformed colonic and breast epithelial cells in a ligand-independent manner. However, the mutant ERBB3 oncogenic activity was dependent on kinase-active ERBB2. Furthermore, we found that anti-ERBB antibodies and small molecule inhibitors effectively blocked mutant ERBB3-mediated oncogenic signaling and disease progression in vivo.

Deregulated HER2 is a target of many approved cancer drugs. We analyzed 111,176 patient tumors and identified recurrent mutations in HER2 transmembrane domain (TMD) and juxtamembrane domain (JMD) that include G660D, R678Q, E693K, and Q709L. Using a saturation mutagenesis screen and testing of patient-derived mutations we found several activating TMD and JMD mutations. Structural modeling and analysis showed that the TMD/JMD mutations function by improving the active dimer interface or stabilizing an activating conformation. Further, we found that HER2 G660D employed asymmetric kinase dimerization for activation and signaling. Importantly, anti-HER2 antibodies and small-molecule kinase inhibitors blocked the activity of TMD/JMD mutants. Consistent with this, a G660D germline mutant lung cancer patient showed remarkable clinical response to HER2 blockade.

Data from the 1000 Genomes project is quite often used as a reference for human genomic analysis. However, its accuracy needs to be assessed to understand the quality of predictions made using this reference. We present here an assessment of the genotyping, phasing, and imputation accuracy data in the 1000 Genomes project. We compare the phased haplotype calls from the 1000 Genomes project to experimentally phased haplotypes for 28 of the same individuals sequenced using the 10X Genomics platform.

ERBB3 is a pseudokinase domain-containing member of the ERBB family of receptor tyrosine kinases (RTKs). Following ligand binding, ERBB receptors homo- or hetero-dimerize, leading to a head-to-tail arrangement of the intracellular kinase domains, where the "receiver" kinase domain of one ERBB is activated by the "activator" domain of the other ERBB in the dimer. In ERBB3, a conserved valine at codon 943 (V943) in the kinase C-terminal domain has been shown to be important for its function as an "activator" kinase in vitro. Here we report a knock-in mouse model where we have modified the endogenous Erbb3 allele to allow for tissue-specific conditional expression of Erbb3 V943R (Erbb3 CKI-V943R ). Additionally, we generated an Erbb3 D850N (Erbb3 CKI-D850N ) conditional knock-in mouse model where the conserved aspartate in the DFG motif of the pseudokinase domain was mutated to abolish any potential residual kinase activity. While Erbb3 D850N/D850N animals developed normally, homozygous Erbb3 V943R/V943R expression during development resulted in embryonic lethality. Further, tissue specific expression of Erbb3 V943R/V943R in the mammary gland epithelium following its activation using MMTV-Cre resulted in delayed elongation of the ductal network during puberty. Single-cell RNA-seq analysis of Erbb3 V943R/V943R mammary glands showed a reduction in a specific subset of fibrinogen-producing luminal epithelial cells.

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the 'venom-ome' and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 'venom-ome-specific toxins' (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery.

Argonaute (AGO) proteins execute microRNA (miRNA)-mediated gene silencing. However, it is unclear whether all 4 mammalian AGO proteins (AGO1, AGO2, AGO3, and AGO4) are required for miRNA activity. We generate Ago1, Ago3, and Ago4-deficient mice (Ago134Δ) and find AGO1/3/4 to be redundant for miRNA biogenesis, homeostasis, or function, a role that is carried out by AGO2. Instead, AGO1/3/4 regulate the expansion of type 2 immunity via precursor mRNA splicing in CD4+ T helper (Th) lymphocytes. Gain- and loss-of-function experiments demonstrate that nuclear AGO3 interacts directly with SF3B3, a component of the U2 spliceosome complex, to aid global mRNA splicing, and in particular the isoforms of the gene Nisch, resulting in a dysregulated Nisch isoform ratio. This work uncouples AGO1, AGO3, and AGO4 from miRNA-mediated RNA interference, identifies an AGO3:SF3B3 complex in the nucleus, and reveals a mechanism by which AGO proteins regulate inflammatory diseases.

Obtaining full-length antibody heavy- and light-chain variable regions from individual B cells at scale remains a challenging problem. Here we use high-throughput single-cell B-cell receptor sequencing (scBCR-seq) to obtain accurately paired full-length variable regions in a massively parallel fashion. We sequenced more than 250,000 B cells from rat, mouse and human repertoires to characterize their lineages and expansion. In addition, we immunized rats with chicken ovalbumin and profiled antigen-reactive B cells from lymph nodes of immunized animals. The scBCR-seq data recovered 81% (n = 56/69) of B-cell lineages identified from hybridomas generated from the same set of B cells subjected to scBCR-seq. Importantly, scBCR-seq identified an additional 710 candidate lineages not recovered as hybridomas. We synthesized, expressed and tested 93 clones from the identified lineages and found that 99% (n = 92/93) of the clones were antigen-reactive. Our results establish scBCR-seq as a powerful tool for antibody discovery.

We sequenced the Hyposidra talaca NPV (HytaNPV) double stranded circular DNA genome using PacBio single molecule sequencing technology. We found that the HytaNPV genome is 139,089 bp long with a GC content of 39.6%. It encodes 141 open reading frames (ORFs) including the 37 baculovirus core genes, 25 genes conserved among lepidopteran baculoviruses, 72 genes known in baculovirus, and 7 genes unique to the HytaNPV genome. It is a group II alphabaculovirus that codes for the F protein and lacks the gp64 gene found in group I alphabaculovirus viruses. Using RNA-seq, we confirmed the expression of the ORFs identified in the HytaNPV genome. Phylogenetic analysis showed HytaNPV to be closest to BusuNPV, SujuNPV and EcobNPV that infect other tea pests, Buzura suppressaria, Sucra jujuba, and Ectropis oblique, respectively. We identified repeat elements and a conserved non-coding baculovirus element in the genome. Analysis of the putative promoter sequences identified motif consistent with the temporal expression of the genes observed in the RNA-seq data.

Tissue homeostasis requires long-term lineage fidelity of somatic stem cells. Whether and how age-related changes in somatic stem cells impact the faithful execution of lineage decisions remains largely unknown. Here, we address this question using genome-wide chromatin accessibility and transcriptome analysis as well as single-cell RNA-seq to explore stem-cell-intrinsic changes in the aging Drosophila intestine. These studies indicate that in stem cells of old flies, promoters of Polycomb (Pc) target genes become differentially accessible, resulting in the increased expression of enteroendocrine (EE) cell specification genes. Consistently, we find age-related changes in the composition of the EE progenitor cell population in aging intestines, as well as a significant increase in the proportion of EE-specified intestinal stem cells (ISCs) and progenitors in aging flies. We further confirm that Pc-mediated chromatin regulation is a critical determinant of EE cell specification in the Drosophila intestine. Pc is required to maintain expression of stem cell genes while ensuring repression of differentiation and specification genes. Our results identify Pc group proteins as central regulators of lineage identity in the intestinal epithelium and highlight the impact of age-related decline in chromatin regulation on tissue homeostasis.

Genomic studies performed in cancer patients and tumor-derived cell lines have identified a high frequency of alterations in components of the mammalian switch/sucrose non-fermentable (mSWI/SNF or BAF) chromatin remodeling complex, including its core catalytic subunit, SMARCA4. Cells exhibiting loss of SMARCA4 rely on its paralog, SMARCA2, making SMARCA2 an attractive therapeutic target. Here we report the genomic profiling of solid tumors from 131,668 cancer patients, identifying 9434 patients with one or more SMARCA4 gene alterations. Homozygous SMARCA4 mutations were highly prevalent in certain tumor types, notably non-small cell lung cancer (NSCLC), and associated with reduced survival. The large sample size revealed previously uncharacterized hotspot missense mutations within the SMARCA4 helicase domain. Functional characterization of these mutations demonstrated markedly reduced remodeling activity. Surprisingly, a few SMARCA4 missense variants partially or fully rescued paralog dependency, underscoring that careful selection criteria must be employed to identify patients with inactivating, homozygous SMARCA4 missense mutations who may benefit from SMARCA2-targeted therapy.

Gallbladder cancer (GBC) is an aggressive gastrointestinal malignancy with no approved targeted therapy. Here, we analyze exomes (n = 160), transcriptomes (n = 115), and low pass whole genomes (n = 146) from 167 gallbladder cancers (GBCs) from patients in Korea, India and Chile. In addition, we also sequence samples from 39 GBC high-risk patients and detect evidence of early cancer-related genomic lesions. Among the several significantly mutated genes not previously linked to GBC are ETS domain genes ELF3 and EHF, CTNNB1, APC, NSD1, KAT8, STK11 and NFE2L2. A majority of ELF3 alterations are frame-shift mutations that result in several cancer-specific neoantigens that activate T-cells indicating that they are cancer vaccine candidates. In addition, we identify recurrent alterations in KEAP1/NFE2L2 and WNT pathway in GBC. Taken together, these define multiple targetable therapeutic interventions opportunities for GBC treatment and management.

Members of the mammalian phosphoinositide-3-OH kinase (PI3K) family of proteins are critical regulators of various cellular process including cell survival, growth, proliferation, and motility. Oncogenic activating mutations in the p110alpha catalytic subunit of the heterodimeric p110/p85 PI3K enzyme are frequent in human cancers. Here we show the presence of frequent mutations in p85alpha in colon cancer, a majority of which occurs in the inter-Src homology-2 (iSH2) domain. These mutations uncouple and retain p85alpha's p110-stabilizing activity, while abrogating its p110-inhibitory activity. The p85alpha mutants promote cell survival, AKT activation, anchorage-independent cell growth, and oncogenesis in a p110-dependent manner.

A propensity for rewiring genetic and epigenetic regulatory networks, thus enabling sustained cell proliferation, suppression of apoptosis, and the ability to evade the immune system, is vital to cancer cell propagation. An increased understanding of how this is achieved is critical for identifying or improving therapeutic interventions. In this study, using acute myeloid leukemia (AML) human cell lines and a custom CRISPR/Cas9 screening platform, we identify the H3K9 methyltransferase SETDB1 as a novel, negative regulator of innate immunity. SETDB1 is overexpressed in many cancers, and loss of this gene in AML cells triggers desilencing of retrotransposable elements that leads to the production of double-stranded RNAs (dsRNAs). This is coincident with induction of a type I interferon response and apoptosis through the dsRNA-sensing pathway. Collectively, our findings establish a unique gene regulatory axis that cancer cells can exploit to circumvent the immune system.

To further understand the molecular distinctions between kidney cancer subtypes, we analyzed exome, transcriptome and copy number alteration data from 167 primary human tumors that included renal oncocytomas and non-clear cell renal cell carcinomas (nccRCCs), consisting of papillary (pRCC), chromophobe (chRCC) and translocation (tRCC) subtypes. We identified ten significantly mutated genes in pRCC, including MET, NF2, SLC5A3, PNKD and CPQ. MET mutations occurred in 15% (10/65) of pRCC samples and included previously unreported recurrent activating mutations. In chRCC, we found TP53, PTEN, FAAH2, PDHB, PDXDC1 and ZNF765 to be significantly mutated. Gene expression analysis identified a five-gene set that enabled the molecular classification of chRCC, renal oncocytoma and pRCC. Using RNA sequencing, we identified previously unreported gene fusions, including ACTG1-MITF fusion. Ectopic expression of the ACTG1-MITF fusion led to cellular transformation and induced the expression of downstream target genes. Finally, we observed upregulation of the anti-apoptotic factor BIRC7 in MiTF-high RCC tumors, suggesting a potential therapeutic role for BIRC7 inhibitors.

Maturity-onset diabetes of the young (MODY) is an early-onset, autosomal dominant form of non-insulin dependent diabetes. Genetic diagnosis of MODY can transform patient management. Earlier data on the genetic predisposition to MODY have come primarily from familial studies in populations of European origin.

Antibodies from B-cell clonal lineages share sequence and structural properties as well as epitope specificity. Clonally unrelated antibodies can similarly share sequence and specificity properties and are said to be convergent. Convergent antibody responses against several antigens have been described in humans and mice and include different classes of shared sequence features. In particular, some antigens and epitopes can induce convergent responses of clonally unrelated antibodies with restricted heavy (VH) and light (VL) chain variable region germline segment usage without similarity in the heavy chain third complementarity-determining region (CDR H3), a critical specificity determinant. Whether these V germline segment-restricted responses reflect a general epitope specificity restriction of antibodies with shared VH/VL pairing is not known. Here, we investigated this question by determining patterns of antigen binding competition between clonally unrelated antigen-specific rat antibodies from paired-chain deep sequencing datasets selected based solely on VH/VL pairing. We found that antibodies with shared VH/VL germline segment pairings but divergent CDR H3 sequences almost invariably have restricted epitope specificity indicated by shared binding competition patterns. This epitope restriction included 82 of 85 clonally unrelated antibodies with 13 different VH/VL pairings binding in 8 epitope groups in 2 antigens. The corollary that antibodies with shared VH/VL pairing and epitope-restricted binding can accommodate widely divergent CDR H3 sequences was confirmed by in vitro selection of variants of anti-human epidermal growth factor receptor 2 antibodies known to mediate critical antigen interactions through CDR H3. Our results show that restricted epitope specificity determined by VH/VL germline segment pairing is a general property of rodent antigen-specific antibodies.