F1000Research F1000Research 2046-1402 F1000 Research Limited London, UK 10.12688/f1000research.44714.1 Data Note Articles Draft genome assembly of the biofuel grass crop Miscanthus sacchariflorus

[version 1; peer review: 2 approved]

 De Vega Jose Conceptualization Formal Analysis Investigation Writing – Original Draft Preparation https://orcid.org/0000-0003-2847-5158 a 1 Donnison Iain Conceptualization Funding Acquisition Project Administration Writing – Review & Editing 2 Dyer Sarah Conceptualization Funding Acquisition Project Administration Writing – Review & Editing 1 Farrar Kerrie Conceptualization Funding Acquisition Project Administration Writing – Review & Editing 2 Earlham Institute, Norwich, NR4 7UZ, UK Institute of Biological, Environmental & Rural Sciences (IBERS) - Aberystwyth University, Aberystwyth, SY23 3EE, UK jose.de-vega@earlham.ac.uk

No competing interests were disclosed.

 18 1 2021 2021 10 29 11 1 2021 Copyright: © 2021 De Vega J et al. 2021 This is an open access article distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Miscanthus sacchariflorus (Maxim.) Hack. is a highly productive C4 perennial rhizomatous biofuel grass crop. M. sacchariflorus is among the most widely distributed species in the genus, particularly at cold northern latitudes, and is one of the progenitor species of the commercial M. × giganteus genotypes. We generated a 2.54 Gb whole-genome assembly of the diploid M. sacchariflorus cv. “Robustus 297” genotype, which represented ~59% of the expected total genome size. We later anchored this assembly using the chromosomes from the M. sinensis genome to generate a second assembly with improved contiguity. We annotated 86,767 and 69,049 protein-coding genes in the unanchored and anchored assemblies, respectively. We estimated our assemblies included ~85% of the M. sacchariflorus genes based on homology and core markers. The utility of the new reference for genomic studies was evidenced by a 99% alignment rate of the RNA-seq reads from the same genotype. The raw data, unanchored and anchored assemblies, and respective gene annotations are publicly available.

 Miscanthus biofuel C4 assembly annotation Biotechnology and Biological Sciences Research Council BBS/E/T/000PR9818 BBS/E/W/10963A01A This work was funded by core strategic funding from the Biotechnology and Biological Sciences Research Council (BBSRC) in projects BBS/E/T/000PR9818 and BBS/E/W/10963A01A. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Introduction

Miscanthus is a genus of C4 perennial rhizomatous grasses native to East Asia and Oceania, and naturally adapted to a wide range of climate zones and land types. Miscanthus sacchariflorus is among the most widely distributed species within the genus. It originated in the Yellow Sea region of China and can be predominantly found in cool latitudes of East Asia with varying ploidy 1 . M. sacchariflorus occurs in both diploid (2n = 38) and tetraploid (2n = 76) forms, where tetraploid M. sacchariflorus genotypes originated by autopolyploidy 2 . M. sacchariflorus probably has the greatest winter hardiness among all the Saccharinae 3 .

Natural interspecific Miscanthus hybrids are commonly observed, even between individuals of different ploidy. For example, introgression of M. sacchariflorus is often found among cultivated European M. sinensis ecotypes 1, 4 . Furthermore, M. x giganteus, a sterile triploid hybrid resulting from the hybridization between M. sinensis and M. sacchariflorus, is the predominant commercially grown species owing to its high biomass productivity and low chemical input requirements. The common occurrence of hybridization events and variable ploidy are challenging to the improvement of these bioenergy grasses and increase the need for genomic resources from different Miscanthus species. A chromosomal-scale reference genome using a doubled-haploid M. sinensis line was recently published 4 .

We assembled, annotated and validated a draft genome from the diploid M. sacchariflorus cv. “Robustus 297” genotype, as well as generating rhizome, stem and leaf RNA-Seq data from the same genotype. This dataset was previously used to verify that both M. sinensis and M. sacchariflorus share the same A/B ancestral tetraploidy 4 . Here, we present the first draft genome of M. sacchariflorus, the second Miscanthus genome available after M. sinensis 4 .

 Methods Plant materials and sequencing

DNA was extracted from leaves from the diploid M. sacchariflorus cv. “Robustus 297” genotype (Biosample SAMN08580354) using the Qiagen DNeasy kit. RNA was also extracted from leaf, stem and root tissues from the same plant. All samples were taken from a plant grown from seed in trays in a glasshouse in 2009. This genotype is established and used in breeding at IBERS (Wales, UK). The RNA-seq libraries were deposited as part of previous work in the BioProject PRJNA639832.

 Whole genome sequencing and assembly

We obtained ~5.86e9 pairs of 100 bp paired-end reads from an Illumina paired-end library with a 560 bp insert-size that was sequenced on Illumina HiSeq 2500 machines in rapid run mode by the Earlham Institute. This represents approximately 50X coverage of the heterozygous content and 100X coverage of the homozygous content of the genome. Read quality was assessed, and contaminants and adaptors removed using Kontaminant 5 . These paired-end short-reads were assembled into 17M contigs with a total length of 3.27 Gb using ABySS 6 version 1.5.1, with default options and a kmer size of 71.

We obtained ~141.1e6 pairs of reads from a Nextera 150 bp mate-pair library with approximately 7 Kb insert-size, which was used for scaffolding the previous contigs together with the paired-end reads, using SSPACE 7 without “extension” step. Nextera mate-pair reads were required to include a fragment of the adaptor to be used in the scaffolding step 5 , and we filtered out sequences shorter than 500 bp. We obtained 589K scaffolds, a total length of 2.54 Gb with an N50 of 10.2 Kb. This whole-genome assembly was denominated “Msac_v2” and is deposited at NCBI in BioProject PRJNA679435.

 Gene model and functional annotations

Our gene structure annotation pipeline 8 used five sources of evidence that were provided to AUGUSTUS 9 (version 2.7) for gene annotation: (1) Repetitive and low complexity regions of the scaffolds identified using RepeatMasker 10 (version open-4.0.5) based on homology with the RepBase 11 public database (Release 20140131) and a new database of repeat elements identified in the assembly with RepeatModeler 12 . The repeats annotation was deposited in Zenodo (See data availability); (2) exon-intron junctions identified by Tophat 13 (version 2.1.0); (3) de novo and genome-guided ab initio transcripts assembled with Trinity 14 (version 2.6.5 )and Cufflinks 15 (version 2.2.1) from RNA-Seq reads obtained from several tissues from the same genotype; (4) ab initio gene models predicted by SNAP 16 (version 29-11-2013) and GeneID 17 (version 1.4.4); and (5) homology-based alignments of transcripts and proteins from Miscanthus sinensis and maize using Exonerate 18 with a minimal identity of 0.7 and coverage of 0.7. Finally, AUGUSTUS 9 was run with the options “genemodel=complete” and “alternatives-from-evidence=true” to ensure that the predicted genes were compatible with all the previous provided evidence.

For the functional annotation of these predicted genes, translated gene sequences were compared with the NCBI non-redundant (nr 20170116) proteins and EBI’s InterPro (version 5.22.61) databases, and the results were imported into Blast2GO 19 to annotate the GO and GO slim terms, enzymatic protein codes and KEGG pathways. A similar GO annotation from translated gene sequences can be done with eggNOG-mapper 20 . These functional descriptors were deposited in Zenodo (See Underlying data).

 Anchoring the whole genome assembly using the Miscanthus sinensis reference

To improve the genome contiguity, we anchored our M. sacchariflorus scaffolds to the Miscanthus sinensis genome 4 . However, no nucleotide content from M. sinensis was incorporated in the M. sacchariflorus assemblies.

Firstly, scaffolds longer than 2 kbps from the whole genome assembly “Msac_v2” were scaffolded again using SSPACE 7 and the M. sinensis mate-pairs reads, the gaps between scaffolds were filled in with Ns. This new whole-genome assembly was denominated “Msac_v3”, and was deposited at NCBI in Bioproject PRJNA435476, under the GenBank accession GCA_002993905. It contains 137,916 scaffolds for a total of 2.074 Gb with an N50 of 25.6 Kbps. The gene annotation was projected to the “Msac_v3” assembly using PASA 21 (version 2.0.1): genes were aligned to the new assembly using GMAP, requiring a minimum identity of 0.85 and coverage of 0.55, and later validated using the default parameters in PASA.

Finally, we obtained the chromosomal position in the M. sinensis chromosomes of the scaffolds from the “Msac_v3” assembly. Using Satsuma2 22 (version untagged-330e3341a1151a978b37), we identified every perfect-identify match between both assemblies (3,635,504 matches in total). The coordinates of these matches in BED 8 format were used as input to the “OrderOrientBySynteny” script from Satsuma2, which identifies the best chromosomal position for each scaffold. These position coordinates are available as an AGP file as part of GCA_002993905, which anchors our final whole-genome assembly to 19 chromosomes (accessions CM00959 to CM009609 in NCBI).

 Completeness assessment

RNA-seq cleaned reads from each tissue were independently aligned to both assembly versions using STAR 23 (version 2.6.0c). BUSCO 24 (version 4.1.4) was used to assess completeness with the single-copy orthologs database for green plants (Viridiplantae, version 2020-09-10). Orthologs were identified using Orthofinder2 25 (version 2.3.12) with default parameters and the option “-msa”, which directly provided comprehensive statistics comparing the provided proteomes. All the proteomes from the other species used ( Table 1) were downloaded from Phytozome (v7.1 DOE-JGI). Genomes were aligned using Minimap2 26 (version 2.17) with the “asm10” parameter for related genomes, secondary alignments (tp:A:S) filtered out, and results visualised using dotPlotly 27 (Github version, latest updated on 4 May 2018).

Completeness statistics of the unanchored and anchored <italic toggle="yes">M. sacchariflorus</italic> whole-genome assemblies in comparison to the <italic toggle="yes">M. sinensis</italic> reference.
Msac_v2 (unanchored) Msac_v3 (anchored by M. sinensis) Reference: M. sinensis 4 .
NCBI bioproject PRJNA679435 PRJNA435476 (GCA_002993905) v.7.1 from Phytozome
Length 2.539 Gb 2.074 Gb 1.68 Gbps
Scaffolds 588,758 scaffolds 137,931 scaffolds * 19 Chrs and 14,414 scaffolds
N20 25.39 Kbps 62.61 Kbps 146.1 Mbps
N50 10.25 Kbps 25.63 Kbps 88.51 Mbps
N80 2.79 Kbps 9.42 Kbps 75.06 Mbps
Max 378.48 Kbps 458.83 Kbps 160.9 Mbps
ANNOTATION Msac_v2 Msac_v3 M. sinensis
Gene models 81,431 68,578 67,967
Proteins 86,767 68,578 ** 67,789
BUSCO Msac_v2 Msac_v3 M. sinensis
Complete 55.5% (48% in single copy) 59.8% (50.4% in single copy) 97.6% (36.2% in single copy)
Fragmented 32.2% 26.4% 1.6%
Missing 12.3% 13.6% 0.8%
RNA MAPPING Msac_v2 Msac_v3 M. sinensis
Reads mapping in the genome once (root, stem and leaf) 76.2% 76.4% 78.8% 75% 76.7% 78.1% 78.8% *** 83.5% 82.5%
Reads mapping in the genome multiple times (root, stem and leaf) 22.5% 23% 20.7% 19.5% 18.8% 17.3% 19.7% *** 15.5% 16.6%

*15 scaffolds from plastids were discarded during the deposit in NCBI resulting in 137,916 scaffolds. ** Only the longest transcript was considered in each projected locus. *** Cross-species alignments.

 Results

We produced two whole-genome assemblies for M. sacchariflorus that we named “Msac_v2” and “Msac_v3”, with total lengths of 2.54 Gbps and 2.074 Gbps, respectively ( Table 1). The difference in size is mainly a result of filtering 402 Mb from sequences under 2 kb in the latter before anchoring to the M. sinensis genome. Our “Msac_v2” assembly covered ~59 % of M. sacchariflorus genome size, which is estimated to be 4.3 Gb 28 . Approximately 40% of the assembly was composed by transposable elements (987.3 Mb; Table 2), including 491 Mb (19.4%) and 154 Mb (6.1%) by copies of the Gypsy and Copia LTRs, respectively; and 180 Mb (7.1%) by several class 2 DNA transposons (MULE, CMC, Harbinger, etc.)

Transposable elements identified in the Miscanthus sacchariflorus genome.
Category Superfamily Coverage(bp) Fraction (2.539 Gb)
Class 1 TEs: retrotransposons (copy and paste) Gypsy LTR 491,915,558 19.37%
Copia LTR 154,244,411 6.08%
Other LTRs 87,661,401 3.45%
SINEs 5,029,476 0.20%
LINEs 25,076,275 0.99%
Other non-LTR retrotransposons 29,192,841 1.15%
Class 2 TEs: DNA transposons (cut and paste) hAT 10,722,644 0.42%
Harbinger/PIF 24,553,614 0.97%
MULE/MuDR 29,733,691 1.17%
Stowaway/TcMar 14,112,359 0.56%
CMC_EnSpm 56,449,907 2.22%
Helitron 10,601,152 0.42%
Other 34,676,501 1.37%
Unclassified TEs Unclassified TEs 5,934,794 0.23%
Non TEs Satellites 5,339,464 0.21%
snRNAs 23,147 0.00%
TOTAL 985,267,235 38.81%

We identified 219,394 primary alignments longer than 2 kb between the unanchored M. sacchariflorus (“Msac_v2”) and M. sinensis. The resulting dotplot ( Figure 1) shows the conserved synteny between both species, which diverged 1.6 Mya 4 . Figure 1 also shows the highly conserved synteny between the pairs of homoeologous chromosomes (e.g. green boxes in chromosomes one and two), and the fusion in chromosome 7 of the chromosome homeolog to chromosome 13; which was also reported in M. sinensis 4 . There are several large inversions between chromosomes 9 and 10, and 3 and 4 (cyan boxes in Figure 1). Our assembly of a heterozygous genotype resulted in multiple heterotigs (heterozygous contigs) containing the alternative or secondary haplotypes (e.g. pink boxes in Figure 1).

Conserved synteny between <italic toggle="yes">M. sacchariflorus</italic> and <italic toggle="yes">M. sinensis</italic> genomes.

The plot shows the primary alignments longer than 2 kbps between both species. The M. sacchariflorus scaffolds (Y-axis) have been sorted by their coordinates in M. sinensis chromosomes (X-axis). Large homoeologous blocks and chromosomal rearrangements are highlighted in boxes.

The utility of our assemblies for genomic studies is evidenced by the proportion of RNA-seq from three different tissues from the same M. sacchariflorus genotype that aligned to the assemblies. On average 99% and 95% of the RNA-seq reads aligned in “Msac_v2” and “Msac_v3”, respectively ( Table 1).

We estimated that we assembled more than 85% of the M. sacchariflorus genes. Furthermore, our assemblies include several alleles of genes in the heterozygous regions of the genome, while the M. sinensis reference was generated from a double-haplotyped genotype. The estimation of the proportion of assembled genes (~85%) was supported by (1) the results from BUSCO, which reported 86.4–87.7% of presented core genes, of which ~2/3rds were complete ( Table 1); and (2) the difference in the number of proteins from related species for which we can identify an ortholog in M. sacchariflorus compared to M. sinensis, as control, using Orthofinder2 ( Table 3).

Number of orthologs between <italic toggle="yes">Miscanthus sinensis</italic> (Msin), <italic toggle="yes">Setaria italica</italic> (Sita; foxtail millet), <italic toggle="yes">Sorghum bicolor</italic> (Sbic; sorghum), <italic toggle="yes">Zea mays</italic> (Zma; maize), and <italic toggle="yes">Panicum virgatum</italic> (Pvi; switchgrass) obtained using Orthofinder 2.
Orthologs Msac_v2 Msac_v3 Msin Sita Sbic Zma Pvi
From Msac_v2 (86,767) - NA 44,151 (50.9%) 36,904 (42.5%) 37,219 (42.9%) 38,478 (44.3%) 45,792 (52.8%)
From Msac_v3 (68,578) NA - 38,122 (55.6%) 32,273 (47.1%) 32,296 (47.1%) 33,395 (48.7%) 38,755 (56.5%)
From Msin (67,789) 43,739 (64.5%) 37,501 (55%) - 41,532 (64.1%) 43,475 (64.1%) 39,986 (58.9%) 45,913 (67.7%)
From Sita (40,599) 26,846 (66.1%) 23,559 (58%) 28,473 (70.1%)
From Sbic (39,441) 27,877 (70.7%) 24,125 (61.2%) 30,907 (78.4%)
From Zma (88,760) 41,530 (46.8%) 35,955 (40.5%) 41,784 (47.1%)
From Pvi (125,439) 63,692 (50.8%) 56,120 (44.7%) 64,271 (51.2%)

Based on the results from Orthofinder2 ( Table 3), we found orthologs in M. sacchariflorus for 64.5% of the M.sinensis annotated proteins, so we estimated ~1/3 rd of the Miscanthus proteins to be specific to each species. On the other hand, we estimated that ~3,000 genes may be missing in the “Msac_v2” annotation based on the number of Sorghum bicolor proteins with orthologues in M. sinensis but absent in M. sacchariflorus. Better estimations were obtained with the other four species, where the genes absent in Msac_v2 compared with M. sinensis were estimated to be 254, 579 and 1627 ( Table 3). Additionally, ~6,000 genes could be missed in “Msac_v3” compared to “Msac_v2” based on the difference in the number of M. sinensis orthologues in each assembly ( Table 3). This is likely from genes in the sequences shorter than 2 Kbps (totalling 402 Mbps) that were filtered out before anchoring. There was a large difference in the proportion of “fragmented” BUSCO genes found in the M. sacchariflorus (32.2%) and M. sinensis (1.6%) assemblies ( Table 1). To assess if that difference had an effect on the quality of the annotation, we compared the number of proteins from M. sacchariflorus and M. sinensis for which we can identify an ortholog in another species ( Table 3); we found the difference between both Miscanthus species ranged between 6,571 proteins when compared to sorghum (43,475 to 37,219; Table 2) to only 121 when compared to maize (39,986 to 38,478, Table 3).

In conclusion, our M. sacchariflorus genome can served as the basis for functional genetic analyses on Miscanthus, one of the main biofuel grass crops used in temperate latitudes. However, there are opportunities to improve it using new approaches, such as long-reads.

 Data availability Underlying data

NCBI BioProject: Miscanthus sacchariflorus cultivar:Robustus 297. Accession number PRJNA435476; https://identifiers.org/bioproject:PRJNA435476.

This BioProject contains the raw paired-end and mate-pair reads.

NCBI BioProject: RNA-seq Miscanthus hybrids with contrasting phenotypes. Accession number PRJNA639832; https://identifiers.org/bioproject:PRJNA639832.

This BioProject contains RNA-seq reads, deposited as part of a previous project 29 .

NCBI BioProject: Miscanthus sacchariflorus cultivar:Robustus 297. Accession number PRJNA679435; https://identifiers.org/bioproject:PRJNA679435.

This Bioproject contains the unanchored “Msac_v2” assemblies and gene annotations under accession JADQCR000000000.

The anchored “Msac_v3” assemblies and gene annotations are deposited under accession accession GCA_002993905 under Bioproject PRJNA435476.

The chromosomal positions in the M. sinensis chromosomes of the scaffolds from the “Msac_v3” assembly are available in an AGP file as part of GCA_002993905, which places the scaffolds in 19 chromosomes (accessions CM009591 to CM009609 in NCBI).

Zenodo: Supplementary dataset to "Draft genome assembly of the biofuel grass crop Miscanthus sacchariflorus". http://doi.org/10.5281/zenodo.4270235.

This project contains the assemblies in FASTA format, gene annotations in GFF3 format, functional annotations in tabulated text format, and AGP file with anchoring information.

Data deposited with Zenodo are available under the terms of the Creative Commons Attribution 4.0 International license (CC-BY 4.0).

 Clark LV Jin X Petersen KK : Population structure of Miscanthus sacchariflorus reveals two major polyploidization events, tetraploid-mediated unidirectional introgression from diploid M. sinensis, and diversity centred around the Yellow Sea. Ann Bot. 2019;124(4):731–48. 30247525 10.1093/aob/mcy161 6821896 Dwiyanti MS Rudolph A Swaminathan K : Genetic analysis of putative triploid Miscanthus hybrids and tetraploid M. sacchariflorus collected from sympatric populations of Kushima, Japan. Bioenergy Res. 2013;6(2):486–93. 10.1007/s12155-012-9274-3 Clark LV Dzyubenko E Dzyubenko N : Ecological characteristics and in situ genetic associations for yield-component traits of wild Miscanthus from eastern Russia. Ann Bot. 2016;118(5):941–55. 27451985 10.1093/aob/mcw137 5055818 Mitros T Session AM James BT : Genome biology of the paleotetraploid perennial biomass crop Miscanthus. Nat Commun. 2020;11(1):5442. 33116128 10.1038/s41467-020-18923-6 7595124 Leggett RM Ramirez-Gonzalez RH Clavijo BJ : Sequencing quality assessment tools to enable data-driven informatics for high throughput genomics. Front Genet. 2013;4:288. 24381581 10.3389/fgene.2013.00288 3865868 Simpson JT Wong K Jackman SD : ABySS: a parallel assembler for short read sequence data. Genome Res. 2009;19(6):1117–23. 19251739 10.1101/gr.089532.108 2694472 Boetzer M Henkel CV Jansen HJ : Scaffolding pre-assembled contigs using SSPACE. Bioinformatics. 2011;27(4):578–9. 21149342 10.1093/bioinformatics/btq683 De Vega JJ Ayling S Hegarty M : Red clover ( Trifolium pratense L.) draft genome provides a platform for trait improvement. Sci Rep. 2015;5:17394. 26617401 10.1038/srep17394 4663792 Stanke M Keller O Gunduz I : AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res. 2006;34(Web Server issue):W435–W9. 16845043 10.1093/nar/gkl200 1538822 Tarailo‐Graovac M Chen N : Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinformatics. 2009;25(1): Chapter 4: Unit 4.10. 19274634 10.1002/0471250953.bi0410s25 Jurka J Kapitonov VV Pavlicek A : Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005;110(1–4):462–7. 16093699 10.1159/000084979 Smit AF Hubley R : RepeatModeler Open-1.0.2008. Reference Source Trapnell C Pachter L Salzberg SL : TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11. 19289445 10.1093/bioinformatics/btp120 2672628 Grabherr MG Haas BJ Yassour M : Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat Biotechnol. 2011;29(7):644–52. 21572440 10.1038/nbt.1883 3571712 Trapnell C Roberts A Goff L : Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78. 22383036 10.1038/nprot.2012.016 3334321 Bromberg Y Rost B : SNAP: predict effect of non-synonymous polymorphisms on function. Nucleic Acids Res. 2007;35(11):3823–35. 17526529 10.1093/nar/gkm238 1920242 Blanco E Parra G Guigó R : Using geneid to identify genes. Curr Protoc Bioinformatics. 2007; Chapter 4: Unit 4.3. 18428791 10.1002/0471250953.bi0403s18 Slater GSC Birney E : Automated generation of heuristics for biological sequence comparison. BMC Bioinformatics. 2005;6(1):31. 15713233 10.1186/1471-2105-6-31 553969 Conesa A Götz S García-Gómez JM : Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6. 16081474 10.1093/bioinformatics/bti610 Huerta-Cepas J Forslund K Coelho LP : Fast Genome-Wide Functional Annotation through Orthology Assignment by eggNOG-Mapper. Mol Biol Evol. 2017;34(8):2115–22. 28460117 10.1093/molbev/msx148 5850834 Haas BJ Delcher AL Mount SM : Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 2003;31(19):5654–66. 14500829 10.1093/nar/gkg770 206470 Clavijo B Wright J Yanes L . Reference Source Dobin A Davis CA Schlesinger F : STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21. 23104886 10.1093/bioinformatics/bts635 3530905 Simão FA Waterhouse RM Ioannidis P : BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31(19):3210–2. 26059717 10.1093/bioinformatics/btv351 Emms DM Kelly S : OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 2019;20(1):238. 31727128 10.1186/s13059-019-1832-y 6857279 Li H : Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018;34(18):3094–100. 29750242 10.1093/bioinformatics/bty191 6137996 Poorten T . Reference Source Rayburn AL Crawford J Rayburn CM : Genome Size of Three Miscanthus Species. Plant Mol Biol Report. 2009;27(2):184. 10.1007/s11105-008-0070-3 De Vega JJ Peel N Purdy SJ : Differential expression of starch and sucrose metabolic genes linked to varying biomass yield in Miscanthus hybrids. BioRxiv. 2020;2020–08. 10.1101/2020.08.04.236885 10.5256/f1000research.47786.r78040 Reviewer response for version 1 Dwiyanti Maria Stefanie 1 Referee https://orcid.org/0000-0001-9042-0612 Research Faculty of Agriculture, Hokkaido University, Sapporo, Japan

Competing interests: No competing interests were disclosed.

 1 3 2021 Copyright: © 2021 Dwiyanti MS 2021 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve

The availability of Miscanthus sacchariflorus genome sequence will be useful for Miscanthus related research, particularly in bioenergy related topics. 

"Better estimations were obtained with the other four species, where the genes absent in Msac_v2 compared with M. sinensis were estimated to be 254, 579 and 1627 (Table 3)."

I found that the difference between number of genes in "Msac_v2" compared to other four species is larger than 254, 579, and 1627; or the way I look into the table is wrong?

Perhaps the sentence above can be reworded so we can easily compare with the Table 3 content.

 Also, what are the predicted functions of genes absent in "Msac_v2" compared to M.sinensis?

This information may provide some clues to trait difference between M. sacchariflorus and M.sinensis.

Are sufficient details of methods and materials provided to allow replication by others?

Yes

Is the rationale for creating the dataset(s) clearly described?

Yes

Are the datasets clearly presented in a useable and accessible format?

Yes

Are the protocols appropriate and is the work technically sound?

Yes

Reviewer Expertise:

Plant genetics and genomics

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

 10.5256/f1000research.47786.r77622 Reviewer response for version 1 Riaño-Pachón Diego Mauricio 1 Referee https://orcid.org/0000-0001-9803-3465 Computational, Evolutionary and Systems Biology Laboratory, Center for Nuclear Energy in Agriculture (CENA), University of São Paulo, Piracicaba, Brazil

Competing interests: No competing interests were disclosed.

 25 2 2021 Copyright: © 2021 Riaño-Pachón DM 2021 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve

The Data Note, "Draft genome assembly of the biofuel grass crop  Miscanthus sacchariflorus", introduces two Miscanthus sacchariflorus genome assemblies, which have been deposited in NCBI under the bioprojects:  PRJNA679435 and  PRJNA435476. Genome sequencing was carried out with Illumina paired end reads and mate-pairs, the assemblies are greatly fragmented, which is expected due to the sequencing technologies used. This is the first Miscanthus sacchariflorus genome assembly, which is of interested for the bioenergy community, and can be used to generate insigths with the genomes of other bioenergy crops, like sorghum and sugarcane.

Suggestions:

Look for contaminant organisms in the final assemblies using BlobPlots.

Provide GenomeScope and Smudgeplots for the clean reads, to generate further statistics prior to assembly.

Are sufficient details of methods and materials provided to allow replication by others?

No

Is the rationale for creating the dataset(s) clearly described?

Yes

Are the datasets clearly presented in a useable and accessible format?

Yes

Are the protocols appropriate and is the work technically sound?

Yes

Reviewer Expertise:

Bioinformatics, genome assembly and annotation.

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.