Alper Kucukural

I-TASSER: a unified platform for automated protein structure and function prediction

The iterative threading assembly refinement (I-TASSER) server is an integrated platform for automated protein structure and function prediction based on the sequence-to-structure-to-function paradigm. Starting from an amino acid sequence, I-TASSER first generates three-dimensional (3D) atomic models from multiple threading alignments and iterative structural assembly simulations. The function of the protein is then inferred by structurally matching the 3D models with other known proteins. The output from a typical server run contains full-length secondary and tertiary structure predictions, and functional annotations on ligand-binding sites, Enzyme Commission numbers and Gene Ontology terms. An estimate of accuracy of the predictions is provided based on the confidence score of the modeling. This protocol provides new insights and guidelines for designing of online server systems for the state-of-the-art protein structure and function predictions. The server is available at http://zhanglab.ccmb.med.umich.edu/I-TASSER.

Biogenesis and function of tRNA fragments during sperm maturation and fertilization in mammals

Offspring affected by sperm small RNAs Paternal dietary conditions in mammals influence the metabolic phenotypes of offspring. Although prior work suggests the involvement of epigenetic pathways, the mechanisms remains unclear. Two studies now show that altered paternal diet affects the level of small RNAs in mouse sperm. Chen et al. injected sperm transfer RNA (tRNA) fragments from males that had been kept on a high-fat diet into normal oocytes. The progeny displayed metabolic disorders and concomitant alteration of genes in metabolic pathways. Sharma et al. observed the biogenesis and function of small tRNA-derived fragments during sperm maturation. Further understanding of the mechanisms by which progeny are affected by parental exposure may affect human diseases such as diet-induced metabolic disorders. Science, this issue p. 397, p. 391 Abundant transfer RNA fragments in maturing mammalian sperm regulate expression of endogenous retroelements in the embryo. Several recent studies link parental environments to phenotypes in subsequent generations. In this work, we investigate the mechanism by which paternal diet affects offspring metabolism. Protein restriction in mice affects small RNA (sRNA) levels in mature sperm, with decreased let-7 levels and increased amounts of 5′ fragments of glycine transfer RNAs (tRNAs). In testicular sperm, tRNA fragments are scarce but increase in abundance as sperm mature in the epididymis. Epididymosomes (vesicles that fuse with sperm during epididymal transit) carry RNA payloads matching those of mature sperm and can deliver RNAs to immature sperm in vitro. Functionally, tRNA-glycine-GCC fragments repress genes associated with the endogenous retroelement MERVL, in both embryonic stem cells and embryos. Our results shed light on sRNA biogenesis and its dietary regulation during posttesticular sperm maturation, and they also link tRNA fragments to regulation of endogenous retroelements active in the preimplantation embryo.

Identification of Neuronal RNA Targets of TDP-43-containing Ribonucleoprotein Complexes

TAR DNA-binding protein 43 (TDP-43) is associated with a spectrum of neurodegenerative diseases. Although TDP-43 resembles heterogeneous nuclear ribonucleoproteins, its RNA targets and physiological protein partners remain unknown. Here we identify RNA targets of TDP-43 from cortical neurons by RNA immunoprecipitation followed by deep sequencing (RIP-seq). The canonical TDP-43 binding site (TG)n is 55.1-fold enriched, and moreover, a variant with adenine in the middle, (TG)nTA(TG)m, is highly abundant among reads in our TDP-43 RIP-seq library. TDP-43 RNA targets can be divided into three different groups: those primarily binding in introns, in exons, and across both introns and exons. TDP-43 RNA targets are particularly enriched for Gene Ontology terms related to synaptic function, RNA metabolism, and neuronal development. Furthermore, TDP-43 binds to a number of RNAs encoding for proteins implicated in neurodegeneration, including TDP-43 itself, FUS/TLS, progranulin, Tau, and ataxin 1 and -2. We also identify 25 proteins that co-purify with TDP-43 from rodent brain nuclear extracts. Prominent among them are nuclear proteins involved in pre-mRNA splicing and RNA stability and transport. Also notable are two neuron-enriched proteins, methyl CpG-binding protein 2 and polypyrimidine tract-binding protein 2 (PTBP2). A PTBP2 consensus RNA binding motif is enriched in the TDP-43 RIP-seq library, suggesting that PTBP2 may co-regulate TDP-43 RNA targets. This work thus reveals the protein and RNA components of the TDP-43-containing ribonucleoprotein complexes and provides a framework for understanding how dysregulation of TDP-43 in RNA metabolism contributes to neurodegeneration.

Novel Observations From Next-Generation RNA Sequencing of Highly Purified Human Adult and Fetal Islet Cell Subsets.

Understanding distinct gene expression patterns of normal adult and developing fetal human pancreatic α- and β-cells is crucial for developing stem cell therapies, islet regeneration strategies, and therapies designed to increase β-cell function in patients with diabetes (type 1 or 2). Toward that end, we have developed methods to highly purify α-, β-, and δ-cells from human fetal and adult pancreata by intracellular staining for the cell-specific hormone content, sorting the subpopulations by flow cytometry, and, using next-generation RNA sequencing, we report the detailed transcriptomes of fetal and adult α- and β-cells. We observed that human islet composition was not influenced by age, sex, or BMI, and transcripts for inflammatory gene products were noted in fetal β-cells. In addition, within highly purified adult glucagon-expressing α-cells, we observed surprisingly high insulin mRNA expression, but not insulin protein expression. This transcriptome analysis from highly purified islet α- and β-cell subsets from fetal and adult pancreata offers clear implications for strategies that seek to increase insulin expression in type 1 and type 2 diabetes.

Simultaneous generation of many RNA-seq libraries in a single reaction

Although RNA-seq is a powerful tool, the considerable time and cost associated with library construction has limited its utilization for various applications. RNAtag-Seq, an approach to generate multiple RNA-seq libraries in a single reaction, lowers time and cost per sample, and it produces data on prokaryotic and eukaryotic samples that are comparable to those generated by traditional strand-specific RNA-seq approaches.

Ebola Virus Glycoprotein with Increased Infectivity Dominated the 2013–2016 Epidemic

Summary The magnitude of the 2013–2016 Ebola virus disease (EVD) epidemic enabled an unprecedented number of viral mutations to occur over successive human-to-human transmission events, increasing the probability that adaptation to the human host occurred during the outbreak. We investigated one nonsynonymous mutation, Ebola virus (EBOV) glycoprotein (GP) mutant A82V, for its effect on viral infectivity. This mutation, located at the NPC1-binding site on EBOV GP, occurred early in the 2013–2016 outbreak and rose to high frequency. We found that GP-A82V had heightened ability to infect primate cells, including human dendritic cells. The increased infectivity was restricted to cells that have primate-specific NPC1 sequences at the EBOV interface, suggesting that this mutation was indeed an adaptation to the human host. GP-A82V was associated with increased mortality, consistent with the hypothesis that the heightened intrinsic infectivity of GP-A82V contributed to disease severity during the EVD epidemic.

ASPeak: an abundance sensitive peak detection algorithm for RIP-Seq

SUMMARY Unlike DNA, RNA abundances can vary over several orders of magnitude. Thus, identification of RNA-protein binding sites from high-throughput sequencing data presents unique challenges. Although peak identification in ChIP-Seq data has been extensively explored, there are few bioinformatics tools tailored for peak calling on analogous datasets for RNA-binding proteins. Here we describe ASPeak (abundance sensitive peak detection algorithm), an implementation of an algorithm that we previously applied to detect peaks in exon junction complex RNA immunoprecipitation in tandem experiments. Our peak detection algorithm yields stringent and robust target sets enabling sensitive motif finding and downstream functional analyses. AVAILABILITY ASPeak is implemented in Perl as a complete pipeline that takes bedGraph files as input. ASPeak implementation is freely available at https://sourceforge.net/projects/as-peak under the GNU General Public License. ASPeak can be run on a personal computer, yet is designed to be easily parallelizable. ASPeak can also run on high performance computing clusters providing efficient speedup. The documentation and user manual can be obtained from http://master.dl.sourceforge.net/project/as-peak/manual.pdf.

Generation and Analysis of Expressed Sequence Tags from Olea europaea L.

Olive (Olea europaea L.) is an important source of edible oil which was originated in Near-East region. In this study, two cDNA libraries were constructed from young olive leaves and immature olive fruits for generation of ESTs to discover the novel genes and search the function of unknown genes of olive. The randomly selected 3840 colonies were sequenced for EST collection from both libraries. Readable 2228 sequences for olive leaf and 1506 sequences for olive fruit were assembled into 205 and 69 contigs, respectively, whereas 2478 were singletons. Putative functions of all 2752 differentially expressed unique sequences were designated by gene homology based on BLAST and annotated using BLAST2GO. While 1339 ESTs show no homology to the database, 2024 ESTs have homology (under 80%) with hypothetical proteins, putative proteins, expressed proteins, and unknown proteins in NCBI-GenBank. 635 ESTs unique genes sequence have been identified by over 80% homology to known function in other species which were not previously described in Olea family. Only 3.1% of total ESTs was shown similarity with olive database existing in NCBI. This generated ESTs data and consensus sequences were submitted to NCBI as valuable source for functional genome studies of olive.

GUIDEseq: a bioconductor package to analyze GUIDE-Seq datasets for CRISPR-Cas nucleases

BackgroundGenome editing technologies developed around the CRISPR-Cas9 nuclease system have facilitated the investigation of a broad range of biological questions. These nucleases also hold tremendous promise for treating a variety of genetic disorders. In the context of their therapeutic application, it is important to identify the spectrum of genomic sequences that are cleaved by a candidate nuclease when programmed with a particular guide RNA, as well as the cleavage efficiency of these sites. Powerful new experimental approaches, such as GUIDE-seq, facilitate the sensitive, unbiased genome-wide detection of nuclease cleavage sites within the genome. Flexible bioinformatics analysis tools for processing GUIDE-seq data are needed.ResultsHere, we describe an open source, open development software suite, GUIDEseq, for GUIDE-seq data analysis and annotation as a Bioconductor package in R. The GUIDEseq package provides a flexible platform with more than 60 adjustable parameters for the analysis of datasets associated with custom nuclease applications. These parameters allow data analysis to be tailored to different nuclease platforms with different length and complexity in their guide and PAM recognition sequences or their DNA cleavage position. They also enable users to customize sequence aggregation criteria, and vary peak calling thresholds that can influence the number of potential off-target sites recovered. GUIDEseq also annotates potential off-target sites that overlap with genes based on genome annotation information, as these may be the most important off-target sites for further characterization. In addition, GUIDEseq enables the comparison and visualization of off-target site overlap between different datasets for a rapid comparison of different nuclease configurations or experimental conditions. For each identified off-target, the GUIDEseq package outputs mapped GUIDE-Seq read count as well as cleavage score from a user specified off-target cleavage score prediction algorithm permitting the identification of genomic sequences with unexpected cleavage activity.ConclusionThe GUIDEseq package enables analysis of GUIDE-data from various nuclease platforms for any species with a defined genomic sequence. This software package has been used successfully to analyze several GUIDE-seq datasets. The software, source code and documentation are freely available at http://www.bioconductor.org/packages/release/bioc/html/GUIDEseq.html.

Comparison of RNA isolation and associated methods for extracellular RNA detection by high-throughput quantitative polymerase chain reaction

MicroRNAs (miRNAs) are small noncoding RNA molecules that function in RNA silencing and posttranscriptional regulation of gene expression. miRNAs in biofluids are being used for clinical diagnosis as well as disease prediction. Efficient and reproducible isolation methods are crucial for extracellular RNA detection. To determine the best methodologies for miRNA detection from plasma, the performance of four RNA extraction kits, including an in-house kit, were determined with miScript miRNA assay technology; all were measured using a high-throughput quantitative polymerase chain reaction (qPCR) platform (BioMark System) with 90 human miRNA assays. In addition, the performances of complementary DNA (cDNA) and preamplification kits for TaqMan miRNA assays and miScript miRNA assays were compared using the same 90 miRNAs on the BioMark System. There were significant quantification cycle (Cq) value differences for the detection of miRNA targets between isolation kits. cDNA, preamplification, and qPCR performances were also varied. In summary, this study demonstrates differences among RNA isolation methods as measured by reverse transcription (RT)-qPCR. Importantly, differences were also noted in cDNA and preamplification performance using TaqMan and miScript. The in-house kit performed better than the other three kits. These findings demonstrate significant variability between isolation and detection methods for low-abundant miRNA detection from biofluids.