@article {447961,
	title = {Genome-scale deletion screening of human long non-coding RNAs using a paired-guide RNA CRISPR{\textendash}Cas9 library},
	journal = {Nature Biotechnology},
	year = {2016},
	month = {2016},
	abstract = {\&nbsp;CRISPR{\textendash}Cas9 screens have been widely adopted to analyze coding-gene functions, but high-throughput screening of non-coding elements using this method is more challenging because indels caused by a single cut in non-coding regions are unlikely to produce a functional knockout. A high-throughput method to produce deletions of non-coding DNA is needed. We report a high-throughput genomic deletion strategy to screen for functional long non-coding RNAs (lncRNAs) that is based on a lentiviral paired-guide RNA (pgRNA) library. Applying our screening method, we identified 51 lncRNAs that can positively or negatively regulate human cancer cell growth. We validated 9 of 51 lncRNA hits using CRISPR{\textendash}Cas9-mediated genomic deletion, functional rescue, CRISPR activation or inhibition and gene-expression profiling. Our high-throughput pgRNA genome deletion method will enable rapid identification of functional mammalian non-coding elements.\&nbsp;},
	url = {http://www.nature.com/nbt/journal/vaop/ncurrent/full/nbt.3715.html},
	author = {Shiyou Zhu* and Li*, Wei and Jingze Liu and Chen, Chen-Hao and Qi Liao and Ping Xu and Han Xu and Xiao, Tengfei and Zhongzheng Cao and Jingyu Peng and Pengfei Yuan and Brown, Myles and Liu, Xiaole Shirley and Wensheng Wei}
}
@article {447931,
	title = {CRISPR-DO for genome-wide CRISPR design and optimization.},
	journal = {Bioinformatics},
	year = {2016},
	month = {2016 Jul 10},
	abstract = {MOTIVATION: Despite the growing popularity in using CRISPR/Cas9 technology for genome editing and gene knockout, its performance still relies on well-designed single guide RNAs (sgRNA). In this study, we propose a web application for the Design and Optimization (CRISPR-DO) of guide sequences that target both coding and non-coding regions in spCas9 CRISPR system across human, mouse, zebrafish, fly and worm genomes. CRISPR-DO uses a computational sequence model to predict sgRNA efficiency, and employs a specificity scoring function to evaluate the potential of off-target effect. It also provides information on functional conservation of target sequences, as well as the overlaps with exons, putative regulatory sequences and single-nucleotide polymorphisms (SNPs). The web application has a user-friendly genome-browser interface to facilitate the selection of the best target DNA sequences for experimental design.

AVAILABILITY AND IMPLEMENTATION: CRISPR-DO is available at http://cistrome.org/crispr/ CONTACT: qiliu@tongji.edu.cn or hanxu@jimmy.harvard.edu or xsliu@jimmy.harvard.eduSupplementary information: Supplementary data are available at Bioinformatics online.},
	issn = {1367-4811},
	doi = {10.1093/bioinformatics/btw476},
	author = {Ma, Jian and K{\"o}ster, Johannes and Qin, Qian and Hu, Shengen and Li, Wei and Chen, Chenhao and Cao, Qingyi and Wang, Jinzeng and Mei, Shenglin and Liu, Qi and Han Xu and Liu, Xiaole Shirley}
}
@article {447866,
	title = {Overcoming Therapeutic Resistance in HER2-Positive Breast Cancers with CDK4/6 Inhibitors.},
	journal = {Cancer Cell},
	volume = {29},
	number = {3},
	year = {2016},
	month = {2016 Mar 14},
	pages = {255-69},
	abstract = {Using transgenic mouse models, cell line-based functional studies, and clinical specimens, we show that cyclin D1/CDK4 mediate resistance to targeted therapy for HER2-positive breast cancer. This is overcome using CDK4/6 inhibitors. Inhibition of CDK4/6 not only suppresses Rb phosphorylation, but also reduces TSC2 phosphorylation and thus partially attenuates mTORC1 activity. This relieves feedback inhibition of upstream EGFR family kinases, resensitizing tumors to EGFR/HER2 blockade. Consequently, dual inhibition of EGFR/HER2 and CDK4/6 invokes a more potent suppression of TSC2 phosphorylation and hence mTORC1/S6K/S6RP activity. The suppression of both Rb and S6RP enhances G1 arrest and a phenotype resembling cellular senescence. In vivo, CDK4/6 inhibitors sensitize patient-derived xenograft tumors to HER2-targeted therapies and delay tumor recurrence in a transgenic model of HER2-positive breast cancer.},
	keywords = {Animals, Breast Neoplasms, Cell Line, Tumor, Cyclin-Dependent Kinase 4, Cyclin-Dependent Kinase 6, Disease Models, Animal, Drug Resistance, Neoplasm, Female, Humans, Mice, Mice, Nude, Mice, Transgenic, Multiprotein Complexes, Neoplasm Recurrence, Local, Phosphorylation, Protein Kinase Inhibitors, Receptor, Epidermal Growth Factor, Receptor, ErbB-2, TOR Serine-Threonine Kinases, Tumor Suppressor Proteins},
	issn = {1878-3686},
	doi = {10.1016/j.ccell.2016.02.006},
	author = {Goel, Shom and Wang, Qi and Watt, April C and Tolaney, Sara M and Dillon, Deborah A and Li, Wei and Susanne Ramm and Palmer, Adam C. and Yuzugullu, Haluk and Varadan, Vinay and Tuck, David and Harris, Lyndsay N and Wong, Kwok-Kin and X. Shirley Liu and Sicinski, Piotr and Winer, Eric P and Krop, Ian E and Zhao, Jean J}
}
@article {447896,
	title = {Accurate inference of isoforms from multiple sample RNA-Seq data.},
	journal = {BMC Genomics},
	volume = {16 Suppl 2},
	year = {2015},
	month = {2015},
	pages = {S15},
	abstract = {BACKGROUND: RNA-Seq based transcriptome assembly has become a fundamental technique for studying expressed mRNAs (i.e., transcripts or isoforms) in a cell using high-throughput sequencing technologies, and is serving as a basis to analyze the structural and quantitative differences of expressed isoforms between samples. However, the current transcriptome assembly algorithms are not specifically designed to handle large amounts of errors that are inherent in real RNA-Seq datasets, especially those involving multiple samples, making downstream differential analysis applications difficult. On the other hand, multiple sample RNA-Seq datasets may provide more information than single sample datasets that can be utilized to improve the performance of transcriptome assembly and abundance estimation, but such information remains overlooked by the existing assembly tools.

RESULTS: We formulate a computational framework of transcriptome assembly that is capable of handling noisy RNA-Seq reads and multiple sample RNA-Seq datasets efficiently. We show that finding an optimal solution under this framework is an NP-hard problem. Instead, we develop an efficient heuristic algorithm, called Iterative Shortest Path (ISP), based on linear programming (LP) and integer linear programming (ILP). Our preliminary experimental results on both simulated and real datasets and comparison with the existing assembly tools demonstrate that (i) the ISP algorithm is able to assemble transcriptomes with a greatly increased precision while keeping the same level of sensitivity, especially when many samples are involved, and (ii) its assembly results help improve downstream differential analysis. The source code of ISP is freely available at http://alumni.cs.ucr.edu/~liw/isp.html.},
	keywords = {Algorithms, Alternative Splicing, Animals, Computational Biology, Computer Simulation, Gene Expression Profiling, Humans, Internet, Models, Genetic, Protein Isoforms, Reproducibility of Results, Sequence Analysis, RNA, Software, Transcriptome},
	issn = {1471-2164},
	doi = {10.1186/1471-2164-16-S2-S15},
	author = {Tasnim, Masruba and Ma, Shining and Yang, Ei-Wen and Jiang, Tao and Li, Wei}
}
@article {447881,
	title = {Network analysis of gene essentiality in functional genomics experiments.},
	journal = {Genome Biol},
	volume = {16},
	year = {2015},
	month = {2015},
	pages = {239},
	abstract = {Many genomic techniques have been developed to study gene essentiality genome-wide, such as CRISPR and shRNA screens. Our analyses of public CRISPR screens suggest protein interaction networks, when integrated with gene expression or histone marks, are highly predictive of gene essentiality. Meanwhile, the quality of CRISPR and shRNA screen results can be significantly enhanced through network neighbor information. We also found network neighbor information to be very informative on prioritizing ChIP-seq target genes and survival indicator genes from tumor profiling. Thus, our study provides a general method for gene essentiality analysis in functional genomic experiments ( http://nest.dfci.harvard.edu ).},
	keywords = {Cell Line, CRISPR-Cas Systems, Gene Expression Profiling, Gene Regulatory Networks, Genes, Essential, Genomics, Humans, Neoplasms, Software, Survival Analysis},
	issn = {1474-760X},
	doi = {10.1186/s13059-015-0808-9},
	author = {Jiang, Peng and Wang, Hongfang and Li, Wei and Zang, Chongzhi and Li, Bo and Wong, Yinling J and Meyer, Cliff and Liu, Jun S and Aster, Jon C and X. Shirley Liu}
}
@article {447886,
	title = {RAS/MAPK Activation Drives Resistance to Smo Inhibition, Metastasis, and Tumor Evolution in Shh Pathway-Dependent Tumors.},
	journal = {Cancer Res},
	volume = {75},
	number = {17},
	year = {2015},
	month = {2015 Sep 1},
	pages = {3623-35},
	abstract = {Aberrant Shh signaling promotes tumor growth in diverse cancers. The importance of Shh signaling is particularly evident in medulloblastoma and basal cell carcinoma (BCC), where inhibitors targeting the Shh pathway component Smoothened (Smo) show great therapeutic promise. However, the emergence of drug resistance limits long-term efficacy, and the mechanisms of resistance remain poorly understood. Using new medulloblastoma models, we identify two distinct paradigms of resistance to Smo inhibition. Sufu mutations lead to maintenance of the Shh pathway in the presence of Smo inhibitors. Alternatively activation of the RAS-MAPK pathway circumvents Shh pathway dependency, drives tumor growth, and enhances metastatic behavior. Strikingly, in BCC patients treated with Smo inhibitor, squamous cell cancers with RAS/MAPK activation emerged from the antecedent BCC tumors. Together, these findings reveal a critical role of the RAS-MAPK pathway in drug resistance and tumor evolution of Shh pathway-dependent tumors.},
	keywords = {Anilides, Animals, Carcinoma, Basal Cell, Cell proliferation, Drug Resistance, Neoplasm, Gene Expression Regulation, Neoplastic, Hedgehog Proteins, Humans, Medulloblastoma, Mice, Mitogen-Activated Protein Kinase Kinases, Pyridines, ras Proteins, Receptors, G-Protein-Coupled, Signal Transduction, Xenograft Model Antitumor Assays},
	issn = {1538-7445},
	doi = {10.1158/0008-5472.CAN-14-2999-T},
	author = {Zhao, Xuesong and Ponomaryov, Tatyana and Ornell, Kimberly J and Zhou, Pengcheng and Dabral, Sukriti K and Pak, Ekaterina and Li, Wei and Atwood, Scott X and Whitson, Ramon J and Chang, Anne Lynn S and Li, Jiang and Oro, Anthony E and Chan, Jennifer A and Kelleher, Joseph F and Segal, Rosalind A}
}
@article {447891,
	title = {Sequence determinants of improved CRISPR sgRNA design.},
	journal = {Genome Res},
	volume = {25},
	number = {8},
	year = {2015},
	month = {2015 Aug},
	pages = {1147-57},
	abstract = {The CRISPR/Cas9 system has revolutionized mammalian somatic cell genetics. Genome-wide functional screens using CRISPR/Cas9-mediated knockout or dCas9 fusion-mediated inhibition/activation (CRISPRi/a) are powerful techniques for discovering phenotype-associated gene function. We systematically assessed the DNA sequence features that contribute to single guide RNA (sgRNA) efficiency in CRISPR-based screens. Leveraging the information from multiple designs, we derived a new sequence model for predicting sgRNA efficiency in CRISPR/Cas9 knockout experiments. Our model confirmed known features and suggested new features including a preference for cytosine at the cleavage site. The model was experimentally validated for sgRNA-mediated mutation rate and protein knockout efficiency. Tested on independent data sets, the model achieved significant results in both positive and negative selection conditions and outperformed existing models. We also found that the sequence preference for CRISPRi/a is substantially different from that for CRISPR/Cas9 knockout and propose a new model for predicting sgRNA efficiency in CRISPRi/a experiments. These results facilitate the genome-wide design of improved sgRNA for both knockout and CRISPRi/a studies.},
	keywords = {Clustered Regularly Interspaced Short Palindromic Repeats, Computational Biology, dna, Gene Knockout Techniques, HL-60 Cells, Humans, Models, Genetic, Mutation Rate, RNA, Guide},
	issn = {1549-5469},
	doi = {10.1101/gr.191452.115},
	author = {Han Xu and Xiao, Tengfei and Chen, Chen-Hao and Li, Wei and Meyer, Clifford A and Wu, Qiu and Wu, Di and Cong, Le and Zhang, Feng and Liu, Jun S and Brown, Myles and X. Shirley Liu}
}
@article {447871,
	title = {Quality control, modeling, and visualization of CRISPR screens with MAGeCK-VISPR.},
	journal = {Genome Biol},
	volume = {16},
	year = {2015},
	month = {2015},
	pages = {281},
	abstract = {High-throughput CRISPR screens have shown great promise in functional genomics. We present MAGeCK-VISPR, a comprehensive quality control (QC), analysis, and visualization workflow for CRISPR screens. MAGeCK-VISPR defines a set of QC measures to assess the quality of an experiment, and includes a maximum-likelihood algorithm to call essential genes simultaneously under multiple conditions. The algorithm uses a generalized linear model to deconvolute different effects, and employs expectation-maximization to iteratively estimate sgRNA knockout efficiency and gene essentiality. MAGeCK-VISPR also includes VISPR, a framework for the interactive visualization and exploration of QC and analysis results. MAGeCK-VISPR is freely available at http://bitbucket.org/liulab/mageck-vispr .},
	keywords = {CRISPR-Cas Systems, Gene Targeting, Genes, Essential, Quality Control, Software},
	issn = {1474-760X},
	doi = {10.1186/s13059-015-0843-6},
	url = {https://genomebiology.biomedcentral.com/articles/10.1186/s13059-015-0843-6},
	author = {Li*, Wei and K{\"o}ster*, Johannes and Han Xu and Chen, Chen-Hao and Xiao, Tengfei and Liu, Jun S and Brown, Myles and X. Shirley Liu}
}
@article {447901,
	title = {MethylPurify: tumor purity deconvolution and differential methylation detection from single tumor DNA methylomes.},
	journal = {Genome Biol},
	volume = {15},
	number = {8},
	year = {2014},
	month = {2014},
	pages = {419},
	abstract = {We propose a statistical algorithm MethylPurify that uses regions with bisulfite reads showing discordant methylation levels to infer tumor purity from tumor samples alone. MethylPurify can identify differentially methylated regions (DMRs) from individual tumor methylome samples, without genomic variation information or prior knowledge from other datasets. In simulations with mixed bisulfite reads from cancer and normal cell lines, MethylPurify correctly inferred tumor purity and identified over 96\% of the DMRs. From patient data, MethylPurify gave satisfactory DMR calls from tumor methylome samples alone, and revealed potential missed DMRs by tumor to normal comparison due to tumor heterogeneity.},
	keywords = {Adenocarcinoma, Algorithms, Breast, Breast Neoplasms, Cell Line, Computational Biology, DNA Methylation, Female, Genetic Heterogeneity, Humans, Lung Neoplasms, Models, Statistical, Sequence Analysis, DNA},
	issn = {1474-760X},
	doi = {10.1186/s13059-014-0419-x},
	author = {Zheng, Xiaoqi and Zhao, Qian and Wu, Hua-Jun and Li, Wei and Wang, Haiyun and Meyer, Clifford A and Qin, Qian Alvin and Han Xu and Zang, Chongzhi and Jiang, Peng and Li, Fuqiang and Hou, Yong and He, Jianxing and Wang, Jun and Wang, Jun and Zhang, Peng and Zhang, Yong and Liu, Xiaole Shirley}
}
@article {447906,
	title = {MAGeCK enables robust identification of essential genes from genome-scale CRISPR/Cas9 knockout screens.},
	journal = {Genome Biol},
	volume = {15},
	number = {12},
	year = {2014},
	month = {2014},
	pages = {554},
	abstract = {We propose the Model-based Analysis of Genome-wide CRISPR/Cas9 Knockout (MAGeCK) method for prioritizing single-guide RNAs, genes and pathways in genome-scale CRISPR/Cas9 knockout screens. MAGeCK demonstrates better performance compared with existing methods, identifies both positively and negatively selected genes simultaneously, and reports robust results across different experimental conditions. Using public datasets, MAGeCK identified novel essential genes and pathways, including EGFR in vemurafenib-treated A375 cells harboring a BRAF mutation. MAGeCK also detected cell type-specific essential genes, including BCR and ABL1, in KBM7 cells bearing a BCR-ABL fusion, and IGF1R in HL-60 cells, which depends on the insulin signaling pathway for proliferation.},
	keywords = {Algorithms, Caspase 9, Cell Line, Tumor, Clustered Regularly Interspaced Short Palindromic Repeats, Computational Biology, Gene Expression Regulation, Neoplastic, Gene Knockout Techniques, Genes, Essential, Genetic Predisposition to Disease, HL-60 Cells, Humans, Neoplasms},
	issn = {1474-760X},
	doi = {10.1186/s13059-014-0554-4},
	author = {Li*, Wei and Xu*, Han and Xiao, Tengfei and Cong, Le and Love, Michael I and Zhang, Feng and Irizarry, Rafael A and Liu, Jun S and Brown, Myles and X. Shirley Liu}
}
@article {447916,
	title = {Improving oligonucleotide fingerprinting of rRNA genes by implementation of polony microarray technology.},
	journal = {J Microbiol Methods},
	volume = {90},
	number = {3},
	year = {2012},
	month = {2012 Sep},
	pages = {235-40},
	abstract = {Improvements to oligonucleotide fingerprinting of rRNA genes (OFRG) were obtained by implementing polony microarray technology. OFRG is an array-based method for analyzing microbial community composition. Polonies are discrete clusters of DNA, produced by solid-phase PCR in hydrogels, and derived from individual, spatially isolated DNA molecules. The advantages of a polony-based OFRG method include higher throughput and reductions in the PCR-induced errors and compositional skew inherent in all other PCR-based community composition methods, including high-throughput sequencing of rRNA genes. Given the similarities between polony microarrays and certain aspects of sequencing methods such as the Illumina platform, we suggest that if concepts presented in this study were implemented in high-throughput sequencing protocols, a reduction of PCR-induced errors and compositional skew may be realized.},
	keywords = {Bacteria, DNA Fingerprinting, Humans, Intestines, Molecular Typing, Oligonucleotide Array Sequence Analysis, phylogeny, RNA, Bacterial, RNA, Ribosomal, Seawater, Water Microbiology},
	issn = {1872-8359},
	doi = {10.1016/j.mimet.2012.05.013},
	author = {Ruegger, Paul M and Bent, Elizabeth and Li, Wei and Jeske, Daniel R and Cui, Xinping and Braun, Jonathan and Jiang, Tao and Borneman, James}
}
@article {447911,
	title = {Transcriptome assembly and isoform expression level estimation from biased RNA-Seq reads.},
	journal = {Bioinformatics},
	volume = {28},
	number = {22},
	year = {2012},
	month = {2012 Nov 15},
	pages = {2914-21},
	abstract = {MOTIVATION: RNA-Seq uses the high-throughput sequencing technology to identify and quantify transcriptome at an unprecedented high resolution and low cost. However, RNA-Seq reads are usually not uniformly distributed and biases in RNA-Seq data post great challenges in many applications including transcriptome assembly and the expression level estimation of genes or isoforms. Much effort has been made in the literature to calibrate the expression level estimation from biased RNA-Seq data, but the effect of biases on transcriptome assembly remains largely unexplored.

RESULTS: Here, we propose a statistical framework for both transcriptome assembly and isoform expression level estimation from biased RNA-Seq data. Using a quasi-multinomial distribution model, our method is able to capture various types of RNA-Seq biases, including positional, sequencing and mappability biases. Our experimental results on simulated and real RNA-Seq datasets exhibit interesting effects of RNA-Seq biases on both transcriptome assembly and isoform expression level estimation. The advantage of our method is clearly shown in the experimental analysis by its high sensitivity and precision in transcriptome assembly and the high concordance of its estimated expression levels with quantitative reverse transcription-polymerase chain reaction data.

AVAILABILITY: CEM is freely available at http://www.cs.ucr.edu/~liw/cem.html.

CONTACT: liw@cs.ucr.edu.

SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.},
	keywords = {Algorithms, Gene Expression Profiling, High-Throughput Nucleotide Sequencing, Humans, Protein Isoforms, Sequence Analysis, RNA, Transcription, Genetic},
	issn = {1367-4811},
	doi = {10.1093/bioinformatics/bts559},
	author = {Li, Wei and Jiang, Tao}
}
@article {447926,
	title = {Inference of isoforms from short sequence reads.},
	journal = {J Comput Biol},
	volume = {18},
	number = {3},
	year = {2011},
	month = {2011 Mar},
	pages = {305-21},
	abstract = {Due to alternative splicing events in eukaryotic species, the identification of mRNA isoforms (or splicing variants) is a difficult problem. Traditional experimental methods for this purpose are time consuming and cost ineffective. The emerging RNA-Seq technology provides a possible effective method to address this problem. Although the advantages of RNA-Seq over traditional methods in transcriptome analysis have been confirmed by many studies, the inference of isoforms from millions of short sequence reads (e.g., Illumina/Solexa reads) has remained computationally challenging. In this work, we propose a method to calculate the expression levels of isoforms and infer isoforms from short RNA-Seq reads using exon-intron boundary, transcription start site (TSS) and poly-A site (PAS) information. We first formulate the relationship among exons, isoforms, and single-end reads as a convex quadratic program, and then use an efficient algorithm (called IsoInfer) to search for isoforms. IsoInfer can calculate the expression levels of isoforms accurately if all the isoforms are known and infer novel isoforms from scratch. Our experimental tests on known mouse isoforms with both simulated expression levels and reads demonstrate that IsoInfer is able to calculate the expression levels of isoforms with an accuracy comparable to the state-of-the-art statistical method and a 60 times faster speed. Moreover, our tests on both simulated and real reads show that it achieves a good precision and sensitivity in inferring isoforms when given accurate exon-intron boundary, TSS, and PAS information, especially for isoforms whose expression levels are significantly high. The software is publicly available for free at http://www.cs.ucr.edu/\~{}jianxing/IsoInfer.html.},
	keywords = {Algorithms, Alternative Splicing, Animals, Computer Simulation, Exons, Gene Expression Profiling, Genomics, High-Throughput Nucleotide Sequencing, Humans, Mice, Models, Genetic, Protein Isoforms, RNA, Time Factors},
	issn = {1557-8666},
	doi = {10.1089/cmb.2010.0243},
	author = {Feng, Jianxing and Li, Wei and Jiang, Tao}
}
@article {447921,
	title = {IsoLasso: a LASSO regression approach to RNA-Seq based transcriptome assembly.},
	journal = {J Comput Biol},
	volume = {18},
	number = {11},
	year = {2011},
	month = {2011 Nov},
	pages = {1693-707},
	abstract = {The new second generation sequencing technology revolutionizes many biology-related research fields and poses various computational biology challenges. One of them is transcriptome assembly based on RNA-Seq data, which aims at reconstructing all full-length mRNA transcripts simultaneously from millions of short reads. In this article, we consider three objectives in transcriptome assembly: the maximization of prediction accuracy, minimization of interpretation, and maximization of completeness. The first objective, the maximization of prediction accuracy, requires that the estimated expression levels based on assembled transcripts should be as close as possible to the observed ones for every expressed region of the genome. The minimization of interpretation follows the parsimony principle to seek as few transcripts in the prediction as possible. The third objective, the maximization of completeness, requires that the maximum number of mapped reads (or ?expressed segments? in gene models) be explained by (i.e., contained in) the predicted transcripts in the solution. Based on the above three objectives, we present IsoLasso, a new RNA-Seq based transcriptome assembly tool. IsoLasso is based on the well-known LASSO algorithm, a multivariate regression method designated to seek a balance between the maximization of prediction accuracy and the minimization of interpretation. By including some additional constraints in the quadratic program involved in LASSO, IsoLasso is able to make the set of assembled transcripts as complete as possible. Experiments on simulated and real RNA-Seq datasets show that IsoLasso achieves, simultaneously, higher sensitivity and precision than the state-of-art transcript assembly tools.},
	keywords = {Algorithms, Animals, Cells, Cultured, Computer Simulation, Exons, Gene Expression, Gene Expression Profiling, Humans, Mice, Models, Genetic, Multivariate Analysis, Protein Isoforms, Regression Analysis, Sequence Analysis, RNA, Software},
	issn = {1557-8666},
	doi = {10.1089/cmb.2011.0171},
	author = {Li, Wei and Feng, Jianxing and Jiang, Tao}
}
