Structural and regulatory evolution of cellular electrophysiological systems

Competition for finite resources as coordination mechanism for morphogenesis: An evolutionary algorithm study of digital embryogeny

The standard view of embryogenesis is one of cooperation driven by the cells' shared genetics and evolutionary interests. However, numerous examples from developmental biology and agriculture reveal a surprising amount of competition among body cells, tissues, and organs for both metabolic and informational resources. To explain the existence of such competition we had hypothesized that evolution uses limiting "reservoirs" of resource molecules as a communication medium - a global scratchpad, to enable tissues across the body to coordinate growth. Here, we test this hypothesis via an evolutionary simulation of embryogeny in silico. Genomes encode state transition rules for cells, such as proliferation, differentiation, and resource use, enabling virtual embryos to develop a specific large-scale morphology. An evolutionary algorithm operates over the genomes, with fitness defined as a function of specific morphological requirements for the final embryo shape. We found that not only does such an algorithm rapidly discover rules for cellular behavior that reliably make embryos with specific anatomical properties, but that it discovers the strategy of using finite resources to coordinate development. Given the option of using finite or infinite reservoirs (which determine cells' ability to carry out specific actions), evolution preferentially uses finite reservoirs, which results in higher fitness and increased consistency (without needing direct selection for morphological invariance). We report aspects of anatomical, physiological/transcriptional, and genomic analysis of evolved virtual embryos that help understand how evolution can use competition among genetically identical subunits within a multicellular body to coordinate reliable, complex morphogenesis. Our results suggest that under some conditions, composite multi-scale systems will promote conflict and artificial scarcity for their components.

OUP accepted manuscript

Animals rely on their sensory systems to inform them of ecologically relevant environmental variation. In the Southern Ocean, the thermal environment has remained between −1.9 and 5 °C for 15 Myr, yet we have no knowledge of how an Antarctic marine organism might sense their thermal habitat as we have yet to discover a thermosensitive ion channel that gates (opens/closes) below 10 °C. Here, we investigate the evolutionary dynamics of transient receptor potential (TRP) channels, which are the primary thermosensors in animals, within cryonotothenioid fishes—the dominant fish fauna of the Southern Ocean. We found cryonotothenioids have a similar complement of TRP channels as other teleosts (∼28 genes). Previous work has shown that thermosensitive gating in a given channel is species specific, and multiple channels act together to sense the thermal environment. Therefore, we combined evidence of changes in selective pressure, gene gain/loss dynamics, and the first sensory ganglion transcriptome in this clade to identify the best candidate TRP channels that might have a functional dynamic range relevant for frigid Antarctic temperatures. We concluded that TRPV1a, TRPA1b, and TRPM4 are the likeliest putative thermosensors, and found evidence of diversifying selection at sites across these proteins. We also put forward hypotheses for molecular mechanisms of other cryonotothenioid adaptations, such as reduced skeletal calcium deposition, sensing oxidative stress, and unusual magnesium homeostasis. By completing a comprehensive and unbiased survey of these genes, we lay the groundwork for functional characterization and answering long-standing thermodynamic questions of thermosensitive gating and protein adaptation to low temperatures.

Transmural gradients in ion channel and auxiliary subunit expression

Evolution has acted to shape the action potential in different regions of the heart in order to produce a maximally stable and efficient pump. This has been achieved by creating regional differences in ion channel expression levels within the heart as well as differences between equivalent cardiac tissues in different species. These region- and species-dependent differences in channel expression are established by regulatory evolution, evolution of the regulatory mechanisms that control channel expression levels. Ion channel auxiliary subunits are obvious targets for regulatory evolution, in order to change channel expression levels and/or modify channel function. This review focuses on the transmural gradients of ion channel expression in the heart and the role that regulation of auxiliary subunit expression plays in generating and shaping these gradients.

Allometric scaling of electrical excitation and propagation in the mammalian heart

Variations in body mass impose constraints on the structure and function of mammalian species, including those of the cardiovascular system. Numerous biological processes, including cardiovascular parameters, have been shown to scale with body mass (BM) according to the law of allometric scaling: Y=Y =a∙BM^b (Y, biological process; a, normalization constant; b, scaling exponent, which in many instances is a multiple of ¼). These parameters include heart and breathing rates, intervals and subintervals of the electrocardiogram (ECG), action potential duration (APD), metabolic rate, and temporal properties of ventricular fibrillation. For instance, the hierarchical branching networks of the vascular system, and of the specialized conduction system in the heart have been proposed to be important determinants of allometric scaling. A global and unifying molecular mechanism of allometric scaling has not been put forth, but changes in gene expression have been proposed to play an important role. Even though it is accepted that differences in body size have cardiovascular effects, the use of scaling in the clinical setting is limited. An increase in the clinical utilization of scaling is thought to lead to improved cardiovascular disease diagnosis and management in patients.

Molecular basis and drug sensitivity of the delayed rectifier (IKr) in the fish heart

Fishes are increasingly used as models for human cardiac diseases, creating a need for a better understanding of the molecular basis of fish cardiac ion currents. To this end we cloned KCNH6 channel of the crucian carp (Carassius carassius) that produces the rapid component of the delayed rectifier K(+) current (IKr), the main repolarising current of the fish heart. KCNH6 (ccErg2) was the main isoform of the Kv11 potassium channel family with relative transcript levels of 98.9% and 99.6% in crucian carp atrium and ventricle, respectively. KCNH2 (ccErg1), an orthologue to human cardiac Erg (Herg) channel, was only slightly expressed in the crucian carp heart. The native atrial IKr and the cloned ccErg2 were inhibited by similar concentrations of verapamil, terfenadine and KB-R7943 (P>0.05), while the atrial IKr was about an order of magnitude more sensitive to E-4031 than ccErg2 (P<0.05) suggesting that some accessory β-subunits may be involved. Sensitivity of the crucian carp atrial IKr to E-4031, terfenadine and KB-R7943 was similar to what has been reported for the Herg channel. In contrast, the sensitivity of the crucian carp IKr to verapamil was approximately 30 times higher than the previously reported values for the Herg current. In conclusion, the cardiac IKr is produced by non-orthologous gene products in fish (Erg2) and mammalian hearts (Erg1) and some marked differences exist in drug sensitivity between fish and mammalian Erg1/2 which need to be taken into account when using fish heart as a model for human heart.

Regulatory evolution and voltage-gated ion channel expression in squid axon: selection-mutation balance and fitness cliffs

It has been suggested that optimization of either axonal conduction velocity or the energy efficiency of action potential conduction predominates in the selection of voltage-gated sodium conductance levels in the squid axon. A population genetics model of channel gene regulatory function was used to examine the role of these and other evolutionary forces on the selection of both sodium and potassium channel expression levels. In this model, the accumulating effects of mutations result in degradation of gene regulatory function, causing channel gene expression to fall to near-zero in the absence of positive selection. In the presence of positive selection, channel expression levels fall to the lowest values consistent with the selection criteria, thereby establishing a selection-mutation balance. Within the parameter space of sodium and potassium conductance values, the physiological performance of the squid axon model showed marked discontinuities associated with conduction failure and excitability. These discontinuities in physiological function may produce fitness cliffs. A fitness cliff associated with conduction failure, combined with the effects of phenotypic noise, can account for the selection of sodium conductance levels, without considering either conduction velocity or metabolic cost. A fitness cliff associated with a transition in axonal excitability, combined with phenotypic noise, can explain the selection of potassium channel expression levels. The results suggest that voltage-gated ion channel expression will fall to low levels, consistent with key functional constraints, even in the absence of positive selection for energy efficiency. Channel expression levels and individual variation in channel expression within the population can be explained by regulatory evolution in combination with genetic variation in regulatory function and phenotypic noise, without resorting to more complex mechanisms, such as activity-dependent homeostasis. Only a relatively small region of the large, nominally isofunctional parameter space for channel expression will normally be occupied, because of the effects of mutation.

Empirically founded genotype-phenotype maps from mammalian cyclic nucleotide-gated ion channels

A major barrier between evolutionary and functional biology is the difficulty of determining appropriate genotype-phenotype-fitness maps, particularly in metazoans. Concrete perspectives towards unifying these approaches are offered by studies on the physiological systems that depend on ion channel dynamics. I focus on the cyclic nucleotide-gated (CNG) channels implicated in the photoreceptor's response to light. From an evolutionary standpoint, sensory systems offers interpretative advantages, as the relation between the sensory response and environment is relatively straightforward. For CNG and other ion channels, extensive data are available about the physiological consequences of scanning mutagenesis on sensitive protein domains, such as the conduction pore. Mutant ion channels can be easily studied in living cells, so that the relation between genotypes and phenotypes is less speculative than usual. By relying on relatively simple theoretical frameworks, I used these data to relate the sequence space with phenotypes at increasing hierarchical levels. These empirical genotype-phenotype and phenotype-phenotype landscapes became smoother at higher integration levels, especially in heterozygous condition. The epistatic interaction between sites was analyzed from double mutant constructs. Magnitude epistasis was common. Moreover, evidence of reciprocal sign epistasis and the presence of permissive mutations were also observed, which suggest how adaptive regions can be connected across maladaptive valleys. The approach I describe suggests a way to better relate the evolutionary dynamics with the underlying physiology.

Electrical excitability of the heart in a Chondrostei fish, the Siberian sturgeon (Acipenser baerii)

Sturgeon (family Acipenseridae) are regarded as living fossils due to their ancient origin and exceptionally slow evolution. To extend our knowledge of fish cardiac excitability to a Chondrostei fish, we examined electrophysiological phenotype of the Siberian sturgeon ( Acipenser baerii) heart with recordings of epicardial ECG, intracellular action potentials (APs), and sarcolemmal ion currents. Epicardial ECG of A. baerii had the typical waveform of the vertebrate ECG with Q-T interval (average duration of ventricular AP) of 650 ± 30 ms and an intrinsic heart rate of 45.5 ± 5 beats min−1 at 20°C. Similar to other fish species, atrial AP was shorter in duration (402 ± 33 ms) than ventricular AP (585 ± 40) ( P < 0.05) at 20°C. Densities of atrial and ventricular Na+ currents were similar (−47.6 ± 4.5 and −53.2 ± 5.1 pA/pF, respectively) and close to the typical values of teleost hearts. Two major K+ currents, the inward rectifier K+ current ( IK1), and the delayed rectifier K+ current ( IKr) were found under basal conditions in sturgeon cardiomyocytes. The atrial IKr (3.3 ± 0.2 pA/pF) was about twice as large as the ventricular IKr (1.3 ± 0.4 pA/pF) ( P < 0.05) conforming to the typical pattern of teleost cardiac IKr. Divergent from other fishes, the ventricular IK1 was remarkably small (−2.5 ± 0.07 pA/pF) and not different from that of the atrial myocytes (−1.9 ± 0.06 pA/pF) ( P > 0.05). Two ligand-gated K+ currents were also found: ACh-activated inward rectifier ( IKACh) was present only in atrial cells, while ATP-sensitive K+ current ( IKATP) was activated by a mitochondrial blocker, CCCP, in both atrial and ventricular cells. The most striking difference to other fishes appeared in Ca2+ currents ( ICa). In atrial myocytes, ICa was predominated by nickel-sensitive and nifedipine-resistant T-type ICa, while ventricular myocytes had mainly nifedipine-sensitive and nickel-resistant L-type ICa. ICaT/ ICaL ratio of the sturgeon atrial myocytes (2.42) is the highest value ever measured for a vertebrate species. In ventricular myocytes, ICaT/ ICaL ratio was 0.09. With the exception of the large atrial ICaT and small ventricular IK1, electrical excitability of A. baerii heart is similar to that of teleost hearts.

Seasonal acclimatization of the cardiac action potential in the Arctic navaga cod (Eleginus navaga, Gadidae)

Freshwater fishes of north-temperate latitudes adjust electrical excitability of the heart to seasonal temperature changes by changing expression levels of ion channel isoforms. However, little is known about thermal responses of action potential (AP) in the hearts of marine polar fishes. To this end, we examined cardiac AP in the atrial myocardium of the Arctic navaga cod (Eleginus navaga) from the White Sea (Russia) acclimatized to winter (March) and summer (September) seasons. Acute increases in temperature from 4 to 10 °C were associated with increases in heart rate, maximum velocity of AP upstroke and negative resting membrane potential, while duration of AP was shortened in both winter-acclimatized and summer-acclimatized navaga hearts. In winter, there was a compensatory shortening (41.1%) of atrial AP duration and this was associated with a strong increase in transcript expression of Erg K(+) channels, known to produce the rapid component of the delayed rectifier K(+) current, I(Kr). Smaller increases were found in the expression of Kir2.1 channels that produce the inward rectifier K(+) current, I(K1). These findings indicate that the heart of navaga cod has a good acclimatory capacity in electrical excitation of cardiac myocytes, which enables cardiac function in the cold-eurythermal waters of the subarctic White Sea.

Evolution of CpG island promoter function underlies changes in KChIP2 potassium channel subunit gene expression in mammalian heart

Scaling of cardiac electrophysiology with body mass requires large changes in the ventricular action potential duration and heart rate in mammals. These changes in cellular electrophysiological function are produced by systematic and coordinated changes in the expression of multiple ion channel and transporter genes. Expression of one important potassium current, the transient outward current (I(to)), changes significantly during mammalian evolution. Changes in I(to) expression are determined, in part, by variation in the expression of an obligatory auxiliary subunit encoded by the KChIP2 gene. The KChIP2 gene is expressed in both cardiac myocytes and neurons and transcription in both cell types is initiated from the same CpG island promoter. Species-dependent variation of KChIP2 expression in heart is mediated by the evolution of the cis-regulatory function of this gene. Surprisingly, the major locus of evolutionary change for KChIP2 gene expression in heart lies within the CpG island core promoter. The results demonstrate that CpG island promoters are not simply permissive for gene expression but can also contribute to tissue-selective expression and, as such, can function as an important locus for the evolution of cis-regulatory function. More generally, evolution of the cis-regulatory function of voltage-gated ion channel genes appears to be an effective and efficient way to modify channel expression levels to optimize electrophysiological function.

Ancient origin of four-domain voltage-gated Na+ channels predates the divergence of animals and fungi

The four-domain voltage-gated Na(+) channels are believed to have arisen in multicellular animals, possibly during the evolution of the nervous system. Recent genomic studies reveal that many ion channels, including Na(+) channels and Ca(2+) channels previously thought to be restricted to animals, can be traced back to one of the unicellular ancestors of animals, Monosiga brevicollis. The eukaryotic supergroup Opisthokonta contains animals, fungi, and a diverse group of their unicellular relatives including M. brevicollis. Here, we demonstrate the presence of a putative voltage-gated Na(+) channel homolog (TtrNa(V)) in the apusozoan protist Thecamonas trahens, which belongs to the unicellular sister group to Opisthokonta. TtrNa(V) displays a unique selectivity motif distinct from most animal voltage-gated Na(+) channels. The identification of TtrNa(V) suggests that voltage-gated Na(+) channels might have evolved before the divergence of animals and fungi. Furthermore, our analyses reveal that Na(V) channels have been lost independently in the amoeboid holozoan Capsaspora owczarzaki of the animal lineage and in several basal fungi. These findings provide novel insights into the evolution of four-domain voltage-gated ion channels, ion selectivity, and membrane excitability in the Opisthokonta lineage.

Late phase of repolarization is autoregenerative and scales linearly with action potential duration in mammals ventricular myocytes: a model study

Scaling of action potential (AP) duration (APD) in mammals of different size is a rather complex phenomenon, dominated by a regulatory type mechanism of ion channels expression. By means of simulations performed on six published mathematical models of cardiac ventricular APs of different mammals, it is shown that AP repolarization is autoregenerative in its later phase (ARRP) and that the duration of such phase scales linearly with APD. For each AP, a 3-D instantaneous time-voltage-current surface is constructed, which has been recently described in a more simplified model. This representation allows us to measure ARRP and to study the contribution to it for different ion currents. It has been found that the existence of an ARRP is not intrinsic to cardiac models formulation; one out of the six models does not show this phase. A linear correlation between ARRP duration and APD in the remaining models is also found. It is shown that ARRP neither simply depend on AP shape nor on APD. Though I(K1) current seems to be the main responsible for determining and modulating this phase, the mechanism by which ARRP scales linearly with APD remains unclear and raises further questions on the scaling strategies of cardiac repolarization in mammals.

Degeneration of an intracellular ion channel in the primate lineage by relaxation of selective constraints

Ion channel genes are highly conserved and are rarely degenerated in the primate lineage leading to humans. So far, the only well-characterized ion channel known to be degenerated in primates is the plasma membrane transient receptor potential channel TRPC2, possibly due to changes in the pheromone signaling. Here, by analyzing the sequence data from ten primate species, we have determined the degeneration process of the TPC3 gene that encodes a member of the two-pore channel (TPC) family recently implicated in Ca(2+) release by nicotinic acid adenine dinucleotide phosphate from intracellular acidic stores in animals. We show that degeneration of TPC3 likely began in the common ancestors of Apes and Old World monkeys through a conserved inactivating mutation, followed by additional deleterious mutations resulting in the generation of a TPC3 pseudogene in the descendant catarrhine lineage. Located at a chromosome recombination hot spot, catarrhine TPC3 pseudogenes underwent a series of lineage-specific rearrangements, including exon deletion and duplication. In contrast, we identify near full-length TPC3 sequences in New World monkeys and Prosimians and show that the gene is subjected to strong purifying selection and therefore likely functional. Our data provide the first evidence for relaxed functional constraints for an intracellular ion channel in primates and shed novel insights into the evolution and regulation of Ca(2+) signaling in the primate lineage.