Feature-by-Feature – Evaluating De Novo Sequence Assembly

The whole-genome sequence assembly (WGSA) problem is among one of the most studied problems in computational biology. Despite the availability of a plethora of tools (i.e., assemblers), all claiming to have solved the WGSA problem, little has been done to systematically compare their accuracy and power. Traditional methods rely on standard metrics and read simulation: while on the one hand, metrics like N50 and number of contigs focus only on size without proportionately emphasizing the information about the correctness of the assembly, comparisons performed on simulated dataset, on the other hand, can be highly biased by the non-realistic assumptions in the underlying read generator. Recently the Feature Response Curve (FRC) method was proposed to assess the overall assembly quality and correctness: FRC transparently captures the trade-offs between contigs' quality against their sizes. Nevertheless, the relationship among the different features and their relative importance remains unknown. In particular, FRC cannot account for the correlation among the different features. We analyzed the correlation among different features in order to better describe their relationships and their importance in gauging assembly quality and correctness. In particular, using multivariate techniques like principal and independent component analysis we were able to estimate the “excess-dimensionality” of the feature space. Moreover, principal component analysis allowed us to show how poorly the acclaimed N50 metric describes the assembly quality. Applying independent component analysis we identified a subset of features that better describe the assemblers performances. We demonstrated that by focusing on a reduced set of highly informative features we can use the FRC curve to better describe and compare the performances of different assemblers. Moreover, as a by-product of our analysis, we discovered how often evaluation based on simulated data, obtained with state of the art simulators, lead to not-so-realistic results

TESLA Technical Design Report Part III: Physics at an e+e- Linear Collider

The TESLA Technical Design Report Part III: Physics at an e+e- Linear ColliderComment: 192 pages, 131 figures. Some figures have reduced quality. Full quality figures can be obtained from http://tesla.desy.de/tdr. Editors - R.-D. Heuer, D.J. Miller, F. Richard, P.M. Zerwa

Search for rare quark-annihilation decays, B --> Ds(*) Phi

We report on searches for B- --> Ds- Phi and B- --> Ds*- Phi. In the context of the Standard Model, these decays are expected to be highly suppressed since they proceed through annihilation of the b and u-bar quarks in the B- meson. Our results are based on 234 million Upsilon(4S) --> B Bbar decays collected with the BABAR detector at SLAC. We find no evidence for these decays, and we set Bayesian 90% confidence level upper limits on the branching fractions BF(B- --> Ds- Phi) Ds*- Phi)<1.2x10^(-5). These results are consistent with Standard Model expectations.Comment: 8 pages, 3 postscript figues, submitted to Phys. Rev. D (Rapid Communications

Transcriptome analysis of haploid male gametophyte development in Arabidopsis

BACKGROUND: The haploid male gametophyte generation of flowering plants consists of two- or three-celled pollen grains. This functional specialization is thought to be a key factor in the evolutionary success of flowering plants. Moreover, pollen ontogeny is also an attractive model in which to dissect cellular networks that control cell growth, asymmetric cell division and cellular differentiation. Our objective, and an essential step towards the detailed understanding of these processes, was to comprehensively define the male haploid transcriptome throughout development. RESULTS: We have developed staged spore isolation procedures for Arabidopsis and used Affymetrix ATH1 genome arrays to identify a total of 13,977 male gametophyte-expressed mRNAs, 9.7% of which were male-gametophyte-specific. The transition from bicellular to tricellular pollen was accompanied by a decline in the number of diverse mRNA species and an increase in the proportion of male gametophyte-specific transcripts. Expression profiles of regulatory proteins and distinct clusters of coexpressed genes were identified that could correspond to components of gametophytic regulatory networks. Moreover, integration of transcriptome and experimental data revealed the early synthesis of translation factors and their requirement to support pollen tube growth. CONCLUSIONS: The progression from proliferating microspores to terminally differentiated pollen is characterized by large-scale repression of early program genes and the activation of a unique late gene-expression program in maturing pollen. These data provide a quantum increase in knowledge concerning gametophytic transcription and lay the foundations for new genomic-led studies of the regulatory networks and cellular functions that operate to specify male gametophyte development

Temporal Dissection of K-rasG12D Mutant In Vitro and In Vivo Using a Regulatable K-rasG12D Mouse Allele

Animal models which allow the temporal regulation of gene activities are valuable for dissecting gene function in tumorigenesis. Here we have constructed a conditional inducible estrogen receptor-K-rasG12D (ER-K-rasG12D) knock-in mice allele that allows us to temporally switch on or off the activity of K-ras oncogenic mutant through tamoxifen administration. In vitro studies using mice embryonic fibroblast (MEF) showed that a dose of tamoxifen at 0.05 µM works optimally for activation of ER-K-rasG12D independent of the gender status. Furthermore, tamoxifen-inducible activation of K-rasG12D promotes cell proliferation, anchor-independent growth, transformation as well as invasion, potentially via activation of downstream MAPK pathway and cell cycle progression. Continuous activation of K-rasG12D in vivo by tamoxifen treatment is sufficient to drive the neoplastic transformation of normal lung epithelial cells in mice. Tamoxifen withdrawal after the tumor formation results in apoptosis and tumor regression in mouse lungs. Taken together, these data have convincingly demonstrated that K-ras mutant is essential for neoplastic transformation and this animal model may provide an ideal platform for further detailed characterization of the role of K-ras oncogenic mutant during different stages of lung tumorigenesis

Dendritic cells : a double-edged sword in immune responses during chagas disease

Dendritic cells (DCs) are the most important member of the antigen presenting cells group due to their ability to recognize antigen at the infection site and their high specialized antigen internalization capacity. These cells have central role in connecting the innate and adaptive immune responses against Trypanosoma cruzi, the causative agent of Chagas disease. These first line defense cells modulate host immune response depending on type, maturation level, cytokine milieu and DC receptor involved in the interactions with T. cruzi, influencing the development of the disease clinic forms. Here, we present a review of DCs–T. cruzi interactions both in human and murine models, pointing out the parasite ability to manipulate DCs activity for the purpose of evading innate immune response and assuring its own survival and persistence

First RNA-seq approach to study fruit set and parthenocarpy in zucchini (Cucurbita pepo L.)

[EN] Background: Zucchini fruit set can be limited due to unfavourable environmental conditions in off-seasons crops that caused ineffective pollination/fertilization. Parthenocarpy, the natural or artificial fruit development without fertilization, has been recognized as an important trait to avoid this problem, and is related to auxin signalling. Nevertheless, differences found in transcriptome analysis during early fruit development of zucchini suggest that other complementary pathways could regulate fruit formation in parthenocarpic cultivars of this species. The development of next-generation sequencing technologies (NGS) as RNA-sequencing (RNA-seq) opens a new horizon for mapping and quantifying transcriptome to understand the molecular basis of pathways that could regulate parthenocarpy in this species. The aim of the current study was to analyze fruit transcriptome of two cultivars of zucchini, a non-parthenocarpic cultivar and a parthenocarpic cultivar, in an attempt to identify key genes involved in parthenocarpy. Results: RNA-seq analysis of six libraries (unpollinated, pollinated and auxin treated fruit in a non-parthenocarpic and parthenocarpic cultivar) was performed mapping to a new version of C. pepo transcriptome, with a mean of 92% success rate of mapping. In the non-parthenocarpic cultivar, 6479 and 2186 genes were differentially expressed (DEGs) in pollinated fruit and auxin treated fruit, respectively. In the parthenocarpic cultivar, 10,497 in pollinated fruit and 5718 in auxin treated fruit. A comparison between transcriptome of the unpollinated fruit for each cultivar has been performed determining that 6120 genes were differentially expressed. Annotation analysis of these DEGs revealed that cell cycle, regulation of transcription, carbohydrate metabolism and coordination between auxin, ethylene and gibberellin were enriched biological processes during pollinated and parthenocarpic fruit set. Conclusion: This analysis revealed the important role of hormones during fruit set, establishing the activating role of auxins and gibberellins against the inhibitory role of ethylene and different candidate genes that could be useful as markers for parthenocarpic selection in the current breeding programs of zucchini.Research worked is supported by the project RTA2014-00078 from the Spanish Institute of Agronomy Research INIA (Instituto Nacional de Investigacion y Tecnologia Agraria y Alimentaria) and also PP.AVA.AVA201601.7, FEDER y FSE (Programa Operativo FSE de Andalucia 2007-2013 "Andalucia se mueve con Europa"). TPV is supported by a FPI scholarship from RTA2011-00044-C02-01/02 project of INIA. The funding agencies were not involved in the design of the study, collection, analysis, and interpretation of data and in writing the manuscript.Pomares-Viciana, T.; Del Rio-Celestino, M.; Roman, B.; Die, J.; Picó Sirvent, MB.; Gómez, P. (2019). First RNA-seq approach to study fruit set and parthenocarpy in zucchini (Cucurbita pepo L.). BMC Plant Biology. 19:1-20. https://doi.org/10.1186/s12870-019-1632-2S12019Varga A, Bruinsma J. Tomato. In: Monselise SP, editor. CRC Handbook of Fruit Set and Development. Boca Raton: CRC Press; 1986. p. 461–80.Nepi M, Cresti L, Guarnieri M, Pacini E. Effect of relative humidity on water content, viability and carbohydrate profile of Petunia hybrid and Cucurbita pepo pollen. Plant Syst Evol. 2010;284:57–64.Gustafson FG. Parthenocarpy: natural and artificial. Bot Rev. 1942;8:599–654.Robinson RW, Reiners S. Parthenocarpy in summer squash. Hortscience. 1999;34:715–7.Pomares-Viciana T, Die J, Del Río-Celestino M, Román B, Gómez P. Auxin signalling regulation during induced and parthenocarpic fruit set in zucchini. Mol Breeding. 2017;37:56.Ozga JA, Reinecke DM. Hormonal interactions in fruit development. J Plant Growth Regul. 2003;22:73–81.Kim IS, Okubo H, Fujieda K. Endogenous levels of IAA in relation to parthenocarpy in cucumber (Cucumis sativus L). Sci Hortic. 1992;52:1–8.Olimpieri I, Siligato F, Caccia R, Mariotti L, Ceccarelli N, Soressi GP, et al. Tomato fruit set driven by pollination or by the parthenocarpic fruit allele are mediated by transcriptionally regulated gibberellin biosynthesis. Planta. 2007;226:877–88.Cui L, Zhang T, Li J, Lou Q, Chen J. Cloning and expression analysis of Cs-TIR1/AFB2: the fruit development-related genes of cucumber (Cucumis sativus L.). Acta Physiol Plant. 2014;36:139–49.De Jong M, Wolters-Arts J, Feron R, Mariani C, Vriezen WH. The Solanum lycopersicum auxin response factor 7 (SlARF7) regulates auxin signalling during tomato fruit set and development. Plant J. 2009;57:160–70.Wang H, Jones B, Li Z, Frasse P, Delalande C, Regad F, Chaabouni S, Latché A, Pech JC, Bouzayen M. The tomato aux/IAA transcription factor IAA9 is involved in fruit development and leaf morphogenesis. Plant Cell. 2005;17(10):2676–92.Goetz M, Vivian-Smith A, Johnson SD, Koltunow AM. AUXIN RESPONSE FACTOR 8 is a negative regulator of fruit initiation in Arabidopsis. Plant Cell. 2006;18(8):1873–86.Mazzucato A, Cellini F, Bouzayen M, Zouine M, Mila I, Minoia S et al. A TILLING allele of the tomato aux/IAA9 gene offers new insights into fruit set mechanisms and perspectives for breeding seedless tomatoes. Mol Breeding. 2015; 35(22):1-15.Blanca J, Cañizares J, Roig C, Ziarsolo P, Nuez F, Picó B. Transcriptome characterization and high throughput SSRs and SNPs discovery in Cucurbita pepo (Cucurbitaceae). BMC Genomics. 2011;12:104.Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.Da Fonseca RR, Albrechtsen A, Themudo GE, Ramos-Madrigal J, Sibbesen JA, Maretty L, et al. Next-generation biology: sequencing and data analysis approaches for non-model organisms. Mar Genomics. 2016;30:3–13.Conesa A, Madrigal P, Tarazona S, Gomez-Cabrero D, Cervera A, McPherson A, et al. A survey of best practices for RNA-seq data analysis. Genome Biol. 2016;17:13.Li J, Cui ZWJ, Zhang T, Guo Q, Xu J, Li J, et al. Transcriptome comparison of global distinctive features between pollination and parthenocarpic fruit set reveals transcriptional phytohormone cross-talk in cucumber (Cucumis sativus L). Plant Cell Physiol. 2014;55(7):1325–42.Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012;28(23):3150–2.Montero-Pau J, Blanca J, Bombarely A, Ziarsolo P, Esteras C, Martí-Gómez C, et al. De novo assembly of the zucchini genome reveals a whole genome duplication associated with the origin of the Cucurbita genus. Plant Biotechnol J. 2017. https://doi.org/10.1111/pbi.12860 .Vriezen WH, Feron R, Maretto F, Keijman J, Mariani C. Changes in tomato ovary transcriptome demonstrate complex hormonal regulation of fruit set. New Phytol. 2008;177:60–76.Tang N, Deng W, Hu G, Hu N, Li Z. Transcriptome profiling reveals the regulatory mechanism underlying pollination dependent and parthenocarpic fruit set mainly mediated by auxin and gibberellin. PLoS One. 2015;10(4):e0125355.Li J, Yan S, Yang W, Li Y, Xia M, Chen Z, et al. Transcriptomic analysis reveals the roles of microtubule-related genes and transcription factors in fruit length regulation in cucumber (Cucumis sativus L.). Sci Rep. 2015;26(5):8031.Mironov V, De Veylder L, Van Montagu M, Inze D. Cyclin-dependent kinases and cell division in plants- the nexus. Plant Cell. 1999;11(4):509–22.Perrot-Rechenmann C. Cellular responses to auxin: division versus expansion. Cold Spring Harb Perspect Biol. 2010;2(5):a001446.De Veylder L, Beeckman T, Beemster GT, Krols L, Terras F, Landrieu I, et al. Functional analysis of cyclin-dependent kinase inhibitors of Arabidopsis. Plant Cell. 2001;13:1653–68.Nieuwland J, Menges M, Murray JAH. The plant cyclins. In: Inze D, editor. Cell cycle control and plant development, vol. 2007. Oxford: Wiley-Blackwell Publishing; 2007. p. 33–61.Menges M, Samland AK, Planchais S, Murray JA. The D-type cyclin CYCD3;1 is limiting for the G1-to-S-phasetransition in Arabidopsis. Plant Cell. 2006;18:893–906.Boruc J, Mylle E, Duda M, De Clercq R, Rombauts S, Geelen D, et al. Systematic localization of the Arabidopsis core cell cycle proteins reveals novel cell division complexes. Plant Physiol. 2010;152(2):553–65.Sampedro J, Cosgrove DJ. The expansin superfamily. Genome Biol. 2005;6:242.Esmon CA, Tinsley AG, Ljung K, Sandberg G, Hearne LB, Liscum E. A gradient of auxin and auxin-dependent transcription precedes tropic growth responses. Proc Natl Acad Sci. 2006;103:236–41.De Folter S, Busscher J, Colombo L, Losa A, Angenent GC. Transcript profiling of transcription factors genes during siliques development in Arabidopsis. Plant Mol Bio. 2004;56:351–3662004.Son O, Cho HY, Kim MR, Lee H, Lee MS, Song E, et al. Induction of a homeodomain-leucine zipper gene by auxin is inhibited by cytokinin in Arabidopsis roots. Biochem Biophys Res Commun. 2005;326(1):203–9.Olsson ASB, Engstroem P, Seoderman E. The homeobox genes ATHB12 and ATHB7 encode potential regulators of growth in response to water deficit in Arabidopsis. Plant Mol Biol. 2004;55:663–77.Merrow SB, Hopp RJ. Storage effects on winter squashes. Associations between the sugar and starch content of and the degree of preference for winter squashes. J Agric Food Chem. 1961;9:321–6.Berg JM, Tymoczko JL, Stryer L. Carbohydrates. In: Freeman WH, editor. Biochemistry. 5th ed. New York: W H Freeman; 2002.Prabhakar V, Löttgert T, Gigolashvili T, Bell K, Flügge UI, Häusler RE. Molecular and functional characterization of the plastid-localized phosphoenolpyruvate enolase (ENO1) from Arabidopsis thaliana. FEBS Lett. 2009;583(6):983–91.Rius SP, Casati P, Iglesias AA, Gomez-Casati DF. Characterization of Arabidopsis lines deficient in GAPC-1, a cytosolic NAD-dependent glyceraldehyde-3-phosphate dehydrogenase. Plant Physiol. 2008;148(3):1655–67.Van der Linde K, Gutsche N, Leffers HM, Lindermayr C, Müller B, Holtgrefe S, et al. Regulation of plant cytosolic aldolase functions by redox-modifications. Plant Physiol Biochem. 2011;49(9):946–57.Lim H, Cho MH, Jeon JS, Bhoo SH, Kwon YK, Hahn TR. Altered expression of pyrophosphate: fructose-6-phosphate 1-phosphotransferase affects the growth of transgenic Arabidopsis plants. Mol Cells. 2009;27(6):641–9.Baud S, Wuillème S, Dubreucq B, De Almeida A, Vuagnat C, Lepiniec L, et al. Function of plastidial pyruvate kinases in seeds of Arabidopsis thaliana. Plant J. 2007;52:405–19.De Jong M, Mariani C, Vriezen WH. The role of auxin and gibberellin in tomato fruit set. J Exp Bot. 2009;60(5):1523–32.Martínez C, Manzano S, Megías Z, Garrido D, Picó B, Jamilena M. Involvement of ethylene biosynthesis and signalling in fruit set and early fruit development in zucchini squash (Cucurbita pepo L.). BMC Plant Biol. 2013;13:139.Serrani JC, Fos M, Atarés A, Garcia-martinez JL. Effect of gibberellin and auxin on parthenocarpic fruit growth induction in the cv. micro-tom of tomato. J Plant Growth Regul. 2007;26:211–21.Mapelli S. Changes in cytokinin in the fruits of parthenocarpic and normal tomatoes. Plant Sci Lett. 1981;22:227–33.Ulmasov T, Hagen G, Guilfoyle TJ. Activation and repression of transcription by auxin-response factors. Proc Natl Acad Sci U S A. 1999;96:5844–9.Ulmasov T, Hagen G, Guilfoyle TJ. Dimerization and DNA binding of auxin response factors. Plant J. 1999;19:309–19.Tiwari SB, Hagen G, Guilfoyle TJ. Aux/IAA proteins contain a potent transcriptional repression domain. Plant Cell. 2004;16:533–43.Switzenberg JA, Beaudry RM, Grumet R. Effect of CRC:: etr1-1 transgene expression on ethylene production, sex expression, fruit set and fruit ripening in transgenic melon (Cucumis melo L.). Transgenic Res. 2015;24(3):497-507.Nitsch LM, Oplaat C, Feron R, Ma Q, Wolters-Arts M, Hedden P, et al. Abscisic acid levels in tomato ovaries are regulated by LeNCED1 and SlCYP707A1. Planta. 2009;229(6):1335–46.Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-seq. Nat Methods. 2008;5(7):621–8.Robinson MD, McCarthy DJ, Smyth GK. Edger: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2008;26(1):139–40.Raza K, Mishra A. A novel anticlustering filtering algorithm for the prediction of genes as a drug target. Am J Bio Engi. 2012;2(5):206–11.Van Iterson M, Boer JM, Menezes RX. Filtering, FDR and power. BMCBioinformatics. 2010;11:450.Conesa A, Götz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.Berardini TZ, Reiser L, Li D, Mezheritsky Y, Muller R, Strait E, Huala E. The Arabidopsis information resource: making and mining the “gold standard” annotated reference plant genome. Genesis. 2015. https://doi.org/10.1002/dvg.22877 .Bairoch A, Apweiler R. The SWISS-PROT protein sequence database and its supplement TrEMBL. Nucleic Acids Res. 2000;28(1):45–8.Johnson M, Zaretskaya I, Raytselis Y, Merezhuk Y, McGinnis S, Madden TL. NCBI BLAST: a better web interface. Nucleic Acids Res. 2008;36:W5–9.Wyatt LE, Strickler SR, Mueller LA, Mazourek M. An acorn squash (Cucurbita pepo ssp. ovifera) fruit and seed transcriptome as a resource for the study of fruit traits in Cucurbita. Hortic Res. 2015;2:14070. https://doi.org/10.1038/hortres.2014.70 .Zhang A, Ren GA, Sun YA, Guo H, Zhang SA, Zhang FA, et al. A high-density genetic map for anchoring genome sequences and identifying QTLs associated with dwarf vine in pumpkin (Cucurbita maxima Duch.). BMC Genomics. 2015;16:1101.Finn RD, Attwood TK, Babbit AB, Bork P, Bridge AJ, Chang HY. InterPro in 2017-beyond protein family and domain annotations. Nucleic Acids Res. https://doi.org/10.1093/nar/gkw1107 .Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Sherlock G. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9.Kanehisa M, Araki M, Goto S, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:480–4