Sc-compReg enables the comparison of gene regulatory networks between conditions using single-cell data.

The comparison of gene regulatory networks between diseased versus healthy individuals or between two different treatments is an important scientific problem. Here, we propose sc-compReg as a method for the comparative analysis of gene expression regulatory networks between two conditions using single cell gene expression (scRNA-seq) and single cell chromatin accessibility data (scATAC-seq). Our software, sc-compReg, can be used as a stand-alone package that provides joint clustering and embedding of the cells from both scRNA-seq and scATAC-seq, and the construction of differential regulatory networks across two conditions. We apply the method to compare the gene regulatory networks of an individual with chronic lymphocytic leukemia (CLL) versus a healthy control. The analysis reveals a tumor-specific B cell subpopulation in the CLL patient and identifies TOX2 as a potential regulator of this subpopulation.

Massively parallel single-cell chromatin landscapes of human immune cell development and intratumoral T cell exhaustion.

Understanding complex tissues requires single-cell deconstruction of gene regulation with precision and scale. Here, we assess the performance of a massively parallel droplet-based method for mapping transposase-accessible chromatin in single cells using sequencing (scATAC-seq). We apply scATAC-seq to obtain chromatin profiles of more than 200,000 single cells in human blood and basal cell carcinoma. In blood, application of scATAC-seq enables marker-free identification of cell type-specific cis- and trans-regulatory elements, mapping of disease-associated enhancer activity and reconstruction of trajectories of cellular differentiation. In basal cell carcinoma, application of scATAC-seq reveals regulatory networks in malignant, stromal and immune cells in the tumor microenvironment. Analysis of scATAC-seq profiles from serial tumor biopsies before and after programmed cell death protein 1 blockade identifies chromatin regulators of therapy-responsive T cell subsets and reveals a shared regulatory program that governs intratumoral CD8+ T cell exhaustion and CD4+ T follicular helper cell development. We anticipate that scATAC-seq will enable the unbiased discovery of gene regulatory factors across diverse biological systems.

Resolving the full spectrum of human genome variation using Linked-Reads.

Large-scale population analyses coupled with advances in technology have demonstrated that the human genome is more diverse than originally thought. To date, this diversity has largely been uncovered using short-read whole-genome sequencing. However, these short-read approaches fail to give a complete picture of a genome. They struggle to identify structural events, cannot access repetitive regions, and fail to resolve the human genome into haplotypes. Here, we describe an approach that retains long range information while maintaining the advantages of short reads. Starting from ∼1 ng of high molecular weight DNA, we produce barcoded short-read libraries. Novel informatic approaches allow for the barcoded short reads to be associated with their original long molecules producing a novel data type known as "Linked-Reads". This approach allows for simultaneous detection of small and large variants from a single library. In this manuscript, we show the advantages of Linked-Reads over standard short-read approaches for reference-based analysis. Linked-Reads allow mapping to 38 Mb of sequence not accessible to short reads, adding sequence in 423 difficult-to-sequence genes including disease-relevant genes STRC, SMN1, and SMN2 Both Linked-Read whole-genome and whole-exome sequencing identify complex structural variations, including balanced events and single exon deletions and duplications. Further, Linked-Reads extend the region of high-confidence calls by 68.9 Mb. The data presented here show that Linked-Reads provide a scalable approach for comprehensive genome analysis that is not possible using short reads alone.

Haplotyping germline and cancer genomes with high-throughput linked-read sequencing.

Haplotyping of human chromosomes is a prerequisite for cataloguing the full repertoire of genetic variation. We present a microfluidics-based, linked-read sequencing technology that can phase and haplotype germline and cancer genomes using nanograms of input DNA. This high-throughput platform prepares barcoded libraries for short-read sequencing and computationally reconstructs long-range haplotype and structural variant information. We generate haplotype blocks in a nuclear trio that are concordant with expected inheritance patterns and phase a set of structural variants. We also resolve the structure of the EML4-ALK gene fusion in the NCI-H2228 cancer cell line using phased exome sequencing. Finally, we assign genetic aberrations to specific megabase-scale haplotypes generated from whole-genome sequencing of a primary colorectal adenocarcinoma. This approach resolves haplotype information using up to 100 times less genomic DNA than some methods and enables the accurate detection of structural variants.

The noncoding RNAs SNORD50A and SNORD50B bind K-Ras and are recurrently deleted in human cancer.

Small nucleolar RNAs (snoRNAs) are conserved noncoding RNAs best studied as ribonucleoprotein (RNP) guides in RNA modification. To explore their role in cancer, we compared 5,473 tumor-normal genome pairs to identify snoRNAs with frequent copy number loss. The SNORD50A-SNORD50B snoRNA locus was deleted in 10-40% of 12 common cancers, where its loss was associated with reduced survival. A human protein microarray screen identified direct SNORD50A and SNORD50B RNA binding to K-Ras. Loss of SNORD50A and SNORD50B increased the amount of GTP-bound, active K-Ras and hyperactivated Ras-ERK1/ERK2 signaling. Loss of these snoRNAs also increased binding by farnesyltransferase to K-Ras and increased K-Ras prenylation, suggesting that KRAS mutation might synergize with SNORD50A and SNORD50B loss in cancer. In agreement with this hypothesis, CRISPR-mediated deletion of SNORD50A and SNORD50B in KRAS-mutant tumor cells enhanced tumorigenesis, and SNORD50A and SNORD50B deletion and oncogenic KRAS mutation co-occurred significantly in multiple human tumor types. SNORD50A and SNORD50B snoRNAs thus directly bind and inhibit K-Ras and are recurrently deleted in human cancer.

Autoproteolytic Activation of a Symbiosis-regulated Truffle Phospholipase A2.

Fungal phospholipases are members of the fungal/bacterial group XIV secreted phospholipases A(2) (sPLA(2)s). TbSP1, the sPLA(2) primarily addressed in this study, is up-regulated by nutrient deprivation and is preferentially expressed in the symbiotic stage of the ectomycorrhizal fungus Tuber borchii. A peculiar feature of this phospholipase and of its ortholog from the black truffle Tuber melanosporum is the presence of a 54-amino acid sequence of unknown functional significance, interposed between the signal peptide and the start of the conserved catalytic core of the enzyme. X-ray diffraction analysis of a recombinant TbSP1 form corresponding to the secreted protein previously identified in T. borchii mycelia revealed a structure comprising the five α-helices that form the phospholipase catalytic module but lacking the N-terminal 54 amino acids. This finding led to a series of functional studies that showed that TbSP1, as well as its T. melanosporum ortholog, is a self-processing pro-phospholipase A(2), whose phospholipase activity increases up to 80-fold following autoproteolytic removal of the N-terminal peptide. Proteolytic cleavage occurs within a serine-rich, intrinsically flexible region of TbSP1, does not involve the phospholipase active site, and proceeds via an intermolecular mechanism. Autoproteolytic activation, which also takes place at the surface of nutrient-starved, sPLA(2) overexpressing hyphae, may strengthen and further control the effects of phospholipase up-regulation in response to nutrient deprivation, also in the context of symbiosis establishment and mycorrhiza formation.

Efficient electron transfer in a protein network lacking specific interactions.

In many biochemical processes, proteins need to bind partners amidst a sea of other molecules. Generally, partner selection is achieved by formation of a single-orientation complex with well-defined, short-range interactions. We describe a protein network that functions effectively in a metabolic electron transfer process but lacks such specific interactions. The soil bacterium Paracoccus denitrificans oxidizes a variety of compounds by channeling electrons into the main respiratory pathway. Upon conversion of methylamine by methylamine dehydrogenase, electrons are transported to the terminal oxidase to reduce molecular oxygen. Steady-state kinetic measurements and NMR experiments demonstrate a remarkable number of possibilities for the electron transfer, involving the cupredoxin amicyanin as well as four c-type cytochromes. The observed interactions appear to be governed exclusively by the electrostatic nature of each of the proteins. It is concluded that Paracoccus provides a pool of cytochromes for efficient electron transfer via weak, ill-defined interactions, in contrast with the view that functional biochemical interactions require well-defined molecular interactions. It is proposed that the lack of requirement for specificity in these interactions might facilitate the integration of new metabolic pathways.

Amicyanin transfers electrons from methylamine dehydrogenase to cytochrome c-551i via a ping-pong mechanism, not a ternary complex.

The first crystal structure of a ternary redox protein complex was comprised of the enzyme methylamine dehydrogenase (MADH) and two electron transfer proteins, amicyanin and cytochrome c-551i from Paracoccus denitrificans [Chen et al. Science 1994, 264, 86-90]. The arrangement of the proteins suggested possible electron transfer from the active site of MADH via the amicyanin copper ion to the cytochrome heme iron, although the distance between the metals is large. We studied the interactions between these proteins in solution. A titration followed by NMR spectroscopy shows that amicyanin binds cytochrome c-551i. The interface comprises the hydrophobic and positive patches of amicyanin, not the binding site observed in the ternary complex. NMR experiments further show that amicyanin binds tightly to MADH with an interface that matches the one observed in the crystal structure and that mostly overlaps with the binding site for cytochrome c-551i. Upon addition of cytochrome c-551i, no changes in the NMR spectrum of MADH-bound amicyanin are observed, suggesting that a possible interaction of the cytochrome with the binary complex must be very weak, with a dissociation constant higher than 2 mM. Reconstitution of the entire redox chain in vitro demonstrates that amicyanin can react rapidly with cytochrome c-551i, but that association of amicyanin with MADH inhibits this reaction. It is concluded that electron transfer from MADH to cytochrome c-551i does not involve a ternary complex but occurs via a ping-pong mechanism in which amicyanin uses the same interface for the reactions with MADH and cytochrome c-551i.