Targeted transcriptome analysis using synthetic long read sequencing uncovers isoform reprograming in the progression of colon cancer.

The characterization of human gene expression is limited by short read lengths, high error rates and large input requirements. Here, we used a synthetic long read (SLR) sequencing approach, LoopSeq, to generate accurate sequencing reads that span full length transcripts using standard short read data. LoopSeq identified isoforms from control samples with 99.4% accuracy and a 0.01% per-base error rate, exceeding the accuracy reported for other long-read technologies. Applied to targeted transcriptome sequencing from colon cancers and their metastatic counterparts, LoopSeq revealed large scale isoform redistributions from benign colon mucosa to primary colon cancer and metastatic cancer and identified several previously unknown fusion isoforms. Strikingly, single nucleotide variants (SNVs) occurred dominantly in specific isoforms and some SNVs underwent isoform switching in cancer progression. The ability to use short reads to generate accurate long-read data as the raw unit of information holds promise as a widely accessible approach in transcriptome sequencing.

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.

Role of cysteine residues in thermal inactivation of fungal Cel6A cellobiohydrolases.

Numerous protein engineering studies have focused on increasing the thermostability of fungal cellulases to improve production of fuels and chemicals from lignocellulosic feedstocks. However, the engineered enzymes still undergo thermal inactivation at temperatures well below the inactivation temperatures of hyperthermophilic cellulases. In this report, we investigated the role of free cysteines in the thermal inactivation of wild-type and engineered fungal family 6 cellobiohydrolases (Cel6A). The mechanism of thermal inactivation of Cel6A is consistent with disulfide bond degradation and thiol-disulfide exchange. Circular dichroism spectroscopy revealed that a thermostable variant lacking free cysteines refolds to a native-like structure and retains activity after heat treatment over the pH range 5-9. Whereas conserved disulfide bonds are essential for retaining activity after heat treatment, free cysteines contribute to irreversible thermal inactivation in engineered thermostable Cel6A as well as Cel6A from Hypocrea jecorina and Humicola insolens.

Engineered thermostable fungal Cel6A and Cel7A cellobiohydrolases hydrolyze cellulose efficiently at elevated temperatures.

Thermostability is an important feature in industrial enzymes: it increases biocatalyst lifetime and enables reactions at higher temperatures, where faster rates and other advantages ultimately reduce the cost of biocatalysis. Here we report the thermostabilization of a chimeric fungal family 6 cellobiohydrolase (HJPlus) by directed evolution using random mutagenesis and recombination of beneficial mutations. Thermostable variant 3C6P has a half-life of 280 min at 75°C and a T(50) of 80.1°C, a ~15°C increase over the thermostable Cel6A from Humicola insolens (HiCel6A) and a ~20°C increase over that from Hypocrea jecorina (HjCel6A). Most of the mutations also stabilize the less-stable HjCel6A, the wild-type Cel6A closest in sequence to 3C6P. During a 60-h Avicel hydrolysis, 3C6P released 2.4 times more cellobiose equivalents at its optimum temperature (T(opt)) of 75°C than HiCel6A at its T(opt) of 60°C. The total cellobiose equivalents released by HiCel6A at 60°C after 60 h is equivalent to the total released by 3C6P at 75°C after ~6 h, a 10-fold reduction in hydrolysis time. A binary mixture of thermostable Cel6A and Cel7A hydrolyzes Avicel synergistically and released 1.8 times more cellobiose equivalents than the wild-type mixture, both mixtures assessed at their respective T(opt). Crystal structures of HJPlus and 3C6P, determined at 1.5 and 1.2 Å resolution, indicate that the stabilization comes from improved hydrophobic interactions and restricted loop conformations by introduced proline residues.

A family of thermostable fungal cellulases created by structure-guided recombination.

SCHEMA structure-guided recombination of 3 fungal class II cellobiohydrolases (CBH II cellulases) has yielded a collection of highly thermostable CBH II chimeras. Twenty-three of 48 genes sampled from the 6,561 possible chimeric sequences were secreted by the Saccharomyces cerevisiae heterologous host in catalytically active form. Five of these chimeras have half-lives of thermal inactivation at 63 degrees C that are greater than the most stable parent, CBH II enzyme from the thermophilic fungus Humicola insolens, which suggests that this chimera collection contains hundreds of highly stable cellulases. Twenty-five new sequences were designed based on mathematical modeling of the thermostabilities for the first set of chimeras. Ten of these sequences were expressed in active form; all 10 retained more activity than H. insolens CBH II after incubation at 63 degrees C. The total of 15 validated thermostable CBH II enzymes have high sequence diversity, differing from their closest natural homologs at up to 63 amino acid positions. Selected purified thermostable chimeras hydrolyzed phosphoric acid swollen cellulose at temperatures 7 to 15 degrees C higher than the parent enzymes. These chimeras also hydrolyzed as much or more cellulose than the parent CBH II enzymes in long-time cellulose hydrolysis assays and had pH/activity profiles as broad, or broader than, the parent enzymes. Generating this group of diverse, thermostable fungal CBH II chimeras is the first step in building an inventory of stable cellulases from which optimized enzyme mixtures for biomass conversion can be formulated.

Integrated regulation of Toll-like receptor responses by Notch and interferon-gamma pathways.

Toll-like receptor (TLR) responses are regulated to avoid toxicity and achieve coordinated responses appropriate for the cell environment. We found that Notch and TLR pathways cooperated to activate canonical Notch target genes, including transcriptional repressors Hes1 and Hey1, and to increase production of canonical TLR-induced cytokines TNF, IL-6, and IL-12. Cooperation by these pathways to increase target gene expression was mediated by the Notch-pathway component and transcription factor RBP-J, which also contributed to lethality after endotoxin injection. TLR- and Notch-induced Hes1 and Hey1 attenuated IL-6 and IL-12 production. This Hes1- and Hey1-mediated feedback inhibitory loop was abrogated by interferon-gamma (IFN-gamma), which blocked TLR-induced activation of canonical Notch target genes by inhibiting Notch2 signaling and downstream transcription. These findings identify new immune functions for RBP-J, Hes, and Hey proteins and provide insights into mechanisms by which Notch, TLR, and IFN-gamma signals are integrated to modulate specific effector functions in macrophages.