ROKET: Associating Somatic Mutation with Clinical Outcomes through Kernel Regression and Optimal Transport

Somatic mutations in cancer patients are inherently sparse and potentially high dimensional. Cancer patients may share the same set of deregulated biological processes perturbed by different sets of somatically mutated genes. Therefore, when assessing the associations between somatic mutations and clinical outcomes, gene-by-gene analyses is often under-powered because it does not capture the complex disease mechanisms shared across cancer patients. Rather than testing genes one by one, an intuitive approach is to aggregate somatic mutation data of multiple genes to assess the joint association. The challenge is how to aggregate such information. Building on the optimal transport method, we propose a principled approach to estimate the similarity of somatic mutation profiles of multiple genes between tumor samples, while accounting for gene-gene similarity defined by gene annotations or empirical mutational patterns. Using such similarities, we can assess the associations between somatic mutations and clinical outcomes by kernel regression. We have applied our method to analyze somatic mutation data of 17 cancer types and identified at least three cancer types harboring associations between somatic mutations and overall survival, progression-free interval or cytolytic activity.

Multiplexed CRISPR-mediated engineering of protein secretory pathway genes in the thermotolerant methylotrophic yeast Ogataea thermomethanolica

CRISPR multiplex gRNA systems have been employed in genome engineering in various industrially relevant yeast species. The thermotolerant methylotrophic yeast Ogataea thermomethanolica TBRC 656 is an alternative host for heterologous protein production. However, the limited secretory capability of this yeast is a bottleneck for protein production. Here, we refined CRISPR-based genome engineering tools for simultaneous mutagenesis and activation of multiple protein secretory pathway genes to improve heterologous protein secretion. We demonstrated that multiplexed CRISPR-Cas9 mutation of up to four genes (SOD1, VPS1, YPT7 and YPT35) in one single cell is practicable. We also developed a multiplexed CRISPR-dCas9 system which allows simultaneous activation of multiple genes in this yeast. 27 multiplexed gRNA combinations were tested for activation of three genes (SOD1, VPS1 and YPT7), three of which were demonstrated to increase the secretion of fungal xylanase and phytase up to 29% and 41%, respectively. Altogether, our study provided a toolkit for mutagenesis and activation of multiple genes in O. thermomethanolica, which could be useful for future strain engineering to improve heterologous protein production in this yeast.

Development of New Rice (Oryza. sativa L.) Breeding Lines through Marker-Assisted Introgression and Pyramiding of Brown Planthopper, Blast, Bacterial Leaf Blight Resistance, and Aroma Genes

Brown planthopper, blast, and bacterial blight are the main biotic stressors of rice and can cause a massive loss in rice production. Aroma is an important character of rice quality. It is of far-reaching significance to breed resistant and high-quality varieties using germplasms with objective genes. In this study, the introgression and pyramiding of brown planthopper (BPH), blast, and bacterial leaf blight (BLB) resistance genes and aroma genes into elite rice maintainers and restorers were conducted through conventional cross-breeding coupled with the marker-assisted selection (MAS) breeding method. Single-plant selection was performed from F2 onwards to select desirable recombinants possessing alleles of interest with suitable phenotypes. Respective linked markers were used in each generation from intercrossing to the F7 generation for tracking the presence of targeted genes. A total of 74 improved lines (ILs) have been developed which possess a combination of 1 to 4 genes for BPH, blast, and BLB resistance and aroma. These ILs showed moderate to high resistance to multiple biotic stresses (BPH, blast and BLB) or aromatic fragrance without obvious negative effects on agronomic traits. As multiple resistance and aromatic traits have become significant objectives in rice breeding, these resistance and/or aroma gene introgressed or pyramided lines have important application prospects. Core ideas: (1) marker-assisted breeding was used to pyramid multiple genes for an elite breeding line; (2) improved lines with the introgression of 1–4 genes were developed to achieve high resistance against various biotic stresses and aroma; (3) new lines were used as donor parents to introgress multiple genes in other genetic backgrounds.

Multiple genes in cis mediate the effects of a single chromatin accessibility variant on aberrant synaptic development and function in human neurons

Despite hundreds of risk loci from genome-wide association studies of neuropsychiatric disorders, causal variants/genes remain largely unknown. Here, in NEUROG2-induced human neurons, we identified 31 risk SNPs in 26 schizophrenia (SZ) risk loci that displayed allele-specific open chromatin (ASoC) and were likely to be functional. Editing the strongest ASoC SNP rs2027349 near vacuolar protein sorting 45 homolog (VPS45) altered the expression of VPS45, lncRNA AC244033.2, and a distal gene, C1orf54, in human neurons. Notably, the global gene expression changes in neurons were enriched for SZ risk and correlated with post-mortem brain gene expression signatures of neuropsychiatric disorders. Neurons carrying the risk allele exhibited increased dendritic complexity, synaptic puncta density, and hyperactivity, which were reversed by knocking-down distinct cis-regulated genes (VPS45, AC244033.2, or C1orf54), suggesting a phenotypic contribution from all three genes. Interestingly, transcriptomic analysis of knockdown cells suggested a non-additive effects of these genes. Our study reveals a compound effect of multiple genes at a single SZ locus on synaptic development and function, providing a mechanistic link between a non-coding SZ risk variant and disease-related cellular phenotypes.

Construction of Multiple Guide RNAs in CRISPR/Cas9 Vector Using Stepwise or Simultaneous Golden Gate Cloning: Case Study for Targeting the FAD2 and FATB Multigene in Soybean

CRISPR/Cas9 is a commonly used technique in reverse-genetics research to knock out a gene of interest. However, when targeting a multigene family or multiple genes, it is necessary to construct a vector with multiple single guide RNAs (sgRNAs) that can navigate the Cas9 protein to the target site. In this protocol, the Golden Gate cloning method was used to generate multiple sgRNAs in the Cas9 vector. The vectors used were pHEE401E_UBQ_Bar and pBAtC_tRNA, which employ a one-promoter/one-sgRNA and a polycistronic-tRNA-gRNA strategy, respectively. Golden Gate cloning was performed with type IIS restriction enzymes to generate gRNA polymers for vector inserts. Four sgRNAs containing the pHEE401E_UBQ_Bar vector and four to six sgRNAs containing the pBAtC_tRNA vector were constructed. In practice, we constructed multiple sgRNAs targeting multiple genes of FAD2 and FATB in soybean using this protocol. These three vectors were transformed into soybeans using the Agrobacterium-mediated method. Using deep sequencing, we confirmed that the T0 generation transgenic soybean was edited at various indel ratios in the predicted target regions of the FAD2 and FATB multigenes. This protocol is a specific guide that allows researchers to easily follow the cloning of multiple sgRNAs into commonly used CRISPR/Cas9 vectors for plants.

Non-Coding RNAs in Response to Drought Stress

Drought stress causes changes in the morphological, physiological, biochemical and molecular characteristics of plants. The response to drought in different plants may vary from avoidance, tolerance and escape to recovery from stress. This response is genetically programmed and regulated in a very complex yet synchronized manner. The crucial genetic regulations mediated by non-coding RNAs (ncRNAs) have emerged as game-changers in modulating the plant responses to drought and other abiotic stresses. The ncRNAs interact with their targets to form potentially subtle regulatory networks that control multiple genes to determine the overall response of plants. Many long and small drought-responsive ncRNAs have been identified and characterized in different plant varieties. The miRNA-based research is better documented, while lncRNA and transposon-derived RNAs are relatively new, and their cellular role is beginning to be understood. In this review, we have compiled the information on the categorization of non-coding RNAs based on their biogenesis and function. We also discuss the available literature on the role of long and small non-coding RNAs in mitigating drought stress in plants.

Evaluating the Potential of Applying Machine Learning Tools to Metabolic Pathway Optimization

Successful engineering of a microbial host for efficient production of a target product from a given substrate can be viewed as an extensive optimization task. Such a task involves the selection of high activity enzymes as well as their gene expression regulatory control elements (i.e., promoters and ribosome binding sites). Finally, there is also the need to tune expression of multiple genes along a heterologous pathway to relieve constraints from rate-limiting step and help reduce metabolic burden on cells from unnecessary over-expression of high activity enzymes. While the aforementioned tasks could be performed through combinatorial experiments, such an approach incurs significant cost, time and effort, which is a handicap that can be relieved by application of modern machine learning tools. Such tools could attempt to predict high activity enzymes from sequence, but they are currently most usefully applied in classifying strong promoters from weaker ones as well as combinatorial tuning of expression of multiple genes. This perspective reviews the application of machine learning tools to aid metabolic pathway optimization through identifying challenges in metabolic engineering that could be overcome with the help of machine learning tools.

687 MDSC gene expression analysis in patients with cancer and the response to inhibition of Bruton's tyrosine kinase

BackgroundMyeloid-derived suppressor cells (MDSC) are an immunosuppressive immature population of myeloid cells that are elevated in cancer patients. Increased levels of MDSC has been linked to dysregulated anti-tumor responses and reduced efficacy of immune checkpoint therapies thus making them an attractive target. MDSC express Bruton's tyrosine kinase (BTK) and can be depleted using ibrutinib, an FDA-approved irreversible inhibitor of BTK. BTK inhibition leads to reduced MDSC expansion/function in murine models and significantly improved activity of anti-PD-1 antibodies. In this study, single cell RNA-seq (scRNA-seq) was used to characterize the gene expression of MDSC from different cancer types and the effect of ibrutinib on MDSC gene expression.MethodsPeripheral blood mononuclear cells were isolated from patients with melanoma (n=2), head & neck (n=1), and breast cancer (n=1). MDSC were isolated via fluorescence activated cell sorting. MDSC isolated from melanoma patients (n=2) were treated in vitro for 4h with 1 uM ibrutinib or DMSO and scRNA-seq was performed using the Chromium 10x Genomics platform. ScRNA-seq samples were analyzed using the standard integrative workflow of Seurat v3, which addresses the sample heterogeneity. Cell clusters were identified using Seurat and annotated using SingleRversion3.12. Identification of gene markers for each cell cluster and cell-cluster-specific differential expression analyses were conducted using Seurat.ResultsBaseline gene expression of MDSC from patients with breast and head & neck cancer revealed similarities among the top expressed genes (S100A8, VCAN, and LYZ). In vitro ibrutinib treatment of MDSC from patients with melanoma resulted in significant changes in gene expression within the MDSC cluster compared to DMSO treatment. GBP1(-1.72 log fold change), IL 1β(-1.27 log fold change), and CXCL8(-0.63 log fold change) were among the top downregulated genes (p<0.001) and RGS2 (0.68 log fold change) and ABHD5(0.52 log fold change) were among the top upregulated genes (p<0.001). MDSC subset (PMN-MDSC, M-MDSC, early-MDSC, and CD14+/CD15+ double positive) gene expression changes mirrored total MDSC gene changes. Ingenuity pathway analysis revealed significant downregulated pathways including TREM1 (p<0.001), nitric oxide signaling (p<0.003), and IL-6 signaling (p<0.004). Multiple genes associated with cellular movement (CXCL8, CXCL10) and activation of macrophages (CXCL10, CCL3) were downregulated (p<0.001). PCR analysis on isolated melanoma MDSC (n=2) treated in vitro with ibrutinib verified downregulation of CXCL8 (0.42 fold decrease, p<0.05) and CXCL10 (0.40 fold decrease, p<0.001).ConclusionsAnalysis via scRNA-seq revealed similar gene expression patterns for MDSC from different cancer patients. There was downregulation of multiple genes and pathways important to MDSC function and migration after BTK inhibition.Ethics ApprovalThe study obtained ethics approval. IRB# 1999C0348

The Functional Association of ACQOS/VICTR with Salt Stress Resistance in Arabidopsis thaliana Was Confirmed by CRISPR-Mediated Mutagenesis

Clustered regularly interspaced palindromic repeat (CRISPR)-mediated mutagenesis has become an important tool in plant research, enabling the characterization of genes via gene knock-out. CRISPR genome editing tools can be applied to generate multi-gene knockout lines. Typically, multiple single-stranded, single guide RNAs (gRNAs) must be expressed in an organism to target multiple genes simultaneously; however, a single gRNA can target multiple genes if the target genes share similar sequences. A gene cluster comprising ACQUIRED OSMOTOLERANCE (ACQOS; AT5G46520) and neighboring nucleotide-binding leucine-rich repeats (NLRs; AT5G46510) is associated with osmotic tolerance. To investigate the role of ACQOS and the tandemly arranged NLR in osmotic tolerance, we introduced small insertion/deletion mutations into two target genes using a single gRNA and obtained transformant plant lines with three different combinations of mutant alleles. We then tested our mutant lines for osmotic tolerance after a salt-stress acclimation period by determining the chlorophyll contents of the mutant seedlings. Our results strongly suggest that ACQOS is directly associated with salt resistance, while the neighboring NLR is not. Here, we confirmed previous findings suggesting the involvement of ACQOS in salt tolerance and demonstrated the usefulness of CRISPR-mediated mutagenesis in validating the functions of genes in a single genetic background.

Von Genen zu Chromosomen: Pflanzenzüchtung mit CRISPR-CAS

Using the CRISPR-Cas system, it has been possible to introduce different kinds of mutations in single or multiple genes for trait improvement in crops. Last year, for the first time, the CRISPR-Cas-mediated induction of different kinds of targeted heritable chromosomal rearrangements has been achieved in plants. This novel application has the potential to revolutionize plant breeding as genetic exchange and linkage drag are now becoming controllable in a targeted manner.