Prunus cerasus cv. 'Montmorency' Whole Genome v1.0 Assembly & Annotation
Goeckeritz, C. Z., Rhoades, K. E., Childs, K. L., Iezzoni, A. F., VanBuren, R., & Hollender, C. A. (2023). Genome of tetraploid sour cherry (Prunus cerasus L.) 'Montmorency' identifies three distinct ancestral Prunus genomes. Horticulture Research, 2023;, uhad097, DOI: https://doi.org/10.1093/hr/uhad097
Three collections of young leaf tissue were made in the spring of 2019 for PacBio sequencing, Illumina sequencing of a gDNA library, and Illumina sequencing of a Hi-C library. Canu was used to assemble approximately 100X PacBio long-reads into contigs. Subsequently, 56X coverage of Illumina 150bp paired-end gDNA reads were aligned to the assembly using bowtie2, and these alignments were used for polishing the initial assembly with Pilon. After polishing, the assembly was scaffolded using Illumina 150bp paired-end Hi-C reads and the Juicer + 3D-DNA pipeline. Misassemblies in the scaffolded chromosome-scale sequences were corrected in JuiceBox Assembly Tools. For annotation, Illumina 150bp paired-end RNA-sequencing was conducted on a variety of tissues from cultivar 'Montmorency' at different developmental time points to aid in gene discovery. In addition to the RNA-seq, Nanopore cDNA-sequencing was conducted on a diverse pool of tissues in an attempt to obtain full-length transcripts for annotation. Manually-curated protein databases were downloaded from Uniprot and arabidopsis.org and aligned to the assembly with Exonerate to make homology-based gene predictions. The RNA-seq and cDNA-sequencing data were aligned to the assembly using STAR and minimap2, respectively, and transcriptomes from both types of data were created using Stringtie. The two transcriptomes and protein alignments were given to MAKER for evidence-based gene predictions. After MAKER's first run, the resulting .gff3 was used to train gene finders Augustus and SNAP. MAKER was run a second time using these evidence-trained gene finders and the subsequent gene predictions were filtered so that each prediction contained at least one known Pfam domain. Then, deFusion was run to identify erroneously fused gene candidates, and more than 2500 gene predictions on the 16 main chromosomes of the assembly were defused as a result of this pipeline. At this point, predictions were again filtered to contain Pfam domains but also to exclude predictions with homology to transposeable elements. Lastly, Apollo was used to manually annotate the intron-exon structure of 67 genes, including all identified DAM (Dormancy Associated MADS-box) genes.
Summary of the ‘Montmorency’ assembly metrics
Homology of the Prunus cerasus Montmorency Genome v1.0 proteins was determined by pairwise sequence comparison using the blastp algorithm against various protein databases. An expectation value cutoff less than 1e-9 was used for the NCBI nr (Release 2021-09) and 1e-6 for the Arabidoposis proteins (Araport11), UniProtKB/SwissProt (Release 2021-09), and UniProtKB/TrEMBL (Release 2021-09) databases. The best hit reports are available for download in Excel format.
The Prunus cerasus Montmorency v1.0 genome gene prediction files are available in FASTA and GFF3 formats.
Functional annotation for the Prunus cerasus Montmorency Genome v1.0 are available for download below. The Prunus cerasus Montmorency Genome v1.0 proteins were analyzed using InterProScan in order to assign InterPro domains and Gene Ontology (GO) terms. Pathways analysis was performed using the KEGG Automatic Annotation Server (KAAS).
Transcript alignments were performed by the GDR Team of Main Bioinformatics Lab at WSU. The alignment tool 'BLAT' was used to map transcripts to the Prunus cerasusMontmorency genome assembly. Alignments with an alignment length of 97% and 97% identify were preserved. The available files are in GFF3 format.