Gene Expression of the EGF System—a Prognostic Model in Non–Small Cell Lung Cancer Patients Without Activating EGFR Mutations

OBJECTIVES: Contradicting results have been demonstrated for the expression of the epidermal growth factor receptor (EGFR) as a prognostic marker in non–small cell lung cancer (NSCLC). The complexity of the EGF system with four interacting receptors and more than a dozen activating ligands is a likely explanation. The aim of this study is to demonstrate that the combined network of receptors and ligands from the EGF system is a prognostic marker. MATERIAL AND METHODS: Gene expression of the receptors EGFR, HER2, HER3, HER4, and the ligands AREG, HB-EGF, EPI, TGF-α, and EGF was measured by quantitative polymerase chain reaction in tumor samples from 100 NSCLC patients without EGFR activating mutations. Results were dichotomized into high or low levels of target expression. Coexpression of the ligands and receptors was observed, and a score was developed based on the summed effect of receptors and ligands. Akaike’s information criteria selected the optimal score. Results were correlated with age, sex, stage, histology, performance status, and overall survival. RESULTS: Patients were randomly split 1:1 to create test and validation cohorts. In multivariate analyses, the only individual prognostic marker was EPI (hazard ratio [HR] 0.38 [0.20-0.72], P = .003). The optimal score in the test cohort was validated as a marker of inferior survival in the validation cohort and by bootstrapping. Multivariate analysis confirmed the combined score as a prognostic marker of inferior survival (HR 3.75 [2.17-6.47], P < .00001). CONCLUSION: Our study has developed a model that takes the complexity of the EGF system into account and shows that this model is a strong prognostic marker in NSCLC patients.Despite advances in the treatment, non–small cell lung cancer (NSCLC) remains the leading cause of cancer-related death in the world [1]. In particular, the overall prognosis is poor for the metastatic stages, with a median overall survival (OS) of only 8 to 10 months. Even in the early nonmetastatic stages, the 5-year survival rate is as low as 50% [2], [3]. Prognostic markers are needed to stratify patients with different risk outcome. Several biomarkers have been evaluated in NSCLC, but only a few have proven to be clinically relevant. An activating mutation in the epidermal growth factor receptor (EGFR) is both a well-described predictive marker of benefit of EGFR-targeted tyrosine kinase inhibitors but also a debated prognostic marker of better OS [4], [5], [6], [7], [8], [9]. As EGFR expression has been associated with OS in head and neck, colorectal, and esophagus cancer [10], [11], [12], attention has been directed toward the use of EGFR expression as a prognostic marker in NSCLC, but contradicting results have been demonstrated [13], [14], [15], [16]. The EGF system is complex, and the effect of ligand-receptor interaction depends on a variety of different factors, which provides a plausible explanation for the divergence observed between studies that only evaluate EGFR expression in general. EGFR is one out of four related receptors from the EGF system and is capable of forming homodimers or heterodimers with one of the three other receptors when activated by a ligand. Several ligands from the EGF system such as amphiregulin (AREG), epidermal growth factor (EGF), and transforming growth factor–α (TGF-α) only activate EGFR, whereas some have the ability to activate several combinations of the four EGF receptors like heparin-binding epidermal growth factor (HB-EGF), epiregulin (EPI), and betacellulin (BCL). Most knowledge on the role of the ligands in NSCLC is from in vitro studies or from smaller clinical studies. In vitro studies have suggested that the biological effect of EGFR activation is dependent on the specific activating ligand as well as the dimerization partner [17]. Yet, no clinical studies have evaluated the effect of the network of receptors and ligands influencing EGFR in NSCLC. Furthermore, the majority of the clinical studies exploring EGFR expression are based on immunohistochemistry which is a semiquantitative method with a great risk of interobserver variability. Quantitative gene expression analyses provide a more accurate measure and are therefore more suitable for studies comparing expression levels. Prospectively, we have collected fresh tumor samples from patients suspected of lung cancer. Accordingly, the aim of this study is to evaluate the gene expression of the network of receptors and ligands of the EGF system affecting EGFR as a prognostic markers in NSCLC.     

Signal Intensities in Preoperative MRI Do Not Reflect Proliferative Activity in Meningioma

BACKGROUND: Identification of high-grade meningiomas in preoperative magnetic resonance imaging (MRI) is important for optimized surgical strategy and best possible resection. Numerous studies investigated subjectively determined morphological features as predictors of tumor biology in meningiomas. The aim of this study was to identify the predictive value of more reliable, quantitatively measured signal intensities in MRI for differentiation of high- and low-grade meningiomas and identification of meningiomas with high proliferation rates, respectively. PATIENTS AND METHODS: Sixty-six patients (56 World Health Organization [WHO] grade I, 9 WHO grade II, and 1 WHO grade I) were included in the study. Preoperative MRI signal intensities (fluid-attenuated inversion recovery [FLAIR], T1 precontrast, and T1 postcontrast as genuine and normalized values) were correlated with Ki-67 expression in tissue sections of resected meningiomas. Differences between the groups (analysis of variance) and Spearman rho correlation were computed using SPSS 22. RESULTS: Mean values of genuine signal intensities of meningiomas in FLAIR, T1 native, and T1 postcontrast were 323.9 ± 59, 332.8 ± 67.9, and 768.5 ± 165.3. Mean values of normalized (to the contralateral white matter) signal intensities of meningiomas in FLAIR, T1 native, and T1 postcontrast were 1.5 ± 0.3, 0.8 ± 0.1, and 1.9 ± 0.4. There was no significant correlation between MRI signal intensities and WHO grade or Ki-67 expression. Signal intensities did not differ significantly between WHO grade I and II/III meningiomas. Ki-67 expression was significantly increased in high-grade meningiomas compared with low-grade meningiomas (P < 0.01). Objectively measured values of MRI signal intensities (FLAIR, T1 precontrast, and T1 postcontrast enhancement) did not distinguish between high-grade and low-grade meningiomas. Furthermore, there was no association between MRI signal intensities and Ki-67 expression representing proliferative activity.Meningiomas are among the most common brain tumors. Their incidence is about 1%, and they account for almost one third of all primary intracranial masses. The majority of meningiomas are very slowly growing and nonsymptomatic or minimally symptomatic entities, discovered as incidental findings on neuroimaging [1]. The World Health Organization (WHO) classification system distinguishes 3 histological grades and 15 subtypes and is a well-accepted tool for prediction of prognosis. Although most meningiomas are benign masses, certain histological subtypes reveal very high recurrence rates despite the tumors’ seemingly total removal. Grade II (atypical) and grade III (anaplastic) meningiomas are associated with an increased risk of recurrence, are more aggressive, and show invasive behavior [2]. Grade I meningiomas are generally considered as benign tumors, but recent studies indicate substantial neurological deficits and impaired long-term survival due to tumor recurrence and stroke despite their low histopathological grading in a considerable proportion of cases [3], [4]. Increased mitotic activity (more than 4 mitoses per 10 high-power fields) and elevated Ki-67 expression (Ki-67 index of more than 5% of nuclei) are reliable histopathological markers for tumor recurrence [2].Because histopathological grading alone does not predict outcome satisfyingly, numerous studies investigated the value of preoperative magnetic resonance imaging (MRI) for prognostics. For example, Liu et al. demonstrated that hyperintensity on diffusion-weighted imaging, heterogeneous gadolinium enhancement, disruption of the arachnoid at brain tumor interface, T2 hyperintense peritumoral edema, and irregular tumor shape were independent predictors of non–grade I meningioma [5]. Other works produced comparable results, although some of these studies underline the importance of positive capsular enhancement [6], [7], whereas others emphasize the predictive value of peritumoral edema [5], [8]. All the above-cited works investigated morphological features of meningiomas summarized in subjective scoring systems, but not one of the studies objectively analyzed values of SIs in commonly used preoperative MRI sequences.Therefore, the aim of this study was to investigate the predictive value of genuine and normalized SIs of standardized preoperative MRI (T1 pre- and postcontrast, T2, and fluid-attenuated inversion recovery [FLAIR]) as in vivo predictors of proliferative activity of meningiomas.     

Systems Dynamic Modeling of a Guard Cell Cl? Channel Mutant Uncovers an Emergent Homeostatic Network Regulating Stomatal Transpiration

Stomata account for much of the 70% of global water usage associated with agriculture and have a profound impact on the water and carbon cycles of the world. Stomata have long been modeled mathematically, but until now, no systems analysis of a plant cell has yielded detail sufficient to guide phenotypic and mutational analysis. Here, we demonstrate the predictive power of a systems dynamic model in Arabidopsis (Arabidopsis thaliana) to explain the paradoxical suppression of channels that facilitate K⁺ uptake, slowing stomatal opening, by mutation of the SLAC1 anion channel, which mediates solute loss for closure. The model showed how anion accumulation in the mutant suppressed the H⁺ load on the cytosol and promoted Ca²⁺ influx to elevate cytosolic pH (pH_i) and free cytosolic Ca²⁺ concentration ([Ca2+]_i), in turn regulating the K⁺ channels. We have confirmed these predictions, measuring pH_i and [Ca2+]_i in vivo, and report that experimental manipulation of pH_i and [Ca2+]_i is sufficient to recover K⁺ channel activities and accelerate stomatal opening in the slac1 mutant. Thus, we uncover a previously unrecognized signaling network that ameliorates the effects of the slac1 mutant on transpiration by regulating the K⁺ channels. Additionally, these findings underscore the importance of H⁺-coupled anion transport for pH_i homeostasis.Guard cells surround stomatal pores in the epidermis of plant leaves and regulate pore aperture to balance the demands for CO₂ in photosynthesis with the need to conserve water by the plant. Transpiration through stomata accounts for much of the 70% of global water usage associated with agriculture, and it has a profound impact on the water and carbon cycles of the world (Gedney et al., 2006; Betts et al., 2007). Guard cells open the pore by transport and accumulation of osmotically active solutes, mainly K⁺ and Cl⁻ and the organic anion malate²⁻ (Mal), to drive water uptake and cell expansion. They close the pore by coordinating the release of these solutes through K⁺ and anion channels at the plasma membrane. The past half-century has generated a wealth of knowledge on guard cell transport, signaling, and homeostasis, resolving the properties of the major transport processes and metabolic pathways for osmotic solute uptake and accumulation, and many of the signaling pathways that control them (Blatt, 2000; Schroeder et al., 2001; McAinsh and Pittman, 2009; Hills et al., 2012). Even so, much of stomatal dynamics remains unresolved, especially how the entire network of transporters in guard cells works to modulate solute flux and how this network is integrated with organic acid metabolism (Wang and Blatt, 2011) to achieve a dynamic range of stomatal apertures.This gap in understanding is most evident in a number of often unexpected observations, many of which have led necessarily to ad hoc interpretations. Among these, recent studies highlighted a diurnal variation in the free cytosolic Ca²⁺ concentration ([Ca2+]_i), high in the daytime despite the activation of primary ion-exporting ATPases, and have been interpreted to require complex levels of regulation (Dodd et al., 2007). Other findings wholly defy intuitive explanation. For example, the tpk1 mutant of Arabidopsis (Arabidopsis thaliana) removes a major pathway for K⁺ flux across the tonoplast and suppresses stomatal closure, yet the mutant has no significant effect on cellular K⁺ content (Gobert et al., 2007). Similarly, the Arabidopsis clcc mutant eliminates the H⁺-Cl⁻ antiporter at the tonoplast; it affects Cl⁻ uptake, reduces vacuolar Cl⁻ content, and slows stomatal opening; however, counterintuitively, it also suppresses stomatal closure (Jossier et al., 2010). In work leading to this study, we observed that the slac1 anion channel mutant of Arabidopsis paradoxically profoundly alters the activities of the two predominant K⁺ channels at the guard cell plasma membrane. The SLAC1 anion channel is a major pathway for anion loss from the guard cells during stomatal closure (Negi et al., 2008; Vahisalu et al., 2008), and its mutation leads to incomplete and slowed closure of stomata in response to physiologically relevant signals of dark, high CO₂, and the water-stress hormone abscisic acid. Guard cells of the slac1 mutant accumulate substantially higher levels of Cl⁻, Mal, and also K⁺ when compared with guard cells of wild-type Arabidopsis (Negi et al., 2008). The latter observation is consistent with additional impacts on K⁺ transport; however, a straightforward explanation for these findings has not been not forthcoming.Quantitative systems analysis offers one approach to such problems. Efforts to model stomatal function generally have been driven by a “top-down” approach (Farquhar and Wong, 1984; Eamus and Shanahan, 2002) and have not incorporated detail essential to understanding the molecular and cellular mechanics that drive stomatal movement. Only recently we elaborated a quantitative systems dynamic approach to modeling the stomatal guard cell that incorporates all of the fundamental properties of the transporters at the plasma membrane and tonoplast, the salient features of osmolite metabolism, and the essential cytosolic pH (pH_i) and [Ca2+]_i buffering characteristics that have been described in the literature (Hills et al., 2012). The model resolved with this approach (Chen et al., 2012b) successfully recapitulated a wide range of known stomatal behaviors, including transport and aperture dependencies on extracellular pH, KCl, and CaCl₂ concentrations, diurnal changes in [Ca2+]_i (Dodd et al., 2007), and oscillations in membrane voltage and [Ca2+]_i thought to facilitate stomatal closure (Blatt, 2000; McAinsh and Pittman, 2009; Chen et al., 2012b). We have used this approach to resolve the mechanism behind the counterintuitive alterations in K⁺ channel activity uncovered in the slac1 mutant of Arabidopsis. Here, we show how anion accumulation in the mutant affects the H⁺ and Ca²⁺ loads on the cytosol, elevating pH_i and [Ca2+]_i, and in turn regulating the K⁺ channels. We have validated the key predictions of the model and, in so doing, have uncovered a previously unrecognized homeostatic network that ameliorates the effects of the slac1 mutant on transpiration from the plant.     

Focus Issue on Calcium Signaling: Analyses of a Gravistimulation-Specific Ca2+ Signature in Arabidopsis using Parabolic Flights

Gravity is a critical environmental factor affecting the morphology and functions of organisms on the Earth. Plants sense changes in the gravity vector (gravistimulation) and regulate their growth direction accordingly. In Arabidopsis (Arabidopsis thaliana) seedlings, gravistimulation, achieved by rotating the specimens under the ambient 1g of the Earth, is known to induce a biphasic (transient and sustained) increase in cytoplasmic calcium concentration ([Ca2+]_c). However, the [Ca2+]_c increase genuinely caused by gravistimulation has not been identified because gravistimulation is generally accompanied by rotation of specimens on the ground (1g), adding an additional mechanical signal to the treatment. Here, we demonstrate a gravistimulation-specific Ca²⁺ response in Arabidopsis seedlings by separating rotation from gravistimulation by using the microgravity (less than 10⁻⁴g) conditions provided by parabolic flights. Gravistimulation without rotating the specimen caused a sustained [Ca2+]_c increase, which corresponds closely to the second sustained [Ca2+]_c increase observed in ground experiments. The [Ca2+]_c increases were analyzed under a variety of gravity intensities (e.g. 0.5g, 1.5g, or 2g) combined with rapid switching between hypergravity and microgravity, demonstrating that Arabidopsis seedlings possess a very rapid gravity-sensing mechanism linearly transducing a wide range of gravitational changes (0.5g–2g) into Ca²⁺ signals on a subsecond time scale.Calcium ion (Ca²⁺) functions as an intracellular second messenger in many signaling pathways in plants (White and Broadley, 2003; Hetherington and Brownlee, 2004; McAinsh and Pittman, 2009; Spalding and Harper, 2011). Endogenous and exogenous signals are spatiotemporally encoded by changing the free cytoplasmic concentration of Ca²⁺ ([Ca2+]_c), which in turn triggers [Ca2+]_c-dependent downstream signaling (Sanders et al., 2002; Dodd et al., 2010). A variety of [Ca2+]_c increases induced by diverse environmental and developmental stimuli are reported, such as phytohormones (Allen et al., 2000), temperature (Plieth et al., 1999; Dodd et al., 2006), and touch (Knight et al., 1991; Monshausen et al., 2009). The [Ca2+]_c increase couples each stimulus and appropriate physiological responses. In the Ca²⁺ signaling pathways, the stimulus-specific [Ca2+]_c pattern (e.g. amplitude and oscillation) provide the critical information for cellular signaling (Scrase-Field and Knight, 2003; Dodd et al., 2010). Therefore, identification of the stimulus-specific [Ca2+]_c signature is crucial for an understanding of the intracellular signaling pathways and physiological responses triggered by each stimulus, as shown in the case of cold acclimation (Knight et al., 1996; Knight and Knight, 2000).Plants often exhibit biphasic [Ca2+]_c increases in response to environmental stimuli. Thus, slow cooling causes a fast [Ca2+]_c transient followed by a second, extended [Ca2+]_c increase in Arabidopsis (Arabidopsis thaliana; Plieth et al., 1999; Knight and Knight, 2000). The Ca²⁺ channel blocker lanthanum (La³⁺) attenuated the fast transient but not the following increase (Knight and Knight, 2000), suggesting that these two [Ca2+]_c peaks have different origins. Similarly, hypoosmotic shock caused a biphasic [Ca2+]_c increase in tobacco (Nicotiana tabacum) suspension-culture cells (Takahashi et al., 1997; Cessna et al., 1998). The first [Ca2+]_c peak was inhibited by gadolinium (Gd³⁺), La³⁺, and the Ca²⁺ chelator EGTA (Takahashi et al., 1997; Cessna et al., 1998), whereas the second [Ca2+]_c increase was inhibited by the intracellular Ca²⁺ store-depleting agent caffeine but not by EGTA (Cessna et al., 1998). The amplitude of the first [Ca2+]_c peak affected the amplitude of the second increase and vice versa (Cessna et al., 1998). These results suggest that even though the two [Ca2+]_c peaks originate from different Ca²⁺ fluxes (e.g. Ca²⁺ influx through the plasma membrane and Ca²⁺ release from subcellular stores, respectively), they are closely interrelated, showing the importance of the kinetic and pharmacological analyses of these [Ca2+]_c increases.Changes in the gravity vector (gravistimulation) could work as crucial environmental stimuli in plants and are generally achieved by rotating the specimens (e.g. +180°) in ground experiments. Use of Arabidopsis seedlings expressing apoaequorin, a Ca²⁺-reporting photoprotein (Plieth and Trewavas, 2002; Toyota et al., 2008a), has revealed that gravistimulation induces a biphasic [Ca2+]_c increase that may be involved in the sensory pathway for gravity perception/response (Pickard, 2007; Toyota and Gilroy, 2013) and the intracellular distribution of auxin transporters (Benjamins et al., 2003; Zhang et al., 2011). These two Ca²⁺ changes have different characteristics. The first transient [Ca2+]_c increase depends on the rotational velocity but not angle, whereas the second sustained [Ca2+]_c increase depends on the rotational angle but not velocity. The first [Ca2+]_c transient was inhibited by Gd³⁺, La³⁺, and the Ca²⁺ chelator 1,2-bis(2-aminophenoxy)ethane-N,N,N′,N′-tetraacetic acid but not by ruthenium red (RR), whereas the second sustained [Ca2+]_c increase was inhibited by all these chemicals. These results suggest that the first transient and second sustained [Ca2+]_c increases are related to the rotational stimulation and the gravistimulation, respectively, and are mediated by distinct molecular mechanisms (Toyota et al., 2008a). However, it has not been demonstrated directly that the second sustained [Ca2+]_c increase is induced solely by gravistimulation; it could be influenced by other factors, such as an interaction with the first transient [Ca2+]_c increase (Cessna et al., 1998), vibration, and/or deformation of plants during the rotation.To elucidate the genuine Ca²⁺ signature in response to gravistimulation in plants, we separated rotation and gravistimulation under microgravity (μg; less than 10⁻⁴g) conditions provided by parabolic flight (PF). Using this approach, we were able to apply rotation and gravistimulation to plants separately (Fig. 1). When Arabidopsis seedlings were rotated +180° under μg conditions, the [Ca2+]_c response to the rotation was transient and almost totally attenuated in a few seconds. Gravistimulation (transition from μg to 1.5g) was then applied to these prerotated specimens at the terminating phase of the PF. This gravistimulation without simultaneous rotation induced a sustained [Ca2+]_c increase. The kinetic properties of this sustained [Ca2+]_c increase were examined under different gravity intensities (0.5g–2g) and sequences of gravity intensity changes (Fig. 2A). This analysis revealed that gravistimulation-specific Ca²⁺ response has an almost linear dependency on gravitational acceleration (0.5g–2g) and an extremely rapid responsiveness of less than 1 s.Open in a separate windowFigure 1.Diagram of the experimental procedures for applying separately rotation and gravistimulation to Arabidopsis seedlings. Rotatory stimulation (green arrow) was applied by rotating the seedlings 180° under μg conditions, and 1.5g 180° rotation gravistimulation (blue arrow) was applied to the prerotated seedlings after μg.Open in a separate windowFigure 2.Acceleration, temperature, humidity, and pressure in an aircraft during flight experiments. A, Accelerations along x, y, and z axes in the aircraft during PF. The direction of flight (FWD) and coordinates (x, y, and z) are indicated in the bottom graph. The inset shows an enlargement of the acceleration along the z axis (gravitational acceleration) during μg conditions lasting for approximately 20 s. B, Temperature, humidity, and pressure in the aircraft during PF. Shaded areas in graphs denote the μg condition.     

An Optimal Frequency in Ca2+ Oscillations for Stomatal Closure Is an Emergent Property of Ion Transport in Guard Cells

Oscillations in cytosolic-free Ca²⁺ concentration ([Ca2+]_i) have been proposed to encode information that controls stomatal closure. [Ca2+]_i oscillations with a period near 10 min were previously shown to be optimal for stomatal closure in Arabidopsis (Arabidopsis thaliana), but the studies offered no insight into their origins or mechanisms of encoding to validate a role in signaling. We have used a proven systems modeling platform to investigate these [Ca2+]_i oscillations and analyze their origins in guard cell homeostasis and membrane transport. The model faithfully reproduced differences in stomatal closure as a function of oscillation frequency with an optimum period near 10 min under standard conditions. Analysis showed that this optimum was one of a range of frequencies that accelerated closure, each arising from a balance of transport and the prevailing ion gradients across the plasma membrane and tonoplast. These interactions emerge from the experimentally derived kinetics encoded in the model for each of the relevant transporters, without the need of any additional signaling component. The resulting frequencies are of sufficient duration to permit substantial changes in [Ca2+]_i and, with the accompanying oscillations in voltage, drive the K⁺ and anion efflux for stomatal closure. Thus, the frequency optima arise from emergent interactions of transport across the membrane system of the guard cell. Rather than encoding information for ion flux, these oscillations are a by-product of the transport activities that determine stomatal aperture.Stomata in the leaf epidermis are the main pathway both for CO₂ entry for photosynthesis and for foliar water loss by transpiration. Guard cells surround the stomatal pore and regulate the aperture, balancing the often conflicting demands for CO₂ and water conservation. Guard cells open and close the pore by expanding and contracting through the uptake and loss, respectively, of osmotic solutes, notably of K⁺, Cl⁻, and malate²⁻ (Mal²⁻; Pandey et al., 2007; Kim et al., 2010; Roelfsema and Hedrich, 2010; Lawson and Blatt, 2014). These transport processes comprise the final effectors of a regulatory network that coordinates transport across the plasma membrane and tonoplast, and maintains the homeostasis of the guard cell. A number of well-defined signals—including light, CO₂, drought and the water stress hormone abscisic acid (ABA)—act on this network, altering transport, solute content, turgor and cell volume, and ultimately stomatal aperture.Much research has focused on stomatal closure, underscoring both Ca²⁺-independent and Ca²⁺-dependent signaling. Of the latter, elevated cytosolic-free Ca²⁺ concentration ([Ca2+]_i) inactivates inward-rectifying K⁺ channels (I_K,in) to prevent K⁺ uptake and activates Cl⁻ (anion) channels (I_Cl) at the plasma membrane to depolarize the membrane and engage K⁺ efflux through outward-rectifying K⁺ channels (I_K,out; Keller et al., 1989; Blatt et al., 1990; Thiel et al., 1992; Lemtiri-Chlieh and MacRobbie, 1994). ABA, and most likely CO₂ (Kim et al., 2010), elevate [Ca2+]_i by facilitating Ca²⁺ entry at the plasma membrane to trigger Ca²⁺ release from endomembrane stores, a process often described as Ca²⁺-induced Ca²⁺ release (Grabov and Blatt, 1998, 1999). The hormone promotes Ca²⁺ influx by activating Ca²⁺ channels (I_Ca) at the plasma membrane, even in isolated membrane patches (Hamilton et al., 2000, 2001), which is linked to reactive oxygen species (Kwak et al., 2003; Wang et al., 2013). In parallel, cADP-ribose and nitric oxide promote endomembrane Ca²⁺ release and [Ca2+]_i elevation (Leckie et al., 1998; Neill et al., 2002; Garcia-Mata et al., 2003; Blatt et al., 2007). Best estimates indicate that endomembrane release accounts for more than 95% of the Ca²⁺ entering the cytosol to raise [Ca2+]_i (Chen et al., 2012; Wang et al., 2012).One feature of stomatal response to ABA, and indeed to a range of stimuli both hormonal as well as external, is its capacity for oscillations both in membrane voltage and [Ca2+]_i. Guard cell [Ca2+]_i at rest is typically around 100 to 200 nm, as it is in virtually all living cells. In response to ABA, [Ca2+]_i can rise above 1 μm—and locally, most likely above 10 μm—often in cyclic transients of tens of seconds to several minutes’ duration in association with oscillations in voltage and stomatal closure (Gradmann et al., 1993; McAinsh et al., 1995; Webb et al., 1996; Grabov and Blatt, 1998, 1999; Staxen et al., 1999; Allen et al., 2001). In principle, cycling in voltage and [Ca2+]_i arises as closure is accelerated with a controlled release of K⁺, Cl⁻, and Mal²⁻ from the guard cell and is subject to extracellular ion concentrations (Gradmann et al., 1993; Chen et al., 2012). However, it has been proposed that these, and similar oscillations in a variety of plant cell models, serve as physiological signals in their own right (McAinsh et al., 1995; Ehrhardt et al., 1996; Taylor et al., 1996). In support of such a signaling role, experiments designed to impose [Ca2+]_i (and voltage) oscillations in guard cells have yielded an optimal frequency for closure with a period near 10 min (Allen et al., 2001). Nonetheless, the studies offer no mechanistic explanation for this optimum that could validate a causal role in signaling, and none has been forthcoming since. Here we address questions of how such optimal frequencies in [Ca2+]_i oscillation arise and their relevance for stomatal closure, using quantitative systems analysis of guard cell transport and homeostasis. Our findings indicate that oscillations in voltage and [Ca2+]_i, and their optima associated with stomatal closure, are most simply explained as emerging from the interactions between ion transporters that drive stomatal closure. Thus, we conclude that these oscillations do not control, but are a by-product of the transport that determines stomatal aperture.     

Age Structure,Changing Demography and Effective Population Size in Atlantic Salmon (Salmo salar)

Effective population size (N_e) is a central evolutionary concept, but its genetic estimation can be significantly complicated by age structure. Here we investigate N_e in Atlantic salmon (Salmo salar) populations that have undergone changes in demography and population dynamics, applying four different genetic estimators. For this purpose we use genetic data (14 microsatellite markers) from archived scale samples collected between 1951 and 2004. Through life table simulations we assess the genetic consequences of life history variation on N_e. Although variation in reproductive contribution by mature parr affects age structure, we find that its effect on N_e estimation may be relatively minor. A comparison of estimator models suggests that even low iteroparity may upwardly bias N_e estimates when ignored (semelparity assumed) and should thus empirically be accounted for. Our results indicate that N_e may have changed over time in relatively small populations, but otherwise remained stable. Our ability to detect changes in N_e in larger populations was, however, likely hindered by sampling limitations. An evaluation of N_e estimates in a demographic context suggests that life history diversity, density-dependent factors, and metapopulation dynamics may all affect the genetic stability of these populations.THE effective size of a population (N_e) is an evolutionary parameter that can be informative on the strength of stochastic evolutionary processes, the relevance of which relative to deterministic forces has been debated for decades (e.g., Lande 1988). Stochastic forces include environmental, demographic, and genetic components, the latter two of which are thought to be more prominent at reduced population size, with potentially detrimental consequences for average individual fitness and population persistence (Newman and Pilson 1997; Saccheri et al. 1998; Frankham 2005). The quantification of N_e in conservation programs is thus frequently advocated (e.g., Luikart and Cornuet 1998; Schwartz et al. 2007), although gene flow deserves equal consideration given its countering effects on genetic stochasticity (Frankham et al. 2003; Palstra and Ruzzante 2008).Effective population size is determined mainly by the lifetime reproductive success of individuals in a population (Wright 1938; Felsenstein 1971). Variance in reproductive success, sex ratio, and population size fluctuations can reduce N_e below census population size (Frankham 1995). Given the difficulty in directly estimating N_e through quantification of these demographic factors (reviewed by Caballero 1994), efforts have been directed at inferring N_e indirectly through measurement of its genetic consequences (see Leberg 2005, Wang 2005, and Palstra and Ruzzante 2008 for reviews). Studies employing this approach have quantified historical levels of genetic diversity and genetic threats to population persistence (e.g., Nielsen et al. 1999b; Miller and Waits 2003; Johnson et al. 2004). N_e has been extensively studied in (commercially important) fish species, due to the common availability of collections of archived samples that facilitate genetic estimation using the temporal method (e.g., Hauser et al. 2002; Shrimpton and Heath 2003; Gomez-Uchida and Banks 2006; Saillant and Gold 2006).Most models relating N_e to a population''s genetic behavior make simplifying assumptions regarding population dynamics. Chiefly among these is the assumption of discrete generations, frequently violated in practice given that most natural populations are age structured with overlapping generations. Here, theoretical predictions still apply, provided that population size and age structure are constant (Felsenstein 1971; Hill 1972). Ignored age structure can introduce bias into temporal genetic methods for the estimation of N_e, especially for samples separated by time spans that are short relative to generation interval (Jorde and Ryman 1995; Waples and Yokota 2007; Palstra and Ruzzante 2008). Moreover, estimation methods that do account for age structure (e.g., Jorde and Ryman 1995) still assume this structure to be constant. Population dynamics will, however, likely be altered as population size changes, thus making precise quantifications of the genetic consequences of acute population declines difficult (Nunney 1993; Engen et al. 2005; Waples and Yokota 2007). This problem may be particularly relevant when declines are driven by anthropogenic impacts, such as selective harvesting regimes, that can affect age structure and N_e simultaneously (Ryman et al. 1981; Allendorf et al. 2008). Demographic changes thus have broad conservation implications, as they can affect a population''s sensitivity to environmental stochasticity and years of poor recruitment (Warner and Chesson 1985; Ellner and Hairston 1994; Gaggiotti and Vetter 1999). Consequently, although there is an urgent need to elucidate the genetic consequences of population declines, relatively little is understood about the behavior of N_e when population dynamics change (but see Engen et al. 2005, 2007).Here we focus on age structure and N_e in Atlantic salmon (Salmo salar) river populations in Newfoundland and Labrador. The freshwater habitat in this part of the species'' distribution range is relatively pristine (Parrish et al. 1998), yet Atlantic salmon in this area have experienced demographic declines, associated with a commercial marine fishery, characterized by high exploitation rates (40–80% of anadromous runs; Dempson et al. 2001). A fishery moratorium was declared in 1992, with rivers displaying differential recovery patterns since then (Dempson et al. 2004b), suggesting a geographically variable impact of deterministic and stochastic factors, possibly including genetics. An evaluation of those genetic consequences thus requires accounting for potential changes in population dynamics as well as in life history. Life history in Atlantic salmon can be highly versatile (Fleming 1996; Hutchings and Jones 1998; Fleming and Reynolds 2004), as exemplified by the high variation in age-at-maturity displayed among and within populations (Hutchings and Jones 1998), partly reflecting high phenotypic plasticity (Hutchings 2004). This diversity is particularly evident in the reproductive biology of males, which can mature as parr during juvenile freshwater stages (Jones and King 1952; Fleming and Reynolds 2004) and/or at various ages as anadromous individuals, when returning to spawn in freshwater from ocean migration. Variability in life history strategies is further augmented by iteroparity, which can be viewed as a bet-hedging strategy to deal with environmental uncertainty (e.g., Orzack and Tuljapurkar 1989; Fleming and Reynolds 2004). Life history diversity and plasticity may allow salmonid fish populations to alter and optimize their life history under changing demography and population dynamics, potentially acting to stabilize N_e. Reduced variance in individual reproductive success at low breeder abundance (genetic compensation) will achieve similar effects and might be a realistic aspect of salmonid breeding systems (Ardren and Kapuscinski 2003; Fraser et al. 2007b). Little is currently known about the relationships between life history plasticity, demographic change and N_e, partly due to scarcity of the multivariate data required for these analyses.Our objective in this article is twofold. First, we use demographic data for rivers in Newfoundland to quantify how life history variation influences age structure in Atlantic salmon and hence N_e and its empirical estimation from genetic data. We find that variation in reproductive contribution by mature parr has a much smaller effect on the estimation of N_e than is often assumed. Second, we use temporal genetic data to estimate N_e and quantify the genetic consequences of demographic changes. We attempt to account for potential sources of bias, associated with (changes in) age structure and life history, by using four different analytical models to estimate N_e: a single-sample estimator using the linkage disequilibrium method (Hill 1981), the temporal model assuming discrete generations (Nei and Tajima 1981; Waples 1989), and two temporal models for species with overlapping generations (Waples 1990a,b; Jorde and Ryman 1995) that differ principally in assumptions regarding iteroparity. A comparison of results from these different estimators suggests that iteroparity may often warrant analytical consideration, even when it is presumably low. Although sometimes limited by statistical power, a quantification and comparison of temporal changes in N_e among river populations suggests a more prominent impact of demographic changes on N_e in relatively small river populations.     

Identification of a Novel Class of Membrane-Bound [NiFe]-Hydrogenases in Thermococcus onnurineus NA1 by In Silico Analysis

In silico analysis of group 4 [NiFe]-hydrogenases from a hyperthermophilic archaeon, Thermococcus onnurineus NA1, revealed a novel tripartite gene cluster consisting of dehydrogenase-hydrogenase-cation/proton antiporter subunits, which may be classified as the new subgroup 4b of [NiFe]-hydrogenases-based on sequence motifs.Hydrogenases are the key enzymes involved in the metabolism of H₂, catalyzing the following chemical reaction: 2H⁺ + 2e⁻ ↔ H₂. Hydrogenases can be classified into [NiFe]-hydrogenases, [FeFe]-hydrogenases, and [Fe]-hydrogenases, based on their distinctive functional core containing the catalytic metal center (11, 17).The genomic analysis of Thermococcus onnurineus NA1, a hyperthermophilic archaeon isolated from a deep-sea hydrothermal vent area, revealed the presence of several distinct gene clusters encoding seven [NiFe]-hydrogenases and one homolog similar to Mbx (membrane-bound oxidoreductase) from Pyrococcus furiosus (1, 6, 8, 12). According to the classification system of hydrogenases by Vignais et al. (17), three hydrogenases (one F₄₂₀-reducing and two NADP-reducing hydrogenases) belong to group 3 [NiFe]-hydrogenases, and four hydrogenases belong to group 4 [NiFe]-hydrogenases. The group 4 hydrogenases are widely distributed among bacteria and archaea (17), with Hyc and Hyf (hydrogenase 3 and 4, respectively) from Escherichia coli (19), Coo (CO-induced hydrogenase) from Rhodospirillum rubrum (4), Ech (energy-converting hydrogenase) from Methanosarcina barkeri (7), and Mbh (membrane-bound hydrogenase) from P. furiosus (6, 10, 12) being relatively well-characterized hydrogenases in this group. One of the four group 4 hydrogenases from T. onnurineus NA1 was found to be similar in sequence to that of P. furiosus Mbh (10).     

Focus Issue on Calcium Signaling: Calcium-Dependent Protein Kinase CPK6 Positively Functions in Induction by Yeast Elicitor of Stomatal Closure and Inhibition by Yeast Elicitor of Light-Induced Stomatal Opening in Arabidopsis

Yeast elicitor (YEL) induces stomatal closure that is mediated by a Ca²⁺-dependent signaling pathway. A Ca²⁺-dependent protein kinase, CPK6, positively regulates activation of ion channels in abscisic acid and methyl jasmonate signaling, leading to stomatal closure in Arabidopsis (Arabidopsis thaliana). YEL also inhibits light-induced stomatal opening. However, it remains unknown whether CPK6 is involved in induction by YEL of stomatal closure or in inhibition by YEL of light-induced stomatal opening. In this study, we investigated the roles of CPK6 in induction by YEL of stomatal closure and inhibition by YEL of light-induced stomatal opening in Arabidopsis. Disruption of CPK6 gene impaired induction by YEL of stomatal closure and inhibition by YEL of light-induced stomatal opening. Activation by YEL of nonselective Ca²⁺-permeable cation channels was impaired in cpk6-2 guard cells, and transient elevations elicited by YEL in cytosolic-free Ca²⁺ concentration were suppressed in cpk6-2 and cpk6-1 guard cells. YEL activated slow anion channels in wild-type guard cells but not in cpk6-2 or cpk6-1 and inhibited inward-rectifying K⁺ channels in wild-type guard cells but not in cpk6-2 or cpk6-1. The cpk6-2 and cpk6-1 mutations inhibited YEL-induced hydrogen peroxide accumulation in guard cells and apoplast of rosette leaves but did not affect YEL-induced hydrogen peroxide production in the apoplast of rosette leaves. These results suggest that CPK6 positively functions in induction by YEL of stomatal closure and inhibition by YEL of light-induced stomatal opening in Arabidopsis and is a convergent point of signaling pathways for stomatal closure in response to abiotic and biotic stress.Stomata, formed by pairs of guard cells, play a critical role in regulation of plant CO₂ uptake and water loss, thus critically influencing plant growth and water stress responsiveness. Guard cells respond to a variety of abiotic and biotic stimuli, such as light, drought, and pathogen attack (Israelsson et al., 2006; Shimazaki et al., 2007; Melotto et al., 2008).Elicitors derived from microbial surface mimic pathogen attack and induce stomatal closure in various plant species such as Solanum lycopersicum (Lee et al., 1999), Commelina communis (Lee et al., 1999), Hordeum vulgare (Koers et al., 2011), and Arabidopsis (Arabidopsis thaliana; Melotto et al., 2006; Khokon et al., 2010). Yeast elicitor (YEL) induces stomatal closure in Arabidopsis (Klüsener et al., 2002; Khokon et al., 2010; Salam et al., 2013). Our recent studies showed that YEL inhibits light-induced stomatal opening and that protein phosphorylation is involved in induction by YEL of stomatal closure and inhibition by YEL of light-induced stomatal opening (Salam et al., 2013).Cytosolic Ca²⁺ has long been recognized as a conserved second messenger in stomatal movement (Shimazaki et al., 2007; Roelfsema and Hedrich 2010; Hubbard et al., 2012). Elevation of cytosolic free Ca²⁺ concentration ([Ca2+]_cyt) is triggered by influx of Ca²⁺ from apoplast and release of Ca²⁺ from intracellular stores in guard cell signaling (Leckie et al., 1998; Hamilton et al., 2000; Pei et al., 2000; Garcia-Mata et al., 2003; Lemtiri-Chlieh et al., 2003). The influx of Ca²⁺ is carried by nonselective Ca²⁺-permeable cation (I_Ca) channels that are activated by plasma membrane hyperpolarization and H₂O₂ (Pei et al., 2000; Murata et al., 2001; Kwak et al., 2003). Elevation of [Ca2+]_cyt activates slow anion (S-type) channels and down-regulates inward-rectifying potassium (K_in) channels in guard cells (Schroeder and Hagiwara, 1989; Grabov and Blatt, 1999). The activation of S-type channels is a hallmark of stomatal closure, and the suppression of K_in channels is favorable to stomatal closure but not to stomatal opening (Pei et al., 1997; Kwak et al., 2001; Xue et al., 2011; Uraji et al., 2012).YEL induces stomatal closure with extracellular H₂O₂ production, intracellular H₂O₂ accumulation, activation of I_Ca channels, and transient [Ca2+]_cyt elevations (Klüsener et al., 2002; Khokon et al., 2010). However, it remains to be clarified whether YEL activates S-type channels and inhibits K_in channels in guard cells.Calcium-dependent protein kinases (CDPKs) are regulators in Ca²⁺-dependent guard cell signaling (Mori et al., 2006; Zhu et al., 2007; Geiger et al., 2010, 2011; Zou et al., 2010; Munemasa et al., 2011; Brandt et al., 2012; Scherzer et al., 2012). In guard cells, CDPKs regulate activation of S-type and I_Ca channels and inhibition of K_in channels (Mori et al., 2006; Zou et al., 2010; Munemasa et al., 2011). A CDPK, CPK6, positively regulates activation of S-type channels and I_Ca channels without affecting H₂O₂ production in abscisic acid (ABA)- and methyl jasmonate (MeJA)-induced stomatal closure (Mori et al., 2006; Munemasa et al., 2011). CPK6 phosphorylates and activates SLOW ANION CHANNEL-ASSOCIATED1 expressed in Xenopus spp. oocyte (Brandt et al., 2012; Scherzer et al., 2012). These findings underline the role of CPK6 in regulation of ion channel activation and stomatal movement, leading us to test whether CPK6 regulates the induction by YEL of stomatal closure and inhibition by YEL of light-induced stomatal opening.In this study, we investigated activation of S-type channels and inhibition of K_in channels by YEL and roles of CPK6 in induction by YEL of stomatal closure and inhibition by YEL of light-induced stomatal opening. For this purpose, we examined the effects of mutation of CPK6 on induction by YEL of stomatal closure and inhibition by YEL of light-induced stomatal opening, activation of I_Ca channels, transient [Ca2+]_cyt elevations, activation of S-type channels, inhibition of K_in channels, H₂O₂ production in leaves, and H₂O₂ accumulation in leaves and guard cells.     

An Increasing Danger of Zoonotic Orthopoxvirus Infections

On May 8, 1980, the World Health Assembly at its 33^rd session solemnly declared that the world and all its peoples had won freedom from smallpox and recommended ceasing the vaccination of the population against smallpox. Currently, a larger part of the world population has no immunity not only against smallpox but also against other zoonotic orthopoxvirus infections. Recently, recorded outbreaks of orthopoxvirus diseases not only of domestic animals but also of humans have become more frequent. All this indicates a new situation in the ecology and evolution of zoonotic orthopoxviruses. Analysis of state-of-the-art data on the phylogenetic relationships, ecology, and host range of orthopoxviruses—etiological agents of smallpox (variola virus, VARV), monkeypox (MPXV), cowpox (CPXV), vaccinia (VACV), and camelpox (CMLV)—as well as the patterns of their evolution suggests that a VARV-like virus could emerge in the course of natural evolution of modern zoonotic orthopoxviruses. Thus, there is an insistent need for organization of the international control over the outbreaks of zoonotic orthopoxvirus infections in various countries to provide a rapid response and prevent them from developing into epidemics.The genus Orthopoxvirus of the family Poxviridae comprises the species variola (smallpox) virus (VARV), with human as its only sensitive host; zoonotic species monkeypox virus (MPXV), cowpox virus (CPXV), vaccinia virus (VACV), and camelpox virus (CMLV); and several others. These orthopoxviruses are immunologically cross-reactive and cross-protective, so that infection with any member of this genus provides protection against infection with any other member of the genus [1], [2]. Traditionally, the species of the Orthopoxvirus genus have been named primarily according to the host animal from which they were isolated and identified based on a range of biological characteristics [1]. Most frequently, zoonotic orthopoxviruses have been initially isolated from animals immediately close to humans being incidental hosts for the virus, the natural carriers of which are, as a rule, wild animals. Correspondingly, the name of an orthopoxvirus species does not reflect the actual animal that is its natural reservoir.With accumulation of the data on complete genome nucleotide sequences for various strains of orthopoxvirus species, it has been found that an interesting feature of the orthopoxvirus genomes is the presence of genes that are intact in one species but fragmented or deleted in another [3]–[8]. These data confirm the concept of a reductive evolution of orthopoxviruses, according to which the gene loss plays an important role in the evolutionary adaptation of progenitor virus to a particular environmental niche (host) and emergence of new virus species [9]. CPXV has the largest genome of all the modern representatives of the genus Orthopoxvirus, and this genome contains all the genes found in the other species of this genus [2], [4], [10]–[12]. Therefore, Cowpox virus was proposed as the closest of all the modern species to the progenitor virus for the genus Orthopoxvirus, while the remaining species, Variola virus included, had appeared as a result of multistage reductive evolution [4], [9], [13].VARV, the most pathogenic species for humans, has the smallest genome of all the orthopoxviruses [2]–[7]. This suggests a potential possibility for emergence of a VARV-like variant from the currently existing zoonotic orthopoxviruses with longer genomes in the course of natural evolution. It is known that although mutational changes are rather a rare event for the poxvirus DNA [13], characteristic of these viruses is the possibility of intermolecular and intramolecular recombinations, as well as genomic insertions and deletions [14], [15]. It has been recently found that duplication/amplification of genomic segments is typical of poxviruses, and in the case of a certain selective pressure (for example, host antiviral defenses), certain genes are able to relatively rapidly accumulate mutations that would provide the virus adaptation to new conditions, including a new host [16].The conducted analysis of the available archive data on smallpox and the history of ancient civilizations as well as the newest data on the evolutionary relationships of orthopoxviruses has allowed me to suggest the hypothesis that smallpox could have repeatedly emerged in the past via evolutionary changes of a zoonotic progenitor virus [17].Because of the cessation of the vaccination against smallpox after its eradication 35 years ago, a tremendous part of the world human population currently has no immunity not only against smallpox, but also against any other zoonotic orthopoxvirus infections. This new situation allows orthopoxviruses to circulate in the human population and, as a consequence, should alter several established concepts on the ecology and range of sensitive hosts for various orthopoxvirus species.The most intricate case is the origin of VACV. For many decades, VACV has been used for vaccinating humans against smallpox, and it was considered that this virus, variolae vaccinae, originates from zoonotic CPXV, introduced to immunization practice by Jenner as early as 1796 [1]. Only in the 20^th century was it found out that the orthopoxvirus strains used for smallpox vaccination significantly differ in their properties from both the natural CPXV isolates recovered from cows and the other orthopoxvirus species examined by that time [18]. Correspondingly, they were regarded as a separate species, Vaccinia virus [19]. Moreover, it was inferred that the VACV natural reservoir was unknown and numerous hypotheses attempted to explain the origin of this virus while passaging progenitor viruses in animals in the process of vaccine production [1], [2], [20].The issue of VACV origin was somewhat clarified after sequencing the complete genome of horsepox virus (HSPV) [21], which appeared to be closely related to the sequenced VACV strains. Only after this was attention paid to the fact that Jenner specified the origin of his vaccine from an infection of the heels of horses (“grease”) and indicated that the vaccine became more suitable for human use after passage through the cow [20]. This suggests that VACV may originate from a zoonotic HSPV, which naturally persisted concurrently with CPXV. Some facts suggest that the infectious materials not only from cow lesions but also from horse lesions were used for smallpox vaccination in the 19^th century. The vaccine lymph from the horse gave the most satisfactory results in inducing an anti-smallpox immunity as well as less side reactions [1]. By all accounts, they gradually commenced using HSPV isolates for smallpox vaccination, the future generations of which recovered decades later were ascribed to the separate species Vaccinia virus [19], rather than CPXV for smallpox vaccination everywhere.Since the 1960s, VACVs have been repeatedly isolated in Brazil [22]. The first VACV isolates were recovered from wild rodents (sentinel mice and rice rat) [23]. Since 1999, an ever-increasing number of exanthematous outbreaks affecting dairy cows and their handlers have been recorded [24]–[27], supplemented recently with outbreaks among horses [28], [29]. Several VACV strains have been isolated during these outbreaks from cows, horses, humans, and rodents [22], [27], [28], [30], [31]. The questions that arise are when and how VACV entered Brazil and the wild nature of the American continent. The more widespread point of view is that VACV strains could be transmitted from vaccinated humans to domestic animals and further to wild ones with subsequent adaptation to the rural environment [22]. My standpoint implies that HSPV/VACV could have been repeatedly accidentally imported from Europe to South America with the infected horses or rodents to be further introduced into wildlife. Possibly, the latter hypothesis more adequately reflects the actual pathway of VACV transmission to the Brazilian environment, since recent phylogenetic studies have suggested an independent origin for South American VACV isolates, distinct from the vaccine strains used on this continent during the WHO smallpox eradication campaign [22], [32]. Presumably, genome-wide sequencing of the viruses will give a more precise answer to the origin of VACV variants isolated in Brazil.In the past, the outbreaks of buffalopox had occurred frequently in various states of India as well as in Pakistan, Bangladesh, Indonesia, Egypt, and other countries [33]. The causative agent, buffalopox virus (BPXV), is closely related to VACV and affiliated with the species Vaccinia virus, genus Orthopoxvirus [2], [34]. Recently, mass outbreaks of buffalopox in domestic buffaloes along with severe zoonotic infection in milk attendants were recorded at various places in India [35], [36]. In several buffalopox outbreaks, the BPXV-caused infections were recorded in cows in the same herds [37]. An increase in BPXV transmission to different species, including buffaloes, cows, and humans, suggests the reemergence of zoonotic buffalopox infection [35], [38]. The buffalopox outbreaks recorded in different distant regions of India are likely to suggest the presence of an abundant natural BPXV reservoir represented by wild animals, most probably rodents. Correspondingly, it is of the paramount importance to perform a large-scale study of the presence of orthopoxviruses in wild animals of India.Thus, yet incomplete data on the modern ecology of VACV and BPXV allow for speculation that the orthopoxviruses belonging to the species Vaccinia virus have a wide host range, are zoonotic, are currently spread over large areas in Eurasia and South America, and that their natural carriers are several rodents.CPXV has relatively low pathogenicity for humans but has a wide range of sensitive animal hosts [2], [39]. Human cowpox is a rare sporadic disease, which develops when CPXV is transmitted from an infected animal to human [2], [40]. This disease is mainly recorded in Europe. In wildlife, CPXV carriers are asymptomatically infected rodents [41], [42]. During the last two decades, reports on an increasing number of CPXV infections in cats, rats, exotic animals, and humans have been published [43]–[47]. Comparative studies of the properties of CPXV isolates recovered from various hosts at different times and in several geographic zones have shown sufficient intraspecific variations [2], [48], [49]. A recent phylogenetic analysis of the complete genomes of 12 CPXV strains recovered from humans and several animal species suggests that they be split into two major Cowpox virus–like and Vaccinia virus–like clades [50]. This means that the criteria of the separation of orthopoxviruses into these two species should be corrected.MPXV is a zoonotic virus causing a human infection similar to smallpox in its clinical manifestations with a lethality rate of 1–8% [51]. The natural reservoir of MPXV is various species of African rodents [8], [10]. The active surveillance data in the same health zone (Democratic Republic of Congo) from the 1980s to 2006–2007 suggest a 20-fold increase in human monkeypox incidence 30 years after the cessation of the smallpox vaccination campaign [52]. This poses the question of whether MPXV can acquire the possibility of a high human-to-human transmission rate, characteristic of VARV, under conditions of a long-term absence of vaccination and considerably higher incidence of human infection. If this occurs, humankind will face a problem considerably more complex than with the smallpox eradication. First and foremost, this is determined by the fact that MPXV, unlike VARV, has its natural reservoir represented by numerous African rodents [2], [53].In its biological properties and according to the data of phylogenetic analysis of the complete virus genomic sequence, CMLV is closest to VARV, the causative agent of smallpox, as compared with the other orthopoxvirus species [1], [8]. Camelpox is recognized as one of the most important viral diseases in camels. This infection was first described in India in 1909. Subsequently, camelpox outbreaks have been reported in many countries of the Middle East, Asia, and Africa [54], [55]. Until recently, it has been commonly accepted that the host range of CMLV is confined to one animal species, camels [1], [55]. However, the first human cases of camelpox have been recently confirmed in India [56]. This suggests that camelpox could be a zoonotic disease. Since camelpox outbreaks occur irregularly in distant regions of the world and the viruses isolated during these outbreaks display different degrees of virulence [55], it is possible to postulate the presence of a wildlife animal reservoir of CMLV other than camels. Since the camelpox outbreaks are usually associated with the rainy season of the year, when rodents are actively reproducing, it is likely that rodents could be the natural carriers of CMLV.It is known that most of the emerging human pathogens originate from zoonotic pathogens [57]–[59]. Many viruses do not cause the disease in their natural reservoir hosts but can be highly pathogenic when transmitted to a new host species. Emerging and reemerging human pathogens more often are those with broad host ranges. The viruses able to infect many animal species are evolutionarily adapted to utilizing different cell mechanisms for their reproduction and, thus, can extend/change their host range with a higher probability [58].There are no fundamental prohibitions for the possible reemergence of smallpox or a similar human disease in the future as a result of natural evolution of the currently existing zoonotic orthopoxviruses. An ever-increasing sensitivity of the human population to zoonotic orthopoxviruses, resulting from cessation of the mass smallpox vaccination, elevates the probability for new variants of these viruses, potentially dangerous for humans, to emerge. However, the current situation is radically different from the ancient one, since many outbreaks of orthopoxvirus infections among domestic animals and humans are recorded and studied.Recently, the efforts of scientists under WHO control are directed to the development of state-of-the-art methods for VARV rapid identification as well as design of new generation safe smallpox vaccines and drugs against VARV and other orthopoxviruses [60]. The designed promising anti-orthopoxvirus drugs display no pronounced virus species specificity. Therefore, they are applicable in the outbreaks caused by any orthopoxvirus species. International acceptance of the designed highly efficient anti-orthopoxvirus drugs ST-246 and CMX001 [60] is of paramount importance.In the areas of high incidence of zoonotic orthopoxviral infections, it would be purposeful to vaccinate domestic and zoo animals as well as the persons closely associated with them using state-of-the-art safe vaccines based on VACV, which has a wide range of sensitive hosts. This would considerably decrease the likelihood for such infections to spread from wildlife into the human environment.In the African region endemic for monkeypox, which also displays a high rate of HIV infection, the population could be vaccinated with the VACV strain MVA, which has been recently demonstrated to be safe even for HIV-infected persons [61].Taking into account the above mentioned increased incidence of outbreaks of animal and human orthopoxvirus infections and their potential danger, it is important to accelerate organization of the international Smallpox Laboratory Network, discussed by the WHO Advisory Committee on Variola Virus Research [62], [63], and orient this network to express diagnosing not only of VARV but also of other zoonotic orthopoxviruses. This will provide constant monitoring of these infections in all parts of the world and make it possible to prevent the development of small outbreaks into expanded epidemics, thereby decreasing the risk of evolutional changes and emergence of an orthopoxvirus highly pathogenic for humans.The international system for clinical sampling and identification of infectious agents has been worked out and optimized while implementing the global smallpox eradication program under the aegis of the WHO as well as anti-epidemic measures and methods for mass vaccination [1]. The accumulated experience is of paramount importance for the establishment of international control not only over currently existing orthopoxvirus infections but also other emerging and reemerging diseases.     

Apamin Does Not Inhibit Human Cardiac Na+ Current,L-type Ca2+ Current or Other Major K+ Currents

Background

Apamin is commonly used as a small-conductance Ca²⁺-activated K⁺ (SK) current inhibitor. However, the specificity of apamin in cardiac tissues remains unclear.

Objective

To test the hypothesis that apamin does not inhibit any major cardiac ion currents.

Methods

We studied human embryonic kidney (HEK) 293 cells that expressed human voltage-gated Na⁺, K⁺ and Ca²⁺ currents and isolated rabbit ventricular myocytes. Whole-cell patch clamp techniques were used to determine ionic current densities before and after apamin administration.

Results

Ca²⁺ currents (CACNA1c+CACNB2b) were not affected by apamin (500 nM) (data are presented as median [25^th percentile;75^th percentile] (from –16 [–20;–10] to –17 [–19;–13] pA/pF, P = NS), but were reduced by nifedipine to –1.6 [–3.2;–1.3] pA/pF (p = 0.008). Na⁺ currents (SCN5A) were not affected by apamin (from –261 [–282;–145] to –268 [–379;–132] pA/pF, P = NS), but were reduced by flecainide to –57 [–70;–47] pA/pF (p = 0.018). None of the major K⁺ currents (I _Ks, I _Kr, I _K1 and I _to) were inhibited by 500 nM of apamin (KCNQ1+KCNE1, from 28 [20]; [37] to 23 [18]; [32] pA/pF; KCNH2+KCNE2, from 28 [24]; [30] to 27 [24]; [29] pA/pF; KCNJ2, from –46 [–48;–40] to –46 [–51;–35] pA/pF; KCND3, from 608 [505;748] to 606 [454;684]). Apamin did not inhibit the I _Na or I _CaL in isolated rabbit ventricular myocytes (I _Na, from –67 [–75;–59] to –68 [–71;–59] pA/pF; I _CaL, from –16 [–17;–14] to –14 [–15;–13] pA/pF, P = NS for both).

Conclusions

Apamin does not inhibit human cardiac Na⁺ currents, L-type Ca²⁺ currents or other major K⁺ currents. These findings indicate that apamin is a specific SK current inhibitor in hearts as well as in other organs.     

Self-Incompatibility-Induced Programmed Cell Death in Field Poppy Pollen Involves Dramatic Acidification of the Incompatible Pollen Tube Cytosol

Self-incompatibility (SI) is an important genetically controlled mechanism to prevent inbreeding in higher plants. SI involves highly specific interactions during pollination, resulting in the rejection of incompatible (self) pollen. Programmed cell death (PCD) is an important mechanism for destroying cells in a precisely regulated manner. SI in field poppy (Papaver rhoeas) triggers PCD in incompatible pollen. During SI-induced PCD, we previously observed a major acidification of the pollen cytosol. Here, we present measurements of temporal alterations in cytosolic pH ([pH]_cyt); they were surprisingly rapid, reaching pH 6.4 within 10 min of SI induction and stabilizing by 60 min at pH 5.5. By manipulating the [pH]_cyt of the pollen tubes in vivo, we show that [pH]_cyt acidification is an integral and essential event for SI-induced PCD. Here, we provide evidence showing the physiological relevance of the cytosolic acidification and identify key targets of this major physiological alteration. A small drop in [pH]_cyt inhibits the activity of a soluble inorganic pyrophosphatase required for pollen tube growth. We also show that [pH]_cyt acidification is necessary and sufficient for triggering several key hallmark features of the SI PCD signaling pathway, notably activation of a DEVDase/caspase-3-like activity and formation of SI-induced punctate actin foci. Importantly, the actin binding proteins Cyclase-Associated Protein and Actin-Depolymerizing Factor are identified as key downstream targets. Thus, we have shown the biological relevance of an extreme but physiologically relevant alteration in [pH]_cyt and its effect on several components in the context of SI-induced events and PCD.Programmed cell death (PCD) in plants is relatively well documented and characterized (Jones and Dangl, 1996; van Doorn, 2011; van Doorn et al., 2011). There is considerable biochemical evidence for the involvement of caspase-like activities in plant PCD (van Doorn and Woltering, 2005). For example, the vacuolar processing enzyme has YVADase (caspase-1-like) activity (Hatsugai et al., 2004; Rojo et al., 2004; Hara-Nishimura et al., 2005), DEVDase (caspase-3-like) and YVADases are associated with PCD in several plant systems (del Pozo and Lam, 1998; Korthout et al., 2000; Danon et al., 2004), and VEIDase (caspase-6-like) is the main caspase-like activity involved in embryonic pattern formation (Bozhkov et al., 2004). However, because plants have no caspase gene homologs (Sanmartín et al., 2005), the nature of their caspase-like enzymes is the subject of considerable debate. Vacuolar cell death is one of two major classes of PCD in plants (van Doorn et al., 2011). It is thought that collapse of the vacuole is a key irreversible step in several plant PCD systems, including during tissue and organ formation, such as the classic differentiation of tracheary elements (Hara-Nishimura and Hatsugai, 2011). Exactly how this is achieved and what processes are involved remain unknown. Until very recently, it was generally thought that the rupturing vacuole releases proteases, hydrolases, and nucleases, allowing cellular disassembly by an autophagy-like process. Some PCD systems cannot be assigned to either class; these include PCD triggered by the hypersensitive response to biotrophic pathogens, PCD in cereal endosperm, and self-incompatibility (SI)-induced PCD (van Doorn et al., 2011).SI is a genetically controlled pollen-pistil cell-cell recognition system. Self-pollen is recognized by the stigma as being genetically identical, resulting in inhibition of pollen tube growth. Most SI systems use tightly linked polymorphic genes: the pollen (male) and pistil (female) S-determinants. In field poppy (Papaver rhoeas), the S-determinants are a 14-kD signaling ligand field poppy stigma S (PrsS) and a unique transmembrane protein field poppy pollen S (PrpS; Foote et al., 1994; Wheeler et al., 2010). These interact in an S-specific manner, and increases in cytosolic free calcium ([Ca2+]_cyt) are triggered in incompatible pollen tubes (Franklin-Tong et al., 1993), resulting in phosphorylation of soluble inorganic pyrophosphatases (sPPases; Rudd et al., 1996; de Graaf et al., 2006), activation of a Mitogen-Activated Protein Kinase (MAPK; Rudd et al., 2003), and increases in reactive oxygen species (ROS) and nitric oxide (Wilkins et al., 2011, 2014). Most of these components are integrated into a signaling network leading to PCD (Bosch et al., 2008; Wilkins et al., 2014). The actin cytoskeleton is a key target in the field poppy SI response, undergoing depolymerization (Snowman et al., 2002) followed by polymerization into highly stable F-actin foci decorated with the actin binding proteins (ABPs) Actin-Depolymerizing Factor (ADF) and Cyclase-Associated Protein (CAP; Poulter et al., 2010, 2011), with both processes being involved in mediating PCD (Thomas et al., 2006). A major player in SI-mediated PCD is a caspase-3-like/DEVDase-like activity (Thomas and Franklin-Tong, 2004; Bosch and Franklin-Tong, 2007). The SI-induced caspase-3-like/DEVDase exhibits maximum substrate cleavage in vitro at pH 5, with peak activity 5 h after SI induction in vivo (Bosch and Franklin-Tong, 2007). The low pH optimum for this caspase-3-like/DEVDase activity is unusual, because most of the cytosolic plant caspase-like activities identified to date have optimal activity close to normal physiological pH (approximate pH, 6.5–7.0; Korthout et al., 2000; Bozhkov et al., 2004; Coffeen and Wolpert, 2004). Because the SI-induced cytosolic-located DEVDase requires a low pH for activity, this suggested that, during SI, the pollen tube cytosol undergoes dramatic acidification. In vivo pH measurements of the cytosol at 1 to 4 h after SI induction confirmed this, when cytosolic pH ([pH]_cyt) had dropped from pH 6.9 to pH 5.5 (Bosch and Franklin-Tong, 2007). This fits the in vitro pH optimum of the caspase-3-like/DEVDase almost exactly, implicating pollen cytosolic acidification as playing a vital role in creating optimal conditions for the activation of the caspase-3-like/DEVDase-like activity and progression of PCD.Under normal cellular conditions, [pH]_cyt is between approximately 6.9 and 7.5 (Kurkdjian and Guern, 1989; Felle, 2001). Pollen tubes, like other tip-growing cells, have [pH]_cyt gradients (Gibbon and Kropf, 1994; Feijó et al., 1999). The [pH]_cyt of the pollen tube shank is an approximate pH of 6.9 to 7.11 (Fricker et al., 1997; Messerli and Robinson, 1998). There has been much debate about the [pH]_cyt gradient, comprising an apical domain with an approximate pH of 6.8 and a subapical alkaline band with an approximate pH of 7.2 to 7.8 in Lilium longiflorum and Lilium formosanum pollen tubes (Fricker et al., 1997; Messerli and Robinson, 1998; Feijó et al., 2001; Lovy-Wheeler et al., 2006). Oscillations of [pH]_cyt between approximate pH values of 6.9 and 7.3 have been linked to tip growth in L. formosanum pollen tubes (Lovy-Wheeler et al., 2006). The vacuole and the apoplast have a highly acidic pH between pH 5 and pH 6 (Katsuhara et al., 1989; Feijó et al., 1999). The majority of studies of pH changes in plant cells reports modest, transient changes in [pH]_cyt of approximately 0.4 and 0.7 pH units during development, gravitropic responses, decreases in light intensity, and addition of elicitors, hormones, and other treatments. For example, during root hair development in Arabidopsis (Arabidopsis thaliana), root [pH]_cyt was elevated from an approximate pH of 7.3 to 7.7 (Bibikova et al., 1998). Root gravitropic responses stimulate small transient [pH]_cyt alterations (Scott and Allen, 1999; Fasano et al., 2001; Johannes et al., 2001). More recently, it has been shown that the [pH]_cyt drops during PCD controlling root cap development; however, exactly how many units the [pH]_cyt decreased was not measured (Fendrych et al., 2014). Other studies investigating [pH]_cyt in response to physiologically relevant signals also report small transient alterations. Light-adapted cells respond to a decrease in light intensity with a rapid transient cytosolic acidification by approximately 0.3 pH units (Felle et al., 1986). Addition of nodulation factors resulted in an increase of 0.2 pH units in root hairs (Felle et al., 1998), and abscisic acid increased the [pH]_cyt of guard cells by 0.3 pH units (Blatt and Armstrong, 1993). Changes in [pH]_cyt are thought to activate stress responses (Felle, 2001). Elicitor treatments resulted in a [pH]_cyt drop of between 0.4 and 0.7 pH units in suspension cells (Mathieu et al., 1996; Kuchitsu et al., 1997), a drop of 0.2 pH units in Nitellopsis obtusa cells treated with salt (Katsuhara et al., 1989), and a drop of 0.3 to 0.7 pH units in Eschscholzia californica (Roos et al., 1998).Here, we investigate SI-induced acidification of the cytosol, providing measurements of physiologically relevant temporal alterations in [pH]_cyt, and identify key targets of this, providing mechanistic insights into these events. The SI-induced acidification plays a pivotal role in the activation of a caspase-3-like/DEVDase activity, the formation of punctate F-actin foci, and ABP localization during SI PCD. We investigate the vacuole as a potential contributor to SI-induced [pH]_cyt acidification.     

Generation and Initial Characterization of FDD Knock In Mice

Background

Mutations in the integral membrane protein 2B [1], also known as BRI₂ [2], a type II trans-membrane domain protein cause two autosomal dominant neurodegenerative diseases, Familial British and Danish Dementia [3]. In these conditions, accumulation of a C-terminal peptide (ABri and ADan) cleaved off from the mutated precursor protein by the pro-protein convertase furin [4], leads to amyloid deposition in the walls of blood vessels and parenchyma of the brain. Recent advances in the understanding of the generation of amyloid in Alzheimer''s disease has lead to the finding that BRI₂ interacts with the Amyloid Precursor Protein (APP), decreasing the efficiency of APP processing to generate Aβ [5], [6], [7]. The interaction between the two precursors, APP and BRI₂, and possibly between Aβ and ABri or ADan, could be important in influencing the rate of amyloid production or the tendency of these peptides to aggregate.

Methodology/Principal Findings

We have generated the first BRI₂ Danish Knock-In (FDD_KI) murine model of FDD, expressing the pathogenic decamer duplication in exon 6 of the BRI₂ gene. FDD_KI mice do not show any evident abnormal phenotype, with normal brain histology and no detectable amyloid deposition in blood vessel walls or parenchyma.

Conclusions/Significance

This new murine mouse model will be important to further understand the interaction between APP and BRI₂, and to provide insights into the molecular basis of FDD.     

Heteromeric AtKC1��AKT1 Channels in Arabidopsis Roots Facilitate Growth under K+-limiting Conditions

Plant growth and development is driven by osmotic processes. Potassium represents the major osmotically active cation in plants cells. The uptake of this inorganic osmolyte from the soil in Arabidopsis involves a root K⁺ uptake module consisting of the two K⁺ channel α-subunits, AKT1 and AtKC1. AKT1-mediated potassium absorption from K⁺-depleted soil was shown to depend on the calcium-sensing proteins CBL1/9 and their interacting kinase CIPK23. Here we show that upon activation by the CBL·CIPK complex in low external potassium homomeric AKT1 channels open at voltages positive of E_K, a condition resulting in cellular K⁺ leakage. Although at submillimolar external potassium an intrinsic K⁺ sensor reduces AKT1 channel cord conductance, loss of cytosolic potassium is not completely abolished under these conditions. Depending on channel activity and the actual potassium gradients, this channel-mediated K⁺ loss results in impaired plant growth in the atkc1 mutant. Incorporation of the AtKC1 subunit into the channel complex, however, modulates the properties of the K⁺ uptake module to prevent K⁺ loss. Upon assembly of AKT1 and AtKC1, the activation threshold of the root inward rectifier voltage gate is shifted negative by approximately −70 mV. Additionally, the channel conductance gains a hypersensitive K⁺ dependence. Together, these two processes appear to represent a safety strategy preventing K⁺ loss through the uptake channels under physiological conditions. Similar growth retardation phenotypes of akt1 and atkc1 loss-of-function mutants in response to limiting K⁺ supply further support such functional interdependence of AKT1 and AtKC1. Taken together, these findings suggest an essential role of AtKC1-like subunits for root K⁺ uptake and K⁺ homeostasis when plants experience conditions of K⁺ limitation.Fundamental plant functions such as control of the membrane potential, osmo-regulation, and turgor-driven growth and movements are based on the availability to gain high cellular potassium concentrations (1). The absorption of this inorganic osmolyte from the soil by the root therefore represents a pivotal process for plant life. Classical experiments by Epstein et al. in 1963 (2) described K⁺ root uptake as a biphasic process mediated by two uptake mechanisms: high affinity potassium transport with apparent affinities of ∼20 μm and a low affinity transport system with K_m values in the millimolar range. During the last decades several molecular components of potassium transport systems have been identified and functionally characterized in plants (3, 4). Mutant analyses, heterologous expression, as well as radiotracer uptake experiments characterized the K⁺ channels AKT1·AtKC1 and members of the HAK·KT·KUP family as major components of the Arabidopsis thaliana root-localized potassium transport system (5–9). In this study we focused on AKT1 and AtKC1, members of the Arabidopsis Shaker-like K⁺ channel family. AKT1 is a voltage-dependent inward-rectifying K⁺ channel mediating potassium uptake over a wide range of external potassium concentrations (10–15). Root cells of the akt1-1 loss-of-function mutant completely lack inward rectifying K⁺ currents (12). As a consequence the growth of akt1-1 seedlings is strongly impaired on low potassium medium (100 μm and less) (11, 12, 15). Rescue of yeast growth on 20 μm K⁺ and patch clamp experiments (16, 17) directly demonstrated that plant inward rectifying K⁺ channels are capable of serving as high affinity potassium uptake transporters. AtKC1 shares its expression pattern with AKT1 (18–20). AtKC1 α-subunits, however, neither form functional channels in akt1-1 knock-out plants nor in heterologous expression systems. In contrast to root cells of akt1-1 loss of function mutants, root protoplasts of AtKC1 null mutants (atkc1-f) still exhibit inward rectifying potassium currents most likely derived from homomeric AKT1 tetramers (20). Inward K⁺ currents in this atkc1-f mutant were characterized by a more positive activation voltage. These data suggested that the AtKC1 α-subunits do not form K⁺ channels per se but modulate the properties of the AKT1·AtKC1 heterocomplex (20–22). Previously, two groups in their ground-breaking studies demonstrated that AKT1 is activated by the CBL²-interacting, serine/threonine kinase, CIPK23, particularly under low K⁺ conditions (23, 24). CIPK23 itself was shown to be activated by the two calcineurin B-like proteins, CBL1 and 9, acting in a Ca²⁺-dependent manner upstream of CIPK23 (25, 26). Genetic disruption of these elements resulted in transgenic plants exhibiting a phenotype comparable with that of the AKT1 loss of function mutant. This regulatory system, based on a calcium sensor, a protein kinase, and a K⁺ channel, was functionally reconstituted in Xenopus oocytes (23, 24, 27), suggesting that these elements are essential and sufficient to operate as a low K⁺-sensitive potassium uptake system. Here we report on the physiological properties of the heteromeric K⁺ uptake module formed by the predominant root potassium uptake channel subunits, AKT1 and AtKC1 and its regulating kinase complex, CBL1 and CIPK23. Our studies show that the physical interaction of the CBL1·CIPK23 complex is specific for AKT1 channels and does not involve the AtKC1 subunit. AKT1 possesses a K⁺ (absence) sensor affecting channel activity at submillimolar K⁺ concentrations by strongly reducing its maximal cord conductance. Despite this K⁺ sensor, upon activation, AKT1 homomeric channels were shown to represent a potassium leak at low external potassium concentrations. Integration of AtKC1 into the K⁺ uptake module, however, prevented potassium loss by modulating both the voltage sensor and conductance in the channel complex. Moreover, activation of the AKT1-like maize channel ZMK1 by CBL1·CIPK23 suggests a conserved interaction and regulation across monocot and dicotyledonous plant species. Our biophysical studies as well as growth assays with plant mutant lines lacking the respective channels underline that acquisition of potassium under limiting K⁺ conditions is mediated via the root AKT1·AtKC1 K⁺ uptake channel complex.     

PKCλ/ι signaling promotes triple-negative breast cancer growth and metastasis

Triple-negative breast cancer (TNBC) is a distinct breast cancer subtype defined by the absence of estrogen receptor (ER), progesterone receptor (PR) and epidermal growth factor receptor 2 (HER2/neu), and the patients with TNBC are often diagnosed with higher rates of recurrence and metastasis. Because of the absence of ER, PR and HER2/neu expressions, TNBC patients are insensitive to HER2-directed and endocrine therapies available for breast cancer treatment. Here, we report that expression of atypical protein kinase C isoform, PKCλ/ι, significantly increased and activated in all invasive breast cancer (invasive ductal carcinoma or IDC) subtypes including the TNBC subtype. Because of the lack of targeted therapies for TNBC, we choose to study PKCλ/ι signaling as a potential therapeutic target for TNBC. Our observations indicated that PKCλ/ι signaling is highly active during breast cancer invasive progression, and metastatic breast cancers, the advanced stages of breast cancer disease that developed more frequently in TNBC patients, are also characterized with high levels of PKCλ/ι expression and activation. Functional analysis in experimental mouse models revealed that depletion of PKCλ/ι significantly reduces TNBC growth as well as lung metastatic colonization. Furthermore, we have identified a PKCλ/ι-regulated gene signature consisting of 110 genes, which are significantly associated with indolent to invasive progression of human breast cancer and poor prognosis. Mechanistically, cytokines such as TGFβ and IL1β could activate PKCλ/ι signaling in TNBC cells and depletion of PKCλ/ι impairs NF-κB p65 (RelA) nuclear localization. We observed that cytokine-PKCλ/ι-RelA signaling axis, at least in part, involved in modulating gene expression to regulate invasion of TNBC cells. Overall, our results indicate that induction and activation of PKCλ/ι promote TNBC growth, invasion and metastasis. Thus, targeting PKCλ/ι signaling could be a therapeutic option for breast cancer, including the TNBC subtype.Breast cancer is a clinically heterogeneous disease and both intra and inter-tumor heterogeneities provide great challenges for developing successful therapies. Expressions (or absence thereof) of estrogen receptor (ER), progesterone receptor (PR) and epidermal growth factor receptor 2 (HER2)/neu are widely used to clinically classify breast tumors into multiple therapeutic groups.¹ The ER/PR-positive and the HER2-positive breast cancer patients could be benefited from endocrine and HER2-targeted therapies.¹ However, triple-negative breast cancers (TNBCs), which represent ∼12–17% of all breast cancer,² lack ER, PR and HER2/neu expressions² and are not responsive to therapies targeting these receptors. Therefore, the only systemic therapy available for TNBC is chemotherapy.³ Furthermore, TNBC is associated with aggressive pathologic features like higher histology grade and mitotic index⁴ and often found to be associated with higher rate of metastasis and recurrence leading to limited clinical outcome.^{5, 6, 7, 8} Recurrence of TNBC tends to recur within a few years after successful initial treatment^{6, 9} and often develops metastasis to the bone, brain and lungs with poor prognosis.^{2, 6} Thus, identification of signaling pathways that regulate malignant progression of breast cancer subtypes, especially TNBCs, would be therapeutically important.In recent years, PKC signaling has been implicated in modulating invasion and metastasis of multiple tumors.^{10, 11} The PKC family consists of multiple serine/threonine kinases and the relative contribution of individual PKC isoforms during cancer progression varies due to pleiotropism.¹² PKC isoforms regulate diverse cellular functions such as cell-cycle regulation, cellular survival, cell–cell communications and apoptosis.¹³ In particular, atypical PKC isoforms, PKCζ and atypical protein kinase C lamda/iota (PKCλ/ι), are known to be important for chemotaxis, cell polarity, migration and wound healing processes.^{14, 15} Aberrations in all these processes are manifested in tumor progression and metastasis.¹⁴ Consistent with these notions, recent studies indicated that atypical PKCs are associated with various human cancers.^{10, 11} Importantly, the PKCλ/ι gene is located at the 3q26.2 genomic region, which is most frequently amplified in human cancer^{16, 17}, and overexpression of PKCλ/ι has been implicated in cancer development in multiple tissues including the lung,^{18, 19} pancreas,²⁰ stomach,²¹ colon,²² esophagus,²³ liver,²⁴ bile duct,²⁵ ovary,¹⁷ prostate²⁶ and brain.²⁷ Recently, few studies have been reported higher expression of PKCλ/ι in ER/PR- and HER-positive breast cancer and also in lymph node metastases.^{28, 29} Kojima et. al.²⁸ showed that PKCλ/ι expression is highly induced in the ER/PR- and HER2-positive IDCs compared with ductal carcinoma in situ (DCIS) and normal breast. PKCλ/ι forms apical-junctional complexes (AJCs) with other polarity proteins such as partitioning defective 3 homolog (PAR3) and partitioning defective 6 homolog (PAR6),^{30, 31, 32, 33} and invasiveness of breast tumor cells was shown to be associated with loss of PKCλ/ι localization from their apical domains.²⁸ In addition, predominant nuclear localization of PKCλ/ι in both normal and atypical ductal hyperplasia (ADH) lesions prompted the concept that PKCλ/ι might be in an inactive state in these lesions.²⁸ However, expression and activation of PKCλ/ι in TNBCs and the functional importance of PKCλ/ι signaling in relation to invasive breast cancer progression and metastasis are very poorly understood.^{10, 11}Here, we studied PKCλ/ι signaling during invasive progression of TNBC. We utilized expression evaluations in triple-negative IDCs as well as metastatic breast cancers of human patients, in vitro and in vivo functional assays, and global gene expression analysis of human patient samples. We concluded that PKCλ/ι signaling is an important regulator for invasion and metastatic progression of human breast cancers including triple-negative subtypes.