cgeb2001's DokuWiki!

You've loaded an old revision of the document! If you save it, you will create a new version with this data.

**What is HSDatabase?**

Backgrounds:

The idea of HSDatabase project comes from several important milestone articles in the psychrophilic algae research. The aims of the database is to develop a platform for comprehensive analysis of the eukaryotic species, particularly those who can survive in harsh environment, to unravel the role of gene duplication in genome adaptation.
  - Draft genome sequence of the Antarctic green alga Chlamydomonas sp. UWO241 (2021)
  - Adaptation to Extreme Antarctic Environments Revealed by the Genome of a Sea Ice Green Alga (2020)
  - A constitutive stress response is a result of low temperature growth in the Antarctic green alga Chlamydomonas sp. UWO241 (2021)
  - Photosynthetic Adaptation to Polar Life: Energy Balance, Photoprotection and Genetic Redundancy (2021)
  - HSDFinder: a BLAST-based strategy for identifying highly similar duplicated genes in eukaryotic genomes (2021)
  - HSDatabase – Identification and functional annotation of highly similar duplicated genes in eukaryotic genomes(2022)
  - Comprehensive analysis of gene duplications in eukaryotic genomes unravelling the convergent adaptive and nonadaptive evolution (2022)

So far, the first step via designing the HSDFinder tool has been reached after so many trails, the selected eukaryotic species have been collected into the HSDatabase. The comprehensive analysis is on the way.

**How to prepare the species files?**

As for the HSDatabase itself, we have offered the species request option which allows users to submit the  species within users' interest.

**How to document that data in HSDatabase?**

There are several necessary files to be documented into the database.
Request a new species
If you wish that a new species would be added in HSDatabase, please use the following form. The new species have to meet the following requirement:

* Peptide sequence (https://github.com/zx0223winner/isoform2one)
  * Blast all-against-all file 
  * InterProscan file
  * KEGG file

The HSDatabase is based on the data provided by the NCBI FTP site. If your species is stored in the FTP site, it will be a valuable help to provide us the FTP links to the peptide database. At least, a link to the species information is required.

**How to analyze the data from HSDFinder?**

Although there is no golden rule to distinguish partial duplicates from more complete ones, it is believed that the relative complete duplicates turn to have at least less than 50% amino acid length difference and same number and function of conserved domain.

* HSD_ Statistics.py (https://github.com/zx0223winner/HSDFinder/tree/master/Tutorial) is to calculate the number of HSDs, number of gene copies, non-redundant gene copies, capturing value, performance score,
     * True HSDs: are those HSDs with gene copies minimum length occupied more than half of the maximum length and have the same function and number of Pfam domains.
     * Incomplete HSDs: are those having different number of conserved domains (Pfam domains) or gene copies encoding the hypothetical proteins have the varied aa length more than 50% of each other.

* HSD_Categories.py is to calculate the gene copies within each group, i.e., 2-group is the HSD group only has two gene copies.
  * HSD_add_on.py is to merge a series of combo thresholds based on the formula:E + (D + (C + (B +A)))

* A = 90%_100aa+(90%_70aa+(90%_50aa+(90%_30aa+90%_10aa)))
     * B = 80%_100aa+(80%_70aa+(80%_50aa+(80%_30aa+80%_10aa)))
     * C = 70%_100aa+(70%_70aa+(70%_50aa+(70%_30aa+70%_10aa)))
     * D = 60%_100aa+(60%_70aa+(60%_50aa+(60%_30aa+60%_10aa)))
     * E = 50%_100aa+(50%_70aa+(50%_50aa+(50%_30aa+50%_10aa)))

<Last updated by Xi Zhang on Oct 6th,2021> Upcoming

Please fill all the letters into the box to prove you're human. P S G J F Please keep this field empty:

Edit summary

Note: By editing this page you agree to license your content under the following license: CC Attribution-Noncommercial-Share Alike 4.0 International