Abstract

Gene duplication is an important evolutionary mechanism capable of providing new genetic material, which in some instances can help organisms adapt to various environmental conditions. Recent studies, for example, have indicated that highly similar duplicate genes (HSDs) are aiding adaptation to extreme conditions via gene dosage. However, for most eukaryotic genomes HSDs remain uncharacterized, partly because they can be hard to identify and categorize efficiently and effectively. Here, we collected and curated HSDs in nuclear genomes from various model animals, land plants and algae and indexed them in an online, open-access sequence repository called HSDatabase. Currently, this database contains 117 864 curated HSDs from 40 distinct genomes; it includes statistics on the total number of HSDs per genome as well as individual HSD copy numbers/lengths and provides sequence alignments of the duplicate gene copies. HSDatabase also allows users to download sequences of gene copies, access genome browsers, and link out to other databases, such as Pfam and Kyoto Encyclopedia of Genes and Genomes. What is more, a built-in Basic Local Alignment Search Tool option is available to conveniently explore potential homologous sequences of interest within and across species. HSDatabase has a user-friendly interface and provides easy access to the source data. It can be used on its own for comparative analyses of gene duplicates or in conjunction with HSDFinder, a newly developed bioinformatics tool for identifying, annotating, categorizing and visualizing HSDs.

Database URL: http://hsdfinder.com/database/

Introduction

Gene duplication is a near-ubiquitous phenomenon throughout the eukaryotic tree of life (1), one that can be advantageous or disadvantageous, depending on the circumstances. For example, under certain conditions, it can be detrimental for an organism to retain highly similar expressed genes (2). Thus, with notable exceptions, it is relatively rare for species to maintain duplicate genes encoding the same functions (3). Nevertheless, it is becoming more apparent that in some situations the generation and maintenance of highly similar duplicate genes (HSDs) is possible, particularly for genes encoding products that are in high demand, such as histones or ribosomal proteins (4). Indeed, there are many examples suggesting that genes involved in stress response, sensory functions, transport and/or metabolism are likely to be fixed as duplicated copies given specific environmental conditions (5).

Recently, Zhang et al. (6) revealed that hundreds of HSDs, involved in diverse cellular processes, are maintained in the psychrophilic Antarctic green alga Chlamydomonas sp. UWO241, which was recently renamed Chlamydomonas priscuii (7). It is believed that these HSDs are aiding its survival via gene dosage (8). Unfortunately, the HSDs from most other eukaryotic genomes, particularly those of algae, remain uncharacterized. This is partly because the experimental methods for identifying HSDs are time-consuming and labor-intensive. Many of the available bioinformatics tools for characterizing homologs are limited by their designs (e.g. they only identify orthologs) or their specificity (e.g. they only identify retrocopies or co-localized duplicates) (9–13). Consequently, we recently developed a web-based tool called HSDFinder that can identify HSDs in eukaryotic genomes with high accuracy and reliability (14). For example, HSDFinder predicted 336 and 265 HSDs in the psychrophilic green algae UWO241 and Chlamydomonas sp. ICE-L (6), respectively, which is consistent with other experimental data (8). By applying HSDFinder to a variety of other species (15), we predicted and cataloged thousands of HSD candidates, which are now curated and documented in a new online repository called HSDatabase. Currently, it houses 117 864 HSDs from 40 eukaryotic species, with a focus on green algae, animals and land plants.

Here, we briefly introduce the general features as well as the procedures and principles for collecting data from HSDatabase. In short, HSDatabase contains information on HSD number, gene copy number and gene copy length. Additionally, the protein functional domains and associated pathways of the HSDs can be retrieved from the Kyoto Encyclopedia of Genes and Genomes (KEGG) and InterProScan (16). A built-in Basic Local Alignment Search Tool (BLAST)-search option is also provided, allowing users to conveniently explore potential homologous sequences of interest within and across species. HSDatabase also provides data on a range of other parameters about gene duplicates, such as the number of HSD per Mb, the most commonly conserved domains among HSDs and the functional categories of HSDs. It is our hope to build a comparative analysis framework across species, especially for best-assembled eukaryotic genomes from species living in extreme environments, to better understand the role of gene duplication in adaptive evolution.

Materials and methods

Database collection

HSDs were identified in 40 well-assembled nuclear genomes from diverse model species, including land plants (e.g. Arabidopsis thaliana and Zea mays), algae (e.g. Chlamydomonas reinhardtii and Fragilariopsis cylindrus) and animals (e.g. Drosophila melanogaster, Homo sapiens and Mus musculus) (Figure 1). We focused on model animal and plant genomes because of their high-quality assemblies and annotations. The genome sequences of the selected species are all retrievable from the National Center for Biotechnology Information (NCBI) (17) (Table 1). The HSDs, which are represented by gene copies with nearly identical lengths and similar gene structures, were identified using HSDFinder (14). The identification method is based on all-against-all BLASTP analyses (18) carried out using uniform homology assessment metrics: E-value cut-off ≤1e−10, amino acid pairwise identity ≥90% and amino acid aligned length variance ≤10. Note, the short form of these parameters is denoted as ‘90%_10aa’. Additionally, putative HSDs were expected to have similar structural information, such as matching protein family (Pfam) domains (19), corresponding InterPro annotations (16) and/or nearly identical conserved residues. The InterProScan tool (16), which is an integrated platform for protein signatures, was used to collect the structural information of the HSDs. The all-against-all BLAST and InterProScan results (tab-delimited files) were fed into HSDFinder to generate HSD candidates in an 8-column tab-delimited file (Figure 2A). These candidates were identified by parsing the BLAST all-against-all protein similarity search results with the homology metrics: amino acid pairwise identity and amino acid aligned length variance. To collect and curate the data in HSDatabase, we performed a series of combo thresholds for filtering putatively functional gene copies (described below at Database curation section).

Taxonomic tree of 40 eukaryotic species in four highlighted categories. Stramenopila, Plantae , Fungi and Animalia are in blue, orange, green and red, respectively. The tree topologies were inferred by Taxonomy Common Tree from NCBI (https://www.ncbi.nlm.nih.gov/Taxonomy/CommonTree/wwwcmt.cgi).
Figure 1.

Taxonomic tree of 40 eukaryotic species in four highlighted categories. Stramenopila, Plantae , Fungi and Animalia are in blue, orange, green and red, respectively. The tree topologies were inferred by Taxonomy Common Tree from NCBI (https://www.ncbi.nlm.nih.gov/Taxonomy/CommonTree/wwwcmt.cgi).

Table 1.

Summary statistics of the curated HSD groups in the selected genomes from HSDatabase

Species name (common name)ClassificationCurated HSD groups #Genome size (Mb)No. of considered genesGene copiesHSDs/genesHSDs/Mb2-group HSDsa3-group HSDs≥4-group HSDsGenome assembly accession numberbRef
Ailuropoda melanoleuca (giant panda)Animalia25342371.822 45080910.1131.0681570440524GCF_002007445.1(23)
Bos taurus (cattle)Animalia22382667.8521 03673990.1060.8391433433372GCF_002263795.1(24)
Canis lupus familiaris (dog)Animalia23722344.0920 94573140.1131.0121541434397GCF_014441545.1(25)
Danio rerio (zebrafish)Animalia37991405.126 43911 9300.1442.7042263689847GCF_000002035.6(26)
Drosophila melanogaster (fruit fly)Animalia782138.9313 73922300.0575.629557118107GCF_000001215.4(27)
Equus caballus (horse)Animalia24032474.9321 12674990.1140.9711509468426GCF_002863925.1(28)
Felis catus (domestic cat)Animalia21672493.1420 45266910.1060.8691409396362GCF_018350175.1(29)
Gadus morhua (Atlantic cod)Animalia3269627.0423 48297570.1395.2131969633667GCF_902167405.1(30)
Gallus gallus (chicken)Animalia20261037.1717 79760920.1141.9531323374329GCF_016699485.2(31)
Gorilla gorilla (western gorilla)Animalia23353063.3620 63266830.1130.7621531434370GCF_008122165.1(32)
Homo sapiens (human)Animalia21782864.1419 53163520.1120.7601442395341GCF_000001405.39(33)
Hypsibius dujardini (waterbear)Animalia2290182.1520 85357100.11012.5721618366306GCA_002082055.1(34)
Loxodonta africana (African savanna elephant)Animalia20313196.7421 09473710.0960.6351327360344GCF_000001905.1(35)
Meleagris gallopavo (turkey)Animalia14861058.6417 97439680.0831.4041061240185GCF_000146605.3(36)
Mus musculus (house mouse)Animalia24022588.6230 73688550.1080.9281480459463GCF_000001635.27(37)
Rattus norvegicus (Norway rat)Animalia24342647.9222 21987570.1100.9191478471485GCF_015227675.2(38)
Saccharomyces cerevisiae (yeast)Fungi39711.83600210060.07333.5593314125GCA_003086655.1(39)
Arabidopsis lyrataPlantae5302202.9729 81716 9010.17826.12231049851213GCF_000004255.2(40)
Arabidopsis thaliana (thale cress)Plantae4428119.7527 56014 2250.16136.97726307931005GCF_000001735.4(41)
Brassica oleracea (wild cabbage)Plantae8918529.9244 38230 5110.20116.829456518412512GCF_000695525.1(42)
Carica papaya (papaya)Plantae2094360.6318 12663110.1165.8071299374421GCF_000150535.2(43)
Chlamydomonas eustigma (green alga)Plantae96666.6314 16121990.06814.49877810979GCA_002335675.1(20)
Chlamydomonas reinhardtii (green alga)Plantae1129111.1119 87031600.06410.161740187202GCF_000002595.2(44)
Chlamydomonas sp. ICE-LPlantae1540541.8617 73138530.0782.8421139224177GCA_013435795.1(21)
Chlamydomonas sp. UWO 241 (green alga)Plantae1112211.6416 01832820.0685.254741139232GCA_016618255.1(6)
Coccomyxa subellipsoidea C-169 (green alga)Plantae36048.83983910150.0377.3732814336GCA_000258705.1(45)
Cucumis sativus (cucumber)Plantae2891240.9920 03895580.14411.9961655532704GCF_000004075.3(46)
Dunaliella salina (green alga)Plantae1589343.718 74038590.0954.6231227194168GCA_002284615.2(47)
Fragilariopsis cylindrus (diatom)Stramenopila112974.7618 11131920.06215.102766172191GCA_001750085.1(48)
Glycine max (soybean)Plantae11 107995.2747 06438 2740.23611.160655912953253GCF_000004515.6(49)
Gonium pectorale (green alga)Plantae1028148.8116 29026690.0636.908719143166GCA_001584585.1(50)
Musa acuminata (dwarf banana)Plantae5489461.5422 17720 9340.17911.893263310831773GCF_000313855.2(51)
Oryza sativa (rice)Plantae4531386.4928 73514 7040.15811.72326418121078GCF_001433935.1(52)
Prunus persica (peach)Plantae3454220.923 13311 2100.14915.6362013611830GCF_000346465.2(53)
Solanum lycopersicum (tomato)Plantae4144809.1825 61213 7110.1625.12123457681031GCF_000188115.4(54)
Solanum tuberosum (potato)Plantae4733768.228 40715 9260.1676.16126478801206GCF_000226075.1(55)
Theobroma cacao (cacao)Plantae3074335.4421 51799330.1439.1641775554745GCF_000208745.1(56)
Vitis vinifera (wine grape)Plantae4039427.2125 83013 6130.1569.45422937221024GCF_000003745.3(57)
Volvox carteri (green alga)Plantae863137.6814 43626000.0616.268509152202GCA_000143455.1(58)
Zea mays (Maize)Plantae68012191.634 32822 4990.1983.103391011461745GCF_902167145.1(59)
Species name (common name)ClassificationCurated HSD groups #Genome size (Mb)No. of considered genesGene copiesHSDs/genesHSDs/Mb2-group HSDsa3-group HSDs≥4-group HSDsGenome assembly accession numberbRef
Ailuropoda melanoleuca (giant panda)Animalia25342371.822 45080910.1131.0681570440524GCF_002007445.1(23)
Bos taurus (cattle)Animalia22382667.8521 03673990.1060.8391433433372GCF_002263795.1(24)
Canis lupus familiaris (dog)Animalia23722344.0920 94573140.1131.0121541434397GCF_014441545.1(25)
Danio rerio (zebrafish)Animalia37991405.126 43911 9300.1442.7042263689847GCF_000002035.6(26)
Drosophila melanogaster (fruit fly)Animalia782138.9313 73922300.0575.629557118107GCF_000001215.4(27)
Equus caballus (horse)Animalia24032474.9321 12674990.1140.9711509468426GCF_002863925.1(28)
Felis catus (domestic cat)Animalia21672493.1420 45266910.1060.8691409396362GCF_018350175.1(29)
Gadus morhua (Atlantic cod)Animalia3269627.0423 48297570.1395.2131969633667GCF_902167405.1(30)
Gallus gallus (chicken)Animalia20261037.1717 79760920.1141.9531323374329GCF_016699485.2(31)
Gorilla gorilla (western gorilla)Animalia23353063.3620 63266830.1130.7621531434370GCF_008122165.1(32)
Homo sapiens (human)Animalia21782864.1419 53163520.1120.7601442395341GCF_000001405.39(33)
Hypsibius dujardini (waterbear)Animalia2290182.1520 85357100.11012.5721618366306GCA_002082055.1(34)
Loxodonta africana (African savanna elephant)Animalia20313196.7421 09473710.0960.6351327360344GCF_000001905.1(35)
Meleagris gallopavo (turkey)Animalia14861058.6417 97439680.0831.4041061240185GCF_000146605.3(36)
Mus musculus (house mouse)Animalia24022588.6230 73688550.1080.9281480459463GCF_000001635.27(37)
Rattus norvegicus (Norway rat)Animalia24342647.9222 21987570.1100.9191478471485GCF_015227675.2(38)
Saccharomyces cerevisiae (yeast)Fungi39711.83600210060.07333.5593314125GCA_003086655.1(39)
Arabidopsis lyrataPlantae5302202.9729 81716 9010.17826.12231049851213GCF_000004255.2(40)
Arabidopsis thaliana (thale cress)Plantae4428119.7527 56014 2250.16136.97726307931005GCF_000001735.4(41)
Brassica oleracea (wild cabbage)Plantae8918529.9244 38230 5110.20116.829456518412512GCF_000695525.1(42)
Carica papaya (papaya)Plantae2094360.6318 12663110.1165.8071299374421GCF_000150535.2(43)
Chlamydomonas eustigma (green alga)Plantae96666.6314 16121990.06814.49877810979GCA_002335675.1(20)
Chlamydomonas reinhardtii (green alga)Plantae1129111.1119 87031600.06410.161740187202GCF_000002595.2(44)
Chlamydomonas sp. ICE-LPlantae1540541.8617 73138530.0782.8421139224177GCA_013435795.1(21)
Chlamydomonas sp. UWO 241 (green alga)Plantae1112211.6416 01832820.0685.254741139232GCA_016618255.1(6)
Coccomyxa subellipsoidea C-169 (green alga)Plantae36048.83983910150.0377.3732814336GCA_000258705.1(45)
Cucumis sativus (cucumber)Plantae2891240.9920 03895580.14411.9961655532704GCF_000004075.3(46)
Dunaliella salina (green alga)Plantae1589343.718 74038590.0954.6231227194168GCA_002284615.2(47)
Fragilariopsis cylindrus (diatom)Stramenopila112974.7618 11131920.06215.102766172191GCA_001750085.1(48)
Glycine max (soybean)Plantae11 107995.2747 06438 2740.23611.160655912953253GCF_000004515.6(49)
Gonium pectorale (green alga)Plantae1028148.8116 29026690.0636.908719143166GCA_001584585.1(50)
Musa acuminata (dwarf banana)Plantae5489461.5422 17720 9340.17911.893263310831773GCF_000313855.2(51)
Oryza sativa (rice)Plantae4531386.4928 73514 7040.15811.72326418121078GCF_001433935.1(52)
Prunus persica (peach)Plantae3454220.923 13311 2100.14915.6362013611830GCF_000346465.2(53)
Solanum lycopersicum (tomato)Plantae4144809.1825 61213 7110.1625.12123457681031GCF_000188115.4(54)
Solanum tuberosum (potato)Plantae4733768.228 40715 9260.1676.16126478801206GCF_000226075.1(55)
Theobroma cacao (cacao)Plantae3074335.4421 51799330.1439.1641775554745GCF_000208745.1(56)
Vitis vinifera (wine grape)Plantae4039427.2125 83013 6130.1569.45422937221024GCF_000003745.3(57)
Volvox carteri (green alga)Plantae863137.6814 43626000.0616.268509152202GCA_000143455.1(58)
Zea mays (Maize)Plantae68012191.634 32822 4990.1983.103391011461745GCF_902167145.1(59)
a

2-group HSDs refers to the number of curated HSD groups with only two gene copies.

b

Accession numbers are from the US NCBI GenBank assembly accession.

Table 1.

Summary statistics of the curated HSD groups in the selected genomes from HSDatabase

Species name (common name)ClassificationCurated HSD groups #Genome size (Mb)No. of considered genesGene copiesHSDs/genesHSDs/Mb2-group HSDsa3-group HSDs≥4-group HSDsGenome assembly accession numberbRef
Ailuropoda melanoleuca (giant panda)Animalia25342371.822 45080910.1131.0681570440524GCF_002007445.1(23)
Bos taurus (cattle)Animalia22382667.8521 03673990.1060.8391433433372GCF_002263795.1(24)
Canis lupus familiaris (dog)Animalia23722344.0920 94573140.1131.0121541434397GCF_014441545.1(25)
Danio rerio (zebrafish)Animalia37991405.126 43911 9300.1442.7042263689847GCF_000002035.6(26)
Drosophila melanogaster (fruit fly)Animalia782138.9313 73922300.0575.629557118107GCF_000001215.4(27)
Equus caballus (horse)Animalia24032474.9321 12674990.1140.9711509468426GCF_002863925.1(28)
Felis catus (domestic cat)Animalia21672493.1420 45266910.1060.8691409396362GCF_018350175.1(29)
Gadus morhua (Atlantic cod)Animalia3269627.0423 48297570.1395.2131969633667GCF_902167405.1(30)
Gallus gallus (chicken)Animalia20261037.1717 79760920.1141.9531323374329GCF_016699485.2(31)
Gorilla gorilla (western gorilla)Animalia23353063.3620 63266830.1130.7621531434370GCF_008122165.1(32)
Homo sapiens (human)Animalia21782864.1419 53163520.1120.7601442395341GCF_000001405.39(33)
Hypsibius dujardini (waterbear)Animalia2290182.1520 85357100.11012.5721618366306GCA_002082055.1(34)
Loxodonta africana (African savanna elephant)Animalia20313196.7421 09473710.0960.6351327360344GCF_000001905.1(35)
Meleagris gallopavo (turkey)Animalia14861058.6417 97439680.0831.4041061240185GCF_000146605.3(36)
Mus musculus (house mouse)Animalia24022588.6230 73688550.1080.9281480459463GCF_000001635.27(37)
Rattus norvegicus (Norway rat)Animalia24342647.9222 21987570.1100.9191478471485GCF_015227675.2(38)
Saccharomyces cerevisiae (yeast)Fungi39711.83600210060.07333.5593314125GCA_003086655.1(39)
Arabidopsis lyrataPlantae5302202.9729 81716 9010.17826.12231049851213GCF_000004255.2(40)
Arabidopsis thaliana (thale cress)Plantae4428119.7527 56014 2250.16136.97726307931005GCF_000001735.4(41)
Brassica oleracea (wild cabbage)Plantae8918529.9244 38230 5110.20116.829456518412512GCF_000695525.1(42)
Carica papaya (papaya)Plantae2094360.6318 12663110.1165.8071299374421GCF_000150535.2(43)
Chlamydomonas eustigma (green alga)Plantae96666.6314 16121990.06814.49877810979GCA_002335675.1(20)
Chlamydomonas reinhardtii (green alga)Plantae1129111.1119 87031600.06410.161740187202GCF_000002595.2(44)
Chlamydomonas sp. ICE-LPlantae1540541.8617 73138530.0782.8421139224177GCA_013435795.1(21)
Chlamydomonas sp. UWO 241 (green alga)Plantae1112211.6416 01832820.0685.254741139232GCA_016618255.1(6)
Coccomyxa subellipsoidea C-169 (green alga)Plantae36048.83983910150.0377.3732814336GCA_000258705.1(45)
Cucumis sativus (cucumber)Plantae2891240.9920 03895580.14411.9961655532704GCF_000004075.3(46)
Dunaliella salina (green alga)Plantae1589343.718 74038590.0954.6231227194168GCA_002284615.2(47)
Fragilariopsis cylindrus (diatom)Stramenopila112974.7618 11131920.06215.102766172191GCA_001750085.1(48)
Glycine max (soybean)Plantae11 107995.2747 06438 2740.23611.160655912953253GCF_000004515.6(49)
Gonium pectorale (green alga)Plantae1028148.8116 29026690.0636.908719143166GCA_001584585.1(50)
Musa acuminata (dwarf banana)Plantae5489461.5422 17720 9340.17911.893263310831773GCF_000313855.2(51)
Oryza sativa (rice)Plantae4531386.4928 73514 7040.15811.72326418121078GCF_001433935.1(52)
Prunus persica (peach)Plantae3454220.923 13311 2100.14915.6362013611830GCF_000346465.2(53)
Solanum lycopersicum (tomato)Plantae4144809.1825 61213 7110.1625.12123457681031GCF_000188115.4(54)
Solanum tuberosum (potato)Plantae4733768.228 40715 9260.1676.16126478801206GCF_000226075.1(55)
Theobroma cacao (cacao)Plantae3074335.4421 51799330.1439.1641775554745GCF_000208745.1(56)
Vitis vinifera (wine grape)Plantae4039427.2125 83013 6130.1569.45422937221024GCF_000003745.3(57)
Volvox carteri (green alga)Plantae863137.6814 43626000.0616.268509152202GCA_000143455.1(58)
Zea mays (Maize)Plantae68012191.634 32822 4990.1983.103391011461745GCF_902167145.1(59)
Species name (common name)ClassificationCurated HSD groups #Genome size (Mb)No. of considered genesGene copiesHSDs/genesHSDs/Mb2-group HSDsa3-group HSDs≥4-group HSDsGenome assembly accession numberbRef
Ailuropoda melanoleuca (giant panda)Animalia25342371.822 45080910.1131.0681570440524GCF_002007445.1(23)
Bos taurus (cattle)Animalia22382667.8521 03673990.1060.8391433433372GCF_002263795.1(24)
Canis lupus familiaris (dog)Animalia23722344.0920 94573140.1131.0121541434397GCF_014441545.1(25)
Danio rerio (zebrafish)Animalia37991405.126 43911 9300.1442.7042263689847GCF_000002035.6(26)
Drosophila melanogaster (fruit fly)Animalia782138.9313 73922300.0575.629557118107GCF_000001215.4(27)
Equus caballus (horse)Animalia24032474.9321 12674990.1140.9711509468426GCF_002863925.1(28)
Felis catus (domestic cat)Animalia21672493.1420 45266910.1060.8691409396362GCF_018350175.1(29)
Gadus morhua (Atlantic cod)Animalia3269627.0423 48297570.1395.2131969633667GCF_902167405.1(30)
Gallus gallus (chicken)Animalia20261037.1717 79760920.1141.9531323374329GCF_016699485.2(31)
Gorilla gorilla (western gorilla)Animalia23353063.3620 63266830.1130.7621531434370GCF_008122165.1(32)
Homo sapiens (human)Animalia21782864.1419 53163520.1120.7601442395341GCF_000001405.39(33)
Hypsibius dujardini (waterbear)Animalia2290182.1520 85357100.11012.5721618366306GCA_002082055.1(34)
Loxodonta africana (African savanna elephant)Animalia20313196.7421 09473710.0960.6351327360344GCF_000001905.1(35)
Meleagris gallopavo (turkey)Animalia14861058.6417 97439680.0831.4041061240185GCF_000146605.3(36)
Mus musculus (house mouse)Animalia24022588.6230 73688550.1080.9281480459463GCF_000001635.27(37)
Rattus norvegicus (Norway rat)Animalia24342647.9222 21987570.1100.9191478471485GCF_015227675.2(38)
Saccharomyces cerevisiae (yeast)Fungi39711.83600210060.07333.5593314125GCA_003086655.1(39)
Arabidopsis lyrataPlantae5302202.9729 81716 9010.17826.12231049851213GCF_000004255.2(40)
Arabidopsis thaliana (thale cress)Plantae4428119.7527 56014 2250.16136.97726307931005GCF_000001735.4(41)
Brassica oleracea (wild cabbage)Plantae8918529.9244 38230 5110.20116.829456518412512GCF_000695525.1(42)
Carica papaya (papaya)Plantae2094360.6318 12663110.1165.8071299374421GCF_000150535.2(43)
Chlamydomonas eustigma (green alga)Plantae96666.6314 16121990.06814.49877810979GCA_002335675.1(20)
Chlamydomonas reinhardtii (green alga)Plantae1129111.1119 87031600.06410.161740187202GCF_000002595.2(44)
Chlamydomonas sp. ICE-LPlantae1540541.8617 73138530.0782.8421139224177GCA_013435795.1(21)
Chlamydomonas sp. UWO 241 (green alga)Plantae1112211.6416 01832820.0685.254741139232GCA_016618255.1(6)
Coccomyxa subellipsoidea C-169 (green alga)Plantae36048.83983910150.0377.3732814336GCA_000258705.1(45)
Cucumis sativus (cucumber)Plantae2891240.9920 03895580.14411.9961655532704GCF_000004075.3(46)
Dunaliella salina (green alga)Plantae1589343.718 74038590.0954.6231227194168GCA_002284615.2(47)
Fragilariopsis cylindrus (diatom)Stramenopila112974.7618 11131920.06215.102766172191GCA_001750085.1(48)
Glycine max (soybean)Plantae11 107995.2747 06438 2740.23611.160655912953253GCF_000004515.6(49)
Gonium pectorale (green alga)Plantae1028148.8116 29026690.0636.908719143166GCA_001584585.1(50)
Musa acuminata (dwarf banana)Plantae5489461.5422 17720 9340.17911.893263310831773GCF_000313855.2(51)
Oryza sativa (rice)Plantae4531386.4928 73514 7040.15811.72326418121078GCF_001433935.1(52)
Prunus persica (peach)Plantae3454220.923 13311 2100.14915.6362013611830GCF_000346465.2(53)
Solanum lycopersicum (tomato)Plantae4144809.1825 61213 7110.1625.12123457681031GCF_000188115.4(54)
Solanum tuberosum (potato)Plantae4733768.228 40715 9260.1676.16126478801206GCF_000226075.1(55)
Theobroma cacao (cacao)Plantae3074335.4421 51799330.1439.1641775554745GCF_000208745.1(56)
Vitis vinifera (wine grape)Plantae4039427.2125 83013 6130.1569.45422937221024GCF_000003745.3(57)
Volvox carteri (green alga)Plantae863137.6814 43626000.0616.268509152202GCA_000143455.1(58)
Zea mays (Maize)Plantae68012191.634 32822 4990.1983.103391011461745GCF_902167145.1(59)
a

2-group HSDs refers to the number of curated HSD groups with only two gene copies.

b

Accession numbers are from the US NCBI GenBank assembly accession.

The workflow of HSDatabase. (A) Steps for using HSDFinder to collect candidate HSDs. (B) Manual curation of HSDs via filtering and adding new HSD candidates prior to being deposited into HSDatabase. (C) Steps of accessing HSD data in HSDatabase, including browsing via organism name, blasting query sequences against the database and searching through the HSD and gene copy IDs.
Figure 2.

The workflow of HSDatabase. (A) Steps for using HSDFinder to collect candidate HSDs. (B) Manual curation of HSDs via filtering and adding new HSD candidates prior to being deposited into HSDatabase. (C) Steps of accessing HSD data in HSDatabase, including browsing via organism name, blasting query sequences against the database and searching through the HSD and gene copy IDs.

Database curation

Prior to uploading data into HSDatabase, we curated HSD candidates by filtering for redundancy and adding the newly curated HSDs (Figure 2B). For genes that have alternative protein products, we selected the longest gene isoform to reduce redundancy. Since highly similar gene copies are grouped together as HSDs based on a simple transitive link between the remaining genes (14), it is possible for some highly duplicated genes to form mega HSD groups with varied gene copy lengths, especially those encoding histones, ribosomal proteins or retro-transcriptases. Moreover, some gene copies might appear multiple times causing redundancy among different HSD groups, which is because the BLAST algorithm limits the maximum target hits by default. In these cases, we manually curated the HSD groups, minimizing redundant gene copies.

Since the similarity of duplicate genes within and among genomes can vary significantly, we added newly curated HSDs to the database using a combination of thresholds to acquire a larger dataset of HSD candidates. We added the HSD candidates one after another at different homology assessment metrics (i.e. HSDs identified at more relaxed thresholds were treated more strictly than those found using more conservative thresholds) (Figure 2B). For example, HSDs identified at a threshold of 90%_30aa were added on to those identified at a threshold of 90%_10aa (denoted as ‘90%_30aa+90%_10aa’); any redundant HSDs candidates picked out at this combo threshold were removed if the more relaxed threshold (i.e. 90%_30aa) had the identical genes or contained the same gene copies from the stricter cut-off (i.e. 90%_10aa). Moreover, any HSD candidates pinpointed at the combo threshold (90%_30aa + 90%_10aa) were removed if the minimum gene copy length was less than half of the maximum gene copy length for each HSD or if HSD candidates had gene copies with incomplete conserved domains (i.e. a different number of Pfam domains). After filtering the combo threshold at (90%_30aa + 90%_10aa), we added on a more relaxed threshold 90%_50aa (i.e. 90%_50aa + (90%_30aa + 90%_10aa)) and then carried out the same HSD candidate removal/filtering process. To minimize redundancy and to acquire a larger dataset of HSD candidates, we processed each selected species with the following combination of thresholds: |${\rm{E + }}\left( {{\rm{D + }}\left( {{\rm{C + }}\left( {{\rm{B + A}}} \right)} \right)} \right)$|⁠.
 
 
 
 

Database implementation

The database was built with the Django 3.0.5 web framework (https://www.djangoproject.com/), and all data were stored in an SQLite 3.36.0 database (https://www.sqlite.org/index.html) on an Amazon web server. Webpage templates used Bootstrap framework (https://getbootstrap.com/), D3.js (https://d3js.org), jQuery (http://jquery.com) and Bootstrap Table (https://bootstrap-table.com/) libraries to establish a user-friendly, front-end interface. On the browse page, NCBI’s Sequence Viewer 3.44.0 (https://www.ncbi.nlm.nih.gov/projects/sviewer/) was employed to build a fast and scalable genome browser.

Results and discussion

Database content and analysis

HSDatabase was built using a relational database (MySQL) allowing the rapid retrieval of data and making resources easily maintainable. One entry corresponds to one eukaryote genome. The genomes can be accessed via the organism table or the taxonomic tree. The genome entry is then split into various subcategories of HSD entries. Database access is via a web interface based on python script and provides various ways to search for HSD entries, including species name, unique HSD IDs and gene copy IDs.

Using HSDFinder (15), we collected and curated 117 864 HSDs (representing 379 844 gene copies) from 40 well-assembled nuclear genomes of diverse model species (Table 1). Various green algae were included because of our specific interest in algal genomics and also because of their relatively modest genome sizes and penchant for gene duplications. For example, the acidophilic green alga Chlamydomonas eustigma is known to have large numbers of gene duplicates in its nuclear genome, including 10 gene copies for arsenate reductase and 20 for glutaredoxin (20). Similarly, the psychrophilic green alga Chlamydomonas sp. ICE-L contains multiple copies of genes encoding carotene biosynthesis-related protein and Lhc-like protein (Lhl4) (21). These data are consistent with our identification of large numbers of HSDs in C. eustigma (276) and ICE-L (265) (Table 1), suggesting a potential adaptative role of gene duplication under different extreme environmental conditions.

Compared to algae, the investigated land plants had higher detected numbers of HSDs as well as larger ratios of HSDs/Mb and HSDs/genes (Table 1). For example, the HSDs/Mb values for Arabidopsis lyrata and A. thaliana are 26 and 37, respectively, whereas the average HSDs/Mb value among selected green algae is 8.2. Compared to algae and land plants, the HSDs/Mb values in animals are generally quite low with the exception of Hypsibius dujardini (13.6) and D. melanogaster (5.6). Two-group HSDs (i.e. HSDs containing two gene copies) represent the majority (>50%) of total HSDs for all explored species.

As for the associated functions of the detected HSDs, three green algal species with relatively large values of HSDs/genes were compared previously. These algae can survive various extreme environmental conditions and include the Antarctic psychrophilic green algae UWO241 (0.068) and ICE-L (0.078) and the acidophilic C. eustigma (0.068) (Table 1). The identified duplicates are involved in a diversity of cellular pathways, including gene expression, cell growth, membrane transport and energy metabolism, but also include ribosomal proteins (6, 14). Although HSDs for protein translation, DNA packaging and photosynthesis are particularly prevalent, around 30% of the HSDs are hypothetical proteins without any Pfam domains.

Database composition and usage

Information about specific HSDs and their associated gene copies for a given species can be obtained through the ‘Browse’ and ‘Search’ tabs, which are located on the menu bar at the top of the page, or using nucleotide/amino acid sequences as queries to search against the database via BLAST (i.e. BLASTP or BLASTX). To categorize duplicated genes into their functional categories, KEGG pathway schematics are available for each species.

Browse

By selecting the ‘Browse’ option from the main menu, users are offered three ways to explore their species of interest. First, they can simply click the organism name on the taxonomic tree containing the 40 species. Secondly, users can select the ‘Plantae and Stramenopila’ or ‘Animalia and Fungi’ tabs (Figure 3A), which contain 23 and 17 species, respectively. Selecting a tab takes users to a summary table that contains the organism names, number of HSDs, species background information, GenBank accession links to genome assemblies, and reference links to PubMed.

Screenshots of the HSDatabase interface. There are four main functions in the menu page: (A) Browse the database via species entries; (B) search the database via the HSDatabase unique ID (e.g. hsd_id_Athaliana_1) or gene ID (e.g. NP_200993.1); (C) use BLAST to search the database via amino acid sequence in FASTA format; and (D) categorize the gene copies and HSDs under the KEGG pathway functional categories.
Figure 3.

Screenshots of the HSDatabase interface. There are four main functions in the menu page: (A) Browse the database via species entries; (B) search the database via the HSDatabase unique ID (e.g. hsd_id_Athaliana_1) or gene ID (e.g. NP_200993.1); (C) use BLAST to search the database via amino acid sequence in FASTA format; and (D) categorize the gene copies and HSDs under the KEGG pathway functional categories.

Selecting a specific species from the browse page leads to the respective HSDs summary page (Figure 4A), which gives data on the total number of HSDs, unique HSD IDs, gene copy GenBank IDs and number of gene copies; it also provides access to the data download function. Choosing one of the HSD ID entries, for example, brings up a page containing information and features of a detected HSD, including the associated gene copies for a unique HSD as well as the GenBank link, the sequence length, the Pfam domain ID/description and the InterPro database ID/description (Figure 4B). Clicking on the ‘genome browser’ tab allows for the visualization of a specific gene copy through the built-in NCBI genome browser (Figure 4C). The ‘FASTA sequence’ tab provides the option to download the sequence data (Figure 4D) and the ‘alignments and identity%’ tab brings up the gene copy alignment and percentage identity matrix created by the built-in Clustal v2.1 tool (22) (Figure 4E).

Summary of database information for a selected species. (A) HSDs collected in a table for a specific species. (B) Basic information of the unique HSD ID, gene copy ID and the associated links to Pfam domains and InterPro databases. (C) Linking gene copies to the genome browser. (D) The FASTA sequence downloads of gene copies. (E) Alignments and percentage sequence identities of gene copies.
Figure 4.

Summary of database information for a selected species. (A) HSDs collected in a table for a specific species. (B) Basic information of the unique HSD ID, gene copy ID and the associated links to Pfam domains and InterPro databases. (C) Linking gene copies to the genome browser. (D) The FASTA sequence downloads of gene copies. (E) Alignments and percentage sequence identities of gene copies.

Search

Through the search option from the main menu users can search unique HSD IDs or gene copy IDs against the database (Figure 3B); they can also set the selection categories to limit the search results, which can improve search efficiency. After activation of the search button, 30 results per page are displayed (Figure 3B) in a four-column table, including HSD name, gene copy name, number of gene copies and the external download link to the output data (tab-delimited file). Users can navigate through the results page or download specific HSD entries. As described in the Browse section, the data file includes various summary statistics on the HSDs (Figure 4B).

BLAST

The BLAST tool bar allows users to input a nucleotide or amino acid sequence (in FASTA format) and carry out a sequence similarity search using BLASTX or BLASTP. Users can specify the species against which the BLAST search will be performed. The E-value and maximum target sequence of results can also be adjusted, but all other parameters remain as default and cannot be changed (Figure 3C). The BLAST search output result is in the standard 13-column tabular format, including the linkable query sequence ID and HSD ID, percentage identity, aligned length and all other BLAST tabular output values. The most similar sequences are arranged at the top.

KEGG

The KEGG page contains details on the associated KEGG pathways of the HSD gene copies for the 40 species. To browse the data for a particular species, users can simply select the organism’s name. The 6-column table lists the gene copies and HSDs under KEGG functional categories (Figure 3D). Gene copies involved in the same KEGG pathway are detailed with the first KEGG category (e.g. Carbohydrate metabolism), then the secondary category (e.g. Glycolysis/Gluconeogenesis) and finally the KEGG pathway function description (e.g. ENO, eno; enolase). The KEGG ID (e.g. K01689) is linked to the external KEGG database, providing more detailed information.

Future direction and limitation

Now that HSDatabase is publicly available, the next step is to analyze duplicate genes across a broader range of species, which we plan on doing in the near future. Currently, the database includes a range of statistics (e.g. number of HSD per Mb), but we hope to add additional data in the coming years, including information on differential expression levels among duplicates, for instance, as well as data on rates of synonymous and nonsynonymous substitutions (dN/dS rates). The biggest challenge moving forward will be determining an appropriate threshold for accurately predicting HSDs. As research on gene duplicates improves, we may need to adjust the metrics (e.g. amino acid pairwise identity and amino acid aligned length variance) to find as many bona fide HSDs as possible.

Presently, there is no standard golden cut-off for identifying HSDs and there might never be one as there a multitude of forces, including lineage/genomic specific ones, that can impact the accuracy of the identification metrics. This is why users can employ different parameters in the HSDFinder tool (from 30 to 100% amino acid pairwise identity and from within 0 to 100 amino acid aligned length variance). In our case, we used a series of combination thresholds to curate the HSDs in HSDatabase. But due to the limitations of this strategy, there are some large groups of HSD candidates in the database that likely diverged in function from one another and, thus, are not inducing a gene dosage benefit. In the database, we have labeled these putatively diverged HSD groups as ‘candidate HSDs’ and have added a warning note that users should proceed with caution when working with these datasets. In the future, our goal is to guide users to species-specific thresholds and deposit more diverse eukaryotic species into the database.

Conclusions

With the decreasing cost of next-generation sequencing, biologists are dealing with ever larger amounts of data. However, many bioinformatics software suites require considerable knowledge of computer scripting and microprogramming. To facilitate the understanding and analysis of gene duplication in nuclear genomes, we developed HSDatabase, which currently contains 117 864 HSDs from 40 well-assembled eukaryotic genomes. In conjunction with HSDatabase, we designed HSDFinder, which can efficiently identify duplicated genes from unannotated genome sequences by integrating the results from InterProScan and KEGG. HSDatabase aims to become a useful platform for the identification and comprehensive analysis of HSDs in eukaryotic genomes, which could aid research into the mechanisms driving genome adaptation. In the future, the database will be updated by incorporating advancements in the field of gene duplication.

Supplementary data

Supplementary data are available at Database Online.

Acknowledgements

We want to especially thank the editors and reviewers for their professional comments that greatly improved this manuscript.

Funding

Discovery Grants from the Natural Sciences and Engineering Research Council of Canada.

Conflict of interest

None declared.

Data availability

The datasets of eukaryotes supporting the conclusions of this article are available from Joint Genome Institute (JGI) (https://phytozome.jgi.doe.gov/pz/portal.html) or National Center for Biotechnology Information (https://www.ncbi.nlm.nih.gov) database.

Author contributions

The study was conceptualized by X.Z. and D.R.S. The data were analyzed by X.Z., and Y.N.H. implemented the HSDatabase website. X.Z. and D.R.S. drafted the manuscript, and all authors commented to produce the manuscript for peer review.

References

1.

Ohno
S.
(
1970
)
Evolution by Gene Duplication
.
Springer
,
Berlin/Heidelberg, Germany
.

2.

Conrad
B.
and
Antonarakis
S.E.
(
2007
)
Gene duplication: a drive for phenotypic diversity and cause of human disease
.
Annu. Rev. Genomics Hum. Genet.
,
8
,
17
35
.

3.

Kubiak
M.R.
and
Makałowska
I.
(
2017
)
Protein-coding genes’ retrocopies and their functions
.
Viruses
,
9
,
1
27
.

4.

Zhang
J.
(
2003
)
Evolution by gene duplication: an update
.
Trends Ecol. Evol. (Amst.)
,
18
,
292
298
.

5.

Kondrashov
F.A.
(
2012
)
Gene duplication as a mechanism of genomic adaptation to a changing environment
.
Proc. Royal Soc. B
,
279
,
5048
5057
.

6.

Zhang
X.
,
Cvetkovska
M.
,
Morgan-Kiss
R.
 et al.  (
2021
)
Draft genome sequence of the Antarctic green alga Chlamydomonas sp. UWO241
.
iScience
,
24
,
1
9
.

7.

Stahl-Rommel
S.
,
Kalra
I.
,
D’Silva
S.
 et al.  (
2022
)
Cyclic electron flow (CEF) and ascorbate pathway activity provide constitutive photoprotection for the photopsychrophile, Chlamydomonas sp. UWO 241 (renamed Chlamydomonas priscuii)
.
Photosyn. Res.
,
151
,
235
250
.

8.

Cvetkovska
M.
,
Szyszka-Mroz
B.
,
Possmayer
M.
 et al.  (
2018
)
Characterization of photosynthetic ferredoxin from the Antarctic alga Chlamydomonas sp. UWO241 reveals novel features of cold adaptation
.
New Phytol.
,
219
,
588
604
.

9.

Rosikiewicz
W.
,
Kabza
M.
,
Kosiński
J.G.
 et al.  (
2017
)
RetrogeneDB—a database of plant and animal retrocopies
.
Database
,
2017
,
1
11
.

10.

Kabza
M.
,
Ciomborowska
J.
and
Makałowska
I.
(
2014
)
RetrogeneDB—a database of animal retrogenes
.
Mol. Biol. Evol.
,
31
,
1646
1648
.

11.

Ouedraogo
M.
,
Bettembourg
C.
,
Bretaudeau
A.
 et al.  (
2012
)
The duplicated genes database: identification and functional annotation of co-localised duplicated genes across genomes
.
PLoS One
,
7
,
1
8
.

12.

Li
L.
,
Stoeckert
C.J.
and
Roos
D.S.
(
2003
)
OrthoMCL: identification of ortholog groups for eukaryotic genomes
.
Genome Res.
,
13
,
2178
2189
.

13.

Zdobnov
E.M.
,
Tegenfeldt
F.
,
Kuznetsov
D.
 et al.  (
2017
)
OrthoDB v9. 1: cataloging evolutionary and functional annotations for animal, fungal, plant, archaeal, bacterial and viral orthologs
.
Nucleic Acids Res.
,
45
,
D744
D749
.

14.

Zhang
X.
,
Hu
Y.
and
Smith
D.R.
(
2021
)
HSDFinder: a BLAST-based strategy for identifying highly similar duplicated genes in eukaryotic genomes
.
Front. Bioinf.
,
1
,
1
12
.

15.

Zhang
X.
,
Hu
Y.
and
Smith
D.R.
(
2021
)
Protocol for HSDFinder: identifying, annotating, categorizing, and visualizing duplicated genes in eukaryotic genomes
.
STAR Protoc.
,
2
,
1
18
.

16.

Quevillon
E.
,
Silventoinen
V.
,
Pillai
S.
 et al.  (
2005
)
InterProScan: protein domains identifier
.
Nucleic Acids Res.
,
33
,
116
120
.

17.

Pruitt
K.D.
,
Tatusova
T.
and
Maglott
D.R.
(
2005
)
NCBI reference sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins
.
Nucleic Acids Res.
,
33
,
D501
D504
.

18.

Altschul
S.F.
 et al.  (
1997
)
Gapped BLAST and PSI-BLAST: a new generation of protein database search programs
.
Nucleic Acids Res.
,
25
,
3389
3402
.

19.

Finn
R.D.
,
Bateman
A.
,
Clements
J.
 et al.  (
2014
)
Pfam: the protein families database
.
Nucleic Acids Res.
,
42
,
222
230
.

20.

Hirooka
S.
,
Hirose
Y.
,
Kanesaki
Y.
 et al.  (
2017
)
Acidophilic green algal genome provides insights into adaptation to an acidic environment
.
Proc. Natl. Acad. Sci.
,
114
,
8304
8313
.

21.

Zhang
Z.
,
Qu
C.
,
Zhang
K.
 et al.  (
2020
)
Adaptation to extreme Antarctic environments revealed by the genome of a sea ice green alga
.
Curr. Biol.
,
30
,
1
12
.

22.

Sievers
F.
,
Wilm
A.
,
Dineen
D.
 et al.  (
2011
)
Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega
.
Mol. Syst. Biol.
,
7
,
1
6
.

23.

Li
R.
,
Fan
W.
,
Tian
G.
 et al.  (
2010
)
The sequence and de novo assembly of the giant panda genome
.
Nature
,
463
,
311
317
.

24.

Zimin
A.V.
,
Delcher
A.L.
,
Florea
L.
 et al.  (
2009
)
A whole-genome assembly of the domestic cow, Bos taurus
.
Genome Biol.
,
10
,
1
10
.

25.

Lindblad-Toh
K.
,
Wade
C.M.
,
Mikkelsen
T.S.
 et al.  (
2005
)
Genome sequence, comparative analysis and haplotype structure of the domestic dog
.
Nature
,
438
,
803
819
.

26.

Howe
K.
,
Clark
M.D.
,
Torroja
C.F.
 et al.  (
2013
)
The zebrafish reference genome sequence and its relationship to the human genome
.
Nature
,
496
,
498
503
.

27.

Hoskins
R.A.
,
Carlson
J.W.
,
Wan
K.H.
 et al.  (
2015
)
The Release 6 reference sequence of the Drosophila melanogaster genome
.
Genome Res.
,
25
,
445
458
.

28.

Wade
C.
,
Giulotto
E.
,
Sigurdsson
S.
 et al.  (
2009
)
Genome sequence, comparative analysis, and population genetics of the domestic horse
.
Science
,
326
,
865
867
.

29.

Lopez
J.V.
,
Cevario
S.
and
O’Brien
S.J.
(
1996
)
Complete nucleotide sequences of the domestic cat (Felis catus) mitochondrial genome and a transposed mtDNA tandem repeat (Numt) in the nuclear genome
.
Genomics
,
33
,
229
246
.

30.

Star
B.
,
Nederbragt
A.J.
,
Jentoft
S.
 et al.  (
2011
)
The genome sequence of Atlantic cod reveals a unique immune system
.
Nature
,
477
,
207
210
.

31.

Viertlboeck
B.C.
,
Habermann
F.A.
,
Schmitt
R.
 et al.  (
2005
)
The chicken leukocyte receptor complex: a highly diverse multigene family encoding at least six structurally distinct receptor types
.
J. Immunol.
,
175
,
385
393
.

32.

Hughes
J.F.
,
Skaletsky
H.
,
Pyntikova
T.
 et al.  (
2005
)
Conservation of Y-linked genes during human evolution revealed by comparative sequencing in chimpanzee
.
Nature
,
437
,
100
103
.

33.

Lander
E.S.
 et al.  (
2001
)
Initial sequencing and analysis of the human genome
.
Nature
,
409
,
860
921
.

34.

Koutsovoulos
G.
,
Kumar
S.
,
Laetsch
D.R.
 et al.  (
2016
)
No evidence for extensive horizontal gene transfer in the genome of the tardigrade Hypsibius dujardini
.
Proc. Natl. Acad. Sci.
,
113
,
5053
5058
.

35.

Guo
Y.
,
Bao
Y.
,
Wang
H.
 et al.  (
2011
)
A preliminary analysis of the immunoglobulin genes in the African elephant (Loxodonta africana)
.
PLoS One
,
6
,
1
14
.

36.

Dalloul
R.A.
,
Long
J.A.
,
Zimin
A.V.
 et al.  (
2010
)
Multi-platform next-generation sequencing of the domestic turkey (Meleagris gallopavo): genome assembly and analysis
.
PLoS Biol.
,
8
,
1
21
.

37.

Church
D.M.
,
Schneider
V.A.
,
Graves
T.
 et al.  (
2011
)
Modernizing reference genome assemblies
.
PLoS Biol.
,
9
,
1
5
.

38.

Gibbs
R.A.
and
Pachter
L.
(
2004
)
Genome sequence of the Brown Norway rat yields insights into mammalian evolution
.
Nature
,
428
,
493
521
.

39.

Shao
Y.
,
Lu
N.
,
Wu
Z.
 et al.  (
2018
)
Creating a functional single-chromosome yeast
.
Nature
,
560
,
331
335
.

40.

Hu
T.T.
,
Pattyn
P.
,
Bakker
E.G.
 et al.  (
2011
)
The Arabidopsis lyrata genome sequence and the basis of rapid genome size change
.
Nat. Genet.
,
43
,
476
481
.

41.

Sloan
D.B.
,
Wu
Z.
and
Sharbrough
J.
(
2018
)
Correction of persistent errors in Arabidopsis reference mitochondrial genomes
.
Plant Cell
,
30
,
525
527
.

42.

Parkin
I.A.
,
Koh
C.
,
Tang
H.
 et al.  (
2014
)
Transcriptome and methylome profiling reveals relics of genome dominance in the mesopolyploid Brassica oleracea
.
Genome Biol.
,
15
,
1
18
.

43.

Ming
R.
,
Hou
S.
,
Feng
Y.
 et al.  (
2008
)
The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus)
.
Nature
,
452
,
991
996
.

44.

Merchant
S.S.
,
Prochnik
S.E.
,
Vallon
O.
 et al.  (
2007
)
The Chlamydomonas genome reveals the evolution of key animal and plant functions
.
Science
,
318
,
245
250
.

45.

Blanc
G.
,
Agarkova
I.
,
Grimwood
J.
 et al.  (
2012
)
The genome of the polar eukaryotic microalga Coccomyxa subellipsoidea reveals traits of cold adaptation
.
Genome Biol.
,
13
,
1
12
.

46.

Li
Q.
,
Li
H.
,
Huang
W.
 et al.  (
2019
)
A chromosome-scale genome assembly of cucumber (Cucumis sativus L.)
.
GigaScience
,
8
,
1
10
.

47.

Polle
J.E.
,
Barry
K.
,
Cushman
J.
 et al.  (
2017
)
Draft nuclear genome sequence of the halophilic and beta-carotene-accumulating green alga Dunaliella salina strain CCAP19/18
.
Genome Announc.
,
5
,
01105
01117
.

48.

Mock
T.
,
Otillar
R.P.
,
Strauss
J.
 et al.  (
2017
)
Evolutionary genomics of the cold-adapted diatom Fragilariopsis cylindrus
.
Nature
,
541
,
536
540
.

49.

Schmutz
J.
,
Cannon
S.B.
,
Schlueter
J.
 et al.  (
2010
)
Genome sequence of the palaeopolyploid soybean
.
Nature
,
463
,
178
183
.

50.

Hanschen
E.R.
,
Marriage
T.N.
,
Ferris
P.J.
 et al.  (
2016
)
The Gonium pectorale genome demonstrates co-option of cell cycle regulation during the evolution of multicellularity
.
Nat. Commun.
,
7
,
1
10
.

51.

Hubert
O.
,
Piral
G.
,
Galas
C.
 et al.  (
2014
)
Changes in ethylene signaling and MADS box gene expression are associated with banana finger drop
.
Plant Sci.
,
223
,
99
108
.

52.

Sakai
H.
,
Lee
S.S.
,
Tanaka
T.
 et al.  (
2013
)
Rice Annotation Project Database (RAP-DB): an integrative and interactive database for rice genomics
.
Plant Cell Physiol.
,
54
,
1
11
.

53.

Verde
I.
,
Abbott
A.G.
,
Scalabrin
S.
 et al.  (
2013
)
The high-quality draft genome of peach (Prunus persica) identifies unique patterns of genetic diversity, domestication and genome evolution
.
Nat. Genet.
,
45
,
487
494
.

54.

Aoki
K.
,
Yano
K.
,
Suzuki
A.
 et al.  (
2010
)
Large-scale analysis of full-length cDNAs from the tomato (Solanum lycopersicum) cultivar Micro-Tom, a reference system for the Solanaceae genomics
.
BMC Genomics
,
11
,
1
16
.

55.

Diambra
L.A.
(
2011
)
Genome sequence and analysis of the tuber crop potato
.
Nature
,
475
,
189
195
.

56.

Argout
X.
,
Salse
J.
,
Aury
J.-M.
 et al.  (
2011
)
The genome of Theobroma cacao
.
Nat. Genet.
,
43
,
101
108
.

57.

Jaillon
O.
,
Aury
J.-M.
,
Noel
B.
 et al.  (
2007
)
The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla
.
Nature
,
449
,
463
467
.

58.

Prochnik
S.E.
,
Umen
J.
,
Nedelcu
A.M.
 et al.  (
2010
)
Genomic analysis of organismal complexity in the multicellular green alga Volvox carteri
.
Science
,
329
,
223
226
.

59.

Soderlund
C.
,
Descour
A.
,
Kudrna
D.
 et al.  (
2009
)
Sequencing, mapping, and analysis of 27,455 maize full-length cDNAs
.
PLoS Genet.
,
5
,
1
13
.

This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial License (https://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact journals.permissions@oup.com

Supplementary data