This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Information like Genome size, gene count, chromosomes numbers, contigs and scaffolds etc using a list of plant species

I have a list of plant species for which I want to get information like:

Genome size; protein coding gene; chromosomes numbers; contigs and scaffolds

Is there any easiest way for this (programmatically or via tools)

genome
• 2,210 views • link
updated by Ram Ram • written by Nelo • Nelo

1 answer

For genome sizes, you might be interested the Plant c-value database https://cvalues.science.kew.org/

When it comes to assembly values I would recommend looking into NCBI and find the representative RefSeq or Genbank genome per species if it exists. This could be automated via Entrez e-utils. Sometimes there is information about the number of contigs also, but this information can be unavailable. If the genome is annotated, you can get information about the number of protein-coding genes from the annotation files, but sometimes such annotation is unavailable. Then you can look for any publications, but don't trust any gene numbers given in publications if you have the annotation file also. In the end, you will have to do this semi-automated and curate manually because the process is error-prone.

• 0 views • link
written by Michael Dondrup Michael Dondrup

Hello Could you please write down the Entrez e-utils command

• 0 views • link
written by Nelo • Nelo

Here is a simple example in e-utils that will fetch simple assembly stats for a species, in this case Medicago truncatula.

The code needs to be run in the linux shell, eg. under bash The first line fetches the accession of the representative assembly for that species, the second line fetches the stats for that assembly.

ACC=$(esearch -db genome -query "Medicago truncatula"[Organism:exp] | efetch -format docsum | xtract -pattern DocumentSummary  -element Assembly_Accession)
esearch -query $ACC -db assembly | efetch -format docsum | xtract -pattern DocumentSummary -block Stats -tab "\n"  -element Stat@category Stat

output:

alt_loci_count  chromosome_count    contig_count    contig_l50  contig_n50  non_chromosome_replicon_count   replicon_count  scaffold_count  scaffold_count  scaffold_count  scaffold_count  scaffold_l50    scaffold_n50    total_length    ungapped_length
0   8   64  7   23305320    2   10  42  10  0   32  4   56236587    430008013   429433753

You may have to modify that script towards your needs.

• 0 views • link
written by Michael Dondrup Michael Dondrup

Log in to answer this question.

More posts like this