The Piezo channel is a mechano-sensitive complex component in the mammalian inner ear hair cell.

The inner ear is the hub where hair cells (HCs) transduce sound, gravity, and head acceleration stimuli to the brain. Hearing and balance rely on mechanosensation, the fastest sensory signals transmitted to the brain. The mechanoelectrical transducer (MET) channel is the entryway for the sound-balance-brain interface, but the channel-complex composition is not entirely known. Here, we report that the mouse utilizes Piezo1 (Pz1) and Piezo2 (Pz2) isoforms as MET-complex components. The Pz channels, expressed in HC stereocilia, and cell lines are co-localized and co-assembled with MET complex partners. Mice expressing non-functional Pz1 and Pz2 at the ROSA26 locus have impaired auditory and vestibular traits that can only be explained if the Pzs are integral to the MET complex. We suggest that Pz subunits constitute part of the MET complex and that interactions with other MET complex components yield functional MET units to generate HC MET currents.

The Piezo channel is central to the mechano-sensitive channel complex in the mammalian inner ear.

The inner ear is the hub where hair cells transduce sound, gravity, and head acceleration stimuli carried by neural codes to the brain. Of all the senses, hearing and balance, which rely on mechanosensation, are the fastest sensory signals transmitted to the central nervous system. The mechanoelectrical transducer (MET) channel in hair cells is the entryway for the sound-balance-brain interface, but the channel's composition has eluded biologists due to its complexity. Here, we report that the mouse utilizes Piezo1 (Pz1) and Piezo2 (Pz2) isoforms as central components of the MET complex. The Pz channel subunits are expressed in hair-cell stereocilia, are co-localized and co-assembled, and are essential components of the MET complex in vitro and in situ, including integration with the transmembrane channel (Tmc1/2) protein. Mice expressing non-functional Pz1 and Pz2, but not functional Pz1 at the ROSA26 locus under the control of hair-cell promoters, have impaired auditory and vestibular traits that can only be explained if Pz channel multimers are integral to the MET complex. We affirm that Pz protein subunits constitute MET channels and that functional interactions with components of the MET complex yield current properties resembling hair-cell MET currents. Our results demonstrate Pz is a MET channel component central to interacting with MET complex proteins. Results account for the MET channel pore and complex.

Label-free quantitative mass spectrometry analysis of differential protein expression in the developing cochlear sensory epithelium.

BACKGROUND: The sensory epithelium of the inner ear converts the mechanical energy of sound to electro-chemical energy recognized by the central nervous system. This process is mediated by receptor cells known as hair cells that express proteins in a timely fashion with the onset of hearing. METHODS: The proteomes of 3, 14, and 30 day-old mice cochlear sensory epithelia were revealed, using label-free quantitative mass spectrometry (LTQ-Orbitrap). Statistical analysis using a one-way ANOVA followed by Bonferroni's post-hoc test was used to show significant differences in protein expression. Ingenuity Pathway Analysis was used to observe networks of differentially expressed proteins, their biological processes, and associated diseases, while Cytoscape software was used to determine putative interactions with select biomarker proteins. These candidate biomarkers were further verified using Western blotting, while coimmunoprecipitation was used to verify putative partners determined using bioinformatics. RESULTS: We show that a comparison across all three proteomes shows that there are 447 differentially expressed proteins, with 387 differentially expressed between postnatal day 3 and 30. Ingenuity Pathway Analysis revealed ~ 62% of postnatal day 3 downregulated proteins are involved in neurological diseases. Several proteins are expressed exclusively on P3, including Parvin α, Drebrin1 (Drb1), Secreted protein acidic and cysteine rich (SPARC), Transmembrane emp24 domain-containing protein 10 (Tmed10). Coimmunoprecipitations showed that Parvin and SPARC interact with integrin-linked protein kinase and the large conductance calcium-activated potassium channel, respectively. CONCLUSIONS: Quantitative mass spectrometry revealed the identification of numerous differentially regulated proteins over three days of postnatal development. These data provide insights into functional pathways regulating normal sensory and supporting cell development in the cochlea that include potential biomarkers. Interacting partners of two of these markers suggest the importance of these complexes in regulating cellular structure and synapse development.

Protein Quantitation of the Developing Cochlea Using Mass Spectrometry.

Mass spectrometry-based proteomics allows for the measurement of hundreds to thousands of proteins in a biological system. Additionally, mass spectrometry can also be used to quantify proteins and peptides. However, observing quantitative differences between biological systems using mass spectrometry-based proteomics can be challenging because it is critical to have a method that is fast, reproducible, and accurate. Therefore, to study differential protein expression in biological samples labeling or label-free quantitative methods can be used. Labeling methods have been widely used in quantitative proteomics, however label-free methods have become equally as popular and more preferred because they produce faster, cleaner, and simpler results. Here, we describe the methods by which proteins are isolated and identified from cochlear sensory epithelia tissues at different ages and quantitatively differentiated using label-free mass spectrometry.

Ultrastructural Identification and Colocalization of Interacting Proteins in the Murine Cochlea by Post-Embedding Immunogold Transmission Electron Microscopy.

Verification of the presence and location of a protein within tissue can be accomplished by western blotting and immunohistochemistry, using either paraffin or frozen sections. Affinity purification by reciprocal coimmunoprecipitations using the tissue of interest can demonstrate the existence of an interacting pair of proteins. Ultimately, the ability to visualize the interaction at the cellular level is desired. Precise location(s) of interacting proteins in situ can be accomplished by ultrastructural localization with high-quality primary antibodies and small-particle-size Au-conjugated secondary antibodies. Visualization can be obtained with a transmission electron microscope fitted with a high-resolution camera permitting magnifications that exceed 2 × 10(5), and, to date, resolution capability of 20+ Mpixels, thus enabling localization of the target protein to within nanometers of the actual location. Here, we report the method by which immunolocalization at the level of the electron microscope is accomplished using the post-embedding technique, i.e., performing antibody labeling of proteins on ultrathin sections of tissue embedded in acrylic resin.

The large conductance calcium-activated potassium channel affects extrinsic and intrinsic mechanisms of apoptosis.

The large-conductance calcium-activated K(+) or BK channel underlies electrical signals in a number of different cell types. Studies show that BK activity can also serve to regulate cellular homeostasis by protecting cells from apoptosis resulting from events such as ischemia. Recent coimmunoprecipitation studies, combined with mass spectrometry, suggest putative protein partners that interact with BK to regulate intrinsic and extrinsic apoptotic pathways. This study tests two of those partners to determine the effects on these two signaling pathways. Through reciprocal coimmunoprecipitation (coIP) experiments, we show that BK interacts with p53 and fas-associated protein with death domain (FADD) in mouse brain and when overexpressed in a heterologous expression system, such as HEK293 cells. Moreover, coIP experiments with N- and C-terminal fragments reveal that FADD interacts with the C-terminus of BK, whereas p53 interacts with either the N- or the C-terminus. Immunolocalization studies show that BK colocalizes with p53 and FADD in the mitochondrion and plasmalemma, respectively. HEK cells that stably express BK are more resistant to apoptosis when p53 or FADD is overexpressed or when their intrinsic and extrinsic pathways are stimulated via mitomycin C or tumor necrosis factor-related apoptosis-inducing ligand (TRAIL), respectively. Moreover, when stimulating with TRAIL, caspase-8 activation decreases in BK-expressing cells. These data suggest that BK is part of a larger complex of proteins that protects against apoptosis by interacting with proapoptotic proteins, such as p53 and FADD.

Intrinsic disorder in the BK channel and its interactome.

The large-conductance Ca2+-activated K+ (BK) channel is broadly expressed in various mammalian cells and tissues such as neurons, skeletal and smooth muscles, exocrine cells, and sensory cells of the inner ear. Previous studies suggest that BK channels are promiscuous binders involved in a multitude of protein-protein interactions. To gain a better understanding of the potential mechanisms underlying BK interactions, we analyzed the abundance, distribution, and potential mechanisms of intrinsic disorder in 27 BK channel variants from mouse cochlea, 104 previously reported BK-associated proteins (BKAPS) from cytoplasmic and membrane/cytoskeletal regions, plus BK ß- and γ-subunits. Disorder was evaluated using the MFDp algorithm, which is a consensus-based predictor that provides a strong and competitive predictive quality and PONDR, which can determine long intrinsically disordered regions (IDRs). Disorder-based binding sites or molecular recognition features (MoRFs) were found using MoRFpred and ANCHOR. BKAP functions were categorized based on Gene Ontology (GO) terms. The analyses revealed that the BK variants contain a number of IDRs. Intrinsic disorder is also common in BKAPs, of which ∼ 5% are completely disordered. However, intrinsic disorder is very differently distributed within BK and its partners. Approximately 65% of the disordered segments in BK channels are long (IDRs) (>50 residues), whereas >60% of the disordered segments in BKAPs are short IDRs that range in length from 4 to 30 residues. Both α and γ subunits showed various amounts of disorder as did hub proteins of the BK interactome. Our analyses suggest that intrinsic disorder is important for the function of BK and its BKAPs. Long IDRs in BK are engaged in protein-protein and protein-ligand interactions, contain multiple post-translational modification sites, and are subjected to alternative splicing. The disordered structure of BK and its BKAPs suggests one of the underlying mechanisms of their interaction.

Bottom-up and shotgun proteomics to identify a comprehensive cochlear proteome.

Proteomics is a commonly used approach that can provide insights into complex biological systems. The cochlear sensory epithelium contains receptors that transduce the mechanical energy of sound into an electro-chemical energy processed by the peripheral and central nervous systems. Several proteomic techniques have been developed to study the cochlear inner ear, such as two-dimensional difference gel electrophoresis (2D-DIGE), antibody microarray, and mass spectrometry (MS). MS is the most comprehensive and versatile tool in proteomics and in conjunction with separation methods can provide an in-depth proteome of biological samples. Separation methods combined with MS has the ability to enrich protein samples, detect low molecular weight and hydrophobic proteins, and identify low abundant proteins by reducing the proteome dynamic range. Different digestion strategies can be applied to whole lysate or to fractionated protein lysate to enhance peptide and protein sequence coverage. Utilization of different separation techniques, including strong cation exchange (SCX), reversed-phase (RP), and gel-eluted liquid fraction entrapment electrophoresis (GELFrEE) can be applied to reduce sample complexity prior to MS analysis for protein identification.

Direct control of Na(+)-K(+)-2Cl(-)-cotransport protein (NKCC1) expression with aldosterone.

Sodium/potassium/chloride cotransporter (NKCC1) proteins play important roles in Na(+) and K(+) concentrations in key physiological systems, including cardiac, vascular, renal, nervous, and sensory systems. NKCC1 levels and functionality are altered in certain disease states, and tend to decline with age. A sensitive, effective way of regulating NKCC1 protein expression has significant biotherapeutic possibilities. The purpose of the present investigation was to determine if the naturally occurring hormone aldosterone (ALD) could regulate NKCC1 protein expression. Application of ALD to a human cell line (HT-29) revealed that ALD can regulate NKCC1 protein expression, quite sensitively and rapidly, independent of mRNA expression changes. Utilization of a specific inhibitor of mineralocorticoid receptors, eplerenone, implicated these receptors as part of the ALD mechanism of action. Further experiments with cycloheximide (protein synthesis inhibitor) and MG132 (proteasome inhibitor) revealed that ALD can upregulate NKCC1 by increasing protein stability, i.e., reducing ubiquitination of NKCC1. Having a procedure for controlling NKCC1 protein expression opens the doors for therapeutic interventions for diseases involving the mis-regulation or depletion of NKCC1 proteins, for example during aging.

In-depth proteomic analysis of mouse cochlear sensory epithelium by mass spectrometry.

Proteomic analysis of sensory organs such as the cochlea is challenging due to its small size and difficulties with membrane protein isolation. Mass spectrometry in conjunction with separation methods can provide a more comprehensive proteome, because of the ability to enrich protein samples, detect hydrophobic proteins, and identify low abundant proteins by reducing the proteome dynamic range. GELFrEE as well as different separation and digestion techniques were combined with FASP and nanoLC-MS/MS to obtain an in-depth proteome analysis of cochlear sensory epithelium from 30-day-old mice. Digestion with LysC/trypsin followed by SCX fractionation and multiple nanoLC-MS/MS analyses identified 3773 proteins with a 1% FDR. Of these, 694 protein IDs were in the plasmalemma. Protein IDs obtained by combining outcomes from GELFrEE/LysC/trypsin with GELFrEE/trypsin/trypsin generated 2779 proteins, of which 606 additional proteins were identified using the GELFrEE/LysC/trypsin approach. Combining results from the different techniques resulted in a total of 4620 IDs, including a number of previously unreported proteins. GO analyses showed high expression of binding and catalytic proteins as well as proteins associated with metabolism. The results show that the application of multiple techniques is needed to provide an exhaustive proteome of the cochlear sensory epithelium that includes many membrane proteins. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium with the data set identifier PXD000231.

The large conductance calcium-activated K(+) channel interacts with the small GTPase Rab11b.

The transduction of sound by the receptor or hair cells of the cochlea leads to the activation of ion channels found in the basal and lateral regions of these cells. Thus, the processing of these transduced signals to the central nervous system is tied to the regulation of baso-lateral ion channels. The large conductance calcium-activated potassium or BK channel was revealed to interact with the small GTPase, Rab11b, which is one of many Rabs found in various endosomal pathways. Immunoelectron microscopy showed the colocalization of these two proteins in receptor cells and auditory neurons. Using Chinese hamster ovary cells as a heterologous expression system, Rab11b increased or decreased BK expression, depending on the overexpression or RNAi knockdown of Rab, respectively. Additional mutation analyses, using a yeast two-hybrid assay, suggested that this GTPase moderately interacts within a region of BK exclusive of the N- or C-terminal tails. These data suggest that this small GTPase regulates BK in a slow recycling process through the endocytic compartment and to the plasmalemma.

Conserved BK channel-protein interactions reveal signals relevant to cell death and survival.

The large-conductance Ca(2+)-activated K(+) (BK) channel and its ß-subunit underlie tuning in non-mammalian sensory or hair cells, whereas in mammals its function is less clear. To gain insights into species differences and to reveal putative BK functions, we undertook a systems analysis of BK and BK-Associated Proteins (BKAPS) in the chicken cochlea and compared these results to other species. We identified 110 putative partners from cytoplasmic and membrane/cytoskeletal fractions, using a combination of coimmunoprecipitation, 2-D gel, and LC-MS/MS. Partners included 14-3-3γ, valosin-containing protein (VCP), stathmin (STMN), cortactin (CTTN), and prohibitin (PHB), of which 16 partners were verified by reciprocal coimmunoprecipitation. Bioinformatics revealed binary partners, the resultant interactome, subcellular localization, and cellular processes. The interactome contained 193 proteins involved in 190 binary interactions in subcellular compartments such as the ER, mitochondria, and nucleus. Comparisons with mice showed shared hub proteins that included N-methyl-D-aspartate receptor (NMDAR) and ATP-synthase. Ortholog analyses across six species revealed conserved interactions involving apoptosis, Ca(2+) binding, and trafficking, in chicks, mice, and humans. Functional studies using recombinant BK and RNAi in a heterologous expression system revealed that proteins important to cell death/survival, such as annexinA5, γ-actin, lamin, superoxide dismutase, and VCP, caused a decrease in BK expression. This revelation led to an examination of specific kinases and their effectors relevant to cell viability. Sequence analyses of the BK C-terminus across 10 species showed putative binding sites for 14-3-3, RAC-α serine/threonine-protein kinase 1 (Akt), glycogen synthase kinase-3ß (GSK3ß) and phosphoinositide-dependent kinase-1 (PDK1). Knockdown of 14-3-3 and Akt caused an increase in BK expression, whereas silencing of GSK3ß and PDK1 had the opposite effect. This comparative systems approach suggests conservation in BK function across different species in addition to novel functions that may include the initiation of signals relevant to cell death/survival.

Identification and quantification of full-length BK channel variants in the developing mouse cochlea.

Maxi-K(+) (BK) channel diversity is attributed to alternative splicing in the kcnma1 gene. The resultant variants manifest themselves in different cell types, tissues, and functions, such as excitation, metabolism, and signaling. Immunoelectron microscopy revealed immunogold particle labeling of BK in apical and basal regions of inner and outer hair cells, respectively. Additional labeling occurs in Deiters' cells and the inner mitochondrial membrane. Identification of full-length sequences reveals 27 BK variants from embryonic and postnatal mouse inner ear, per classification by tail motif, VYR, DEC, and ERL, and by exon usage. Three predicted start codons are found encoding MAN, MSS, and MDA, of which MDA shows the greatest expression through all stages in development, whereas MAN is undetectable. Complex splice sites occur between exons 9 and 10 and between 21 and 23. Spliced-in/out exons between 8 and 10 reveal a short fragment composed of exons 8 + 10, detectable on postnatal day (PD) 14 and PD30, and a longer fragment composed of exons 8 + 9 + 10 that is upregulated on embryonic day (ED) 14. Spliced-in exons 22 or 23 are expressed on ED14 but decrease over time; however, exon 22 increases again on PD34. Using tail-specific primers, qRT-PCR from ED14, PD4, -14, and -30 shows that BK-VYR and -ERL dominate expression on ED14, whereas DEC dominates after birth in all cochlear regions. The localization of BK and the changes in expression of its exons and tail types, by alternative splicing during development, may contribute to cochlear organization, acquisition of hearing, and intracellular signaling.

The large-conductance Ca(2+)-activated K(+) channel interacts with the apolipoprotein ApoA1.

Owing to the multifaceted functions of the large conductance Ca(2+)-activated K(+) channel (BK), identification of protein-protein interactions is essential in determining BK regulation. A yeast two-hybrid screening of a cochlear cDNA library revealed a BK-ApoA1 interaction. Patch clamp recordings of excised membrane patches from transfected HEK293 cells showed that ApoA1 inhibits the BK alpha-subunit by significantly increasing activation and deactivation times, and shifting half-activation voltage to more positive potentials. Reciprocal coimmunoprecipitations verified the BK-ApoA1 interaction using excised sensory epithelium and ganglia. Additionally, immunocolocalization studies revealed BK and ApoA1 expression in both receptor cells and auditory neurons. These data suggest new avenues of investigation, given the importance of apolipoproteins in neurological diseases.

A protein interaction network for the large conductance Ca(2+)-activated K(+) channel in the mouse cochlea.

The large conductance Ca(2+)-activated K(+) or BK channel has a role in sensory/neuronal excitation, intracellular signaling, and metabolism. In the non-mammalian cochlea, the onset of BK during development correlates with increased hearing sensitivity and underlies frequency tuning in non-mammals, whereas its role is less clear in mammalian hearing. To gain insights into BK function in mammals, coimmunoprecipitation and two-dimensional PAGE, combined with mass spectrometry, were used to reveal 174 putative BKAPs from cytoplasmic and membrane/cytoskeletal fractions of mouse cochlea. Eleven BKAPs were verified using reciprocal coimmunoprecipitation, including annexin, apolipoprotein, calmodulin, hippocalcin, and myelin P0, among others. These proteins were immunocolocalized with BK in sensory and neuronal cells. A bioinformatics approach was used to mine databases to reveal binary partners and the resultant protein network, as well as to determine previous ion channel affiliations, subcellular localization, and cellular processes. The search for binary partners using the IntAct molecular interaction database produced a putative global network of 160 nodes connected with 188 edges that contained 12 major hubs. Additional mining of databases revealed that more than 50% of primary BKAPs had prior affiliations with K(+) and Ca(2+) channels. Although a majority of BKAPs are found in either the cytoplasm or membrane and contribute to cellular processes that primarily involve metabolism (30.5%) and trafficking/scaffolding (23.6%), at least 20% are mitochondrial-related. Among the BKAPs are chaperonins such as calreticulin, GRP78, and HSP60 that, when reduced with siRNAs, alter BKalpha expression in CHO cells. Studies of BKalpha in mitochondria revealed compartmentalization in sensory cells, whereas heterologous expression of a BK-DEC splice variant cloned from cochlea revealed a BK mitochondrial candidate. The studies described herein provide insights into BK-related functions that include not only cell excitation, but also cell signaling and apoptosis, and involve proteins concerned with Ca(2+) regulation, structure, and hearing loss.

PPTX, a pentraxin domain-containing protein, interacts with the T1 domain of K v 4.

Voltage-gated K(+) (K(v)) channels reside as tetramers in the membrane. The events that coordinate folding, trafficking, and tetramerization are mediated by an array of associated proteins and phospholipids whose identification is vital to understanding the dynamic nature of channel expression and activity. An interaction between an A-type K(+) channel, K(v)4.2, and a protein containing a pentraxin domain (PPTX) was demonstrated in the cochlea (Duzhyy et al. [ 2005] J. Biol. Chem. 280:15165-15172). Here, we present results based on fold recognition and homology modeling that revealed the tetramerization (T1) domain of K(v)4.2 as a potential docking site for interacting proteins such as PPTX. By using this model, putative sites were experimentally tested with the yeast two-hybrid system to assay interactions between PPTX and the T1 domain of K(v)4.2 wild type (wt) and mutants (mut). Results showed that amino acid residues 86 and 118 in the T1 domain are essential for interaction, because replacing these negatively charged with neutrally charged amino acids inhibits interactions. Cotransfections of Chinese hamster ovary cells with PPTX and K(v)4.2wt further revealed that PPTX increases K(v)4.2 wt expression in vitro when analyzing total lysates, whereas interactions with K(v)4.2 microt resulted in a decrease. These studies suggest that portions of the T1 domain can act as docking sites for proteins such as PPTX, further underscoring the significance of this domain.

The use of 2-D gels to identify novel protein-protein interactions in the cochlea.

Functional proteomics comprises a wide range of technologies for the identification of novel protein-protein interactions and biological markers. Studies of protein-protein interactions have gained from the development of techniques and technologies such as immunoprecipitation, preparative two-dimensional (2-D) gel electrophoresis for peptide mass fingerprinting (PMF), using matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry, and liquid chromatography-tandem mass spectrometry (LC-MS/MS). These applications enabled the discovery of putative protein partners without a priori knowledge of which one(s) might be relevant. Here, we report the methods by which membrane proteins are isolated from cochlear tissues and prepared for identification by mass spectrometry techniques.

In vivo verification of protein interactions in the inner ear by coimmunoprecipitation.

Genomics has provided us with vast amounts of data and thus, the challenge to identify and characterize gene products. Proteomics analysis, using methods such as yeast two-hybrid screenings, isoelectric focusing, and mass spectroscopy, generate potentially useful information. To determine functional relationships between and among proteins, however, the initial data for putative protein interactions must first be validated. One technique, which is considered the gold standard, is coimmunoprecipitation.

Development and regeneration of hair cells share common functional features.

The structural phenotype of neural connections in the auditory brainstem is sculpted by spontaneous and stimulus-induced neural activities during development. However, functional and molecular mechanisms of spontaneous action potentials (SAPs) in the developing cochlea are unknown. Additionally, it is unclear how regenerating hair cells establish their neural ranking in the constellation of neurons in the brainstem. We have demonstrated that a transient Ca(2+) current produced by the Ca(v)3.1 channel is expressed early in development to initiate spontaneous Ca(2+) spikes. Ca(v)1.3 currents, typical of mature hair cells, appeared later in development. Moreover, there is a surprising disappearance of the Ca(v)3.1 current that coincides with the attenuation of the transient Ca(2+) current as the electrical properties of hair cells transition to the mature phenotype. Remarkably, this process is recapitulated during hair-cell regeneration, suggesting that the transient expression of Ca(v)3.1 and the ensuing SAPs are signatures of hair cell development and regeneration.

Ion channel gene expression in the inner ear.

The ion channel genome is still being defined despite numerous publications on the subject. The ion channel transcriptome is even more difficult to assess. Using high-throughput computational tools, we surveyed all available inner ear cDNA libraries to identify genes coding for ion channels. We mapped over 100,000 expressed sequence tags (ESTs) derived from human cochlea, mouse organ of Corti, mouse and zebrafish inner ear, and rat vestibular end organs to Homo sapiens, Mus musculus, Danio rerio, and Rattus norvegicus genomes. A survey of EST data alone reveals that at least a third of the ion channel genome is expressed in the inner ear, with highest expression occurring in hair cell-enriched mouse organ of Corti and rat vestibule. Our data and comparisons with other experimental techniques that measure gene expression show that every method has its limitations and does not per se provide a complete coverage of the inner ear ion channelome. In addition, the data show that most genes produce alternative transcripts with the same spectrum across multiple organisms, no ion channel gene variants are unique to the inner ear, and many splice variants have yet to be annotated. Our high-throughput approach offers a qualitative computational and experimental analysis of ion channel genes in inner ear cDNA collections. A lack of data and incomplete gene annotations prevent both rigorous statistical analyses and comparisons of entire ion channelomes derived from different tissues and organisms.