1
|
Zang H, Guo S, Dong S, Song Y, Li K, Fan X, Qiu J, Zheng Y, Jiang H, Wu Y, Lü Y, Chen D, Guo R. Construction of a Full-Length Transcriptome of Western Honeybee Midgut Tissue and Improved Genome Annotation. Genes (Basel) 2024; 15:728. [PMID: 38927663 PMCID: PMC11202838 DOI: 10.3390/genes15060728] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Revised: 05/22/2024] [Accepted: 05/26/2024] [Indexed: 06/28/2024] Open
Abstract
Honeybees are an indispensable pollinator in nature with pivotal ecological, economic, and scientific value. However, a full-length transcriptome for Apis mellifera, assembled with the advanced third-generation nanopore sequencing technology, has yet to be reported. Here, nanopore sequencing of the midgut tissues of uninoculated and Nosema ceranae-inoculated A. mellifera workers was conducted, and the full-length transcriptome was then constructed and annotated based on high-quality long reads. Next followed improvement of sequences and annotations of the current reference genome of A. mellifera. A total of 5,942,745 and 6,664,923 raw reads were produced from midguts of workers at 7 days post-inoculation (dpi) with N. ceranae and 10 dpi, while 7,100,161 and 6,506,665 raw reads were generated from the midguts of corresponding uninoculated workers. After strict quality control, 6,928,170, 6,353,066, 5,745,048, and 6,416,987 clean reads were obtained, with a length distribution ranging from 1 kb to 10 kb. Additionally, 16,824, 17,708, 15,744, and 18,246 full-length transcripts were respectively detected, including 28,019 nonredundant ones. Among these, 43,666, 30,945, 41,771, 26,442, and 24,532 full-length transcripts could be annotated to the Nr, KOG, eggNOG, GO, and KEGG databases, respectively. Additionally, 501 novel genes (20,326 novel transcripts) were identified for the first time, among which 401 (20,255), 193 (13,365), 414 (19,186), 228 (12,093), and 202 (11,703) were respectively annotated to each of the aforementioned five databases. The expression and sequences of three randomly selected novel transcripts were confirmed by RT-PCR and Sanger sequencing. The 5' UTR of 2082 genes, the 3' UTR of 2029 genes, and both the 5' and 3' UTRs of 730 genes were extended. Moreover, 17,345 SSRs, 14,789 complete ORFs, 1224 long non-coding RNAs (lncRNAs), and 650 transcription factors (TFs) from 37 families were detected. Findings from this work not only refine the annotation of the A. mellifera reference genome, but also provide a valuable resource and basis for relevant molecular and -omics studies.
Collapse
Affiliation(s)
- He Zang
- College of Bee Science and Biomedicine, Fujian Agriculture and Forestry University, Fuzhou 350002, China; (H.Z.); (S.G.); (S.D.); (Y.S.); (K.L.); (X.F.); (J.Q.); (Y.Z.); (D.C.)
- National & Local United Engineering Laboratory of Natural Biotoxin, Fuzhou 350002, China
- Apitherapy Research Institute of Fujian Province, Fuzhou 350002, China
| | - Sijia Guo
- College of Bee Science and Biomedicine, Fujian Agriculture and Forestry University, Fuzhou 350002, China; (H.Z.); (S.G.); (S.D.); (Y.S.); (K.L.); (X.F.); (J.Q.); (Y.Z.); (D.C.)
| | - Shunan Dong
- College of Bee Science and Biomedicine, Fujian Agriculture and Forestry University, Fuzhou 350002, China; (H.Z.); (S.G.); (S.D.); (Y.S.); (K.L.); (X.F.); (J.Q.); (Y.Z.); (D.C.)
| | - Yuxuan Song
- College of Bee Science and Biomedicine, Fujian Agriculture and Forestry University, Fuzhou 350002, China; (H.Z.); (S.G.); (S.D.); (Y.S.); (K.L.); (X.F.); (J.Q.); (Y.Z.); (D.C.)
| | - Kunze Li
- College of Bee Science and Biomedicine, Fujian Agriculture and Forestry University, Fuzhou 350002, China; (H.Z.); (S.G.); (S.D.); (Y.S.); (K.L.); (X.F.); (J.Q.); (Y.Z.); (D.C.)
| | - Xiaoxue Fan
- College of Bee Science and Biomedicine, Fujian Agriculture and Forestry University, Fuzhou 350002, China; (H.Z.); (S.G.); (S.D.); (Y.S.); (K.L.); (X.F.); (J.Q.); (Y.Z.); (D.C.)
- National & Local United Engineering Laboratory of Natural Biotoxin, Fuzhou 350002, China
- Apitherapy Research Institute of Fujian Province, Fuzhou 350002, China
| | - Jianfeng Qiu
- College of Bee Science and Biomedicine, Fujian Agriculture and Forestry University, Fuzhou 350002, China; (H.Z.); (S.G.); (S.D.); (Y.S.); (K.L.); (X.F.); (J.Q.); (Y.Z.); (D.C.)
- National & Local United Engineering Laboratory of Natural Biotoxin, Fuzhou 350002, China
- Apitherapy Research Institute of Fujian Province, Fuzhou 350002, China
| | - Yidi Zheng
- College of Bee Science and Biomedicine, Fujian Agriculture and Forestry University, Fuzhou 350002, China; (H.Z.); (S.G.); (S.D.); (Y.S.); (K.L.); (X.F.); (J.Q.); (Y.Z.); (D.C.)
| | - Haibin Jiang
- Apiculture Science Institute of Jilin Province, Jilin 132000, China; (H.J.); (Y.W.)
| | - Ying Wu
- Apiculture Science Institute of Jilin Province, Jilin 132000, China; (H.J.); (Y.W.)
| | - Yang Lü
- Mudanjiang Branch of Heilongjiang Academy of Agricultural Sciences, Mudanjiang 157000, China;
| | - Dafu Chen
- College of Bee Science and Biomedicine, Fujian Agriculture and Forestry University, Fuzhou 350002, China; (H.Z.); (S.G.); (S.D.); (Y.S.); (K.L.); (X.F.); (J.Q.); (Y.Z.); (D.C.)
- National & Local United Engineering Laboratory of Natural Biotoxin, Fuzhou 350002, China
- Apitherapy Research Institute of Fujian Province, Fuzhou 350002, China
| | - Rui Guo
- College of Bee Science and Biomedicine, Fujian Agriculture and Forestry University, Fuzhou 350002, China; (H.Z.); (S.G.); (S.D.); (Y.S.); (K.L.); (X.F.); (J.Q.); (Y.Z.); (D.C.)
- National & Local United Engineering Laboratory of Natural Biotoxin, Fuzhou 350002, China
- Apitherapy Research Institute of Fujian Province, Fuzhou 350002, China
| |
Collapse
|
2
|
Madrigal G, Minhas BF, Catchen J. Klumpy: A tool to evaluate the integrity of long-read genome assemblies and illusive sequence motifs. Mol Ecol Resour 2024:e13982. [PMID: 38800997 DOI: 10.1111/1755-0998.13982] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2024] [Accepted: 05/13/2024] [Indexed: 05/29/2024]
Abstract
The improvement and decreasing costs of third-generation sequencing technologies has widened the scope of biological questions researchers can address with de novo genome assemblies. With the increasing number of reference genomes, validating their integrity with minimal overhead is vital for establishing confident results in their applications. Here, we present Klumpy, a tool for detecting and visualizing both misassembled regions in a genome assembly and genetic elements (e.g. genes) of interest in a set of sequences. By leveraging the initial raw reads in combination with their respective genome assembly, we illustrate Klumpy's utility by investigating antifreeze glycoprotein (afgp) loci across two icefishes, by searching for a reported absent gene in the northern snakehead fish, and by scanning the reference genomes of a mudskipper and bumblebee for misassembled regions. In the two former cases, we were able to provide support for the noncanonical placement of an afgp locus in the icefishes and locate the missing snakehead gene. Furthermore, our genome scans were able identify an unmappable locus in the mudskipper reference genome and identify a putative repetitive element shared among several species of bees.
Collapse
Affiliation(s)
- Giovanni Madrigal
- Department of Evolution, Ecology, and Behavior, University of Illinois at Urbana-Champaign, Urbana, Illinois, USA
| | - Bushra Fazal Minhas
- Informatics Program, University of Illinois at Urbana-Champaign, Urbana, Illinois, USA
| | - Julian Catchen
- Department of Evolution, Ecology, and Behavior, University of Illinois at Urbana-Champaign, Urbana, Illinois, USA
- Informatics Program, University of Illinois at Urbana-Champaign, Urbana, Illinois, USA
| |
Collapse
|
3
|
Chen G, Stepanenko A, Borisjuk N. Contrasting patterns of 5S rDNA repeats in European and Asian ecotypes of greater duckweed, Spirodela polyrhiza (Lemnaceae). FRONTIERS IN PLANT SCIENCE 2024; 15:1378683. [PMID: 38711607 PMCID: PMC11070557 DOI: 10.3389/fpls.2024.1378683] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Accepted: 04/09/2024] [Indexed: 05/08/2024]
Abstract
Ribosomal DNA (rDNA) contains highly conserved, specifically organized sequences encoding ribosomal RNAs (rRNAs) separated by variable non-transcribed intergenic spacers (NTSs) and is abundant in eukaryotic genomes. These characteristics make the rDNA an informative molecular target to study genome organization, molecular evolution, and phylogenetics. In this study, we characterized the 5S rDNA repeats in the greater duckweed Spiroldela polyrhiza, a species known for its small size, rapid growth, highly conserved genome organization, and low mutation rate. Sequence analysis of at least 12 individually cloned PCR fragments containing the 5S rDNA units for each of six ecotypes that originated from Europe (Ukraine) and Asia (China) revealed two distinct types of 5S rDNA repeats containing NTSs of different lengths and nucleotide compositions. The shorter 5S rDNA repeat units had a highly homogeneous 400-bp NTS, with few ecotype- or region-specific single-nucleotide polymorphisms (SNPs). The longer 5S rDNA units had NTSs of 1056-1084 bp with characteristic intra- and inter-genomic variants due to specific SNPs and insertions/deletions of 4-15-bp DNA elements. We also detected significant variability in the ratio of short/long 5S rDNA variants between ecotypes of S. polyrhiza. The contrasting dynamics of the two types of 5S rDNA units, combined with the unusually low repeat copy number (for plants) in S. polyrhiza (46-220 copies per genome), shows that this species could serve as an excellent model for examining the mechanisms of concerted evolution and functional significance of rDNA variability.
Collapse
Affiliation(s)
- Guimin Chen
- School of Life Sciences, Huaiyin Normal University, Huai’an, China
| | - Anton Stepanenko
- School of Life Sciences, Huaiyin Normal University, Huai’an, China
- Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), Gatersleben, Germany
- Department of Molecular Genetics, Institute of Cell Biology and Genetic Engineering, Kyiv, Ukraine
| | - Nikolai Borisjuk
- School of Life Sciences, Huaiyin Normal University, Huai’an, China
| |
Collapse
|
4
|
Badger JH, Giordano R, Zimin A, Wappel R, Eskipehlivan SM, Muller S, Donthu R, Soto-Adames F, Vieira P, Zasada I, Goodwin S. Direct sequencing of insect symbionts via nanopore adaptive sampling. CURRENT OPINION IN INSECT SCIENCE 2024; 61:101135. [PMID: 37926187 DOI: 10.1016/j.cois.2023.101135] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/13/2023] [Revised: 10/16/2023] [Accepted: 10/17/2023] [Indexed: 11/07/2023]
Abstract
Insect symbionts can alter their host phenotype and their effects can range from beneficial to pathogenic. Moreover, many insects exhibit co-infections, making their study more challenging. Less than 1% of insect species have high-quality referenced genomes available and fewer still also have their symbionts sequenced. Two methods are commonly used to sequence symbionts: whole-genome sequencing to concomitantly capture the host and bacterial genomes, or isolation of the symbiont's genome before sequencing. These methods are limited when dealing with rare or poorly characterized symbionts. Long-read technology is an important tool to generate high-quality genomes as they can overcome high levels of heterozygosity, repeat content, and transposable elements that confound short-read methods. Oxford Nanopore (ONT) adaptive sampling allows a sequencing instrument to select or reject sequences in real time. We describe a method based on ONT adaptive sampling (subtractive) approach that readily permitted the sequencing of the complete genomes of mitochondria, Buchnera and its plasmids (pLeu, pTrp), and Wolbachia genomes in two aphid species, Aphis glycines and Pentalonia nigronervosa. Adaptive sampling is able to retrieve organelles such as mitochondria and symbionts that have high representation in their hosts such as Buchnera and Wolbachia, but is less successful at retrieving symbionts in low concentrations.
Collapse
Affiliation(s)
- Jonathan H Badger
- Genetics and Microbiome Core, Laboratory of Integrative Cancer Immunology, Center for Cancer Research, National Cancer Institute, Bethesda, MD, USA
| | - Rosanna Giordano
- Institute of Environment, Florida International University, Miami, FL, USA
| | - Aleksey Zimin
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Robert Wappel
- Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | | | | | - Ravikiran Donthu
- Centre for Life Sciences, Mahindra University, Bahadurpally, Hyderabad 500043, India
| | - Felipe Soto-Adames
- Florida Department of Agriculture and Consumer Services, Department of Plant Industry, Gainesville, FL 32614, USA
| | - Paulo Vieira
- USDA-ARS Agricultural Research Center, Mycology & Nematology Genetic Diversity & Biology Laboratory, Beltsville, MD, USA
| | - Inga Zasada
- USDA-ARS Horticultural Crops Research Laboratory, Corvallis, OR, USA
| | - Sara Goodwin
- Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA.
| |
Collapse
|
5
|
Duarte VDS, Porcellato D. Host DNA depletion methods and genome-centric metagenomics of bovine hindmilk microbiome. mSphere 2024; 9:e0047023. [PMID: 38054728 PMCID: PMC10826364 DOI: 10.1128/msphere.00470-23] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2023] [Accepted: 10/20/2023] [Indexed: 12/07/2023] Open
Abstract
Bovine mastitis is a multi-etiological and complex disease, resulting in serious economic consequences for dairy farmers and industry. In recent years, the microbiological evaluation of raw milk has been investigated in-depth using next-generation sequencing approaches such as metataxonomic analysis. Despite this, host DNA is a major concern in the shotgun metagenomic sequencing of microbial communities in milk samples, and it represents a big challenge. In this study, we aimed to evaluate different methods for host DNA depletion and/or microbial DNA enrichment and assess the use of PCR-based whole genome amplification in milk samples with high somatic cell count (SCC) by using short- and long-read sequencing technologies. Our results evidenced that DNA extraction performed differently in terms of host DNA removal, impacting metagenome composition and functional profiles.. Moreover, the ratio of SCC/bacteria ultimately impacts microbial DNA yield, and samples with low SCC (SCC below 100,000 cells/mL) are the most problematic. When milk samples with high SCC (SCC above 200,000 cells/mL) underwent multiple-displacement amplification (MDA), we successfully recovered high-quality metagenome-assembled genomes (MAGs), and long-read sequencing was feasible even for samples with low DNA concentration. By associating MDA and short-read sequencing, we recovered two times more MAGs than in untreated samples, and an ongoing co-infection not reported by traditional methods was detected for mastitis pathogen. Overall, this new approach will improve the detection of mastitis-associated microorganisms and make it possible to examine host-microbiome interactions in bovine mastitis.IMPORTANCENext-generation sequencing technologies have been widely used to gain new insights into the diversity of the microbial community of milk samples and dairy products for different purposes such as microbial safety, profiling of starter cultures, and host-microbiome interactions. Milk is a complex food matrix, and additionally, the presence of host nucleic acid sequences is considered a contaminant in untargeted high-throughput sequencing studies. Therefore, genomic-centric metagenomic studies of milk samples focusing on the health-disease status in dairy cattle are still scarce, which makes it difficult to evaluate the microbial ecophysiology of bovine hindmilk. This study provides an alternative method for genome-centric metagenome studies applied to hindmilk samples with high somatic cell content, which is indispensable to examining host-microbiome interactions in bovine mastitis.
Collapse
Affiliation(s)
- Vinícius da Silva Duarte
- Faculty of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences, Ås, Norway
| | - Davide Porcellato
- Faculty of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences, Ås, Norway
| |
Collapse
|
6
|
Heckenhauer J, Stewart RJ, Ríos-Touma B, Powell A, Dorji T, Frandsen PB, Pauls SU. Characterization of the primary structure of the major silk gene, h-fibroin, across caddisfly (Trichoptera) suborders. iScience 2023; 26:107253. [PMID: 37529107 PMCID: PMC10387566 DOI: 10.1016/j.isci.2023.107253] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2023] [Revised: 04/05/2023] [Accepted: 06/27/2023] [Indexed: 08/03/2023] Open
Abstract
Larvae of caddisflies (Trichoptera) produce silk to build various underwater structures allowing them to exploit a wide range of aquatic environments. The silk adheres to various substrates underwater and has high tensile strength, extensibility, and toughness and is of interest as a model for biomimetic adhesives. As a step toward understanding how the properties of underwater silk evolved in Trichoptera, we used genomic data to identify full-length sequences and characterize the primary structure of the major silk protein, h-fibroin, across the order. The h-fibroins have conserved termini and basic motif structure with high variation in repeating modules and variation in the percentage of amino acids, mainly proline. This finding might be linked to differences in mechanical properties related to the different silk usage and sets a starting point for future studies to screen and correlate amino acid motifs and other sequence features with quantifiable silk properties.
Collapse
Affiliation(s)
- Jacqueline Heckenhauer
- LOEWE Centre for Translational Biodiversity Genomics (LOEWE-TBG), Frankfurt, Hesse 60325, Germany
- Department of Terrestrial Zoology, Senckenberg Research Institute and Natural History Museum Frankfurt, Frankfurt, Hesse 60325, Germany
| | - Russell J. Stewart
- Department of Biomedical Engineering, University of Utah, Salt Lake City, UT 84112, USA
| | - Blanca Ríos-Touma
- Facultad de Ingenierías y Ciencias Aplicadas, Ingeniería Ambiental, Grupo de Investigación en Biodiversidad, Medio Ambiente y Salud (BIOMAS), Universidad de Las Américas, Quito, EC 170124, Ecuador
| | - Ashlyn Powell
- Department of Plant and Wildlife Science, Brigham Young University, Provo, UT 84602, USA
| | - Tshering Dorji
- Department of Environment and Climate Studies, Royal University of Bhutan, Punakha 13001, Bhutan
| | - Paul B. Frandsen
- LOEWE Centre for Translational Biodiversity Genomics (LOEWE-TBG), Frankfurt, Hesse 60325, Germany
- Department of Plant and Wildlife Science, Brigham Young University, Provo, UT 84602, USA
- Data Science Lab, Smithsonian Institution, Washington, DC 20560, USA
| | - Steffen U. Pauls
- LOEWE Centre for Translational Biodiversity Genomics (LOEWE-TBG), Frankfurt, Hesse 60325, Germany
- Department of Terrestrial Zoology, Senckenberg Research Institute and Natural History Museum Frankfurt, Frankfurt, Hesse 60325, Germany
- Institute for Insect Biotechnology, Justus-Liebig-University, Gießen, Hesse 35392; Germany
| |
Collapse
|
7
|
Wang Z, Kim W, Wang YW, Yakubovich E, Dong C, Trail F, Townsend JP, Yarden O. The Sordariomycetes: an expanding resource with Big Data for mining in evolutionary genomics and transcriptomics. FRONTIERS IN FUNGAL BIOLOGY 2023; 4:1214537. [PMID: 37746130 PMCID: PMC10512317 DOI: 10.3389/ffunb.2023.1214537] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/30/2023] [Accepted: 06/06/2023] [Indexed: 09/26/2023]
Abstract
Advances in genomics and transcriptomics accompanying the rapid accumulation of omics data have provided new tools that have transformed and expanded the traditional concepts of model fungi. Evolutionary genomics and transcriptomics have flourished with the use of classical and newer fungal models that facilitate the study of diverse topics encompassing fungal biology and development. Technological advances have also created the opportunity to obtain and mine large datasets. One such continuously growing dataset is that of the Sordariomycetes, which exhibit a richness of species, ecological diversity, economic importance, and a profound research history on amenable models. Currently, 3,574 species of this class have been sequenced, comprising nearly one-third of the available ascomycete genomes. Among these genomes, multiple representatives of the model genera Fusarium, Neurospora, and Trichoderma are present. In this review, we examine recently published studies and data on the Sordariomycetes that have contributed novel insights to the field of fungal evolution via integrative analyses of the genetic, pathogenic, and other biological characteristics of the fungi. Some of these studies applied ancestral state analysis of gene expression among divergent lineages to infer regulatory network models, identify key genetic elements in fungal sexual development, and investigate the regulation of conidial germination and secondary metabolism. Such multispecies investigations address challenges in the study of fungal evolutionary genomics derived from studies that are often based on limited model genomes and that primarily focus on the aspects of biology driven by knowledge drawn from a few model species. Rapidly accumulating information and expanding capabilities for systems biological analysis of Big Data are setting the stage for the expansion of the concept of model systems from unitary taxonomic species/genera to inclusive clusters of well-studied models that can facilitate both the in-depth study of specific lineages and also investigation of trait diversity across lineages. The Sordariomycetes class, in particular, offers abundant omics data and a large and active global research community. As such, the Sordariomycetes can form a core omics clade, providing a blueprint for the expansion of our knowledge of evolution at the genomic scale in the exciting era of Big Data and artificial intelligence, and serving as a reference for the future analysis of different taxonomic levels within the fungal kingdom.
Collapse
Affiliation(s)
- Zheng Wang
- Department of Biostatistics, Yale School of Public Health, New Haven, CT, United States
| | - Wonyong Kim
- Korean Lichen Research Institute, Sunchon National University, Suncheon, Republic of Korea
| | - Yen-Wen Wang
- Department of Biostatistics, Yale School of Public Health, New Haven, CT, United States
| | - Elizabeta Yakubovich
- Department of Plant Pathology and Microbiology, The Robert H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University of Jerusalem, Rehovot, Israel
| | - Caihong Dong
- Institute of Microbiology, Chinese Academy of Sciences, Beijing, China
| | - Frances Trail
- Department of Plant Biology, Michigan State University, East Lansing, MI, United States
- Department of Plant, Soil and Microbial Sciences, Michigan State University, East Lansing, MI, United States
| | - Jeffrey P. Townsend
- Department of Biostatistics, Yale School of Public Health, New Haven, CT, United States
- Department of Ecology and Evolutionary Biology, Program in Microbiology, and Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, United States
| | - Oded Yarden
- Department of Plant Pathology and Microbiology, The Robert H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University of Jerusalem, Rehovot, Israel
| |
Collapse
|