ncbi_utils.inc

Go to the documentation of this file.

<?php

function tpps_ncbi_get_taxonomy($organism, $rank) {
  // Subkingdom is a little different.
  if ($rank == 'subkingdom') {
    return tpps_ncbi_get_subkingdom($organism);
  }

  $id = tpps_ncbi_get_taxon_id($organism);
  $fetch = new EFetch('taxonomy');
  $fetch->addParam('id', "$id");
  $response = $fetch->get()->xml();
  if (empty($response->Taxon->LineageEx->Taxon)) {
    drupal_set_message(t('Unable to load Lineage information from NCBI for organism: @org', array('@org' => $organism)));
    return FALSE;
  }

  $lineage = $response->Taxon->LineageEx->Taxon;
  foreach ($lineage as $item) {
    if ($item->Rank == $rank) {
      return $item->ScientificName;
    }
  }
  drupal_set_message("Unable to determine $rank from NCBI for organism: $organism");
  return FALSE;
}

function tpps_ncbi_get_subkingdom($organism) {
  $id = tpps_ncbi_get_taxon_id($organism);
  $fetch = new EFetch('taxonomy');
  $fetch->addParam('id', "$id");
  $response = $fetch->get()->xml();
  if (empty($response->Taxon->Lineage)) {
    drupal_set_message(t('Unable to load Lineage information from NCBI for organism: @org', array('@org' => $organism)));
    return FALSE;
  }

  $lineage = $response->Taxon->Lineage;
  if (preg_match('/(angiosperm|gymnosperm)/', $lineage, $matches)) {
    if (!empty($matches[1])) {
      return $matches[1];
    }
  }
  drupal_set_message(t('Unable to determine subkingdom from NCBI for organism: @org', array('@org' => $organism)));
  return FALSE;
}

function tpps_ncbi_get_taxon_id($organism, $list = FALSE) {
  $formatted_org = implode('+', explode(' ', $organism));
  $search = new ESearch('taxonomy');
  $search->addParam('term', "{$formatted_org}[Scientific+Name]");

  $err_msg = 'Unable to load Taxon ID from NCBI for organism named: ' . $organism;
  try {
    $response = $search->get()->xml();
    if (empty($response->IdList->Id)) {
      drupal_set_message($err_msg);
      return FALSE;
    }
    if (!$list) {
      return $response->IdList->Id;
    }
    return $response->IdList;
  }
  catch (\Exception $e) {
    drupal_set_message($err_msg . ' EUtils Error: ' . $e->getMessage());
  }
}