Gossypium gossypioides (D6) genome ISU_v1

Overview
Analysis NameGossypium gossypioides (D6) genome ISU_v1
MethodAbySS & Chromosomer (v. 2.0.1 & v August 2017)
Source (v1)
Date performed2021-02-17

Materials and Methods

DNA was extracted from leaves of D6-5. Samples were submitted to BGI. Genomics (Hong Kong) for Illumina library preparation and 2 x100 bp sequencing. Reads were trimmed and filtered with Trimmomatic v0.32 (Bolger et al. 2014) with the following options: 1) sequence adapter removal, 2) removal of leading and/or trailing bases when the quality score (Q) <28, 3) removal of bases after average Q < 28 (8 nt window) or single base quality <10, and 4) removal of reads <85 nt. Detailed parameters can be found at https://github.com/IGBB/D_Cottons_USDA, last accessed December 18, 2018.

About the assembly

Trimmed data were independently assembled for each species via ABySS v2.0.1 (Simpson et al. 2009), using every 5th kmer value from 40 through 100. A single assembly with the highest E-size (Salzberg et al. 2012) was selected for each species and subsequently annotated with MAKER v2.31.6 (Holt and Yandell 2011) using evidence from: 1) the NCBI G. raimondii EST database (Udall et al. 2006), 2) G. raimondii reference genome predicted proteins, as hosted by CottonGen.org (Paterson et al. 2012), and 3) three ab initio gene prediction programs, that is Genemark v4.30 (Borodovsky et al. 2003), SNAP v2013-11-29 (Korf 2004), and Augustus v3.0.3 (Stanke et al. 2006). Both the SNAP and Augustus models were trained using BUSCO v2.0 (Simão et al. 2015). Chromosomer version 0.1.3 (Tamazian et al. 2016), a reference-assisted scaffolder, was used to scaffold the selected assemblies against the gold standard G. raimondii genome. MAKER v2.31.6 (Holt and Yandell 2011) was used to transfer the previous annotations to the Chromosomer-based scaffolds by rerunning MAKER and using the transcripts from the original annotation as evidence. Assemblies are also available under NCBI PRJNA488266.

Summary of the Updated Assembly Numbers
Total sequence length 664,724,314
Total ungapped length 640,396,013
Gaps between scaffolds 0
Number of scaffolds 344,899
Scaffold N50 32,963,687
Scaffold L50 10
Number of contigs 590,898
Contig N50 3,979
Contig L50 41,251
Total number of chromosomes and plasmids 13
Number of component sequences (WGS or clone) 344,899

Publication

Grover, et al. Insights into the evolution of the New World diploid cottons (Gossypium, subgenus Houzingenia) based on genome sequencing. Genome Biol Evol. 2019 Jan 1;11(1):53-71.

 

Assembly

The Gossypium gossypioides D6 ISU Genome v1.0 assembly file is available in FASTA format.

Downloads

Chromosomes (FASTA file) Gossypium_gossypioides_D6_ISU_v1.0.fasta.gz

 

Functional Analysis

Functional annotation files for the Gossypium gossypioides ISU Genome v1.0 are available for download below. The Gossypium gossypioides ISU 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).

Downloads

GO assignments from InterProScan D6_ISU_v1_genes2GO.xlsx.gz
IPR assignments from InterProScan D6_ISU_v1_genes2IPR.xlsx.gz
Proteins mapped to KEGG Orthologs D6_ISU_v1_KEGG-orthologis.xlsx.gz
Proteins mapped to KEGG Pathways D6_ISU_v1_KEGG-pathways.xlsx.gz

 

Gene Predictions

The Gossypium gossypioides D6 ISU v1.0 genome gene prediction files are available in FASTA and GFF3 formats.

Downloads

Protein sequences  (FASTA file) Gossypium_gossypioides_D6_ISU_v1.0.proteins.fasta.gz
CDS  (FASTA file) Gossypium_gossypioides_D6_ISU_v1.0.cds.fasta.gz
RNA sequences  (FASTA file) Gossypium_gossypioides_D6_ISU_v1.0.rnas.fasta.gz
Genes (GFF3 file) Gossypium_gossypioides_D6_ISU_v1.0.genes.gff3.gz

 

Homology

Homology of the Gossypium gossypioides ISU 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 2018-05) and 1e-6  for the Arabidoposis proteins (Araport11), UniProtKB/SwissProt (Release 2019-01), and UniProtKB/TrEMBL (Release 2019-01) databases. The best hit reports are available for download in Excel format. 

 

Protein Homologs

G.gossypioides ISU Genome v1.0 proteins with NCBI nr homologs (EXCEL file) D6_ISU_v1_vs_nr.xlsx.gz
G.gossypioides ISU Genome v1.0 proteins with NCBI nr (FASTA file) D6_ISU_v1_vs_nr_hit.fasta.gz
G.gossypioides ISU Genome v1.0 proteins without NCBI nr (FASTA file) D6_ISU_v1_vs_nr_noHit.fasta.gz
G.gossypioides ISU Genome v1.0 proteins with arabidopsis (Araport11) homologs (EXCEL file) D6_ISU_v1_vs_tair.xlsx.gz
G.gossypioides ISU Genome v1.0 proteins with arabidopsis (Araport11) (FASTA file) D6_ISU_v1_vs_tair_hit.fasta.gz
G.gossypioides ISU Genome v1.0 proteins without arabidopsis (Araport11) (FASTA file) D6_ISU_v1_vs_tair_noHit.fasta.gz
G.gossypioides ISU Genome v1.0 proteins with SwissProt homologs (EXCEL file) D6_ISU_v1_vs_swissprot.xlsx.gz
G.gossypioides ISU Genome v1.0 proteins with SwissProt (FASTA file) D6_ISU_v1_vs_swissprot_hit.fasta.gz
G.gossypioides ISU Genome v1.0 proteins without SwissProt (FASTA file) D6_ISU_v1_vs_swissprot_noHit.fasta.gz
G.gossypioides ISU Genome v1.0 proteins with TrEMBL homologs (EXCEL file) D6_ISU_v1_vs_trembl.xlsx.gz
G.gossypioides ISU Genome v1.0 proteins with TrEMBL (FASTA file) D6_ISU_v1_vs_trembl_hit.fasta.gz
G.gossypioides ISU Genome v1.0 proteins without TrEMBL (FASTA file) D6_ISU_v1_vs_trembl_noHit.fasta.gz

 

Markers
Marker alignments were performed by the CottonGen Team of Main Bioinformatics Lab at WSU. The alignment tool 'BLAT' was used to map marker sequences from CottonGen to the Gossypium gossypioides ISU me assembly. Markers required 90% identity over 97% of their length. For SSRs & RFLPs, gap size was restricted to 1000bp or less with less than 2 gaps. For dbSNPs and Indels gap size was restricted to 2bp with less than 2 gaps. The available files are in GFF3 format. Markers available in CottonGen are linked to JBrowse.
 
CottonGen SNP markers mapped to genome G.gossypioides_ISU-D6_SNP
CottonGen RFLP markers mapped to genome G.gossypioides_ISU-D6_RFLP
CottonGen SSR markers mapped to genome G.gossypioides_ISU-D6_SSR
CottonGen InDel markers mapped to genome G.gossypioides_ISU-D6_InDel

 

Transcript Alignments
Transcript alignments were performed by the CottonGen Team of Main Bioinformatics Lab at WSU. The alignment tool 'BLAT' was used to map transcripts to the G. gossypioides genome assembly. Alignments with an alignment length of 97% and 90% identify were preserved. The available files are in GFF3 format.

 

G. arboreum CottonGen RefTrans v1 G.gossypioides_ISU-D6_g.arboreum_cottongen_reftransV1
G. barbadense CottonGen RefTrans v1 G.gossypioides_ISU-D6_G.barbadense_cottongen_reftransV1
G. hirsutum CottonGen RefTrans v1 G.gossypioides_ISU-D6_g.hirsutum_cottongen_reftransV1
G. raimondii CottonGen RefTrans v1 G.gossypioides_ISU-D6_g.raimondii_cottongen_reftransV1