Search

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., Changes in cell state underlie the difference between health and disease. Here, the authors propose a computational framework for the integration of gene expression and chromatin-accessibility data from single cells to identify differences in gene regulation in cell types across two conditions.

Understanding complex tissues requires single-cell deconstruction of gene regulation with precision and scale. Here we present a massively parallel droplet-based platform for mapping transposase-accessible chromatin in tens of thousands of single cells per sample (scATAC-seq). We obtain and analyze chromatin profiles of over 200,000 single cells in two primary human systems. In blood, scATAC-seq allows marker-free identification of cell type-specificcis- andtrans-regulatory elements, mapping of disease-associated enhancer activity, and reconstruction of trajectories of differentiation from progenitors to diverse and rare immune cell types. In basal cell carcinoma, scATAC-seq reveals regulatory landscapes of malignant, stromal, and immune cell types in the tumor microenvironment. Moreover, scATAC-seq of serial tumor biopsies before and after PD-1 blockade allows identification of chromatin regulators and differentiation trajectories of therapy-responsive intratumoral T cell subsets, revealing a shared regulatory program driving CD8+T cell exhaustion and CD4+T follicular helper cell development. We anticipate that droplet-based single-cell chromatin accessibility will provide a broadly applicable means of identifying regulatory factors and elements that underlie cell type and function.

Small nucleolar RNAs (snoRNAs) are conserved noncoding RNAs best studied as ribonucleoprotein (RNP) guides in RNA modification1,2. 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.

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.

Large-scale population based 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, standard short-read approaches, used primarily due to accuracy, throughput and costs, fail to give a complete picture of a genome. They struggle to identify large, balanced structural events, cannot access repetitive regions of the genome and fail to resolve the human genome into its two haplotypes. Here we describe an approach that retains long range information while harnessing the advantages of short reads. Starting from only ∼1ng of DNA, we produce barcoded short read libraries. The use of novel informatic approaches allows for the barcoded short reads to be associated with the long molecules of origin producing a novel datatype known as ‘Linked-Reads’. This approach allows for simultaneous detection of small and large variants from a single Linked-Read library. We have previously demonstrated the utility of whole genome Linked-Reads (lrWGS) for performing diploid, de novo assembly of individual genomes (Weisenfeld et al. 2017). In this manuscript, we show the advantages of Linked-Reads over standard short read approaches for reference based analysis. We demonstrate the ability of Linked-Reads to reconstruct megabase scale haplotypes and to recover parts of the genome that are typically inaccessible to short reads, including phenotypically important genes such as STRC, SMN1 and SMN2. We demonstrate the ability of both lrWGS and Linked-Read Whole Exome Sequencing (lrWES) to identify complex structural variations, including balanced events, single exon deletions, and single exon duplications. 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 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.

Fungal phospholipases are members of the fungal/bacterial group XIV secreted phospholipases A2 (sPLA2s). TbSP1, the sPLA2 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 A2, 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, sPLA2 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. Background: TbSP1 is a phospholipase A2 strongly up-regulated during the symbiotic phase of the truffle Tuber borchii. Results: An activated enzyme species composed of five α-helices is generated by self-proteolysis through an intermolecular reaction. Conclusion: TbSP1 autoproteolysis is a site-specific post-translational modification not involving the phospholipase active site. Significance: Autoproteolytic activation is described for the first time for a microbial PLA2, with possible implications for symbiosis establishment.

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.

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.