Analysis of chromosome karyotype and genome size in echiuran Urechis unicinctus Drasche, 1880 (Polychaeta, Urechidae)

Comparative Cytogenetics, Mar 2019

Karyotype and genome size are two primary cytogenetic characteristics of species, which are of great significance to the study of cytogenetics, taxonomy, phylogenesis, evolution as well as molecular biology. However, this basic cytogenetic information in echiurans is lacking. Therefore, we analyzed characteristics of karyotype and genome size in the echiuran worm Urechis unicinctus Drasche, 1880. In this study, coelomic cells of U. unicinctus were used for analyzing the genome size by a flow cytometry with chicken erythrocytes as DNA standard, and the 2C DNA content was determined to be 1.85 pg, which was corresponded to the genome size of 904.58 Mbp approximately. Furthermore, trochophores of U. unicinctus were dissociated and cells were utilized for preparing the chromosomes stained with DAPI, and the karyotype was determined as 2n = 30 (10m + 6sm + 6st + 8t), FN=52. Our data provided the basic cytogenetic information of U. unicinctus, which could be utilized in taxonomic study and whole-genome sequencing in future.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://compcytogen.pensoft.net/article/31448/download/pdf/

Analysis of chromosome karyotype and genome size in echiuran Urechis unicinctus Drasche, 1880 (Polychaeta, Urechidae)

Analysis of chromosome karyotype and genome size in echiuran Urechis unicinctus Drasche Analysis of chromosome karyotype and genome size in echiuran Urechis unicinctus Drasche, 1880 (Polychaeta, Urechidae) Zhenkui Qin 0 Xueyu Li 0 Danwen Liu 0 Qing Wang 0 Li Lu 0 Zhifeng Zhang 0 Corresponding author: Zhifeng Zhang (zzfp @ouc.edu.cn) 0 Ministry of Education Key Laboratory of Marine Genetics and Breeding, College of Marine Life Sciences, Ocean University of China , Qingdao 266003 , China Karyotype and genome size are two primary cytogenetic characteristics of species, which are of great significance to the study of cytogenetics, taxonomy, phylogenesis, evolution as well as molecular biology. However, this basic cytogenetic information in echiurans is lacking. Therefore, we analyzed characteristics of karyotype and genome size in the echiuran worm Urechis unicinctus Drasche, 1880. In this study, coelomic cells of U. unicinctus were used for analyzing the genome size by a flow cytometry with chicken erythrocytes as DNA standard, and the 2C DNA content was determined to be 1.85 pg, which was corresponded to the genome size of 904.58 Mbp approximately. Furthermore, trochophores of U. unicinctus were dissociated and cells were utilized for preparing the chromosomes stained with DAPI, and the karyotype was determined as 2n = 30 (10m + 6sm + 6st + 8t), FN=52. Our data provided the basic cytogenetic information of U. unicinctus, which could be utilized in taxonomic study and whole-genome sequencing in future. eol>Urechis unicinctus; karyotype; genome size; flow cytometry Introduction Echiurans (spoon worms) are a group of marine worms which are unsegmented, coelomate, bilaterally symmetrical and soft-bodied (Goto 2016) . Traditionally, echiurans have been excluded from Annelida because of their non-segmented characteristics (Fisher and MacGinitie 1928) . Recently, based on molecular phylogenetic data (Struck et al. 2007, Wu et al. 2009, Struck et al. 2011, Weigert et al. 2014, Andrade et al. 2015, Goto 2016; ) the species in Echiura have often been considered as a group of derived annelids that secondarily lost segmentation. Thus, a controversial issue of whether echiurans belong to Annelida or a separate Echiura phylum has emerged and this needs to be elucidated from different research scopes. Chromosome karyotype and genome size are two important cytogenetic characteristics and have been applied widely in taxonomic, phylogenetic and evolutionary studies (Ipucha et al. 2007, Kashmenskaya and Polyakov 2008, Leit?o et al. 2010, Palomina et al. 2017) . Regrettably, only very few studies about chromosomes in the echiuran worms are available (Griffin 1899, Lefevre 1907, Singhal and Dattagupta 1980) , which were reported several decades ago, and no report was related to their genome size. Therefore, these basic cytogenetic characteristics of echiurans need to be revealed urgently. Karyotype, including chromosome number and composition, could reflect the taxonomic relationship between species and be used as a tool to explore biological diversity (Dobigny et al. 2004, Leit?o et al. 2010, Cioffi et al. 2012) . Gallardo-Esc?rateand Del R?o-Portilla (2007) analyzed the cytogenetical relationships of three abalone species Haliotis corrugata W. Wood, 1828, H. fulgens Philippi, 1845, and H. rufescens Swainson, 1822 based on their chromosomal morphology, and proposed that H. rufescens and H. corrugata are cytogenetically more similar to each other than to H. fulgens. Ipucha et al. (2007) discussed the phylogenetic relationship by analyzing the karyotypes of seven species from Nereididae, and suggested the karyotypes are relatively similar and stable in nereidid species at the family level, while the main mechanism of chromosomal evolution could be pericentric inversions. Genome size is the total DNA content within a single copy genome and is also referred to C-value, which is specific in every species and ranges from 0.02 pg (Pratylenchus coffeae Zimmermann, 1898, a plantparasitic nematode) to 132.83 pg (Protopterus aethiopicus Heckel, 1851, a marbled lungfish) in animals (Gregory 2018) . C-value estimation is important for genomic sequencing and analysis (Gregory 2005) . Nevertheless, variation of genome size in different species is rarely used as a direct or single factor in evolution analysis due to the C-value paradox (the phenomenon that C-value is inconsistent with the complexity of biological structure or composition), which means the genome size among organisms were diverse and possessed no relationship to organismal complexity (Gregory 2010) . Urechis unicinctus, a commercial echiuran worm inhabiting the U-shaped burrows in the coastal mud flats, has unique roles in animal evolution, coastal sediment improvement and marine drug development (Liu et al. 2015) . In this study, chromosome counting and composition? analysis of U. unicinctus were carried out in the well-developed trochophore for the first time, and the genome size of U. unicinctus coelomic cells was also determined using flow cytometry. We aim to reveal the basic cytogenetic characteristics of U. unicinctus, which provide useful information for taxonomic and genomic studies. Material and methods Animals U. unicinctus adults with 9.96 ? 0.42 cm in body length were purchased from an aquatic product market, which were collected from a coastal intertidal flat in Yantai, China. Chromosome preparation and karyotype analysis Sampling The mature sperms and oocytes were obtained by dissecting the nephridia of the healthy U. unicinctus, respectively. Artificial fertilization was conducted by mixing sperms and oocytes at a ratio of 10:1 in filter seawater (FSW), and then these fertilized eggs were cultivated until hatched in FSW (19.7 ? 0.3 ?C, salinity 29 PSU, pH 8.29 ? 0.02). The hatched trochophores were collected using a 500 mesh sieve. Chromosome preparation and karyotype analysis Chromosomes of U. unicinctus trochophores were prepared as described by Earley (1975) with some modifications. The larvae were incubated in FSW containing 0.02% colchicine for 2.5 h at room temperature, and then transferred into Ca2+/Mg2+-free Dulbecco?s phosphate buffered saline (CMF-DPBS, 137 mM NaCl, 2.7 mM KCl, 8.1 mM Na2HPO4?7H2O, 1.1 mM KH2PO4) for 30 min with continuous mild shaking to obtain the dissociated larval cells. Successively, the cells were treated with 0.075 mol/L KCl for 30 min, collected by centrifuging at 1200 g for 5 min, and then fixed three times with cold Carnoy?s fixative (3 ethanol: 1 glacial acetic acid) for 15 min each. After centrifugation, the cells were re-suspended in 50% glacial acetic acid aqueous solution, and then dropped onto preheated clean glass slides at 56 ?C, and air-dried. Finally, these samples were stained with 10 ?g/ml 4?, 6-diamidino-2-phenylindole (DAPI) (Solarbio, China) for 15 min, and were examined and photographed under a Nikon Eclipse 80i fluorescence microscope (Nikon, Tokyo, Japan). The chromosome lengths, chromosome relative lengths and arm ratios from wellformed chromosomes in metaphase were measured and calculated using MICROMEASURE 3.3 software. Data were presented as mean ? SEM (n = 5). Chromosomes were classified according to the description of Levan et al. (1964) , and the homologous chromosome was assigned based on the similarities in length and centromere position using PHOTOSHOP CS6 software. The idiogram was constructed according to the arm ratio and relative length of the chromosomes. The karyotype was classified as described by Stebbins (1971) . Estimation of DNA content Preparation of single cell samples Coelomic fluids from 11 healthy worms were collected by syringe puncturing U. unicinctus body wall, respectively, and three duplicate samples were obtained from each individual. The coelomic cells were pelleted by centrifugation at 1000 g (4 ?C) for 5 min, washed three times with PBS (pH 7.2), and then resuspended with PBS (pH 7.2). The suspension was added dropwise to the precooled 70% ethanol and fixed overnight at 4 ?C. The next day, cell samples were collected through a 50 ?m nylon mesh filtration, adjusted to 5?105 cells/ml, digested with 20 ?g/?l RNase A for 10 min and then stained with 1 ?g/?l propidium iodide (PI) for 30 min in the dark at room temperature. Chicken erythrocytes (2C = 25 pg DNA) were used as an internal standard (Darzynkiewicz and Juan 1997) . Fresh blood was acquired by puncturing the heart and mixed with 5% sodium citrate to prevent coagulation, and the subsequential processing procedures were as described in above. Flow cytometric analysis Twelve samples, including a chicken erythrocyte, a U. unicinctus coelomic cell, and ten mixed samples containing 500 ?l erythrocytes and 500 ?l coelomic cells, were analyzed using a Coulter Cytomics FC500-MPL flow cytometer (Beckman, California, USA) equipped with a 488 nm laser source to detect the DNA content, and the output was processed in the software FLOWJO 7.6.1. Coefficients of variation (CV) were adjusted below 5% to ensure the reliability. The DNA content of U. unicinctus was then calculated according to Dole?el et al. (2007) using the equation: Y = N/M ? X (Y means the 2C DNA content of U. unicinctus, N means the fluorescence mean values of U. unicinctus samples, M means the fluorescence mean values of internal standard (chicken erythrocytes), and X means the 2C DNA content of internal standard). The genome size of U. unicinctus was calculated using the equation: genome size (bp) = (0.978 ? 109) ? DNA content (pg). Results Chromosome characteristics of U. unicinctus The well dispersed metaphase chromosomes from U. unicinctus trochophore cells are shown in Fig. 1. The diploid chromosome number of U. unicinctus was 30, and fundamental number (FN) was 52. Fifteen pairs of the homologous chromosomes were matched, including 5 pairs of metacentric (m), 3 pairs of submetacentric (sm), 3 pairs of subtelocentric (st) and 4 pairs of telocentric chromosomes (t). No secondary constriction or satellite was found here. According to the measurement data (Table 1), the karyotype formula was deduced as 2n = 30 (10m + 6sm + 6st + 8t) (Fig. 1). Moreover, the idiogram was drawn based on the results above (Fig. 2). The index of karyotypic asymmetry (AsK) was 74.4%, which was the ratio of total length of long arms to that of? all chromosomes. The chromosome size was 2.36?5.93 ?m, with the longest to shortest chromosome length ratio (L/S) 2.51, and the percentage of chromosomes with arm ratio greater than 2:1 was 66.7%. Therefore, the karyotype of U. unicinctus was classified as 3B. Genome size of U. unicinctus The frequency histograms of DNA content from chicken erythrocytes and U. unicinctus coelomic cells were presented based on the flow cytometric analyses (Fig. 3). No overlap between the two peaks (Fig. 3C) indicated that chicken erythrocytes as the internal standard was suitable for DNA content determination of U. unicinctus coelomic cells. The 2C mean values of chicken erythrocytes (M) and U. unicinctus coelomic cells (N) in ten mixed samples and their ratios were presented in Table 2. The results showed p, length of short arm; q, length of long arm; Total, total length of chromosome; RL, relative length of chromosome; AR, arm ratio of long arm to short arm from metaphase chromosomes. m, metacentric chromosome; sm, submetacentric chromosome; st, subtelocentric chromosome; t, telocentric chromosome. that the average ratio of U. unicinctus coelomic cells to chicken erythrocytes was 0.74, therefore the 2C DNA content of U. unicinctus was calculated to be 1.85 pg, and its genome size was 904.58 Mb. Discussion Studies on the chromosomes in echiurans are very limited, and all of them were conducted several decades ago (Griffin 1899, Lefevre 1907, Singhal and Dattagupta 1980) . The only one study of the karyotype was performed in a Bonellidae species, Achaetobonellia maculata Fisher, 1953, (Singhal and Dattagupta 1980) , which has 2n? =? 20? (10m), whereas others were focusing on the status and motion of chromosomes during cell division. Until now, no more karyotypic information of echiurans has been investigated. In this study, we obtained the clear karyotype of U. unicinctus using fluorescent staining technique and estimated its genome size as well. The karyotype of U. unicinctus was 2n = 30 (10m + 6sm + 6st + 8t), the 2C DNA content was 1.85 pg, and the genome size was 904.58 Mb approximately. This is the first study conducted in an Urechidae animal. Karyotypic information could be utilized to study the taxonomic relationships of species and biological diversity (Dobigny et al. 2004, Ipucha et al. 2007, Leit?o et al. 2010, Cioffi et al. 2012) . In this study, we collected the karyotypic data of multiple 1A 2A 3A 4B 2B 2B 3B 3B 2B Reference This study Singhal and Dattagupta 1980 Wang et al. 2008 Shi et al. 2013 Ipucha et al. 2007 Ipucha et al. 2007 Leit?o et al. 2010 Anisimov et al. 2015 Kashmenskaya and Polyakov?2008 Kashmenskaya and Polyakov?2008 echiurans, sipunculids and annelids, and made some comparisons. It appeared that U. unicinctus possesses similar number and morphology of chromosomes with annelids (Table 3), and its karyotypic asymmetry was closely concentrated with annelids (Fig.?4). However, more data and analysis were required to determine the phylogenetic relationship between Echiurans and Annelids in future. In general, genetic information of higher organisms is more complex than that of lower organisms, so the genomic size of higher organisms is relatively greater. However, there is no inevitable correlation between genome size and organismal complexity, because genome often contains a large number of highly repetitive DNA sequences, resulting in the conflict of DNA content and its evolutionary level. Gregory and Hebert estimated genome sizes from 12 species of freshwater oligochaetes ranging from 0.8 to 7.6 pg, and 15 species of earthworms varied from 0.4 to 1.2 pg (Gregory and Hebert 2002) , suggesting that there is such a wide variation in the DNA content even between related species. Variation in genome size of polychaete taxa is not evenly distributed, as species inhabiting interstitial environments have smaller size (0.06?1.1?pg), whereas macrobenthic species are larger (0.4?7.2 pg), and the difference has been considered to adaptation of different environments (Soldi et al. 1994, Gambi et al. 1997, Gregory 2018) . In addition, the DNA content among different species was also found to be independent of chromosome numbers, which was also concluded by EI-Shehawi and Elseehy (2017) that no correlation between genome size and chromosome number after the comparison of more than 6000 records. In echiurans, there has been no report of the nuclear genome size up to now. In the present study, the genome size of U.?unicinctus was estimated as 1.85 pg, which is relatively small and could also be resulted from the adaptation to the hash and variable intertidal environment. The determination of U. unicinctus genome size maybe of little significance for the study of the evolutionary status of Echiura, but it could provide effective data support for large-scale whole-genome sequencing of U. unicinctus in the near future. Conclusion In the present study, the karyotype of an Urechidae animal, U. unicinctus, was discovered for the first time as 2n = 30 (10m + 6sm + 6st + 8t), FN=52. Meanwhile, the 2C DNA content was detected to be 1.85 pg and its genome size was estimated as 904.58?Mb. Our study provided effective cytogenetic information for taxonomic study and whole-genome sequencing of U. unicinctus. Acknowledgement This work was supported by the National Natural Science Foundation of China (31572601). Andrade SC , Novo M , Kawauchi GY , Worsaae K , Pleijel F , Giribet G , Rouse GW ( 2015 ) Articulating ?archiannelids?: Phylogenomics and annelid relationships, with emphasis on meiofaunal taxa . Molecular Biology and Evolution 32 : 2860 - 2875 . https://doi.org/10.1093/ molbev/msv157 Anisimov AP , Roslik GV , Ganin GN ( 2015 ) Cytogenetic description of the earthworm Drawida ghilarovi Gates, 1969 (Oligochaeta, Moniligastridae) from the southern Russian Far East . Comparative Cytogenetics 9 : 565 - 577 . https://doi.org/10.3897/CompCytogen.v9i4. 5741 Cioffi MdB , Molina WF , Antonio RF , Bertollo LAC ( 2012 ) Chromosomes as tools for discovering biodiversity-The case of Erythrinidae fish family . In: Tirunilai P (Ed.) Recent Trends in Cytogenetic Studies-Methodologies and Applications . Intech, Rijeka, 125 - 146 . Darzynkiewicz Z , Juan G ( 1997 ) DNA content measurement for DNA ploidy and cell cycle analysis . Current Protocols in Cytometry 7 .5: 1 - 24 . https://doi.org/10.1002/0471142956.cy0705s00 Dobigny G , Ducroz J-F , Robinson TJ , Volobouev V , Sites J ( 2004 ) Cytogenetics and cladistics . Systematic Biology 53 : 470 - 484 . https://doi.org/10.1080/10635150490445698 Dole?el J , Greilhuber J , Suda J ( 2007 ) Estimation of nuclear DNA content in plants using flow cytometry . Nature Protocols 2 : 2233 - 2244 . https://doi.org/10.1038/nprot. 2007 .310 Earley EM ( 1975 ) Chromosome preparations from monolayer cell cultures . TCA manual / Tissue Culture Association 1 : 31 - 35 . https://doi.org/10.1007/BF00914434 EI-Shehawi AM , Elseehy MM ( 2017 ) Genome size and chromosome number relationship contradicts the principle of darwinian evolution from common ancestor . Journal of Phylogenetics and Evolutionary Biology 5 : 179 . Fisher WK , MacGinitie GE ( 1928 ) XVII-The natural history of an echiuroid worm . Journal of Natural History 1 : 204 - 213 . https://doi.org/10.1080/00222932808672765 Gallardo-Esc?rate C , Del R?o-Portilla M ( 2007 ) Karyotype composition in three California abalones and their relationship with genome size . Journal of Shellfish Research 26 : 825 - 832 . https://doi.org/10.2983/ 0730 - 8000 ( 2007 ) 26 [825:KCITCA] 2.0 .CO; 2 Gambi MC , Ramella L , Sella G , Protto P , Aldieri E ( 1997 ) Variation in genome size of benthic polychaetes: systematic and ecological relationship . Journal of the Marine Biological Association of the United Kingdom 77 : 1045 - 1057 . https://doi.org/10.1017/ S0025315400038613 Goto R ( 2016 ) A comprehensive molecular phylogeny of spoon worms (Echiura, Annelida): Implications for morphological evolution, the origin of dwarf males, and habitat shifts . Molecular Phylogenetics and Evolution 99 : 247 - 260 . https://doi.org/10.1016/j. ympev. 2016 . 03 .003 Gregory TR ( 2005 ) Genome Size Evolution in Animals . In: Gregory TR (Ed.): The Evolution of the Genome . Elsevier, San Diego, 4- 71 . https://doi.org/10.1016/B978-012301463-4/ 50003 - 6 Gregory TR ( 2010 ) Coincidence, coevolution, or causation? DNA content, cell size, and the C-value enigma . Biological Reviews 76 : 65 - 101 . Gregory TR ( 2018 ) Animal Genome Size Database . http://www.genomesize. com [accessed 7 November 2018 ] Gregory TR , Hebert PDN ( 2002 ) Genome size estimates for some oligochaete annelids . Canadian Journal of Zoology 80 : 1485 - 1489 . https://doi.org/10.1139/z02- 145 Griffin BB ( 1899 ) Studies on the maturation, fertilization, and cleavage of Thalassema and Zirphaea . Journal of Morphology 15 : 583 - 634 . https://doi.org/10.1002/jmor.1050150302 Ipucha MC , Santos GC , Lana PDC , Sbalqueiro I ( 2007 ) Cytogenetic characterization of seven South American species of Nereididae (Annelida: Polychaeta): implications for the karyotypic evolution . Journal of Basic and Applied Genetics 18 : 27 - 38 . Kashmenskaya MN , Polyakov AV ( 2008 ) Karyotype analysis of five species of earthworms (Oligochaeta: Lumbricidae) . Comparative Cytogenetics 2 : 121 - 125 . Leit?o A , Carvalho S , Ben-Hamadou R , Gaspar MB ( 2010 ) Cytogenetics of Hediste diversicolor (Annelida: Polychaeta) and comparative karyological analysis within Nereididae . Aquatic Biology 10 : 193 - 201 . https://doi.org/10.3354/ab00278 Levan A , Fredga K , Sandberg AA ( 1964 ) Nomenclature for centromeric position on chromosomes . Hereditas 52 : 201 - 220 . https://doi.org/10.1111/j.1601- 5223 . 1964 .tb01953.x Lefevre G ( 1907 ) Artificial parthenogenesis in Thalassema mellita . Journal of Experimental Zoology 4 : 91 - 149 . https://doi.org/10.1002/jez.1400040105 Li L , Qian J ( 1993 ) A contribution to the karyotypes of two endemic species of pinus in China with a discussion on the systematic position of the genus . Acta Botanica Yunnanica 15 : 47 - 56 . Liu X , Zhang L , Zhang Z , Ma X , Liu J ( 2015 ) Transcriptional response to sulfide in the Echiuran Worm Urechis unicinctus by digital gene expression analysis . BMC Genomics 16 : 829 . https://doi.org/10.1186/s12864-015-2094-z Shi S , Wang Q , Deng Y ( 2013 ) Karyotype analysis of Phasolosma esculenta . Guangdong Agricultural Sciences 24 : 132 - 133 . Singhal RN , Dattagupta AK ( 1980 ) Cytological studies of Achaetobonellia maculata Fisher (Echiura: Bonellidae) . Caryologia 33 : 449 - 459 . https://doi.org/10.1080/00087114. 1980 .10796859 Soldi R , Ramella L , Gambi MC , Sordino P , Sella G ( 1994 ) Genome size in polychaetes: relationship with body length and life habit . M?moires du Mus?um National d'Histoire Naturelle 162 : 129 - 135 . Stebbins GL ( 1971 ) Chromosomal evolution in higher plants . Edward Arnold , London, 216 pp. Struck TH , Schult N , Kusen T , Hickman E , Bleidorn C , McHugh D , Halanych KM ( 2007 ) Annelid phylogeny and the status of Sipuncula and Echiura . BMC Evolutionary Biology 7 : 57 . https://doi.org/10.1186/ 1471 -2148-7-57 Struck TH , Paul C , Hill N , Hartmann S , Hosel C , Kube M , Lieb B , Meyer A , Tiedemann R , Purschke G , Bleidorn C ( 2011 ) Phylogenomic analyses unravel annelid evolution . Nature 471 : 95 - 98 . https://doi.org/10.1038/nature09864 Wang Q , Du Y , Lin S , Du X , Deng Y ( 2008 ) Karyotype analysis of Sipunculus nudus . Chinese Journal of Zoology 43 : 95 - 98 . Weigert A , Helm C , Meyer M , Nickel B , Arendt D , Hausdorf B , Santos SR , Halanych KM , Purschke G , Bleidorn C ( 2014 ) Illuminating the base of the annelid tree using transcriptomics . Molecular Biology and Evolution 31 : 1391 - 1401 . https://doi.org/10.1093/molbev/msu080 Wu Z , Shen X , Sun M , Ren J , Wang Y , Huang Y , Liu B ( 2009 ) Phylogenetic analyses of complete mitochondrial genome of Urechis unicinctus (Echiura) support that echiurans are derived annelids . Molecular Phylogenetics and Evolution 52 : 558 - 562 . https://doi. org/10.1016/j.ympev. 2009 . 03 .009


This is a preview of a remote PDF: https://compcytogen.pensoft.net/article/31448/download/pdf/

Zhenkui Qin, Xueyu Li, Danwen Liu, Qing Wang, Li Lu, Zhifeng Zhang. Analysis of chromosome karyotype and genome size in echiuran Urechis unicinctus Drasche, 1880 (Polychaeta, Urechidae), Comparative Cytogenetics, 2019, 75-85, DOI: doi:10.3897/CompCytogen.v13i1.31448