ngs_backbone: a pipeline for read cleaning, mapping and SNP calling using Next Generation Sequence

Background: The possibilities offered by next generation sequencing (NGS) platforms are revolutionizing biotechnological laboratories. Moreover, the combination of NGS sequencing and affordable high-throughput genotyping technologies is facilitating the rapid discovery and use of SNPs in non-model species. However, this abundance of sequences and polymorphisms creates new software needs. To fulfill these needs, we have developed a powerful, yet easy-to-use application. Results: The ngs_backbone software is a parallel pipeline capable of analyzing Sanger, 454, Illumina and SOLiD (Sequencing by Oligonucleotide Ligation and Detection) sequence reads. Its main supported analyses are: read cleaning, transcriptome assembly and annotation, read mapping and single nucleotide polymorphism (SNP) calling and selection. In order to build a truly useful tool, the software development was paired with a laboratory experiment. All public tomato Sanger EST reads plus 14.2 million Illumina reads were employed to test the tool and predict polymorphism in tomato. The cleaned reads were mapped to the SGN tomato transcriptome obtaining a coverage of 4.2 for Sanger and 8.5 for Illumina. 23,360 single nucleotide variations (SNVs) were predicted. A total of 76 SNVs were experimentally validated, and 85% were found to be real. Conclusions: ngs_backbone is a new software package capable of analyzing sequences produced by NGS technologies and predicting SNVs with great accuracy. In our tomato example, we created a highly polymorphic collection of SNVs that will be a useful resource for tomato researchers and breeders. The software developed along with its documentation is freely available under the AGPL license and can be downloaded from http://bioinf. comav.upv.es/ngs_backbone/ or http://github.com/JoseBlanca/franklin.Blanca Postigo, JM.; Pascual Bañuls, L.; Ziarsolo Areitioaurtena, P.; Nuez Viñals, F.; Cañizares Sales, J. (2011). Ngs_backbone: a pipeline for read cleaning, mapping and SNP calling using Next Generation Sequence. BMC Genomics. 12:1-8. doi:10.1186/1471-2164-12-285S1812Metzker ML: Sequencing technologies - the next generation. Nature Reviews Genetics. 2010, 11 (1): 31-46. 10.1038/nrg2626.454 sequencing. [ http://www.454.com/ ]Illumina Inc. [ http://www.illumina.com/ ]Flicek P, Birney E: Sense from sequence reads: methods for alignment and assembly (vol 6, pg S6, 2009). Nature Methods. 2010, 7 (6): 479-479.Chevreux B, Pfisterer T, Drescher B, Driesel AJ, Muller WEG, Wetter T, Suhai S: Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs. Genome Research. 2004, 14 (6): 1147-1159. 10.1101/gr.1917404.Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25 (14): 1754-1760. 10.1093/bioinformatics/btp324.Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biology. 2009, 10 (3):Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Genome Project Data P: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.1000 Genomes. A deep Catalog of Human Genetic Variation. [ http://1000genomes.org/wiki/doku.php?id=1000_genomes:analysis:vcf4.0 ]The seqanswers internet forum. [ http://seqanswers.com/ ]Blankenberg D, Taylor J, Schenck I, He JB, Zhang Y, Ghent M, Veeraraghavan N, Albert I, Miller W, Makova KD, Ross CH, Nekrutenko A: A framework for collaborative analysis of ENCODE data: Making large-scale analyses biologist-friendly. Genome Research. 2007, 17 (6): 960-964. 10.1101/gr.5578007.CloVR Automated Sequence Analysis from Your Desktop. [ http://clovr.org/ ]Papanicolaou A, Stierli R, Ffrench-Constant RH, Heckel DG: Next generation transcriptomes for next generation genomes using est2assembly. Bmc Bioinformatics. 2009, 10:Applied Biosystems by life technologies. [ http://www.appliedbiosystems.com/absite/us/en/home/applications-technologies/solid-next-generation-sequencing.html ]Wall PK, Leebens-Mack J, Chanderbali AS, Barakat A, Wolcott E, Liang HY, Landherr L, Tomsho LP, Hu Y, Carlson JE, Ma H, Schuster SC, Soltis DE, Soltis PS, Altman N, dePamphilis CW: Comparison of next generation sequencing technologies for transcriptome characterization. Bmc Genomics. 2009, 10:Murchison EP, Tovar C, Hsu A, Bender HS, Kheradpour P, Rebbeck CA, Obendorf D, Conlan C, Bahlo M, Blizzard CA, Pyecroft S, Kreiss A, Kellis M, Stark A, Harkins TT, Marshall Graves JA, Woods GM, Hanon GJ, Papenfuss AT: The Tasmanian Devil Transcriptome Reveals Schwann Cell Origins of a Clonally Transmissible Cancer. Science. 2010, 327 (5961): 84-87. 10.1126/science.1180616.Parchman TL, Geist KS, Grahnen JA, Benkman CW, Buerkle CA: Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery. Bmc Genomics. 2010, 11:Babik W, Stuglik M, Qi W, Kuenzli M, Kuduk K, Koteja P, Radwan J: Heart transcriptome of the bank vole (Myodes glareolus): towards understanding the evolutionary variation in metabolic rate. BMC Genomics. 2010, 11: 390-10.1186/1471-2164-11-390.Miller JC, Tanksley SD: RFLP analysis of phylogenetic-relationships and genetic-variation in the genus Lycopersicon. Theoretical and Applied Genetics. 1990, 80 (4): 437-448.Williams CE, Stclair DA: Phenetic relationships and levels of variability detected by restriction-fragment-length-polymorphism and random amplified polymorphic DNA analysis of cultivated and wild accessions of Lycopersicon-esculentum. Genome. 1993, 36 (3): 619-630. 10.1139/g93-083.Rick CM: Tomato, Lycopersicon esculentum (Solanaceae). Evolution of crop plants. Edited by: Simmonds NW. 1976, London: Longman Group, 268-273.Labate JA, Baldo AM: Tomato SNP discovery by EST mining and resequencing. Molecular Breeding. 2005, 16 (4): 343-349. 10.1007/s11032-005-1911-5.Yano K, Watanabe M, Yamamoto N, Maeda F, Tsugane T, Shibata D: Expressed sequence tags (EST) database of a miniature tomato cultivar, Micro-Tom. Plant and Cell Physiology. 2005, 46: S139-S139.Jimenez-Gomez JM, Maloof JN: Sequence diversity in three tomato species: SNPs, markers, and molecular evolution. Bmc Plant Biology. 2009, 9:Yang WC, Bai XD, Kabelka E, Eaton C, Kamoun S, van der Knaap E, Francis D: Discovery of single nucleotide polymorphisms in Lycopersicon esculentum by computer aided analysis of expressed sequence tags. Molecular Breeding. 2004, 14 (1): 21-34.Van Deynze A, Stoffel K, Buell CR, Kozik A, Liu J, van der Knaap E, Francis D: Diversity in conserved genes in tomato. Bmc Genomics. 2007, 8:Sim SC, Robbins MD, Chilcott C, Zhu T, Francis DM: Oligonucleotide array discovery of polymorphisms in cultivated tomato (Solanum lycopersicum L.) reveals patterns of SNP variation associated with breeding. Bmc Genomics. 2009, 10:Bioinformatics at COMAV. [ http://bioinf.comav.upv.es/ngs_backbone/index.html ]Broad institute. [ http://www.broadinstitute.org/igv ]Bioinformatics at COMAV. [ http://bioinf.comav.upv.es/ngs_backbone/install.html ]Github social coding. [ http://github.com/JoseBlanca/franklin ]Chou HH, Holmes MH: DNA sequence quality trimming and vector removal. Bioinformatics. 2001, 17 (12): 1093-1104. 10.1093/bioinformatics/17.12.1093.Picard. [ http://picard.sourceforge.net/index.shtml ]McKenna A, Hanna M, Banks E, Sivachenko A, Citulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA: The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Research. 2010, 20: 1297-1303. 10.1101/gr.107524.110.Sol Genomics Network. [ ftp://ftp.solgenomics.net/ ]NCBI Genbank. [ http://www.ncbi.nlm.nih.gov/genbank/ ]Gundry CN, Vandersteen JG, Reed GH, Pryor RJ, Chen J, Wittwer CT: Amplicon melting analysis with labeled primers: A closed-tube method for differentiating homozygotes and heterozygotes. Clinical Chemistry. 2003, 49 (3): 396-406. 10.1373/49.3.396

Total Elbow Arthroplasty

Total elbow arthroplasty has continued to evolve over time. Elbow implants may be linked or unlinked. Unlinked implants are attractive for patients with relatively well preserved bone stock and ligaments, but many favor linked implants, since they prevent instability and allow replacement for a wider spectrum of indications. Inflammatory arthropathies such as rheumatoid arthritis represent the classic indication for elbow arthroplasty. Indications have been expanded to include posttraumatic osteoarthritis, acute distal humerus fractures, distal humerus nonunions and reconstruction after tumor resection. Elbow arthroplasty is very successful in terms of pain relief, motion and function. However, its complication rate remains higher than arthroplasty of other joints. The overall success rate is best for patients with inflammatory arthritis and elderly patients with acute distal humerus fractures, worse for patients with posttraumatic osteoarthritis. The most common complications of elbow arthroplasty include infection, loosening, wear, triceps weakness and ulnar neuropathy. When revision surgery becomes necessary, bone augmentation techniques provide a reasonable outcome

M2 Macrophages Activate WNT Signaling Pathway in Epithelial Cells: Relevance in Ulcerative Colitis

Macrophages, which exhibit great plasticity, are important components of the inflamed tissue and constitute an essential element of regenerative responses. Epithelial Wnt signalling is involved in mechanisms of proliferation and differentiation and expression of Wnt ligands by macrophages has been reported. We aim to determine whether the macrophage phenotype determines the expression of Wnt ligands, the influence of the macrophage phenotype in epithelial activation of Wnt signalling and the relevance of this pathway in ulcerative colitis. Human monocyte-derived macrophages and U937-derived macrophages were polarized towards M1 or M2 phenotypes and the expression of Wnt1 and Wnt3a was analyzed by qPCR. The effects of macrophages and the role of Wnt1 were analyzed on the expression of β-catenin, Tcf-4, c-Myc and markers of cell differentiation in a co-culture system with Caco-2 cells. Immunohistochemical staining of CD68, CD206, CD86, Wnt1, β-catenin and c-Myc were evaluated in the damaged and non-damaged mucosa of patients with UC. We also determined the mRNA expression of Lgr5 and c-Myc by qPCR and protein levels of β-catenin by western blot. Results show that M2, and no M1, activated the Wnt signaling pathway in co-culture epithelial cells through Wnt1 which impaired enterocyte differentiation. A significant increase in the number of CD206+ macrophages was observed in the damaged mucosa of chronic vs newly diagnosed patients. CD206 immunostaining co-localized with Wnt1 in the mucosa and these cells were associated with activation of canonical Wnt signalling pathway in epithelial cells and diminution of alkaline phosphatase activity. Our results show that M2 macrophages, and not M1, activate Wnt signalling pathways and decrease enterocyte differentiation in co-cultured epithelial cells. In the mucosa of UC patients, M2 macrophages increase with chronicity and are associated with activation of epithelial Wnt signalling and diminution in enterocyte differentiation

The Chemotherapeutic Drug 5-Fluorouracil Promotes PKR-Mediated Apoptosis in a p53- Independent Manner in Colon and Breast Cancer Cells

The chemotherapeutic drug 5-FU is widely used in the treatment of a range of cancers, but resistance to the drug remains a major clinical problem. Since defects in the mediators of apoptosis may account for chemo-resistance, the identification of new targets involved in 5-FU-induced apoptosis is of main clinical interest. We have identified the ds-RNA-dependent protein kinase (PKR) as a key molecular target of 5-FU involved in apoptosis induction in human colon and breast cancer cell lines. PKR distribution and activation, apoptosis induction and cytotoxic effects were analyzed during 5-FU and 5-FU/IFNα treatment in several colon and breast cancer cell lines with different p53 status. PKR protein was activated by 5-FU treatment in a p53-independent manner, inducing phosphorylation of the protein synthesis translation initiation factor eIF-2α and cell death by apoptosis. Furthermore, PKR interference promoted a decreased response to 5-FU treatment and those cells were not affected by the synergistic antitumor activity of 5-FU/IFNα combination. These results, taken together, provide evidence that PKR is a key molecular target of 5-FU with potential relevance in the clinical use of this drug

Evidence of human occupation in Mexico around the Last Glacial Maximum.

The initial colonization of the Americas remains a highly debated topic1, and the exact timing of the first arrivals is unknown. The earliest archaeological record of Mexico-which holds a key geographical position in the Americas-is poorly known and understudied. Historically, the region has remained on the periphery of research focused on the first American populations2. However, recent investigations provide reliable evidence of a human presence in the northwest region of Mexico3,4, the Chiapas Highlands5, Central Mexico6 and the Caribbean coast7-9 during the Late Pleistocene and Early Holocene epochs. Here we present results of recent excavations at Chiquihuite Cave-a high-altitude site in central-northern Mexico-that corroborate previous findings in the Americas10-17of cultural evidence that dates to the Last Glacial Maximum (26,500-19,000 years ago)18, and which push back dates for human dispersal to the region possibly as early as 33,000-31,000 years ago. The site yielded about 1,900 stone artefacts within a 3-m-deep stratified sequence, revealing a previously unknown lithic industry that underwent only minor changes over millennia. More than 50 radiocarbon and luminescence dates provide chronological control, and genetic, palaeoenvironmental and chemical data document the changing environments in which the occupants lived. Our results provide new evidence for the antiquity of humans in the Americas, illustrate the cultural diversity of the earliest dispersal groups (which predate those of the Clovis culture) and open new directions of research