1
|
Molitor C, Kurowski TJ, Fidalgo de Almeida PM, Kevei Z, Spindlow DJ, Chacko Kaitholil SR, Iheanyichi JU, Prasanna HC, Thompson AJ, Mohareb FR. A chromosome-level genome assembly of Solanum chilense, a tomato wild relative associated with resistance to salinity and drought. Front Plant Sci 2024; 15:1342739. [PMID: 38525148 PMCID: PMC10957597 DOI: 10.3389/fpls.2024.1342739] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/22/2023] [Accepted: 02/12/2024] [Indexed: 03/26/2024]
Abstract
Introduction Solanum chilense is a wild relative of tomato reported to exhibit resistance to biotic and abiotic stresses. There is potential to improve tomato cultivars via breeding with wild relatives, a process greatly accelerated by suitable genomic and genetic resources. Methods In this study we generated a high-quality, chromosome-level, de novo assembly for the S. chilense accession LA1972 using a hybrid assembly strategy with ~180 Gbp of Illumina short reads and ~50 Gbp long PacBio reads. Further scaffolding was performed using Bionano optical maps and 10x Chromium reads. Results The resulting sequences were arranged into 12 pseudomolecules using Hi-C sequencing. This resulted in a 901 Mbp assembly, with a completeness of 95%, as determined by Benchmarking with Universal Single-Copy Orthologs (BUSCO). Sequencing of RNA from multiple tissues resulting in ~219 Gbp of reads was used to annotate the genome assembly with an RNA-Seq guided gene prediction, and for a de novo transcriptome assembly. This chromosome-level, high-quality reference genome for S. chilense accession LA1972 will support future breeding efforts for more sustainable tomato production. Discussion Gene sequences related to drought and salt resistance were compared between S. chilense and S. lycopersicum to identify amino acid variations with high potential for functional impact. These variants were subsequently analysed in 84 resequenced tomato lines across 12 different related species to explore the variant distributions. We identified a set of 7 putative impactful amino acid variants some of which may also impact on fruit development for example the ethylene-responsive transcription factor WIN1 and ethylene-insensitive protein 2. These variants could be tested for their ability to confer functional phenotypes to cultivars that have lost these variants.
Collapse
Affiliation(s)
- Corentin Molitor
- The Bioinformatics Group, School of Water, Energy and Environment, Cranfield University, Wharley End, United Kingdom
| | - Tomasz J. Kurowski
- The Bioinformatics Group, School of Water, Energy and Environment, Cranfield University, Wharley End, United Kingdom
| | | | - Zoltan Kevei
- Soil, Agrifood and Biosciences, Cranfield University, Wharley End, United Kingdom
| | - Daniel J. Spindlow
- The Bioinformatics Group, School of Water, Energy and Environment, Cranfield University, Wharley End, United Kingdom
| | - Steffimol R. Chacko Kaitholil
- The Bioinformatics Group, School of Water, Energy and Environment, Cranfield University, Wharley End, United Kingdom
| | - Justice U. Iheanyichi
- The Bioinformatics Group, School of Water, Energy and Environment, Cranfield University, Wharley End, United Kingdom
| | - H. C. Prasanna
- Division of Vegetable Crops, ICAR-Indian Institute of Horticultural Research, Bangalore, India
| | - Andrew J. Thompson
- Soil, Agrifood and Biosciences, Cranfield University, Wharley End, United Kingdom
| | - Fady R. Mohareb
- The Bioinformatics Group, School of Water, Energy and Environment, Cranfield University, Wharley End, United Kingdom
| |
Collapse
|
2
|
Traut W, Sahara K, ffrench-Constant RH. Lepidopteran Synteny Units reveal deep chromosomal conservation in butterflies and moths. G3 (Bethesda) 2023; 13:jkad134. [PMID: 37310934 PMCID: PMC10411566 DOI: 10.1093/g3journal/jkad134] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/23/2023] [Revised: 05/31/2023] [Accepted: 06/05/2023] [Indexed: 06/15/2023]
Abstract
DNA is compacted into individual particles or chromosomes that form the basic units of inheritance. However, different animals and plants have widely different numbers of chromosomes. This means that we cannot readily tell which chromosomes are related to which. Here, we describe a simple technique that looks at the similarity of genes on each chromosome and thus gives us a true picture of their homology or similarity through evolutionary time. We use this new system to look at the chromosomes of butterflies and moths or Lepidoptera. We term the associated synteny units, Lepidopteran Synteny Units (LSUs). Using a sample of butterfly and moth genomes from across evolutionary time, we show that LSUs form a simple and reliable method of tracing chromosomal homology back through time. Surprisingly, this technique reveals that butterfly and moth chromosomes show conserved blocks dating back to their sister group the Trichoptera. As Lepidoptera have holocentric chromosomes, it will be interesting to see if similar levels of synteny are shown in groups of animals with monocentric chromosomes. The ability to define homology via LSU analysis makes it considerably easier to approach many questions in chromosomal evolution.
Collapse
Affiliation(s)
- Walther Traut
- Institut für Biologie, Zentrum für Medizinische Struktur- und Zellbiologie, Universität zu Lübeck, Ratzeburger Allee 160, D-23562 Lübeck, Germany
| | - Ken Sahara
- Laboratory of Molecular Entomology, Faculty of Agriculture, Iwate University, 3-18-8, Ueda, Morioka 020-8550, Japan
| | | |
Collapse
|
3
|
Arora D, Hernandez AG, Walden KKO, Fields CJ, Yan G. First Draft Genome Assembly of Root-Lesion Nematode Pratylenchus scribneri Generated Using Long-Read Sequencing. Int J Mol Sci 2023; 24:ijms24087311. [PMID: 37108472 PMCID: PMC10138993 DOI: 10.3390/ijms24087311] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2023] [Revised: 04/06/2023] [Accepted: 04/13/2023] [Indexed: 04/29/2023] Open
Abstract
Root-lesion nematodes (genus Pratylenchus) belong to a diverse group of plant-parasitic nematodes (PPN) with a worldwide distribution. Despite being an economically important PPN group of more than 100 species, genome information related to Pratylenchus genus is scarcely available. Here, we report the draft genome assembly of Pratylenchus scribneri generated on the PacBio Sequel IIe System using the ultra-low DNA input HiFi sequencing workflow. The final assembly created using 500 nematodes consisted of 276 decontaminated contigs, with an average contig N50 of 1.72 Mb and an assembled draft genome size of 227.24 Mb consisting of 51,146 predicted protein sequences. The benchmarking universal single-copy ortholog (BUSCO) analysis with 3131 nematode BUSCO groups indicated that 65.4% of the BUSCOs were complete, whereas 24.0%, 41.4%, and 1.8% were single-copy, duplicated, and fragmented, respectively, and 32.8% were missing. The outputs from GenomeScope2 and Smudgeplots converged towards a diploid genome for P. scribneri. The data provided here will facilitate future studies on host plant-nematode interactions and crop protection at the molecular level.
Collapse
Affiliation(s)
- Deepika Arora
- Department of Plant Pathology, North Dakota State University, Fargo, ND 58105, USA
| | - Alvaro G Hernandez
- Roy J. Carver Biotechnology Center, University of Illinois at Urbana-Champaign, Urbana, IL 61820, USA
| | - Kimberly K O Walden
- Roy J. Carver Biotechnology Center, University of Illinois at Urbana-Champaign, Urbana, IL 61820, USA
| | - Christopher J Fields
- Roy J. Carver Biotechnology Center, University of Illinois at Urbana-Champaign, Urbana, IL 61820, USA
| | - Guiping Yan
- Department of Plant Pathology, North Dakota State University, Fargo, ND 58105, USA
| |
Collapse
|
4
|
Ye YX, Li DT, Zhang SY, Shen ZC, Zhang CX. Chromosome-level Genome Assembly and Sex-specific Differential Transcriptome of the White-backed Planthopper, Sogatella furcifera. Curr Genomics 2023; 23:400-411. [PMID: 37920557 PMCID: PMC10173415 DOI: 10.2174/1389202924666230102092822] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Revised: 11/14/2022] [Accepted: 12/16/2022] [Indexed: 01/04/2023] Open
Abstract
Background The white-backed planthopper (WBPH), Sogatella furcifera, causes great damage to many crops (mainly rice) by direct feeding or transmitting plant viruses. The previous genome assembly was generated by second-generation sequencing technologies, with a contig N50 of only 51.5 kb, and contained a lot of heterozygous sequences. Methods We utilized third-generation sequencing technologies and Hi-C data to generate a high-quality chromosome-level assembly. We also provide a large amount of transcriptome data for full-length transcriptome analysis and gender differential expression analysis. Results The final assembly comprised 56.38 Mb, with a contig N50 of 2.20 Mb and a scaffold N50 of 45.25 Mb. Fourteen autosomes and one X chromosome were identified. More than 99.5% of the assembled bases located on the 15 chromosomes. 95.9% of the complete BUSCO Hemiptera genes were detected in the final assembly and 16,880 genes were annotated. 722 genes were relatively highly expressed in males, while 60 in the females. Conclusion The integrated genome, definite sex chromosomes, comprehensive transcriptome profiles, high efficiency of RNA interference and short life cycle substantially made WBPH an efficient research object for functional genomics.
Collapse
Affiliation(s)
- Yu-Xuan Ye
- Institute of Insect Sciences, Zhejiang University, Hangzhou, 310058, China
- The rural development academy, Zhejiang University, Hangzhou, 310058, China
| | - Dan-Ting Li
- Zhejiang Provincial Key Laboratory of Biometrology and Inspection and Quarantine, College of Life Science, China Jiliang University, Hangzhou, 310018, China
| | - Si-Yu Zhang
- Institute of Insect Sciences, Zhejiang University, Hangzhou, 310058, China
| | - Zhi-Cheng Shen
- Institute of Insect Sciences, Zhejiang University, Hangzhou, 310058, China
| | - Chuan-Xi Zhang
- Institute of Insect Sciences, Zhejiang University, Hangzhou, 310058, China
- State Key Laboratory for Managing Biotic and Chemical Threats to the Quality and Safety of Agro-products, Institute of Plant Virology, Ningbo University, Ningbo, 315211, China
| |
Collapse
|
5
|
Tyagi A, Sharma S, Srivastava H, Singh A, Kaila T, Ali S, Gaikwad AB, Singh NK, Gaikwad K. Transcriptome profiling of two contrasting pigeon pea ( Cajanus cajan) genotypes in response to waterlogging stress. Front Genet 2023; 13:1048476. [PMID: 36704340 PMCID: PMC9871235 DOI: 10.3389/fgene.2022.1048476] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2022] [Accepted: 12/22/2022] [Indexed: 01/12/2023] Open
Affiliation(s)
- Anshika Tyagi
- National Institute for Plant Biotechnology, New Delhi, India,Department of Biotechnology Yeungnam University, Gyeongsan, South Korea
| | - Sandhya Sharma
- National Institute for Plant Biotechnology, New Delhi, India
| | | | - Anuradha Singh
- National Institute for Plant Biotechnology, New Delhi, India
| | - Tanvi Kaila
- National Institute for Plant Biotechnology, New Delhi, India
| | - Sajad Ali
- Department of Biotechnology Yeungnam University, Gyeongsan, South Korea
| | | | - N. K Singh
- National Institute for Plant Biotechnology, New Delhi, India
| | - Kishor Gaikwad
- National Institute for Plant Biotechnology, New Delhi, India,*Correspondence: Kishor Gaikwad,
| |
Collapse
|
6
|
Timilsena PR, Barrett CF, Piñeyro-Nelson A, Wafula EK, Ayyampalayam S, McNeal JR, Yukawa T, Givnish TJ, Graham SW, Pires JC, Davis JI, Ané C, Stevenson DW, Leebens-Mack J, Martínez-Salas E, Álvarez-Buylla ER, dePamphilis CW. Phylotranscriptomic Analyses of Mycoheterotrophic Monocots Show a Continuum of Convergent Evolutionary Changes in Expressed Nuclear Genes From Three Independent Nonphotosynthetic Lineages. Genome Biol Evol 2022; 15:6965378. [PMID: 36582124 PMCID: PMC9887272 DOI: 10.1093/gbe/evac183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2021] [Revised: 12/13/2022] [Accepted: 12/18/2022] [Indexed: 12/31/2022] Open
Abstract
Mycoheterotrophy is an alternative nutritional strategy whereby plants obtain sugars and other nutrients from soil fungi. Mycoheterotrophy and associated loss of photosynthesis have evolved repeatedly in plants, particularly in monocots. Although reductive evolution of plastomes in mycoheterotrophs is well documented, the dynamics of nuclear genome evolution remains largely unknown. Transcriptome datasets were generated from four mycoheterotrophs in three families (Orchidaceae, Burmanniaceae, Triuridaceae) and related green plants and used for phylogenomic analyses to resolve relationships among the mycoheterotrophs, their relatives, and representatives across the monocots. Phylogenetic trees based on 602 genes were mostly congruent with plastome phylogenies, except for an Asparagales + Liliales clade inferred in the nuclear trees. Reduction and loss of chlorophyll synthesis and photosynthetic gene expression and relaxation of purifying selection on retained genes were progressive, with greater loss in older nonphotosynthetic lineages. One hundred seventy-four of 1375 plant benchmark universally conserved orthologous genes were undetected in any mycoheterotroph transcriptome or the genome of the mycoheterotrophic orchid Gastrodia but were expressed in green relatives, providing evidence for massively convergent gene loss in nonphotosynthetic lineages. We designate this set of deleted or undetected genes Missing in Mycoheterotrophs (MIM). MIM genes encode not only mainly photosynthetic or plastid membrane proteins but also a diverse set of plastid processes, genes of unknown function, mitochondrial, and cellular processes. Transcription of a photosystem II gene (psb29) in all lineages implies a nonphotosynthetic function for this and other genes retained in mycoheterotrophs. Nonphotosynthetic plants enable novel insights into gene function as well as gene expression shifts, gene loss, and convergence in nuclear genomes.
Collapse
Affiliation(s)
- Prakash Raj Timilsena
- Department of Biology, The Pennsylvania State University, University Park, Pennsylvania
| | - Craig F Barrett
- Department of Biology, West Virginia University, Morgantown, West Virginia
| | - Alma Piñeyro-Nelson
- Departamento de Producción Agrícola y Animal, Universidad Autónoma Metropolitana-Xochimilco, Mexico City, Mexico,Centro de Ciencias de la Complejidad, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | - Eric K Wafula
- Department of Biology, The Pennsylvania State University, University Park, Pennsylvania
| | | | - Joel R McNeal
- Department of Ecology, Evolution, and Organismal Biology, Kennesaw State University, Georgia
| | - Tomohisa Yukawa
- Tsukuba Botanical Garden, National Museum of Nature and Science, 1-1, Amakubo 4, Tsukuba, 305-0005, Japan
| | - Thomas J Givnish
- Department of Botany, University of Wisconsin-Madison, Madison, Wisconsin
| | - Sean W Graham
- Department of Botany, University of British Columbia, Vancouver, British Columbia, V6T 1Z4Canada
| | - J Chris Pires
- Division of Biological Sciences, University of Missouri–Columbia, Columbia, Missouri
| | - Jerrold I Davis
- School of Integrative Plant Sciences and L.H. Bailey Hortorium, Cornell University, Ithaca, New York, 1485
| | - Cécile Ané
- Department of Botany, University of Wisconsin-Madison, Madison, Wisconsin,Department of Statistics, University of Wisconsin–Madison, Madison, Wisconsin
| | | | - Jim Leebens-Mack
- Department of Plant Biology, University of Georgia, Athens, Georgia, 3060
| | - Esteban Martínez-Salas
- Departmento de Botánica, Instituto de Biología, Universidad Nacional Autónoma de México, Mexico City, México
| | - Elena R Álvarez-Buylla
- Departamento de Ecología Funcional, Instituto de Ecología, Universidad Nacional Autónoma de México, Mexico City, Mexico,Centro de Ciencias de la Complejidad, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | | |
Collapse
|
7
|
Timilsena PR, Wafula EK, Barrett CF, Ayyampalayam S, McNeal JR, Rentsch JD, McKain MR, Heyduk K, Harkess A, Villegente M, Conran JG, Illing N, Fogliani B, Ané C, Pires JC, Davis JI, Zomlefer WB, Stevenson DW, Graham SW, Givnish TJ, Leebens-Mack J, dePamphilis CW. Phylogenomic resolution of order- and family-level monocot relationships using 602 single-copy nuclear genes and 1375 BUSCO genes. Front Plant Sci 2022; 13:876779. [PMID: 36483967 PMCID: PMC9723157 DOI: 10.3389/fpls.2022.876779] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2022] [Accepted: 09/29/2022] [Indexed: 05/26/2023]
Abstract
We assess relationships among 192 species in all 12 monocot orders and 72 of 77 families, using 602 conserved single-copy (CSC) genes and 1375 benchmarking single-copy ortholog (BUSCO) genes extracted from genomic and transcriptomic datasets. Phylogenomic inferences based on these data, using both coalescent-based and supermatrix analyses, are largely congruent with the most comprehensive plastome-based analysis, and nuclear-gene phylogenomic analyses with less comprehensive taxon sampling. The strongest discordance between the plastome and nuclear gene analyses is the monophyly of a clade comprising Asparagales and Liliales in our nuclear gene analyses, versus the placement of Asparagales and Liliales as successive sister clades to the commelinids in the plastome tree. Within orders, around six of 72 families shifted positions relative to the recent plastome analysis, but four of these involve poorly supported inferred relationships in the plastome-based tree. In Poales, the nuclear data place a clade comprising Ecdeiocoleaceae+Joinvilleaceae as sister to the grasses (Poaceae); Typhaceae, (rather than Bromeliaceae) are resolved as sister to all other Poales. In Commelinales, nuclear data place Philydraceae sister to all other families rather than to a clade comprising Haemodoraceae+Pontederiaceae as seen in the plastome tree. In Liliales, nuclear data place Liliaceae sister to Smilacaceae, and Melanthiaceae are placed sister to all other Liliales except Campynemataceae. Finally, in Alismatales, nuclear data strongly place Tofieldiaceae, rather than Araceae, as sister to all the other families, providing an alternative resolution of what has been the most problematic node to resolve using plastid data, outside of those involving achlorophyllous mycoheterotrophs. As seen in numerous prior studies, the placement of orders Acorales and Alismatales as successive sister lineages to all other extant monocots. Only 21.2% of BUSCO genes were demonstrably single-copy, yet phylogenomic inferences based on BUSCO and CSC genes did not differ, and overall functional annotations of the two sets were very similar. Our analyses also reveal significant gene tree-species tree discordance despite high support values, as expected given incomplete lineage sorting (ILS) related to rapid diversification. Our study advances understanding of monocot relationships and the robustness of phylogenetic inferences based on large numbers of nuclear single-copy genes that can be obtained from transcriptomes and genomes.
Collapse
Affiliation(s)
- Prakash Raj Timilsena
- Department of Biology and Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, PA, United States
| | - Eric K. Wafula
- Department of Biology and Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, PA, United States
| | - Craig F. Barrett
- Department of Biology, West Virginia University, Morgantown, WV, United States
| | - Saravanaraj Ayyampalayam
- Georgia Advanced Computing Resource Center, University of Georgia, Athens, GA, United States
- Department of Plant Biology, University of Georgia, Athens, GA, United States
| | - Joel R. McNeal
- Department of Ecology, Evolution, and Organismal Biology, Biology Kennesaw State University, Kennesaw, GA, United States
| | - Jeremy D. Rentsch
- Department of Biology, Francis Marion University, Florence, SC, United States
| | - Michael R. McKain
- Department of Biological Sciences, University of Alabama, Tuscaloosa, AL, United States
| | - Karolina Heyduk
- School of Life Sciences, University of Hawai’i at Mānoa, Honolulu, HI, United States
| | - Alex Harkess
- HudsonAlpha Institute for Biotechnology, Huntsville, AL, United States
| | - Matthieu Villegente
- Institut des Sciences Exactes et Appliquees (ISEA), University of New Caledonia, Noumea, New Caledonia
| | - John G. Conran
- Australian Centre for Evolutionary Biology and Biodiversity & Sprigg Geobiology Centre, School of Biological Sciences, University of Adelaide, Adelaide, SA, Australia
| | - Nicola Illing
- Department of Molecular and Cell Biology, University of Cape Town, Cape Town, South Africa
| | - Bruno Fogliani
- Institut des Sciences Exactes et Appliquees (ISEA), University of New Caledonia, Noumea, New Caledonia
| | - Cécile Ané
- Department of Botany, University of Wisconsin-Madison, Madison, WI, United States
- Department of Statistics, University of Wisconsin–Madison, Madison, WI, United States
| | - J. Chris Pires
- Division of Biological Sciences and Bond Life Sciences Center, University of Missouri, Columbia, MO, United States
| | - Jerrold I. Davis
- School of Integrative Plant Sciences and L.H. Bailey Hortorium, Cornell University, Ithaca, NY, United States
| | - Wendy B. Zomlefer
- Department of Plant Biology, University of Georgia, Athens, GA, United States
| | | | | | - Thomas J. Givnish
- Department of Botany, University of Wisconsin-Madison, Madison, WI, United States
| | - James Leebens-Mack
- Department of Plant Biology, University of Georgia, Athens, GA, United States
| | - Claude W. dePamphilis
- Department of Biology and Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, PA, United States
| |
Collapse
|
8
|
Ward CM, Perry KD, Baker G, Powis K, Heckel DG, Baxter SW. A haploid diamondback moth (Plutella xylostella L.) genome assembly resolves 31 chromosomes and identifies a diamide resistance mutation. Insect Biochem Mol Biol 2021; 138:103622. [PMID: 34252570 DOI: 10.1016/j.ibmb.2021.103622] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/25/2021] [Revised: 07/04/2021] [Accepted: 07/04/2021] [Indexed: 05/21/2023]
Abstract
The diamondback moth, Plutella xylostella (L.), is a highly mobile brassica crop pest with worldwide distribution and can rapidly evolve resistance to insecticides, including group 28 diamides. Reference genomes assembled using Illumina sequencing technology have provided valuable resources to advance our knowledge regarding the biology, origin and movement of diamondback moth, and more recently with its sister species, Plutella australiana. Here we apply a trio binning approach to sequence and annotate a chromosome level reference genome of P. xylostella using PacBio Sequel and Dovetail Hi-C sequencing technology and identify a point mutation that causes resistance to commercial diamides. A P. xylostella population collected from brassica crops in the Lockyer Valley, Australia (LV-R), was reselected for chlorantraniliprole resistance then a single male was crossed to a P. australiana female and a hybrid pupa sequenced. A chromosome level 328 Mb P. xylostella genome was assembled with 98.1% assigned to 30 autosomes and the Z chromosome. The genome was highly complete with 98.4% of BUSCO Insecta genes identified and RNAseq informed protein prediction annotated 19,002 coding genes. The LV-R strain survived recommended field application doses of chlorantraniliprole, flubendiamide and cyclaniliprole. Some hybrids also survived these doses, indicating significant departure from recessivity, which has not been previously documented for diamides. Diamide chemicals modulate insect Ryanodine Receptors (RyR), disrupting calcium homeostasis, and we identified an amino acid substitution (I4790K) recently reported to cause diamide resistance in a strain from Japan. This chromosome level assembly provides a new resource for insect comparative genomics and highlights the emergence of diamide resistance in Australia. Resistance management plans need to account for the fact that resistance is not completely recessive.
Collapse
Affiliation(s)
- C M Ward
- School of Biological Sciences, University of Adelaide, 5005, Australia
| | - K D Perry
- South Australian Research and Development Institute, Urrbrae, 5064, Australia
| | - G Baker
- South Australian Research and Development Institute, Urrbrae, 5064, Australia
| | - K Powis
- South Australian Research and Development Institute, Urrbrae, 5064, Australia
| | - D G Heckel
- Department of Entomology, Max Planck Institute for Chemical Ecology, 07745, Jena, Germany
| | - S W Baxter
- Bio21 Institute, School of BioSciences, University of Melbourne, 3052, Australia.
| |
Collapse
|
9
|
Yamaguchi K, Kadota M, Nishimura O, Ohishi Y, Naito Y, Kuraku S. Technical considerations in Hi-C scaffolding and evaluation of chromosome-scale genome assemblies. Mol Ecol 2021; 30:5923-5934. [PMID: 34432923 PMCID: PMC9292758 DOI: 10.1111/mec.16146] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2020] [Revised: 07/28/2021] [Accepted: 08/18/2021] [Indexed: 12/15/2022]
Abstract
The recent development of ecological studies has been fueled by the introduction of massive information based on chromosome‐scale genome sequences, even for species for which genetic linkage is not accessible. This was enabled mainly by the application of Hi‐C, a method for genome‐wide chromosome conformation capture that was originally developed for investigating the long‐range interaction of chromatins. Performing genomic scaffolding using Hi‐C data is highly resource‐demanding and employs elaborate laboratory steps for sample preparation. It starts with building a primary genome sequence assembly as an input, which is followed by computation for genome scaffolding using Hi‐C data, requiring careful validation. This article presents technical considerations for obtaining optimal Hi‐C scaffolding results and provides a test case of its application to a reptile species, the Madagascar ground gecko (Paroedura picta). Among the metrics that are frequently used for evaluating scaffolding results, we investigate the validity of the completeness assessment of chromosome‐scale genome assemblies using single‐copy reference orthologues.
Collapse
Affiliation(s)
- Kazuaki Yamaguchi
- Laboratory for Phyloinformatics, RIKEN Center for Biosystems Dynamics Research, Kobe, Japan
| | - Mitsutaka Kadota
- Laboratory for Phyloinformatics, RIKEN Center for Biosystems Dynamics Research, Kobe, Japan
| | - Osamu Nishimura
- Laboratory for Phyloinformatics, RIKEN Center for Biosystems Dynamics Research, Kobe, Japan
| | - Yuta Ohishi
- Laboratory for Phyloinformatics, RIKEN Center for Biosystems Dynamics Research, Kobe, Japan
| | - Yuki Naito
- Database Center for Life Science (DBCLS), Mishima, Japan
| | - Shigehiro Kuraku
- Laboratory for Phyloinformatics, RIKEN Center for Biosystems Dynamics Research, Kobe, Japan.,Molecular Life History Laboratory, National Institute of Genetics, Mishima, Japan.,Department of Genetics, Sokendai (Graduate University for Advanced Studies), Mishima, Japan
| |
Collapse
|
10
|
Rödelsperger C. The community-curated Pristionchus pacificus genome facilitates automated gene annotation improvement in related nematodes. BMC Genomics 2021; 22:216. [PMID: 33765927 PMCID: PMC7992802 DOI: 10.1186/s12864-021-07529-x] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2020] [Accepted: 03/12/2021] [Indexed: 01/30/2023] Open
Abstract
Background The nematode Pristionchus pacificus is an established model organism for comparative studies with Caenorhabditis elegans. Over the past years, it developed into an independent animal model organism for elucidating the genetic basis of phenotypic plasticity. Community-based curations were employed recently to improve the quality of gene annotations of P. pacificus and to more easily facilitate reverse genetic studies using candidate genes from C. elegans. Results Here, I demonstrate that the reannotation of phylogenomic data from nine related nematode species using the community-curated P. pacificus gene set as homology data substantially improves the quality of gene annotations. Benchmarking of universal single copy orthologs (BUSCO) estimates a median completeness of 84% which corresponds to a 9% increase over previous annotations. Nevertheless, the ability to infer gene models based on homology already drops beyond the genus level reflecting the rapid evolution of nematode lineages. This also indicates that the highly curated C. elegans genome is not optimally suited for annotating non-Caenorhabditis genomes based on homology. Furthermore, comparative genomic analysis of apparently missing BUSCO genes indicates a failure of ortholog detection by the BUSCO pipeline due to the insufficient sample size and phylogenetic breadth of the underlying OrthoDB data set. As a consequence, the quality of multiple divergent nematode genomes might be underestimated. Conclusions This study highlights the need for optimizing gene annotation protocols and it demonstrates the benefit of a high quality genome for phylogenomic data of related species. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07529-x.
Collapse
Affiliation(s)
- Christian Rödelsperger
- Department for Integrative Evolutionary Biology, Max Planck Institute for Developmental Biology, Max-Planck-Ring 9, 72076, Tübingen, Germany.
| |
Collapse
|
11
|
Jauhal AA, Newcomb RD. Assessing genome assembly quality prior to downstream analysis: N50 versus BUSCO. Mol Ecol Resour 2021; 21:1416-1421. [PMID: 33629477 DOI: 10.1111/1755-0998.13364] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2020] [Accepted: 02/16/2021] [Indexed: 12/21/2022]
Abstract
With the ever-increasing number of publicly available eukaryotic genome assemblies and user-friendly bioinformatics tools, there are increasing opportunities for researchers to use genomic resources in their research. While there are multiple dimensions to genome quality, it is often reduced to a single score that may not be correlated with other metrics, or appropriate for all applications of an assembly. To assess whether the commonly reported N50 value could reliably predict a separate dimension of genome quality, gene space completeness, we performed a meta-analysis of 611 published articles on eukaryotic genomes that used BUSCO scores, in addition to the typical N50 score. We found that although assemblies with relatively high contig and scaffold N50 values consistently had high BUSCO scores, a high BUSCO score could also be obtained from assemblies with a low N50. This reinforces that despite its ubiquity, N50 is not a perfect proxy for all measures of genome accuracy. Our data also suggests that variations in BUSCO scores among assemblies with poor N50 scores may be related to the number of introns in conserved eukaryotic genes. We stress the importance of screening and evaluating assembly quality based on the appropriate tools and urge increased reporting of additional genome assessment metrics in addition to N50. We also discuss the potential limitations of BUSCO and suggest improvements for assessing gene space within genome assemblies.
Collapse
Affiliation(s)
- April A Jauhal
- School of Biological Sciences, University of Auckland, Auckland, New Zealand.,The New Zealand Institute for Plant & Food Research, Auckland, New Zealand
| | - Richard D Newcomb
- The New Zealand Institute for Plant & Food Research, Auckland, New Zealand
| |
Collapse
|
12
|
Carneiro CM, Noble JD, Pietras A, Moler P, Austin JD. Iso-seq analysis and functional annotation of the Santa Fe cave crayfish (Procambarus erythrops) transcriptome. Mar Genomics 2021; 58:100842. [PMID: 34217485 DOI: 10.1016/j.margen.2021.100842] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2020] [Revised: 01/08/2021] [Accepted: 01/09/2021] [Indexed: 10/22/2022]
Abstract
The genus Procambarus represents a diverse genus of freshwater crayfish that includes epigean species, stygobitic species, and at least one parthenogenic species. Despite its evolutionary, ecological, and economic importance, most genomic and transcriptomic resources for this genus are limited to a couple of model species. We sequenced the transcriptome of a non-model species, P. erythrops, a geographically restricted stygobitic species from Florida. RNA isolated from gill, muscle and eye tissue was pooled to create a de novo transcriptome assembly using Single Molecule Real-Time sequencing (PacBio), resulting in 19,442 full-length isoforms. The assembly has been deposited in the NCBI (BioProject PRJNA657230). These data will make an important contribution to the comparative study of transcriptome evolution in crayfish and crustaceans.
Collapse
Affiliation(s)
- Celine M Carneiro
- Department of Wildlife Ecology and Conservation, University of Florida, Gainesville, FL 32611, United States of America; School of Natural Resources and Environment, University of Florida, Gainesville, FL 32611, United States of America
| | - Jerald D Noble
- Plant Molecular and Cellular Biology Graduate Program, University of Florida, Gainesville, FL 32611, United States of America
| | - Adele Pietras
- Department of Wildlife Ecology and Conservation, University of Florida, Gainesville, FL 32611, United States of America
| | - Paul Moler
- Florida Fish and Wildlife Conservation Commission, Fish and Wildlife Research Insituute, 1105 SW Williston Road, Gainesville FL32601, United States of America
| | - James D Austin
- Department of Wildlife Ecology and Conservation, University of Florida, Gainesville, FL 32611, United States of America; School of Natural Resources and Environment, University of Florida, Gainesville, FL 32611, United States of America.
| |
Collapse
|
13
|
Dias GR, Dupim EG, Vanderlinde T, Mello B, Carvalho AB. A phylogenomic study of Steganinae fruit flies (Diptera: Drosophilidae): strong gene tree heterogeneity and evidence for monophyly. BMC Evol Biol 2020; 20:141. [PMID: 33138771 DOI: 10.1186/s12862-020-01703-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2020] [Accepted: 10/19/2020] [Indexed: 12/18/2022] Open
Abstract
BACKGROUND The Drosophilidae family is traditionally divided into two subfamilies: Drosophilinae and Steganinae. This division is based on morphological characters, and the two subfamilies have been treated as monophyletic in most of the literature, but some molecular phylogenies have suggested Steganinae to be paraphyletic. To test the paraphyletic-Steganinae hypothesis, here, we used genomic sequences of eight Drosophilidae (three Steganinae and five Drosophilinae) and two Ephydridae (outgroup) species and inferred the phylogeny for the group based on a dataset of 1,028 orthologous genes present in all species (> 1,000,000 bp). This dataset includes three genera that broke the monophyly of the subfamilies in previous works. To investigate possible biases introduced by small sample sizes and automatic gene annotation, we used the same methods to infer species trees from a set of 10 manually annotated genes that are commonly used in phylogenetics. RESULTS Most of the 1,028 gene trees depicted Steganinae as paraphyletic with distinct topologies, but the most common topology depicted it as monophyletic (43.7% of the gene trees). Despite the high levels of gene tree heterogeneity observed, species tree inference in ASTRAL, in PhyloNet, and with the concatenation approach strongly supported the monophyly of both subfamilies for the 1,028-gene dataset. However, when using the concatenation approach to infer a species tree from the smaller set of 10 genes, we recovered Steganinae as a paraphyletic group. The pattern of gene tree heterogeneity was asymmetrical and thus could not be explained solely by incomplete lineage sorting (ILS). CONCLUSIONS Steganinae was clearly a monophyletic group in the dataset that we analyzed. In addition to ILS, gene tree discordance was possibly the result of introgression, suggesting complex branching processes during the early evolution of Drosophilidae with short speciation intervals and gene flow. Our study highlights the importance of genomic data in elucidating contentious phylogenetic relationships and suggests that phylogenetic inference for drosophilids based on small molecular datasets should be performed cautiously. Finally, we suggest an approach for the correction and cleaning of BUSCO-derived genomic datasets that will be useful to other researchers planning to use this tool for phylogenomic studies.
Collapse
|
14
|
Abstract
In daily practice of de novo genome assembly and gene prediction, it would be a natural urge to evaluate their products. Different programs and parameter settings give rise to variable outputs, which leaves a decision of which output to adopt for downstream analysis for addressing biological questions. Instead of superficial assessment of length-based statistics of output sequences (e.g., N50 scaffold length), completeness assessment by means of scoring the coverage of reference orthologs has been increasingly utilized.We previously launched a web service, gVolante ( https://gvolante.riken.jp /), to provide a user-friendly interface and a uniform environment for completeness assessment with the pipelines CEGMA and BUSCO. Completeness assessments performed on gVolante report scores based on not just the coverage of reference genes but also on sequence lengths, allowing quality control in multiple aspects. This chapter focuses on the procedure for such assessment and provides technical tips for higher accuracy.
Collapse
|
15
|
Abstract
Genomics drives the current progress in molecular biology, generating unprecedented volumes of data. The scientific value of these sequences depends on the ability to evaluate their completeness using a biologically meaningful approach. Here, we describe the use of the BUSCO tool suite to assess the completeness of genomes, gene sets, and transcriptomes, using their gene content as a complementary method to common technical metrics. The chapter introduces the concept of universal single-copy genes, which underlies the BUSCO methodology, covers the basic requirements to set up the tool, and provides guidelines to properly design the analyses, run the assessments, and interpret and utilize the results.
Collapse
Affiliation(s)
- Mathieu Seppey
- Department of Genetic Medicine and Development, Swiss Institute of Bioinformatics, University of Geneva Medical School, Geneva, Switzerland
| | - Mosè Manni
- Department of Genetic Medicine and Development, Swiss Institute of Bioinformatics, University of Geneva Medical School, Geneva, Switzerland
| | - Evgeny M Zdobnov
- Department of Genetic Medicine and Development, Swiss Institute of Bioinformatics, University of Geneva Medical School, Geneva, Switzerland.
| |
Collapse
|
16
|
Briones C, Nuñez JJ, Pérez M, Espinoza-Rojas D, Molina-Quiroz C, Guiñez R. De novo male gonad transcriptome draft for the marine mussel Perumytilus purpuratus with a focus on its reproductive-related proteins. J Genomics 2018; 6:127-132. [PMID: 30510598 PMCID: PMC6275399 DOI: 10.7150/jgen.27864] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2018] [Accepted: 09/04/2018] [Indexed: 12/17/2022] Open
Abstract
Perumytilus purpuratus is a marine mussel considered a bioengineer species with a broad distribution in the Pacific and Atlantic coast of South America. Studies have shown two geographically and genetically differentiated subpopulations at molecular level and in sperm morphological traits. To open avenues for molecular research on P. purpuratus, a global de novo transcriptome from gonadal tissue of mature males was sequenced using the Illumina platform. From a total of 126.38 million reads, 37,765 transcripts were successfully annotated. BUSCO analysis determined a level of 89% completeness for the assembled transcriptome. The functional gene ontology (GO) annotation indicated that, in terms of abundance, the transcripts related with molecular function were the most represented, followed by those related with biological process and cellular components. Additionally, a subset of GO annotations generated using the "sperm" term resulted in a total of 1,294 sequences where the biological process category was the more represented, with transcripts strongly associated to sperm-processes required for fertilization, and with processes where the sperm-egg interaction could be implicated. Our work will contribute to the evolutionary understanding of the molecular mechanisms related to tissue-specific functions. This work reports the first male gonad transcriptome for the mussel P. purpuratus, generating a useful transcriptomic resource for this species and other closely related mytilids.
Collapse
Affiliation(s)
- Carolina Briones
- Instituto de Ciencias Marinas y Limnológicas, Universidad Austral de Chile, Casilla 567, Valdivia, Chile
| | - José J Nuñez
- Instituto de Ciencias Marinas y Limnológicas, Universidad Austral de Chile, Casilla 567, Valdivia, Chile
| | - Montse Pérez
- Instituto Español de Oceanografía, Centro Oceanográfico de Vigo, 36200 Vigo, España
| | | | | | - Ricardo Guiñez
- Instituto de Ciencias Naturales Alexander von Humboldt, Facultad de Ciencias del Mar y Recursos Biológicos, Universidad de Antofagasta, Casilla 170, Antofagasta, Chile
| |
Collapse
|
17
|
Carruthers M, Yurchenko AA, Augley JJ, Adams CE, Herzyk P, Elmer KR. De novo transcriptome assembly, annotation and comparison of four ecological and evolutionary model salmonid fish species. BMC Genomics 2018; 19:32. [PMID: 29310597 PMCID: PMC5759245 DOI: 10.1186/s12864-017-4379-x] [Citation(s) in RCA: 52] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2017] [Accepted: 12/11/2017] [Indexed: 12/21/2022] Open
Abstract
Background Salmonid fishes exhibit high levels of phenotypic and ecological variation and are thus ideal model systems for studying evolutionary processes of adaptive divergence and speciation. Furthermore, salmonids are of major interest in fisheries, aquaculture, and conservation research. Improving understanding of the genetic mechanisms underlying traits in these species would significantly progress research in these fields. Here we generate high quality de novo transcriptomes for four salmonid species: Atlantic salmon (Salmo salar), brown trout (Salmo trutta), Arctic charr (Salvelinus alpinus), and European whitefish (Coregonus lavaretus). All species except Atlantic salmon have no reference genome publicly available and few if any genomic studies to date. Results We used paired-end RNA-seq on Illumina to generate high coverage sequencing of multiple individuals, yielding between 180 and 210 M reads per species. After initial assembly, strict filtering was used to remove duplicated, redundant, and low confidence transcripts. The final assemblies consisted of 36,505 protein-coding transcripts for Atlantic salmon, 35,736 for brown trout, 33,126 for Arctic charr, and 33,697 for European whitefish and are made publicly available. Assembly completeness was assessed using three approaches, all of which supported high quality of the assemblies: 1) ~78% of Actinopterygian single-copy orthologs were successfully captured in our assemblies, 2) orthogroup inference identified high overlap in the protein sequences present across all four species (40% shared across all four and 84% shared by at least two), and 3) comparison with the published Atlantic salmon genome suggests that our assemblies represent well covered (~98%) protein-coding transcriptomes. Thorough comparison of the generated assemblies found that 84-90% of transcripts in each assembly were orthologous with at least one of the other three species. We also identified 34-37% of transcripts in each assembly as paralogs. We further compare completeness and annotation statistics of our new assemblies to available related species. Conclusion New, high-confidence protein-coding transcriptomes were generated for four ecologically and economically important species of salmonids. This offers a high quality pipeline for such complex genomes, represents a valuable contribution to the existing genomic resources for these species and provides robust tools for future investigation of gene expression and sequence evolution in these and other salmonid species. Electronic supplementary material The online version of this article (10.1186/s12864-017-4379-x) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Madeleine Carruthers
- Institute of Biodiversity, Animal Health & Comparative Medicine, College of Medical, Veterinary & Life Sciences, University of Glasgow, G12 8QQ, Glasgow, UK
| | - Andrey A Yurchenko
- Institute of Biodiversity, Animal Health & Comparative Medicine, College of Medical, Veterinary & Life Sciences, University of Glasgow, G12 8QQ, Glasgow, UK
| | - Julian J Augley
- Glasgow Polyomics, Wolfson Wohl Cancer Research Centre, University of Glasgow, G61 1QH, Glasgow, UK.,Present Address: Fios Genomics Ltd., Nine Edinburgh Bioquarter, 9 Little France Road, Edinburgh, EH16 4UX, UK
| | - Colin E Adams
- Institute of Biodiversity, Animal Health & Comparative Medicine, College of Medical, Veterinary & Life Sciences, University of Glasgow, G12 8QQ, Glasgow, UK.,Scottish Centre for Ecology and the Natural Environment, University of Glasgow, Rowardennan, G63 0AW, UK
| | - Pawel Herzyk
- Glasgow Polyomics, Wolfson Wohl Cancer Research Centre, University of Glasgow, G61 1QH, Glasgow, UK.,Institute of Molecular, Cell & Systems Biology, College of Medical, Veterinary & Life Sciences, University of Glasgow, G12 8QQ, Glasgow, UK
| | - Kathryn R Elmer
- Institute of Biodiversity, Animal Health & Comparative Medicine, College of Medical, Veterinary & Life Sciences, University of Glasgow, G12 8QQ, Glasgow, UK.
| |
Collapse
|
18
|
Theissinger K, Falckenhayn C, Blande D, Toljamo A, Gutekunst J, Makkonen J, Jussila J, Lyko F, Schrimpf A, Schulz R, Kokko H. De Novo assembly and annotation of the freshwater crayfish Astacus astacus transcriptome. Mar Genomics 2016; 28:7-10. [PMID: 26988698 DOI: 10.1016/j.margen.2016.02.006] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2016] [Revised: 02/01/2016] [Accepted: 02/22/2016] [Indexed: 02/07/2023]
Abstract
We generated RNA-seq data to assemble the transcriptome of the noble crayfish (Astacus astacus) from four combined tissues (abdominal muscle, hepatopancreas, ovaries, green glands). A total of 194 million read pairs with a length of 100 bp were generated. The transcriptome was assembled de novo using Trinity software, producing 158,649 non-redundant transcripts. Lowly expressed transcripts were filtered out leaving 45,415 transcripts of which 14,559 were found to contain open reading frames with predicted gene function. The Transrate software revealed that 91% of the total reads were realigned to the assembly. Furthermore, BUSCO analysis indicated that our assembly is 64% complete. A total of 13,770 transcripts were assigned at least one GO term. This first de novo transcriptome assembly is an important foundation for future genomic research on the noble crayfish and adds to the general knowledge and further characterization of transcriptomes of non-model organisms.
Collapse
Affiliation(s)
- Kathrin Theissinger
- University of Koblenz-Landau, Institute for Environmental Sciences, Fortstrasse 7, 76829 Landau, Germany.
| | - Cassandra Falckenhayn
- Division of Epigenetics, DKFZ-ZMBH Alliance, German Cancer Research Center, Im Neuenheimer Feld 580, 69120 Heidelberg, Germany
| | - Daniel Blande
- The University of Eastern Finland, Department of Environmental and Biological Sciences, P.O. Box 1627, 70211 Kuopio, Finland
| | - Anna Toljamo
- The University of Eastern Finland, Department of Environmental and Biological Sciences, P.O. Box 1627, 70211 Kuopio, Finland
| | - Julian Gutekunst
- Division of Epigenetics, DKFZ-ZMBH Alliance, German Cancer Research Center, Im Neuenheimer Feld 580, 69120 Heidelberg, Germany
| | - Jenny Makkonen
- The University of Eastern Finland, Department of Environmental and Biological Sciences, P.O. Box 1627, 70211 Kuopio, Finland
| | - Japo Jussila
- The University of Eastern Finland, Department of Environmental and Biological Sciences, P.O. Box 1627, 70211 Kuopio, Finland
| | - Frank Lyko
- Division of Epigenetics, DKFZ-ZMBH Alliance, German Cancer Research Center, Im Neuenheimer Feld 580, 69120 Heidelberg, Germany
| | - Anne Schrimpf
- University of Koblenz-Landau, Institute for Environmental Sciences, Fortstrasse 7, 76829 Landau, Germany
| | - Ralf Schulz
- University of Koblenz-Landau, Institute for Environmental Sciences, Fortstrasse 7, 76829 Landau, Germany
| | - Harri Kokko
- The University of Eastern Finland, Department of Environmental and Biological Sciences, P.O. Box 1627, 70211 Kuopio, Finland
| |
Collapse
|