1
|
Alvarez RV, Landsman D. GTax: improving de novo transcriptome assembly by removing foreign RNA contamination. Genome Biol 2024; 25:12. [PMID: 38191464 PMCID: PMC10773103 DOI: 10.1186/s13059-023-03141-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2022] [Accepted: 12/08/2023] [Indexed: 01/10/2024] Open
Abstract
The cost and complexity of generating a complete reference genome means that many organisms lack an annotated reference. An alternative is to use a de novo reference transcriptome. This technology is cost-effective but is susceptible to off-target RNA contamination. In this manuscript, we present GTax, a taxonomy-structured database of genomic sequences that can be used with BLAST to detect and remove foreign contamination in RNA sequencing samples before assembly. In addition, we use a de novo transcriptome assembly of Solanum lycopersicum (tomato) to demonstrate that removing foreign contamination in sequencing samples reduces the number of assembled chimeric transcripts.
Collapse
Affiliation(s)
- Roberto Vera Alvarez
- Computational Biology Branch, National Center for Biotechnology Information, Intramural Research Program, National Library of Medicine, NIH, Bethesda, MD, USA
| | - David Landsman
- Computational Biology Branch, National Center for Biotechnology Information, Intramural Research Program, National Library of Medicine, NIH, Bethesda, MD, USA.
| |
Collapse
|
2
|
Shi X, Li J, Liu T, Zhao H, Leng H, Sun K, Feng J. Divergence of cochlear transcriptomics between reference‑based and reference‑free transcriptome analyses among Rhinolophus ferrumequinum populations. PLoS One 2023; 18:e0288404. [PMID: 37432940 DOI: 10.1371/journal.pone.0288404] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Accepted: 06/26/2023] [Indexed: 07/13/2023] Open
Abstract
Differences in gene expression within tissues can lead to differences in tissue function. Understanding the transcriptome of a species helps elucidate the molecular mechanisms underlying phenotypic divergence. According to the presence or absence of a reference genome of for a studied species, transcriptome analyses can be divided into reference‑based and reference‑free methods, respectively. Presently, comparisons of complete transcriptome analysis results between those two methods are still rare. In this study, we compared the cochlear transcriptome analysis results of greater horseshoe bats (Rhinolophus ferrumequinum) from three lineages in China with different acoustic phenotypes using reference‑based and reference‑free methods to explore their differences in subsequent analysis. The results gained by reference-based results had lower false-positive rates and were more accurate because differentially expressed genes among the three populations obtained by this method had greater reliability and a higher annotation rate. Some phenotype-related enrichment terms, including those related to inorganic molecules and proton transmembrane channels, were also obtained only by the reference-based method. However, the reference‑based method might have the limitation of incomplete information acquisition. Thus, we believe that a combination of reference‑free and reference‑based methods is ideal for transcriptome analyses. The results of our study provided a reference for the selection of transcriptome analysis methods in the future.
Collapse
Affiliation(s)
- Xiaoxiao Shi
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun, Jilin, China
| | - Jun Li
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun, Jilin, China
| | - Tong Liu
- Department of Life Science, Jilin Agricultural University, Changchun, Jilin, China
| | - Hanbo Zhao
- Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural, Shenzhen, China
| | - Haixia Leng
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun, Jilin, China
| | - Keping Sun
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun, Jilin, China
- Key Laboratory of Vegetation Ecology, Ministry of Education, Changchun, Jilin, China
| | - Jiang Feng
- Jilin Provincial Key Laboratory of Animal Resource Conservation and Utilization, Northeast Normal University, Changchun, Jilin, China
- Department of Life Science, Jilin Agricultural University, Changchun, Jilin, China
| |
Collapse
|
3
|
De Novo Transcriptome Dataset Generation of the Swamp Buffalo Brain and Non-Brain Tissues. BIOMED RESEARCH INTERNATIONAL 2022; 2022:4472940. [PMID: 36408285 PMCID: PMC9668446 DOI: 10.1155/2022/4472940] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/07/2022] [Accepted: 09/23/2022] [Indexed: 12/02/2022]
Abstract
The sequenced data availability opened new horizons related to buffalo genetic control of economic traits and genomic diversity. The visceral organs (brain, liver, etc.) significantly involved in energy metabolism, docility, or social interactions. We performed swamp buffalo transcriptomic profiling of 24 different tissues (brain and non-brain) to identify novel transcripts and analyzed the differentially expressed genes (DEGs) of brain vs. non-brain tissues with their functional annotation. We obtained 178.57 Gb clean transcriptomic data with GC contents 52.77%, reference genome alignment 95.36%, exonic coverage 88.49%. Totally, 26363 mRNAs transcripts including 5574 novel genes were obtained. Further, 7194 transcripts were detected as DEGs by comparing brain vs. non-brain tissues group, of which 3,999 were upregulated and 3,195 downregulated. These DEGs were functionally associated with cellular metabolic activities, signal transduction, cytoprotection, and structural and binding activities. The related functional pathways included cancer pathway, PI3k-Akt signaling, axon guidance, JAK-STAT signaling, basic cellular metabolism, thermogenesis, and oxidative phosphorylation. Our study provides an in-depth understanding of swamp buffalo transcriptomic data including DEGs potentially involved in basic cellular activities and development that helped to maintain their working capacity and social interaction with humans, and also, helpful to disclose the genetic architecture of different phenotypic traits and their gene expression regulation.
Collapse
|
4
|
De Novo Transcriptomic Characterization Enables Novel Microsatellite Identification and Marker Development in Betta splendens. Life (Basel) 2021; 11:life11080803. [PMID: 34440547 PMCID: PMC8400612 DOI: 10.3390/life11080803] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2021] [Revised: 08/05/2021] [Accepted: 08/07/2021] [Indexed: 01/10/2023] Open
Abstract
The wild populations of the commercially valuable ornamental fish species, Betta splendens, and its germplasm resources have long been threatened by habitat degradation and contamination with artificially bred fish. Because of the lack of effective marker resources, population genetics research projects are severely hampered. To generate genetic data for developing polymorphic simple sequence repeat (SSR) markers and identifying functional genes, transcriptomic analysis was performed. Illumina paired-end sequencing yielded 105,505,486 clean reads, which were then de novo assembled into 69,836 unigenes. Of these, 35,751 were annotated in the non-redundant, EuKaryotic Orthologous Group, Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes and Gene Ontology databases. A total of 12,751 SSR loci were identified from the transcripts and 7970 primer pairs were designed. One hundred primer pairs were randomly selected for PCR validation and 53 successfully generated target amplification products. Further validation demonstrated that 36% (n = 19) of the 53 amplified loci were polymorphic. These data could not only enrich the genetic information for the identification of functional genes but also effectively facilitate the development of SSR markers. Such knowledge would accelerate further studies on the genetic variation and evolution, comparative genomics, linkage mapping and molecular breeding in B. splendens.
Collapse
|
5
|
Holtze S, Gorshkova E, Braude S, Cellerino A, Dammann P, Hildebrandt TB, Hoeflich A, Hoffmann S, Koch P, Terzibasi Tozzini E, Skulachev M, Skulachev VP, Sahm A. Alternative Animal Models of Aging Research. Front Mol Biosci 2021; 8:660959. [PMID: 34079817 PMCID: PMC8166319 DOI: 10.3389/fmolb.2021.660959] [Citation(s) in RCA: 36] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2021] [Accepted: 04/08/2021] [Indexed: 12/23/2022] Open
Abstract
Most research on mechanisms of aging is being conducted in a very limited number of classical model species, i.e., laboratory mouse (Mus musculus), rat (Rattus norvegicus domestica), the common fruit fly (Drosophila melanogaster) and roundworm (Caenorhabditis elegans). The obvious advantages of using these models are access to resources such as strains with known genetic properties, high-quality genomic and transcriptomic sequencing data, versatile experimental manipulation capabilities including well-established genome editing tools, as well as extensive experience in husbandry. However, this approach may introduce interpretation biases due to the specific characteristics of the investigated species, which may lead to inappropriate, or even false, generalization. For example, it is still unclear to what extent knowledge of aging mechanisms gained in short-lived model organisms is transferable to long-lived species such as humans. In addition, other specific adaptations favoring a long and healthy life from the immense evolutionary toolbox may be entirely missed. In this review, we summarize the specific characteristics of emerging animal models that have attracted the attention of gerontologists, we provide an overview of the available data and resources related to these models, and we summarize important insights gained from them in recent years. The models presented include short-lived ones such as killifish (Nothobranchius furzeri), long-lived ones such as primates (Callithrix jacchus, Cebus imitator, Macaca mulatta), bathyergid mole-rats (Heterocephalus glaber, Fukomys spp.), bats (Myotis spp.), birds, olms (Proteus anguinus), turtles, greenland sharks, bivalves (Arctica islandica), and potentially non-aging ones such as Hydra and Planaria.
Collapse
Affiliation(s)
- Susanne Holtze
- Department of Reproduction Management, Leibniz Institute for Zoo and Wildlife Research, Berlin, Germany
| | - Ekaterina Gorshkova
- Center for Precision Genome Editing and Genetic Technologies for Biomedicine, Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, Moscow, Russia
- Faculty of Biology, Lomonosov Moscow State University, Moscow, Russia
| | - Stan Braude
- Department of Biology, Washington University in St. Louis, St. Louis, MO, United States
| | - Alessandro Cellerino
- Biology Laboratory, Scuola Normale Superiore, Pisa, Italy
- Leibniz Institute on Aging – Fritz Lipmann Institute, Jena, Germany
| | - Philip Dammann
- Department of General Zoology, Faculty of Biology, University of Duisburg-Essen, Essen, Germany
- Central Animal Laboratory, University Hospital Essen, Essen, Germany
| | - Thomas B. Hildebrandt
- Department of Reproduction Management, Leibniz Institute for Zoo and Wildlife Research, Berlin, Germany
- Faculty of Veterinary Medicine, Free University of Berlin, Berlin, Germany
| | - Andreas Hoeflich
- Division Signal Transduction, Institute for Genome Biology, Leibniz Institute for Farm Animal Biology, Dummerstorf, Germany
| | - Steve Hoffmann
- Computational Biology Group, Leibniz Institute on Aging – Fritz Lipmann Institute, Jena, Germany
| | - Philipp Koch
- Core Facility Life Science Computing, Leibniz Institute on Aging – Fritz Lipmann Institute, Jena, Germany
| | - Eva Terzibasi Tozzini
- Department of Biology and Evolution of Marine Organisms, Stazione Zoologica Anton Dohrn, Naples, Italy
| | - Maxim Skulachev
- Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Moscow, Russia
| | - Vladimir P. Skulachev
- Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Moscow, Russia
| | - Arne Sahm
- Computational Biology Group, Leibniz Institute on Aging – Fritz Lipmann Institute, Jena, Germany
| |
Collapse
|
6
|
Wang Q, Lu W, Yang J, Jiang L, Zhang Q, Kan X, Yang X. Comparative transcriptomics in three Passerida species provides insights into the evolution of avian mitochondrial complex I. COMPARATIVE BIOCHEMISTRY AND PHYSIOLOGY D-GENOMICS & PROTEOMICS 2018; 28:27-36. [DOI: 10.1016/j.cbd.2018.06.002] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/12/2017] [Revised: 06/04/2018] [Accepted: 06/13/2018] [Indexed: 02/02/2023]
|
7
|
Frias-Soler RC, Villarín Pildaín L, Hotz-Wagenblatt A, Kolibius J, Bairlein F, Wink M. De novo annotation of the transcriptome of the Northern Wheatear ( Oenanthe oenanthe). PeerJ 2018; 6:e5860. [PMID: 30498627 PMCID: PMC6251345 DOI: 10.7717/peerj.5860] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2017] [Accepted: 10/02/2018] [Indexed: 11/20/2022] Open
Abstract
We have sequenced a partial transcriptome of the Northern Wheatear (Oenanthe oenanthe), a species with one of the longest migrations on Earth. The transcriptome was constructed de novo using RNA-Seq sequence data from the pooled mRNA of six different tissues: brain, muscle, intestine, liver, adipose tissue and skin. The samples came from nine captive-bred wheatears collected at three different stages of the endogenous autumn migratory period: (1) lean birds prior the onset of migration, (2) during the fattening stage and (3) individuals at their migratory body mass plateau, when they have almost doubled their lean body mass. The sample structure used to build up the transcriptome of the Northern Wheatears concerning tissue composition and time guarantees the future survey of the regulatory genes involved in the development of the migratory phenotype. Through the pre-migratory period, birds accomplish outstanding physical and behavioural changes that involve all organ systems. Nevertheless, the molecular mechanisms through which birds synchronize and control hyperphagia, fattening, restlessness increase, immunity boosting and tuning the muscles for such endurance flight are still largely unknown. The use of RNA-Seq has emerged as a powerful tool to analyse complex traits on a broad scale, and we believe it can help to characterize the migratory phenotype of wheatears at an unprecedented level. The primary challenge to conduct quantitative transcriptomic studies in non-model species is the availability of a reference transcriptome, which we have constructed and described in this paper. The cDNA was sequenced by pyrosequencing using the Genome Sequencer Roche GS FLX System; with single paired-end reads of about 400 bp. We estimate the total number of genes at 15,640, of which 67% could be annotated using Turkey and Zebra Finch genomes, or protein sequence information from SwissProt and NCBI databases. With our study, we have made a first step towards understanding the migratory phenotype regarding gene expression of a species that has become a model to study birds long-distance migrations.
Collapse
Affiliation(s)
- Roberto Carlos Frias-Soler
- Institute of Pharmacy and Molecular Biotechnology, Heidelberg University, Heidelberg, Baden Württemberg, Germany.,Institute of Avian Research, Wilhelmshaven, Germany
| | - Lilian Villarín Pildaín
- Institute of Pharmacy and Molecular Biotechnology, Heidelberg University, Heidelberg, Baden Württemberg, Germany
| | - Agnes Hotz-Wagenblatt
- Bioinformatics Group, Core Facility Genomics and Proteomics, German Cancer Research Center, Heidelberg University, Heidelberg, Baden Württemberg, Germany
| | - Jonas Kolibius
- Institute of Pharmacy and Molecular Biotechnology, Heidelberg University, Heidelberg, Baden Württemberg, Germany
| | | | - Michael Wink
- Institute of Pharmacy and Molecular Biotechnology, Heidelberg University, Heidelberg, Baden Württemberg, Germany
| |
Collapse
|
8
|
MacManes MD. The Oyster River Protocol: a multi-assembler and kmer approach for de novo transcriptome assembly. PeerJ 2018; 6:e5428. [PMID: 30083482 PMCID: PMC6078068 DOI: 10.7717/peerj.5428] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2017] [Accepted: 07/21/2018] [Indexed: 11/24/2022] Open
Abstract
Characterizing transcriptomes in non-model organisms has resulted in a massive increase in our understanding of biological phenomena. This boon, largely made possible via high-throughput sequencing, means that studies of functional, evolutionary, and population genomics are now being done by hundreds or even thousands of labs around the world. For many, these studies begin with a de novo transcriptome assembly, which is a technically complicated process involving several discrete steps. The Oyster River Protocol (ORP), described here, implements a standardized and benchmarked set of bioinformatic processes, resulting in an assembly with enhanced qualities over other standard assembly methods. Specifically, ORP produced assemblies have higher Detonate and TransRate scores and mapping rates, which is largely a product of the fact that it leverages a multi-assembler and kmer assembly process, thereby bypassing the shortcomings of any one approach. These improvements are important, as previously unassembled transcripts are included in ORP assemblies, resulting in a significant enhancement of the power of downstream analysis. Further, as part of this study, I show that assembly quality is unrelated with the number of reads generated, above 30 million reads. Code Availability: The version controlled open-source code is available at https://github.com/macmanes-lab/Oyster_River_Protocol. Instructions for software installation and use, and other details are available at http://oyster-river-protocol.rtfd.org/.
Collapse
Affiliation(s)
- Matthew D MacManes
- Department of Molecular, Cellular and Biomedical Sciences, University of New Hampshire, Durham, NH, USA
| |
Collapse
|
9
|
Genes Integral to the Reproductive Function of Male Reproductive Tissues Drive Heterogeneity in Evolutionary Rates in Japanese Quail. G3-GENES GENOMES GENETICS 2018; 8:39-51. [PMID: 29158338 PMCID: PMC5765365 DOI: 10.1534/g3.117.300095] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Early comparative genomics studies originally uncovered a nonintuitive pattern; genes involved in reproduction appeared to evolve more rapidly than other classes of genes. Currently, the emerging consensus is that genes encoding reproductive proteins evolve under variable selective pressures, producing more heterogeneous divergence patterns than previously appreciated. Here, we investigate a facet of that heterogeneity and explore the factors that drive male reproductive tissue-based heterogeneity in evolutionary rates. In Japanese quail (Coturnix japonica), genes with enriched expression in the testes evolve much more rapidly than those enriched in the foam gland (FG), a novel gland that secretes an airy foam that males transfer to females during mating. We compared molecular evolutionary patterns among (1) genes with induced expression in breeding vs. wintering conditions for both tissues and (2) genes that encode foam proteins (FPs) vs. those with varying degrees of expression specificity in the FG. We report two major findings. First, genes upregulated in breeding condition testes evolve exceptionally rapidly, while those induced in breeding condition FGs evolve slowly. These differences hold even after correcting for hormonally-dependent gene expression and chromosomal location. Second, genes encoding FPs are extremely conserved in terms of gene identity and sequence. Together, these finding suggest that genes involved in the reproductive function of each tissue drive the marked rate of heterogeneity.
Collapse
|
10
|
Huang Y, Xiong JL, Gao XC, Sun XH. Transcriptome analysis of the Chinese giant salamander ( Andrias davidianus) using RNA-sequencing. GENOMICS DATA 2017; 14:126-131. [PMID: 29159068 PMCID: PMC5675895 DOI: 10.1016/j.gdata.2017.10.005] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/21/2017] [Revised: 10/21/2017] [Accepted: 10/21/2017] [Indexed: 12/21/2022]
Abstract
The Chinese giant salamander (Andrias davidianus) is an economically important animal on academic value. However, the genomic information of this species has been less studied. In our study, the transcripts of A. davidianus were obtained by RNA-seq to conduct a transcriptomic analysis. In total 132,912 unigenes were generated with an average length of 690 bp and N50 of 1263 bp by de novo assembly using Trinity software. Using a sequence similarity search against the nine public databases (CDD, KOG, NR, NT, PFAM, Swiss-prot, TrEMBL, GO and KEGG databases), a total of 24,049, 18,406, 36,711, 15,858, 20,500, 27,515, 36,705, 28,879 and 10,958 unigenes were annotated in databases, respectively. Of these, 6323 unigenes were annotated in all database and 39,672 unigenes were annotated in at least one database. Blasted with KEGG pathway, 10,958 unigenes were annotated, and it was divided into 343 categories according to different pathways. In addition, we also identified 29,790 SSRs. This study provided a valuable resource for understanding transcriptomic information of A. davidianus and laid a foundation for further research on functional gene cloning, genomics, genetic diversity analysis and molecular marker exploitation in A. davidianus.
Collapse
Affiliation(s)
- Yong Huang
- College of Animal Science and Technology, Henan University of Science and Technology, Luoyang 471023, China
| | - Jian Li Xiong
- College of Animal Science and Technology, Henan University of Science and Technology, Luoyang 471023, China
| | - Xiao Chan Gao
- College of Animal Science and Technology, Henan University of Science and Technology, Luoyang 471023, China
| | - Xi Hong Sun
- College of Animal Science and Technology, Henan University of Science and Technology, Luoyang 471023, China
| |
Collapse
|
11
|
Moazzzam Jazi M, Seyedi SM, Ebrahimie E, Ebrahimi M, De Moro G, Botanga C. A genome-wide transcriptome map of pistachio (Pistacia vera L.) provides novel insights into salinity-related genes and marker discovery. BMC Genomics 2017; 18:627. [PMID: 28814265 PMCID: PMC5559799 DOI: 10.1186/s12864-017-3989-7] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2017] [Accepted: 08/01/2017] [Indexed: 12/18/2022] Open
Abstract
Background Pistachio (Pistacia vera L.) is one of the most important commercial nut crops worldwide. It is a salt-tolerant and long-lived tree, with the largest cultivation area in Iran. Climate change and subsequent increased soil salt content have adversely affected the pistachio yield in recent years. However, the lack of genomic/global transcriptomic sequences on P. vera impedes comprehensive researches at the molecular level. Hence, whole transcriptome sequencing is required to gain insight into functional genes and pathways in response to salt stress. Results RNA sequencing of a pooled sample representing 24 different tissues of two pistachio cultivars with contrasting salinity tolerance under control and salt treatment by Illumina Hiseq 2000 platform resulted in 368,953,262 clean 100 bp paired-ends reads (90 Gb). Following creating several assemblies and assessing their quality from multiple perspectives, we found that using the annotation-based metrics together with the length-based parameters allows an improved assessment of the transcriptome assembly quality, compared to the solely use of the length-based parameters. The generated assembly by Trinity was adopted for functional annotation and subsequent analyses. In total, 29,119 contigs annotated against all of five public databases, including NR, UniProt, TAIR10, KOG and InterProScan. Among 279 KEGG pathways supported by our assembly, we further examined the pathways involved in the plant hormone biosynthesis and signaling as well as those to be contributed to secondary metabolite biosynthesis due to their importance under salinity stress. In total, 11,337 SSRs were also identified, which the most abundant being dinucleotide repeats. Besides, 13,097 transcripts as candidate stress-responsive genes were identified. Expression of some of these genes experimentally validated through quantitative real-time PCR (qRT-PCR) that further confirmed the accuracy of the assembly. From this analysis, the contrasting expression pattern of NCED3 and SOS1 genes were observed between salt-sensitive and salt-tolerant cultivars. Conclusion This study, as the first report on the whole transcriptome survey of P. vera, provides important resources and paves the way for functional and comparative genomic studies on this major tree to discover the salinity tolerance-related markers and stress response mechanisms for breeding of new pistachio cultivars with more salinity tolerance. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3989-7) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Maryam Moazzzam Jazi
- Plant Biotechnology Department, National Institute of Genetic Engineering and Biotechnology, Tehran, Iran
| | - Seyed Mahdi Seyedi
- Plant Biotechnology Department, National Institute of Genetic Engineering and Biotechnology, Tehran, Iran.
| | - Esmaeil Ebrahimie
- School of Medicine, The University of Adelaide, Adelaide, Australia.,Institute of Biotechnology, Shiraz University, Shiraz, Iran.,Division of Information Technology, Engineering and the Environment, School of Information Technology and Mathematical Sciences, University of South Australia, Adelaide, Australia.,School of Biological Sciences, Faculty of Science and Engineering, Flinders University, Adelaide, Australia
| | | | - Gianluca De Moro
- Center of Marine Sciences (CCMAR), University of Algarve, Faro, Portugal
| | - Christopher Botanga
- Department of Biological Sciences, Chicago State University, Chicago, IL, USA
| |
Collapse
|
12
|
Chen D, Chen F, Chen C, Chen X, Mao Y. Transcriptome analysis of three cotton pests reveals features of gene expressions in the mesophyll feeder Apolygus lucorum. SCIENCE CHINA-LIFE SCIENCES 2017; 60:826-838. [PMID: 28730342 DOI: 10.1007/s11427-017-9065-3] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/12/2017] [Accepted: 04/05/2017] [Indexed: 11/27/2022]
Abstract
The green mirid bug Apolygus lucorum is an agricultural pest that is known to cause damage to more than 150 plant species. Here, we report the transcriptomes of A. lucorum at three different developmental stages (the second and fifth instar nymphs and adults). A total of 98,236 unigenes with an average length of 1,335 nt was obtained, of which 50,640 were annotated, including those encoding digestive enzymes and cytochrome P450s. Comparisons with cotton bollworm and cotton aphid transcriptomes revealed distinct features of A. lucorum as a mesophyll feeder. The gene expression dynamics varied during development from young nymphs to adults. The high-quality transcriptome data and the gene expression dynamics reported here provide valuable data for a more comprehensive understanding of the physiology and development of mirid bugs, and for mining targets for their control.
Collapse
Affiliation(s)
- Dianyang Chen
- National Key Laboratory of Plant Molecular Genetics, CAS Center for Excellence in Molecular Plant Sciences, Institute of Plant Physiology and Ecology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai, 200032, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Fangyan Chen
- National Key Laboratory of Plant Molecular Genetics, CAS Center for Excellence in Molecular Plant Sciences, Institute of Plant Physiology and Ecology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai, 200032, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Chunyu Chen
- National Key Laboratory of Plant Molecular Genetics, CAS Center for Excellence in Molecular Plant Sciences, Institute of Plant Physiology and Ecology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai, 200032, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Xiaoya Chen
- National Key Laboratory of Plant Molecular Genetics, CAS Center for Excellence in Molecular Plant Sciences, Institute of Plant Physiology and Ecology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai, 200032, China.
| | - Yingbo Mao
- National Key Laboratory of Plant Molecular Genetics, CAS Center for Excellence in Molecular Plant Sciences, Institute of Plant Physiology and Ecology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai, 200032, China.
| |
Collapse
|
13
|
Cerveau N, Jackson DJ. Combining independent de novo assemblies optimizes the coding transcriptome for nonconventional model eukaryotic organisms. BMC Bioinformatics 2016; 17:525. [PMID: 27938328 PMCID: PMC5148890 DOI: 10.1186/s12859-016-1406-x] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2016] [Accepted: 12/01/2016] [Indexed: 12/01/2022] Open
Abstract
Background Next-generation sequencing (NGS) technologies are arguably the most revolutionary technical development to join the list of tools available to molecular biologists since PCR. For researchers working with nonconventional model organisms one major problem with the currently dominant NGS platform (Illumina) stems from the obligatory fragmentation of nucleic acid material that occurs prior to sequencing during library preparation. This step creates a significant bioinformatic challenge for accurate de novo assembly of novel transcriptome data. This challenge becomes apparent when a variety of modern assembly tools (of which there is no shortage) are applied to the same raw NGS dataset. With the same assembly parameters these tools can generate markedly different assembly outputs. Results In this study we present an approach that generates an optimized consensus de novo assembly of eukaryotic coding transcriptomes. This approach does not represent a new assembler, rather it combines the outputs of a variety of established assembly packages, and removes redundancy via a series of clustering steps. We test and validate our approach using Illumina datasets from six phylogenetically diverse eukaryotes (three metazoans, two plants and a yeast) and two simulated datasets derived from metazoan reference genome annotations. All of these datasets were assembled using three currently popular assembly packages (CLC, Trinity and IDBA-tran). In addition, we experimentally demonstrate that transcripts unique to one particular assembly package are likely to be bioinformatic artefacts. For all eight datasets our pipeline generates more concise transcriptomes that in fact possess more unique annotatable protein domains than any of the three individual assemblers we employed. Another measure of assembly completeness (using the purpose built BUSCO databases) also confirmed that our approach yields more information. Conclusions Our approach yields coding transcriptome assemblies that are more likely to be closer to biological reality than any of the three individual assembly packages we investigated. This approach (freely available as a simple perl script) will be of use to researchers working with species for which there is little or no reference data against which the assembly of a transcriptome can be performed. Electronic supplementary material The online version of this article (doi:10.1186/s12859-016-1406-x) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Nicolas Cerveau
- Department of Geobiology, Goldschmidtstr.3, Georg-August University of Göttingen, 37077, Göttingen, Germany
| | - Daniel J Jackson
- Department of Geobiology, Goldschmidtstr.3, Georg-August University of Göttingen, 37077, Göttingen, Germany.
| |
Collapse
|
14
|
Kang SW, Patnaik BB, Hwang HJ, Park SY, Chung JM, Song DK, Patnaik HH, Lee JB, Kim C, Kim S, Park HS, Park SH, Park YS, Han YS, Lee JS, Lee YS. Sequencing and de novo assembly of visceral mass transcriptome of the critically endangered land snail Satsuma myomphala: Annotation and SSR discovery. COMPARATIVE BIOCHEMISTRY AND PHYSIOLOGY D-GENOMICS & PROTEOMICS 2016; 21:77-89. [PMID: 28107688 DOI: 10.1016/j.cbd.2016.10.004] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/24/2016] [Revised: 10/24/2016] [Accepted: 10/26/2016] [Indexed: 12/12/2022]
Abstract
Satsuma myomphala is critically endangered through loss of natural habitats, predation by natural enemies, and indiscriminate collection. It is a protected species in Korea but lacks genomic resources for an understanding of varied functional processes attributable to evolutionary success under natural habitats. For assessing the genetic information of S. myomphala, we performed for the first time, de novo transcriptome sequencing and functional annotation of expressed sequences using Illumina Next-Generation Sequencing (NGS) platform and bioinformatics analysis. We identified 103,774 unigenes of which 37,959, 12,890, and 17,699 were annotated in the PANM (Protostome DB), Unigene, and COG (Clusters of Orthologous Groups) databases, respectively. In addition, 14,451 unigenes were predicted under Gene Ontology functional categories, with 4581 assigned to a single category. Furthermore, 3369 sequences with 646 having Enzyme Commission (EC) numbers were mapped to 122 pathways in the Kyoto Encyclopedia of Genes and Genomes Pathway database. The prominent protein domains included the Zinc finger (C2H2-like), Reverse Transcriptase, Thioredoxin-like fold, and RNA recognition motif domain. Many unigenes with homology to immunity, defense, and reproduction-related genes were screened in the transcriptome. We also detected 3120 putative simple sequence repeats (SSRs) encompassing dinucleotide to hexanucleotide repeat motifs from >1kb unigene sequences. A list of PCR primers of SSR loci have been identified to study the genetic polymorphisms. The transcriptome data represents a valuable resource for further investigations on the species genome structure and biology. The unigenes information and microsatellites would provide an indispensable tool for conservation of the species in natural and adaptive environments.
Collapse
Affiliation(s)
- Se Won Kang
- Department of Life Science and Biotechnology, College of Natural Sciences, Soonchunhyang University, 22 Soonchunhyangro, Shinchang-myeon, Asan, Chungcheongnam-do 31538, Republic of Korea
| | - Bharat Bhusan Patnaik
- Department of Life Science and Biotechnology, College of Natural Sciences, Soonchunhyang University, 22 Soonchunhyangro, Shinchang-myeon, Asan, Chungcheongnam-do 31538, Republic of Korea; Trident School of Biotech Sciences, Trident Academy of Creative Technology (TACT), Bhubaneswar, Odisha, 751024, India
| | - Hee-Ju Hwang
- Department of Life Science and Biotechnology, College of Natural Sciences, Soonchunhyang University, 22 Soonchunhyangro, Shinchang-myeon, Asan, Chungcheongnam-do 31538, Republic of Korea
| | - So Young Park
- Biodiversity Conservation & Change Research Division, Nakdonggang National Institute of Biological Resources, 137, Donam 2-gil, Sangju-si, Gyeongsangbuk-do, 37242, Republic of Korea
| | - Jong Min Chung
- Department of Life Science and Biotechnology, College of Natural Sciences, Soonchunhyang University, 22 Soonchunhyangro, Shinchang-myeon, Asan, Chungcheongnam-do 31538, Republic of Korea
| | - Dae Kwon Song
- Department of Life Science and Biotechnology, College of Natural Sciences, Soonchunhyang University, 22 Soonchunhyangro, Shinchang-myeon, Asan, Chungcheongnam-do 31538, Republic of Korea
| | - Hongray Howrelia Patnaik
- Department of Life Science and Biotechnology, College of Natural Sciences, Soonchunhyang University, 22 Soonchunhyangro, Shinchang-myeon, Asan, Chungcheongnam-do 31538, Republic of Korea
| | - Jae Bong Lee
- Korea Zoonosis Research Institute (KOZRI), Chonbuk National University, 820-120 Hana-ro, Iksan, Jeollabuk-do 54528, Republic of Korea
| | - Changmu Kim
- National Institute of Biological Resources, 42, Hwangyeong-ro, Seo-gu, Incheon 22689, Republic of Korea
| | - Soonok Kim
- National Institute of Biological Resources, 42, Hwangyeong-ro, Seo-gu, Incheon 22689, Republic of Korea
| | - Hong Seog Park
- Research Institute, GnC BIO Co., Ltd., 621-6 Banseok-dong, Yuseong-gu, Daejeon 34069, Republic of Korea
| | - Seung-Hwan Park
- Biological Resource Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), 181 Ipsin-gil (Shinjeong0dong), Jungeup-si, Jeollabuk-do,56212, Republic of Korea
| | - Young-Su Park
- Department of Nursing, College of Medicine, Soonchunhyang University, 22 Soonchunhyangro, Shinchang-myeon, Asan, Chungcheongnam-do 31538, Republic of Korea
| | - Yeon Soo Han
- College of Agriculture and Life Science, Chonnam National University, 77 Yongbong-ro, Buk-gu, Gwangju 61186, Republic of Korea
| | - Jun Sang Lee
- Institute of Environmental Research, Kangwon National University, 1 Kangwondaehak-gil, Chuncheon, Gangwon-do-si 243341, Republic of Korea
| | - Yong Seok Lee
- Department of Life Science and Biotechnology, College of Natural Sciences, Soonchunhyang University, 22 Soonchunhyangro, Shinchang-myeon, Asan, Chungcheongnam-do 31538, Republic of Korea.
| |
Collapse
|
15
|
Zhao LP, Bolouri H. Object-oriented regression for building predictive models with high dimensional omics data from translational studies. J Biomed Inform 2016; 60:431-45. [PMID: 26972839 PMCID: PMC5097461 DOI: 10.1016/j.jbi.2016.03.001] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2015] [Revised: 02/23/2016] [Accepted: 03/01/2016] [Indexed: 12/31/2022]
Abstract
Maturing omics technologies enable researchers to generate high dimension omics data (HDOD) routinely in translational clinical studies. In the field of oncology, The Cancer Genome Atlas (TCGA) provided funding support to researchers to generate different types of omics data on a common set of biospecimens with accompanying clinical data and has made the data available for the research community to mine. One important application, and the focus of this manuscript, is to build predictive models for prognostic outcomes based on HDOD. To complement prevailing regression-based approaches, we propose to use an object-oriented regression (OOR) methodology to identify exemplars specified by HDOD patterns and to assess their associations with prognostic outcome. Through computing patient's similarities to these exemplars, the OOR-based predictive model produces a risk estimate using a patient's HDOD. The primary advantages of OOR are twofold: reducing the penalty of high dimensionality and retaining the interpretability to clinical practitioners. To illustrate its utility, we apply OOR to gene expression data from non-small cell lung cancer patients in TCGA and build a predictive model for prognostic survivorship among stage I patients, i.e., we stratify these patients by their prognostic survival risks beyond histological classifications. Identification of these high-risk patients helps oncologists to develop effective treatment protocols and post-treatment disease management plans. Using the TCGA data, the total sample is divided into training and validation data sets. After building up a predictive model in the training set, we compute risk scores from the predictive model, and validate associations of risk scores with prognostic outcome in the validation data (P-value=0.015).
Collapse
Affiliation(s)
- Lue Ping Zhao
- Division of Public Health Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA, United States; Department of Biostatistics and Epidemiology, University of Washington School of Public Health, Seattle, WA, United States.
| | - Hamid Bolouri
- Division of Human Biology, Fred Hutchinson Cancer Research Center, Seattle, WA, United States
| |
Collapse
|
16
|
Rispe C, Legeai F, Papura D, Bretaudeau A, Hudaverdian S, Le Trionnaire G, Tagu D, Jaquiéry J, Delmotte F. De novo transcriptome assembly of the grapevine phylloxera allows identification of genes differentially expressed between leaf- and root-feeding forms. BMC Genomics 2016; 17:219. [PMID: 26968158 PMCID: PMC4787006 DOI: 10.1186/s12864-016-2530-8] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2015] [Accepted: 02/24/2016] [Indexed: 11/28/2022] Open
Abstract
Background Grapevine phylloxera, an insect related to true aphids, is a major historic pest of viticulture only controlled through the selection of resistant rootstocks or through quarantine regulations where grapevine is cultivated own-rooted. Transcriptomic data could help understand the bases of its original life-traits, including a striking case of polyphenism, with forms feeding on roots and forms feeding in leaf-galls. Comparisons with true aphids (for which complete genomes have been sequenced) should also allow to link differences in life-traits of the two groups with changes in gene repertoires or shifts in patterns of expression. Results We sequenced transcriptomes of the grapevine phylloxera (Illumina technology), choosing three life-stages (adults on roots or on leaf galls, and eggs) to cover a large catalogue of transcripts, and performed a de novo assembly. This resulted in 105,697 contigs, which were annotated: most contigs had a best blastx hit to the pea aphid (phylogenetically closest complete genome), while very few bacterial hits were recorded (except for Probionibacterium acnes). Coding sequences were predicted from this data set (17,372 sequences), revealing an extremely high AT-bias (at the third codon position). Differential expression (DE) analysis among root-feeding and gall-feeding showed that i) the root-feeding form displayed a much larger number of differentially expressed transcripts ii) root-feeding biased genes were enriched in some categories, for example cuticular proteins and genes associated with cell-cell signaling iii) leaf-galling-biased genes were enriched in genes associated with the nucleus and DNA-replication, suggesting a metabolism more oriented towards fast and active multiplication. We also identified a gene family with a very high expression level (copies totaling nearly 10 % of the reads) in the grapevine phylloxera (both in root and leaf galling forms), but usually expressed at very low levels in true aphids (except in sexual oviparous females). These transcripts thus appear to be associated with oviparity. Conclusions Our study illustrated major intraspecific changes in transcriptome profiles, related with different life-styles (and the feeding on roots versus in leaf-galls). At a different scale, we could also illustrate one major shift in expression levels associated with changes in life-traits that occurred along evolution and that respectively characterize (strictly oviparous) grapevine phylloxera and (mostly viviparous) true aphids. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2530-8) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Claude Rispe
- Present Address: BIOEPAR, INRA, Oniris, La Chantrerie, F-44307, Nantes, France. .,IGEPP, INRA, F-35653, Le Rheu cedex, France.
| | - Fabrice Legeai
- IGEPP, BIPAA, INRA, Campus Beaulieu, Rennes, France.,Institut National de Recherche en Informatique et en Automatique, Institut de Recherche en Informatique et Systèmes Aléatoires, Genscale, Campus Beaulieu, Rennes, France
| | | | - Anthony Bretaudeau
- IGEPP, BIPAA, INRA, Campus Beaulieu, Rennes, France.,Institut National de Recherche en Informatique et en Automatique, Institut de Recherche en Informatique et Systèmes Aléatoires, Genouest, Campus Beaulieu, Rennes, France
| | | | | | - Denis Tagu
- IGEPP, INRA, F-35653, Le Rheu cedex, France
| | - Julie Jaquiéry
- IGEPP, INRA, F-35653, Le Rheu cedex, France.,Present address: University of Rennes 1, UMR CNRS 6553 EcoBio, 35042, Rennes, France
| | | |
Collapse
|
17
|
Marchant A, Mougel F, Mendonça V, Quartier M, Jacquin-Joly E, da Rosa JA, Petit E, Harry M. Comparing de novo and reference-based transcriptome assembly strategies by applying them to the blood-sucking bug Rhodnius prolixus. INSECT BIOCHEMISTRY AND MOLECULAR BIOLOGY 2016; 69:25-33. [PMID: 26005117 DOI: 10.1016/j.ibmb.2015.05.009] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/21/2014] [Revised: 04/08/2015] [Accepted: 05/15/2015] [Indexed: 06/04/2023]
Abstract
High Throughput Sequencing capabilities have made the process of assembling a transcriptome easier, whether or not there is a reference genome. But the quality of a transcriptome assembly must be good enough to capture the most comprehensive catalog of transcripts and their variations, and to carry out further experiments on transcriptomics. There is currently no consensus on which of the many sequencing technologies and assembly tools are the most effective. Many non-model organisms lack a reference genome to guide the transcriptome assembly. One question, therefore, is whether or not a reference-based genome assembly gives better results than de novo assembly. The blood-sucking insect Rhodnius prolixus-a vector for Chagas disease-has a reference genome. It is therefore a good model on which to compare reference-based and de novo transcriptome assemblies. In this study, we compared de novo and reference-based genome assembly strategies using three datasets (454, Illumina, 454 combined with Illumina) and various assembly software. We developed criteria to compare the resulting assemblies: the size distribution and number of transcripts, the proportion of potentially chimeric transcripts, how complete the assembly was (completeness evaluated both through CEGMA software and R. prolixus proteome fraction retrieved). Moreover, we looked for the presence of two chemosensory gene families (Odorant-Binding Proteins and Chemosensory Proteins) to validate the assembly quality. The reference-based assemblies after genome annotation were clearly better than those generated using de novo strategies alone. Reference-based strategies revealed new transcripts, including new isoforms unpredicted by automatic genome annotation. However, a combination of both de novo and reference-based strategies gave the best result, and allowed us to assemble fragmented transcripts.
Collapse
Affiliation(s)
- A Marchant
- UMR EGCE (Laboratoire Evolution, Génomes, Comportement, Ecologie), Univ. Paris-Sud, CNRS, IRD, IDEEV, Univ. Paris-Saclay, Avenue de la Terrasse, Bâtiment 13, BP1 - 91198 Gif-sur-Yvette, France.
| | - F Mougel
- UMR EGCE (Laboratoire Evolution, Génomes, Comportement, Ecologie), Univ. Paris-Sud, CNRS, IRD, IDEEV, Univ. Paris-Saclay, Avenue de la Terrasse, Bâtiment 13, BP1 - 91198 Gif-sur-Yvette, France; Université Paris Sud, Orsay, France
| | - V Mendonça
- UMR EGCE (Laboratoire Evolution, Génomes, Comportement, Ecologie), Univ. Paris-Sud, CNRS, IRD, IDEEV, Univ. Paris-Saclay, Avenue de la Terrasse, Bâtiment 13, BP1 - 91198 Gif-sur-Yvette, France; Faculdade de Ciências Farmacêuticas, Universidade Estadual Paulista (UNESP), Araraquara, Brazil
| | - M Quartier
- UMR EGCE (Laboratoire Evolution, Génomes, Comportement, Ecologie), Univ. Paris-Sud, CNRS, IRD, IDEEV, Univ. Paris-Saclay, Avenue de la Terrasse, Bâtiment 13, BP1 - 91198 Gif-sur-Yvette, France; University of Neuchâtel, Institute of Biology, Neuchâtel, Switzerland
| | - E Jacquin-Joly
- INRA, UMR 1392, Institut d'Ecologie et des Sciences de l'Environnement de Paris (iEES-Paris), Versailles, France
| | - J A da Rosa
- Faculdade de Ciências Farmacêuticas, Universidade Estadual Paulista (UNESP), Araraquara, Brazil
| | - E Petit
- UMR EGCE (Laboratoire Evolution, Génomes, Comportement, Ecologie), Univ. Paris-Sud, CNRS, IRD, IDEEV, Univ. Paris-Saclay, Avenue de la Terrasse, Bâtiment 13, BP1 - 91198 Gif-sur-Yvette, France; Université Paris Sud, Orsay, France
| | - M Harry
- UMR EGCE (Laboratoire Evolution, Génomes, Comportement, Ecologie), Univ. Paris-Sud, CNRS, IRD, IDEEV, Univ. Paris-Saclay, Avenue de la Terrasse, Bâtiment 13, BP1 - 91198 Gif-sur-Yvette, France; Université Paris Sud, Orsay, France.
| |
Collapse
|
18
|
Deng T, Pang C, Lu X, Zhu P, Duan A, Tan Z, Huang J, Li H, Chen M, Liang X. De Novo Transcriptome Assembly of the Chinese Swamp Buffalo by RNA Sequencing and SSR Marker Discovery. PLoS One 2016; 11:e0147132. [PMID: 26766209 PMCID: PMC4713091 DOI: 10.1371/journal.pone.0147132] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2015] [Accepted: 12/29/2015] [Indexed: 01/11/2023] Open
Abstract
The Chinese swamp buffalo (Bubalis bubalis) is vital to the lives of small farmers and has tremendous economic importance. However, a lack of genomic information has hampered research on augmenting marker assisted breeding programs in this species. Thus, a high-throughput transcriptomic sequencing of B. bubalis was conducted to generate transcriptomic sequence dataset for gene discovery and molecular marker development. Illumina paired-end sequencing generated a total of 54,109,173 raw reads. After trimming, de novo assembly was performed, which yielded 86,017 unigenes, with an average length of 972.41 bp, an N50 of 1,505 bp, and an average GC content of 49.92%. A total of 62,337 unigenes were successfully annotated. Among the annotated unigenes, 27,025 (43.35%) and 23,232 (37.27%) unigenes showed significant similarity to known proteins in NCBI non-redundant protein and Swiss-Prot databases (E-value < 1.0E-5), respectively. Of these annotated unigenes, 14,439 and 15,813 unigenes were assigned to the Gene Ontology (GO) categories and EuKaryotic Ortholog Group (KOG) cluster, respectively. In addition, a total of 14,167 unigenes were assigned to 331 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Furthermore, 17,401 simple sequence repeats (SSRs) were identified as potential molecular markers. One hundred and fifteen primer pairs were randomly selected for amplification to detect polymorphisms. The results revealed that 110 primer pairs (95.65%) yielded PCR amplicons and 69 primer pairs (60.00%) presented polymorphisms in 35 individual buffaloes. A phylogenetic analysis showed that the five swamp buffalo populations were clustered together, whereas two river buffalo breeds clustered separately. In the present study, the Illumina RNA-seq technology was utilized to perform transcriptome analysis and SSR marker discovery in the swamp buffalo without using a reference genome. Our findings will enrich the current SSR markers resources and help spearhead molecular genetic research studies on the swamp buffalo.
Collapse
Affiliation(s)
- Tingxian Deng
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Chunying Pang
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Xingrong Lu
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Peng Zhu
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Anqin Duan
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Zhengzhun Tan
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Jian Huang
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Hui Li
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Mingtan Chen
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
| | - Xianwei Liang
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction technology, Ministry of Agriculture, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, Guangxi, P. R. China
- * E-mail:
| |
Collapse
|
19
|
Matsuzaki M, Mizushima S, Hiyama G, Hirohashi N, Shiba K, Inaba K, Suzuki T, Dohra H, Ohnishi T, Sato Y, Kohsaka T, Ichikawa Y, Atsumi Y, Yoshimura T, Sasanami T. Lactic acid is a sperm motility inactivation factor in the sperm storage tubules. Sci Rep 2015; 5:17643. [PMID: 26619826 PMCID: PMC4664960 DOI: 10.1038/srep17643] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2015] [Accepted: 11/03/2015] [Indexed: 12/23/2022] Open
Abstract
Although successful fertilization depends on timely encounters between sperm and egg, the decoupling of mating and fertilization often confers reproductive advantages to internally fertilizing animals. In several vertebrate groups, postcopulatory sperm viability is prolonged by storage in specialized organs within the female reproductive tract. In birds, ejaculated sperm can be stored in a quiescent state within oviductal sperm storage tubules (SSTs), thereby retaining fertilizability for up to 15 weeks at body temperature (41 °C); however, the mechanism by which motile sperm become quiescent within SSTs is unknown. Here, we show that low oxygen and high lactic acid concentrations are established in quail SSTs. Flagellar quiescence was induced by lactic acid in the concentration range found in SSTs through flagellar dynein ATPase inactivation following cytoplasmic acidification (<pH 6.0). The long-term preservation of sperm morphology under hypoxic and high temperature conditions indicates that a combination of these factors enables sperm cells to survive during the ovulation cycles. Our findings suggested a novel physiological role for lactic acid in promoting sperm quiescence in SSTs and opened up a new opportunity for technological improvement in prolonging sperm longevity at ambient or body temperature.
Collapse
Affiliation(s)
- Mei Matsuzaki
- Department of Applied Biological Chemistry, Faculty of Agriculture, Shizuoka University, 836 Ohya, Shizuoka, Shizuoka 422-8529, Japan.,United Graduate School of Agricultural Science, Gifu University, 1-1 Yanagido, Gifu 501-1193, Japan
| | - Shusei Mizushima
- Department of Applied Biological Chemistry, Faculty of Agriculture, Shizuoka University, 836 Ohya, Shizuoka, Shizuoka 422-8529, Japan
| | - Gen Hiyama
- Department of Applied Biological Chemistry, Faculty of Agriculture, Shizuoka University, 836 Ohya, Shizuoka, Shizuoka 422-8529, Japan
| | - Noritaka Hirohashi
- Oki Marine Biological Station, Education and Research Center for Biological Resources, Faculty of Life and Environmental Science, Shimane University, 194 Kamo, Okinoshima-cho, Oki, Shimane 685-0024, Japan
| | - Kogiku Shiba
- Shimoda Marine Research Center, University of Tsukuba, 5-10-1 Shimoda, Shizuoka 415-0025, Japan
| | - Kazuo Inaba
- Shimoda Marine Research Center, University of Tsukuba, 5-10-1 Shimoda, Shizuoka 415-0025, Japan
| | - Tomohiro Suzuki
- Research Institute of Green Science and Technology, Shizuoka University, 836 Ohya, Shizuoka, Shizuoka 422-8529, Japan
| | - Hideo Dohra
- Research Institute of Green Science and Technology, Shizuoka University, 836 Ohya, Shizuoka, Shizuoka 422-8529, Japan
| | - Toshiyuki Ohnishi
- Department of Applied Biological Chemistry, Faculty of Agriculture, Shizuoka University, 836 Ohya, Shizuoka, Shizuoka 422-8529, Japan.,Research Institute of Green Science and Technology, Shizuoka University, 836 Ohya, Shizuoka, Shizuoka 422-8529, Japan
| | - Yoshikatsu Sato
- Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan
| | - Tetsuya Kohsaka
- Department of Applied Biological Chemistry, Faculty of Agriculture, Shizuoka University, 836 Ohya, Shizuoka, Shizuoka 422-8529, Japan
| | - Yoshinobu Ichikawa
- Department of Applied Biological Chemistry, Faculty of Agriculture, Shizuoka University, 836 Ohya, Shizuoka, Shizuoka 422-8529, Japan
| | - Yusuke Atsumi
- Graduate School of Bioagricultural Sciences, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan
| | - Takashi Yoshimura
- Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan.,Graduate School of Bioagricultural Sciences, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan
| | - Tomohiro Sasanami
- Department of Applied Biological Chemistry, Faculty of Agriculture, Shizuoka University, 836 Ohya, Shizuoka, Shizuoka 422-8529, Japan
| |
Collapse
|
20
|
Sharpe ML, Dearden PK, Gimenez G, Krause KL. Comparative RNA seq analysis of the New Zealand glowworm Arachnocampa luminosa reveals bioluminescence-related genes. BMC Genomics 2015; 16:825. [PMID: 26486607 PMCID: PMC4617951 DOI: 10.1186/s12864-015-2006-2] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2015] [Accepted: 10/03/2015] [Indexed: 01/12/2023] Open
Abstract
Background The New Zealand glowworm is the larva of a carnivorous fungus gnat that produces bioluminescence to attract prey. The bioluminescent system of the glowworm is evolutionarily distinct from other well-characterised systems, especially that of the fireflies, and the molecules involved have not yet been identified. We have used high throughput sequencing technology to produce a transcriptome for the glowworm and identify transcripts encoding proteins that are likely to be involved in glowworm bioluminescence. Results Here we report the sequencing and annotation of the first transcriptome of the glowworm, and a differential analysis of expression from the glowworm light organ compared with non-light organ tissue. The analysis identified six transcripts encoding proteins that are potentially involved in glowworm bioluminescence. Three of these proteins are members of the ANL superfamily of adenylating enzymes, with similar amino acid sequences to that of the luciferase enzyme found in fireflies (31 to 37 % identical), and are candidate luciferases for the glowworm bioluminescent system. The remaining three transcripts encode putative aminoacylase, phosphatidylethanolamine-binding and glutathione S-transferase proteins. Conclusions This research provides a basis for further biochemical studies into how the glowworm produces light, and a source of genetic information to aid future ecological and evolutionary studies of the glowworm. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2006-2) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Miriam L Sharpe
- Department of Biochemistry, University of Otago, Dunedin, New Zealand.
| | - Peter K Dearden
- Department of Biochemistry, University of Otago, Dunedin, New Zealand.
| | - Gregory Gimenez
- Otago Genomics & Bioinformatics Facility, University of Otago, Dunedin, New Zealand.
| | - Kurt L Krause
- Department of Biochemistry, University of Otago, Dunedin, New Zealand.
| |
Collapse
|
21
|
Morandin C, Dhaygude K, Paviala J, Trontti K, Wheat C, Helanterä H. Caste-biases in gene expression are specific to developmental stage in the ant Formica exsecta. J Evol Biol 2015; 28:1705-18. [PMID: 26172873 DOI: 10.1111/jeb.12691] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2015] [Revised: 06/30/2015] [Accepted: 06/30/2015] [Indexed: 02/02/2023]
Abstract
Understanding how a single genome creates and maintains distinct phenotypes is a central goal in evolutionary biology. Social insects are a striking example of co-opted genetic backgrounds giving rise to dramatically different phenotypes, such as queen and worker castes. A conserved set of molecular pathways, previously envisioned as a set of 'toolkit' genes, has been hypothesized to underlie queen and worker phenotypes in independently evolved social insect lineages. Here, we investigated the toolkit from a developmental point of view, using RNA-Seq to compare caste-biased gene expression patterns across three life stages (pupae, emerging adult and old adult) and two female castes (queens and workers) in the ant Formica exsecta. We found that the number of genes with caste-biased expression increases dramatically from pupal to old adult stages. This result suggests that phenotypic differences between queens and workers at the pupal stage may derive from a relatively low number of caste-biased genes, compared to higher number of genes required to maintain caste differences at the adult stage. Gene expression patterns were more similar among castes within developmental stages than within castes despite the extensive phenotypic differences between queens and workers. Caste-biased expression was highly variable among life stages at the level of single genes, but more consistent when gene functions (gene ontology terms) were investigated. Finally, we found that a large part of putative toolkit genes were caste-biased at least in some life stages in F. exsecta, and the caste-biases, but not their direction, were more often shared between F. exsecta and other ant species than between F. exsecta and bees. Our results indicate that gene expression should be examined across several developmental stages to fully reveal the genetic basis of polyphenisms.
Collapse
Affiliation(s)
- C Morandin
- Centre of Excellence in Biological Interactions, Department of Biosciences, Helsinki University, Helsinki, Finland.,Tvärminne Zoological Station, University of Helsinki, Hanko, Finland
| | - K Dhaygude
- Centre of Excellence in Biological Interactions, Department of Biosciences, Helsinki University, Helsinki, Finland
| | - J Paviala
- Centre of Excellence in Biological Interactions, Department of Biosciences, Helsinki University, Helsinki, Finland
| | - K Trontti
- Centre of Excellence in Biological Interactions, Department of Biosciences, Helsinki University, Helsinki, Finland
| | - C Wheat
- Department of Zoology, Population Genetics, Stockholm University, Stockholm, Sweden
| | - H Helanterä
- Centre of Excellence in Biological Interactions, Department of Biosciences, Helsinki University, Helsinki, Finland.,Tvärminne Zoological Station, University of Helsinki, Hanko, Finland
| |
Collapse
|
22
|
Pelin A, Selman M, Aris-Brosou S, Farinelli L, Corradi N. Genome analyses suggest the presence of polyploidy and recent human-driven expansions in eight global populations of the honeybee pathogen Nosema ceranae. Environ Microbiol 2015; 17:4443-58. [PMID: 25914091 DOI: 10.1111/1462-2920.12883] [Citation(s) in RCA: 54] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2015] [Revised: 04/13/2015] [Accepted: 04/15/2015] [Indexed: 12/23/2022]
Abstract
Nosema ceranae is a microsporidian pathogen whose infections have been associated with recent global declines in the populations of western honeybees (Apis mellifera). Despite the outstanding economic and ecological threat that N. ceranae may represent for honeybees worldwide, many aspects of its biology, including its mode of reproduction, propagation and ploidy, are either very unclear or unknown. In the present study, we set to gain knowledge in these biological aspects by re-sequencing the genome of eight isolates (i.e. a population of spores isolated from one single beehive) of this species harvested from eight geographically distant beehives, and by investigating their level of polymorphism. Consistent with previous analyses performed using single gene sequences, our analyses uncovered the presence of very high genetic diversity within each isolate, but also very little hive-specific polymorphism. Surprisingly, the nature, location and distribution of this genetic variation suggest that beehives around the globe are infected by a population of N. ceranae cells that may be polyploid (4n or more), and possibly clonal. Lastly, phylogenetic analyses based on genome-wide single-nucleotide polymorphism data extracted from these parasites and mitochondrial sequences from their hosts all failed to support the current geographical structure of our isolates.
Collapse
Affiliation(s)
- Adrian Pelin
- Canadian Institute for Advanced Research, Department of Biology; University of Ottawa, Ottawa, ON, Canada
| | - Mohammed Selman
- Canadian Institute for Advanced Research, Department of Biology; University of Ottawa, Ottawa, ON, Canada
| | - Stéphane Aris-Brosou
- Departments of Biology and of Mathematics & Statistics, University of Ottawa, Ottawa, ON, Canada
| | - Laurent Farinelli
- FASTERIS S.A., Ch. du Pont-du-Centenaire 109, P.O. Box 28, Plan-les-Ouates, CH-1228, Geneva, Switzerland
| | - Nicolas Corradi
- Canadian Institute for Advanced Research, Department of Biology; University of Ottawa, Ottawa, ON, Canada
| |
Collapse
|
23
|
Powell D, Knibb W, Remilton C, Elizur A. De-novo transcriptome analysis of the banana shrimp (Fenneropenaeus merguiensis) and identification of genes associated with reproduction and development. Mar Genomics 2015; 22:71-8. [PMID: 25936497 DOI: 10.1016/j.margen.2015.04.006] [Citation(s) in RCA: 48] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2015] [Revised: 04/20/2015] [Accepted: 04/20/2015] [Indexed: 01/07/2023]
Abstract
The banana shrimp Fenneropenaeus merguiensis is a commercially important marine crustacean for world aquaculture and fisheries. Despite this, limited genetic information is available for it and many other penaeid shrimp species. Here we present the first in-depth analysis of the transcriptional content of 8 different tissues from the banana shrimp using RNA-Seq technologies. A total of over 1 million single-end and over 49 million paired-end reads were obtained from Roche 454FLX and illumina sequencing platforms, respectively, resulting in an assembly of 124,631 transcripts with an N50 of 1,332 and mean length of 514 nt. A total of 59,179 putative protein sequences obtained from the assembled transcripts were annotated using public protein sequence databases and assigned 20,430 BLAST hits, 16,866 GO terms and 13,304 KOG categories. Further analysis revealed a rich set of transcript sequences exhibiting homology with genes associated with reproduction, sex determination and development and distinguished the tissues responsible for this expression. This report adds a substantial contribution to the sequence data currently available for F. merguiensis, providing valuable resources for further research.
Collapse
Affiliation(s)
- Daniel Powell
- Faculty of Science, Health, Education and Engineering, University of the Sunshine Coast, Maroochydore, QLD 4558, Australia.
| | - Wayne Knibb
- Faculty of Science, Health, Education and Engineering, University of the Sunshine Coast, Maroochydore, QLD 4558, Australia.
| | | | - Abigail Elizur
- Faculty of Science, Health, Education and Engineering, University of the Sunshine Coast, Maroochydore, QLD 4558, Australia.
| |
Collapse
|
24
|
In search of pathogens: transcriptome-based identification of viral sequences from the pine processionary moth (Thaumetopoea pityocampa). Viruses 2015; 7:456-79. [PMID: 25626148 PMCID: PMC4353898 DOI: 10.3390/v7020456] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2014] [Revised: 12/29/2014] [Accepted: 01/13/2015] [Indexed: 01/06/2023] Open
Abstract
Thaumetopoea pityocampa (pine processionary moth) is one of the most important pine pests in the forests of Mediterranean countries, Central Europe, the Middle East and North Africa. Apart from causing significant damage to pinewoods, T. pityocampa occurrence is also an issue for public and animal health, as it is responsible for dermatological reactions in humans and animals by contact with its irritating hairs. High throughput sequencing technologies have allowed the fast and cost-effective generation of genetic information of interest to understand different biological aspects of non-model organisms as well as the identification of potential pathogens. Using these technologies, we have obtained and characterized the transcriptome of T. pityocampa larvae collected in 12 different geographical locations in Turkey. cDNA libraries for Illumina sequencing were prepared from four larval tissues, head, gut, fat body and integument. By pooling the sequences from Illumina platform with those previously published using the Roche 454-FLX and Sanger methods we generated the largest reference transcriptome of T. pityocampa. In addition, this study has also allowed identification of possible viral pathogens with potential application in future biocontrol strategies.
Collapse
|
25
|
Koning-Boucoiran CFS, Esselink GD, Vukosavljev M, van 't Westende WPC, Gitonga VW, Krens FA, Voorrips RE, van de Weg WE, Schulz D, Debener T, Maliepaard C, Arens P, Smulders MJM. Using RNA-Seq to assemble a rose transcriptome with more than 13,000 full-length expressed genes and to develop the WagRhSNP 68k Axiom SNP array for rose (Rosa L.). FRONTIERS IN PLANT SCIENCE 2015; 6:249. [PMID: 25954285 PMCID: PMC4404716 DOI: 10.3389/fpls.2015.00249] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/28/2015] [Accepted: 03/27/2015] [Indexed: 05/18/2023]
Abstract
In order to develop a versatile and large SNP array for rose, we set out to mine ESTs from diverse sets of rose germplasm. For this RNA-Seq libraries containing about 700 million reads were generated from tetraploid cut and garden roses using Illumina paired-end sequencing, and from diploid Rosa multiflora using 454 sequencing. Separate de novo assemblies were performed in order to identify single nucleotide polymorphisms (SNPs) within and between rose varieties. SNPs among tetraploid roses were selected for constructing a genotyping array that can be employed for genetic mapping and marker-trait association discovery in breeding programs based on tetraploid germplasm, both from cut roses and from garden roses. In total 68,893 SNPs were included on the WagRhSNP Axiom array. Next, an orthology-guided assembly was performed for the construction of a non-redundant rose transcriptome database. A total of 21,740 transcripts had significant hits with orthologous genes in the strawberry (Fragaria vesca L.) genome. Of these 13,390 appeared to contain the full-length coding regions. This newly established transcriptome resource adds considerably to the currently available sequence resources for the Rosaceae family in general and the genus Rosa in particular.
Collapse
Affiliation(s)
| | - G. Danny Esselink
- Wageningen UR Plant Breeding, Wageningen University and Research CentreWageningen, Netherlands
| | - Mirjana Vukosavljev
- Wageningen UR Plant Breeding, Wageningen University and Research CentreWageningen, Netherlands
| | | | - Virginia W. Gitonga
- Wageningen UR Plant Breeding, Wageningen University and Research CentreWageningen, Netherlands
| | - Frans A. Krens
- Wageningen UR Plant Breeding, Wageningen University and Research CentreWageningen, Netherlands
| | - Roeland E. Voorrips
- Wageningen UR Plant Breeding, Wageningen University and Research CentreWageningen, Netherlands
| | - W. Eric van de Weg
- Wageningen UR Plant Breeding, Wageningen University and Research CentreWageningen, Netherlands
| | - Dietmar Schulz
- Abteilung Molekulare Pflanzenzüchtung, Institute for Plant Genetics, Leibnitz University HannoverHannover, Germany
| | - Thomas Debener
- Abteilung Molekulare Pflanzenzüchtung, Institute for Plant Genetics, Leibnitz University HannoverHannover, Germany
| | - Chris Maliepaard
- Wageningen UR Plant Breeding, Wageningen University and Research CentreWageningen, Netherlands
| | - Paul Arens
- Wageningen UR Plant Breeding, Wageningen University and Research CentreWageningen, Netherlands
| | - Marinus J. M. Smulders
- Wageningen UR Plant Breeding, Wageningen University and Research CentreWageningen, Netherlands
- *Correspondence: Marinus J. M. Smulders, Wageningen UR Plant Breeding, Wageningen University and Research Centre, PO Box 386, NL-6708 PB Wageningen, Netherlands
| |
Collapse
|