1
|
González-Prendes R, Pena RN, Richart C, Nadal J, Ros-Freixedes R. Long-read de novo assembly of the red-legged partridge (Alectoris rufa) genome. Sci Data 2024; 11:908. [PMID: 39191744 PMCID: PMC11349902 DOI: 10.1038/s41597-024-03659-5] [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: 01/30/2024] [Accepted: 07/18/2024] [Indexed: 08/29/2024] Open
Abstract
The red-legged partridge (Alectoris rufa) is a popular game bird species that is in decline in several regions of southwestern Europe. The introduction of farm-reared individuals of a distinct genetic make-up in hunting reserves can result in genetic swamping of wild populations. Here we present a de novo genome assembly for the red-legged partridge based on long-read sequencing technology. The assembled genome size is 1.14 Gb, with scaffold N50 of 37.6 Mb and contig N50 of 29.5 Mb. Our genome is highly contiguous and contains 97.06% of complete avian core genes. Overall, the quality of this genome assembly is equivalent to those available for other close relatives such as the Japanese quail or the chicken. This genome assembly will contribute to the understanding of genetic dynamics of wild populations of red-legged partridges with releases of farm-reared reinforcements and to appropriate management decisions of such populations.
Collapse
Affiliation(s)
- Rayner González-Prendes
- Animal Breeding and Genomics, Wageningen University & Research, 6708PB, Wageningen, The Netherlands
| | - Ramona Natacha Pena
- Departament de Ciència Animal, Universitat de Lleida, Lleida, Spain
- Agrotecnio-CERCA Center, Lleida, Spain
| | - Cristóbal Richart
- Departament de Medicina i Cirurgia, Universitat Rovira i Virgili, Tarragona, Spain
| | - Jesús Nadal
- Departament de Ciència Animal, Universitat de Lleida, Lleida, Spain.
| | - Roger Ros-Freixedes
- Departament de Ciència Animal, Universitat de Lleida, Lleida, Spain.
- Agrotecnio-CERCA Center, Lleida, Spain.
| |
Collapse
|
2
|
Bhowmik O, Rahman T, Kalyanaraman A. Maptcha: an efficient parallel workflow for hybrid genome scaffolding. BMC Bioinformatics 2024; 25:263. [PMID: 39118013 PMCID: PMC11313021 DOI: 10.1186/s12859-024-05878-4] [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: 03/01/2024] [Accepted: 07/22/2024] [Indexed: 08/10/2024] Open
Abstract
BACKGROUND Genome assembly, which involves reconstructing a target genome, relies on scaffolding methods to organize and link partially assembled fragments. The rapid evolution of long read sequencing technologies toward more accurate long reads, coupled with the continued use of short read technologies, has created a unique need for hybrid assembly workflows. The construction of accurate genomic scaffolds in hybrid workflows is complicated due to scale, sequencing technology diversity (e.g., short vs. long reads, contigs or partial assemblies), and repetitive regions within a target genome. RESULTS In this paper, we present a new parallel workflow for hybrid genome scaffolding that would allow combining pre-constructed partial assemblies with newly sequenced long reads toward an improved assembly. More specifically, the workflow, called Maptcha, is aimed at generating long scaffolds of a target genome, from two sets of input sequences-an already constructed partial assembly of contigs, and a set of newly sequenced long reads. Our scaffolding approach internally uses an alignment-free mapping step to build a ⟨ contig,contig ⟩ graph using long reads as linking information. Subsequently, this graph is used to generate scaffolds. We present and evaluate a graph-theoretic "wiring" heuristic to perform this scaffolding step. To enable efficient workload management in a parallel setting, we use a batching technique that partitions the scaffolding tasks so that the more expensive alignment-based assembly step at the end can be efficiently parallelized. This step also allows the use of any standalone assembler for generating the final scaffolds. CONCLUSIONS Our experiments with Maptcha on a variety of input genomes, and comparison against two state-of-the-art hybrid scaffolders demonstrate that Maptcha is able to generate longer and more accurate scaffolds substantially faster. In almost all cases, the scaffolds produced by Maptcha are at least an order of magnitude longer (in some cases two orders) than the scaffolds produced by state-of-the-art tools. Maptcha runs significantly faster too, reducing time-to-solution from hours to minutes for most input cases. We also performed a coverage experiment by varying the sequencing coverage depth for long reads, which demonstrated the potential of Maptcha to generate significantly longer scaffolds in low coverage settings ( 1 × - 10 × ).
Collapse
Affiliation(s)
- Oieswarya Bhowmik
- School of Electrical Engineering and Computer Science, Washington State University, Pullman, WA, 99164, USA.
| | - Tazin Rahman
- School of Electrical Engineering and Computer Science, Washington State University, Pullman, WA, 99164, USA
| | - Ananth Kalyanaraman
- School of Electrical Engineering and Computer Science, Washington State University, Pullman, WA, 99164, USA
| |
Collapse
|
3
|
Fava S, Sollitto M, Racaku M, Iannucci A, Benazzo A, Ancona L, Gratton P, Florian F, Pallavicini A, Ciofi C, Cesaroni D, Gerdol M, Sbordoni V, Bertorelle G, Trucchi E. Chromosome-Level Reference Genome of the Ponza Grayling (Hipparchia sbordonii), an Italian Endemic and Endangered Butterfly. Genome Biol Evol 2024; 16:evae136. [PMID: 39023104 PMCID: PMC11255612 DOI: 10.1093/gbe/evae136] [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] [Accepted: 06/13/2024] [Indexed: 07/20/2024] Open
Abstract
Islands are crucial evolutionary hotspots, providing unique opportunities for differentiation of novel biodiversity and long-term segregation of endemic species. Islands are also fragile ecosystems, where biodiversity is more exposed to environmental and anthropogenic pressures than on continents. The Ponza grayling, Hipparchia sbordonii, is an endemic butterfly species that is currently found only in two tiny islands of the Pontine archipelago, off the coast of Italy, occupying an area smaller than 10 km2. It has been classified as Endangered (IUCN) because of the extremely limited area of occurrence, population fragmentation, and the recent demographic decline. Thanks to a combination of different assemblers of long and short genomic reads, bulk transcriptome RNAseq, and synteny analysis with phylogenetically close butterflies, we produced a highly contiguous, chromosome-scale annotated reference genome for the Ponza grayling, including 28 autosomes and the Z sexual chromosomes. The final assembly spanned 388.61 Gb with a contig N50 of 14.5 Mb and a BUSCO completeness score of 98.5%. Synteny analysis using four other butterfly species revealed high collinearity with Hipparchia semele and highlighted 10 intrachromosomal inversions longer than 10 kb, of which two appeared on the lineage leading to H. sbordonii. Our results show that a chromosome-scale reference genome is attainable also when chromatin conformation data may be impractical or present specific technical challenges. The high-quality genomic resource for H. sbordonii opens up new opportunities for the accurate assessment of genetic diversity and genetic load and for the investigations of the genomic novelties characterizing the evolutionary path of this endemic island species.
Collapse
Affiliation(s)
- Sebastiano Fava
- Department of Life and Environmental Sciences, Marche Polytechnic University, Ancona, Italy
| | - Marco Sollitto
- Department of Life Sciences, University of Trieste, Trieste, Italy
| | - Mbarsid Racaku
- Department of Life Sciences, University of Trieste, Trieste, Italy
| | | | - Andrea Benazzo
- Department of Life Sciences and Biotechnology, University of Ferrara, Ferrara, Italy
| | - Lorena Ancona
- Department of Life and Environmental Sciences, Marche Polytechnic University, Ancona, Italy
| | - Paolo Gratton
- Department of Biology, University of Rome “Tor Vergata”, Rome, Italy
| | - Fiorella Florian
- Department of Life Sciences, University of Trieste, Trieste, Italy
| | | | - Claudio Ciofi
- Department of Biology, University of Florence, Florence, Italy
| | | | - Marco Gerdol
- Department of Life Sciences, University of Trieste, Trieste, Italy
| | - Valerio Sbordoni
- Department of Biology, University of Rome “Tor Vergata”, Rome, Italy
| | - Giorgio Bertorelle
- Department of Life Sciences and Biotechnology, University of Ferrara, Ferrara, Italy
| | - Emiliano Trucchi
- Department of Life and Environmental Sciences, Marche Polytechnic University, Ancona, Italy
| |
Collapse
|
4
|
Bhattarai UR, Poulin R, Gemmell NJ, Dowle E. Genome assembly and annotation of the mermithid nematode Mermis nigrescens. G3 (BETHESDA, MD.) 2024; 14:jkae023. [PMID: 38301266 PMCID: PMC10989877 DOI: 10.1093/g3journal/jkae023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/29/2023] [Revised: 01/21/2024] [Accepted: 01/22/2024] [Indexed: 02/03/2024]
Abstract
Genetic studies of nematodes have been dominated by Caenorhabditis elegans as a model species. A lack of genomic resources has limited the expansion of genetic research to other groups of nematodes. Here, we report a draft genome assembly of a mermithid nematode, Mermis nigrescens. Mermithidae are insect parasitic nematodes with hosts including a wide range of terrestrial arthropods. We sequenced, assembled, and annotated the whole genome of M. nigrescens using nanopore long reads and 10X Chromium link reads. The assembly is 524 Mb in size consisting of 867 scaffolds. The N50 value is 2.42 Mb, and half of the assembly is in the 30 longest scaffolds. The assembly BUSCO score from the eukaryotic database (eukaryota_odb10) indicates that the genome is 86.7% complete and 5.1% partial. The genome has a high level of heterozygosity (6.6%) with a repeat content of 83.98%. mRNA-seq reads from different sized nematodes (≤2 cm, 3.5-7 cm, and >7 cm body length) representing different developmental stages were also generated and used for the genome annotation. Using ab initio and evidence-based gene model predictions, 12,313 protein-coding genes and 24,186 mRNAs were annotated. These genomic resources will help researchers investigate the various aspects of the biology and host-parasite interactions of mermithid nematodes.
Collapse
Affiliation(s)
- Upendra R Bhattarai
- Department of Anatomy, University of Otago, Dunedin 9016, New Zealand
- Department of Organismic & Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
| | - Robert Poulin
- Department of Zoology, University of Otago, Dunedin 9016, New Zealand
| | - Neil J Gemmell
- Department of Anatomy, University of Otago, Dunedin 9016, New Zealand
| | - Eddy Dowle
- Department of Anatomy, University of Otago, Dunedin 9016, New Zealand
| |
Collapse
|
5
|
Hoge C, de Manuel M, Mahgoub M, Okami N, Fuller Z, Banerjee S, Baker Z, McNulty M, Andolfatto P, Macfarlan TS, Schumer M, Tzika AC, Przeworski M. Patterns of recombination in snakes reveal a tug-of-war between PRDM9 and promoter-like features. Science 2024; 383:eadj7026. [PMID: 38386752 DOI: 10.1126/science.adj7026] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2023] [Accepted: 01/04/2024] [Indexed: 02/24/2024]
Abstract
In some mammals, notably humans, recombination occurs almost exclusively where the protein PRDM9 binds, whereas in vertebrates lacking an intact PRDM9, such as birds and canids, recombination rates are elevated near promoter-like features. To determine whether PRDM9 directs recombination in nonmammalian vertebrates, we focused on an exemplar species with a single, intact PRDM9 ortholog, the corn snake (Pantherophis guttatus). Analyzing historical recombination rates along the genome and crossovers in pedigrees, we found evidence that PRDM9 specifies the location of recombination events, but we also detected a separable effect of promoter-like features. These findings reveal that the uses of PRDM9 and promoter-like features need not be mutually exclusive and instead reflect a tug-of-war that is more even in some species than others.
Collapse
Affiliation(s)
- Carla Hoge
- Department of Biological Sciences, Columbia University, New York, NY, USA
| | - Marc de Manuel
- Department of Biological Sciences, Columbia University, New York, NY, USA
| | - Mohamed Mahgoub
- The Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, USA
| | - Naima Okami
- Department of Biological Sciences, Columbia University, New York, NY, USA
| | - Zachary Fuller
- Department of Biological Sciences, Columbia University, New York, NY, USA
| | - Shreya Banerjee
- Department of Biology, Stanford University, Stanford, CA, USA
| | - Zachary Baker
- Department of Systems Biology, Columbia University, New York, NY, USA
| | - Morgan McNulty
- Department of Biological Sciences, Columbia University, New York, NY, USA
| | - Peter Andolfatto
- Department of Biological Sciences, Columbia University, New York, NY, USA
| | - Todd S Macfarlan
- The Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, USA
| | - Molly Schumer
- Department of Biology, Stanford University, Stanford, CA, USA
- Howard Hughes Medical Institute, Stanford, CA, USA
| | - Athanasia C Tzika
- Laboratory of Artificial & Natural Evolution (LANE), Department of Genetics & Evolution, University of Geneva, Geneva, Switzerland
| | - Molly Przeworski
- Department of Biological Sciences, Columbia University, New York, NY, USA
- Department of Systems Biology, Columbia University, New York, NY, USA
| |
Collapse
|
6
|
Dimens PV, Jones KL, Margulies D, Scholey V, Cusatti S, McPeak B, Hildahl TE, Saillant EAE. Genomic resources for the Yellowfin tuna Thunnus albacares. Mol Biol Rep 2024; 51:232. [PMID: 38281308 DOI: 10.1007/s11033-023-09117-6] [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: 07/08/2023] [Accepted: 12/06/2023] [Indexed: 01/30/2024]
Abstract
BACKGROUND The Yellowfin tuna (Thunnus albacares) is a large tuna exploited by major fisheries in tropical and subtropical waters of all oceans except the Mediterranean Sea. Genomic studies of population structure, adaptive variation or of the genetic basis of phenotypic traits are needed to inform fisheries management but are currently limited by the lack of a reference genome for this species. Here we report a draft genome assembly and a linkage map for use in genomic studies of T. albacares. METHODS AND RESULTS Illumina and PacBio SMRT sequencing were used in combination to generate a hybrid assembly that comprises 743,073,847 base pairs contained in 2,661 scaffolds. The assembly has a N50 of 351,587 and complete and partial BUSCO scores of 86.47% and 3.63%, respectively. Double-digest restriction associated DNA (ddRAD) was used to genotype the 2 parents and 164 of their F1 offspring resulting from a controlled breeding cross, retaining 19,469 biallelic single nucleotide polymorphism (SNP) loci. The SNP loci were used to construct a linkage map that features 24 linkage groups that represent the 24 chromosomes of yellowfin tuna. The male and female maps span 1,243.8 cM and 1,222.9 cM, respectively. The map was used to anchor the assembly in 24 super-scaffolds that contain 79% of the yellowfin tuna genome. Gene prediction identified 46,992 putative genes 20,203 of which could be annotated via gene ontology. CONCLUSIONS The draft reference will be valuable to interpret studies of genome wide variation in T. albacares and other Scombroid species.
Collapse
Affiliation(s)
- Pavel V Dimens
- School of Ocean Science and Engineering, The University of Southern Mississippi, Ocean Springs, MS, 39564, USA
| | | | - Daniel Margulies
- Inter-American Tropical Tuna Commission, 8901 La Jolla Shores Drive, La Jolla, CA, 92037, USA
| | - Vernon Scholey
- Inter-American Tropical Tuna Commission, 8901 La Jolla Shores Drive, La Jolla, CA, 92037, USA
| | - Susana Cusatti
- Inter-American Tropical Tuna Commission, 8901 La Jolla Shores Drive, La Jolla, CA, 92037, USA
| | - Brooke McPeak
- School of Ocean Science and Engineering, The University of Southern Mississippi, Ocean Springs, MS, 39564, USA
| | - Tami E Hildahl
- School of Ocean Science and Engineering, The University of Southern Mississippi, Ocean Springs, MS, 39564, USA
| | - Eric A E Saillant
- School of Ocean Science and Engineering, The University of Southern Mississippi, Ocean Springs, MS, 39564, USA.
| |
Collapse
|
7
|
Cloutier S, Edwards T, Zheng C, Booker HM, Islam T, Nabetani K, Kutcher HR, Molina O, You FM. Fine-mapping of a major locus for Fusarium wilt resistance in flax (Linum usitatissimum L.). TAG. THEORETICAL AND APPLIED GENETICS. THEORETISCHE UND ANGEWANDTE GENETIK 2024; 137:27. [PMID: 38245903 PMCID: PMC10800302 DOI: 10.1007/s00122-023-04528-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/04/2023] [Accepted: 12/11/2023] [Indexed: 01/23/2024]
Abstract
KEY MESSAGE Fine-mapping of a locus on chromosome 1 of flax identified an S-lectin receptor-like kinase (SRLK) as the most likely candidate for a major Fusarium wilt resistance gene. Fusarium wilt, caused by the soil-borne fungal pathogen Fusarium oxysporum f. sp. lini, is a devastating disease in flax. Genetic resistance can counteract this disease and limit its spread. To map major genes for Fusarium wilt resistance, a recombinant inbred line population of more than 700 individuals derived from a cross between resistant cultivar 'Bison' and susceptible cultivar 'Novelty' was phenotyped in Fusarium wilt nurseries at two sites for two and three years, respectively. The population was genotyped with 4487 single nucleotide polymorphism (SNP) markers. Twenty-four QTLs were identified with IciMapping, 18 quantitative trait nucleotides with 3VmrMLM and 108 linkage disequilibrium blocks with RTM-GWAS. All models identified a major QTL on chromosome 1 that explained 20-48% of the genetic variance for Fusarium wilt resistance. The locus was estimated to span ~ 867 Kb but included a ~ 400 Kb unresolved region. Whole-genome sequencing of 'CDC Bethune', 'Bison' and 'Novelty' produced ~ 450 Kb continuous sequences of the locus. Annotation revealed 110 genes, of which six were considered candidate genes. Fine-mapping with 12 SNPs and 15 Kompetitive allele-specific PCR (KASP) markers narrowed down the interval to ~ 69 Kb, which comprised the candidate genes Lus10025882 and Lus10025891. The latter, a G-type S-lectin receptor-like kinase (SRLK) is the most likely resistance gene because it is the only polymorphic one. In addition, Fusarium wilt resistance genes previously isolated in tomato and Arabidopsis belonged to the SRLK class. The robust KASP markers can be used in marker-assisted breeding to select for this major Fusarium wilt resistance locus.
Collapse
Affiliation(s)
- S Cloutier
- Ottawa Research and Development Centre, Agriculture and Agri-Food Canada, 960 Carling Avenue, Ottawa, ON, K1A 0C6, Canada.
| | - T Edwards
- Ottawa Research and Development Centre, Agriculture and Agri-Food Canada, 960 Carling Avenue, Ottawa, ON, K1A 0C6, Canada
| | - C Zheng
- Ottawa Research and Development Centre, Agriculture and Agri-Food Canada, 960 Carling Avenue, Ottawa, ON, K1A 0C6, Canada
| | - H M Booker
- Crop Development Centre, University of Saskatchewan, 51 Campus Drive, Saskatoon, SK, S7N 5A8, Canada
- Department of Plant Agriculture, Ontario Agricultural College, University of Guelph, 50 Stone Road E, Guelph, ON, N1G 2W1, Canada
| | - T Islam
- Crop Development Centre, University of Saskatchewan, 51 Campus Drive, Saskatoon, SK, S7N 5A8, Canada
| | - K Nabetani
- Crop Development Centre, University of Saskatchewan, 51 Campus Drive, Saskatoon, SK, S7N 5A8, Canada
| | - H R Kutcher
- Crop Development Centre, University of Saskatchewan, 51 Campus Drive, Saskatoon, SK, S7N 5A8, Canada
| | - O Molina
- Morden Research and Development Centre, Agriculture and Agri-Food Canada, 101 Route 100, Morden, MB, R6M 1Y5, Canada
| | - F M You
- Ottawa Research and Development Centre, Agriculture and Agri-Food Canada, 960 Carling Avenue, Ottawa, ON, K1A 0C6, Canada.
| |
Collapse
|
8
|
Xiong W, van Workum DJM, Berke L, Bakker LV, Schijlen E, Becker FFM, van de Geest H, Peters S, Michelmore R, van Treuren R, Jeuken M, Smit S, Schranz ME. Genome assembly and analysis of Lactuca virosa: implications for lettuce breeding. G3 (BETHESDA, MD.) 2023; 13:jkad204. [PMID: 37740775 PMCID: PMC10627274 DOI: 10.1093/g3journal/jkad204] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Revised: 07/16/2023] [Accepted: 07/19/2023] [Indexed: 09/25/2023]
Abstract
Lettuce (Lactuca sativa L.) is a leafy vegetable crop with ongoing breeding efforts related to quality, resilience, and innovative production systems. To breed resilient and resistant lettuce in the future, valuable genetic variation found in close relatives could be further exploited. Lactuca virosa (2x = 2n = 18), a wild relative assigned to the tertiary lettuce gene pool, has a much larger genome (3.7 Gbp) than Lactuca sativa (2.5 Gbp). It has been used in interspecific crosses and is a donor to modern crisphead lettuce cultivars. Here, we present a de novo reference assembly of L. virosa with high continuity and complete gene space. This assembly facilitated comparisons to the genome of L. sativa and to that of the wild species L. saligna, a representative of the secondary lettuce gene pool. To assess the diversity in gene content, we classified the genes of the 3 Lactuca species as core, accessory, and unique. In addition, we identified 3 interspecific chromosomal inversions compared to L. sativa, which each may cause recombination suppression and thus hamper future introgression breeding. Using 3-way comparisons in both reference-based and reference-free manners, we show that the proliferation of long-terminal repeat elements has driven the genome expansion of L. virosa. Further, we performed a genome-wide comparison of immune genes, nucleotide-binding leucine-rich repeat, and receptor-like kinases among Lactuca spp. and indicated the evolutionary patterns and mechanisms behind their expansions. These genome analyses greatly facilitate the understanding of genetic variation in L. virosa, which is beneficial for the breeding of improved lettuce varieties.
Collapse
Affiliation(s)
- Wei Xiong
- Biosystematics Group, Wageningen University & Research, P.O. Box 16, Wageningen, 6700 AA, The Netherlands
| | - Dirk-Jan M van Workum
- Bioinformatics Group, Wageningen University & Research, P.O. Box 633, Wageningen, 6700 AP, The Netherlands
| | - Lidija Berke
- Biosystematics Group, Wageningen University & Research, P.O. Box 16, Wageningen, 6700 AA, The Netherlands
| | - Linda V Bakker
- Bioscience, Wageningen University & Research, P.O. Box 16, Wageningen, 6700 AA, The Netherlands
| | - Elio Schijlen
- Bioscience, Wageningen University & Research, P.O. Box 16, Wageningen, 6700 AA, The Netherlands
| | - Frank F M Becker
- Biosystematics Group, Wageningen University & Research, P.O. Box 16, Wageningen, 6700 AA, The Netherlands
- Laboratory of Genetics, Wageningen University & Research, P.O. Box 16, Wageningen, 6700 AA, The Netherlands
| | - Henri van de Geest
- Bioscience, Wageningen University & Research, P.O. Box 16, Wageningen, 6700 AA, The Netherlands
| | - Sander Peters
- Bioscience, Wageningen University & Research, P.O. Box 16, Wageningen, 6700 AA, The Netherlands
| | - Richard Michelmore
- The Genome Center, Genome & Biomedical Sciences Facility, University of California, Davis, 451 East Health Sciences Drive, Davis, CA 95616-8816, USA
| | - Rob van Treuren
- Centre for Genetic Resources, the Netherlands (CGN), Wageningen University & Research, P.O. Box 16, Wageningen, 6700 AA, The Netherlands
| | - Marieke Jeuken
- Plant Breeding Group, Wageningen University & Research, P.O. Box 386, Wageningen, 6700 AJ, The Netherlands
| | - Sandra Smit
- Bioinformatics Group, Wageningen University & Research, P.O. Box 633, Wageningen, 6700 AP, The Netherlands
| | - M Eric Schranz
- Biosystematics Group, Wageningen University & Research, P.O. Box 16, Wageningen, 6700 AA, The Netherlands
| |
Collapse
|
9
|
Zheng JW, Lu Y, Yang YF, Huang D, Li DW, Wang X, Gao Y, Yang WD, Guan Y, Li HY. Systematic dissection of genomic features determining the vast diversity of conotoxins. BMC Genomics 2023; 24:598. [PMID: 37814244 PMCID: PMC10561478 DOI: 10.1186/s12864-023-09689-4] [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: 02/13/2023] [Accepted: 09/20/2023] [Indexed: 10/11/2023] Open
Abstract
BACKGROUND Conus, a highly diverse species of venomous predators, has attracted significant attention in neuroscience and new drug development due to their rich collection of neuroactive peptides called conotoxins. Recent advancements in transcriptome, proteome, and genome analyses have facilitated the identification of conotoxins within Conus' venom glands, providing insights into the genetic features and evolutionary patterns of conotoxin genes. However, the underlying mechanism behind the extraordinary hypervariability of conotoxins remains largely unknown. RESULTS We analyzed the transcriptomes of 34 Conus species, examining various tissues such as the venom duct, venom bulb, and salivary gland, leading to the identification of conotoxin genes. Genetic variation analysis revealed that a subset of these genes (15.78% of the total) in Conus species underwent positive selection (Ka/Ks > 1, p < 0.01). Additionally, we reassembled and annotated the genome of C. betulinus, uncovering 221 conotoxin-encoding genes. These genes primarily consisted of three exons, with a significant portion showing high transcriptional activity in the venom ducts. Importantly, the flanking regions and adjacent introns of conotoxin genes exhibited a higher prevalence of transposon elements, suggesting their potential contribution to the extensive variability observed in conotoxins. Furthermore, we detected genome duplication in C. betulinus, which likely contributed to the expansion of conotoxin gene numbers. Interestingly, our study also provided evidence of introgression among Conus species, indicating that interspecies hybridization may have played a role in shaping the evolution of diverse conotoxin genes. CONCLUSIONS This study highlights the impact of adaptive evolution and introgressive hybridization on the genetic diversity of conotoxin genes and the evolution of Conus. We also propose a hypothesis suggesting that transposable elements might significantly contribute to the remarkable diversity observed in conotoxins. These findings not only enhance our understanding of peptide genetic diversity but also present a novel approach for peptide bioengineering.
Collapse
Affiliation(s)
- Jian-Wei Zheng
- Key Laboratory of Aquatic Eutrophication and Control of Harmful Algal Blooms of Guangdong Higher Education Institute, College of Life Science and Technology, Jinan University, Guangzhou, 510632, China
- College of Food Science and Engineering, Foshan University of Science and Technology, Foshan, 528231, China
| | - Yang Lu
- Key Laboratory of Aquatic Eutrophication and Control of Harmful Algal Blooms of Guangdong Higher Education Institute, College of Life Science and Technology, Jinan University, Guangzhou, 510632, China
| | - Yu-Feng Yang
- Key Laboratory of Aquatic Eutrophication and Control of Harmful Algal Blooms of Guangdong Higher Education Institute, College of Life Science and Technology, Jinan University, Guangzhou, 510632, China
| | - Dan Huang
- Key Laboratory of Aquatic Eutrophication and Control of Harmful Algal Blooms of Guangdong Higher Education Institute, College of Life Science and Technology, Jinan University, Guangzhou, 510632, China
| | - Da-Wei Li
- Key Laboratory of Aquatic Eutrophication and Control of Harmful Algal Blooms of Guangdong Higher Education Institute, College of Life Science and Technology, Jinan University, Guangzhou, 510632, China
| | - Xiang Wang
- Key Laboratory of Aquatic Eutrophication and Control of Harmful Algal Blooms of Guangdong Higher Education Institute, College of Life Science and Technology, Jinan University, Guangzhou, 510632, China
| | - Yang Gao
- Gulou Hospital, Nanjing University, Nanjing, China
| | - Wei-Dong Yang
- Key Laboratory of Aquatic Eutrophication and Control of Harmful Algal Blooms of Guangdong Higher Education Institute, College of Life Science and Technology, Jinan University, Guangzhou, 510632, China
| | - Yuanfang Guan
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, USA
| | - Hong-Ye Li
- Key Laboratory of Aquatic Eutrophication and Control of Harmful Algal Blooms of Guangdong Higher Education Institute, College of Life Science and Technology, Jinan University, Guangzhou, 510632, China.
| |
Collapse
|
10
|
Cicconardi F, Milanetti E, Pinheiro de Castro EC, Mazo-Vargas A, Van Belleghem SM, Ruggieri AA, Rastas P, Hanly J, Evans E, Jiggins CD, Owen McMillan W, Papa R, Di Marino D, Martin A, Montgomery SH. Evolutionary dynamics of genome size and content during the adaptive radiation of Heliconiini butterflies. Nat Commun 2023; 14:5620. [PMID: 37699868 PMCID: PMC10497600 DOI: 10.1038/s41467-023-41412-5] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2022] [Accepted: 08/30/2023] [Indexed: 09/14/2023] Open
Abstract
Heliconius butterflies, a speciose genus of Müllerian mimics, represent a classic example of an adaptive radiation that includes a range of derived dietary, life history, physiological and neural traits. However, key lineages within the genus, and across the broader Heliconiini tribe, lack genomic resources, limiting our understanding of how adaptive and neutral processes shaped genome evolution during their radiation. Here, we generate highly contiguous genome assemblies for nine Heliconiini, 29 additional reference-assembled genomes, and improve 10 existing assemblies. Altogether, we provide a dataset of annotated genomes for a total of 63 species, including 58 species within the Heliconiini tribe. We use this extensive dataset to generate a robust and dated heliconiine phylogeny, describe major patterns of introgression, explore the evolution of genome architecture, and the genomic basis of key innovations in this enigmatic group, including an assessment of the evolution of putative regulatory regions at the Heliconius stem. Our work illustrates how the increased resolution provided by such dense genomic sampling improves our power to generate and test gene-phenotype hypotheses, and precisely characterize how genomes evolve.
Collapse
Affiliation(s)
- Francesco Cicconardi
- School of Biological Sciences, Bristol University, Bristol, United Kingdom.
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom.
| | - Edoardo Milanetti
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, 00185, Rome, Italy
- Center for Life Nano- & Neuro-Science, Italian Institute of Technology, Viale Regina Elena 291, 00161, Rome, Italy
| | | | - Anyi Mazo-Vargas
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, 14853, USA
| | - Steven M Van Belleghem
- Department of Biology, University of Puerto Rico, Rio Piedras, PR, Puerto Rico
- Ecology, Evolution and Conservation Biology, Biology Department, KU Leuven, Leuven, Belgium
| | | | - Pasi Rastas
- Institute of Biotechnology, University of Helsinki, Helsinki, Finland
| | - Joseph Hanly
- Department of Biological Sciences, The George Washington University, Washington DC, WA, 20052, USA
- Smithsonian Tropical Research Institute, Panama City, Panama
| | - Elizabeth Evans
- Department of Biology, University of Puerto Rico, Rio Piedras, PR, Puerto Rico
| | - Chris D Jiggins
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - W Owen McMillan
- Smithsonian Tropical Research Institute, Panama City, Panama
| | - Riccardo Papa
- Department of Biology, University of Puerto Rico, Rio Piedras, PR, Puerto Rico
- Molecular Sciences and Research Center, University of Puerto Rico, San Juan, PR, Puerto Rico
- Comprehensive Cancer Center, University of Puerto Rico, San Juan, PR, Puerto Rico
| | - Daniele Di Marino
- Department of Life and Environmental Sciences, New York-Marche Structural Biology Center (NY-MaSBiC), Polytechnic University of Marche, Via Brecce Bianche, 60131, Ancona, Italy
- Neuronal Death and Neuroprotection Unit, Department of Neuroscience, Mario Negri Institute for Pharmacological Research-IRCCS, Via Mario Negri 2, 20156, Milano, Italy
- National Biodiversity Future Center (NBFC), Palermo, Italy
| | - Arnaud Martin
- Department of Biological Sciences, The George Washington University, Washington DC, WA, 20052, USA
| | - Stephen H Montgomery
- School of Biological Sciences, Bristol University, Bristol, United Kingdom.
- Smithsonian Tropical Research Institute, Panama City, Panama.
| |
Collapse
|
11
|
Schmeing S, Robinson MD. Gapless provides combined scaffolding, gap filling, and assembly correction with long reads. Life Sci Alliance 2023; 6:e202201471. [PMID: 37142439 PMCID: PMC10166144 DOI: 10.26508/lsa.202201471] [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: 04/03/2022] [Revised: 04/17/2023] [Accepted: 04/18/2023] [Indexed: 05/06/2023] Open
Abstract
Continuity, correctness, and completeness of genome assemblies are important for many biological projects. Long reads represent a major driver towards delivering high-quality genomes, but not everybody can achieve the necessary coverage for good long read-only assemblies. Therefore, improving existing assemblies with low-coverage long reads is a promising alternative. The improvements include correction, scaffolding, and gap filling. However, most tools perform only one of these tasks and the useful information of reads that supported the scaffolding is lost when running separate programs successively. Therefore, we propose a new tool for combined execution of all three tasks using PacBio or Oxford Nanopore reads. gapless is available at: https://github.com/schmeing/gapless.
Collapse
Affiliation(s)
- Stephan Schmeing
- Department of Molecular Life Sciences, University of Zurich, Zurich, Switzerland
- SIB Swiss Institute of Bioinformatics, University of Zurich, Zurich, Switzerland
| | - Mark D Robinson
- Department of Molecular Life Sciences, University of Zurich, Zurich, Switzerland
- SIB Swiss Institute of Bioinformatics, University of Zurich, Zurich, Switzerland
| |
Collapse
|
12
|
Wong J, Coombe L, Nikolić V, Zhang E, Nip KM, Sidhu P, Warren RL, Birol I. Linear time complexity de novo long read genome assembly with GoldRush. Nat Commun 2023; 14:2906. [PMID: 37217507 DOI: 10.1038/s41467-023-38716-x] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2022] [Accepted: 05/11/2023] [Indexed: 05/24/2023] Open
Abstract
Current state-of-the-art de novo long read genome assemblers follow the Overlap-Layout-Consensus paradigm. While read-to-read overlap - its most costly step - was improved in modern long read genome assemblers, these tools still often require excessive RAM when assembling a typical human dataset. Our work departs from this paradigm, foregoing all-vs-all sequence alignments in favor of a dynamic data structure implemented in GoldRush, a de novo long read genome assembly algorithm with linear time complexity. We tested GoldRush on Oxford Nanopore Technologies long sequencing read datasets with different base error profiles sourced from three human cell lines, rice, and tomato. Here, we show that GoldRush achieves assembly scaffold NGA50 lengths of 18.3-22.2, 0.3 and 2.6 Mbp, for the genomes of human, rice, and tomato, respectively, and assembles each genome within a day, using at most 54.5 GB of random-access memory, demonstrating the scalability of our genome assembly paradigm and its implementation.
Collapse
Affiliation(s)
- Johnathan Wong
- Canada's Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC, V5Z 4S6, Canada.
| | - Lauren Coombe
- Canada's Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC, V5Z 4S6, Canada
| | - Vladimir Nikolić
- Canada's Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC, V5Z 4S6, Canada
| | - Emily Zhang
- Canada's Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC, V5Z 4S6, Canada
| | - Ka Ming Nip
- Canada's Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC, V5Z 4S6, Canada
| | - Puneet Sidhu
- Canada's Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC, V5Z 4S6, Canada
| | - René L Warren
- Canada's Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC, V5Z 4S6, Canada
| | - Inanç Birol
- Canada's Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, BC, V5Z 4S6, Canada.
| |
Collapse
|
13
|
Luo J, Guan T, Chen G, Yu Z, Zhai H, Yan C, Luo H. SLHSD: hybrid scaffolding method based on short and long reads. Brief Bioinform 2023; 24:7152317. [PMID: 37141142 DOI: 10.1093/bib/bbad169] [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: 07/03/2022] [Revised: 01/08/2023] [Accepted: 04/12/2023] [Indexed: 05/05/2023] Open
Abstract
In genome assembly, scaffolding can obtain more complete and continuous scaffolds. Current scaffolding methods usually adopt one type of read to construct a scaffold graph and then orient and order contigs. However, scaffolding with the strengths of two or more types of reads seems to be a better solution to some tricky problems. Combining the advantages of different types of data is significant for scaffolding. Here, a hybrid scaffolding method (SLHSD) is present that simultaneously leverages the precision of short reads and the length advantage of long reads. Building an optimal scaffold graph is an important foundation for getting scaffolds. SLHSD uses a new algorithm that combines long and short read alignment information to determine whether to add an edge and how to calculate the edge weight in a scaffold graph. In addition, SLHSD develops a strategy to ensure that edges with high confidence can be added to the graph with priority. Then, a linear programming model is used to detect and remove remaining false edges in the graph. We compared SLHSD with other scaffolding methods on five datasets. Experimental results show that SLHSD outperforms other methods. The open-source code of SLHSD is available at https://github.com/luojunwei/SLHSD.
Collapse
Affiliation(s)
- Junwei Luo
- School of Software, Henan Polytechnic University, Jiaozuo 454003, China
| | - Ting Guan
- School of Software, Henan Polytechnic University, Jiaozuo 454003, China
| | - Guolin Chen
- School of Software, Henan Polytechnic University, Jiaozuo 454003, China
| | - Zhonghua Yu
- School of Software, Henan Polytechnic University, Jiaozuo 454003, China
| | - Haixia Zhai
- School of Software, Henan Polytechnic University, Jiaozuo 454003, China
| | - Chaokun Yan
- School of Computer and Information Engineering, Henan University, Kaifeng 475001, China
| | - Huimin Luo
- School of Computer and Information Engineering, Henan University, Kaifeng 475001, China
| |
Collapse
|
14
|
Coombe L, Warren RL, Wong J, Nikolic V, Birol I. ntLink: A Toolkit for De Novo Genome Assembly Scaffolding and Mapping Using Long Reads. Curr Protoc 2023; 3:e733. [PMID: 37039735 PMCID: PMC10091225 DOI: 10.1002/cpz1.733] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
Abstract
With the increasing affordability and accessibility of genome sequencing data, de novo genome assembly is an important first step to a wide variety of downstream studies and analyses. Therefore, bioinformatics tools that enable the generation of high-quality genome assemblies in a computationally efficient manner are essential. Recent developments in long-read sequencing technologies have greatly benefited genome assembly work, including scaffolding, by providing long-range evidence that can aid in resolving the challenging repetitive regions of complex genomes. ntLink is a flexible and resource-efficient genome scaffolding tool that utilizes long-read sequencing data to improve upon draft genome assemblies built from any sequencing technologies, including the same long reads. Instead of using read alignments to identify candidate joins, ntLink utilizes minimizer-based mappings to infer how input sequences should be ordered and oriented into scaffolds. Recent improvements to ntLink have added important features such as overlap detection, gap-filling, and in-code scaffolding iterations. Here, we present three basic protocols demonstrating how to use each of these new features to yield highly contiguous genome assemblies, while still maintaining ntLink's proven computational efficiency. Further, as we illustrate in the alternate protocols, the lightweight minimizer-based mappings that enable ntLink scaffolding can also be utilized for other downstream applications, such as misassembly detection. With its modularity and multiple modes of execution, ntLink has broad benefit to the genomics community, from genome scaffolding and beyond. ntLink is an open-source project and is freely available from https://github.com/bcgsc/ntLink. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: ntLink scaffolding using overlap detection Basic Protocol 2: ntLink scaffolding with gap-filling Basic Protocol 3: Running in-code iterations of ntLink scaffolding Alternate Protocol 1: Generating long-read to contig mappings with ntLink Alternate Protocol 2: Using ntLink mappings for genome assembly correction with Tigmint-long Support Protocol: Installing ntLink.
Collapse
Affiliation(s)
- Lauren Coombe
- Canada’s Michael Smith Genome Sciences Centre at BC Cancer, 570 W 7 Ave, Vancouver, BC V5Z 4S6, 604-707-5900
| | - René L. Warren
- Canada’s Michael Smith Genome Sciences Centre at BC Cancer, 570 W 7 Ave, Vancouver, BC V5Z 4S6, 604-707-5900
| | - Johnathan Wong
- Canada’s Michael Smith Genome Sciences Centre at BC Cancer, 570 W 7 Ave, Vancouver, BC V5Z 4S6, 604-707-5900
| | - Vladimir Nikolic
- Canada’s Michael Smith Genome Sciences Centre at BC Cancer, 570 W 7 Ave, Vancouver, BC V5Z 4S6, 604-707-5900
| | - Inanc Birol
- Canada’s Michael Smith Genome Sciences Centre at BC Cancer, 570 W 7 Ave, Vancouver, BC V5Z 4S6, 604-707-5900
| |
Collapse
|
15
|
Laso-Pérez R, Wu F, Crémière A, Speth DR, Magyar JS, Zhao K, Krupovic M, Orphan VJ. Evolutionary diversification of methanotrophic ANME-1 archaea and their expansive virome. Nat Microbiol 2023; 8:231-245. [PMID: 36658397 PMCID: PMC9894754 DOI: 10.1038/s41564-022-01297-4] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2022] [Accepted: 11/29/2022] [Indexed: 01/21/2023]
Abstract
'Candidatus Methanophagales' (ANME-1) is an order-level clade of archaea responsible for anaerobic methane oxidation in deep-sea sediments. The diversity, ecology and evolution of ANME-1 remain poorly understood. In this study, we use metagenomics on deep-sea hydrothermal samples to expand ANME-1 diversity and uncover the effect of virus-host dynamics. Phylogenetic analyses reveal a deep-branching, thermophilic family, 'Candidatus Methanospirareceae', closely related to short-chain alkane oxidizers. Global phylogeny and near-complete genomes show that hydrogen metabolism within ANME-1 is an ancient trait that was vertically inherited but differentially lost during lineage diversification. Metagenomics also uncovered 16 undescribed virus families so far exclusively targeting ANME-1 archaea, showing unique structural and replicative signatures. The expansive ANME-1 virome contains a metabolic gene repertoire that can influence host ecology and evolution through virus-mediated gene displacement. Our results suggest an evolutionary continuum between anaerobic methane and short-chain alkane oxidizers and underscore the effects of viruses on the dynamics and evolution of methane-driven ecosystems.
Collapse
Affiliation(s)
- Rafael Laso-Pérez
- MARUM, Center for Marine Environmental Science, and Department of Geosciences, University of Bremen, Bremen, Germany.
- Systems Biology Department, Centro Nacional de Biotecnología (CNB-CSIC), Madrid, Spain.
| | - Fabai Wu
- ZJU-Hangzhou Global Scientific and Technological Innovation Center, Hangzhou, China.
- Ocean College, Zhejiang University, Zhoushan, China.
- Donghai Laboratory, Zhoushan, China.
- Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA, USA.
| | - Antoine Crémière
- Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA, USA
| | - Daan R Speth
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, USA
- Max-Planck Institute for Marine Microbiology, Bremen, Germany
| | - John S Magyar
- Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA, USA
| | - Kehan Zhao
- ZJU-Hangzhou Global Scientific and Technological Innovation Center, Hangzhou, China
| | - Mart Krupovic
- Institut Pasteur, Université Paris Cité, CNRS UMR6047, Archaeal Virology Unit, Paris, France.
| | - Victoria J Orphan
- Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA, USA.
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, USA.
| |
Collapse
|
16
|
Merges D, Dal Grande F, Valim H, Singh G, Schmitt I. Gene abundance linked to climate zone: Parallel evolution of gene content along elevation gradients in lichenized fungi. Front Microbiol 2023; 14:1097787. [PMID: 37032854 PMCID: PMC10073550 DOI: 10.3389/fmicb.2023.1097787] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2022] [Accepted: 02/23/2023] [Indexed: 04/11/2023] Open
Abstract
Introduction Intraspecific genomic variability affects a species' adaptive potential toward climatic conditions. Variation in gene content across populations and environments may point at genomic adaptations to specific environments. The lichen symbiosis, a stable association of fungal and photobiont partners, offers an excellent system to study environmentally driven gene content variation. Many of these species have remarkable environmental tolerances, and often form populations across different climate zones. Here, we combine comparative and population genomics to assess the presence and absence of genes in high and low elevation genomes of two lichenized fungi of the genus Umbilicaria. Methods The two species have non-overlapping ranges, but occupy similar climatic niches in North America (U. phaea) and Europe (U. pustulata): high elevation populations are located in the cold temperate zone and low elevation populations in the Mediterranean zone. We assessed gene content variation along replicated elevation gradients in each of the two species, based on a total of 2050 individuals across 26 populations. Specifically, we assessed shared orthologs across species within the same climate zone, and tracked, which genes increase or decrease in abundance within populations along elevation. Results In total, we found 16 orthogroups with shared orthologous genes in genomes at low elevation and 13 at high elevation. Coverage analysis revealed one ortholog that is exclusive to genomes at low elevation. Conserved domain search revealed domains common to the protein kinase superfamily. We traced the discovered ortholog in populations along five replicated elevation gradients on both continents and found that the number of this protein kinase gene linearly declined in abundance with increasing elevation, and was absent in the highest populations. Discussion We consider the parallel loss of an ortholog in two species and in two geographic settings a rare find, and a step forward in understanding the genomic underpinnings of climatic tolerances in lichenized fungi. In addition, the tracking of gene content variation provides a widely applicable framework for retrieving biogeographical determinants of gene presence/absence patterns. Our work provides insights into gene content variation of lichenized fungi in relation to climatic gradients, suggesting a new research direction with implications for understanding evolutionary trajectories of complex symbioses in relation to climatic change.
Collapse
Affiliation(s)
- Dominik Merges
- Senckenberg Biodiversity and Climate Research Centre, Frankfurt am Main, Germany
- LOEWE Centre for Translational Biodiversity Genomics (LOEWE-TBG), Frankfurt am Main, Germany
- Department of Forest Mycology and Plant Pathology, Swedish University of Agricultural Sciences, Uppsala, Sweden
- *Correspondence: Dominik Merges,
| | - Francesco Dal Grande
- Senckenberg Biodiversity and Climate Research Centre, Frankfurt am Main, Germany
- LOEWE Centre for Translational Biodiversity Genomics (LOEWE-TBG), Frankfurt am Main, Germany
- Department of Biology, University of Padova, Padua, Italy
- National Biodiversity Future Center (NBFC), Palermo, Italy
| | - Henrique Valim
- Senckenberg Biodiversity and Climate Research Centre, Frankfurt am Main, Germany
- LOEWE Centre for Translational Biodiversity Genomics (LOEWE-TBG), Frankfurt am Main, Germany
| | - Garima Singh
- Senckenberg Biodiversity and Climate Research Centre, Frankfurt am Main, Germany
- LOEWE Centre for Translational Biodiversity Genomics (LOEWE-TBG), Frankfurt am Main, Germany
- Department of Biology, University of Padova, Padua, Italy
| | - Imke Schmitt
- Senckenberg Biodiversity and Climate Research Centre, Frankfurt am Main, Germany
- LOEWE Centre for Translational Biodiversity Genomics (LOEWE-TBG), Frankfurt am Main, Germany
- Goethe University Frankfurt, Institute of Ecology, Evolution and Diversity, Frankfurt am Main, Germany
| |
Collapse
|
17
|
Barros CP, Derks MFL, Mohr J, Wood BJ, Crooijmans RPMA, Megens HJ, Bink MCAM, Groenen MAM. A new haplotype-resolved turkey genome to enable turkey genetics and genomics research. Gigascience 2022; 12:giad051. [PMID: 37489751 PMCID: PMC10360393 DOI: 10.1093/gigascience/giad051] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2022] [Revised: 12/12/2022] [Accepted: 06/27/2023] [Indexed: 07/26/2023] Open
Abstract
BACKGROUND The domesticated turkey (Meleagris gallopavo) is a species of significant agricultural importance and is the second largest contributor, behind broiler chickens, to world poultry meat production. The previous genome is of draft quality and partly based on the chicken (Gallus gallus) genome. A high-quality reference genome of M. gallopavo is essential for turkey genomics and genetics research and the breeding industry. RESULTS By adopting the trio-binning approach, we were able to assemble a high-quality chromosome-level F1 assembly and 2 parental haplotype assemblies, leveraging long-read technologies and genome-wide chromatin interaction data (Hi-C). From a total of 40 chromosomes (2n = 80), we captured 35 chromosomes in a single scaffold, showing much improved genome completeness and continuity compared to the old assembly build. The 3 assemblies are of higher quality than the previous draft quality assembly and comparable to the chicken assemblies (GRCg7) shown by the largest contig N50 (26.6 Mb) and comparable BUSCO gene set completeness scores (96-97%). Comparative analyses confirm a previously identified large inversion of around 19 Mbp on the Z chromosome not found in other Galliformes. Structural variation between the parent haplotypes was identified, which poses potential new target genes for breeding. CONCLUSIONS We contribute a new high-quality turkey genome at the chromosome level, benefiting turkey genetics and other avian genomics research as well as the turkey breeding industry.
Collapse
Affiliation(s)
- Carolina P Barros
- Wageningen University and Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| | - Martijn F L Derks
- Wageningen University and Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| | - Jeff Mohr
- Hybrid Turkeys, 650 Riverbend Drive Suite C, Kitchener, ON N2K 3S2, Canada
| | - Benjamin J Wood
- Hybrid Turkeys, 650 Riverbend Drive Suite C, Kitchener, ON N2K 3S2, Canada
- School of Veterinary Science, University of Queensland, Gatton, QLD 4343, Australia
| | | | - Hendrik-Jan Megens
- Wageningen University and Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| | - Marco C A M Bink
- Hendrix Genetics Research, Technology & Services, Boxmeer, AC 5830, The Netherlands
| | - Martien A M Groenen
- Wageningen University and Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| |
Collapse
|
18
|
Bentz EJ, Ophir AG. Chromosome-scale genome assembly of the African giant pouched rat (Cricetomys ansorgei) and evolutionary analysis reveals evidence of olfactory specialization. Genomics 2022; 114:110521. [PMID: 36351561 DOI: 10.1016/j.ygeno.2022.110521] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2022] [Revised: 10/28/2022] [Accepted: 11/04/2022] [Indexed: 11/07/2022]
Abstract
The Southern giant pouched rat, Cricetomys ansorgei, is a large rodent best known for its ability to detect landmines using its impressive sense of smell. Their powerful chemosensory abilities enable subtle discrimination of chemical social signals, and female pouched rats demonstrate a unique reproductive physiology hypothesized to be mediated by pheromonal mechanisms. Thus, C. ansorgei represents a novel mammalian model for chemosensory physiology, social behavior, and pheromonal control of reproductive physiology. We present the first chromosome-scale genomic sequence of the pouched rat encoding 22,671 protein coding genes, including 1571 olfactory receptors, and provide a glance into the evolutionary history of this species. Functional enrichment analysis reveals genetic expansions specific to the pouched rat are enriched for functions related to olfactory specialization. Overall, this assembly is of reference-quality, and will serve as a useful and informative genomic sequence on which we can confidently base future molecular research involving the pouched rat.
Collapse
Affiliation(s)
- Ehren J Bentz
- Department of Psychology, Cornell University, Ithaca, NY, USA.
| | | |
Collapse
|
19
|
Singh G, Dal Grande F, Schmitt I. Genome mining as a biotechnological tool for the discovery of novel biosynthetic genes in lichens. FRONTIERS IN FUNGAL BIOLOGY 2022; 3:993171. [PMID: 37746187 PMCID: PMC10512267 DOI: 10.3389/ffunb.2022.993171] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/13/2022] [Accepted: 08/30/2022] [Indexed: 09/26/2023]
Abstract
Natural products (NPs) and their derivatives are a major contributor to modern medicine. Historically, microorganisms such as bacteria and fungi have been instrumental in generating drugs and lead compounds because of the ease of culturing and genetically manipulating them. However, the ever-increasing demand for novel drugs highlights the need to bioprospect previously unexplored taxa for their biosynthetic potential. Next-generation sequencing technologies have expanded the range of organisms that can be explored for their biosynthetic content, as these technologies can provide a glimpse of an organism's entire biosynthetic landscape, without the need for cultivation. The entirety of biosynthetic genes can be compared to the genes of known function to identify the gene clusters potentially coding for novel products. In this study, we mine the genomes of nine lichen-forming fungal species of the genus Umbilicaria for biosynthetic genes, and categorize the biosynthetic gene clusters (BGCs) as "associated product structurally known" or "associated product putatively novel". Although lichen-forming fungi have been suggested to be a rich source of NPs, it is not known how their biosynthetic diversity compares to that of bacteria and non-lichenized fungi. We found that 25%-30% of biosynthetic genes are divergent as compared to the global database of BGCs, which comprises 1,200,000 characterized biosynthetic genes from plants, bacteria, and fungi. Out of 217 BGCs, 43 were highly divergant suggesting that they potentially encode structurally and functionally novel NPs. Clusters encoding the putatively novel metabolic diversity comprise polyketide synthases (30), non-ribosomal peptide synthetases (12), and terpenes (1). Our study emphasizes the utility of genomic data in bioprospecting microorganisms for their biosynthetic potential and in advancing the industrial application of unexplored taxa. We highlight the untapped structural metabolic diversity encoded in the lichenized fungal genomes. To the best of our knowledge, this is the first investigation identifying genes coding for NPs with potentially novel properties in lichenized fungi.
Collapse
Affiliation(s)
- Garima Singh
- Senckenberg Biodiversity and Climate Research Centre (SBiK-F), Frankfurt am Main, Germany
- LOEWE Center for Translational Biodiversity Genomics (TBG), Frankfurt am Main, Germany
- Department of Biology, University of Padova, Padova, Italy
| | - Francesco Dal Grande
- Senckenberg Biodiversity and Climate Research Centre (SBiK-F), Frankfurt am Main, Germany
- LOEWE Center for Translational Biodiversity Genomics (TBG), Frankfurt am Main, Germany
- Department of Biology, University of Padova, Padova, Italy
| | - Imke Schmitt
- Senckenberg Biodiversity and Climate Research Centre (SBiK-F), Frankfurt am Main, Germany
- LOEWE Center for Translational Biodiversity Genomics (TBG), Frankfurt am Main, Germany
- Institute of Ecology, Diversity and Evolution, Goethe University, Frankfurt am Main, Germany
| |
Collapse
|
20
|
Lan T, Li H, Yang S, Shi M, Han L, Sahu SK, Lu Y, Wang J, Zhou M, Liu H, Huang J, Wang Q, Zhu Y, Wang L, Xu Y, Lin C, Liu H, Hou Z. The chromosome-scale genome of the raccoon dog: Insights into its evolutionary characteristics. iScience 2022; 25:105117. [PMID: 36185367 PMCID: PMC9523411 DOI: 10.1016/j.isci.2022.105117] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2022] [Revised: 07/07/2022] [Accepted: 09/08/2022] [Indexed: 11/28/2022] Open
Affiliation(s)
- Tianming Lan
- BGI Life Science Joint Research Center, Northeast Forestry University, Harbin 150040, China
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen 518083, China
| | - Haimeng Li
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen 518083, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Shangchen Yang
- College of Life Sciences, Zhejiang University, Hangzhou 310058, China
| | - Minhui Shi
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen 518083, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Lei Han
- College of Wildlife and Protected Area, Northeast Forestry University, Harbin 150040, China
| | - Sunil Kumar Sahu
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen 518083, China
| | - Yaxian Lu
- College of Wildlife and Protected Area, Northeast Forestry University, Harbin 150040, China
| | - Jiangang Wang
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen 518083, China
| | - Mengchao Zhou
- College of Wildlife and Protected Area, Northeast Forestry University, Harbin 150040, China
| | - Hui Liu
- Key Laboratory of Genetics and Germplasm Innovation of Tropical Special Forest Trees and Ornamental Plants (Ministry of Education), College of Forestry, Hainan University, Haikou 570228, China
| | - Junxuan Huang
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen 518083, China
| | - Qing Wang
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen 518083, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yixin Zhu
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen 518083, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Li Wang
- College of Life Sciences, Zhejiang University, Hangzhou 310058, China
| | - Yanchun Xu
- BGI Life Science Joint Research Center, Northeast Forestry University, Harbin 150040, China
- College of Wildlife and Protected Area, Northeast Forestry University, Harbin 150040, China
- Corresponding author
| | - Chuyu Lin
- Shenzhen Zhong Nong Jing Yue Biotech Company Limited, Shenzhen 518120, China
- Corresponding author
| | - Huan Liu
- BGI Life Science Joint Research Center, Northeast Forestry University, Harbin 150040, China
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen 518083, China
- Guangdong Provincial Key Laboratory of Genome Read and Write, BGI-Shenzhen, Shenzhen 518120, China
- Corresponding author
| | - Zhijun Hou
- BGI Life Science Joint Research Center, Northeast Forestry University, Harbin 150040, China
- College of Wildlife and Protected Area, Northeast Forestry University, Harbin 150040, China
- Corresponding author
| |
Collapse
|
21
|
Bhattarai UR, Katuwal M, Poulin R, Gemmell NJ, Dowle E. Genome assembly and annotation of the European earwig Forficula auricularia (subspecies B). G3 (BETHESDA, MD.) 2022; 12:jkac199. [PMID: 35972389 PMCID: PMC9526046 DOI: 10.1093/g3journal/jkac199] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Accepted: 07/06/2022] [Indexed: 11/14/2022]
Abstract
The European earwig Forficula auricularia is an important model for studies of maternal care, sexual selection, sociality, and host-parasite interactions. However, detailed genetic investigations of this species are hindered by a lack of genomic resources. Here, we present a high-quality hybrid genome assembly for Forficula auricularia using Nanopore long-reads and 10× linked-reads. The final assembly is 1.06 Gb in length with 31.03% GC content. It consists of 919 scaffolds with an N50 of 12.55 Mb. Half of the genome is present in only 20 scaffolds. Benchmarking Universal Single-Copy Orthologs scores are ∼90% from 3 sets of single-copy orthologs (eukaryotic, insect, and arthropod). The total repeat elements in the genome are 64.62%. The MAKER2 pipeline annotated 12,876 protein-coding genes and 21,031 mRNAs. Phylogenetic analysis revealed the assembled genome as that of species B, one of the 2 known genetic subspecies of Forficula auricularia. The genome assembly, annotation, and associated resources will be of high value to a large and diverse group of researchers working on dermapterans.
Collapse
Affiliation(s)
| | - Mandira Katuwal
- Department of Anatomy, University of Otago, Dunedin 9016, New Zealand
| | - Robert Poulin
- Department of Zoology, University of Otago, Dunedin 9016, New Zealand
| | - Neil J Gemmell
- Department of Anatomy, University of Otago, Dunedin 9016, New Zealand
| | - Eddy Dowle
- Department of Anatomy, University of Otago, Dunedin 9016, New Zealand
| |
Collapse
|
22
|
Ogaji YO, Lee RC, Sawbridge TI, Cocks BG, Daetwyler HD, Kaur S. De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea. J Fungi (Basel) 2022; 8:884. [PMID: 36012871 PMCID: PMC9410150 DOI: 10.3390/jof8080884] [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: 06/22/2022] [Revised: 08/12/2022] [Accepted: 08/15/2022] [Indexed: 11/16/2022] Open
Abstract
Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga's genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host-pathogen interactions.
Collapse
Affiliation(s)
- Yvonne O. Ogaji
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, 5 Ring Road, Melbourne, VIC 3083, Australia
- School of Applied Systems Biology, La Trobe University, Melbourne, VIC 3086, Australia
| | - Robert C. Lee
- Centre for Crop and Disease Management, School of Molecular and Life Sciences, Curtin University, Perth, WA 6102, Australia
| | - Tim I. Sawbridge
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, 5 Ring Road, Melbourne, VIC 3083, Australia
- School of Applied Systems Biology, La Trobe University, Melbourne, VIC 3086, Australia
| | - Benjamin G. Cocks
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, 5 Ring Road, Melbourne, VIC 3083, Australia
- School of Applied Systems Biology, La Trobe University, Melbourne, VIC 3086, Australia
| | - Hans D. Daetwyler
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, 5 Ring Road, Melbourne, VIC 3083, Australia
- School of Applied Systems Biology, La Trobe University, Melbourne, VIC 3086, Australia
| | - Sukhjiwan Kaur
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, 5 Ring Road, Melbourne, VIC 3083, Australia
| |
Collapse
|
23
|
A Candidate Gene Cluster for the Bioactive Natural Product Gyrophoric Acid in Lichen-Forming Fungi. Microbiol Spectr 2022; 10:e0010922. [PMID: 35867425 PMCID: PMC9430680 DOI: 10.1128/spectrum.00109-22] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Natural products of lichen-forming fungi are structurally diverse and have a variety of medicinal properties. Despite this, they have limited implementation in industry mostly because the corresponding genes are unknown for most of their natural products. Here, we implement a long-read sequencing and bioinformatic approach to identify the putative biosynthetic gene cluster of the bioactive natural product gyrophoric acid (GA). Using 15 high-quality genomes representing nine GA-producing species of the lichen-forming fungal genus Umbilicaria, we identify the most likely GA cluster and investigate the cluster gene organization and composition across the nine species. Our results show that GA clusters are promiscuous within Umbilicaria, and only three genes are conserved across species, including the polyketide synthase (PKS) gene. In addition, our results suggest that the same cluster codes for different, but structurally similar compounds, namely, GA, umbilicaric-, and hiascic acid, bringing new evidence that lichen metabolite diversity is also generated through regulatory mechanisms at the molecular level. Ours is the first study to identify the most likely GA cluster and, thus, provides essential information to open new avenues for biotechnological approaches to producing and modifying GA and similar lichen-derived compounds. GA PKS is the first tridepside PKS to be identified. IMPORTANCE The implementation of natural products in the pharmaceutical industry relies on the possibility of modifying the natural product (NP) pathway to optimize yields and pharmacological effects. Characterization of genes and pathways underlying natural product biosynthesis is a major bottleneck for exploiting the medicinal properties of the natural products. Genome mining is a promising and relatively cost- and time-effective approach to utilize unexplored NP resources for drug discovery. In this study, we identify the most likely gene cluster for the lichen-forming fungal depside gyrophoric acid in nine Umbilicaria species. This compound shows cytotoxic and antiproliferative properties against several cancer cell lines and is also a broad-spectrum antimicrobial agent. This information paves the way for generating GA analogs with modified properties by selective activation/deactivation of genes.
Collapse
|
24
|
The role of zinc in the adaptive evolution of polar phytoplankton. Nat Ecol Evol 2022; 6:965-978. [PMID: 35654896 DOI: 10.1038/s41559-022-01750-x] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2021] [Accepted: 03/28/2022] [Indexed: 12/20/2022]
Abstract
Zinc is an essential trace metal for oceanic primary producers with the highest concentrations in polar oceans. However, its role in the biological functioning and adaptive evolution of polar phytoplankton remains enigmatic. Here, we have applied a combination of evolutionary genomics, quantitative proteomics, co-expression analyses and cellular physiology to suggest that model polar phytoplankton species have a higher demand for zinc because of elevated cellular levels of zinc-binding proteins. We propose that adaptive expansion of regulatory zinc-finger protein families, co-expanded and co-expressed zinc-binding proteins families involved in photosynthesis and growth in these microalgal species and their natural communities were identified to be responsible for the higher zinc demand. The expression of their encoding genes in eukaryotic phytoplankton metatranscriptomes from pole-to-pole was identified to correlate not only with dissolved zinc concentrations in the upper ocean but also with temperature, suggesting that environmental conditions of polar oceans are responsible for an increased demand of zinc. These results suggest that zinc plays an important role in supporting photosynthetic growth in eukaryotic polar phytoplankton and that this has been critical for algal colonization of low-temperature polar oceans.
Collapse
|
25
|
Li M, Li LM. RegScaf: a regression approach to scaffolding. Bioinformatics 2022; 38:2675-2682. [PMID: 35561180 PMCID: PMC9326850 DOI: 10.1093/bioinformatics/btac174] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2021] [Revised: 02/19/2022] [Accepted: 03/23/2022] [Indexed: 11/13/2022] Open
Abstract
MOTIVATION Crucial to the correctness of a genome assembly is the accuracy of the underlying scaffolds that specify the orders and orientations of contigs together with the gap distances between contigs. The current methods construct scaffolds based on the alignments of 'linking' reads against contigs. We found that some 'optimal' alignments are mistaken due to factors such as the contig boundary effect, particularly in the presence of repeats. Occasionally, the incorrect alignments can even overwhelm the correct ones. The detection of the incorrect linking information is challenging in any existing methods. RESULTS In this study, we present a novel scaffolding method RegScaf. It first examines the distribution of distances between contigs from read alignment by the kernel density. When multiple modes are shown in a density, orientation-supported links are grouped into clusters, each of which defines a linking distance corresponding to a mode. The linear model parameterizes contigs by their positions on the genome; then each linking distance between a pair of contigs is taken as an observation on the difference of their positions. The parameters are estimated by minimizing a global loss function, which is a version of trimmed sum of squares. The least trimmed squares estimate has such a high breakdown value that it can automatically remove the mistaken linking distances. The results on both synthetic and real datasets demonstrate that RegScaf outperforms some popular scaffolders, especially in the accuracy of gap estimates by substantially reducing extremely abnormal errors. Its strength in resolving repeat regions is exemplified by a real case. Its adaptability to large genomes and TGS long reads is validated as well. AVAILABILITY AND IMPLEMENTATION RegScaf is publicly available at https://github.com/lemontealala/RegScaf.git. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Mengtian Li
- National Center of Mathematics and Interdisciplinary Sciences, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Lei M Li
- To whom correspondence should be addressed.
| |
Collapse
|
26
|
Hoopes G, Meng X, Hamilton JP, Achakkagari SR, de Alves Freitas Guesdes F, Bolger ME, Coombs JJ, Esselink D, Kaiser NR, Kodde L, Kyriakidou M, Lavrijssen B, van Lieshout N, Shereda R, Tuttle HK, Vaillancourt B, Wood JC, de Boer JM, Bornowski N, Bourke P, Douches D, van Eck HJ, Ellis D, Feldman MJ, Gardner KM, Hopman JCP, Jiang J, De Jong WS, Kuhl JC, Novy RG, Oome S, Sathuvalli V, Tan EH, Ursum RA, Vales MI, Vining K, Visser RGF, Vossen J, Yencho GC, Anglin NL, Bachem CWB, Endelman JB, Shannon LM, Strömvik MV, Tai HH, Usadel B, Buell CR, Finkers R. Phased, chromosome-scale genome assemblies of tetraploid potato reveal a complex genome, transcriptome, and predicted proteome landscape underpinning genetic diversity. MOLECULAR PLANT 2022; 15:520-536. [PMID: 35026436 DOI: 10.1016/j.molp.2022.01.003] [Citation(s) in RCA: 65] [Impact Index Per Article: 32.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/19/2021] [Revised: 10/19/2021] [Accepted: 01/07/2022] [Indexed: 05/11/2023]
Abstract
Cultivated potato is a clonally propagated autotetraploid species with a highly heterogeneous genome. Phased assemblies of six cultivars including two chromosome-scale phased genome assemblies revealed extensive allelic diversity, including altered coding and transcript sequences, preferential allele expression, and structural variation that collectively result in a highly complex transcriptome and predicted proteome, which are distributed across the homologous chromosomes. Wild species contribute to the extensive allelic diversity in tetraploid cultivars, demonstrating ancestral introgressions predating modern breeding efforts. As a clonally propagated autotetraploid that undergoes limited meiosis, dysfunctional and deleterious alleles are not purged in tetraploid potato. Nearly a quarter of the loci bore mutations are predicted to have a high negative impact on protein function, complicating breeder's efforts to reduce genetic load. The StCDF1 locus controls maturity, and analysis of six tetraploid genomes revealed that 12 allelic variants of StCDF1 are correlated with maturity in a dosage-dependent manner. Knowledge of the complexity of the tetraploid potato genome with its rampant structural variation and embedded deleterious and dysfunctional alleles will be key not only to implementing precision breeding of tetraploid cultivars but also to the construction of homozygous, diploid potato germplasm containing favorable alleles to capitalize on heterosis in F1 hybrids.
Collapse
Affiliation(s)
- Genevieve Hoopes
- Department of Plant Biology, Michigan State University, East Lansing, MI 48824, USA
| | - Xiaoxi Meng
- Department of Horticultural Science, University of Minnesota, St. Paul, MN 55108, USA
| | - John P Hamilton
- Department of Plant Biology, Michigan State University, East Lansing, MI 48824, USA
| | - Sai Reddy Achakkagari
- Department of Plant Science, McGill University, Sainte-Anne-de-Bellevue, QC H9X 3V9, Canada
| | | | - Marie E Bolger
- IBG-4 Bioinformatics, Forschungszentrum Jülich, Wilhelm Johnen Str, 52428 Jülich, Germany
| | - Joseph J Coombs
- Department of Plant, Soil, and Microbial Sciences, Michigan State University, East Lansing, MI 48824, USA
| | - Danny Esselink
- Plant Breeding, Wageningen University & Research, Plant Breeding, 6708 PB Wageningen, the Netherlands
| | - Natalie R Kaiser
- Department of Plant, Soil, and Microbial Sciences, Michigan State University, East Lansing, MI 48824, USA; Bayer Crop Science, Woodland, CA 95695, USA
| | - Linda Kodde
- Plant Breeding, Wageningen University & Research, Plant Breeding, 6708 PB Wageningen, the Netherlands
| | - Maria Kyriakidou
- Department of Plant Science, McGill University, Sainte-Anne-de-Bellevue, QC H9X 3V9, Canada
| | - Brian Lavrijssen
- Plant Breeding, Wageningen University & Research, Plant Breeding, 6708 PB Wageningen, the Netherlands
| | - Natascha van Lieshout
- Plant Breeding, Wageningen University & Research, Plant Breeding, 6708 PB Wageningen, the Netherlands
| | - Rachel Shereda
- Department of Plant Biology, Michigan State University, East Lansing, MI 48824, USA
| | - Heather K Tuttle
- Department of Horticultural Science, University of Minnesota, St. Paul, MN 55108, USA
| | | | - Joshua C Wood
- Department of Plant Biology, Michigan State University, East Lansing, MI 48824, USA
| | | | - Nolan Bornowski
- Department of Plant Biology, Michigan State University, East Lansing, MI 48824, USA
| | - Peter Bourke
- Plant Breeding, Wageningen University & Research, Plant Breeding, 6708 PB Wageningen, the Netherlands
| | - David Douches
- Department of Plant, Soil, and Microbial Sciences, Michigan State University, East Lansing, MI 48824, USA
| | - Herman J van Eck
- Plant Breeding, Wageningen University & Research, Plant Breeding, 6708 PB Wageningen, the Netherlands
| | - Dave Ellis
- International Potato Center, 1895 Avenida La Molina, Lima, Peru
| | | | - Kyle M Gardner
- Agriculture and Agri-Food Canada Fredericton Research and Development Centre, Fredericton, NB E3B 4Z7, Canada
| | | | - Jiming Jiang
- Department of Plant Biology, Michigan State University, East Lansing, MI 48824, USA; Department of Horticulture, Michigan State University, East Lansing, MI 48824, USA
| | - Walter S De Jong
- School of Integrative Plant Science, Cornell University, Ithaca, NY 14853-1901, USA
| | - Joseph C Kuhl
- Department of Plant Sciences, University of Idaho, Moscow, ID 83844, USA
| | - Richard G Novy
- USDA-ARS, Small Grains and Potato Germplasm Research, Aberdeen, ID 83210, USA
| | - Stan Oome
- HZPC Research B.V., Edisonweg 5, 8501 XG Joure, the Netherlands
| | - Vidyasagar Sathuvalli
- Department of Crop and Soil Science, Oregon State University, Hermiston, OR 97838, USA
| | - Ek Han Tan
- School of Biology and Ecology, University of Maine, 5735 Hitchner Hall Orono, ME 04469, USA
| | - Remco A Ursum
- HZPC Research B.V., Edisonweg 5, 8501 XG Joure, the Netherlands
| | - M Isabel Vales
- Department of Horticultural Sciences, Texas A&M University, College Station, TX 77843-2133, USA
| | - Kelly Vining
- Department of Horticulture, Oregon State University, Corvallis, OR 97331, USA
| | - Richard G F Visser
- Plant Breeding, Wageningen University & Research, Plant Breeding, 6708 PB Wageningen, the Netherlands
| | - Jack Vossen
- Plant Breeding, Wageningen University & Research, Plant Breeding, 6708 PB Wageningen, the Netherlands
| | - G Craig Yencho
- Department of Horticultural Science, North Carolina State University, Raleigh, NC 27695-7609, USA
| | - Noelle L Anglin
- International Potato Center, 1895 Avenida La Molina, Lima, Peru; USDA-ARS, Small Grains and Potato Germplasm Research, Aberdeen, ID 83210, USA
| | - Christian W B Bachem
- Plant Breeding, Wageningen University & Research, Plant Breeding, 6708 PB Wageningen, the Netherlands
| | - Jeffrey B Endelman
- Department of Horticulture, University of Wisconsin-Madison, Madison, WI 53706, USA
| | - Laura M Shannon
- Department of Horticultural Science, University of Minnesota, St. Paul, MN 55108, USA
| | - Martina V Strömvik
- Department of Plant Science, McGill University, Sainte-Anne-de-Bellevue, QC H9X 3V9, Canada
| | - Helen H Tai
- Agriculture and Agri-Food Canada Fredericton Research and Development Centre, Fredericton, NB E3B 4Z7, Canada
| | - Björn Usadel
- IBG-4 Bioinformatics, Forschungszentrum Jülich, Wilhelm Johnen Str, 52428 Jülich, Germany; Institute for Biological Data Science, Heinrich Heine University, Düsseldorf, 40225 Düsseldorf, Germany
| | - C Robin Buell
- Department of Plant Biology, Michigan State University, East Lansing, MI 48824, USA; Plant Resilience Institute, Michigan State University, East Lansing, MI 48824, USA.
| | - Richard Finkers
- Plant Breeding, Wageningen University & Research, Plant Breeding, 6708 PB Wageningen, the Netherlands.
| |
Collapse
|
27
|
Delorme Q, Costa R, Mansour Y, Fiston-Lavier AS, Chateau A. Involving repetitive regions in scaffolding improvement. J Bioinform Comput Biol 2021; 19:2140016. [PMID: 34923926 DOI: 10.1142/s0219720021400163] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
In this paper, we investigate througth a premilinary study the influence of repeat elements during the assembly process. We analyze the link between the presence and the nature of one type of repeat element, called transposable element (TE) and misassembly events in genome assemblies. We propose to improve assemblies by taking into account the presence of repeat elements, including TEs, during the scaffolding step. We analyze the results and relate the misassemblies to TEs before and after correction.
Collapse
Affiliation(s)
- Quentin Delorme
- LIRMM, Univ Montpellier, CNRS, Montpellier, France.,Laboratoire MIVEGEC (Université de Montpellier, CNRS 5290, IRD 229), Centre de Recherche en Écologie et Évolution de la Santé (CREES), Institut de Recherche pour le Développement (IRD), F-34394, Montpellier, France
| | - Rémy Costa
- LIRMM, Univ Montpellier, CNRS, Montpellier, France.,IGH-UMR9002, Univ Montpellier, CNRS, Montpellier, France
| | - Yasmine Mansour
- LIRMM, Univ Montpellier, CNRS, Montpellier, France.,ISEM, Univ Montpellier, CNRS, IRD, Montpellier, France
| | - Anna-Sophie Fiston-Lavier
- ISEM, Univ Montpellier, CNRS, IRD, Montpellier, France.,Institut Universitaire de France (IUF), France
| | | |
Collapse
|
28
|
Coombe L, Li JX, Lo T, Wong J, Nikolic V, Warren RL, Birol I. LongStitch: high-quality genome assembly correction and scaffolding using long reads. BMC Bioinformatics 2021; 22:534. [PMID: 34717540 PMCID: PMC8557608 DOI: 10.1186/s12859-021-04451-7] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2021] [Accepted: 10/19/2021] [Indexed: 12/12/2022] Open
Abstract
BACKGROUND Generating high-quality de novo genome assemblies is foundational to the genomics study of model and non-model organisms. In recent years, long-read sequencing has greatly benefited genome assembly and scaffolding, a process by which assembled sequences are ordered and oriented through the use of long-range information. Long reads are better able to span repetitive genomic regions compared to short reads, and thus have tremendous utility for resolving problematic regions and helping generate more complete draft assemblies. Here, we present LongStitch, a scalable pipeline that corrects and scaffolds draft genome assemblies exclusively using long reads. RESULTS LongStitch incorporates multiple tools developed by our group and runs in up to three stages, which includes initial assembly correction (Tigmint-long), followed by two incremental scaffolding stages (ntLink and ARKS-long). Tigmint-long and ARKS-long are misassembly correction and scaffolding utilities, respectively, previously developed for linked reads, that we adapted for long reads. Here, we describe the LongStitch pipeline and introduce our new long-read scaffolder, ntLink, which utilizes lightweight minimizer mappings to join contigs. LongStitch was tested on short and long-read assemblies of Caenorhabditis elegans, Oryza sativa, and three different human individuals using corresponding nanopore long-read data, and improves the contiguity of each assembly from 1.2-fold up to 304.6-fold (as measured by NGA50 length). Furthermore, LongStitch generates more contiguous and correct assemblies compared to state-of-the-art long-read scaffolder LRScaf in most tests, and consistently improves upon human assemblies in under five hours using less than 23 GB of RAM. CONCLUSIONS Due to its effectiveness and efficiency in improving draft assemblies using long reads, we expect LongStitch to benefit a wide variety of de novo genome assembly projects. The LongStitch pipeline is freely available at https://github.com/bcgsc/longstitch .
Collapse
Affiliation(s)
- Lauren Coombe
- Canada's Michael Smith Genome Sciences Centre, BC Cancer Research, 100-570 West 7th Avenue, Vancouver, BC, V5Z 4S6, Canada.
| | - Janet X Li
- Canada's Michael Smith Genome Sciences Centre, BC Cancer Research, 100-570 West 7th Avenue, Vancouver, BC, V5Z 4S6, Canada
| | - Theodora Lo
- Canada's Michael Smith Genome Sciences Centre, BC Cancer Research, 100-570 West 7th Avenue, Vancouver, BC, V5Z 4S6, Canada
| | - Johnathan Wong
- Canada's Michael Smith Genome Sciences Centre, BC Cancer Research, 100-570 West 7th Avenue, Vancouver, BC, V5Z 4S6, Canada
| | - Vladimir Nikolic
- Canada's Michael Smith Genome Sciences Centre, BC Cancer Research, 100-570 West 7th Avenue, Vancouver, BC, V5Z 4S6, Canada
| | - René L Warren
- Canada's Michael Smith Genome Sciences Centre, BC Cancer Research, 100-570 West 7th Avenue, Vancouver, BC, V5Z 4S6, Canada
| | - Inanc Birol
- Canada's Michael Smith Genome Sciences Centre, BC Cancer Research, 100-570 West 7th Avenue, Vancouver, BC, V5Z 4S6, Canada
| |
Collapse
|
29
|
Kawato S, Nishitsuji K, Arimoto A, Hisata K, Kawamitsu M, Nozaki R, Kondo H, Shinzato C, Ohira T, Satoh N, Shoguchi E, Hirono I. Genome and transcriptome assemblies of the kuruma shrimp, Marsupenaeus japonicus. G3 (BETHESDA, MD.) 2021; 11:jkab268. [PMID: 34515781 PMCID: PMC8527471 DOI: 10.1093/g3journal/jkab268] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/07/2021] [Accepted: 07/18/2021] [Indexed: 11/25/2022]
Abstract
The kuruma shrimp Marsupenaeus japonicus (order Decapoda, family Penaeidae) is an economically important crustacean that occurs in shallow, warm seas across the Indo-Pacific. Here, using a combination of Illumina and Oxford Nanopore Technologies platforms, we produced a draft genome assembly of M. japonicus (1.70 Gbp; 18,210 scaffolds; scaffold N50 = 234.9 kbp; 34.38% GC, 93.4% BUSCO completeness) and a complete mitochondrial genome sequence (15,969 bp). As with other penaeid shrimp genomes, the M. japonicus genome is extremely rich in simple repeats, which occupies 27.4% of the assembly. A total of 26,381 protein-coding gene models (94.7% BUSCO completeness) were predicted, of which 18,005 genes (68.2%) were assigned functional description by at least one method. We also produced an Illumina-based transcriptome shotgun assembly (40,991 entries; 93.0% BUSCO completeness) and a PacBio Iso-Seq transcriptome assembly (25,415 entries; 67.5% BUSCO completeness). We envision that the M. japonicus genome and transcriptome assemblies will serve as useful resources for the basic research, fisheries management, and breeding programs of M. japonicus.
Collapse
Affiliation(s)
- Satoshi Kawato
- Laboratory of Genome Science, Tokyo University of Marine Science and Technology, Tokyo 108-8477, Japan
| | - Koki Nishitsuji
- Marine Genomics Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Asuka Arimoto
- Marine Genomics Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Kanako Hisata
- Marine Genomics Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Mayumi Kawamitsu
- DNA Sequencing Section, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Reiko Nozaki
- Laboratory of Genome Science, Tokyo University of Marine Science and Technology, Tokyo 108-8477, Japan
| | - Hidehiro Kondo
- Laboratory of Genome Science, Tokyo University of Marine Science and Technology, Tokyo 108-8477, Japan
| | - Chuya Shinzato
- Atmosphere and Ocean Research Institute, The University of Tokyo, Chiba 277-0882, Japan
| | - Tsuyoshi Ohira
- Faculty of Science, Department of Biological Sciences, Kanagawa University, Kanagawa 221-8686, Japan
| | - Noriyuki Satoh
- Marine Genomics Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Eiichi Shoguchi
- Marine Genomics Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa 904-0495, Japan
| | - Ikuo Hirono
- Laboratory of Genome Science, Tokyo University of Marine Science and Technology, Tokyo 108-8477, Japan
| |
Collapse
|
30
|
Duong TA, Aylward J, Ametrano CG, Poudel B, Santana QC, Wilken PM, Martin A, Arun-Chinnappa KS, de Vos L, DiStefano I, Grewe F, Huhndorf S, Lumbsch HT, Rakoma JR, Poudel B, Steenkamp ET, Sun Y, van der Nest MA, Wingfield MJ, Yilmaz N, Wingfield BD. IMA Genome - F15 : Draft genome assembly of Fusarium pilosicola, Meredithiella fracta, Niebla homalea, Pyrenophora teres hybrid WAC10721, and Teratosphaeria viscida. IMA Fungus 2021; 12:30. [PMID: 34645521 PMCID: PMC8513234 DOI: 10.1186/s43008-021-00077-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/05/2021] [Indexed: 12/01/2022] Open
Affiliation(s)
- Tuan Anh Duong
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Private Bag X20, Hatfield, 0028, South Africa
| | - Janneke Aylward
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Private Bag X20, Hatfield, 0028, South Africa
- Department of Conservation Ecology and Entomology, Stellenbosch University, Private Bag X1, Matieland, 7602, South Africa
| | - Claudio Gennaro Ametrano
- Field Museum, Department of Science and Education, Grainger Bioinformatics Center, 1400 S Lake Shore Drive, Chicago, IL 60605, USA
| | - Barsha Poudel
- Centre for Crop Health, University of Southern Queensland, Toowoomba, QLD, 4350, Australia
| | - Quentin Carlo Santana
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Private Bag X20, Hatfield, 0028, South Africa
| | - Pieter Markus Wilken
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Private Bag X20, Hatfield, 0028, South Africa
| | - Anke Martin
- Centre for Crop Health, University of Southern Queensland, Toowoomba, QLD, 4350, Australia
| | - Kiruba Shankari Arun-Chinnappa
- Centre for Crop Health, University of Southern Queensland, Toowoomba, QLD, 4350, Australia
- PerkinElmer Pty LTD., Level 2, Building 5, Brandon Business Park 530-540, Springvale Road, Glen Waverley, VIC, 3150, Australia
| | - Lieschen de Vos
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Private Bag X20, Hatfield, 0028, South Africa
| | - Isabel DiStefano
- Field Museum, Department of Science and Education, Grainger Bioinformatics Center, 1400 S Lake Shore Drive, Chicago, IL 60605, USA
| | - Felix Grewe
- Field Museum, Department of Science and Education, Grainger Bioinformatics Center, 1400 S Lake Shore Drive, Chicago, IL 60605, USA
| | - Sabine Huhndorf
- Field Museum, Department of Science and Education, Grainger Bioinformatics Center, 1400 S Lake Shore Drive, Chicago, IL 60605, USA
| | - Helge Thorsten Lumbsch
- Field Museum, Department of Science and Education, Grainger Bioinformatics Center, 1400 S Lake Shore Drive, Chicago, IL 60605, USA
| | - Jostina Raesetsa Rakoma
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Private Bag X20, Hatfield, 0028, South Africa
| | - Barsha Poudel
- Centre for Crop Health, University of Southern Queensland, Toowoomba, QLD, 4350, Australia
| | - Emma Theodora Steenkamp
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Private Bag X20, Hatfield, 0028, South Africa
| | - Yukun Sun
- Field Museum, Department of Science and Education, Grainger Bioinformatics Center, 1400 S Lake Shore Drive, Chicago, IL 60605, USA
| | - Magriet A van der Nest
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Private Bag X20, Hatfield, 0028, South Africa
- Biotechnology Platform, Agricultural Research Council, Onderstepoort, Pretoria, 0110, South Africa
| | - Michael John Wingfield
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Private Bag X20, Hatfield, 0028, South Africa
| | - Neriman Yilmaz
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Private Bag X20, Hatfield, 0028, South Africa
| | - Brenda Diana Wingfield
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Private Bag X20, Hatfield, 0028, South Africa.
| |
Collapse
|
31
|
El-Bondkly EAM, El-Bondkly AAM, El-Bondkly AAM. Marine endophytic fungal metabolites: A whole new world of pharmaceutical therapy exploration. Heliyon 2021; 7:e06362. [PMID: 33869822 PMCID: PMC8035529 DOI: 10.1016/j.heliyon.2021.e06362] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2020] [Revised: 04/16/2020] [Accepted: 02/22/2021] [Indexed: 12/19/2022] Open
Abstract
The growing threat arises due to diseases such as cancer and the infections around the world leading to a critical requirement for novel and constructive compounds with unique ways of action capable of combating these deadly diseases. At present, it is evident that endophytic fungi constitute an enormous as well as comparatively untapped source of great biodiversity that can be considered as a wellspring of effective novel natural products for medical, agricultural and industrial use. Marine endophytic fungi have been found in every marine plants (algae, seagrass, driftwood, mangrove plants), marine vertebrates (mainly, fish) or marine invertebrates (mainly, sponge and coral) inter- and intra-cellular without causing any palpable symptoms of illness. Since evolution of microbes and eukaryotes to a higher level, coevolution has resulted in specific interaction mechanisms. Endophytic fungi are known to influence the life cycle and are necessary for the homeostasis of their eukaryotic hosts and the chemical signals of their host have been shown to activate gene expression in endophytes to induce expression of endophytic secondary metabolites. Marine endophytic fungi are receiving increasing attention by chemists because of their varied and structurally unmatched compounds that have strong biological roles in life as lead pharmaceutical compounds, including anticancer, antiviral, insulin mimetic, antineurodegenerative, antimicrobial, antioxidant and immuno-suppressant compounds. Moreover, fungal endophytes proved to have different biological activities for exploitation in the environmental and agricultural sustainability.
Collapse
|
32
|
Luo J, Wei Y, Lyu M, Wu Z, Liu X, Luo H, Yan C. A comprehensive review of scaffolding methods in genome assembly. Brief Bioinform 2021; 22:6149347. [PMID: 33634311 DOI: 10.1093/bib/bbab033] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2020] [Revised: 01/21/2021] [Accepted: 01/22/2021] [Indexed: 12/20/2022] Open
Abstract
In the field of genome assembly, scaffolding methods make it possible to obtain a more complete and contiguous reference genome, which is the cornerstone of genomic research. Scaffolding methods typically utilize the alignments between contigs and sequencing data (reads) to determine the orientation and order among contigs and to produce longer scaffolds, which are helpful for genomic downstream analysis. With the rapid development of high-throughput sequencing technologies, diverse types of reads have emerged over the past decade, especially in long-range sequencing, which have greatly enhanced the assembly quality of scaffolding methods. As the number of scaffolding methods increases, biology and bioinformatics researchers need to perform in-depth analyses of state-of-the-art scaffolding methods. In this article, we focus on the difficulties in scaffolding, the differences in characteristics among various kinds of reads, the methods by which current scaffolding methods address these difficulties, and future research opportunities. We hope this work will benefit the design of new scaffolding methods and the selection of appropriate scaffolding methods for specific biological studies.
Collapse
Affiliation(s)
- Junwei Luo
- College of Computer Science and Technology, Henan Polytechnic University, Jiaozuo, China
| | - Yawei Wei
- College of Computer Science and Technology, Henan Polytechnic University, Jiaozuo, China
| | - Mengna Lyu
- College of Computer Science and Technology, Henan Polytechnic University, Jiaozuo, China
| | - Zhengjiang Wu
- College of Computer Science and Technology, Henan Polytechnic University, Jiaozuo, China
| | - Xiaoyan Liu
- College of Computer Science and Technology, Henan Polytechnic University, Jiaozuo, China
| | - Huimin Luo
- School of Computer and Information Engineering, Henan University, Kaifeng, China
| | - Chaokun Yan
- School of Computer and Information Engineering, Henan University, Kaifeng, China
| |
Collapse
|
33
|
Zhao M, Wang W, Zhang F, Ma C, Liu Z, Yang MH, Chen W, Li Q, Cui M, Jiang K, Feng C, Li JT, Ma L. A chromosome-level genome of the mud crab (Scylla paramamosain estampador) provides insights into the evolution of chemical and light perception in this crustacean. Mol Ecol Resour 2021; 21:1299-1317. [PMID: 33464679 DOI: 10.1111/1755-0998.13332] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2020] [Revised: 01/09/2021] [Accepted: 01/13/2021] [Indexed: 02/06/2023]
Abstract
Mud crabs, found throughout the Indo-Pacific region, are coastal species that are important fisheries resources in many tropical and subtropical Asian countries. Here, we present a chromosome-level genome assembly of a mud crab (Scylla paramamosain). The genome is 1.55 Gb (contig N50 191 kb) in length and encodes 17,821 proteins. The heterozygosity of the assembled genome was estimated to be 0.47%. Effective population size analysis suggested that an initial large population size of this species was maintained until 200 thousand years ago. The contraction of cuticle protein and opsin genes compared with Litopenaeus vannamei is assumed to be correlated with shell hardness and light perception ability, respectively. Furthermore, the analysis of three chemoreceptor gene families, the odorant receptor (OR), gustatory receptor (GR) and ionotropic receptor (IR) families, suggested that the mud crab has no OR genes and shows a contraction of GR genes and expansion of IR genes. The numbers of the three gene families were similar to those in three other decapods but different from those in two nondecapods and insects. In addition, IRs were more diversified in decapods than in nondecapod crustaceans, and most of the expanded IRs in the mud crab genome were clustered with the antennal IR clades. These findings suggested that IRs might exhibit more diverse functions in decapods than in nondecapods, which may compensate for the smaller number of GR genes. Decoding the S. paramamosain genome not only provides insight into the genetic changes underpinning ecological traits but also provides valuable information for improving the breeding and aquaculture of this species.
Collapse
Affiliation(s)
- Ming Zhao
- Key Laboratory of East China Sea Fishery Resources Exploitation, Ministry of Agriculture and Rural Affairs, East China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Shanghai, China.,College of Fisheries and Life Science, Shanghai Ocean University, Shanghai, China
| | - Wei Wang
- Key Laboratory of East China Sea Fishery Resources Exploitation, Ministry of Agriculture and Rural Affairs, East China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Shanghai, China
| | - Fengying Zhang
- Key Laboratory of East China Sea Fishery Resources Exploitation, Ministry of Agriculture and Rural Affairs, East China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Shanghai, China
| | - Chunyan Ma
- Key Laboratory of East China Sea Fishery Resources Exploitation, Ministry of Agriculture and Rural Affairs, East China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Shanghai, China
| | - Zhiqiang Liu
- Key Laboratory of East China Sea Fishery Resources Exploitation, Ministry of Agriculture and Rural Affairs, East China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Shanghai, China
| | - Meidi-Huang Yang
- Key Laboratory of Aquatic Genomics, Ministry of Agriculture and Rural Affairs, Beijing Key Laboratory of Fishery Biotechnology, Chinese Academy of Fishery Sciences, Beijing, China
| | - Wei Chen
- Key Laboratory of East China Sea Fishery Resources Exploitation, Ministry of Agriculture and Rural Affairs, East China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Shanghai, China
| | - Qingsong Li
- Key Laboratory of Aquatic Genomics, Ministry of Agriculture and Rural Affairs, Beijing Key Laboratory of Fishery Biotechnology, Chinese Academy of Fishery Sciences, Beijing, China
| | - Mingshu Cui
- Key Laboratory of Aquatic Genomics, Ministry of Agriculture and Rural Affairs, Beijing Key Laboratory of Fishery Biotechnology, Chinese Academy of Fishery Sciences, Beijing, China
| | - Keji Jiang
- Key Laboratory of East China Sea Fishery Resources Exploitation, Ministry of Agriculture and Rural Affairs, East China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Shanghai, China
| | - Chunlei Feng
- Key Laboratory of East China Sea Fishery Resources Exploitation, Ministry of Agriculture and Rural Affairs, East China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Shanghai, China
| | - Jiong Tang Li
- Key Laboratory of Aquatic Genomics, Ministry of Agriculture and Rural Affairs, Beijing Key Laboratory of Fishery Biotechnology, Chinese Academy of Fishery Sciences, Beijing, China
| | - Lingbo Ma
- Key Laboratory of East China Sea Fishery Resources Exploitation, Ministry of Agriculture and Rural Affairs, East China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Shanghai, China.,College of Fisheries and Life Science, Shanghai Ocean University, Shanghai, China
| |
Collapse
|
34
|
Nguyen SH, Cao MD, Coin LJM. Real-time resolution of short-read assembly graph using ONT long reads. PLoS Comput Biol 2021; 17:e1008586. [PMID: 33471816 PMCID: PMC7850483 DOI: 10.1371/journal.pcbi.1008586] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2020] [Revised: 02/01/2021] [Accepted: 11/30/2020] [Indexed: 11/19/2022] Open
Abstract
A streaming assembly pipeline utilising real-time Oxford Nanopore Technology (ONT) sequencing data is important for saving sequencing resources and reducing time-to-result. A previous approach implemented in npScarf provided an efficient streaming algorithm for hybrid assembly but was relatively prone to mis-assemblies compared to other graph-based methods. Here we present npGraph, a streaming hybrid assembly tool using the assembly graph instead of the separated pre-assembly contigs. It is able to produce more complete genome assembly by resolving the path finding problem on the assembly graph using long reads as the traversing guide. Application to synthetic and real data from bacterial isolate genomes show improved accuracy while still maintaining a low computational cost. npGraph also provides a graphical user interface (GUI) which provides a real-time visualisation of the progress of assembly. The tool and source code is available at https://github.com/hsnguyen/assembly.
Collapse
Affiliation(s)
- Son Hoang Nguyen
- Institute for Molecular Bioscience, the University of Queensland, St Lucia, Brisbane, Australia
- * E-mail: (SHN); (LC)
| | - Minh Duc Cao
- Institute for Molecular Bioscience, the University of Queensland, St Lucia, Brisbane, Australia
| | - Lachlan J. M. Coin
- Institute for Molecular Bioscience, the University of Queensland, St Lucia, Brisbane, Australia
- Department of Microbiology and Immunology, The University of Melbourne, Parkville, Australia
- Department of Clinical Pathology, The University of Melbourne, Parkville, Australia
- Department of Infectious Disease, Imperial College London, London, UK
- * E-mail: (SHN); (LC)
| |
Collapse
|
35
|
Ubiquitous Selfish Toxin-Antidote Elements in Caenorhabditis Species. Curr Biol 2021; 31:990-1001.e5. [PMID: 33417886 DOI: 10.1016/j.cub.2020.12.013] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2020] [Revised: 10/28/2020] [Accepted: 12/09/2020] [Indexed: 11/22/2022]
Abstract
Toxin-antidote elements (TAs) are selfish genetic dyads that spread in populations by selectively killing non-carriers. TAs are common in prokaryotes, but very few examples are known in animals. Here, we report the discovery of maternal-effect TAs in both C. tropicalis and C. briggsae, two distant relatives of C. elegans. In C. tropicalis, multiple TAs combine to cause a striking degree of intraspecific incompatibility: five elements reduce the fitness of >70% of the F2 hybrid progeny of two Caribbean isolates. We identified the genes underlying one of the novel TAs, slow-1/grow-1, and found that its toxin, slow-1, is homologous to nuclear hormone receptors. Remarkably, although previously known TAs act during embryonic development, maternal loading of slow-1 in oocytes specifically slows down larval development, delaying the onset of reproduction by several days. Finally, we found that balancing selection acting on linked, conflicting TAs hampers their ability to spread in populations, leading to more stable genetic incompatibilities. Our findings indicate that TAs are widespread in Caenorhabditis species and target a wide range of developmental processes and that antagonism between them may cause lasting incompatibilities in natural populations. We expect that similar phenomena exist in other animal species.
Collapse
|
36
|
Batra SS, Levy-Sakin M, Robinson J, Guillory J, Durinck S, Vilgalys TP, Kwok PY, Cox LA, Seshagiri S, Song YS, Wall JD. Accurate assembly of the olive baboon (Papio anubis) genome using long-read and Hi-C data. Gigascience 2020; 9:giaa134. [PMID: 33283855 PMCID: PMC7719865 DOI: 10.1093/gigascience/giaa134] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2020] [Revised: 07/14/2020] [Accepted: 11/04/2020] [Indexed: 11/14/2022] Open
Abstract
BACKGROUND Baboons are a widely used nonhuman primate model for biomedical, evolutionary, and basic genetics research. Despite this importance, the genomic resources for baboons are limited. In particular, the current baboon reference genome Panu_3.0 is a highly fragmented, reference-guided (i.e., not fully de novo) assembly, and its poor quality inhibits our ability to conduct downstream genomic analyses. FINDINGS Here we present a de novo genome assembly of the olive baboon (Papio anubis) that uses data from several recently developed single-molecule technologies. Our assembly, Panubis1.0, has an N50 contig size of ∼1.46 Mb (as opposed to 139 kb for Panu_3.0) and has single scaffolds that span each of the 20 autosomes and the X chromosome. CONCLUSIONS We highlight multiple lines of evidence (including Bionano Genomics data, pedigree linkage information, and linkage disequilibrium data) suggesting that there are several large assembly errors in Panu_3.0, which have been corrected in Panubis1.0.
Collapse
Affiliation(s)
- Sanjit Singh Batra
- Computer Science Division, University of California Berkeley, 2626 Hearst Avenue, Berkeley, CA 94720, USA
| | - Michal Levy-Sakin
- Cardiovascular Research Institute, University of California San Francisco, 555 Mission Bay Boulevard South, San Francisco, CA 94158, USA
| | - Jacqueline Robinson
- Institute for Human Genetics, University of California San Francisco, 513 Parnassus Avenue, San Francisco, CA 94143, USA
| | - Joseph Guillory
- Department of Molecular Biology, Genentech, Inc., 1 DNA Way, South San Francisco, CA 94080, USA
| | - Steffen Durinck
- Department of Molecular Biology, Genentech, Inc., 1 DNA Way, South San Francisco, CA 94080, USA
- Bioinformatics and Computational Biology Department, Genentech, Inc., 1 DNA Way, South San Francisco, CA 94080, USA
| | - Tauras P Vilgalys
- Department of Evolutionary Anthropology, Duke University, 130 Science Drive, Durham, NC 27705, USA
| | - Pui-Yan Kwok
- Cardiovascular Research Institute, University of California San Francisco, 555 Mission Bay Boulevard South, San Francisco, CA 94158, USA
- Institute for Human Genetics, University of California San Francisco, 513 Parnassus Avenue, San Francisco, CA 94143, USA
| | - Laura A Cox
- Center for Precision Medicine, Department of Internal Medicine, Section of Molecular Medicine, Wake Forest School of Medicine, 475 Vine Drive, Winston-Salem, NC 27101, USA
- Southwest National Primate Research Center, Texas Biomedical Research Institute, 8715 W. Military Drive, San Antonio, TX 78227, USA
| | - Somasekar Seshagiri
- Department of Molecular Biology, Genentech, Inc., 1 DNA Way, South San Francisco, CA 94080, USA
| | - Yun S Song
- Computer Science Division, University of California Berkeley, 2626 Hearst Avenue, Berkeley, CA 94720, USA
- Department of Statistics, University of California Berkeley, 367 Evans Hall, Berkeley, CA 94720, USA
- Chan Zuckerberg Biohub, Mission Bay, San Francisco, CA 94158, USA
| | - Jeffrey D Wall
- Institute for Human Genetics, University of California San Francisco, 513 Parnassus Avenue, San Francisco, CA 94143, USA
| |
Collapse
|
37
|
Reva ON, Larisa SA, Mwakilili AD, Tibuhwa D, Lyantagaye S, Chan WY, Lutz S, Ahrens CH, Vater J, Borriss R. Complete genome sequence and epigenetic profile of Bacillus velezensis UCMB5140 used for plant and crop protection in comparison with other plant-associated Bacillus strains. Appl Microbiol Biotechnol 2020; 104:7643-7656. [PMID: 32651600 DOI: 10.1007/s00253-020-10767-w] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2020] [Revised: 06/11/2020] [Accepted: 07/02/2020] [Indexed: 01/29/2023]
Abstract
The application of biocontrol biopesticides based on plant growth-promoting rhizobacteria (PGPR), particularly members of the genus Bacillus, is considered a promising perspective to make agricultural practices sustainable and ecologically safe. Recent advances in genome sequencing by third-generation sequencing technologies, e.g., Pacific Biosciences' Single Molecule Real-Time (PacBio SMRT) platform, have allowed researchers to gain deeper insights into the molecular and genetic mechanisms of PGPR activities, and to compare whole genome sequences and global patterns of epigenetic modifications. In the current work, this approach was used to sequence and compare four Bacillus strains that exhibited various PGPR activities including the strain UCMB5140, which is used in the commercial biopesticide Phytosubtil. Whole genome comparison and phylogenomic inference assigned the strain UCMB5140 to the species Bacillus velezensis. Strong biocontrol activities of this strain were confirmed in several bioassays. Several factors that affect the evolution of active PGPR B. velezensis strains were identified: (1) horizontal acquisition of novel non-ribosomal peptide synthetases (NRPS) and adhesion genes; (2) rearrangements of functional modules of NRPS genes leading to strain specific combinations of their encoded products; (3) gain and loss of methyltransferases that can cause global alterations in DNA methylation patterns, which eventually may affect gene expression and regulate transcription. Notably, we identified a horizontally transferred NRPS operon encoding an uncharacterized polypeptide antibiotic in B. velezensis UCMB5140. Other horizontally acquired genes comprised a possible adhesin and a methyltransferase, which may explain the strain-specific methylation pattern of the chromosomal DNA of UCMB5140. KEY POINTS: • Whole genome sequence of the active PGPR Bacillus velezensis UCMB5140. • Identification of genetic determinants responsible for PGPR activities. • Role of methyltransferases and epigenetic mechanisms in evolution of bacteria.
Collapse
Affiliation(s)
- Oleg N Reva
- Centre for Bioinformatics and Computational Biology, Department of Biochemistry, Genetics and Microbiology, University of Pretoria, Hillcrest, Lynnwood Rd., Pretoria, South Africa.
| | - Safronova A Larisa
- Innovation and Technology Transfer Laboratory, DK Zabolotny Institute of Microbiology and Virology, 154 Zabolotnogo Str, Kyiv, 03143, Ukraine
| | - Aneth D Mwakilili
- Department of Molecular Biology and Biotechnology, University of Dar es Salaam, Dar es Salaam, Tanzania.,Plant Protection Department, Swedish University of Agricultural Sciences (SLU), Alnarp, Sweden
| | - Donatha Tibuhwa
- Department of Molecular Biology and Biotechnology, University of Dar es Salaam, Dar es Salaam, Tanzania
| | - Sylvester Lyantagaye
- Department of Molecular Biology and Biotechnology, University of Dar es Salaam, Dar es Salaam, Tanzania
| | - Wai Yin Chan
- Biotechnology Platform (BTP), Agricultural Research Council, Onderstepoort Veterinary Research Campus, Old Soutpan Rd, Onderstepoort, South Africa.,Department of Biochemistry, Genetics and Microbiology, University of Pretoria, Pretoria, South Africa.,Forestry and Agricultural Biotechnology Institute (FABI), DST-NRF Centre of Excellence in Tree Health Biotechnology (CTHB), University of Pretoria, Pretoria, South Africa
| | - Stefanie Lutz
- Agroscope, Molecular Diagnostics, Genomics and Bioinformatics & SIB Swiss Institute of Bioinformatics, Müller-Thurgau-Str. 29, 8820, Wädenswil, Switzerland
| | - Christian H Ahrens
- Agroscope, Molecular Diagnostics, Genomics and Bioinformatics & SIB Swiss Institute of Bioinformatics, Müller-Thurgau-Str. 29, 8820, Wädenswil, Switzerland
| | | | - Rainer Borriss
- Institut für Biologie, Humboldt Universität Berlin, Berlin, Germany
| |
Collapse
|