I would like to fetch 1:1 orthologs for GBRA2_HUMAN, a cannonical GABAA receptor alpha-2 subunit sequence as defined by Uniprot (Uniprot ID: P47869), however it seems that this sequence is not available in the OMA database. Could you please let me know if there is an alternative way to obtain the orthologs for this gene and what might be the reason for the missing entry?
I tried using the OMAdb R package:
getHOG(id ='GBRA2_HUMAN', members = TRUE)
THE OMA REST API request failed:https://omabrowser.org/api/hog/GBRA2_HUMAN/members/ Here's the original error message: Not Found (HTTP 404).
If I use
mapSequence(sequence) using the cannonical sequence from Uniprot, I get three matched targets, neither of which is the cannonical human GABAA alpha2 subunit:
 "A0A2R8ZZ28" "H2QPE5" "G3QVC6"
On the other hand, other human GABAA subunits seem to be present in the database (e.g. GBRA1_HUMAN, GBRA3_HUMAN).
The GBRA2_HUMAN sequence I used for searches is the following:
sequence <- 'MKTKLNIYNMQFLLFVFLVWDPARLVLANIQEDEAKNNITIFTRILDRLLDGYDNRLRPG LGDSITEVFTNIYVTSFGPVSDTDMEYTIDVFFRQKWKDERLKFKGPMNILRLNNLMASK IWTPDTFFHNGKKSVAHNMTMPNKLLRIQDDGTLLYTMRLTVQAECPMHLEDFPMDAHSC PLKFGSYAYTTSEVTYIWTYNASDSVQVAPDGSRLNQYDLLGQSIGKETIKSSTGEYTVM TAHFHLKRKIGYFVIQTYLPCIMTVILSQVSFWLNRESVPARTVFGVTTVLTMTTLSISA RNSLPKVAYATAMDWFIAVCYAFVFSALIEFATVNYFTKRGWAWDGKSVVNDKKKEKASV MIQNNAYAVAVANYAPNLSKDPVLSTISKSATTPEPNKKPENKPAEAKKTFNSVSKIDRM SRIVFPVLFGTFNLVYWATYLNREPVLGVSP'