Prunus ferganensis Whole Genome v2.0 Assembly & Annotation
The study has not been published yet. The initial title is named as below: Pan-genome analyses of peach and its wild relatives provide insights into the genetics of disease resistance and species adaptation. And the manuscript has been submitted to bioRxiv.
Peach (Prunus persica) is the third most produced fruit crop, and is widely cultivated in temperate and subtropical regions. Due to its small genome size, peach has been used as a model plant for comparative and functional genomic researches of the Rosaceae family. In 2013, a high-quality reference genome sequence of peach constructed with the Sanger whole-genome shotgun approach was released by International Peach Genome Initiative. It is well known that wild germplasm contributes a significant proportion of the genetic resources of major crop species, and significant phenotypic differences in fruit size, flavor, and stress tolerance were found among P. persica and its wild relatives, P. mira, P. davidiana, P. kansuensis, and P. ferganensis. It is necessary to study genetic variations of peach and its wild relatives from a broader perspective, such as pan-genome analyses. P. ferganensis is an attractive model for studying the differentiation of peach because it has a close relationship with P. persica.
Genome facts and statistics
The genome of P. davidiana was assembled using variety “Zhouxing Shan Tao” through Illumina platforms. Illumina reads from the wild species of P. persica were assembled using ALLPATHS-LG, and gaps in the assemblies were filled using GapCloser V1.12. Mate-paired reads were then used to generate scaffolds using SSPACE. Illumina reads from the wild specie were assembled using ALLPATHS-LG, and gaps in the assemblies were filled using GapCloser V1.12. Mate-paired reads were then used to generate scaffolds using SSPACE. In additon, PacBio SMRT reads were de novo assembled using FALCON. Approximately 26.54 Gb of PacBio SMRT reads were first pairwise compared, and the longest 60 coverage of subreads were selected as seeds to do error correction. All PacBio SMRT reads were mapped back to the assembled contigs with Blast and the Arrow program implemented in SMRT Link (PacBio) was used for error correction with default parameters. The Illumina paired-end reads were then mapped to the corrected contigs to perform the second round of error correction. To further improve the continuity of the assembly, SSPACE (v3.0) was used to build scaffolds using reads from all the mate pair libraries. FragScaff v1-1 was further applied to build superscaffolds using the barcoded sequencing reads. Finally, Hi-C data were used to correct superscaffolds and cluster the scaffolds into pseudochromosomes. A total of 90.53% of them were allocated to eight pseudochromosomes. The contig N50 and scaffold N50 sizes of the final assembly were 26.47 Mb and 28.25 Mb, respectively. Gene prediction was performed using a combination of homology, ab initio and transcriptome based approaches resulting in 28,587 protein-coding genes in P. ferganensis.
Homology of the Prunus ferganensis genome v2.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.
Functional annotation for the Prunus ferganensis genome v2.0 are available for download below. The Prunus ferganensis genome v2.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 ferganensis genome assembly. Alignments with an alignment length of 97% and 97% identify were preserved. The available files are in GFF3 format.