The BioMart Community Portal (www.biomart.org) is a community-driven effort to provide a unified interface to biomedical databases that are distributed worldwide. The portal provides access to numerous database projects supported by 30 scientific organizations. It includes over 800 different biological datasets spanning genomics, proteomics, model organisms, cancer data, ontology information and more. All resources available through the portal are independently administered and funded by their host organizations. The BioMart data federation technology provides a unified interface to all the available data. The latest version of the portal comes with many new databases that have been created by our ever-growing community. It also comes with better support and extensibility for data analysis and visualization tools. A new addition to our toolbox, the enrichment analysis tool is now accessible through graphical and web service interface. The BioMart community portal averages over one million requests per day. Building on this level of service and the wealth of information that has become available, the BioMart Community Portal has introduced a new, more scalable and cheaper alternative to the large data stores maintained by specialized organizations.

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.

Lung cancer is a highly heterogeneous disease in terms of both underlying genetic lesions and response to therapeutic treatments. We performed deep whole-genome sequencing and transcriptome sequencing on 19 lung cancer cell lines and three lung tumor/normal pairs. Overall, our data show that cell line models exhibit similar mutation spectra to human tumor samples. Smoker and never-smoker cancer samples exhibit distinguishable patterns of mutations. A number of epigenetic regulators, including KDM6A, ASH1L, SMARCA4, and ATAD2, are frequently altered by mutations or copy number changes. A systematic survey of splice-site mutations identified 106 splice site mutations associated with cancer specific aberrant splicing, including mutations in several known cancer-related genes. RAC1b, an isoform of the RAC1 GTPase that includes one additional exon, was found to be preferentially up-regulated in lung cancer. We further show that its expression is significantly associated with sensitivity to a MAP2K (MEK) inhibitor PD-0325901. Taken together, these data present a comprehensive genomic landscape of a large number of lung cancer samples and further demonstrate that cancer-specific alternative splicing is a widespread phenomenon that has potential utility as therapeutic biomarkers. The detailed characterizations of the lung cancer cell lines also provide genomic context to the vast amount of experimental data gathered for these lines over the decades, and represent highly valuable resources for cancer biology.

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.

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.

Many subtypes of cortical interneurons (CINs) are found in adult mouse cortices, but the mechanism generating their diversity remains elusive. We performed single-cell RNA sequencing on the mouse embryonic medial ganglionic eminence (MGE), the major birthplace for CINs, and on MGE-like cells differentiated from embryonic stem cells. Two distinct cell types were identified as proliferating neural progenitors and immature neurons, both of which comprised sub-populations. Although lineage development of MGE progenitors was reconstructed and immature neurons were characterized as GABAergic, cells that might correspond to precursors of different CINs were not identified. A few non-neuronal cell types were detected, including microglia. In vitro MGE-like cells resembled bona fide MGE cells but expressed lower levels of Foxg1 and Epha4. Together, our data provide detailed understanding of the embryonic MGE developmental program and suggest how CINs are specified.

Baboons are a widely used nonhuman primate model for biomedical, evolutionary, and basic genetics research. Despite this importance, the genomic resources for baboons are limited. In particular, the current baboon reference genome Panu_3.0 is a highly fragmented, reference-guided (i.e., not fully de novo) assembly, and its poor quality inhibits our ability to conduct downstream genomic analyses.

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.

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.

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.

Given the frequency with which MAP kinase signaling is dysregulated in cancer, much effort has been focused on inhibiting RAS signaling for therapeutic benefit. KSR1, a pseudokinase that interacts with RAF, is a potential target; it was originally cloned in screens for suppressors of constitutively active RAS, and its deletion prevents RAS-mediated transformation of mouse embryonic fibroblasts. In this work, we used a genetically engineered mouse model of pancreatic cancer to assess whether KSR1 deletion would influence tumor development in the setting of oncogenic RAS. We found that Ksr1-/- mice on this background had a modest but significant improvement in all-cause morbidity compared to Ksr1+/+ and Ksr1+/- cohorts. Ksr1-/- mice, however, still developed tumors, and precursor pancreatic intraepithelial neoplastic (PanIN) lesions were detected within a similar timeframe compared to Ksr1+/+ mice. No significant differences in pERK expression or in proliferation were noted. RNA sequencing also did not reveal any unique genetic signature in Ksr1-/- tumors. Further studies will be needed to determine whether and in what settings KSR inhibition may be clinically useful.

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.

Next-generation sequencings platforms coupled with advanced bioinformatic tools enable re-sequencing of the human genome at high-speed and large cost savings. We compare sequencing platforms from Roche/454(GS FLX), Illumina/HiSeq (HiSeq 2000), and Life Technologies/SOLiD (SOLiD 3 ECC) for their ability to identify single nucleotide substitutions in whole genome sequences from the same human sample. We report on significant GC-related bias observed in the data sequenced on Illumina and SOLiD platforms. The differences in the variant calls were investigated with regards to coverage, and sequencing error. Some of the variants called by only one or two of the platforms were experimentally tested using mass spectrometry; a method that is independent of DNA sequencing. We establish several causes why variants remained unreported, specific to each platform. We report the indel called using the three sequencing technologies and from the obtained results we conclude that sequencing human genomes with more than a single platform and multiple libraries is beneficial when high level of accuracy is required.

Biological studies involve a growing number of distinct high-throughput experiments to characterize samples of interest. There is a lack of methods to visualize these different genomic datasets in a versatile manner. In addition, genomic data analysis requires integrated visualization of experimental data along with constantly changing genomic annotation and statistical analyses.

Gastric cancer is the second leading cause of worldwide cancer mortality, yet the underlying genomic alterations remain poorly understood. Here we perform exome and transcriptome sequencing and SNP array assays to characterize 51 primary gastric tumours and 32 cell lines. Meta-analysis of exome data and previously published data sets reveals 24 significantly mutated genes in microsatellite stable (MSS) tumours and 16 in microsatellite instable (MSI) tumours. Over half the patients in our collection could potentially benefit from targeted therapies. We identify 55 splice site mutations accompanied by aberrant splicing products, in addition to mutation-independent differential isoform usage in tumours. ZAK kinase isoform TV1 is preferentially upregulated in gastric tumours and cell lines relative to normal samples. This pattern is also observed in colorectal, bladder and breast cancers. Overexpression of this particular isoform activates multiple cancer-related transcription factor reporters, while depletion of ZAK in gastric cell lines inhibits proliferation. These results reveal the spectrum of genomic and transcriptomic alterations in gastric cancer, and identify isoform-specific oncogenic properties of ZAK.

Whole-genome sequencing (WGS) can identify variants that cause genetic disease, but the time required for sequencing and analysis has been a barrier to its use in acutely ill patients. In the present study, we develop an approach for ultra-rapid nanopore WGS that combines an optimized sample preparation protocol, distributing sequencing over 48 flow cells, near real-time base calling and alignment, accelerated variant calling and fast variant filtration for efficient manual review. Application to two example clinical cases identified a candidate variant in <8 h from sample preparation to variant identification. We show that this framework provides accurate variant calls and efficient prioritization, and accelerates diagnostic clinical genome sequencing twofold compared with previous approaches.

Metastasis is the principal cause of cancer related deaths. Tumor invasion is essential for metastatic spread. However, determinants of invasion are poorly understood. We addressed this knowledge gap by leveraging a unique attribute of kidney cancer. Renal tumors invade into large vessels forming tumor thrombi (TT) that migrate extending sometimes into the heart. Over a decade, we prospectively enrolled 83 ethnically-diverse patients undergoing surgical resection for grossly invasive tumors at UT Southwestern Kidney Cancer Program. In this study, we perform comprehensive histological analyses, integrate multi-region genomic studies, generate in vivo models, and execute functional studies to define tumor invasion and metastatic competence. We find that invasion is not always associated with the most aggressive clone. Driven by immediate early genes, invasion appears to be an opportunistic trait attained by subclones with diverse oncogenomic status in geospatial proximity to vasculature. We show that not all invasive tumors metastasize and identify determinants of metastatic competency. TT associated with metastases are characterized by higher grade, mTOR activation and a particular immune contexture. Moreover, TT grade is a better predictor of metastasis than overall tumor grade, which may have implications for clinical practice.

Tandem gene amplification is a frequent and dynamic source of antibiotic resistance in bacteria. Ongoing expansions and contractions of repeat arrays during population growth are expected to manifest as cell-to-cell differences in copy number (CN). As a result, a clonal bacterial culture could comprise subpopulations of cells with different levels of antibiotic sensitivity that result from variable gene dosage. Despite the high potential for misclassification of heterogenous cell populations as either antibiotic-susceptible or fully resistant in clinical settings, and the concomitant risk of inappropriate treatment, CN distribution among cells has defied analysis. Here, we use the MinION single-molecule nanopore sequencer to uncover CN heterogeneity in clonal populations of Escherichia coli and Acinetobacter baumannii grown from single cells isolated while selecting for resistance to an optimized arylomycin, a member of a recently discovered class of Gram-negative antibiotic. We found that gene amplification of the arylomycin target, bacterial type I signal peptidase LepB, is a mechanism of unstable arylomycin resistance and demonstrate in E. coli that amplification instability is independent of RecA. This instability drives the emergence of a nonuniform distribution of lepB CN among cells with a range of 1 to at least 50 copies of lepB identified in a single clonal population. In sum, this remarkable heterogeneity, and the evolutionary plasticity it fuels, illustrates how gene amplification can enable bacterial populations to respond rapidly to novel antibiotics. This study establishes a rationale for further nanopore-sequencing studies of heterogeneous cell populations to uncover CN variability at single-molecule resolution.

The cellular origin of sporadic pancreatic neuroendocrine tumors (PNETs) is obscure. Hormone expression suggests that these tumors arise from glucagon-producing alpha cells or insulin-producing β cells, but instability in hormone expression prevents linage determination. We utilize loss of hepatic glucagon receptor (GCGR) signaling to drive alpha cell hyperproliferation and tumor formation to identify a cell of origin and dissect mechanisms that drive progression. Using a combination of genetically engineered Gcgr knockout mice and GCGR-inhibiting antibodies, we show that elevated plasma amino acids drive the appearance of a proliferative population of SLC38A5+ embryonic progenitor-like alpha cells in mice. Further, we characterize tumors from patients with rare bi-allelic germline GCGR loss-of-function variants and find prominent tumor-cell-associated expression of the SLC38A5 paralog SLC7A8 as well as markers of active mTOR signaling. Thus, progenitor cells arise from adult alpha cells in response to metabolic signals and, when inductive signals are chronically present, drive tumor initiation.