Analysis of genomes, sequence alignments and topology analysis were conducted as reported previously
[58],
[154]. BLASTP and TBLASTN analyses to identify homologues of Ca
2+, Na
+ and non-selective cation channel subunits were carried out using the following human sequences (protein accession number in parentheses): full-length or pore sequences of IP
3R1 (Q14643.2; pore region residues 2536–2608) or RyR1 (P21817.3; pore region residues 4877–4948), and sequences of human TrpA1 (NP_015628; N-truncated sequence residues 765-end), TrpV1 (NP_061197; N-truncated sequence residues 430-end), TrpC1 (P48995; N-truncated sequence residues 350-end), CNGA1 (EAW93049; transmembrane sequence residues 200–420), CNGB1 (NP_001288), HCN2 (NP_001185.3; full-length, and TMD residues 200–470), NMDA receptor NR1 (Q05586), NMDA receptor N2 (Q12879), AMPA receptor GRIA1 (P42261.2), kainate receptor GRIK1 (P39086), nAChR-alpha1 (ABR09427), purinergic receptor P2X4 (NP_002551.2), pannexin-1 (AAH16931), Orai1 (NP_116179.2), STIM1 (AAH21300), TPC1 (NP_001137291.1), TPC2 (NP_620714.2), TrpP1 (NP_001009944), TrpP2 (NP_000288), TrpM1 (NP_002411), TrpML1 (NP_065394), CatSper1 (Q8NEC5.3), acid-sensing ion channel-1 (ASIC1) (P78348.3), mitochondrial uniporter (NP_612366.1), Ca
v1.2 (NP_955630.2), Na
v1.1 (NP_001189364), Piezo-1 (NP_001136336), Piezo-2 (NP_071351) and NALCN (AAH64343). Sequences of the
S. cerevisiae Ca
2+ channel Cch1 (CAA97244), Mid1 (NP_014108) and TrpY1 (NP_014730), as well as
Arabidopsis thaliana TPC1 (AAK39554) were also used to search for fungal homologues. The sequence of the MCU auxiliary subunit MICU1 (NP_006068.2) was also used. Searches to identify K
+ channel homologues were carried out using the following sequences of diverse human K
+ channels (protein accession number in parentheses): K
v1.2 (NP_004965.1), K
v7.1 (NP_000209.2) and K
v11.1 (hERG1) (Q12809.1); K
ir1.1 (ROMK1) (NP_000211.1), K
ir2.1 (IRK1) (NP_000882.1), K
ir3.1 (GIRK1) (NP_002230.1), K
ir4.1 (P78508.1), K
ir5.1 (Q9NPI9.1), K
ir6.1 (K
ATP1) (Q15842.1), K
ir6.2 (NP_000516.3) and K
ir7.1 (CAA06878.1); K
2P1.1 (TWIK1) (NP_002236.1), K
2P2.1 (TREK1) (NP_001017425.2), K
2P3.1 (TASK1) (NP_002237.1), K
2P13.1 (THIK1) (NP_071337.2), K
2P16.1 (TALK1) (NP_001128577.1) and K
2P18.1 (TRESK2) (NP_862823.1); K
Ca1.1 (BK) (NP_001154824.1), K
Ca2.1 (SK1) (NP_002239.2), K
Ca2.2 (SK2) (NP_067627), K
Ca3.1 (IK/SK4) (NP_002241.1) and K
Ca4.1 (SLACK/K
Na) (NP_065873.2). Other K
+ channel sequences were also used to search for fungal homologues, including: bacterial KcsA (P0A334), bacterial cyclic nucleotide-gated MlotiK1 (Q98GN8.1), archaeal depolarization-activated K
vAP (Q9YDF8.1), archaeal hyperpolarization-activated MVP (Q57603.1), archaeal Ca
2+-activated MthK (O27564.1), and TOK1 from
S. cerevisiae (CAA89386.1). Plant K
+ channel sequences were also used, including: the vacuolar outwardly rectifying, Ca
2+-regulated vacuolar two-pore TPK1 channel (NP_200374.1); vacuolar KCO3 (NP_001190480.1); the pollen plasma membrane TPK4 (NP_171752.1), the inward rectifier KAT1 (NP_199436.1), the outward rectifier SKOR (pore region of NP_186934.1, residues 271–340 to avoid ankyrin hits), and AKT1 (NP_180233.1). We also searched for homologues of Hv1 proton channel subunits (NP_115745.2). Default BLAST parameters for assessing statistical significance and for filtering were used (
ie. an Expect threshold of 10, and SEG filtering).
Several procedures ensured that hits were likely homologues of cation channel subunits. Firstly, the presence of multiple transmembrane domains was confirmed using TOPCONS
[155]. Secondly, reciprocal BLASTP searches (non-redundant protein database at NCBI) were made, using the identified fungal hits as bait, and only proteins that gave the original mammalian protein family as hits were analyzed further. Thirdly, the presence of conserved domains was confirmed using the Conserved Domains Database (NCBI). In addition, for homologues of K
+ channel subunits, only hits with regions of sequence similarity that encompassed the selectivity filter sequence of the K
+ channel subunit used as bait were acknowledged. Also, where possible, pore homology was confirmed by sequence alignment using ClustalW2.1 (European Bioinformatics Institute). Multiple sequence alignments were made using ClustalW2.1 and physiochemical residue colours are shown. Where shown, asterisks below the alignment indicate positions that have a single fully conserved residue, while colons indicate positions that have residues with highly similar properties (scoring >0.5 in the Gonnet PAM 250 matrix, ClustalW2). For phylogenetic analysis, multiple sequence alignments were made with MUSCLE v3.7 using default parameters. After using GBLOCKS at high stringency to remove regions of low confidence, and removing gaps, Maximum Likelihood analysis was carried out using PhyML v3.0 (WAG substitution model; 4 substitution rate categories; default estimated gamma distribution parameters; default estimated proportions of invariable sites; 100 bootstrapped data sets). Phylogenetic trees were depicted using TreeDyn (v198.3). MUSCLE, GBLOCKS, PhyML and TreeDyn were all functions of Phylogeny.fr (
http://www.phylogeny.fr/)
[156].