Feeds:
Posts
Comments

Archive for the ‘Evolutionary cognition’ Category

Therapeutic Implications for Targeted Therapy from the Resurgence of Warburg ‘Hypothesis’

Writer and Curator: Larry H. Bernstein, MD, FCAP 

(Note that each portion of the discussion is followed by a reference)

It is now a time to pause after almost a century of a biological scientific discoveries that have transformed the practice of medicine and impacted the lives of several generations of young minds determined to probe the limits of our knowledge.  In the century that we have entered into the scientific framework of medicine has brought together a difficult to grasp evolution of the emergence of human existence from wars, famine, droughts, storms, infectious diseases, and insect born pestilence with betterment of human lives, only unevenly divided among societal classes that have existed since time immemorial. In this short time span there have emerged several generations of physicians who have benefited from a far better medical education that their forebears could have known. In this expansive volume on cancer, we follow an incomplete and continuing challenge to understand cancer, a disease that has become associated with longer life spans in developed nations.

While there are significant improvements in the diagnosis and treatment of cancers, there is still a personal as well as locality factor in the occurrence of this group of diseases, which has been viewed incorrectly as a “dedifferentiation” of mature tissue types and the emergence of a cell phenotype that is dependent on glucose, reverts to a cancer “stem cell type” (loss of stemness), loses cell to cell adhesion, loses orderly maturation, and metastasizes to distant sites. At the same time, physician and nurses are stressed in the care of patients by balancing their daily lives and maintaining a perspective.

The conceptual challenge of cancer diagnosis and management has seemed insurmountable, but owes much to the post World War I activities of Otto Heinrich Warburg. It was Warburg who made the observation that cancer cells metabolize glucose by fermentation in much the way Pasteur 60 years earlier observed fermentation of yeast cells. This metabolic phenomenon occurs even in the presence of an oxygen supply, which would provide a huge deficit in ATP production compared with respiration. The cancer cell is “addicted to glucose” and produced lactic acid. Warburg was awarded the Nobel Prize in Medicine for this work in 1931.

In the last 15 years there has been a resurgence of work on the Warburg effect that sheds much new light on the process that was not previously possible, with significant therapeutic implications.  In the first place, the metabolic mechanism for the Warburg effect was incomplete even at the beginning of the 21st century.  This has been partly rectified with the enlightening elucidation of genome modifications, cellular metabolic regulation, and signaling pathways.

The following developments have become central to furthering our understanding of malignant transformation.

  1. There is usually an identifiable risk factor, such as, H. pylori, or of a chronic inflammatory state, as in the case of Barrett’s esophagus.
  2. There are certain changes in glucose metabolism that have been unquestionably been found in the evolution of this disease. The changes are associated with major changes in metabolic pathways, miRN signaling, and the metabolism geared to synthesis of cells with an impairment of the cell death cycle. In these changes, mitochondrial function is central to both the impaired respiration and the autophagy geared to the synthesis of cancer cells.

The emergence of this cell prototype is characterized by the following, again related to the Warburg effect:

  1. Cancer cells oxidize a decreased fraction of the pyruvate generated from glycolysis
  2. The mitochondrial pyruvate carrier (MPC), composed of the products of the MPC1 and MPC2 genes, modulates fractional pyruvate oxidation. MPC1 is deleted or underexpressed in multiple cancers and correlates with poor prognosis.
  3. Cancer cells tend to express a partially inhibited splice variant of pyruvate kinase (PK-M2), leading to decreased pyruvate production.
  4. The two proteins that mediate pyruvate conversion to lactate and its export, M-type lactate dehydrogenase and the monocarboxylate transporter MCT-4, are commonly upregulated in cancer cells leading to decreased pyruvate oxidation.
  5. The enzymatic step following mitochondrial entry is the conversion of pyruvate to acetyl-CoA by the pyruvate dehydrogenase (PDH) complex. Cancer cells frequently exhibit increased expression of the PDH kinase PDK1, which phosphorylates and inactivates PDH. This PDH regulatory mechanism is required for oncogene induced transformation and reversed in oncogene-induced senescence.
  6. The PDK inhibitor dichloroacetate has shown some clinical efficacy, which correlates with increased pyruvate oxidation. One of the simplest mechanisms to explain decreased mitochondrial pyruvate oxidation in cancer cells, a loss of mitochondrial pyruvate import, has been observed repeatedly over the past 40 years. This process has been impossible to study at a molecular level until recently, however, as the identities of the protein(s) that mediate mitochondrial pyruvate uptake were unknown.
  7. The mitochondrial pyruvate carrier (MPC) as a multimeric complex that is necessary for efficient mitochondrial pyruvate uptake. The MPC contains two distinct proteins, MPC1 and MPC2; the absence of either leads to a loss of mitochondrial pyruvate uptake and utilization in yeast, flies, and mammalian cells.

A Role for the Mitochondrial Pyruvate Carrier as a Repressor of the Warburg Effect and Colon Cancer Cell Growth

John C. Schell, Kristofor A. Olson, Lei Jiang, Amy J. Hawkins, et al.
Molecular Cell Nov 6, 2014; 56: 400–413.
http://dx.doi.org/10.1016/j.molcel.2014.09.026

In addition to the above, the following study has therapeutic importance:

Glycolysis has become a target of anticancer strategies. Glucose deprivation is sufficient to induce growth inhibition and cell death in cancer cells. The increased glucose transport in cancer cells has been attributed primarily to the upregulation of glucose transporter 1 (Glut1),  1 of the more than 10 glucose transporters that are responsible for basal glucose transport in almost all cell types. Glut1 has not been targeted until very recently due to the lack of potent and selective inhibitors.

First, Glut1 antibodies were shown to inhibit cancer cell growth. Other Glut1 inhibitors and glucose transport inhibitors, such as fasentin and phloretin, were also shown to be effective in reducing cancer cell growth. A group of inhibitors of glucose transporters has been recently identified with IC50 values lower than 20mmol/L for inhibiting cancer cell growth. However, no animal or detailed mechanism studies have been reported with these inhibitors.

Recently, a small molecule named STF-31 was identified that selectively targets the von Hippel-Lindau (VHL) deficient kidney cancer cells. STF-31 inhibits VHL deficient cancer cells by inhibiting Glut1. It was further shown that daily intraperitoneal injection of a soluble analogue of STF-31 effectively reduced the growth of tumors of VHL-deficient cancer cells grafted on nude mice. On the other hand, STF-31 appears to be an inhibitor with a narrow cell target spectrum.

These investigators recently reported the identification of a group of novel small compounds that inhibit basal glucose transport and reduce cancer cell growth by a glucose deprivation–like mechanism. These compounds target Glut1 and are efficacious in vivo as anticancer agents. A novel representative compound WZB117 not only inhibited cell growth in cancer cell lines but also inhibited cancer growth in a nude mouse model. Daily intraperitoneal injection of WZB117 resulted in a more than 70% reduction in the size of human lung cancer of A549 cell origin. Mechanism studies showed that WZB117 inhibited glucose transport in human red blood cells (RBC), which express Glut1 as their sole glucose transporter. Cancer cell treatment with WZB117 led to decreases in levels of Glut1 protein, intracellular ATP, and glycolytic enzymes. All these changes were followed by increase in ATP sensing enzyme AMP-activated protein kinase (AMPK) and declines in cyclin E2 as well as phosphorylated retinoblastoma, resulting in cell-cycle arrest, senescence, and necrosis. Addition of extracellular ATP rescued compound-treated cancer cells, suggesting that the reduction of intracellular ATP plays an important role in the anticancer mechanism of the molecule.

A Small-Molecule Inhibitor of Glucose Transporter 1 Downregulates Glycolysis, Induces Cell-Cycle Arrest, and Inhibits Cancer Cell Growth In Vitro and In Vivo

Yi Liu, Yanyan Cao, Weihe Zhang, Stephen Bergmeier, et al.
Mol Cancer Ther Aug 2012; 11(8): 1672–82
http://dx.doi.org://10.1158/1535-7163.MCT-12-0131

Alterations in cellular metabolism are among the most consistent hallmarks of cancer. These investigators have studied the relationship between increased aerobic lactate production and mitochondrial physiology in tumor cells. To diminish the ability of malignant cells to metabolize pyruvate to lactate, M-type lactate dehydrogenase levels were knocked down by means of LDH-A short hairpin RNAs. Reduction in LDH-A activity resulted in stimulation of mitochondrial respiration and decrease of mitochondrial membrane potential. It also compromised the ability of these tumor cells to proliferate under hypoxia. The tumorigenicity of the LDH-A-deficient cells was severely diminished, and this phenotype was reversed by complementation with the human ortholog LDH-A protein. These results demonstrate that LDH-A plays a key role in tumor maintenance.

The results are consistent with a functional connection between alterations in glucose metabolism and mitochondrial physiology in cancer. The data also reflect that the dependency of tumor cells on glucose metabolism is a liability for these cells under limited-oxygen conditions. Interfering with LDH-A activity as a means of blocking pyruvate to lactate conversion could be exploited therapeutically. Because individuals with complete deficiency of LDH-A do not show any symptoms under ordinary circumstances, the genetic data suggest that inhibition of LDH-A activity may represent a relatively nontoxic approach to interfere with tumor growth.

Attenuation of LDH-A expression uncovers a link between glycolysis, mitochondrial physiology, and tumor maintenance

Valeria R. Fantin Julie St-Pierre and Philip Leder
Cancer Cell Jun 2006; 9: 425–434.
http://dx.doi.org:/10.1016/j.ccr.2006.04.02

The widespread clinical use of positron-emission tomography (PET) for the detection of aerobic glycolysis in tumors and recent findings have rekindled interest in Warburg’s theory. Studies on the physiological changes in malignant conversion provided a metabolic signature for the different stages of tumorigenesis; during tumorigenesis, an increase in glucose uptake and lactate production have been detected. The fully transformed state is most dependent on aerobic glycolysis and least dependent on the mitochondrial machinery for ATP synthesis.

Tumors ferment glucose to lactate even in the presence of oxygen (aerobic glycolysis; Warburg effect). The pentose phosphate pathway (PPP) allows glucose conversion to ribose for nucleic acid synthesis and glucose degradation to lactate. The nonoxidative part of the PPP is controlled by transketolase enzyme reactions. We have detected upregulation of a mutated transketolase transcript (TKTL1) in human malignancies, whereas transketolase (TKT) and transketolase-like-2 (TKTL2) transcripts were not upregulated. Strong TKTL1 protein expression was correlated to invasive colon and urothelial tumors and to poor patients outcome. TKTL1 encodes a transketolase with unusual enzymatic properties, which are likely to be caused by the internal deletion of conserved residues. We propose that TKTL1 upregulation in tumors leads to enhanced, oxygen-independent glucose usage and a lactate based matrix degradation. As inhibition of transketolase enzyme reactions suppresses tumor growth and metastasis, TKTL1 could be the relevant target for novel anti-transketolase cancer therapies. We suggest an individualized cancer therapy based on the determination of metabolic changes in tumors that might enable the targeted inhibition of invasion and metastasis.

Other important links between cancer-causing genes and glucose metabolism have been already identified. Activation of the oncogenic kinase Akt has been shown to stimulate glucose uptake and metabolism in cancer cells and renders these cells susceptible to death in response to glucose withdrawal. Such tumor cells have been shown to be dependent on glucose because the ability to induce fatty acid oxidation in response to glucose deprivation is impaired by activated Akt. In addition, AMP-activated protein kinase (AMPK) has been identified as a link between glucose metabolism and the cell cycle, thereby implicating p53 as an essential component of metabolic cell-cycle control.

Expression of transketolase TKTL1 predicts colon and urothelial cancer patient survival: Warburg effect reinterpreted

S Langbein, M Zerilli, A zur Hausen, W Staiger, et al.
British Journal of Cancer (2006) 94, 578–585.
http://dx.doi.org:/10.1038/sj.bjc.6602962

The unique metabolic profile of cancer (aerobic glycolysis) might confer apoptosis resistance and be therapeutically targeted. Compared to normal cells, several human cancers have high mitochondrial membrane potential (DJm) and low expression of the K+ channel Kv1.5, both contributing toapoptosis resistance. Dichloroacetate (DCA) inhibits mitochondrial pyruvate dehydrogenase kinase (PDK), shifts metabolism from glycolysis to glucose oxidation, decreases DJm, increases mitochondrial H2O2, and activates Kv channels in all cancer, but not normal, cells; DCA upregulates Kv1.5 by an NFAT1-dependent mechanism. DCA induces apoptosis, decreases proliferation, and inhibits tumor growth, without apparent toxicity. Molecular inhibition of PDK2 by siRNA mimics DCA. The mitochondria-NFAT-Kv axis and PDK are important therapeutic targets in cancer; the orally available DCA is a promising selective anticancer agent.

Cancer progression and its resistance to treatment depend, at least in part, on suppression of apoptosis. Although mitochondria are recognized as regulators of apoptosis, their importance as targets for cancer therapy has not been adequately explored or clinically exploited. In 1930, Warburg suggested that mitochondrial dysfunction in cancer results in a characteristic metabolic phenotype, that is, aerobic glycolysis (Warburg, 1930). Positron emission tomography (PET) imaging has now confirmed that most malignant tumors have increased glucose uptake and metabolism. This bioenergetic feature is a good marker of cancer but has not been therapeutically pursued..

The small molecule DCA is a metabolic modulator that has been used in humans for decades in the treatment of lactic acidosis and inherited mitochondrial diseases. Without affecting normal cells, DCA reverses the metabolic electrical remodeling that we describe in several cancer lines (hyperpolarized mitochondria, activated NFAT1, downregulated Kv1.5), inducing apoptosis and decreasing tumor growth. DCA in the drinking water at clinically relevant doses for up to 3 months prevents and reverses tumor growth in vivo, without apparent toxicity and without affecting hemoglobin, transaminases, or creatinine levels. The ease of delivery, selectivity, and effectiveness  make DCA an attractive candidate for proapoptotic cancer therapy which can be rapidly translated into phase II–III clinical trials.

A Mitochondria-K+ Channel Axis Is Suppressed in Cancer and Its Normalization Promotes Apoptosis and Inhibits Cancer Growth

Sebastien Bonnet, Stephen L. Archer, Joan Allalunis-Turner, et al.

Cancer Cell Jan 2007; 11: 37–51.
http://dx.doi.org:/10.1016/j.ccr.2006.10.020

Tumor cells, just as other living cells, possess the potential for proliferation, differentiation, cell cycle arrest, and apoptosis. There is a specific metabolic phenotype associated with each of these conditions, characterized by the production of both energy and special substrates necessary for the cells to function in that particular state. Unlike that of normal living cells, the metabolic phenotype of tumor cells supports the proliferative state. Aim: To present the metabolic hypothesis that (1) cell transformation and tumor growth are associated with the activation of metabolic enzymes that increase glucose carbon utilization for nucleic acid synthesis, while enzymes of the lipid and amino acid synthesis pathways are activated in tumor growth inhibition, and (2) phosphorylation and allosteric and transcriptional regulation of intermediary metabolic enzymes and their substrate availability together mediate and sustain cell transformation from one condition to another. Conclusion: Evidence is presented that demonstrates opposite changes in metabolic phenotypes induced by TGF-β, a cell transforming agent, and tumor growth-inhibiting phytochemicals such as genistein and Avemar, or novel synthetic antileukemic drugs such as STI571 (Gleevec).  Intermediary metabolic enzymes that mediate the growth signaling pathways and promote malignant cell transformation may serve as high efficacy nongenetic novel targets for cancer therapies.

A Metabolic Hypothesis of Cell Growth and Death in Pancreatic Cancer

Laszlo G. Boros, Wai-Nang Paul Lee, and Vay Liang W. Go
Pancreas 2002; 24(1):26–33

Clear cell renal cell carcinoma (ccRCC) is the most common pathological subtype of kidney cancer. Here, we integrated an unbiased genome-wide RNA interference screen for ccRCC survival regulators with an analysis of recurrently overexpressed genes in ccRCC to identify new therapeutic targets in this disease. One of the most potent survival regulators, the monocarboxylate transporter MCT4 (SLC16A3), impaired ccRCC viability in all eight ccRCC lines tested and was the seventh most overexpressed gene in a meta-analysis of five ccRCC expression datasets.

MCT4 silencing impaired secretion of lactate generated through glycolysis and induced cell cycle arrest and apoptosis. Silencing MCT4 resulted in intracellular acidosis, and reduction in intracellular ATP production together with partial reversion of the Warburg effect in ccRCC cell lines. Intra-tumoral heterogeneity in the intensity of MCT4 protein expression was observed in primary ccRCCs.

MCT4 protein expression analysis based on the highest intensity of expression in primary ccRCCs was associated with poorer relapse-free survival, whereas modal intensity correlated with Fuhrman nuclear grade. Consistent with the potential selection of subclones enriched for MCT4 expression during disease progression, MCT4 expression was greater at sites of metastatic disease. These data suggest that MCT4 may serve as a novel metabolic target to reverse the Warburg effect and limit disease progression in ccRCC.

Clear cell carcinoma (ccRCC) is the commonest subtype of renal cell carcinoma, accounting for 80% of cases. These tumors are highly resistant to cytotoxic chemotherapy and until recently, systemic treatment options for advanced ccRCC were limited to cytokine based therapies, such as interleukin-2 and interferon-α. Recently, anti-angiogenic drugs and mTOR inhibitors, all targeting the HIF–VEGF axis which is activated in up to 91% of ccRCCs through loss of the VHL tumor suppressor gene [1], have been shown to be effective in metastatic ccRCC [2–5]. Although these drugs increase overall survival to more than 2 years [6], resistance invariably occurs, making the identification of new molecular targets a major clinical need to improve outcomes in patients with metastatic ccRCC.

Genome-wide RNA interference analysis of renal carcinoma survival regulators identifies MCT4 as a Warburg effect metabolic target

Marco Gerlinger, Claudio R Santos, Bradley Spencer-Dene, et al.
J Pathol 2012; 227: 146–156
http://dx.doi.org:/10.1002/path.4006

Hypoxia-inducible factor 1 (HIF-1) plays a key role in the reprogramming of cancer metabolism by activating transcription of genes encoding glucose transporters and glycolytic enzymes, which take up glucose and convert it to lactate; pyruvate dehydrogenase kinase 1, which shunts pyruvate away from the mitochondria; and BNIP3, which triggers selective mitochondrial autophagy. The shift from oxidative to glycolytic metabolism allows maintenance of redox homeostasis and cell survival under conditions of prolonged hypoxia. Many metabolic abnormalities in cancer cells increase HIF-1 activity. As a result, a feed-forward mechanism can be activated that drives HIF-1 activation and may promote tumor progression.

Metastatic cancer is characterized by reprogramming of cellular metabolism leading to increased uptake of glucose for use as both an anabolic and a catabolic substrate. Increased glucose uptake is such a reliable feature that it is utilized clinically to detect metastases by positron emission tomography using 18F-fluorodeoxyglucose (FDG-PET) with a sensitivity of >90% [1]. As with all aspects of cancer biology, the details of metabolic reprogramming differ widely among individual tumors. However, the role of specific signaling pathways and transcription factors in this process is now understood in considerable detail. This review will focus on the involvement of hypoxia-inducible factor 1 (HIF-1) in both mediating metabolic reprogramming and responding to metabolic alterations. The placement of HIF-1 both upstream and downstream of cancer metabolism results in a feed-forward mechanism that may play a major role in the development of the invasive, metastatic, and lethal cancer phenotype.

O2 concentrations are significantly reduced in many human cancers compared with the surrounding normal tissue. The median PO2 in breast cancers is 10 mmHg, as compared with65 mmHg in normal breast tissue. Reduced O2 availability induces HIF-1, which regulates the transcription of hundreds of genes that encode proteins involved in every aspect of cancer biology, including: cell immortalization and stem cell maintenance; genetic instability; glucose and energy metabolism; vascularization; autocrine growth factor signaling; invasion and metastasis; immune evasion; and resistance to chemotherapy and radiation therapy.

HIF-1 is a transcription factor that consists of an O2 regulated HIF-1a and a constitutively expressed HIF-1b subunit. In well-oxygenated cells, HIF-1a is hydroxylated on proline residue 402 (Pro-402) and/or Pro-564 by prolyl hydroxylase domain protein 2 (PHD2), which uses O2 and a-ketoglutarate as substrates in a reaction that generates CO2 and succinate as byproducts. Prolylhydroxylated HIF-1a is bound by the von Hippel–Lindau tumor suppressor protein (VHL), which recruits an E3-ubiquitin ligase that targets HIF-1a for proteasomal degradation (Figure 1a). Asparagine 803 in the transactivation domain is hydroxylated in well-oxygenated cells by factor inhibiting HIF-1 (FIH-1), which blocks the binding of the coactivators p300 and CBP. Under hypoxic conditions, the prolyl and asparaginyl hydroxylation reactions are inhibited by substrate (O2) deprivation and/or the mitochondrial generation of reactive oxygen species (ROS), which may oxidize Fe(II) present in the catalytic center of the hydroxylases.

The finding that acute changes in PO2 increase mitochondrial ROS production suggests that cellular respiration is optimized at physiological PO2 to limit ROS generation and that any deviation in PO2 – up or down – results in increased ROS generation. If hypoxia persists, induction of HIF-1 leads to adaptive mechanisms to reduce ROS and re-establish homeostasis, as described below. Prolyl and asparaginyl hydroxylation provide a molecular mechanism by which changes in cellular oxygenation can be transduced to the nucleus as changes in HIF-1 activity.

HIF-1: upstream and downstream of cancer metabolism

Gregg L Semenza
Current Opinion in Genetics & Development 2010, 20:51–56

This review comes from a themed issue on Genetic and cellular mechanisms of oncogenesis Edited by Tony Hunter and Richard Marais

http://dx.doi.org:/10.1016/j.gde.2009.10.009

Hypoxia-inducible factor 1 (HIF-1) regulates the transcription of many genes involved in key aspects of cancer biology, including immortalization, maintenance of stem cell pools, cellular dedifferentiation, genetic instability, vascularization, metabolic reprogramming, autocrine growth factor signaling, invasion/metastasis, and treatment failure. In animal models, HIF-1 overexpression is associated with increased tumor growth, vascularization, and metastasis, whereas HIF-1 loss-of-function has the opposite effect, thus validating HIF-1 as a target. In further support of this conclusion, immunohistochemical detection of HIF-1a overexpression in biopsy sections is a prognostic factor in many cancers. A growing number of novel anticancer agents have been shown to inhibit HIF-1 through a  variety of molecular mechanisms. Determining which combination of drugs to administer to any given patient remains a major obstacle to improving cancer treatment outcomes.

Intratumoral hypoxia The majority of locally advanced solid tumors contain regions of reduced oxygen availability. Intratumoral hypoxia results when cells are located too far from a functional blood vessel for diffusion of adequate amounts of O2 as a result of rapid cancer cell proliferation and the formation of blood vessels that are structurally and functionally abnormal. In the most extreme case, O2 concentrations are below those required for survival, resulting in cell death and establishing a selection for cancer cells in which apoptotic pathways are inactivated, anti-apoptotic pathways are activated, or invasion/metastasis pathways that promote escape from the hypoxic microenvironment are activated. This hypoxic adaptation may arise by alterations in gene expression or by mutations in the genome or both and is associated with reduced patient survival.

Hypoxia-inducible factor 1 (HIF-1) The expression of hundreds of genes is altered in each cell exposed to hypoxia. Many of these genes are regulated by HIF-1. HIF-1 is a heterodimer formed by the association of an O2-regulated HIF1a subunit with a constitutively expressed HIF-1b subunit. The structurally and functionally related HIF-2a protein also dimerizes with HIF-1b and regulates an overlapping battery of target genes. Under nonhypoxic conditions, HIF-1a (as well as HIF-2a) is subject to O2-dependent prolyl hydroxylation and this modification is required for binding of the von Hippel–Lindau tumor suppressor protein (VHL), which also binds to Elongin C and thereby recruits a ubiquitin ligase complex that targets HIF-1a for ubiquitination and proteasomal degradation. Under hypoxic conditions, the rate of hydroxylation and ubiquitination declines, resulting in accumulation of HIF-1a. Immunohistochemical analysis of tumor biopsies has revealed high levels of HIF-1a in hypoxic but viable tumor cells surrounding areas of necrosis.

Genetic alterations in cancer cells increase HIF-1 activity In the majority of clear-cell renal carcinomas, VHL function is lost, resulting in constitutive activation of HIF-1. After re-introduction of functional VHL, renal carcinoma cell lines are no longer tumorigenic, but can be made tumorigenic by expression of HIF2a in which the prolyl residues that are subject to hydroxylation have been mutated. In addition to VHL loss-of-function, many other genetic alterations that inactivate tumor suppressors

Evaluation of HIF-1 inhibitors as anticancer agents

Gregg L. Semenza
Drug Discovery Today Oct 2007; 12(19/20).
http://dx.doi.org:/10.1016/j.drudis.2007.08.006

Hypoxia-inducible factor-1 (HIF-1), which is present at high levels in human tumors, plays crucial roles in tumor promotion by upregulating its target genes, which are involved in anaerobic energy metabolism, angiogenesis, cell survival, cell invasion, and drug resistance. Therefore, it is apparent that the inhibition of HIF-1 activity may be a strategy for treating cancer. Recently, many efforts to develop new HIF-1-targeting agents have been made by both academic and pharmaceutical industry laboratories. The future success of these efforts will be a new class of HIF-1-targeting anticancer agents, which would improve the prognoses of many cancer patients. This review focuses on the potential of HIF-1 as a target molecule for anticancer therapy, and on possible strategies to inhibit HIF-1 activity. In addition, we introduce YC-1 as a new anti-HIF-1, anticancer agent. Although YC-1 was originally developed as a potential therapeutic agent for thrombosis and hypertension, recent studies demonstrated that YC-1 suppressed HIF-1 activity and vascular endothelial growth factor expression in cancer cells. Moreover, it halted tumor growth in immunodeficient mice without serious toxicity during the treatment period. Thus, we propose that YC-1 is a good lead compound for the development of new anti-HIF-1, anticancer agents.

Although many anticancer regimens have been introduced to date, their survival benefits are negligible, which is the reason that a more innovative treatment is required. Basically, the identification of the specific molecular features of tumor promotion has allowed for rational drug discovery in cancer treatment, and drugs have been screened based upon the modulation of specific molecular targets in tumor cells. Target-based drugs should satisfy the following two conditions.

First, they must act by a described mechanism.

Second, they must reduce tumor growth in vivo, associated with this mechanism.

Many key factors have been found to be involved in the multiple steps of cell growth signal-transduction pathways. Targeting these factors offers a strategy for preventing tumor growth; for example, competitors or antibodies blocking ligand–receptor interaction, and receptor tyrosine kinase inhibitors, downstream pathway inhibitors (i.e., RAS farnesyl transferase inhibitors, mitogen-activated protein kinase and mTOR inhibitors), and cell-cycle arresters (i.e., cyclin-dependent kinase inhibitors) could all be used to inhibit tumor growth.

In addition to the intracellular events, tumor environmental factors should be considered to treat solid tumors. Of these, hypoxia is an important cancer-aggravating factor because it contributes to the progression of a more malignant phenotype, and to the acquisition of resistance to radiotherapy and chemotherapy. Thus, transcription factors that regulate these hypoxic events are good targets for anticancer therapy and in particular HIF-1 is one of most compelling targets. In this paper, we introduce the roles of HIF-1 in tumor promotion and provide a summary of new anticancer strategies designed to inhibit HIF-1 activity.

New anticancer strategies targeting HIF-1

Eun-Jin Yeo, Yang-Sook Chun, Jong-Wan Park
Biochemical Pharmacology 68 (2004) 1061–1069
http://dx.doi.org:/10.1016/j.bcp.2004.02.040

Classical work in tumor cell metabolism focused on bioenergetics, particularly enhanced glycolysis and suppressed oxidative phosphorylation (the ‘Warburg effect’). But the biosynthetic activities required to create daughter cells are equally important for tumor growth, and recent studies are now bringing these pathways into focus. In this review, we discuss how tumor cells achieve high rates of nucleotide and fatty acid synthesis, how oncogenes and tumor suppressors influence these activities, and how glutamine metabolism enables macromolecular synthesis in proliferating cells.

Otto Warburg’s demonstration that tumor cells rapidly use glucose and convert the majority of it to lactate is still the most fundamental and enduring observation in tumor metabolism. His work, which ushered in an era of study on tumor metabolism focused on the relationship between glycolysis and cellular bioenergetics, has been revisited and expanded by generations of tumor biologists. It is now accepted that a high rate of glucose metabolism, exploited clinically by 18FDGPET scanning, is a metabolic hallmark of rapidly dividing cells, correlates closely with transformation, and accounts for a significant percentage of ATP generated during cell proliferation. A ‘metabolic transformation’ is required for tumorigenesis. Research over the past few years has reinforced this idea, revealing the conservation of metabolic activities among diverse tumor types, and proving that oncogenic mutations can promote metabolic autonomy by driving nutrient uptake to levels that often exceed those required for cell growth and proliferation.

In order to engage in replicative division, a cell must duplicate its genome, proteins, and lipids and assemble the components into daughter cells; in short, it must become a factory for macromolecular biosynthesis. These activities require that cells take up extracellular nutrients like glucose and glutamine and allocate them into metabolic pathways that convert them into biosynthetic precursors (Figure 1). Tumor cells can achieve this phenotype through changes in the expression of enzymes that determine metabolic flux rates, including nutrient transporters and enzymes [8– 10]. Current studies in tumor metabolism are revealing novel mechanisms for metabolic control, establishing which enzyme isoforms facilitate the tumor metabolic phenotype, and suggesting new targets for cancer therapy.

The ongoing challenge in tumor cell metabolism is to understand how individual pathways fit together into the global metabolic phenotype of cell growth. Here we discuss two biosynthetic activities required by proliferating tumor cells: production of ribose-5 phosphate for nucleotide biosynthesis and production of fatty acids for lipid biosynthesis. Nucleotide and lipid biosynthesis share three important characteristics.

  • First, both use glucose as a carbon source.
  • Second, both consume TCA cycle intermediates, imposing the need for a mechanism to replenish the cycle.
  • Third, both require reductive power in the form of NADPH.

In this Essay, we discuss the possible drivers, advantages, and potential liabilities of the altered metabolism of cancer cells (Figure 1, not shown). Although our emphasis on the Warburg effect reflects the focus of the field, we would also like to encourage a broader approach to the study of cancer metabolism that takes into account the contributions of all interconnected small molecule pathways of the cell.

The Tumor Microenvironment Selects for Altered Metabolism One compelling idea to explain the Warburg effect is that the altered metabolism of cancer cells confers a selective advantage for survival and proliferation in the unique tumor microenvironment. As the early tumor expands, it outgrows the diffusion limits of its local blood supply, leading to hypoxia and stabilization of the hypoxia-inducible transcription factor, HIF. HIF initiates a transcriptional program that provides multiple solutions to hypoxic stress (reviewed in Kaelin and Ratcliffe, 2008). Because a decreased dependence on aerobic respiration becomes advantageous, cell metabolism is shifted toward glycolysis by the increased expression of glycolytic enzymes, glucose transporters, and inhibitors of mitochondrial metabolism. In addition, HIF stimulates angiogenesis (the formation of new blood vessels) by upregulating several factors, including most prominently vascular endothelial growth factor (VEGF).

Blood vessels recruited to the tumor microenvironment, however, are disorganized, may not deliver blood effectively, and therefore do not completely alleviate hypoxia (reviewed in Gatenby and Gillies, 2004). The oxygen levels within a tumor vary both spatially and temporally, and the resulting rounds of fluctuating oxygen levels potentially select for tumors that constitutively upregulate glycolysis. Interestingly, with the possible exception of tumors that have lost the von Hippel-Lindau protein (VHL), which normally mediates degradation of HIF, HIF is still coupled to oxygen levels, as evident from the heterogeneity of HIF expression within the tumor microenvironment. Therefore, the Warburg effect—that is, an uncoupling of glycolysis from oxygen levels—cannot be explained solely by upregulation of HIF. Other molecular mechanisms are likely to be important, such as the metabolic changes induced by oncogene activation and tumor suppressor loss.

Oncogene Activation Drives Changes in Metabolism Not only may the tumor microenvironment select for a deranged metabolism, but oncogene status can also drive metabolic changes. Since Warburg’s time, the biochemical study of cancer metabolism has been overshadowed by efforts to identify the mutations that contribute to cancer initiation and progression. Recent work, however, has demonstrated that the key components of the Warburg effect—

  • increased glucose consumption,
  • decreased oxidative phosphorylation, and
  • accompanying lactate production—
  • are also distinguishing features of oncogene activation.

The signaling molecule Ras, a powerful oncogene when mutated, promotes glycolysis (reviewed in Dang and Semenza, 1999; Ramanathan et al., 2005). Akt kinase, a well-characterized downstream effector of insulin signaling, reprises its role in glucose uptake and utilization in the cancer setting (reviewed in Manning and Cantley, 2007), whereas the Myc transcription factor upregulates the expression of various metabolic genes (reviewed in Gordan et al., 2007). The most parsimonious route to tumorigenesis may be activation of key oncogenic nodes that execute a proliferative program, of which metabolism may be one important arm. Moreover, regulation of metabolism is not exclusive to oncogenes.

Cancer Cell Metabolism: Warburg & Beyond

Hsu PP & Sabatini DM
Cell  Sep 5, 2008; 134, 703-705
http://dx.doi.org:/10.1016/j.cell.2008.08.021

Tumor cells respond to growth signals by the activation of protein kinases, altered gene expression and significant modifications in substrate flow and redistribution among biosynthetic pathways. This results in a proliferating phenotype with altered cellular function. These transformed cells exhibit unique anabolic characteristics, which includes increased and preferential utilization of glucose through the non-oxidative steps of the pentose cycle for nucleic acid synthesis but limited de novo fatty  acid   synthesis   and   TCA   cycle   glucose   oxidation. This  primarily nonoxidative anabolic profile reflects an undifferentiated highly proliferative aneuploid cell phenotype and serves as a reliable metabolic biomarker to determine cell proliferation rate and the level of cell transformation/differentiation in response to drug treatment.

Novel drugs effective in particular cancers exert their anti-proliferative effects by inducing significant reversions of a few specific non-oxidative anabolic pathways. Here we present evidence that cell transformation of various mechanisms is sustained by a unique disproportional substrate distribution between the two branches of the pentose cycle for nucleic acid synthesis, glycolysis and the TCA cycle for fatty acid synthesis and glucose oxidation. This can be demonstrated by the broad labeling and unique specificity of [1,2-13C2]glucose to trace a large number of metabolites in the metabolome. Stable isotope-based dynamic metabolic profiles (SIDMAP) serve the drug discovery process by providing a powerful new tool that integrates the metabolome into a functional genomics approach to developing new drugs. It can be used in screening kinases and their metabolic targets, which can therefore be more efficiently characterized, speeding up and improving drug testing, approval and labeling processes by saving trial and error type study costs in drug testing.

Metabolic Biomarker and Kinase Drug Target Discovery in Cancer Using Stable Isotope-Based Dynamic Metabolic Profiling (SIDMAP)

László G. Boros, Daniel J. Brackett and George G. Harrigan
Current Cancer Drug Targets, 2003, 3, 447-455 447

Pyruvate constitutes a critical branch point in cellular carbon metabolism. We have identified two proteins, Mpc1 and Mpc2, as essential for mitochondrial pyruvate transport in yeast, Drosophila, and humans. Mpc1 and Mpc2 associate to form an ~150 kilodalton complex in the inner mitochondrial membrane. Yeast and Drosophila mutants lacking MPC1 display impaired pyruvate metabolism, with an accumulation of upstream metabolites and a depletion of tricarboxylic acid cycle intermediates. Loss of yeast Mpc1 results in defective mitochondrial pyruvate uptake, while silencing of MPC1 or MPC2 in mammalian cells impairs pyruvate oxidation. A point mutation in MPC1 provides resistance to a known inhibitor of the mitochondrial pyruvate carrier. Human genetic studies of three families with children suffering from lactic acidosis and hyperpyruvatemia revealed a causal locus that mapped to MPC1, changing single amino acids that are conserved throughout eukaryotes. These data demonstrate that Mpc1 and Mpc2 form an essential part of the mitochondrial pyruvate carrier.

A Mitochondrial Pyruvate Carrier Required for Pyruvate Uptake in Yeast, Drosophila , and Humans

Daniel K. Bricker, Eric B. Taylor, John C. Schell, Thomas Orsak, et al.
Science Express 24 May 2012
http://dx.doi.org:/10.1126/science.1218099

Adenosine deaminase acting on RNA (ADAR) enzymes convert adenosine (A) to inosine (I) in double-stranded (ds) RNAs. Since Inosine is read as Guanosine, the biological consequence of ADAR enzyme activity is an A/G conversion within RNA molecules. A-to-I editing events can occur on both coding and non-coding RNAs, including microRNAs (miRNAs), which are small regulatory RNAs of ~20–23 nucleotides that regulate several cell processes by annealing to target mRNAs and inhibiting their translation. Both miRNA precursors and mature miRNAs undergo A-to-I RNA editing, affecting the miRNA maturation process and activity. ADARs can also edit 3′ UTR of mRNAs, further increasing the interplay between mRNA targets and miRNAs. In this review, we provide a general overview of the ADAR enzymes and their mechanisms of action as well as miRNA processing and function. We then review the more recent findings about the impact of ADAR-mediated activity on the miRNA pathway in terms of biogenesis, target recognition, and gene expression regulation.

Review ADAR Enzyme and miRNA Story: A Nucleotide that Can Make the Difference 

Sara Tomaselli, Barbara Bonamassa, Anna Alisi, Valerio Nobili, Franco Locatelli and Angela Gallo
Int. J. Mol. Sci. 19 Nov 2013; 14, 22796-22816 http://dx.doi.org:/10.3390/ijms141122796

The fermented wheat germ extract (FWGE) nutraceutical (Avemar™), manufactured under “good manufacturing practice” conditions and, fulfilling the self-affirmed “generally recognized as safe” status in the United States, has been approved as a “dietary food for special medical purposes for cancer patients” in Europe. In this paper, we report the adjuvant use of this nutraceutical in the treatment of high-risk skin melanoma patients. Methods: In a randomized, pilot, phase II clinical trial, the efficacy of dacarbazine (DTIC)-based adjuvant chemotherapy on survival parameters of melanoma patients was compared to that of the same treatment supplemented with a 1-year long administration of FWGE. Results: At the end of an additional 7-year-long follow-up period, log-rank analyses (Kaplan-Meier estimates) showed significant differences in both progression-free (PFS) and overall survival (OS) in favor of the FWGE group. Mean PFS: 55.8 months (FWGE group) versus 29.9 months (control group), p  0.0137. Mean OS: 66.2 months (FWGE group) versus 44.7 months (control group), p < 0.0298. Conclusions: The inclusion of Avemar into the adjuvant protocols of high-risk skin melanoma patients is highly recommended.

Adjuvant Fermented Wheat Germ Extract (Avemar™) Nutraceutical Improves Survival of High-Risk Skin Melanoma Patients: A Randomized, Pilot, Phase II Clinical Study with a 7-Year Follow-Up

LV Demidov, LV Manziuk, GY Kharkevitch, NA Pirogova, and EV Artamonova
Cancer Biotherapy & Radiopharmaceuticals 2008; 23(4)
http://dx.doi.org:/10.1089/cbr.2008.0486

Cancer cells possess unique metabolic signatures compared to normal cells, including shifts in aerobic glycolysis, glutaminolysis, and de novo biosynthesis of macromolecules. Targeting these changes with agents (drugs and dietary components) has been employed as strategies to reduce the complications associated with tumorigenesis. This paper highlights the ability of several food components to suppress tumor-specific metabolic pathways, including increased expression of glucose transporters, oncogenic tyrosine kinase, tumor-specific M2-type pyruvate kinase, and fatty acid synthase, and the detection of such effects using various metabonomic technologies, including liquid chromatography/mass spectrometry (LC/MS) and stable isotope-labeled MS. Stable isotope-mediated tracing technologies offer exciting opportunities for defining specific target(s) for food components. Exposures, especially during the early transition phase from normal to cancer, are critical for the translation of knowledge about food components into effective prevention strategies. Although appropriate dietary exposures needed to alter cellular metabolism remain inconsistent and/or ill-defined, validated metabonomic biomarkers for dietary components hold promise for establishing effective strategies for cancer prevention.

Bioactive Food Components and Cancer-Specific Metabonomic Profiles

Young S. Kim and John A. Milner
Journal of Biomedicine and Biotechnology 2011, Art ID 721213, 9 pages
http://dx.doi.org:/10.1155/2011/721213

This reviewer poses the following observation.  The importance of the pyridine nucleotide reduced/oxidized ratio has not been alluded to here, but the importance cannot be understated. It has relevance to the metabolic functions of anabolism and catabolism of the visceral organs.  The importance of this has ties to the pentose monophosphate pathway. The importance of the pyridine nucleotide transhydrogenase reaction remains largely unexplored.  In reference to the NAD-redox state, the observation was made by Nathan O. Kaplan that the organs may be viewed with respect to their primary functions in anabolic or high energy catabolic activities. Thus we find that the endocrine organs are largely tied to anabolic functioning, and to NADP, whereas cardiac and skeletal muscle are highly dependent on NAD. The consequence of this observed phenomenon appears to be related to a difference in the susceptibility to malignant transformation.  In the case of the gastrointestinal tract, the rate of turnover of the epithelium is very high. However, with the exception of the liver, there is no major activity other than cell turnover. In the case of the liver, there is a major commitment to synthesis of lipids, storage of fuel, and synthesis of proteins, which is largely anabolic, but there is also a major activity in detoxification, which is not.  In addition, the liver has a double circulation. As a result, a Zahn infarct is uncommon.  Now we might also consider the heart.  The heart is a muscle syncytium with a high need for oxygen.  Cutting of the oxygen supply makes the myocytes vulnerable to ischemic insult and abberant rhythm abnormalities.  In addition, the cardiomyocyte can take up lactic acid from the circulation for fuel, which is tied to the utilization of lactate from vigorous skeletal muscle activity.  The skeletal muscle is tied to glycolysis in normal function, which has a poor generation of ATP, so that the recycling of excess lactic acid is required by cardiac muscle and hepatocytes.  This has not been a part of the discussion, but this reviewer considers it important to remember in considering the organ-specific tendencies to malignant transformation.

Comment (Aurelian Udristioiu):

Otto Warburg observed that many cancers lose their capacity for mitochondrial respiration, limiting ATP production to anaerobic glycolytic pathways. The phenomenon is particularly prevalent in aggressive malignancies, most of which are also hypoxic [1].
Hypoxia induces a stochastic imbalance between the numbers of reduced mitochondrial species vs. available oxygen, resulting in increased reactive oxygen species (ROS) whose toxicity can lead to apoptotic cell death.
Mechanism involves inhibition of glycolytic ATP production via a Randle-like cycle while increased uncoupling renders cancers unable to produce compensatory ATP from respiration-.generation in the presence of intact tricarboxylic acid (TCA) enzyme.
One mitochondrial adaptation to increased ROS is over-expression of the uncoupling protein 2 (UCP2) that has been reported in multiple human cancer cell lines [2-3]. Increased UCP2 expression was also associated with reduced ATP production in malignant oxyphilic mouse leukemia and human lymphoma cell lines [4].
Hypoxia reduces the ability of cells to maintain their energy levels, because less ATP is obtained from glycolysis than from oxidative phosphorylation. Cells adapt to hypoxia by activating the expression of mutant genes in glycolysis.
-Severe hypoxia causes a high mutation rate, resulting in point mutations that may be explained by reduced DNA mismatch repairing activity.
The most direct induction of apoptosis caused by hypoxia is determined by the inhibition of the electron carrier chain from the inner membrane of the mitochondria. The lack of oxygen inhibits the transport of protons and thereby causes a decrease in membrane potential. Cell survival under conditions of mild hypoxia is mediated by phosphoinositide-3 kinase (PIK3) using severe hypoxia or anoxia, and then cells initiate a cascade of events that lead to apoptosis [5].
After DNA damage, a very important regulator of apoptosis is the p53 protein. This tumor suppressor gene has mutations in over 60% of human tumors and acts as a suppressor of cell division. The growth-suppressive effects of p53 are considered to be mediated through the transcriptional trans-activation activity of the protein. In addition to the maturational state of the clonal tumor, the prognosis of patients with CLL is dependent of genetic changes within the neoplastic cell population.

1.Warburg O. On the origin of cancer cells. Science 1956; 123 (3191):309-314
PubMed Abstract ; Publisher Full Text

2.Giardina TM, Steer JH, Lo SZ, Joyce DA. Uncoupling protein-2 accumulates rapidly in the inner mitochondrial membrane during mitochondrial reactive oxygen stress in macrophages. Biochim Biophys Acta 2008, 1777(2):118-129. PubMed Abstract | Publisher Full Text

3. Horimoto M, Resnick MB, Konkin TA, Routhier J, Wands JR, Baffy G. Expression of uncoupling protein-2 in human colon cancer. Clin Cancer Res 2004; 10 (18 Pt1):6203-6207. PubMed Abstract | Publisher Full Text

4. Randle PJ, England PJ, Denton RM. Control of the tricarboxylate cycle and it interactions with glycolysis during acetate utilization in rat heart. Biochem J 1970; 117(4):677-695. PubMed Abstract | PubMed Central Full Text

5. Gillies RJ, Robey I, Gatenby RA. Causes and consequences of increased glucose metabolism of cancers. J Nucl Med 2008; 49(Suppl 2):24S-42S. PubMed Abstract | Publisher Full Text

Shortened version of Comment –

Hypoxia induces a stochastic imbalance between the numbers of reduced mitochondrial species vs. available oxygen, resulting in increased reactive oxygen species (ROS) whose toxicity can lead to apoptotic cell death.
Mechanism involves inhibition of glycolytic ATP production via a Randle-like cycle while increased uncoupling renders cancers unable to produce compensatory ATP from respiration-.generation in the presence of intact tricarboxylic acid (TCA) enzyme.
One mitochondrial adaptation to increased ROS is over-expression of the uncoupling protein 2 (UCP2) that has been reported in multiple human cancer cell lines. Increased UCP2 expression was also associated with reduced ATP production in malignant oxyphilic mouse leukemia and human lymphoma cell lines.
Severe hypoxia causes a high mutation rate, resulting in point mutations that may be explained by reduced DNA mismatch repairing activity.

Read Full Post »

Artificial Intelligence Versus the Scientist: Who Will Win?

Will DARPA Replace the Human Scientist: Not So Fast, My Friend!

Writer, Curator: Stephen J. Williams, Ph.D.

Article ID #168: Artificial Intelligence Versus the Scientist: Who Will Win?. Published on 3/2/2015

WordCloud Image Produced by Adam Tubman

scientistboxingwithcomputer

Last month’s issue of Science article by Jia You “DARPA Sets Out to Automate Research”[1] gave a glimpse of how science could be conducted in the future: without scientists. The article focused on the U.S. Defense Advanced Research Projects Agency (DARPA) program called ‘Big Mechanism”, a $45 million effort to develop computer algorithms which read scientific journal papers with ultimate goal of extracting enough information to design hypotheses and the next set of experiments,

all without human input.

The head of the project, artificial intelligence expert Paul Cohen, says the overall goal is to help scientists cope with the complexity with massive amounts of information. As Paul Cohen stated for the article:

“‘

Just when we need to understand highly connected systems as systems,

our research methods force us to focus on little parts.

                                                                                                                                                                                                               ”

The Big Mechanisms project aims to design computer algorithms to critically read journal articles, much as scientists will, to determine what and how the information contributes to the knowledge base.

As a proof of concept DARPA is attempting to model Ras-mutation driven cancers using previously published literature in three main steps:

  1. Natural Language Processing: Machines read literature on cancer pathways and convert information to computational semantics and meaning

One team is focused on extracting details on experimental procedures, using the mining of certain phraseology to determine the paper’s worth (for example using phrases like ‘we suggest’ or ‘suggests a role in’ might be considered weak versus ‘we prove’ or ‘provide evidence’ might be identified by the program as worthwhile articles to curate). Another team led by a computational linguistics expert will design systems to map the meanings of sentences.

  1. Integrate each piece of knowledge into a computational model to represent the Ras pathway on oncogenesis.
  2. Produce hypotheses and propose experiments based on knowledge base which can be experimentally verified in the laboratory.

The Human no Longer Needed?: Not So Fast, my Friend!

The problems the DARPA research teams are encountering namely:

  • Need for data verification
  • Text mining and curation strategies
  • Incomplete knowledge base (past, current and future)
  • Molecular biology not necessarily “requires casual inference” as other fields do

Verification

Notice this verification step (step 3) requires physical lab work as does all other ‘omics strategies and other computational biology projects. As with high-throughput microarray screens, a verification is needed usually in the form of conducting qPCR or interesting genes are validated in a phenotypical (expression) system. In addition, there has been an ongoing issue surrounding the validity and reproducibility of some research studies and data.

See Importance of Funding Replication Studies: NIH on Credibility of Basic Biomedical Studies

Therefore as DARPA attempts to recreate the Ras pathway from published literature and suggest new pathways/interactions, it will be necessary to experimentally validate certain points (protein interactions or modification events, signaling events) in order to validate their computer model.

Text-Mining and Curation Strategies

The Big Mechanism Project is starting very small; this reflects some of the challenges in scale of this project. Researchers were only given six paragraph long passages and a rudimentary model of the Ras pathway in cancer and then asked to automate a text mining strategy to extract as much useful information. Unfortunately this strategy could be fraught with issues frequently occurred in the biocuration community namely:

Manual or automated curation of scientific literature?

Biocurators, the scientists who painstakingly sort through the voluminous scientific journal to extract and then organize relevant data into accessible databases, have debated whether manual, automated, or a combination of both curation methods [2] achieves the highest accuracy for extracting the information needed to enter in a database. Abigail Cabunoc, a lead developer for Ontario Institute for Cancer Research’s WormBase (a database of nematode genetics and biology) and Lead Developer at Mozilla Science Lab, noted, on her blog, on the lively debate on biocuration methodology at the Seventh International Biocuration Conference (#ISB2014) that the massive amounts of information will require a Herculaneum effort regardless of the methodology.

Although I will have a future post on the advantages/disadvantages and tools/methodologies of manual vs. automated curation, there is a great article on researchinformation.infoExtracting More Information from Scientific Literature” and also see “The Methodology of Curation for Scientific Research Findings” and “Power of Analogy: Curation in Music, Music Critique as a Curation and Curation of Medical Research Findings – A Comparison” for manual curation methodologies and A MOD(ern) perspective on literature curation for a nice workflow paper on the International Society for Biocuration site.

The Big Mechanism team decided on a full automated approach to text-mine their limited literature set for relevant information however was able to extract only 40% of relevant information from these six paragraphs to the given model. Although the investigators were happy with this percentage most biocurators, whether using a manual or automated method to extract information, would consider 40% a low success rate. Biocurators, regardless of method, have reported ability to extract 70-90% of relevant information from the whole literature (for example for Comparative Toxicogenomics Database)[3-5].

Incomplete Knowledge Base

In an earlier posting (actually was a press release for our first e-book) I had discussed the problem with the “data deluge” we are experiencing in scientific literature as well as the plethora of ‘omics experimental data which needs to be curated.

Tackling the problem of scientific and medical information overload

pubmedpapersoveryears

Figure. The number of papers listed in PubMed (disregarding reviews) during ten year periods have steadily increased from 1970.

Analyzing and sharing the vast amounts of scientific knowledge has never been so crucial to innovation in the medical field. The publication rate has steadily increased from the 70’s, with a 50% increase in the number of original research articles published from the 1990’s to the previous decade. This massive amount of biomedical and scientific information has presented the unique problem of an information overload, and the critical need for methodology and expertise to organize, curate, and disseminate this diverse information for scientists and clinicians. Dr. Larry Bernstein, President of Triplex Consulting and previously chief of pathology at New York’s Methodist Hospital, concurs that “the academic pressures to publish, and the breakdown of knowledge into “silos”, has contributed to this knowledge explosion and although the literature is now online and edited, much of this information is out of reach to the very brightest clinicians.”

Traditionally, organization of biomedical information has been the realm of the literature review, but most reviews are performed years after discoveries are made and, given the rapid pace of new discoveries, this is appearing to be an outdated model. In addition, most medical searches are dependent on keywords, hence adding more complexity to the investigator in finding the material they require. Third, medical researchers and professionals are recognizing the need to converse with each other, in real-time, on the impact new discoveries may have on their research and clinical practice.

These issues require a people-based strategy, having expertise in a diverse and cross-integrative number of medical topics to provide the in-depth understanding of the current research and challenges in each field as well as providing a more conceptual-based search platform. To address this need, human intermediaries, known as scientific curators, are needed to narrow down the information and provide critical context and analysis of medical and scientific information in an interactive manner powered by web 2.0 with curators referred to as the “researcher 2.0”. This curation offers better organization and visibility to the critical information useful for the next innovations in academic, clinical, and industrial research by providing these hybrid networks.

Yaneer Bar-Yam of the New England Complex Systems Institute was not confident that using details from past knowledge could produce adequate roadmaps for future experimentation and noted for the article, “ “The expectation that the accumulation of details will tell us what we want to know is not well justified.”

In a recent post I had curated findings from four lung cancer omics studies and presented some graphic on bioinformatic analysis of the novel genetic mutations resulting from these studies (see link below)

Multiple Lung Cancer Genomic Projects Suggest New Targets, Research Directions for

Non-Small Cell Lung Cancer

which showed, that while multiple genetic mutations and related pathway ontologies were well documented in the lung cancer literature there existed many significant genetic mutations and pathways identified in the genomic studies but little literature attributed to these lung cancer-relevant mutations.

KEGGinliteroanalysislungcancer

  This ‘literomics’ analysis reveals a large gap between our knowledge base and the data resulting from large translational ‘omic’ studies.

Different Literature Analyses Approach Yeilding

A ‘literomics’ approach focuses on what we don NOT know about genes, proteins, and their associated pathways while a text-mining machine learning algorithm focuses on building a knowledge base to determine the next line of research or what needs to be measured. Using each approach can give us different perspectives on ‘omics data.

Deriving Casual Inference

Ras is one of the best studied and characterized oncogenes and the mechanisms behind Ras-driven oncogenenis is highly understood.   This, according to computational biologist Larry Hunt of Smart Information Flow Technologies makes Ras a great starting point for the Big Mechanism project. As he states,” Molecular biology is a good place to try (developing a machine learning algorithm) because it’s an area in which common sense plays a minor role”.

Even though some may think the project wouldn’t be able to tackle on other mechanisms which involve epigenetic factors UCLA’s expert in causality Judea Pearl, Ph.D. (head of UCLA Cognitive Systems Lab) feels it is possible for machine learning to bridge this gap. As summarized from his lecture at Microsoft:

“The development of graphical models and the logic of counterfactuals have had a marked effect on the way scientists treat problems involving cause-effect relationships. Practical problems requiring causal information, which long were regarded as either metaphysical or unmanageable can now be solved using elementary mathematics. Moreover, problems that were thought to be purely statistical, are beginning to benefit from analyzing their causal roots.”

According to him first

1) articulate assumptions

2) define research question in counter-inference terms

Then it is possible to design an inference system using calculus that tells the investigator what they need to measure.

To watch a video of Dr. Judea Pearl’s April 2013 lecture at Microsoft Research Machine Learning Summit 2013 (“The Mathematics of Causal Inference: with Reflections on Machine Learning”), click here.

The key for the Big Mechansism Project may me be in correcting for the variables among studies, in essence building a models system which may not rely on fully controlled conditions. Dr. Peter Spirtes from Carnegie Mellon University in Pittsburgh, PA is developing a project called the TETRAD project with two goals: 1) to specify and prove under what conditions it is possible to reliably infer causal relationships from background knowledge and statistical data not obtained under fully controlled conditions 2) develop, analyze, implement, test and apply practical, provably correct computer programs for inferring causal structure under conditions where this is possible.

In summary such projects and algorithms will provide investigators the what, and possibly the how should be measured.

So for now it seems we are still needed.

References

  1. You J: Artificial intelligence. DARPA sets out to automate research. Science 2015, 347(6221):465.
  2. Biocuration 2014: Battle of the New Curation Methods [http://blog.abigailcabunoc.com/biocuration-2014-battle-of-the-new-curation-methods]
  3. Davis AP, Johnson RJ, Lennon-Hopkins K, Sciaky D, Rosenstein MC, Wiegers TC, Mattingly CJ: Targeted journal curation as a method to improve data currency at the Comparative Toxicogenomics Database. Database : the journal of biological databases and curation 2012, 2012:bas051.
  4. Wu CH, Arighi CN, Cohen KB, Hirschman L, Krallinger M, Lu Z, Mattingly C, Valencia A, Wiegers TC, John Wilbur W: BioCreative-2012 virtual issue. Database : the journal of biological databases and curation 2012, 2012:bas049.
  5. Wiegers TC, Davis AP, Mattingly CJ: Collaborative biocuration–text-mining development task for document prioritization for curation. Database : the journal of biological databases and curation 2012, 2012:bas037.

Other posts on this site on include: Artificial Intelligence, Curation Methodology, Philosophy of Science

Inevitability of Curation: Scientific Publishing moves to embrace Open Data, Libraries and Researchers are trying to keep up

A Brief Curation of Proteomics, Metabolomics, and Metabolism

The Methodology of Curation for Scientific Research Findings

Scientific Curation Fostering Expert Networks and Open Innovation: Lessons from Clive Thompson and others

The growing importance of content curation

Data Curation is for Big Data what Data Integration is for Small Data

Cardiovascular Original Research: Cases in Methodology Design for Content Co-Curation The Art of Scientific & Medical Curation

Exploring the Impact of Content Curation on Business Goals in 2013

Power of Analogy: Curation in Music, Music Critique as a Curation and Curation of Medical Research Findings – A Comparison

conceived: NEW Definition for Co-Curation in Medical Research

Reconstructed Science Communication for Open Access Online Scientific Curation

Search Results for ‘artificial intelligence’

 The Simple Pictures Artificial Intelligence Still Can’t Recognize

Data Scientist on a Quest to Turn Computers Into Doctors

Vinod Khosla: “20% doctor included”: speculations & musings of a technology optimist or “Technology will replace 80% of what doctors do”

Where has reason gone?

Read Full Post »

The Reconstruction of Life Processes requires both Genomics and Metabolomics to explain Phenotypes and Phylogenetics

Writer and Curator: Larry H. Bernstein, MD, FCAP 

 

phylogenetics

phylogenetics

http://upload.wikimedia.org/wikipedia/commons/thumb/1/12/CollapsedtreeLabels-simplified.svg/200px-CollapsedtreeLabels-simplified.svg.png

 

This discussion that completes and is an epicrisis (summary and critical evaluation) of the series of discussions that preceded it.

  1. Innervation of Heart and Heart Rate
  2. Action of hormones on the circulation
  3. Allogeneic Transfusion Reactions
  4. Graft-versus Host reaction
  5. Unique problems of perinatal period
  6. High altitude sickness
  7. Deep water adaptation
  8. Heart-Lung-and Kidney
  9. Acute Lung Injury

The concept inherent in this series is that the genetic code is an imprint that is translated into a message.  It is much the same as a blueprint, or a darkroom photographic image that has to be converted to a print. It is biologically an innovation of evolutionary nature because it establishes a simple and reproducible standard for the transcription of the message through the transcription of the message using strings of nucleotides (oligonucleotides) that systematically transfer the message through ribonucleotides that communicate in the cytoplasm with the cytoskeleton based endoplasmic reticulum (ER), composing a primary amino acid sequence.  This process is a quite simple and convenient method of biological activity.  However, the simplicity ends at this step.  The metabolic components of the cell are organelles consisting of lipoprotein membranes and a cytosol which have particularly aligned active proteins, as in the inner membrane of the mitochondrion, or as in the liposome or phagosome, or the structure of the  ER, each of which is critical for energy transduction and respiration, in particular, for the mitochondria, cellular remodeling or cell death, with respect to the phagosome, and construction of proteins with respect to the ER, and anaerobic glycolysis and the hexose monophosphate shunt in the cytoplasmic domain.  All of this refers to structure and function, not to leave out the membrane assigned transport of inorganic, and organic ions (electrolytes and metabolites).

I have identified a specific role of the ER, the organelles, and cellular transactions within and between cells that is orchestrated.  But what I have outlined is a somewhat limited and rigid model that does not reach into the dynamics of cellular transactions.  The DNA has expression that may be old, no longer used messages, and this is perhaps only part of a significant portion of “dark matter”.  There is also nuclear DNA that is enmeshed with protein, mRNA that is a copy of DNA, and mDNA  is copied to ribosomal RNA (rRNA).  There is also rDNA. The classic model is DNA to RNA to protein.  However, there is also noncoding RNA, which plays an important role in regulation of transcription.

This has been discussed in other articles.  But the important point is that proteins have secondary structure through disulfide bonds, which is determined by position of sulfur amino acids, and by van der Waal forces, attraction and repulsion. They have tertiary structure, which is critical for 3-D structure.  When like subunits associate, or dissimilar oligomers, then you have heterodimers and oligomers.  These constructs that have emerged over time interact with metabolites within the cell, and also have an important interaction with the extracellular environment.

When you take this into consideration then a more complete picture emerges. The primitive cell or the multicellular organism lives in an environment that has the following characteristics – air composition, water and salinity, natural habitat, temperature, exposure to radiation, availability of nutrients, and exposure to chemical toxins or to predators.  In addition, there is a time dimension that proceeds from embryonic stage to birth in mammals, a rapid growth phase, a tapering, and a decline.  The time span is determined by body size, fluidity of adaptation, and environmental factors.  This is covered in great detail in this work.  The last two pieces are in the writing stage that completes the series. Much content has already be presented in previous articles.

The function of the heart, kidneys and metabolism of stressful conditions have already been extensively covered in http://pharmaceuticalintelligence.com  in the following and more:

The Amazing Structure and Adaptive Functioning of the Kidneys: Nitric Oxide – Part I

http://pharmaceuticalintelligence.com/2012/11/26/the-amazing-structure-and-adaptive-functioning-of-the-kidneys/

Nitric Oxide and iNOS have Key Roles in Kidney Diseases – Part II

http://pharmaceuticalintelligence.com/2012/11/26/nitric-oxide-and-inos-have-key-roles-in-kidney-diseases/

The pathological role of IL-18Rα in renal ischemia/reperfusion injury – Nature.com

http://pharmaceuticalintelligence.com/2014/10/24/the-pathological-role-of-il-18r%CE%B1-in-renal-ischemiareperfusion-injury-nature-com/

Summary, Metabolic Pathways

http://pharmaceuticalintelligence.com/2014/10/23/summary-metabolic-pathways/

 

Read Full Post »

Depth Underwater and Underground

Writer and Curator: Larry H. Bernstein, MD, FCAP 

 

Introduction

Deep diving for mammals is dangerous for humans and land based animals for too long, and it has dangerous consequences, most notable in nitrogen emboli  with very deep underwater diving. Other mammals live in water and have adapted to a water habitat.  This is another topic that needs further exploration.

Deep diving has different meanings depending on the context. Even in recreational diving the meaning may vary:

In recreational diving, a depth below about 30 metres (98 ft), where nitrogen narcosis becomes a significant hazard for most divers, may be considered a “deep dive”

In technical diving, a depth below about 60 metres (200 ft) where hypoxic breathing gas becomes necessary to avoid oxygen toxicity may be considered a “deep dive”.

Early experiments carried out by Comex S.A. (Compagnie maritime d’expertises) using hydrox and trimix attained far greater depths than any recreational technical diving. One example being the Comex Janus IV open-sea dive to 501 metres (1,644 ft) in 1977. The open-sea diving depth record was achieved in 1988 by a team of Comex divers who performed pipe line connection exercises at a depth of 534 metres (1,752 ft) in the Mediterranean Sea as part of the Hydra 8 program. These divers needed to breathe special gas mixtures because they were exposed to very high ambient pressure (more than 50 times atmospheric pressure).

Then there is the adaptation to the water habitat as a living environment. The two main types of aquatic ecosystems are marine ecosystems and freshwater ecosystems.

http://en.wikipedia.org/wiki/Deep_diving

Marine ecosystems are part of the earth’s aquatic ecosystem. The habitats that make up this vast system range from the productive nearshore regions to the barren ocean floor. The marine waters may be fully saline, brackish or nearly fresh. The saline waters have a salinity of 35-50 ppt (= parts per thousand). The freshwater has a salinity of less than 0.5 ppt. The brackish water lies in between these 2. Marine habitats are situated from the coasts, over the continental shelf to the open ocean and deep sea. The ecosystems are sometimes linked with each other and are sometimes replacing each other in other geographical regions. The reason why habitats differ from another is because of the physical factors that influence the functioning and diversity of the habitats. These factors are temperature, salinity, tides, currents, wind, wave action, light and substrate.

Marine ecosystems are home to a host of different species ranging from planktonic organisms that form the base of the marine food web to large marine mammals. Many species rely on marine ecosystems for both food and shelter from predators. They are very important to the overall health of both marine and terrestrial environments. Coastal habitats are those above the spring high tide limit or above the mean water level in non-tidal waters.  They are close to the sea and include habitats such as coastal dunes and sandy shores, beaches , cliffs and supralittoral habitats. Coastal habitats alone account for approximately 30% of all marine biological productivity.

http://www.marbef.org/wiki/marine_habitats_and_ecosystems

All plant and animal life forms are included from the microscopic picoplankton all the way to the majestic blue whale, the largest creature in the sea—and for that matter in the world. It wasn’t until the writings of Aristotle from 384-322 BC that specific references to marine life were recorded. Aristotle identified a variety of species including crustaceans, echinoderms, mollusks, and fish.
Today’s classification system was developed by Carl Linnaeus external link as an important tool for use in the study of biology and for use in the protection of biodiversity. Without very specific classification information and a naming system to identify species’ relationships, scientists would be limited in attempts to accurately describe the relationships among species. Understanding these relationships helps predict how ecosystems can be altered by human or natural factors.

Preserving biodiversity is facilitated by taxonomy. Species data can be better analyzed to determine the number of different species in a community and to determine how they might be affected by environmental stresses. Family, or phylogenetic, trees for species help predict environmental impacts on individual species and their relatives.

http://marinebio.org/oceans/marine-taxonomy/

For generations, whales and other marine mammals have intrigued humans. 2,400 years ago, Aristotle, a Greek scientist and philosopher, recognized that whales are mammals, not fish, because they nurse their young and breathe air like other mammals. There are numerous myths and legends surrounding marine mammals. The Greeks believed that killing a dolphin was as bad as murdering a human. An Amazon legend said that river dolphins came to shore dressed as men to woo pretty girls during fiestas. During the Middle Ages, there were numerous legends surrounding the narwhals’ amazing tusk, which was thought to have come from the unicorn.

Insert movie

Marine mammals evolved from their land dwelling ancestors over time by developing adaptations to life in the water. To aid swimming, the body has become streamlined and the number of body projections has been reduced. The ears have shrunk to small holes in size and shape. Mammary glands and sex organs are not part of the external physiology, and posterior (hind) limbs are no longer present.

Mechanisms to prevent heat loss have also been developed. The cylindrical body shape with small appendages reduces the surface area to volume ratio of the body, which reduces heat loss. Marine mammals also have a counter current heat exchange mechanism created by convergent evolution external link where the heat from the arteries is transferred to the veins as they pass each other before getting to extremities, thus reducing heat loss. Some marine mammals also have a thick layer of fur with a water repellent undercoat and/or a thick layer of blubber that can’t be compressed. The blubber provides insulation, a food reserve, and aids with buoyancy. These heat loss adaptations can also lead to overheating for animals that spend time out of the water. To prevent overheating, seals or sea lions will swim close to the surface with their front flippers waving in the air. They also flick sand onto themselves to keep the sun from directly hitting their skin. Blood vessels can also be expanded to act as a sort of radiator.

One of the major behavioral adaptations of marine mammals is their ability to swim and dive. Pinnipeds swim by paddling their flippers while sirenians and cetaceans move their tails or flukes up and down.

Some marine mammals can swim at relatively high speeds. Sea lions swim up to 35 kph and orcas can reach 50 kph. The fastest marine mammal, however, is the common dolphin, which reaches speeds up to 64 kph. While swimming, these animals take very quick breaths. For example, fin whales can empty and refill their huge lungs in less than 2 seconds. Marine mammals’ larynx and esophagus close automatically when they open their mouths to catch prey during dives. Oxygen is stored in hemoglobin in the blood and in myoglobin in the muscles. The lungs are also collapsible so that air is pushed into the windpipe preventing excess nitrogen from being absorbed into the tissues. Decreasing pressure can cause excess nitrogen to expand in the tissues as animals ascend to shallower depths, which can lead to decompression sickness,  aka “the bends.” Bradycardia, the reduction of heart rate by 10 to 20%, also takes place to aid with slowing respiration during dives and the blood flow to non-essential body parts. These adaptations allow sea otters to stay submerged for 4 to 5 minutes and dive to depths up to 55 m. Pinnipeds can often stay down for 30 minutes and reach average depths of 150-250 m. One marine mammal with exceptional diving skills is the Weddell seal, which can stay submerged for at least 73 minutes at a time at depths up to 600 m. The length and depth of whale dives depends on the species. Baleen whales feed on plankton near the surface of the water and have no need to dive deeply so they are rarely seen diving deeper than 100 m external link. Toothed whales seek larger prey at deeper depths and some can stay down for hours at depths of up to 2,250 m external link.

http://marinebio.org/oceans/marine-mammals/

Human Experience

Albert Behnke: Nitrogen Narcosis

Casey A. Grover and David H. Grover
The Journal of Emergency Medicine, 2014; 46(2):225–227
http://dx.doi.org/10.1016/j.jemermed.2013.08.080

As early as 1826, divers diving to great depths noted that descent often resulted in a phenomenon of intoxication and euphoria. In 1935, Albert Behnke discovered nitrogen as the cause of this clinical syndrome, a condition now known as nitrogen narcosis. Nitrogen narcosis consists of the development of euphoria, a false sense of security, and impaired judgment upon underwater descent using compressed air below 34 atmospheres (99 to 132 feet). At greater depths, symptoms can progress to loss of consciousness. The syndrome remains relatively unchanged in modern diving when compressed air is used. Behnke’s use of non-nitrogencontaining gas mixtures subsequent to his discovery during the 1939 rescue of the wrecked submarine USS Squalus pioneered the use of non-nitrogencontaining gas mixtures, which are used by modern divers when working at great depth to avoid the effects of nitrogen narcosis.

Behnke’s first duty station as a licensed physician was as assistant medical officer for Submarine Division 20 in San Diego, which was then commanded by one of the Navy’s rising stars, Captain Chester W. Nimitz of World War II fame.
In this setting, Dr. Behnke spent his free time constructively by learning to dive, using the traditional ‘‘hard-hat’’ gear aboard the USS Ortalon, a submarine rescue vessel to which he also rotated. Diving was not a notable specialty of the Navy at the time, and the service was slow in developing the infrastructure for it. Dr. Behnke devoted his efforts to research on the topic of diving medicine, as well as developing a more sound understanding of the biophysics of diving. In 1932, he wrote a letter to the Surgeon General describing some of his observations on arterial gas embolism, which earned him some accolades from the Navy and resulted in his transfer to Harvard’s School of Public Health as a graduate fellow. After 2 years at Harvard, the Navy assigned duty to Dr. Behnke at the Navy’s submarine escape training tower at Pearl Harbor. He worked extensively here on developing techniques for rescuing personnel from disabled submarines on the sea floor. In 1937, he was one of three Navy physicians assigned to the Navy’s Experimental Diving Unit. This team worked on improving the rescue system, plus updating the diving recompression tables originally developed by the British in 1908.

The intoxicating effects of diving were first described by a French physician named Colladon in 1826, who reported that descent in a diving bell resulted in his feeling a ‘‘state of excitement as though I had drunk some alcoholic liquor’’.
The etiology of this phenomenon remained largely unknown until the 1930s, when the British military researcher Damant again highlighted the issue, and reported very unpredictable behavior in his divers during descents as deep as 320 feet during the British Admiralty Deep Sea diving trials. Two initial theories arose as to the etiology for this effect, the first being from psychological causes by Hill and Phillip in 1932, and the second being from oxygen toxicity by Haldane in 1935.

Dr. Behnke and his colleagues at the Harvard School of Public Health had another idea as to the etiology of this phenomenon. In 1935, based on observation of individuals in experiments with a pressure chamber, Dr. Behnke published an article in the American Journal of Physiology in which he posited that nitrogen was the etiology of the intoxicating effects of diving.

Nitrogen narcosis, described as ‘‘rapture of the deep’’ by Jacques Cousteau, still remains a relatively common occurrence in modern diving, despite major advances in diving technology since Behnke’s initial description of the pathophysiologic cause of the condition in 1935. The development of symptoms of this condition varies from diver to diver, but usually begins when a depth of 4 atmospheres (132 feet) is reached in divers using compressed air. More sensitive divers can develop symptoms at only 3 atmospheres (99 feet), and other divers may not be affected up to depths as high as 6 atmospheres (198 feet). Interestingly, tolerance to nitrogen narcosis can be developed by frequent diving and exposure to the effects of compressed air at depth.

  1. Acott C. A brief history of diving and decompression illness. SPUMS J 1999;29:98–109.
    2. Bornmann R. Dr. Behnke, founder of UHMS, dies. Pressure 1992; 21:14.
    3. Behnke AR, Thomson RM, Motley P. The psychologic effects from breathing air at 4 atmospheric pressures. Am J Physiol 1935; 112:554–8.
    4. Behnke AR, Johnson FS, Poppen JR, Motley P. The effect of oxygen on man at pressures from 1 to 4 atmospheres. AmJ Physiol 1934; 110:565–72.

Exhaled nitric oxide concentration and decompression-induced bubble formation: An index of decompression severity in humans?

J.-M. Pontier, Buzzacott, J. Nastorg, A.T. Dinh-Xuan, K. Lambrechts
Nitric Oxide 39 (2014) 29–34
http://dx.doi.org/10.1016/j.niox.2014.04.005

Introduction: Previous studies have highlighted a decreased exhaled nitric oxide concentration (FE NO) in divers after hyperbaric exposure in a dry chamber or following a wet dive. The underlying mechanisms of this decrease remain however unknown. The aim of this study was to quantify the separate effects of submersion, hyperbaric hyperoxia exposure and decompression-induced bubble formation on FE NO after a wet dive.
Methods: Healthy experienced divers (n = 31) were assigned to either

  • a group making a scuba-air dive (Air dive),
  • a group with a shallow oxygen dive protocol (Oxygen dive) or

a group making a deep dive breathing a trimix gas mixture (deep-dive).
Bubble signals were graded with the KISS score. Before and after each dive FE NO values were measured using a hand-held electrochemical analyzer.
Results: There was no change in post-dive values of FE NO values (expressed in ppb = parts per billion) in the Air dive group (15.1 ± 3.6 ppb vs. 14.3 ± 4.7 ppb, n = 9, p = 0.32). There was a significant decrease in post-dive values of FE NO in the Oxygen dive group (15.6 ± 6 ppb vs. 11.7 ± 4.7 ppb, n = 9, p = 0.009). There was an even more pronounced decrease in the deep dive group (16.4 ± 6.6 ppb vs. 9.4 ± 3.5 ppb, n = 13, p < 0.001) and a significant correlation between KISS bubble score >0 (n = 13) and percentage decrease in post-dive FE NO values (r = -0.53, p = 0.03). Discussion: Submersion and hyperbaric hyperoxia exposure cannot account entirely for these results suggesting the possibility that, in combination, one effect magnifies the other. A main finding of the present study is a significant relationship between reduction in exhaled NO concentration and dive-induced bubble formation. We postulate that exhaled NO concentration could be a useful index of decompression severity in healthy human divers.

Brain Damage in Commercial Breath-Hold Divers

Kiyotaka Kohshi, H Tamaki, F Lemaıtre, T Okudera, T Ishitake, PJ Denoble
PLoS ONE 9(8): e105006 http://dx.doi.org:/10.1371/journal.pone.0105006

Background: Acute decompression illness (DCI) involving the brain (Cerebral DCI) is one of the most serious forms of diving related injuries which may leave residual brain damage. Cerebral DCI occurs in compressed air and in breath-hold divers, likewise. We conducted this study to investigate whether long-term breath-hold divers who may be exposed to repeated symptomatic and asymptomatic brain injuries, show brain damage on magnetic resonance imaging (MRI).
Subjects and Methods: Our study subjects were 12 commercial breath-hold divers (Ama) with long histories of diving work in a district of Japan. We obtained information on their diving practices and the presence or absence of medical problems, especially DCI events. All participants were examined with MRI to determine the prevalence of brain lesions.
Results: Out of 12 Ama divers (mean age: 54.965.1 years), four had histories of cerebral DCI events, and 11 divers demonstrated ischemic lesions of the brain on MRI studies. The lesions were situated in the cortical and/or subcortical area (9 cases), white matters (4 cases), the basal ganglia (4 cases), and the thalamus (1 case). Subdural fluid collections were seen in 2 cases. Conclusion: These results suggest that commercial breath-hold divers are at a risk of clinical or subclinical brain injury which may affect the long-term neuropsychological health of divers.

Decompression illness

Richard D Vann, Frank K Butler, Simon J Mitchell, Richard E Moon
Lancet 2010; 377: 153–64

Decompression illness is caused by intravascular or extravascular bubbles that are formed as a result of reduction in environmental pressure (decompression). The term covers both arterial gas embolism, in which alveolar gas or venous gas emboli (via cardiac shunts or via pulmonary vessels) are introduced into the arterial circulation, and decompression sickness, which is caused by in-situ bubble formation from dissolved inert gas. Both syndromes can occur in divers, compressed air workers, aviators, and astronauts, but arterial gas embolism also arises from iatrogenic causes unrelated to decompression. Risk of decompression illness is
affected by immersion, exercise, and heat or cold. Manifestations range from itching and minor pain to neurological symptoms, cardiac collapse, and death. First aid treatment is 100% oxygen and definitive treatment is recompression to increased pressure, breathing 100% oxygen. Adjunctive treatment, including fluid administration and prophylaxis against venous thromboembolism in paralyzed patients, is also recommended. Treatment is, in most cases, effective although residual deficits can remain in serious cases, even after several recompressions.

Bubbles can have mechanical, embolic, and biochemical effects with manifestations ranging from trivial to fatal. Clinical manifestations can be caused by direct effects from extravascular (autochthonous) bubbles such as mechanical distortion of tissues causing pain, or vascular obstruction causing stroke-like signs and symptoms. Secondary effects can cause delayed symptom onset up to 24 h after surfacing. Endothelial damage by intravascular bubbles can cause capillary leak, extravasation of plasma, and haemoconcentration. Impaired endothelial function, as measured by decreased effects of vasoactive compounds, has been reported in animals and might occur in man. Hypotension can occur in severe cases. Other effects include platelet activation and deposition, leucocyte-endothelial adhesion, and possibly consequences of vascular occlusion believed to occur in thromboembolic stroke such as ischaemia-reperfusion injury, and apoptosis.

Classification of initial and of all eventual manifestations of decompression illness in 2346 recreational diving accidents reported to the Divers Alert Network from 1998 to 2004 For all instances of pain, 58% consisted of joint pain, 35% muscle pain, and 7% girdle pain. Girdle pain often portends spinal cord involvement. Constitutional symptoms included headache, lightheadedness, inappropriate fatigue, malaise, nausea or vomiting, and anorexia. Muscular discomfort included stiffness, pressure, cramps, and spasm but excluded pain. Pulmonary manifestations included dyspnoea and cough.

Other than depth and time, risk of decompression sickness is affected by other factors that affect inert gas exchange and bubble formation, such as immersion (vs dry hyperbaric chamber exposure), exercise, and temperature. Immersion decreases venous pooling and increases venous return and cardiac output. Warm environments improve peripheral perfusion by promoting vasodilation, whereas cool temperatures decrease perfusion through vasoconstriction. Exercise increases both peripheral perfusion and temperature. The effect of environmental conditions on risk of decompression sickness is dependent on the phase of the pressure exposure. Pressure, exercise, immersion, or a hot environment increase inert gas uptake and risk of decompression sickness. During decom-pression these factors increase inert gas elimination and therefore decrease the risk of decompression sickness. Conversely, uptake is reduced during rest or in a cold environment, hence a diver resting in a cold environment on the bottom has decreased risk of decompression sickness. Rest or low temperatures during decompression increase the risk. If exercise occurs after decompression when super-saturation is present, bubble formation increases and risk of decompression sickness rises.

Exercise at specific times before a dive can decrease the risk of serious decompression sickness in animals and incidence of venous gas emboli in both animals and man. The mechanisms of these effects are unknown but might involve modulation of nitric oxide production and effects on endothelium. Venous gas emboli and risk of decompression sickness increase slightly with age and body-mass index.

Arterial gas embolism should be suspected if a diver has a new onset of altered consciousness, confusion, focal cortical signs, or seizure during ascent or within a few minutes after surfacing from a compressed gas dive.

If the diver spends much time at depth and might have absorbed substantial inert gas before surfacing, arterial gas embolism and serious decompression sickness can coexist, and in such cases, spinal cord manifestations can predominate. Other organ systems, such as the heart, can also be affected, but the clinical diagnosis of gas embolism is not reliable without CNS manifestations. Arterial gas embolism is rare in altitude exposure; if cerebral symptoms occur after altitude exposure, the cause is usually decompression sickness.

Nondermatomal hypoaesthesia and truncal ataxia are common in neurological decompression sickness and can be missed by cursory examination. Pertinent information includes level of consciousness and mental status, cranial nerve function, and motor strength. Coordination can be affected disproportionately, and abnormalities can be detected by assessment of finger-nose movement, and, with eyes open and closed, ability to stand and walk and do heel-toe walking backwards and forwards. Many of these simple tests can be done on the scene by untrained companions.

Panel: Differential diagnosis of decompression illness
Inner-ear barotrauma
Middle-ear or maxillary sinus overinfl ation
Contaminated diving gas and oxygen toxic effects
Musculoskeletal strains or trauma sustained before, during, or after diving
Seafood toxin ingestion (ciguatera, pufferfish, paralytic shellfish poisoning)
Immersion pulmonary edema
Water aspiration
Decompression chamber

Decompression chamber

Decompression chamber. fluidic or pneumatic ventilator is shown at the left. The infusion pump is contained within a plastic cover, in which 100% nitrogen is used to decrease the fi re risk in the event of an electrical problem. The monitor screen is outside the chamber and can be seen through the viewing port. Photo from Duke University Medical Center, with permission.

Long-term outcomes of 69 divers with spinal cord decompressionsickness, by manifestation
n %
No residual symptoms 34 49·3
Any residual symptom 35 50·7
Mild paraesthesias, weakness, or pain 14 20·3
Some impairment of daily activities 21 30·4
Difficulty walking 11 15·9
Impaired micturition 13 18·8
Impaired defecation 15 21·7
Impaired sexual function 15 21·7

Decompression illness occurs in a small population but is an international problem that few physicians are trained to recognise or manage. Although its manifestations are often mild, the potential for permanent injury exists in severe cases, especially if unrecognised or inadequately treated. Emergency medical personnel should be aware of manifestations of decompression illness in the setting of a patient with a history of recent diving or other exposure to substantial pressure change, and should contact an appropriate consultation service for advice.

Diving Medicine: Contemporary Topics and Their Controversies

Michael B. Strauss and Robert C. Borer, Jr
Am J Emerg Med 2001; 19:232-238
http://dx.doi.org:/10.1053/ajem.2001.22654

SCUBA diving is a popular recreational sport. Although serious injuries occur infrequently, when they do knowledge of diving medicine and/or where to obtain appropriate consultation is essential. The emergency physician is likely to be the first physician contact the injured diver has. We discuss 8 subjects
in diving medicine which are contemporary, yet may have controversies associated with them. From this information the physician dealing primarily with the injured diver will have a basis for understanding and managing, as
well as where to find additional help, for his/her patients’ diving injuries.

Over the past 10 years, new knowledge and equipment improvements have made diving safer and more enjoyable. Estimates of actively participating sports divers show a striking increase over this time interval while the number of SCUBA diving deaths annually has remained nearly level at approximately 100. A further indicator of recreational diving safety is that reflected in the nearly constant number of diving injuries (1000 per annum) over the most recent 5 reported years, or approximately 0.53 to 3.4 incidents/10,000 dives.

Divers Alert Network.
The Divers Alert Network (DAN) is a nonprofit organization directed and staffed by experts in the specialty of diving medicine.6 DAN provides immediate consultation for both divers and physicians in the diagnosis and initial management of diving injuries. This 24-hour service is available free world-wide through a dedicated emergency telephone line: 1-919-684-4326. The DAN staff will also identify the nearest appropriate recompression treatment facility and knowledgeable physicians for an expedient referral. General diving medical inquiries can be answered during normal weekday hours either through an information telephone line: 1-919-684-2948 or through an interactive web site http://www.diversalertnetwork.org.

Use of 100% Oxygen for Initial, on the Scene, Management of Diving Accidents
The breathing of pure oxygen is crucial for the initial management of the diving related problems of arterial gas embolism (AGE), decompression sickness (DCS), pulmonary barotrauma (thoracic squeeze), aspiration pneumonitis, and hypoxic encephalopathy associated with near drowning. In 1985, Dick reported that in many cases the neurologic symptoms of AGE and DCS were resolved with the immediate breathing of pure oxygen on the surface. The breathing of pure oxygen reduces bubble size by increasing the differential pressure for the inert gas to diffuse out of the bubble and it also speeds the washout of inert gas from body tissues. The early elimination of the bubble prevents hypoxia and the interaction of the bubble with the blood vessel lining. This interaction leads to secondary problems of capillary leak, bleeding, inflammation, ischemia, and cell death. These secondary problems are the reasons not all DCS symptoms resolve with recompression chamber treatment. The immediate use of pure oxygen for the medical management of these diving problems is analogous to the use of cardiopulmonary resuscitation for the witnessed cardiac arrest; the sooner initiated the better the results.

Diving Education

Medical Fitness for Diving

Asthma has the potential risk for AGE. Neuman reviewed the subject of asthma and diving. He and his coauthors recommend that asthmatics who are asymptomatic, not on medications and have no exercised induced abnormality on pulmonary function studies be allowed to dive.

Conditions leading to loss of consciousness, such as insulin dependent diabetes and epilepsy, can result in drowning. Carefully controlled diving studies in diabetics, who are free from complications, are now defining the safe requirements for diving. Epilepsy remains as a disqualification except in individuals with a history of febrile seizures ending prior to 5 years of age.

Availability of Hyperbaric Oxygen Treatment Facilities

The availability of these chambers makes it possible for divers who become symptomatic after SCUBA diving to readily receive recompression treatment. This is important because the closer the initiation of recompression treatment to the onset of DCS (and AGE) signs and symptoms, the greater the likelihood of full recovery.

Improved Diving Equipment

Mixed and Rebreather Gas Diving
Mixed gas diving involves changing the breathing gas from air which has 20% oxygen to higher oxygen percentages (nitrox). As the amount of oxygen is increased in the gas mixture, the amount of the inert gas (nitrogen) is reduced. With oxygen enriched air there is less tissue deposition of inert gas per unit of time under water for any given depth. However, because of increased oxygen partial pressures, the seizure threshold for oxygen toxicity is lowered. For normal sports diving activities, oxygen toxicity with mixed gas diving is only a theoretical concern.

Decompression Illness is More Than Bubbles

When AGE occurs, DCS symptoms may be concurrent or appear during or after recompression treatment even though the decompression tables were not violated on the dive. When DCS occurs in this situation it appears resistant to recompression treatment (Neuman) perhaps because of the inflammatory reaction generated by the bubble-blood vessel interaction from the AGE. In cases of DCI where components of both DCS and AGE are suspected, the diver should be observed for a minimum of 24 hours after the recompression treatment is completed for the delayed onset of DCS.

No theory of DCS discounts the primary role of bubbles in this condition. However, new information suggests that there are precursors to bubble formation and post-bubbling events that occur as a consequence of the bubbles. As mentioned earlier, venous gas emboli are a common occurrence diving ascent and ordinarily are filtered out harmlessly by the lungs. Precursors to DCS include stasis, dehydration and too rapid of ascents. These conditions allow the ubiquitous VGE to enlarge, coalesce and occlude the venous side of the circulation. Massive venous bubbling to the lungs can cause pulmonary vessel obstruction described as the chokes. If right to left shunts occur in the heart, VGE can become AGE to the brain. If the arterial flow is slow enough and/or the gradients large enough, autochthonus (ie, spontaneous) bubbles can form in the arterial circulation and lead to any of the consequences of AGE. In such situations it could be difficult to determine whether the DCI event was from AGE or DCS even after careful analysis of the dive profile. Hollenbeck’s model for diving paraplegia includes the setting of venous stasis (Batson’s plexus of veins) in the spinal canal, bubble formation, bubble enlargement possibly from off gassing of the spinal cord, blood vessel occlusion, and venous side infarctions of the spinal cord.
Contemporary Management of DCS

Problem Intervention Effect
Bubble Recompression
with HBO
Reduce bubble size
1. Washout inert gas.
2. Change bubble composition by diffusion.
Stasis and dehydration Hydration: oral fluids if alert, IV fluids otherwise. Improve blood flow.
InflammationCell Ischemia ? Anti-inflammatory medicationsHBO Reduce interaction between bubble and blood vessel endothelium.
Improve oxygen availability to hypoxic tissues, reduce edema and also reduces the interaction between bubble and blood vessel endothelium.

.

Conclusions

We anticipate that in the future there will be further improvements for the safety and enjoyment of the recreational SCUBA diver. For example, the dive computer of the future will be able to individualize dive profiles for different personal medical parameters such as age, body composition and fitness level. Diver locators could quickly target a missing diver and save time and gas consumption as well as prevent serious diving mishaps. Drugs may be developed that would minimize the effect of bubbles interacting with body tissues and prevent DCS and AGE.

Extracorporeal membrane oxygenation therapy for pulmonary decompression illness

Yutaka Kondo, Masataka Fukami and Ichiro Kukita
Kondo et al. Critical Care 2014; 18:438 http://ccforum.com/content/18/3/438/10.1186/cc13935

Pulmonary decompression illness is rarely observed in clinical settings, and most patients die prior to hospitalization. We administered ECMO therapy to rescue a patient, even though this therapy has rarely been reported with good outcome in patients with decompression illness. In addition, we had to select venovenous ECMO even with the patient showing right ventricular failure. A lot of physicians may select venoarterial ECMO if the patient shows right ventricular failure, but the important physiological mechanism of pulmonary decompression illness is massive air embolism in the pulmonary arteries, and the bubbles diminish within the first 24 hours. The management of decompression illness therefore differs substantially from the usual right-sided heart failure.

Extremes of barometric pressure

Jane E Risdall, David P Gradwell
Anaesthesia and Intensive Care Medicine 16:2
Ascent to elevated altitude, commonly achieved through flight, by climbing or by residence in highland regions, exposes the individual to reduced ambient pressure. Although there are physical manifestations of this exposure as a consequence of Boyle’s law, the primary physiological challenge is of hypobaric hypoxia. The acute physiological and longer-term adaptive responses of the cardiovascular, respiratory, hematological and neurological systems to altitude are described, together with an outline of the presentation and management of acute mountain sickness, high-altitude pulmonary edema and high-altitude cerebral edema. While many millions experience modest exposure to altitude as a result of flight in pressurized aircraft, fewer individuals are exposed to increased ambient pressure. The pressure changes during diving and hyperbaric exposures result in greater changes in gas load and gas toxicity. Physiological effects include the consequences of increased work of breathing and redistribution of circulating volume. Neurological manifestations may be the direct result of pressure or a consequence of gas toxicity at depth. Increased tissue gas loads may result in decompression illness on return to surface or subsequent ascent in flight.

  • understand the physical effects of changes in ambient pressure and the physiological consequences on the cardiovascular respiratory and neurological systems
  • gain an awareness that exposure to reduced ambient pressure produces both acute and more chronic effects, with differing signs, symptoms and time to onset at various altitudes
  • develop an awareness of the toxic effects of ‘inert’ gases at increased ambient pressures and the pathogenesis and management of decompression illness

Decompression illness According to Henry’s law, at a constant temperature the amount of gas which dissolves in a liquid is proportional to the pressure of that gas or its partial pressure, if it is part of a mixture of gases. Breathing gases at increased ambient pressure will increase the amount of each gas dissolved in the fluid phases of body tissues. On ascent this excess gas has to be given up. If the ascent is controlled at a sufficiently slow rate, elimination will be via the respiratory system. If the ascent is too fast, excess gas may come out of solution and form free bubbles in the tissues or circulation. Bubbles may contain any of the gases in the breathing mixture, but it is the presence of inert gas bubbles (nitrogen or helium) that are thought most likely to give rise to problems, since the elimination of excess oxygen is achieved by metabolism as well as ventilation. These bubbles may act as venous emboli or may trigger inflammatory tissue responses giving rise to symptoms of decompression illness (DCI). Signs and symptoms of DCI may appear up to 48 hours after exposure to increased ambient pressure and include joint pains, motor and sensory deficits, dyspnoea, cough and skin rashes.

Neurological effects of deep diving

Marit Grønning, Johan A. Aarli
Journal of the Neurological Sciences 304 (2011) 17–21
http://dx.doi.org:/10.1016/j.jns.2011.01.021

Deep diving is defined as diving to depths more than 50 m of seawater (msw), and is mainly used for occupational and military purposes. A deep dive is characterized by the compression phase, the bottom time and the decompression phase. Neurological and neurophysiologic effects are demonstrated in divers during the compression phase and the bottom time. Immediate and transient neurological effects after deep dives have been shown in some divers. However, the results from the epidemiological studies regarding long term neurological effects from deep diving are conflicting and still not conclusive.

Possible immediate neurological effects of deep diving
Syndrome Pressure
Hyperoxia/oxygen seizures >152 kPa (5 msw)
HypoxiaHypercapnia
Nitrogen narcosis >354 kPa (25 msw)
High pressure nervous syndrome >1.6 MPa (150 msw)
Neurological decompression sickness

Neurological effects have been demonstrated, both clinically and neurophysiologically in divers during the compression phase and the bottom time. Studies of divers before and after deep dives have shown immediate and transient neurological effects in some divers. However, the results from the epidemiological and clinical studies regarding long term neurological effects from deep diving are conflicting and still not conclusive. Prospective clinical studies with sufficient power and sensitivity are needed to solve this important issue.

Today deep diving to more than 100 msw is routinely performed globally in the oil- and gas industry. In the North Sea remote underwater intervention and maintenance is performed by the use of remotely operated vehicles (ROV), both in conjunction to and as an alternative to manned underwater operations. There will, however, always be a need for human divers in the technically more advanced underwater operations and for contingency repair operations.

P300 latency indexes nitrogen narcosis

Barry Fowler, Janice Pogue and Gerry Porlier
Electroencephalography, and clinical Neurophysiology, 1990, 75:221-229

This experiment investigated the effects of nitrogen narcosis on reaction time (RT) and P300 latency and amplitude, Ten subjects breathed either air or a non-narcotic 20% oxygen-80% helium (heliox) mixture in a hyperbaric chamber at 6.5, 8.3 and 10 atmospheres absolute (ATA), The subjects responded under controlled accuracy conditions to visually presented male or female names in an oddball paradigm. Single-trial analysis revealed a strong relationship between RT and P300 latency, both of which were slowed in a dose-related manner by hyperbaric air but not by heliox. A clear-cut dose-response relationship could not be established for P300 amplitude. These results indicate that P300 latency indexes nitrogen narcosis and are interpreted as support for the slowed processing model of inert gas narcosis.

Adaptation to Deep Water Habitat

Effects of hypoxia on ionic regulation, glycogen utilization and antioxidative ability in the gills and liver of the aquatic air-breathing fish Trichogaster microlepis

Chun-Yen Huang, Hui-Chen Lina, Cheng-Huang Lin
Comparative Biochemistry and Physiology, Part A 179 (2015) 25–34
http://dx.doi.org/10.1016/j.cbpa.2014.09.001

We examined the hypothesis that Trichogaster microlepis, a fish with an accessory air-breathing organ, uses a compensatory strategy involving changes in both behavior and protein levels to enhance its gas exchange ability. This compensatory strategy enables the gill ion-regulatory metabolism to maintain homeostasis during exposure to hypoxia. The present study aimed to determinewhether ionic regulation, glycogen utilization and antioxidant activity differ in terms of expression under hypoxic stresses; fish were sampled after being subjected to 3 or 12 h of hypoxia and 12 h of recovery under normoxia. The air-breathing behavior of the fish increased under hypoxia. No morphological modification of the gills was observed. The expression of carbonic anhydrase II did not vary among the treatments. The Na+/K+-ATPase enzyme activity did not decrease, but increases in Na+/K+-ATPase protein expression and ionocyte levels were observed. The glycogen utilization increased under hypoxia as measured by glycogen phosphorylase protein expression and blood glucose level, whereas the glycogen content decreased. The enzyme activity of several components of the antioxidant system in the gills, including catalase, glutathione peroxidase, and superoxidase dismutase, increased in enzyme activity. Based on the above data, we concluded that T. microlepis is a hypoxia-tolerant species that does not exhibit ion-regulatory suppression but uses glycogen to maintain energy utilization in the gills under hypoxic stress. Components of the antioxidant system showed increased expression under the applied experimental treatments.

Divergence date estimation and a comprehensive molecular tree of extant cetaceans

Michael R. McGowen , Michelle Spaulding, John Gatesy
Molecular Phylogenetics and Evolution 53 (2009) 891–906
http://dx.doi.org/10.1016/j.ympev.2009.08.018

Cetaceans are remarkable among mammals for their numerous adaptations to an entirely aquatic existence, yet many aspects of their phylogeny remain unresolved. Here we merged 37 new sequences from the nuclear genes RAG1 and PRM1 with most published molecular data for the group (45 nuclear loci, transposons, mitochondrial genomes), and generated a supermatrix consisting of 42,335 characters. The great majority of these data have never been combined. Model-based analyses of the supermatrix produced a solid, consistent phylogenetic hypothesis for 87 cetacean species. Bayesian analyses corroborated odontocete (toothed whale) monophyly, stabilized basal odontocete relationships, and completely resolved branching events within Mysticeti (baleen whales) as well as the problematic speciose clade Delphinidae (oceanic dolphins). Only limited conflicts relative to maximum likelihood results were recorded, and discrepancies found in parsimony trees were very weakly supported. We utilized the Bayesian supermatrix tree to estimate divergence dates among lineages using relaxed-clock methods. Divergence estimates revealed rapid branching of basal odontocete lineages near the Eocene–Oligocene boundary, the antiquity of river dolphin lineages, a Late Miocene radiation of balaenopteroid mysticetes, and a recent rapid radiation of Delphinidae beginning [1]10 million years ago. Our comprehensive,  time calibrated tree provides a powerful evolutionary tool for broad-scale comparative studies of Cetacea.

Mitogenomic analyses provide new insights into cetacean origin and evolution

Ulfur Arnason, Anette Gullberg, Axel Janke
Gene 333 (2004) 27–34
http://dx.doi.org:/10.1016/j.gene.2004.02.010

The evolution of the order Cetacea (whales, dolphins, porpoises) has, for a long time, attracted the attention of evolutionary biologists. Here we examine cetacean phylogenetic relationships on the basis of analyses of complete mitochondrial genomes that represent all extant cetacean families. The results suggest that the ancestors of recent cetaceans had an explosive evolutionary radiation 30–35 million years before present. During this period, extant cetaceans divided into the two primary groups, Mysticeti (baleen whales) and Odontoceti (toothed whales). Soon after this basal split, the Odontoceti diverged into the four extant lineages, sperm whales, beaked whales, Indian river dolphins and delphinoids (iniid river dolphins, narwhals/belugas, porpoises and true dolphins). The current data set has allowed test of two recent morphological hypotheses on cetacean origin. One of these hypotheses posits that Artiodactyla and Cetacea originated from the extinct group Mesonychia, and the other that Mesonychia/Cetacea constitutes a sister group to Artiodactyla. The current results are inconsistent with both these hypotheses. The findings suggest that the claimed morphological similarities between Mesonychia and Cetacea are the result of evolutionary convergence rather than common ancestry.

The order Cetacea traditionally includes three suborders: the extinct Archaeoceti and the recent Odontoceti and Mysticeti. It is commonly believed that the evolution of ancestral cetaceans from terrestrial to marine (aquatic) life was accompanied by a fast and radical morphological adaptation. Such a scenario may explain why it was, for a long time, difficult to morphologically establish the position of Cetacea in the mammalian tree and even to settle whether Cetacea constituted a monophyletic group.

Biochemical analyses in the 1950s  and 1960s had shown a closer relationship between cetaceans and artiodactyls (even-toed hoofed mammals) than between cetaceans and any other eutherian order and karyological studies in the late 1960s and early 1970s unequivocally supported cetacean monophyly (Arnason, 1969, 1974). The nature of the relationship between cetaceans and artiodactyls was resolved in phylogenetic studies of mitochondrial (mt) cytochrome b (cytb) genes (Irwin and Arnason, 1994; Arnason and Gullberg, 1996) that placed Cetacea within the order Artiodactyla itself as the sister group of the Hippopotamidae (see also Sarich, 1993). The Hippopotamidae/ Cetacea relationship was subsequently supported in studies of nuclear data (Gatesy et al., 1996; Gatesy, 1997) and statistically established in analysis of complete mt genomes (Ursing and Arnason, 1998). The relationship has also been confirmed in analyses of combined nuclear and mt sequences (Gatesy et al., 1999; Cassens et al., 2000) and in studies of short interspersed repetitive elements (SINEs). Artiodactyla and Cetacea are now commonly referred to as Cetartiodactyla.

Previous analyses of the complete cytb gene of more than 30 cetacean species (Arnason and Gullberg, 1996) identified five primary lineages of recent cetaceans, viz., Mysticeti and the four odontocete lineages Physeteridae (sperm whales), Platanistidae (Indian river dolphins), Ziphiidae (beaked whales) and Delphinoidea (iniid river dolphins, porpoises, narwhals and dolphins). However, these studies left unresolved the relationships of the five lineages as well as those between the three delphinoid families Monodontidae (narwhals, belugas), Phocoenidae (porpoises) and Delphinidae (dolphins). Similarly, the relationships between the four mysticete families Balaenidae (right whales), Neobalaenidae (pygmy right whales), Eschrichtiidae (gray whales) and Balaenopteridae (rorquals) were not conclusively resolved in analyses of cytb genes.

Fig. (not shown). Cetartiodactyl relationships and the estimated times of their divergences. The tree was established on the basis of maximum likelihood analysis of the concatenated amino acid (aa) sequences of 12 mt protein-coding genes. Length of alignment 3610 aa. Support values for branches A–H are shown in the insert.
Cetruminantia (branch A) receives moderate support and Cetancodonta (B) strong support. Cetacea (C) splits into monophyletic Mysticeti (baleen whales) and monophyletic Odontoceti (toothed whales). Odontoceti has four basal lineages, Physeteridae (sperm whales: represented by the sperm and pygmy sperm whales), Ziphiidae (beaked whales: bottlenose and Baird’s beaked whales), Platanistidae (Indian river dolphins: Indian river dolphin) and Delphinoidea. Delphinoidea encompasses the families Iniidae (iniid river dolphins: Amazon river dolphin, La Plata dolphin), Monodontidae (narwhals/belugas: narwhal), Phocoenidae (porpoises: harbour porpoise) and Delphinidae (dolphins: white-beaked dolphin). The common odontocete branch and the branches separating the four cetacean lineages are short. These relationships are therefore somewhat unstable (cf. Section 3.1 and Table 1). Iniid river dolphins (F) are solidly nested within the Delphinoidea (E). Thus, traditional river dolphins (Platanistidae + Iniidae) do not form a monophyletic unit. Molecular estimates of divergence times (Sanderson 2002) were based on two calibration points, A/C-60 and O/M-35 (cf. Section 3.2). Due to the short lengths of internal branches, some estimates for these divergences overlap. NJ: neighbor joining; MP: maximum parsimony; LBP: local bootstrap probability; QP: quartet puzzling. The bar shows the number of aa substitutions per site.

The limited molecular resolution among basal cetacean lineages has been known for some time. Studies of hemoglobin and myoglobin (Goodman, 1989; Czelusniak et al., 1990) have either joined Physeteridae and Mysticeti to the exclusion of Delphinoidea (myoglobin data) or Mysticeti and Delphinoidea to the exclusion of Physeteridae (hemoglobin data). Thus, neither of the data sets identified monophyletic Odontoceti by joining the two odontocete lineages (Physeteridae and Delphinoidea) to the exclusion of Mysticeti. A similar instability was recognized and cautioned against in analyses of some mt data, notably, sequences of rRNA genes (Arnason et al., 1993b). The suggestion (Milinkovitch et al., 1993) of a sister group relationship between Physeteridae and the mysticete family Balaenopteridae (rorquals) was based on a myoglobin data set (which joins Physeteridae and Mysticeti to the exclusion of Delphinoidea) that was complemented with partial data of the mt 16S rRNA gene.

The cetancodont divergence times calculated using A/C-60 and O/M-35 as references have been included in Fig. 1. As a result of the short branches separating several cetacean lineages, the estimates of these divergences overlap. The same observation has been made in calculations based on SINE flanking sequences (Nikaido et al., 2001). There is a general consistency between the current and the flanking sequence datings, except for those involving the Balaenopteridae, which are somewhat younger in our analysis than in the SINEs study. The currently estimated age of the divergence between Hippopotamus and Cetacea (c53.5 MYBP) is consistent with the age (>50 MY) of the oldest archaeocete fossils identified so far (Bajpai and Gingerich, 1998). This suggests that the ages allocated to the two references, A/C-60 (the divergence between ruminant artiodactyls and cetancodonts) and O/M-35 (the divergence between odontocetes and mysticetes) are reasonably accurate.

The dating of the divergence between the blue and fin whales is of interest regarding hybridization between closely related mammalian species. Previous molecular analyses (Arnason et al., 1991b; Spilliaert et al., 1991) demonstrated the occurrence of hybridization between these two species. These studies, which were based on three hybrids (one female and two males), showed that either species could be the mother or father in these hybridizations. The two male hybrids had rudimentary testes, whereas the female hybrid was in her second pregnancy. This suggests that the blue and fin whales may be close to the limit for permissible species hybridization among mammals.

The current data set has allowed examination of the coherence between the molecular results and two prevalent morphological hypotheses related to cetacean evolution. The first hypothesis, which in essence originates from Van Valen (1966, 1968), postulates that monophyletic Artiodactyla and monophyletic Cetacea evolved separately from the extinct Palaeocene group Mesonychia. This hypothesis was recently reinforced in a morphological study (Thewissen et al., 2001) that included mesonychians, two archaeocete taxa (Ambuloocetus and Pakicetus) and some extant and fossil artiodactyls. The study of Thewissen et al. (2001) showed a sister group relationship between monophyletic Artiodactyla and monophyletic Cetacea, with Mesonychia as the basal sister group of Artiodactyla/Cetacea, a conclusion consistent with the palaeontological age of Mesonychia relative to that of Artiodactyla and Cetacea. The second hypothesis favours a sister group relationship between Mesonychia and Cetacea with the Mesonychia/Cetacea clade as the sister group of monophyletic Artiodactyla (O’Leary and Geisler, 1999; see also Gatesy and O’Leary, 2001).

Although the position of Mesonychia differs in the two morphological hypotheses, both correspond to a sister group relationship between Cetacea and monophyletic Artiodactyla among extant cetartiodactyls. Thus, both hypotheses can be tested against the current data set. The result of such a test has been included in Table 1, topology (m)(not shown). As evident, both these morphological hypotheses are incongruent with the mitogenomic findings.

Morphological studies have not provided an answer to the question whether mysticetes and odontocetes had separate origins among the archaeocetes (Fordyce and de Muizon, 2001). However, the long common cetacean branch and the short branches separating the five extant cetacean lineages strongly suggest an origin of modern cetaceans from the same archaeocete group (probably the Dorudontidae).

The limbs of Ambulocetus constitute somewhat of an evolutionary enigma. As evident in Thewissen et al.’s (1994) paper, Ambulocetus has very large hind limbs compared to its forelimbs, a difference that is less pronounced in later silhouette drawings of the animal. It is nevertheless evident that evolution from the powerful hindlimbs of Ambulocetus to their rudimentation in archaeocetes constitutes a remarkable morphological reversal if Ambulocetus is connected to the cetacean branch after the separation of the hippopotamid and cetacean lineages.

For natural reasons, systematic schemes have traditionally been based on external morphological characteristics. The rates of morphological and molecular evolution are rarely (if ever) strictly correlated, however, and this may give rise to inconsistency between traditional systematics and molecular findings. The emerging consensus that the order Cetacea resides within another traditional order, Artiodactyla, makes apparent the incongruity in cetartiodactyl nomenclature (Graur and Higgins, 1994). In this instance, a possible solution for maintaining reasonable consistency between nomenclature and phylogeny would be to recognize Cetartiodactyla as an order with three suborders: Suina, Tylopoda and Cetruminantia. According to such a scheme, Cetacea would (together with the Hippopotamidae) constitute a parvorder within the infraorder Cetancodonta.

Cytochrome b and Bayesian inference of whale phylogeny

Laura May-Collado, Ingi Agnarsson
Molecular Phylogenetics and Evolution 38 (2006) 344–354
http://dx.doi.org//10.1016/j.ympev.2005.09.019

In the mid 1990s cytochrome b and other mitochondrial DNA data reinvigorated cetacean phylogenetics by proposing many novel

and provocative hypotheses of cetacean relationships. These results sparked a revision and reanalysis of morphological datasets, and the collection of new nuclear DNA data from numerous loci. Some of the most controversial mitochondrial hypotheses have now become benchmark clades, corroborated with nuclear DNA and morphological data; others have been resolved in favor of more traditional views. That major conflicts in cetacean phylogeny are disappearing is encouraging. However, most recent papers aim specifically to resolve higher-level conflicts by adding characters, at the cost of densely sampling taxa to resolve lower-level relationships. No molecular study to date has included more than 33 cetaceans. More detailed molecular phylogenies will provide better tools for evolutionary studies. Until more genes are available for a high number of taxa, can we rely on readily available single gene mitochondrial data? Here, we estimate the phylogeny of 66 cetacean taxa and 24 outgroups based on Cytb sequences. We judge the reliability of our phylogeny based on the recovery of several deep-level benchmark clades. A Bayesian phylogenetic analysis recovered all benchmark clades and for the Wrst time supported Odontoceti monophyly based exclusively on analysis of a single mitochondrial gene. The results recover the monophyly, with the exception of only one taxa within Cetacea, and the most recently proposed super- and subfamilies. In contrast, parsimony never recovered all benchmark clades and was sensitive to a priori weighting decisions. These results provide the most detailed phylogeny of Cetacea to date and highlight the utility of both Bayesian methodology in general, and of Cytb in cetacean phylogenetics. They furthermore suggest that dense taxon sampling, like dense character sampling, can overcome problems in phylogenetic reconstruction.

Some long standing debates are all but resolved: our understanding of deeper level cetacean phylogeny has grown strong. However, the strong focus of most recent studies, aiming specifically to resolve these higher level conflicts by adding mostly characters rather than taxa, has left our understanding of lower level relationships among whale species lagging behind. Mitogenomic data, for example, is available only for 16 cetacean species, and no molecular study to date has included more than 33 cetaceans. It seems timely to focus on more detailed (genus, and species level) molecular phylogenies. These will provide better tools for detailed evolutionary studies, and are necessary to test existing morphological phylogenetic hypotheses, and current cetacean classification.

We judge the reliability of our phylogeny based on the recovery of the previously mentioned benchmark clades, in addition to the less controversial clades Perissodactyla, Euungulata (sensu Waddell et al., 2001; Perissodactyla+ Cetartiodactyla), Cetacea, and Mysticeti. Because Cytb is thought to be most reliable at lower taxonomic levels (due to high substitution rates), recovering ‘known’ deeper clades gives credibility to these new findings which have not been addressed by studies using few taxa. We compare the performance of Bayesian analyses versus parsimony under four different models, and briefly examine the sensitivity of the results to taxon sampling. We use our results to discuss agreement and remaining conflict in cetacean phylogenetics, and provide comments on current classification.

The Bayesian analysis recovered all seven benchmark clades. Support for five of the benchmark clades is high (100 posterior probabilities) but rather low for Cetancodonta (79) and marginal for the monophyly of Odontoceti. The analysis also recovered all but one family level, and most sub- and superfamily level cetacean taxa. The results broadly corroborate current cetacean classiffcation, while also pointing to some lower-level groups that may need redefinition.

Many recent cetacean phylogenetic studies include relatively few taxa, in part due to a focus on generating more characters to resolve higher level phylogenetics. While addressing crucial questions and providing the backbone for lower level phylogenies, such studies have limited utility for classification, and for comparative evolutionary studies. In some cases sparse taxon sampling may also confound the results. Of course, taxon sampling is usually simply constrained by the availability of character data, but for some reason many studies have opted to include only one, or a few outgroup taxa, even if many are available.

We find that as long as outgroup taxon sampling was extensive, Bayesian analyses of Cytb recovered all the a priori identified benchmark clades. When only a few outgroups were chosen, however, the Bayesian analysis negated Odontoceti monophyly, as have many previous parsimony analyses of mitochondrial DNA. Furthermore, in almost every detailed comparison possible our results mirror the findings O’Leary et al. (2004), the most ‘character-complete’ (but including relatively few cetacean taxa) analysis to date (37,000 characters from morphology, SINE, and 51 gene fragments). This result gives credibility to our findings, including previously untested lower level clades.

  • Monophyly and placement of Mysticeti (baleen whales).
  • Monophyly of Odontoceti (toothed whales)
  • Delphinoids
  • River Dolphins
  • Beaked and sperm whales

A major goal of phylogenetics is a phylogeny of life (i.e., many taxa), based on multiple lines of evidence (many characters of many types). However, when phylogenies based on relatively few characters can be judged reliable based on external evidence (taxonomic congruence with other phylogenies using many characters, but few taxa), they seem like very promising and useful ‘first guess’ hypotheses. The evolution of sexual dimorphism, echolocation, social behavior, and whistles and other communicative signals, and major ecological shifts (e.g., transition to fresh water) are among the numerous interesting questions in cetacean biology that this phylogeny can help answer.

Deep-diving sea lions exhibit extreme bradycardia in long duration dives

Birgitte I. McDonald1, and Paul J. Ponganis
The Journal of Experimental Biology (2014) 217, 1525-1534 http://dx.doi.org:/10.1242/jeb.098558

Heart rate and peripheral blood flow distribution are the primary determinants of the rate and pattern of oxygen store utilization and ultimately breath-hold duration in marine endotherms. Despite this, little is known about how otariids (sea lions and fur seals) regulate heart rate (fH) while diving. We investigated dive fH in five adult female California sea lions (Zalophus californianus) during foraging trips by instrumenting them with digital electrocardiogram (ECG) loggers and time depth recorders. In all dives, dive fH (number of beats/duration; 50±9 beats min−1) decreased compared with surface rates (113±5 beats min−1), with all dives exhibiting an instantaneous fH below resting (<54 beats min−1) at some point during the dive. Both dive fH and minimum instantaneous fH significantly decreased with increasing dive duration. Typical instantaneous fH profiles of deep dives (>100 m) consisted of:

(1) an initial rapid decline in fH resulting in the lowest instantaneous fH of the dive at the end of descent, often below 10 beats min−1 in dives longer than 6 min in duration;
(2) a slight increase in fH to ~10–40 beats min−1 during the bottom portion of the dive; and
(3) a gradual increase in fH during ascent with a rapid increase prior to surfacing.

Thus, fH regulation in deep-diving sea lions is not simply a progressive bradycardia. Extreme bradycardia and the presumed associated reductions in pulmonary and peripheral blood flow during late descent of deep dives should

(a) contribute to preservation of the lung oxygen store,
(b) increase dependence of muscle on the myoglobin-bound oxygen store,
(c) conserve the blood oxygen store and
(d) help limit the absorption of nitrogen at depth.

This fH profile during deep dives of sea lions may be characteristic of deep-diving marine endotherms that dive on inspiration as similar fH profiles have been recently documented in the emperor penguin, another deep diver that dives on inspiration.

The resting ƒH measured in this study (54±6 beats min−1) was lower than predicted for an animal of similar size (~80 beats min−1 for an 80 kg mammal). In part, this may be due to the fact that the sea lions were probably sleeping. The resting ƒH in our study was also lower than previous measurements in captive juvenile California sea lions (87±17 beats min−1, average mass 30 kg)  and wild Antarctic fur seals (78±5 beats min−1, body mass 30–50 kg). However, we found a significant negative relationship between mass and resting ƒH even with our small sample size of five sea lions (resting ƒH = –0.58 Mb +100.26, r2=0.81, F1,3=12.37, P=0.039). For a 30 kg sea lion, this equation predicts a resting ƒH of 83 beats min−1, which is similar to what was measured previously in juvenile sea lions, suggesting this equation may be useful in estimating resting ƒH in sea lions.

The sea lions exhibited a distinct sinus arrhythmia fluctuating between a minimum of 42±9 and a maximum of 87±12 beats min−1, comparable to the sinus arrhythmias described in other diving birds and mammals, including sea lions. The minimum instantaneous ƒH during the sinus arrhythmia was similar to the mean minimum ƒH in dives less than 3 min (37±7 beats min−1), indicating that in dives less than 3 min (estimated cADL), ƒH only decreased to levels observed during exhalation at rest. This is consistent with observations in emperor penguins and elephant seals, where it was proposed that in dives shorter than the aerobic dive limit (ADL) the reduction in ƒH is regulated by a mechanism of cardiorespiratory control similar to that governing the respiratory sinus arrhythmia, with a further reduction only occurring in dives longer than the ADL.

Fig. 3. (not shown) Instantaneous fH and dive depth profiles of a California sea lion (CSL12_2). Data are from (A) a short, shallow dive (1.3 min, 45 m), (B) a mid-duration dive (4.8 min, 239 m) and (C) a long-duration dive (8.5 min, 305 m). Minimum instantaneous fH reached 37 beats min−1 in the short dive
(A) 19 beats min−1 in the mid-duration dive
(B) and 7 beats min−1 in the long duration dive
(C) Prominent features typical of mid- and long-duration dives include

  • a surface interval tachycardia (pre- and post-dive);
  • a steady rapid decrease in fH during initial descent;
  • a gradual decline in fH towards the end of descent with the lowest fH of the dive at the end of descent;
  • a slight increase and sometimes variable fH during the bottom portion of the dive; and
  • a slow increase in fH during ascent,
  • often ending in a rapid increase just before surfacing.

We obtained the first diving ƒH data from wild sea lions on natural foraging trips, demonstrating how they regulate ƒH over a range of dive durations. Sea lions always decreased dive ƒH from surface ƒH values; however, individual sea lions exhibited different dive ƒH, accounting for a significant amount of the variation in the relationship between dive duration and ƒH (intra-individual correlation: 75–81%)). The individual differences in dive ƒH exhibited in this study suggest that different dive capacities of individual sea lions may partially account for the range of dive strategies exhibited in a previous study (Villegas-Amtmann et al., 2011). Despite the individual differences in ƒH, the pattern of the dive ƒH response was similar in all the sea lions. As predicted, sea lions only consistently displayed a true bradycardia on mid- to long- duration dives (>4 min) (Fig. 5A). Additionally, as seen in freely diving phocids, dive ƒH and minimum ƒH were negatively related to dive duration, with the longest duration dives having the lowest dive ƒH and displaying the most intense bradycardia, often below 10 beats min−1 (Fig. 5A,B).

Profiles of mean fH at 10 s intervals of dives

Profiles of mean fH at 10 s intervals of dives

Fig 4.  Profiles of mean fH at 10 s intervals of dives for (A) six duration categories and (B) five depth categories. Standard error bars are shown. Data were pooled from 461 dives performed by five sea lions. The number of dives in each category and the number of sea lions performing the dives in each category are provided in the keys.

The mild bradycardia and the dive ƒH profiles observed in the shorter duration dives (<3 min) were similar to those observed in trained juvenile California sea lions and adult Stellar sea lions, but much more intense than ƒH observed in freely diving Antarctic fur seals. Surprisingly, although dive ƒH of trained Steller sea lions was similar, Steller sea lions regularly exhibited lower minimum ƒH, with minimum ƒH almost always less than 20 beats min−1 in dives less than 2 min in duration. In the wild, California sea lions rarely exhibited a minimum ƒH less than 20 beats min−1 in similar duration dives (Fig. 5B), suggesting greater blood oxygen transport during these natural short-duration dives.

Fig. 5. (not shown)  fH decreases with increasing dive duration. Dive duration versus (A) dive fH (total number of beats/dive duration), (B) minimum instantaneous fH and (C) bottom fH (total beats at bottom of dive/bottom time) for California sea lions (461 dives from five sea lions).

Although California sea lions are not usually considered exceptional divers, they exhibited extreme bradycardia, comparable to that of the best diving phocids, during their deep dives. In dives greater than 6 min in duration, minimum ƒH was usually less than 10 beats min−1 and sometimes as low as 6 beats mins−1 (Fig. 5B), which is similar to extreme divers such as emperor penguins (3 beats min−1), elephant seals (3 beats min−1), grey seals (2 beats min−1) and Weddell seals (<10 beats min−1), and even as low as what was observed in forced submersion studies. Thus, similar to phocids, the extreme bradycardia exhibited during forced submersions is also a routine component of the sea lion’s physiological repertoire, allowing them to perform long-duration dives.

While the degree of bradycardia observed in long dives of California sea lions was similar to the extreme bradycardia observed in phocids, the ƒH profiles were quite different. In general, phocid ƒH decreases abruptly upon submergence. The intensity of the initial phocid bradycardia either remains relatively stable or intensifies as the dive progresses, and does not start to increase until the seal begins its ascent. In contrast, the ƒH profiles of sea lions were more complex, showing a more gradual decrease during descent, with the minimum ƒH of the dive usually towards the end of descent (Figs 3, 6). There was often a slight increase in ƒH during the bottom portion of the dive, and as soon as the sea lions started to ascend, the ƒH slowly started to increase, often becoming irregular during the middle of ascent, before increasing rapidly as the sea lion approached the surface.

Fig. 6. (not shown) Instantaneous fH and dive depth profiles of the longest dive (10.0 min, 385 m) from a California sea lion (CSL12_1). During this dive, instantaneous fH reached 7 beats min−1 and was less than 20 beats min−1 for over 5.5 min. Post-dive fH was high in the first 0.5–1 min after surfacing, but then declined to ~100 beats min−1 towards the end of the surface interval.

Implications for pulmonary gas exchange

The moderate dive ƒH in short, shallow dives compared with the much slower ƒH of deep long-duration dives suggests more pulmonary blood flow and greater potential for reliance on lung O2. Most of these dives were to depths of less than 100 m (well below the estimated depth of lung collapse near 200 m), so maintenance of a moderate ƒH during these dives may allow sea lions to maximise use of the potentially significant lung O2 stores (~16% of total body O2 stores) throughout the dive. This is supported by venous blood O2 profiles, where, occasionally, there was no decrease in venous blood O2 between the beginning and end of the dive; this can only occur if pulmonary gas exchange continues throughout the dive. Greater utilization of the lung O2 store in sea lions is consistent with higher dive ƒH in other species that both dive on inspiration and typically perform shallow dives (dolphins, porpoises, some penguin species), and in deeper diving species when they perform shallow dives (emperor penguins).

In deeper dives of sea lions, although ƒH was lower and bradycardia more extreme, the diving ƒH profiles suggest that pulmonary gas exchange is also important. In long-duration dives, even though ƒH started to decrease upon or shortly after submergence, the decrease was not as abrupt as in phocids. Additionally, in long deep dives, despite having overall low dive ƒH, there were more heart beats before resting ƒH was reached compared with short, shallow dives. In dives less than 3 min in duration, there were ~10–15 beats until instantaneous ƒH reached resting values. In longer duration dives (>3 min), there were usually ~30–40 beats before instantaneous ƒH reached resting values. We suggest the greater number of heart beats early in these deeper dives enables more gas exchange and blood O2 uptake at shallow depths, thus allowing utilisation of the postulated larger respiratory O2 stores in deeper dives The less abrupt decline in ƒH we observed in sea lions is similar to the more gradual declines documented in emperor penguins and porpoises, where it has also been proposed that the gradual decrease in ƒH allows them to maximise pulmonary gas exchange at shallower depths. However, as sea lions swam deeper, ƒH decreased further (Figs 3, 6), and by 200 m depth (the approximate depth of lung collapse, instantaneous ƒH was 14 beats min−1. Such an extreme decline in ƒH in conjunction with increased pulmonary shunting due to lung compression at greater depths will result in minimization of both O2 and N2 uptake by blood, even before the depth of full lung collapse (100% pulmonary shunt) is reached.

Implications for blood flow

ƒH is often used as a proxy to estimate blood flow and perfusion during diving because of the relative ease of its measurement. This is based on the assumption that stroke volume does not change during diving in sea lions, and, hence, changes in ƒH directly reflect changes in cardiac output. As breath-hold divers maintain arterial pressure while diving, changes in cardiac output should be associated with changes in peripheral vascular resistance and changes in blood flow to tissues. In Weddell seals, a decrease in cardiac output of ~85% during forced submersions resulted in an 80–100% decrease in tissue perfusion in all tissues excluding the brain, adrenal glands and lung. Sea lions exhibited extremely low instantaneous ƒH values that often remained low for significant portions of the dive (Figs 4, 6), suggesting severe decreases in tissue perfusion in dives greater than 5 min in duration. In almost all dives greater than 6 min in duration, instantaneous ƒH reached 10 beats min−1, and stayed below 20 beats min−1 for more than a minute. At a ƒH of 20 beats min−1, cardiac output will be ~36% of resting cardiac output and only about 18% of average surface cardiac output. At these levels of cardiac suppression, most of this flow should be directed towards the brain and heart.

Conclusions

We successfully obtained diving ƒH profiles from a deep-diving otariid during natural foraging trips. We found that

(1) ƒH decreases during all dives, but true and more intense bradycardia only occurred in longer duration dives and
(2) in the longest duration dives, ƒH and presumed cardiac output were as low as 20% of resting values.

We conclude that, although initial high ƒH promotes gas exchange early in deep dives, the extremely low ƒH in late descent of deep dives (a) preserves lung O2, (b) conserves blood O2, (c) increases the dependence of muscle on myoglobin-bound O2 and (d) limits N2 absorption at depth. This ƒH profile, especially during the late descent/early bottom phase of deep dives is similar to that of deep-diving emperor penguins, and may be characteristic of deep diving endotherms that dive on inspiration.

Dive duration was the fixed effect in all models, and to account for the lack of independence caused by having many dives from the same individual, individual (sea lion ID) was included as a random effect. Covariance and random effect structures of the full models were evaluated using Akaike’s information criterion (AIC) and examination of residual plots. AICs from all the tested models are presented with the best model in bold.

Additionally, dives were classified as short-duration (less than 3 min, minimum cADL), mid-duration (3–5 min, range of cADLs) or long-duration (>5 min) dives. Differences in pre-dive ƒH, dive ƒH, minimum ƒH, post-dive ƒH, and heart beats to resting between the categories were investigated using mixed effects ANOVA, followed by post hoc Tukey tests. In all models, dive duration category was the fixed effect and individual (sea lion ID) was included as a random effect. Model fit was accessed by examination of the residuals. All means are expressed ±s.d. and results of the Tukey tests were considered significant at P<0.05. Statistical analysis was performed in R.

Investigating Annual Diving Behaviour by Hooded Seals (Cystophora cristata) within the Northwest Atlantic Ocean

Julie M. Andersen, Mette Skern-Mauritzen, Lars Boehme
PLoS ONE 8(11): e80438. http://dx.doi.org:/10.1371/journal.pone.0080438

With the exception of relatively brief periods when they reproduce and molt, hooded seals, Cystophora cristata, spend most of the year in the open ocean where they undergo feeding migrations to either recover or prepare for the next fasting period. Valuable insights into habitat use and diving behavior during these periods have been obtained by attaching Satellite Relay Data Loggers (SRDLs) to 51 Northwest (NW) Atlantic hooded seals (33 females and 18 males) during icebound fasting periods (200422008). Using General Additive Models (GAMs) we describe habitat use in terms of First Passage Time (FPT) and analyze how bathymetry, seasonality and FPT influence the hooded seals’ diving behavior described by maximum dive depth, dive duration and surface duration. Adult NW Atlantic hooded seals exhibit a change in diving activity in areas where they spend .20 h by increasing maximum dive depth, dive duration and surface duration, indicating a restricted search behavior. We found that male and female hooded seals are spatially segregated and that diving behavior varies between sexes in relation to habitat properties and seasonality. Migration periods are described by increased dive duration for both sexes with a peak in May, October and January. Males demonstrated an increase in dive depth and dive duration towards May (post-breeding/pre-molt) and August–October (post-molt/pre-breeding) but did not show any pronounced increase in surface duration. Females dived deepest and had the highest surface duration between December and January (post-molt/pre-breeding). Our results suggest that the smaller females may have a greater need to recover from dives than that of the larger males. Horizontal segregation could have evolved as a result of a resource partitioning strategy to avoid sexual competition or that the energy requirements of males and females are different due to different energy expenditure during fasting periods.

Novel locomotor muscle design in extreme deep-diving whales

P. Velten, R. M. Dillaman, S. T. Kinsey, W. A. McLellan and D. A. Pabst
The Journal of Experimental Biology 216, 1862-1871
http://dx.doi.org:/10.1242/jeb.081323

Most marine mammals are hypothesized to routinely dive within their aerobic dive limit (ADL). Mammals that regularly perform deep, long-duration dives have locomotor muscles with elevated myoglobin concentrations that are composed of predominantly large, slow-twitch (Type I) fibers with low mitochondrial volume densities (Vmt). These features contribute to extending ADL by increasing oxygen stores and decreasing metabolic rate. Recent tagging studies, however, have challenged the view that two groups of extreme deep-diving cetaceans dive within their ADLs. Beaked whales (including Ziphius cavirostris and Mesoplodon densirostris) routinely perform the deepest and longest average dives of any air-breathing vertebrate, and short-finned pilot whales (Globicephala macrorhynchus) perform high-speed sprints at depth. We investigated the locomotor muscle morphology and estimated total body oxygen stores of several species within these two groups of cetaceans to determine whether they

(1) shared muscle design features with other deep divers and
(2) performed dives within their calculated ADLs.

Muscle of both cetaceans displayed high myoglobin concentrations and large fibers, as predicted, but novel fiber profiles for diving mammals. Beaked whales possessed a sprinterʼs fiber-type profile, composed of ~80% fast-twitch (Type II) fibers with low Vmt. Approximately one-third of the muscle fibers of short-finned pilot whales were slow-twitch, oxidative, glycolytic fibers, a rare fiber type for any mammal. The muscle morphology of beaked whales likely decreases the energetic cost of diving, while that of short-finned pilot whales supports high activity events. Calculated ADLs indicate that, at low metabolic rates, both beaked and short-finned pilot whales carry sufficient onboard oxygen to aerobically support their dives.

Serial cross-sections of the m. longissimus dorsi of Mesoplodon densirostris

Serial cross-sections of the m. longissimus dorsi of Mesoplodon densirostris

Fig. Serial cross-sections of the m. longissimus dorsi of Mesoplodon densirostris (A–D) and Globicephala macrorhynchus (E–H). Scale bars, 50μm. Muscle sections stained for the alkaline (A,E) and acidic (B,F) preincubations of myosin ATPase were used to distinguish Type I and II fibers. Muscle sections stained for succinate dehydrogenase (C,G) and α-glycerophosphate dehydrogenase (D,H) were used to distinguish glycolytic (gl), oxidative (o) and intermediate (i) fibers.

Previous studies of the locomotor muscles of deep-diving marine mammals have demonstrated that these species share a suite of adaptations that increase onboard oxygen stores while slowing the rate at which these stores are utilized, thus extending ADL. Their locomotor muscles display elevated myoglobin concentrations and are composed predominantly of large Type I fibers. Vmt are also lower in deep divers than in shallow divers or athletic terrestrial species. The results of this study indicate that beaked whales and short-finned pilot whales do not uniformly display these characteristics and that each possesses a novel fiber profile compared with those of other deep divers.

The phylogeny of Cetartiodactyla: The importance of dense taxon sampling, missing data, and the remarkable promise of cytochrome b to provide reliable species-level phylogenies

Ingi Agnarsson, Laura J. May-Collado
Molecular Phylogenetics and Evolution 48 (2008) 964–985
http://dx.doi.org:/10.1016/j.ympev.2008.05.046

We perform Bayesian phylogenetic analyses on cytochrome b sequences from 264 of the 290 extant cetartiodactyl mammals (whales plus even-toed ungulates) and two recently extinct species, the ‘Mouse Goat’ and the ‘Irish Elk’. Previous primary analyses have included only a small portion of the species diversity within Cetartiodactyla, while a complete supertree analysis lacks resolution and branch lengths limiting its utility for comparative studies. The benefits of using a single-gene approach include rapid phylogenetic estimates for a large number of species. However, single-gene phylogenies often differ dramatically from studies involving multiple datasets suggesting that they often are unreliable. However, based on recovery of benchmark clades—clades supported in prior studies based on multiple independent datasets—and recovery of undisputed traditional taxonomic groups, Cytb performs extraordinarily well in resolving cetartiodactyl phylogeny when taxon sampling is dense. Missing data, however, (taxa with partial sequences) can compromise phylogenetic accuracy, suggesting a tradeoff between the benefits of adding taxa and introducing question marks. In the full data, a few species with a short sequences appear misplaced, however, sequence length alone seems a poor predictor of this phenomenon as other taxa.

The mammalian superorder Cetartiodactyla (whales and eventoed ungulates) contains nearly 300 species including many of immense commercial importance (cow, pig, and sheep) and of conservation interest and aesthetic value (antelopes, deer, giraffe, dolphins, and whales) (MacDonald, 2006). Certain members of this superorder count among the best studied organisms on earth, whether speaking morphologically, behaviorally, physiologically or genetically. Understanding the interrelationships among cetartiodactyl species, therefore, is of obvious importance with equally short sequences were not conspicuously misplaced. Although we recommend awaiting a better supported phylogeny based on more character data to reconsider classification and taxonomy within Cetartiodactyla, the new phylogenetic hypotheses provided here represent the currently best available tool for comparative species-level studies within this group. Cytb has been sequenced for a large percentage of mammals and appears to be a reliable phylogenetic marker as long as taxon sampling is dense. Therefore, an opportunity exists now to reconstruct detailed phylogenies of most of the major mammalian clades to rapidly provide much needed tools for species-level comparative studies.

Our results support the following relationship among the four major cetartiodactylan lineages (((Tylopoda ((Cetancodonta (Ruminantia + Suina))), with variable support. This arrangement has not been suggested previously, to our knowledge (see review in O’Leary and Gatesy, 2008 and discussion).

Relationships among clades within Cetancodonta are identical to those found by May-Collado and Agnarsson (2006).

Within Ruminantia all our analyzes suggest the following relationships among families: (((((Tragulidae((((Antilocapridae(((Giraffidae(( Cervidae(Moschidae + Bovidae))))) with relatively high support, supporting the subdivision of Ruminantia into Tragulina and Pecora.
In the rare cases where our results are inconsistent with benchmark clades, ad hoc explanations seem reasonable. The placement of M. meminna (Tragulidae) within Bovidae is likely an artifact of missing data, although remarkably it is the only conspicuous misplacement of a species across the whole phylogeny at the family level (while three species appear to be misplaced at the subfamily level within Cervidae in the full analysis, see Fig. 5a). This is supported by the fact that the placement of Moschiola receives low support, and the removal of Moschiola prior to analysis increases dramatically the support for clades close to where it nested (not shown, analysis available from authors), suggesting it had a tendency to ‘jump around’. Two other possibilities cannot be ruled out, however. One, that possibly the available sequence in Genbank may be mislabeled. And second, it should be kept in mind that the validity of Tragulidae has never been tested with molecular data including more than two species.

Oxygen and carbon dioxide fluctuations in burrows of subterranean blind mole rats indicate tolerance to hypoxic–hypercapnic stresses

Imad Shams, Aaron Avivi, Eviatar Nevo
Comparative Biochemistry and Physiology, Part A 142 (2005) 376 – 382
http://dx.doi.org:/10.1016/j.cbpa.2005.09.003

The composition of oxygen (O2), carbon dioxide (CO2), and soil humidity in the underground burrows from three species of the Israeli subterranean mole rat Spalax ehrenbergi superspecies were studied in their natural habitat. Two geographically close populations of each species from contrasting soil types were probed. Maximal CO2 levels (6.1%) and minimal O2 levels (7.2%) were recorded in northern Israel in the breeding mounds of S. carmeli in a flooded, poor drained field of heavy clay soil with very high volumetric water content. The patterns of gas fluctuations during the measurement period among the different Spalax species studied were similar. The more significant differentiation in gas levels was not among species, but between neighboring populations inhabiting heavy soils or light soils: O2 was lower and CO2 was higher in the heavy soils (clay and basaltic) compared to the relatively light soils (terra rossa and rendzina). The extreme values of gas concentration, which occurred during the rainy season, seemed to fluctuate with partial flooding of the tunnels, animal digging activity, and over-crowded breeding mounds inhabited by a nursing female and her offspring. The gas composition and soil water content in neighboring sites with different soil types indicated large differences in the levels of hypoxic–hypercapnic stress in different populations of the same species. A growing number of genes associated with hypoxic stress have been shown to exhibit structural and functional differences between the subterranean Spalax and the aboveground rat (Rattus norvegicus), probably reflecting the molecular adaptations that Spalax went through during 40 million years of evolution to survive efficiently in the severe fluctuations in gas composition in the underground habitat.

map of the studied sites

map of the studied sites

Schematic map of the studied sites: S. galili (2n =52): 1— Rehania (chalk); 2— Dalton (basaltic); S. golani (2n =54): 3— Majdal Shams (terra tossa); 4—Masa’ada (basaltic soils); S. carmeli (2n =58): 5— Al-Maker (heavy clay); 6— Muhraqa (terra rossa).

Comparison of gas composition (O2 and CO2) and water content between light and heavy soils inhabited by S. carmeli

Comparison of gas composition (O2 and CO2) and water content between light and heavy soils inhabited by S. carmeli

Comparison of gas composition (O2 and CO2) and water content between light and heavy soils inhabited by S. carmeli, Al-Maker (heavy soil) and Muhraqa (light soil). AverageTSD of measurements in the burrows of approximately 10 animals at a given date is presented. **p <0.01, T-test and Mann– Whitney test).

Subterranean mammals, which live in closed underground burrow systems, experience an atmosphere that is different from the atmosphere above-ground. Gas exchange between these two atmospheres depends on diffusion through the soil, which in turn, depends on soil particle size, water content, and burrow depth. Heavy soils (clay and basaltic), hold water and have little air space for gas diffusion. A large deviation from external gas composition is found in the burrows of Spalax living in these soil types. The maximal measured concentration of CO2 was 6.1% in Spalax breeding mounds, which is one of the highest concentrations among studied mammals in natural conditions. At the same time 7.2% O2 was measured in water saturated heavy clay soil

seasonal variation from August to March in mean O2, CO2, and soil water content

seasonal variation from August to March in mean O2, CO2, and soil water content

Example of seasonal variation from August to March in mean O2, CO2, and soil water content (VWC) in the Al-Maker population (2n =58, heavy soil). Values are presented as mean TSD.

In this study new data were presented for a wild mammal that survives in an extreme hypoxic–hypercapnic environment. Interestingly, the very low concentrations of O2 experienced by Spalax are correlated with the expression pattern of hypoxia related genes.  So far, we have shown higher and longer-term mRNA expression of erythropoietin, the main factor that regulates the level of circulating red blood cells, in subterranean Spalax compared to the above-ground rat in response to hypoxic stress, as well as differences in the response of erythropoietin to hypoxia in different populations of Spalax experiencing different hypoxic stress in nature. We also demonstrated that erythropoietin pattern of expression is different in Spalax than in Rattus throughout development, a pattern suggesting more efficient hypoxic tolerance in Spalax starting as early as in the embryonic stages. Furthermore, vascular endothelial growth factor (VEGF), which is a critical angiogenic factor that responds to hypoxia, is constitutively expressed at maximal levels in Spalax muscles, the most energy consuming tissue during digging. This level is 1.6-fold higher than in Rattus muscles and is correlated with significantly higher blood vessel concentration in the Spalax muscles compared to the Rattus muscles. Likewise, myoglobin the globin involved in oxygen homeostasis in skeletal muscles, exhibits different expression pattern under normoxia and in response to hypoxia in Spalax muscles compared to rat muscles as well as between different populations of Spalax exposed to different hypoxic stress in nature (unpublished results). Similarly, neuroglobin, a brain-specific globin involved in reversible oxygen binding, i.e., presumably in cellular homeostasis, is expressed differently in the Spalax brain compared to Rattus brain. Like erythropoietin and myoglobin also neuroglobin is expressed differently in Spalax populations experiencing different oxygen supply (unpublished results). Furthermore, Spalax p53 harbors two amino acid substitutions in its binding domain, which are identical to mutations found in p53 of human cancer cells. These substitutions endow Spalax p53 with several-fold higher activation of cell arrest and DNA repair genes compared to human p53 and favor activation of DNA repair genes over apoptotic genes. The study of specific tumoral variants indicates that such preference of growth arrest over apoptosis possibly results as a response to the hypoxic environmental stress known in tumors. Differences in the structure of other molecules related to homeostasis, namely, hemoglobin, haptoglobin (Nevo, 1999), and cytoglobin (unpublished) were also observed in Spalax.

Stress, adaptation, and speciation in the evolution of the blind mole rat, Spalax, in Israel

Eviatar Nevo
Molecular Phylogenetics and Evolution 66 (2013) 515–525
http://dx.doi.org/10.1016/j.ympev.2012.09.008

Environmental stress played a major role in the evolution of the blind mole rat superspecies Spalax ehrenbergi, affecting its adaptive evolution and ecological speciation underground. Spalax is safeguarded all of its life underground from aboveground climatic fluctuations and predators. However, it encounters multiple stresses in its underground burrows including darkness, energetics, hypoxia, hypercapnia, food scarcity, and pathogenicity. Consequently, it evolved adaptive genomic, proteomic, and phenomic complexes to cope with those stresses. Here I describe some of these adaptive complexes, and their theoretical and applied perspectives. Spalax mosaic molecular and organismal evolution involves reductions or regressions coupled with expansions or progressions caused by evolutionary tinkering and natural genetic engineering. Speciation of Spalax in Israel occurred in the Pleistocene, during the last 2.00–2.35 Mya, generating four species associated intimately with four climatic regimes with increasing aridity stress southwards and eastwards representing an ecological speciational adaptive trend: (Spalax golani, 2n = 54?S. galili, 2n = 52?S. carmeli, 2n = 58?S. judaei, 2n = 60). Darwinian ecological speciation occurred gradually with relatively little genetic change by Robertsonian chromosomal and genic mutations. Spalax genome sequencing has just been completed. It involves multiple adaptive complexes to life underground and is an evolutionary model to a few hundred underground mammals. It involves great promise in the future for medicine, space flight, and deep-sea diving.

Stress is a major driving force of evolution (Parsons, 2005; Nevo, 2011). Parsons defined stress as the ‘‘environmental factor causing potential injurious changes to biological systems with a potential for impacts on evolutionary processes’’. The global climatic transition from the middle Eocene to the early Oligocene (45–35 Ma = Million years ago) led to extensive convergent evolution underground of small subterranean mammals across the planet (Nevo, 1999; Lacey et al., 2000; Bennett and Faulkes, 2000; Begall et al., 2007). The subterranean ecotope provided small mammals with shelter from predators and extreme aboveground climatic stressful fluctuations of temperature and humidity. However, they had to evolve genomic adaptive complexes for the immense underground stresses of darkness, energy for burrowing in solid soil, low productivity and food scarcity, hypoxia, hypercapnia, and high infectivity. These stresses have been described in Nevo (1999, 2011) and Nevo et al. (2001); and Nevo list of Spalax publication at http://evolution.haifa.ac.il with many cited references relevant to these stresses).

blind subterranean mole rat of the Spalax ehrenbergi superspecies

blind subterranean mole rat of the Spalax ehrenbergi superspecies

The blind subterranean mole rat of the Spalax ehrenbergi superspecies in Israel. An extreme example of adaptation to life underground

Circadian rhythm and genes

adaptive circadian genes. We identified the circadian rhythm of Spalax
(Nevo et al., 1982) and described, cloned, sequenced, and expressed several circadian genes in Spalax. These include Clock, MOP3, three Period (Per), and cryptochromes (Avivi et al., 2001, 2002, 2003). The Spalax circadian genes are differentially conserved, yet characterized by a significant number of amino acid substitutions. The glutamine-rich area of Clock, which is assumed to function in circadian rhythmicity, is expanded in Spalax compared with that of mice and humans and is different in amino acid composition from that of rats. All three Per genes of Spalax oscillate with a periodicity of 24 h in the suprachaismatic nucleus, eye, and Harderian gland and are expressed in peripheral organs. Per genes are involved in clock resetting. Spalax Per 3 is unique in mammals though its function is still unresolved. The Spalax Per genes contribute to the unique adaptive circadian rhythm to life underground. The cryptochrome (Cry) genes, found in animals and plants, act both as photoreceptors and as ingredients of the negative feedback mechanism of the biological Clock. The CRY 1 protein is significantly closer to the human homolog than to that of mice, as was also shown in parts of the immunogenetic system. Both Cry 1 and Cry 2 mRNAs were found in the SCN, eye, harderian gland, and in peripheral tissues. Remarkably, the distinctly hypertrophied harderian gland is central in Spalax’s unique underground circadian rhythmicity (Pevet et al., 1984).

  • Spalax eye mosaic evolution
  • Gene expression in the eye of Spalax
  • Brain evolution in Spalax to underground stresses
  • Spalax: four species in Israel

The morphological, physiological, and behavioral Spalax eye patterns are underlain by gene expression representing regressive and progressive associated transcripts. Regressive transcripts involve B-2 microglobulin, transketolase, four keratins, alpha enolase, and different heat shock proteins. Several proteins may be involved in eye degeneration. These include heat shock protein 90alpha (hsp90alpha), found also in the blind fish Astyanax mexicanus, two transcripts of programmed cell death proteins, oculospanin, and peripherin 2, both belonging to the Tetraspanin family, in which 60 different mutations cause eye degeneration in humans. Several progressive transcripts in the Spalax eye are found in the retina of many mammals involving gluthatione, peroxidase 4, B spectrin, and Ankyrin; the last two characterize rod cells in the retina. Some transcripts are involved in metabolic processing of retinal, a vertebrate key component in phototransduction, and a relative of vitamin A.

cross section of the developing eye of the mole rat

cross section of the developing eye of the mole rat

Light micrographs showing cross section of the developing eye of the mole rat Spalax ehrenbergi. (A) Optic cup and lens vesicle initially develop normally (x100). (B) Eye at a later embryonic stage. Note appearance of iris-ciliary body rudiment (arrows), and development of the lens nucleus (L). ON, optic nerve (x100). (C) Eye at a still later fetal stage. Note massive growth of the iris-ciliary body complex colobomatous opening (arrow) (x100). (D) Early postnatal stage. The iris-ciliary body complex completely fills the chamber. The lens is vascularized and vacuolated (x100). (E) Adult eye. Eyelids are completely closed and pupil is absent. Note atrophic appearance of the optic disc region (arrow) (x65). (F) Higher magnification of the adult retina. The different retinal layers are retained: PE, pigment epithelium: RE, receptor layer; ON, outer nuclear layer: IN, inner nuclear layer; GC, ganglion cell layer (x500) (from Sanyal et al., 1990, Fig. 1).

The brains of subterranean mammals underwent dramatic evolution in accordance with underground stresses for digging and photoperiodic perception associated with vibrational, tactile, vocal, olfactory, and magnetic communication systems replacing sight, as is seen in Spalax. The brain of Spalax is twice as large as that of the laboratory rat of the same body size. The somatosensory region in the isocortex of Spalax is 1.7 times, the thalamic nuclei 1.3 times, and the motor cortex 3.1 times larger than in the sighted laboratory rat Rattus norvegicus matched to body size.

The ecological stress determinant in Spalax brain evolution is highlighted by the four species of the Spalax ehrenbergi superspecies in Israel. They differentiated chromosomally (by means of Robertsonian mutations and fission), allopatrically, and clinally southwards into four species associated with different climatic regimes, following the gradient of increasing aridity stress and decreasing predictability southwards towards the desert: Spalax galili (2n = 52) ->S. golani (2n = 54)->S. carmeli (2n = 58)->S. judaei (2n = 60), and eastwards S. galili ->S. golani (2n = 52–>54) (Fig. 2). This chromosomal speciation trend southwards is associated with the regional aridity stress southwards (and eastwards) in Israel, budding new species adapted genomically, proteomically, and phenomically (i.e., in morphology, physiology, and behavior) to increasing stresses of higher solar radiation, temperature, and drought southwards (Nevo, 1999; Nevo et al., 2001; Nevo
list of Spalax at http://evolution.haifa.ac.il). A uniquely recent discovery of incipient sympatric ecological speciation at a microscale in Spalax triggered by local stresses occurs within Spalax galili.

retinal input to primary visual structures in Spalax

retinal input to primary visual structures in Spalax

Relative degree of retinal input to primary visual structures in Spalax, hamster, rat, and Spalacopus cyanus (South American Octodontidae, ‘‘coruro’’). These rodents are of similar body size (120–140 g). B. Relative degree of change in the proportions of retinal input to different primary visual structures in Spalax compared with measures obtained in other rodents. A relative progressive development in Spalax is seen in structures involved in photoperiodic and neuroendocrine functions (SCN, BNST).The main regressive feature is the drastic relative reduction of retinal input to the superior colliculus. The main regressive feature is the drastic reduction of retinal input to the superior colliculus. The relative size of other visual structures in Spalax is modified compared to that of the other species. c. Comparison of the absolute size (volume, mm3 x 10-4) of visual structures in Spalax and other rodents. The size of the SCN is equivalent in all species. The vLGN and dLGN are reduced by 87–93% in Spalax. The retino-recipient layers of the superior colliculus are reduced by 97%. Abbreviations: SCN: suprachiasmatic nucleus; BNST: bed nucleus of the stria terminalis; dLGN: dorsal lateral geniculate nucleus; SC: superior colliculus [From Cooper et al., 1993 (Fig 3)].

Subterranean life has a high energetic cost if an animal has to burrow in order to obtain its food. For a 150 g Thomomys bottae, burrowing 1 m may be 360–3400 times more expensive energetically than moving the same distance on the surface (Vleck, 1979). Mean rates of oxygen consumption during burrowing at 22 oC are from 2.8 to 7.1 times the RMR. Vleck developed a model examining the energetics of foraging by burrowing and found that, in the desert, Thomomys adjusts the burrow segment length to minimize the cost of burrowing. Since burrowing becomes less economic as body size increases, Vleck (1981) predicted that the maximum possible body size that a subterranean mammal can attain depends on a balance between habitat productivity and the cost of burrowing in local soils. Vleck’s cost of burrowing hypothesis has been verified in multiple cases. Heth (1989) demonstrated longer burrows in the rendzina soil and shorter ones in the terra rossa soil, associating lower productivity in the former for Spalax.

Food is a limiting factor for subterranean mammals. The abundance and distribution of food explain some of the ecological, physiological, and behavioral characteristics of subterranean mammals. In a field test of Spalax foraging strategy, we concluded that Spalax was a generalist due to the constraints of the subterranean ecotope. Restricted foraging time primarily during the winter when soil is wet, and the high energetic investment of tunneling to get to food items is significantly reduced than in summertime.
We also identified a decrease in the basic metabolic rate towards the desert, i.e., economizing energetics. The maintenance of adequate O2 transport in a subterranean mammal confronting hypoxia requires adaptation along the O2 transport system, achieved by increasing the flow of O2 in the convection systems (ventilation and perfusion) and by reduction of oxygen pressure (PO2) gradients at the diffusion barriers (lung blood, blood-tissue (Arieli, 1990). The PO2 gradient between blood capillaries and respiring mitochondria capillaries is large, and any adaptation at this level could be significant for O2 transport. Reduction of diffusion distance in a muscle can be achieved, like in Spalax, by increasing the number of capillaries that surround muscle fiber or by reducing fiber areas.

Geographic distribution in Israel of the four chromosomal species belonging to the S. ehrenbergi superspecies

Geographic distribution in Israel of the four chromosomal species belonging to the S. ehrenbergi superspecies

Geographic distribution in Israel of the four chromosomal species belonging to the S. ehrenbergi superspecies that are separated by narrow hybrid zones (2n = 52, 54, 58, and 60, now named as S. galili, S. golani, S. carmeli, and S. judaei, respectively; see Nevo et al., 2001).

Spalacid evolution, based on mtDNA, is driven by climatic oscillations and stresses. The underground ecotope provided subterranean mammals with shelter from extreme climate (temperature and humidity) fluctuations, and predators. However, they had to extensively and intensively adapt to the multiple underground stresses (darkness, energetic, low productivity and
food scarcity, hypoxia, hypercapnia, and high infectivity). All subterranean mammals, including spalacids as an extreme case, share convergent molecular and organismal adaptations to their shared unique underground ecotope. Evolution underground, as exemplified here in spalacids, led to mosaic molecular and organismal evolutionary syndromes to cope with multiple stresses.

Speciation involves all rates – from gradual to rapid. Subterranean mammals, with the spalacid example discussed above, provide uniquely rich evolutionary global tests of speciation and adaptation, convergence, regression, progression, and mosaic evolutionary processes. Adaptation and speciation underground was one of the most dramatic natural experiments verifying Darwinian evolution.

The Spalax genome sequencing has just been completed. It is being analyzed and will soon be published in 2012. This will be a milestone in understanding how numerous mammals across the globe, who found underground shelter from climatic fluctuations and stresses above ground, cope with the new suite of stresses they encountered underground, demanding a new engineering overhaul on all organizational levels, selecting for adaptive complexes to cope with the new underground stresses. The main current and future challenges are to compare and contrast genome sequences and identify the genomic basis of adaptation and speciation.

This global Cenozoic experiment could answer the following open questions: How heterozygous is the whole genome? How prevalent are retrotransposons and what is their functional role? How many genes are involved in the Spalax genome and how are they regulated? What are the genic and regulatory networks resisting the multiple stresses underground? How much of the Spalax genome is conserved and how much is reorganized to cope with the underground stresses? How is the solitary blind mole rat, Spalax, different from the social naked mole rat Heterocephalus? How are the processes of reduction, expansion, and genetic tinkering and engineering reflected across the genome? How effective is copy number variation in regulation? Is there similarity in the transcriptomes of subterranean mammals? How could we harness the rich genome repertoire of Spalax to revolutionize medicine, especially in the realm of hypoxia tolerance and the related major diseases of the western world, e.g., cancer, stroke, and cardiovascular diseases? What is the phylogenetic origin of Spalax? How much of the Spalax genome represents its phylogenetic roots and how much of coding and noncoding genomic regions are shared with other subterranean mammals across the globe in adapting to life underground?

The Atmospheric Environment of the Fossorial Mole Rat (Spalax Ehrenbergi): Effects of Season, Soil Texture, Rain, Temperature and Activity

  1. Arieli
    Comp Biochen Physiol. 1978; 63A:569-5151. The fossorial mole rat (Spalax ehrenbergi) may inhabit heavy soil with low gas permeability.
  2. Air composition in burrows in heavy soil deviates from atmospheric air more than that of burrows in light soil.
  3. In winter and spring O2 and CO2 concentrations in breeding mounds were 16.5% O2 and 2.5-3x CO2 and the extreme values measured were 14.0% O2 and 4.8% Cot.
  4. Hypoxia and hypercapnia in the burrow develop shortly after rain and when ambient temperature drops.
  5. Composition of the burrows air is influenced by the solubility of CO2 in soil water and by faster penetration of oxygen than outflowing of CO2.

Hypo-osmotic stress-induced physiological and ion-osmoregulatory responses in European sea bass (Dicentrarchus labrax) are modulated differentially by nutritional status

Amit Kumar Sinha, AF Dasan, R Rasoloniriana, N Pipralia, R Blust, G De Boeck
Comparative Biochemistry and Physiology, Part A 181 (2015) 87–99
http://dx.doi.org/10.1016/j.cbpa.2014.11.024

We investigated the impact of nutritional status on the physiological, metabolic and ion-osmoregulatory performance of European sea bass (Dicentrarchus labrax)when acclimated to seawater (32 ppt), brackishwater (20 and 10 ppt) and hyposaline water (2.5 ppt) for 2 weeks. Following acclimation to different salinities, fish were either fed or fasted (unfed for 14 days). Plasma osmolality, [Na+], [Cl−] and muscle water contentwere severely altered in fasted fish acclimated to 10 and 2.5 ppt in comparison to normal seawater-acclimated fish, suggesting ion regulation and acid–base balance disturbances. In contrast to feed-deprived fish, fed fish were able to avoid osmotic perturbation more effectively. This was accompanied by an increase in Na+/K+-ATPase expression and activity, transitory activation of H+-ATPase (only at 2.5 ppt) and down-regulation of Na+/K+/2Cl− gene expression. Ammonia excretion rate was inhibited to a larger extent in fasted fish acclimated to low salinities while fed fish were able to excrete efficiently. Consequently, the build-up of ammonia in the plasma of fed fish was relatively lower. Energy stores, especially glycogen and lipid, dropped in the fasted fish at low salinities and progression towards the anaerobic metabolic pathway became evident by an increase in plasma lactate level. Overall, the results indicate no osmotic stress in both feeding treatments within the salinity range of 32 to 20 ppt. However, at lower salinities (10–2.5 ppt) feed deprivation tends to reduce physiological, metabolic, ion-osmo-regulatory and molecular compensatory mechanisms and thus limits the fish’s abilities to adapt to a hypo-osmotic environment.

The absence of ion-regulatory suppression in the gills of the aquatic air-breathing fish Trichogaster lalius during oxygen stress

Chun-Yen Huang, Hsueh-Hsi Lin, Cheng-Huang Lin, Hui-Chen Lin
Comparative Biochemistry and Physiology, Part A 179 (2015) 7–16
http://dx.doi.org/10.1016/j.cbpa.2014.08.017

The strategy for most teleost to survive in hypoxic or anoxic conditions is to conserve energy expenditure, which can be achieved by suppressing energy-consuming activities such as ion regulation. However, an air-breathing fish can cope with hypoxic stress using a similar adjustment or by enhancing gas exchange ability, both behaviorally and physiologically. This study examined Trichogaster lalius, an air-breathing fish without apparent gill modification, for their gill ion-regulatory abilities and glycogen utilization under a hypoxic  treatment. We recorded air-breathing frequency, branchial morphology, and the expression of ion-regulatory proteins (Na+/K+-ATPase and vacuolar-type H+-ATPase) in the 1st and 4th gills and labyrinth organ (LO), and the expression of glycogen utilization (GP, glycogen phosphorylase protein expression and glycogen content) and other protein responses (catalase, CAT; carbonic anhydrase II, CAII; heat shock protein 70, HSP70; hypoxia-inducible factor-1α, HIF-1α; proliferating cell nuclear antigen, PCNA; superoxidase dismutase, SOD) in the gills of T. lalius after 3 days in hypoxic and restricted conditions. No morphological modification of the 1st and 4th gills was observed. The air breathing behavior of the fish and CAII protein expression both increased under hypoxia. Ion-regulatory abilities were not suppressed in the hypoxic or restricted groups, but glycogen utilization was enhanced within the groups. The expression of HIF-1α, HSP70 and PCNA did not vary among the treatments. Regarding the antioxidant system, decreased CAT enzyme activity was observed among the groups. In conclusion, during hypoxic stress, T. lalius did not significantly reduce energy consumption but enhanced gas exchange ability and glycogen expenditure.

The combined effect of hypoxia and nutritional status on metabolic and ionoregulatory responses of common carp (Cyprinus carpio)

Sofie Moyson, HJ Liew, M Diricx, AK Sinha, R Blusta, G De Boeck
Comparative Biochemistry and Physiology, Part A 179 (2015) 133–143
http://dx.doi.org/10.1016/j.cbpa.2014.09.017

In the present study, the combined effects of hypoxia and nutritional status were examined in common carp (Cyprinus carpio), a relatively hypoxia tolerant cyprinid. Fish were either fed or fasted and were exposed to hypoxia (1.5–1.8mgO2 L−1) at or slightly above their critical oxygen concentration during 1, 3 or 7 days followed by a 7 day recovery period. Ventilation initially increased during hypoxia, but fasted fish had lower ventilation frequencies than fed fish. In fed fish, ventilation returned to control levels during hypoxia, while in fasted fish recovery only occurred after reoxygenation. Due to this, C. carpio managed, at least in part, to maintain aerobic metabolism during hypoxia: muscle and plasma lactate levels remained relatively stable although they tended to be higher in fed fish (despite higher ventilation rates). However, during recovery, compensatory responses differed greatly between both feeding regimes: plasma lactate in fed fish increased with a simultaneous breakdown of liver glycogen indicating increased energy use, while fasted fish seemed to economize energy and recycle decreasing plasma lactate levels into increasing liver glycogen levels. Protein was used under both feeding regimes during hypoxia and subsequent recovery: protein levels reduced mainly in liver for fed fish and in muscle for fasted fish. Overall, nutritional status had a greater impact on energy reserves than the lack of oxygen with a lower hepatosomatic index and lower glycogen stores in fasted fish. Fasted fish transiently increased Na+/K+-ATPase activity under hypoxia, but in general ionoregulatory balance proved to be only slightly disturbed, showing that sufficient energy was left for ion regulation.

The effect of temperature and body size on metabolic scope of activity in juvenile Atlantic cod Gadus morhua L.

Bjørn Tirsgaard, Jane W. Behrens, John F. Steffensen
Comparative Biochemistry and Physiology, Part A 179 (2015) 89–94
http://dx.doi.org/10.1016/j.cbpa.2014.09.033

Changes in ambient temperature affect the physiology and metabolism and thus the distribution of fish. In this study we used intermittent flow respirometry to determine the effect of temperature (2, 5, 10, 15 and 20 °C) and wet body mass (BM) (~30–460 g) on standard metabolic rate (SMR, mg O2 h−1), maximum metabolic rate (MMR, mg O2 h−1) and metabolic scope (MS, mg O2 h−1) of juvenile Atlantic cod. SMR increased with BM irrespectively of temperature, resulting in an average scaling exponent of 0.87 (0.82–0.92). Q10 values were 1.8–2.1 at temperatures between 5 and 15 °C but higher (2.6–4.3) between 2 and 5 °C and lower (1.6–1.4) between 15 and 20 °C in 200 and 450 g cod. MMR increased with temperature in the smallest cod (50 g) but in the larger cod MMR plateaued between 10, 15 and 20 °C. This resulted in a negative correlation between the optimal temperature for MS (Topt) and BM, Topt being respectively 14.5, 11.8 and 10.9 °C in a 50, 200 and 450 g cod. Irrespective of BM cold water temperatures resulted in a reduction (30–35%) of MS whereas the reduction of MS at warm temperatures was only evident for larger fish (200 and 450 g), caused by plateauing of MMR at 10 °C and above. Warm temperatures thus seem favorable for smaller (50 g) juvenile cod, but not for larger conspecifics (200 and 450 g).

Read Full Post »

Evolution and Medicine

Reporter and Curator: Larry H. Bernstein, MD, FCAP 

 

http://paleoaerie.org/2015/01/21/what-has-evolution-done-for-me-lately/

Excerpt of article

Cancer is an inescapable fact of life. All of us with either die from it or know someone who will. Cancer is so prevalent because it isn’t a disease in the way a flu or a cold is. No outside force or germ is needed to cause cancer (although it can). It arises from the very way we are put together.  Most of the genes that are needed for multicellular life have been found to be associated with cancer. Cancer is a result of our natural genetic machinery that has been built up over billions of years breaking down over time.

CLONAL EVOLUTION OF CANCER. MEL GREAVES.HTTP://WWW.SCIENCE-CONNECTIONS.COM/TRENDS/SCIENCE_CONTENT/EVOLUTION_6.HTM

Cancer is not only a result of evolutionary processes, cancer itself follows evolutionary theory as it grows. The immune system places a selective pressure on cancer cells, keeping it in check until the cancer evolves a way to avoid it and surpass it in a process known as immunoediting. Cancers face selective pressures in the microenvironments in which they grow. Due to the fast growth of cancer cells, they suck up oxygen in the tissues, causing wildly fluctuating oxygen levels as the body tries to get oxygen to the tissues. This sort of situation is bad for normal tissues and so it is for cancer, at least until they evolve and adapt. At some point, some cancer cells will develop the ability to use what is called aerobic glycolysis to make the ATP we use for energy. Ordinarily, our cells only use glycolysis when they run out of oxygen because aerobic respiration (aka oxidative phosphorylation) is far more efficient. Cancer cells, on the other hand, learn to use glycolysis all the time, even in the presence of abundant oxygen. They may not grow as quickly when there is plenty of oxygen, but they are far better than normal cells at hypoxic, or low oxygen, conditions, which they create by virtue of their metabolism. Moreover, they are better at taking up nutrients because many of the metabolic pathways for aerobic respiration also influence nutrient uptake, so shifting those pathways to nutrient uptake rather than metabolism ensures cancer cells get first pick of any nutrients in the area. The Warburg Effect, as this is called, works by selective pressures hindering those cells that can’t do so and favoring those that can. Because cancer cells have loose genetic controls and they are constantly dividing, the cancer population can evolve, whereas the normal cells cannot.

Evolutionary theory can also be used to track cancer as it metastasizes. If a person has several tumors, it is possible to take biopsies of each one and use standard cladistic programs that are normally used to determine evolutionary relationships between organisms to find which tumor is the original tumor. If the original tumor is not one of those biopsied, it will tell you where the cancer originated within the body. You can thus track the progression of cancer throughout a person’s body. Expanding on this, one can even track the effect of cancer through its effects on how organisms interact within ecosystems, creating its own evolutionary stamp on the environment as its effects radiate throughout the ecosystem.

I’ve talked about cancer at decent length (although I could easily go one for many more pages) because it is less well publicly known than some of the other ways that evolutionary theory helps us out in medicine. The increasing resistance of bacteria and viruses to antibiotics is well known. Antibiotic resistance follows standard evolutionary processes, with the result that antibiotic resistant bacteria are expected to kill 10 million people a year by 2050.  People have to get a new flu shot every year because the flu viruses are legion and they evolve rapidly to bypass old vaccinations.  If we are to accurately predict how the viruses may adapt and properly prepare vaccines for the coming year, evolutionary theory must be taken into account. Without it, the vaccines are much less likely to be effective. Evolutionary studies have pointed out important changes in the Ebola virus and how those changes areaffecting its lethality, which will need to be taken into account for effective treatments. Tracking the origins of viruses, like the avian flu or swine flu, gives us information that will be useful in combating them or even stopping them at their source before they become a problem.

HTTP://WWW.MEDSCAPE.COM/VIEWARTICLE/756378

 

 

Read Full Post »

Highlights of a Green Evolution

Reporter and Curator: Larry H Bernstein, MD, FCAP 

 

 

Chlorophyll

chlorophyll coloration to leaves

chlorophyll coloration to leaves

Paul May
School of Chemistry, University of Bristol
VRML, Jmol, and Chime versions

Chlorophyll is the molecule that absorbs sunlight and uses its energy to
synthesize carbohydrates from CO2 and water. This process is known as
photosynthesis. Animals and humans obtain their food supply by eating plants.

In 1780, the famous English chemist Joseph Priestley found that plants could “restore air which has been injured by the burning of candles.” He placed a mint
plant into a vessel of water for several days, then found that “the air would neither extinguish a candle, nor was it all inconvenient to a mouse which I put into it”.
He discovered that plants produce oxygen. Then, in 1794,  Antoine Lavoisier
discovered oxidation.  It fell to a Dutchman, Jan Ingenhousz,  to make the next
major contribution to the mechanism of photosynthesis.
Having heard of Priestley’s experiments, he  spent a summer near London doing
over 500 experiments, to discover that light plays a major role in photosynthesis.
He noted that plants not only have the faculty to correct bad air in six to ten days,
but they perform this in a few hours; owing to the influence of light of the sun
upon the plant.

Very soon after, more pieces of the puzzle were found by two chemists working
in Geneva. Jean Senebier, found that “fixed air” (CO2) was taken up during photosynthesis, and Theodore de Saussure discovered that the other reactant
necessary was water. The final contribution came from a German surgeon,
Julius Robert Mayer ,

Julius Robert Mayer

Julius Robert Mayer

who recognised that plants convert solar energy into chemical energy. He said:
“Nature has put itself the problem of how to catch in flight light streaming to
the Earth and to store the most elusive of all powers in rigid form. The plants
take in one form of power, light; and produce another power, chemical
difference.” The actual chemical equation which takes place is the reaction
between carbon dioxide and water, catalyzed by sunlight, to produce glucose
and a waste product, oxygen. The glucose sugar is either directly used as an
energy source by the plant for metabolism or growth, or is polymerized to form
starch, so it can be stored until needed. The waste oxygen is excreted into the
atmosphere, where it is made use of by plants and animals for respiration.

http://www.chm.bris.ac.uk/motm/chlorophyll/photosth.gif

Chlorophyll as a Photoreceptor

Chlorophyll is the molecule that traps this ‘most elusive of all powers’ – and is
called a photoreceptor. It is found in the chloroplasts of green plants,
and is what makes green plants, green. The basic structure of a chlorophyll
molecule is a porphyrin ring, co-ordinated to a central atom. This is very
similar in structure to the heme group found in hemoglobin, except that in
heme the central atom is iron, whereas in chlorophyll it is magnesium.

chphyll

http://www.chm.bris.ac.uk/motm/chlorophyll/chphyll.gif

Click for 3D structure file

Click for 3D structure file

There are actually 2 main types of chlorophyll, named a and b. They differ only
slightly, in the composition of a sidechain (in a it is – H3, in b it is CHO). Both of these
two chlorophylls are very effective photoreceptors because they contain a network of
alternating single and double bonds, and the orbitals can delocalize stabilizing the
structure. Such delocalised polyenes have very strong absorption bands in the visible
regions of the spectrum, allowing the plant to absorb the energy from sunlight.

chloroabs

http://www.chm.bris.ac.uk/motm/chlorophyll/chloroabs.gif

The different side groups in the 2 chlorophylls ‘tune’ the absorption spectrum to
slightly different wavelengths, so that light that is not significantly absorbed by
chlorophyll a, at, say, 460nm, will instead be captured by chlorophyll b, which
absorbs strongly at that wavelength. Thus these two kinds of chlorophyll
complement each other in absorbing sunlight. Plants can obtain all their energy
requirements from the blue and red parts of the spectrum, however, there is still
a large spectral region, between 500-600nm, where very little light is absorbed.

This light is in the green region of the spectrum, and since it is reflected, this
is the reason plants appear green. Chlorophyll absorbs so strongly that it can
mask other less intense colours. Some of these more delicate colours (from
molecules such as carotene and quercetin) are revealed when the chlorophyll
molecule decays in the Autumn, and the woodlands turn red, orange,and
golden brown. Chlorophyll can also be damaged when vegetation is cooked,
since the central Mg atom is replaced by hydrogen ions. This affects the energy
levels within the molecule, causing its absorbance spectrum to alter. Thus cooked
leaves change colour – often becoming a paler, insipid yellowy green.

As the chlorophyll in leaves decays in the autumn, the green colour fades and is
replaced by the oranges and reds of carotenoids.

Chlorophyll in Plants

The chlorophyll molecule is the active part that absorbs the sunlight, but just as with
hemoglobin, in order to do its job (synthesising carbohydrates) it needs to be attached
to the backbone of a very complicated protein. This protein may look haphazard in
design, but it has exactly the correct structure to orient the chlorophyll molecules in
the optimal position to enable them to react with nearby CO2 and H2O molecules in
a very efficient manner. Several chlorophyll molecules are lurking inside this bacterial
photoreceptor protein (right).

References:

Introduction to Organic Chemistry, Streitweiser and Heathcock (MacMillan, New York,
1981).

Biochemistry, L. Stryer (W.H. Freeman and Co, San Francisco, 1975).

Wikipedia – Chlorophyll

Chlorophyll (also chlorophyl) is a green pigment found in cyanobacteria and the
chloroplasts of algae and plants.  Its name is derived from the Greek words χλωρός,
chloros (“green”) and φύλλον, phyllon (“leaf”).  Chlorophyll is an extremely important
biomolecule, critical in photosynthesis, which allows plants to absorb energy from light. Chlorophyll absorbs light most strongly in the blue portion of the
electromagnetic spectrum, followed by the red portion. Conversely, it is a poor
absorber of green and near-green portions of the spectrum, hence the green
color of chlorophyll-containing tissues. chlorophyll was first isolated by
Joseph Bienaimé Caventou and Pierre Joseph Pelletier in 1817.

Absorption maxima of chlorophylls against the spectrum of white light

Chlorofilab.svg

Chlorophyll is found in high concentrations in chloroplasts of plant cells.

Clorofila_3
http://upload.wikimedia.org/wikipedia/commons/thumb/0/05/Clorofila_3.jpg/
120px-Clorofila_3.jpg

These chlorophyll maps show milligrams of chlorophyll per cubic meter of seawater
each month. Places where chlorophyll amounts were very low, indicating very low
numbers of phytoplankton, are blue. Places where chlorophyll concentrations were
high, meaning many phytoplankton were growing, are yellow.

chlophyll world map

chlophyll world map

http://upload.wikimedia.org/wikipedia/commons/thumb/e/e3/
MY1DMM_CHLORA.ogv/220px–MY1DMM_CHLORA.ogv.jpg

Chlorophyll and photosynthesis

Chlorophyll is vital for photosynthesis, which allows plants to absorb energy from light.

Chlorophyll molecules are specifically arranged in and around photosystems that are
embedded in the thylakoid membranes of chloroplasts. In these complexes,
chlorophyll serves two primary functions. The function of the vast majority of
chlorophyll (up to several hundred molecules per photosystem) is to absorb light and
transfer that light energy by resonance energy transfer to a specific chlorophyll pair
in the reaction center of the photosystems.

The two currently accepted photosystem units are Photosystem II and Photosystem I,
which have their own distinct reaction center chlorophylls, named P680 and P700,
respectively. These pigments are named after the wavelength (in nanometers) of their
red-peak absorption maximum. The identity, function and spectral properties of the
types of chlorophyll in each photosystem are distinct and determined by each other
and the protein structure surrounding them. Once extracted from the protein into a
solvent (such as acetone or methanol), these chlorophyll pigments can be separated
in a simple paper chromatography experiment and, based on the number of polar
groups between chlorophyll a and chlorophyll b, will chemically separate out on the
paper.

The function of the reaction center chlorophyll is to use the energy absorbed by and
transferred to it from the other chlorophyll pigments in the photosystems to undergo
a charge separation, a specific redox reaction in which the chlorophyll donates an
electron into a series of molecular intermediates called an electron transport chain.
The charged reaction center chlorophyll (P680+) is then reduced back to its ground
state by accepting an electron. In Photosystem II, the electron that reduces P680+
ultimately comes from the oxidation of water into O2 and H+ through several
intermediates.

This reaction is how photosynthetic organisms such as plants produce O2 gas, and
is the source for practically all the O2 in Earth’s atmosphere. Photosystem I typically
works in series with Photosystem II; thus the P700+ of Photosystem I is usually
reduced, via many intermediates in the thylakoid membrane, by electrons ultimately
from Photosystem II. Electron transfer reactions in the thylakoid membranes are
complex, however, and the source of electrons used to reduce P700+ can vary.

The electron flow produced by the reaction center chlorophyll pigments is used to
shuttle H+ ions across the thylakoid membrane, setting up a chemiosmotic potential
used mainly to produce ATP chemical energy; and those electrons ultimately reduce
NADP+ to NADPH, a universal reductant used to reduce CO2 into sugars as well as
for other biosynthetic reductions.

Reaction center chlorophyll–protein complexes are capable of directly absorbing light
and performing charge separation events without other chlorophyll pigments, but the
absorption cross section (the likelihood of absorbing a photon under a given light
intensity) is small. Thus, the remaining chlorophylls in the photosystem and antenna
pigment protein complexes associated with the photosystems all cooperatively absorb
and funnel light energy to the reaction center. Besides chlorophyll a, there are other
pigments, called accessory pigments, which occur in these pigment–protein
antenna complexes.

Chemical structure

Chlorophyll is a chlorin pigment, which is structurally similar to and produced through the same metabolic pathway as other porphyrin pigments such as heme. At the center
of the chlorin ring is a magnesium ion. This was discovered in 1906, and was the first
time that magnesium had been detected in living tissue. or the structures depicted in

this article, some of the ligands attached to the Mg2+ center are omitted for clarity.
The chlorin ring can have several different side chains, usually including a long
phytol chain. There are a few different forms that occur naturally, but the most
widely distributed form in terrestrial plants is chlorophyll a.

Chlorophyll-a-3D

Chlorophyll-a-3D

http://upload.wikimedia.org/wikipedia/commons/thumb/9/92/
Chlorophyll-a-3D-vdW.png/220px-Chlorophyll-a-3D-vdW.png

Space-filling model of the chlorophyll a molecule

After initial work done by German chemist Richard Willstätter spanning from 1905 to
1915, the general structure of chlorophyll a was elucidated by Hans Fischer in 1940.
By 1960, when most of the stereochemistry of chlorophyll a was known, Robert Burns
Woodward published a total synthesis of the molecule. In 1967, the last remaining
stereochemical elucidation was completed by Ian Fleming, and in 1990 Woodward
and co-authors published an updated synthesis. Chlorophyll was announced to be
present in cyanobacteria and other oxygenic microorganisms that form stromatolites
in 2010; a molecular formula of C55H70O6N4Mg and a structure of (2-formyl)-chlorophyll a were deduced based on NMR, optical and mass spectra.

When leaves degreen in the process of plant senescence, chlorophyll is converted
to a group of colourless tetrapyrroles known as nonfluorescent chlorophyll catabolites
(NCC’s) with the general structure:

These compounds have also been identified in several ripening fruits

Nonfluorescentchlorophilcatabolites.svg

http://upload.wikimedia.org/wikipedia/commons/thumb/c/c7/Nonfluorescent
chlorophilcatabolites.svg/241px-Nonfluorescentchlorophilcatabolites.svg.png

Absorbance spectra of free chlorophyll a (blue) and b (red) in a solvent. The spectra
of chlorophyll molecules are slightly modified in vivo depending on specific pigment-
protein interactions.

Chlorophyll_ab_spectra

Chlorophyll_ab_spectra

http://upload.wikimedia.org/wikipedia/commons/thumb/2/23/Chlorophyll_ab_
spectra-en.svg/220px-Chlorophyll_ab_spectra-en.svg.png

Complementary light absorbance of anthocyanins with chlorophylls

Anthocyanins are other plant pigments. The absorbance pattern responsible for the
red color of anthocyanins may be complementary to that of green chlorophyll in
photosynthetically active tissues such as young Quercus coccifera leaves. It may
protect the leaves from attacks by plant eaters that may be attracted by green color.

Superposition of spectra of chlorophyll a and b with oenin (malvidin 3O glucoside),
a typical anthocyanidin, showing that, while chlorophylls absorb in the blue and
yellow/red parts of the visible spectrum, oenin absorbs mainly in the green part
of the spectrum, where chlorophylls don’t absorb at all.

Superposition of spectra of chlorophyll a and b with oenin

Superposition of spectra of chlorophyll a and b with oenin

http://upload.wikimedia.org/wikipedia/commons/thumb/f/f0/Spectra_Chlorophyll_
ab_oenin_%281%29.PNG/220px-Spectra_Chlorophyll_ab_oenin_%281%29.PNG

Many important natural substances are chelates. In chelates a central metal ion is
bonded to a large organic molecule, a molecule composed of carbon, hydrogen, and
other elements such as oxygen and nitrogen. One such chelate is chlorophyll, the
green pigment of plants. In chlorophyll the central ion is magnesium, and the large
organic molecule is a porphyrin. The porphyrin contains four nitrogen atoms that form
bonds to magnesium in a square planar arrangement. There are several forms of
chlorophyll. The structure of one form, chlorophyll a, is shown.

chlrphyl

http://scifun.chem.wisc.edu/chemweek/chlrphyl/chlrphyl.gif

(As you can see from the molecular structure, the “chloro” in chlorophyll does not
mean that it contains the element chlorine. The chloro portion of the word is from
the Greek chloros, which means yellowish green. The name of the element chlorine
comes from the same source. Chlorine is a yellowish green gas.)

Chlorophyll is one of the most important chelates in nature. It is capable of
channeling the energy of sunlight into chemical energy through the process of
photosynthesis. In photosynthesis, the energy absorbed by chlorophyll transforms
carbon dioxide and
water into carbohydrates and oxygen.

CO2 + H2O ——- (CH2O) + O2

(In this equation, (CH2O) is the empirical formula of carbohydrates.) The chemical
energy stored by photosynthesis in carbohydrates drives biochemical reactions in
nearly all living organisms.

In the photosynthetic reaction, carbon dioxide is reduced by water; in other words,
electrons are transferred from water to carbon dioxide. Chlorophyll assists this
transfer. When chlorophyll absorbs light energy, an electron in chlorophyll is excited
from a lower energy state to a higher energy state. In this higher energy state, this
electron is more readily transferred to another molecule. This starts a chain of
electron-transfer steps, which ends with an electron transferred to carbon dioxide.

Meanwhile, the chlorophyll which gave up an electron can accept an electron from
another molecule. This is the end of a process which starts with the removal of an
electron from water. Thus, chlorophyll is at the center of the photosynthetic
oxidation-reduction reaction between carbon dioxide and water.

Other molecules with structures similar to that of chlorophyll play important roles in
other biochemical electron-transfer (oxidation-reduction) reactions. Heme consists
of a porphyrin similar to that in chlorophyll and an iron(II) ion in the center of the
porphyrin. Heme is bright red. In the red blood cells of vertebrates, heme is bound
to proteins forming hemoglobin. Hemoglobin combines with oxygen in the lungs, gills,
or other respiratory surfaces and releases it in the tissues. In muscle cells, myoglobin,
the name given to hemoglobin in muscles, stores oxygen as an electron source for
energy-releasing oxidation-reduction reactions.

Another relative of chlorophyll is vitamin B12. Vitamin B12 contains a cobalt ion at
the center of the porphyrin. Like heme, vitamin B12 is bright red. It is essential to
digestion and nutritional absorption in animals. The exact way it functions is not
known. Because vitamin B12 is not produced by higher plants, a strictly vegetarian
diet can lead to vitamin B12 deficiency. However, it is produced by molds and
bacteria which grow on most foods.

The intense color of chlorophyll suggests that it may be useful as a commercial
pigment. In fact, chlorophyll a is a green dye (Natural Green 3) used in soaps and
cosmetics. The absorption spectrum of chlorophyll (below) shows that it absorbs
strongly in the red and blue-violet regions of the visible spectrum. Because it absorbs
red and blue-violet light, the light it reflects and transmits appears green. Commercial
pigments with structures similar to chlorophyll have been produced in a range of colors.
Some of these have slightly modified porphyrins, such as having hydrogen atoms
replaced with chlorine atoms. Others have different metal ions. For example, one
bright blue pigment has a copper(I) ion at the center of the porphyrin and is used
primarily in coloring fabrics.

http://scifun.chem.wisc.edu/chemweek/chlrphyl/clrphlsp.gif

Read Full Post »

The Colors of Life Function

Writer and Curator: Larry H. Bernstein, MD, FCAP 

2.5.1 Type 1 Copper Proteins

The Cu(II) state of this category has an intense blue color due to a thiolate ligand
to Cu(II) charge transfer, and unusual EPR properties arising from the asymmetrical
Cu site (distorted trigonal-pyramidal). The proteins all have a low molecular
mass and have, so far, rather arbitrarily been divided into sub-groups, such as
azurins, plastocyanins, pseudoazurins, amicyanins and various other blue
proteins. Of these the azurins, amicyanins, pseudo-azurins and plastocyanins
apparently have similar copper coordination by two histidine, one cysteine and
one methionine residue. Where the function of Type I copper proteins is known,
it is invariably electron transfer. As yet the names for these proteins are all trivial
and are often derived from source, function or color. The different classes are
usually discerned on the basis of their primary and tertiary structure.

The first bacterial blue proteins to be described were called azurins. Rusticyanin is
another example of a bacterial protein. It has unusual properties with a reduction
potential of 680 mV, and is functional at pH 2. The azurins have well-defined electron
-transfer functions.

The so-called pseudo-azurins differ from the azurins in the N-terminal amino acid
sequence and the optical spectra, which resemble those of plastocyanins.

The blue proteins known as plastocyanins occur in plants, blue-green and green
algae. Their electron transfer role is well defined, i.e. from the bc1 complex
(EC 1.10.2.2) to the photooxidized P-700.

Amicyanins are electron carriers between methylamine dehydrogenase and
cytochrome c, with a characteristic amino acid sequence.

Of the remaining blue proteins stellacyanin is a well- known example. Umecyanin,
plantacyanin and mavicyanin are also considered to belong to this group.
Although these proteins undergo redox reactions in vitro, their true biological
function remains unknown. Most of these proteins exhibit an unusual EPR signal
in which the copper hyperfine splitting pattern is poorly resolved. There is good
evidence that at least for stellacyanin, methionine does not function as a ligand
for copper.

2.5.2 Type 2 Copper Proteins

The copper centres in these proteins are spectroscopically consistent with square
planar or pyramidal coordination, containing oxygen and/or nitrogen ligation.
The Cu(II) is EPR active, with a ‘normal’ signal. There is no intense blue color.
This group includes the copper/zinc superoxide dismutase (EC 1.15.1.1),
dopamine b-monooxygenase (EC 1.14.17.1), galactose oxidase (EC 1.1.3.9)
and the various copper-containing amine oxidases. Some members of this last
group may also contain an organic prosthetic group, such as PQQ
(see section 10), or a modified amino-acid residue.

2.5.3 Type 3 Copper Proteins

In this group a pair of copper atoms comprise a dinuclear centre, with no EPR
activity as for single Cu’s. The best known example of an enzyme containing a
single Type 3 centre is tyrosinase (catechol oxidase, EC 1.10.3.1). This protein
contains a metal center which is a structural analogue of the dinuclear copper
center in hemocyanin (ref 31).

2.5.4 Multi-Copper Oxidases

In addition to the above, there are several proteins with catalytic activity that
contain Types 1, 2 and 3 centres in various stoichiometric ratios. These
include L-ascorbate oxidase (EC 1.10.3.3), laccase (EC 1.10.3.2) and
ceruloplasmin (ferro-oxidase, EC 1.16.3.1), the latter two having aromatic diamine
and diphenol oxidase activity. There is growing evidence that in these proteins
the Type 2 and Type 3 copper centres are juxtaposed. Recently it has been
shown that in L-ascorbate oxidase, a trinuclear copper site is present, consisting
of a type 3 copper site, very close (3.9 Å) and possibly bridged to a type 2 copper
site (ref 32). There is a view that ceruloplasmin functions as a ferro-oxidase
and the Fe(III) produced in this reaction can then oxidize the same substrates
as laccase.

2.5.5 Copper Centres in Cytochrome Oxidase

There are two copper centres that appear to be unique. Both are present in
cytochrome-c oxidase (EC 1.9.3.1). The first appears to be an isolated metal ion
and has been referred to as Cud and CuA. The second appears to be part
of a dinuclear centre with cytochrome a3. It has been referred to as Cuu,
Cua3 and CuB. At the moment the ascriptions CuA and CuB are most frequently
used; however, the recent discovery (ref 33) of a cytochrome oxidase in which
cytochrome a has been replaced by cytochrome b, leads to the recommendation
that CuB shall be referred to as Cua3.

There is a striking similarity between two of the Cu centres of N2O reductase
and CuA (ref 34, 35).

2.5.6 Molybdenum enzymes (general)

Molybdenum enzymes contain molybdenum at the catalytic center responsible
for reaction with substrate. They may be divided into those that contain
the iron-molybdenum cofactor and those that contain the pterin-molybdenum
cofactor.

2.5.7 Additional centers

If a molybdenum enzyme contains flavin, it may be called either a molybdenum
flavoprotein or a flavomolybdenum protein, as indicated above. Other centers
should be treated similarly, e.g. an iron-sulfur molybdenum protein.

2.5.8 Molybdenum enzymes containing the iron-molybdenum cofactor

The only enzymes at present known to belong to this group are the nitrogenases
(EC 1.18.6.1; and EC 1.19.6.1): see pp 89-116 in (ref 36) and pp 91-100 in (ref 37).

2.5.9 Molybdenum enzymes containing the pterin-molybdenum cofactor

These enzymes [see pp 411-415 in (ref 36) and (ref 38)] may be divided
into those in which the molybdenum bears a cyanide-labile sulfido (or thio
– see Note 1) ligand i.e. containing the S2- ligand as Mo=S) and those
lacking this ligand. The former group includes xanthine oxidase (EC 1.1.3.22),
xanthine dehydrogenase (EC 1.1.1.204), aldehyde oxidase (EC 1.2.3.1) and
purine hydroxylase (EC: see Note 2 and 3). These may be called ‘molybdenum-
containing hydroxylase’ as is widely done. Molybdenum enzymes lacking the
sulfide (thio) ligand include sulfite oxidase (EC 1.8.3.1), NAD(P)+-independent
aldehyde dehydrogenase and nitrate reductases (assimilatory and dissimilatory)
(EC 1.6.6.1-3).

2.5.10 Molybdenum enzymes containing the pterin-molybdenum cofactor

These enzymes [see pp 411-415 in (ref 36) and (ref 38)] may be divided into those
in which the molybdenum bears a cyanide-labile sulfido (or thio – see Note 1)
ligand i.e. containing the S2- ligand as Mo=S) and those lacking this ligand. The
former group includes xanthine oxidase (EC 1.1.3.22), xanthine dehydrogenase
(EC 1.1.1.204), aldehyde oxidase (EC 1.2.3.1) and purine hydroxylase. These
may be called ‘molybdenum-containing hydroxylase’ as is widely done.
Molybdenum enzymes lacking the sulfide (thio) ligand include sulfite oxidase
(EC 1.8.3.1), NAD(P)+-independent aldehyde dehydrogenase and nitrate
reductases (assimilatory and dissimilatory) (EC 1.6.6.1-3).

2.5.11 Metal-Substituted Metalloproteins

Scientists from several areas, dealing with spectroscopy and electron-transfer
mechanisms, often use metalloproteins in which a metal at the active site has
been substituted by another metal ion, like Co, Zn, Hg, Cd. Examples are zinc-
substituted cytochromes and cobalt-substituted ferredoxins.

The names for such modified proteins are easily given by using indications
like: ‘zinc-substituted ….’. In case of multi-metal proteins, where ambiguity might
arise about which metal has been substituted, one could easily add in parentheses
the name of the metal that has been replaced, such as: cobalt- substituted [Fe]
nitrogenase.

In formulae fragments or short names one could use the following notation:
[3Fe1Co-4S]2+, cytochrome c'[Fe[arrow right]CoFe], plastocyanin[Cu
[arrow right]Hg].

Ambler, R.P. (1980) in From Cyclotrons to Cytochromes (Kaplan, N.O. &
Robinson, A., eds) Academic Press, New York

Moore, G. & Pettigrew, F.(1987) Cytochromes c, Springer-Verlag, Berlin

Bartsch, R.G. (1963) in Bacterial Photosynthesis (Gest, H., San Pietro, A. &
Vernon, L.P., ed.) p. 315, Antioch Press, Yellow Springs, Ohio.

Stiefel, E.I. & Cramer, S.P. (1985) in Molybdenum Enzymes (Spiro, T.G., ed.),
Wiley-Interscience, New York, 89-116.

Smith B.E. et al. (1988), in Nitrogen Fixation Hundred Years After (Bothe,
H., de Bruijn, F.J. & Newton, W.E., ed.), Gustav Fischer, Stuttgart, New York,
91-100

Type-2 copper-containing enzymes.
MacPherson IS1, Murphy ME.
Cell Mol Life Sci. 2007 Nov;64(22):2887-99.

Type-2  Cu sites are found in all the major branches of life and are often
involved in the catalysis of oxygen species. Four type-2 Cu protein
families are selected as model systems for review: amine oxidases,
Cu monooxygenases, nitrite reductase/multicopper oxidase, and
CuZn superoxide dismutase. For each model protein, the availability
of multiple crystal structures and detailed enzymological studies provides
a detailed molecular view of the type-2 Cu site and delineation of the
mechanistic role of the Cu in biological function. Comparison of these
model proteins leads to the identification of common properties of the
Cu sites and insight into the evolution of the trinuclear active site found
in multicopper oxidases.

Copper proteins and copper enzymes.
Cass AE, Hill HA.
Ciba Found Symp. 1980;79:71-91.
http://www.chm.bris.ac.uk/motm/caeruloplasmin/copper_proteins/t1.htm

The copper proteins that function in homeostasis, electron transport, dioxygen
transport and oxidation are discussed. Particular emphasis is placed on the
role of the ligands, their type and disposition which, in conjunction with other
residues in the active site, determine the role of the copper ion. It is proposed that
copper proteins can be considered in four groups. Those in Group I contain a
single copper ion in an approximately tetrahedral environment with nitrogen and
sulphur-containing ligands. Group II proteins have a single copper ion in a square-
planar-like arrangement. Group III proteins have two copper ions in close
proximity. Group IV consists of multi-opper proteins, composed of sites
representative of the other three groups.

Such centers owe their name to the intense blue coloration of the corresponding
Cu(II) proteins. The color is particularly distinctive since the metal centers are
so optically diluted in these metalloenzymes that only intense absorption in the
visible region, resulting from symmetry allowed electronic transitions, can give
rise to conspicuous colors. In contrast, the comparatively pale blue color of normal
Cu(II)) is the result of forbidden electronic transitions between d-orbitals of
different symmetry; in Cu2+(aq) this gives a molar extinction coefficient of
10 M-1cm-1 from a broad absorption between 10,000 cm-1 and 15,000 cm-1
compared to about 3000 M-1cm-1 observed for blue Cu(II) centers.  For the
T1 centers the intense absorption is attributed to a ligand-to-metal charge
transfer between the Cu2+ and a bonded cysteinate ligand. Typically, as in
azurin or plastocyanin this occurs around 16,000 cm-1. Ceruloplasmin has
three T1 centers, and the blue absorption is at 16,400 cm-1 (610nm).

Plastocyanine geometry

around the copper Crystal structures show a very irregular ‘tetrahedral’ coordination
with two sulphurs from methionine and cysteinate, and two histidine nitrogens.
However a comparison of azurin with plastocyanin shows that the geometry
is in some ways closer to a trigonal bipyramid, with and without one extra apical
ligand, so that azurin has a weakly bound glutamine oxygen, and plastocyanine
does not. The T1 coppers in caruloplasmin are in plastocyanine-type domains.
Each of these are coordinated to two histidines and a cysteine, in two of the T1
domains there is also a methionine residue, the third T1 domain has a leucine
residue which may only have a van der Waals type contact with the copper.

T1 copper centers are functional in the reversible electron transfer:

Cu2+ + e-   =   Cu+

The strongly distorted geometry represents a compromise (entactic-state
situation) between d10 Cu(I), with its preferred tetrahedral or trigonal
coordination through soft sulfur ligands, and d9 Cu(II) with preferential
square planar or square pyramidal geometry and nitrogen ligand
coordination.   This irregular, high energy arrangement at the metal
center resembles the transition-state geometry between the tetrahedral
and square planar equilibrium configurations of the two oxidation states
involved and permits enhanced rates of electron transfer. The potential
range for proteins with T1 copper centers runs from 180 mV in
stellacyanin to 680 mV in rusticyanin.

Zinc proteins: enzymes, storage proteins, transcription factors, and replication
proteins.
Coleman JE.
Annu Rev Biochem. 1992;61:897-946.

In the past five years there has been a great expansion in our knowledge of
the role of zinc in the structure and function of proteins. Not only is zinc
required for essential catalytic functions in enzymes (more than 300 are known
at present), but also it stabilizes and even induces the folding of protein
subdomains. The latter functions have been most dramatically illustrated
by the discovery of the essential role of zinc in the folding of the DNA-binding
domains of eukaryotic transcription factors, including the zinc
finger transcription factors, the large family of hormone receptor proteins,
and the zinc cluster transcription factors from yeasts. Similar functions are
highly probable for the zinc found in the RNA polymerases and the zinc-
containing accessory proteins involved in nucleic acid replication. The rapid
increase in the number and nature of the proteins in which zinc functions
is not unexpected since zinc is the second most abundant trace metal found in
eukaryotic organisms, second only to iron. If one subtracts the amount of iron
found in hemoglobin, zinc becomes the most abundant trace metal found
in the human body.

Zinc Coordination Spheres in Protein Structures
ACS ChemWorx
Mikko Laitaoja , Jarkko Valjakka , and Janne Jänis
Inorg. Chem., 2013, 52 (19), pp 10983–10991
http://dx.doi.org:/10.1021/ic401072d
Sept 23, 2013

Synopsis
A statistical analysis in terms of zinc coordinating amino acids, metal-to-ligand
bond lengths, coordination number, and structural classification was performed,
revealing coordination spheres from classical tetrahedral cysteine/histidine binding
sites to more complex binuclear sites with carboxylated lysine residues. According
to the results, coordination spheres of hundreds of crystal structures in the PDB
could be misinterpreted due to symmetry-related molecules or missing electron
densities for ligands.

Protein-folding location can regulate manganese-binding versus copper- or
zinc-binding.
Tottey S, Waldron KJ, Firbank SJ, Reale B, Bessant C, Sato K, Cheek TR, et al.
Nature. 2008 Oct 23;455(7216):1138-42. http://dx.doi.org:/10.1038/nature07340

Metals are needed by at least one-quarter of all proteins. Although metallo-
chaperones insert the correct metal into some proteins, they have not been
found for the vast majority, and the view is that most metalloproteins acquire
their metals directly from cellular pools. However, some metals form more
stable complexes with proteins than do others. For instance, as described
in the Irving-Williams series, Cu(2+) and Zn(2+) typically form more stable
complexes than Mn(2+). Thus it is unclear what cellular mechanisms manage
metal acquisition by most nascent proteins. To investigate this question, we
identified the most abundant Cu(2+)-protein, CucA (Cu(2+)-cupin A), and the
most abundant Mn(2+)-protein, MncA (Mn(2+)-cupin A), in the periplasm of
the cyanobacterium Synechocystis PCC 6803. Each of these newly identified
proteins binds its respective metal via identical  ligands within a cupin fold.
Consistent with the Irving-Williams series, MncA only binds Mn(2+) after
folding in solutions containing at least a 10(4) times molar excess of Mn(2+)
over Cu(2+) or Zn(2+). However once MncA has bound Mn(2+), the metal
does not exchange with Cu(2+). MncA and CucA have signal peptides for
different export pathways into the periplasm, Tat and Sec respectively. Export
by the Tat pathway allows MncA to fold in the cytoplasm, which contains only
tightly bound copper or Zn(2+) (refs 10-12) but micromolar Mn(2+) (ref. 13). In
contrast, CucA folds in the periplasm to acquire Cu(2+). These results reveal
a mechanism whereby the compartment in which a protein folds overrides its
binding preference to control its metal content. They explain why the cytoplasm
must contain only tightly bound and buffered copper and Zn(2+).

Predicting copper-, iron-, and zinc-binding proteins in pathogenic species of the
Paracoccidioides genus
GB Tristão, L do Prado Assunção, LPA dos Santos, CL Borges, MG Silva-Bailão,
CM de Almeida Soares, G Cavallaro and AM Bailão*
Front. Microbiol., 9 Jan 2015 http://dx.doi.org:/10.3389/fmicb.2014.00761

Approximately one-third of all proteins have been estimated to contain at least
one metal cofactor, and these proteins are referred to as metalloproteins. These
represent one of the most diverse classes of proteins, containing metal ions that
bind to specific sites to perform catalytic, regulatory and structural functions.
Bioinformatic tools have been developed to predict metalloproteins encoded by
an organism based only on its genome sequence. Its function and the type of
metal binder can also be predicted via a bioinformatics approach.  Paracoccidioides
complex includes termodimorphic pathogenic fungi that are found as saprobic
mycelia in the environment and as yeast, the parasitic form, in host tissues. They
are the etiologic agents of Paracoccidioidomycosis, a prevalent systemic mycosis
in Latin America. Many metalloproteins are important for the virulence of several
pathogenic microorganisms. Accordingly, the present work aimed to predict the
copper, iron and zinc proteins encoded by the genomes of three phylogenetic species
of Paracoccidioides (Pb01, Pb03, andPb18). The metalloproteins were identified
using bioinformatics approaches based on structure, annotation and domains. Cu-,
Fe-, and Zn-binding proteins represent 7% of the total proteins encoded by
Paracoccidioides spp. genomes. Zinc proteins were the most abundant metallo-
proteins, representing 5.7% of the fungus proteome, whereas copper and iron
proteins represent 0.3 and 1.2%, respectively. Functional classification revealed that
metalloproteins are related to many cellular processes. Furthermore, it was observed
that many of these metalloproteins serve as virulence factors in the biology of the
fungus. Thus, it is concluded that the Cu, Fe, and Zn metalloproteomes of the
Paracoccidioides spp. are of the utmost importance for the biology and virulence
of these particular human pathogens.

Zinc finger proteins: new insights into structural and functional diversity
John H Laity, Brian M Lee, Peter E Wright
Current Opinion in Structural Biology Feb 2001; 11(1): 39–46
http://epigenie.com/key-epigenetic-players/chromatin-modifying-and-dna-
binding-proteins/zinc-finger-proteins/

Zinc finger proteins are among the most abundant proteins in eukaryotic genomes.
Their functions are extraordinarily diverse and include DNA recognition, RNA
packaging, transcriptional activation, regulation of apoptosis, protein folding
and assembly, and lipid binding. Zinc finger structures are as diverse as their
functions. Structures have recently been reported for many new zinc finger
domains with novel topologies, providing important insights into structure/function
relationships. In addition, new structural studies of proteins containing the
classical Cys2His2 zinc finger motif have led to novel insights into mechanisms
of DNA binding and to a better understanding of their broader functions in
transcriptional regulation.

Zinc Finger Proteins

Zinc finger (ZnF) proteins are a massive, diverse family of proteins that serve a
wide variety of biological functions. Due to their diversity, it is difficult to come up
with a simple definition of what unites all ZnF proteins; however, the most common
approach is to define them as all small, functional domains that require coordination
by at least one zinc ion (Laity et al., 2001). The zinc ion serves to stabilize the
integration of the protein itself, and is generally not involved in binding targets.
The “finger” refers to the secondary structures (α-helix and β-sheet) that are
held together by the Zn ion. Zinc finger containing domains typically serve
as interactors, binding DNA, RNA, proteins or small molecules (Laity et al., 2001).

ZnF Protein Families

Cys2His2 was the first domain discovered (also known as Krüppel-type). It was
initially discovered as a repeating domain in the IIIA transcription factor in
Xenopus laevis (Brown et al., 1985; Miller et al., 1985). IIIA has nine repeats
of the 30 amino acids that make up the Cys2His2 domain. Each domain forms
a left-handed ββα secondary structure, and coordinates a Zn ion between
two cysteines on the β-sheet hairpin and two histidines in the α-helix, hence
the name Cys2His2 (Lee et al., 1989). These resides are highly conserved,
as well as a general hydrophobic core that allows the helix to form. The other
residues can show great sequence diversity (Michael et al., 1992). Cys2His2
zinc fingers that bind DNA tend to have 2-4 tandem domains as part of a
larger protein. The residues of the alpha helices form specific contacts with a
specific DNA sequence motif by “reading” the nucleotides in major groove
of DNA (Elrod-Erickson et al., 1996; Pavletich and Pabo, 1991). Cys2His2
proteins are the biggest group of transcription factors in most species. Non-
DNA binding proteins can have much more flexible tertiary structure.
Examples of Cys2His2 proteins include the Inhibitor of Apoptosis (IAP) family
of proteins and the CTFC transcription factor.

Treble clef fingers are a very diverse group of ZnF protiens both in terms of
structure and function. What makes them a family is a shared fold at their core
that looks a little like a musical treble clef, especially if you squint (Grishin,
2001). Most treble clef finger motifs have a β hairpin, a variable loop region,
a β hairpin, and an α helix. The “knuckle” of the β hairpin and the α helix contain
the Cys-x-x-Cys sequence necessary to coordinate the Zn ion. Treble clef
fingers often form the core of protein structures, for example the L24E and
S14 ribosomal proteins and the RING finger family.

Zinc ribbons are a little less structurally complex than the other two major groups.
Zinc ribbons contain two zinc knuckles, often β hairpins, coordinating a zinc ion via
a two Cys residures separated by 2-4 other residues on one knuckle, and a Cys-x-x-
Cys on the other (Hahn and Roberts, 2000). Examples of zinc ribbon-containing
proteins include the basal transcription factors TFIIS and TFIIB that for a complex
with RNAPII to bind DNA, and the Npl4 nuclear core protein that uses a zinc ribbon
to bind ubiquitin (Alam et al., 2004). Cys2His2, treble clef fingers, and zinc ribbons
form the majority of zinc fingers, but there are several other smaller groups that
don’t fit neatly into these three. Green fluorescent protein as a marker for gene
expression.

Metallothionein proteins expression, copper and zinc concentrations, and lipid
peroxidation level in a rodent model for amyotrophic lateral sclerosis
E Tokuda, Shin-Ichi Ono,  K Ishige, A Naganuma, Y Ito, T Suzuki
Toxicology Jan 2007; 229(1–2): 33–41

It has been hypothesized that copper-mediated oxidative stress contributes to the
pathogenesis of familial amyotrophic lateral sclerosis (ALS), a fatal motor neuron
disease in humans. To verify this hypothesis, we examined the copper and zinc
concentrations and the amounts of lipid peroxides, together with that of the
expression of metallothionein (MT) isoforms in a mouse model [superoxide
dismutase1 transgenic (SOD1 Tg) mouse] of ALS. The expression of MT-I and
MT-II (MT-I/II) isoforms were measured together with Western blotting, copper
level, and lipid peroxides amounts increased in an age-dependent manner in the
spinal cord, the region responsible for motor paralysis. A significant increase was
already seen as early as 8-week-old SOD1 Tg mice, at which time the mice had not
yet exhibited motor paralysis, and showed a further increase at 16 weeks of age,
when paralysis was evident. Inversely, the spinal zinc level had significantly
decreased at both 8 and 16 weeks of age. The third isoform, the MT-III level,
remained at the same level as an 8-week-old wild-type mouse, finally increasing
to a significant level at 16 weeks of age. It has been believed that a mutant SOD1
protein, encoded by a mutant SOD1, gains a novel cytotoxic function while
maintaining its original enzymatic activity, and causes motor neuron death
(gain-of-toxic function). Copper-mediated oxidative stress seems to be a probable
underlying pathogenesis of gain-of-toxic function. Taking the above current
concepts and the classic functions of MT into account, MTs could have a disease
modifying property: the MT-I/II isoform for attenuating the gain-of-toxic function
at the early stage of the disease, and the MT-III isoform at an advanced stage.

Prion protein expression level alters regional copper, iron and zinc content in
the mouse brain
MJ Pushie,  IJ Pickering, GR Martin, S Tsutsui, FR Jirik and GN George
Metallomics, 2011,3, 206-214 http://dx.doi.org:/10.1039/C0MT00037J

The central role of the prion protein (PrP) in a family of fatal neurodegenerate
diseases has garnered considerable research interest over the past two decades.
Moreover, the role of PrP in neuronal development, as well as its apparent role
in metal homeostasis, is increasingly of interest. The host-encoded form of the
prion protein (PrPC) binds multiple copper atoms via its N-terminal domain
and can influence brain copper and iron levels. The importance of PrPC to the
regulation of brain metal homeostasis and metal distribution, however, is not
fully understood. We therefore employed synchrotron-based X-ray fluorescence
imaging to map the level and distributions of several key metals in the brains of
mice that express different levels of PrPC. Brain sections from wild-type, prion
gene knockout (Prnp−/−) and PrPC over-expressing mice revealed striking
variation in the levels of iron, copper, and even zinc in specific brain regions as
a function of PrPC expression. Our results indicate that one important function
of PrPC may be to regulate the amount and distribution of specific metals within
the central nervous system. This raises the possibility that PrPC levels, or its
activity, might regulate the progression of diseases in which altered metal
homeostasis is thought to play a pathogenic role such as Alzheimer’s,
Parkinson’s and Wilson’s diseases and disorders such as hemochromatosis.

Zinc & Copper Imbalances: Immense Biochemical Implications
Mar 27, 2013 by Michael McEvoy
http://metabolichealing.com/zinc-copper-imbalances-immense-biochemical-
implications/

The status of zinc and copper levels may have profound implications for
many people. Much has been written about the significance of these two
trace elements for many, many years. Many health conditions may be
directly caused by abnormal zinc and copper levels.

With all of the recent attention given to methylation status, gene mutations,
MTHFR, and the associated neurological and mental/behavioral disorders
that may ensue, zinc and copper status remains a pivotal ratio in these regards.

While zinc toxicity and copper deficiency are possible, the subject of this
article is on the more common imbalance: copper toxicity and zinc deficiency.

The Physiological Roles Of Zinc & Copper

Zinc and copper are antagonists. The balance between these two trace
elements is an example of the effects of biological dualism. While zinc
toxicity is possible, far more common is zinc deficiency and copper toxicity.
Both zinc and copper play essential roles in the body, and there can be a
number of causes for why imbalances ensue.

It may be easier to identify the roles that zinc doesn’t play in the body,
than the roles it does play. Zinc is an essential trace element that activates
several hundred enzymatic reactions. These reactions are fundamental
to life and biological activity. Some of the activities that zinc are involved in:

  • DNA & RNA synthesis
  • Gene expression
  • Nervous system function
  • Immune function & immune signaling such as cell
    apoptosis
  • Neuronal transmission
  • Brain function
  • Zinc possesses powerful anabolic activities in the cells
  • Formation of zinc proteins known as “zinc fingers”
  • Zinc is essential for blood clotting and platelet formation
  • Zinc is involved in Vitamin A synthesis
  • Folate is made available through zinc enzyme reactions
  • Along with copper, Zinc makes up the antioxidant
    enzyme
    system, ZnCu superoxide dismutase
  • Steroidal hormone synthesis
  • Growth & development of children
  • Testosterone and semen formation
  • The highest concentration of zinc is found in the
    male prostate gland

Copper is an essential trace element serving many important functions
as well. However, copper is well documented to induce several toxic effects
in the body, when elevated. Because copper is a pro-oxidant when free and
unbound, it can quickly generate free radicals.

The major sources for copper toxicity are: exposure to industrial forms
of copper such as copper pipes, copper cookware, birth control, exposure
to copper-based fungicides. Diets high in copper and low in zinc may play
a role in copper toxicity. Pyrrole disorder, which causes depletion of zinc,
may result in elevated levels of copper.

Some of the essential roles copper plays in the body:

  • Connective tissue formation
  • ATP synthesis
  • Iron metabolism
  • Brain health via neurotransmitter synthesis
  • Gene transcription
  • Synthesis of the antioxidant superoxide dismutase
  • Skin pigmentation
  • Nerve tissue: myelin sheath formation
  • Copper tends to rise when estrogen is dominant

Perhaps one of the first reports that zinc and copper imbalances play
a role in human health and disease was their detection in mental
disorders made by Carl Pfeiffer, MD, PhD. Dr. Pfeiffer identified a
condition known as pyrrole disorder, sometimes referred to as
pyrroluria or “mauve factor”.

As it turns out, pyrrole disorder is a major biochemical imbalance
in many people with chronic illnesses such as chronic Lyme disease,
autism, schizophrenia, depression, bi-polar, and chronic fatigue
syndrome. Pyrroles are a byproduct of hemoglobin synthesis.
Apparently, some individuals are more predisposed towards producing
higher amounts of pyrroles. When pyrroles are excessive, they irreversibly
bind to zinc and vitamin B6, causing their excretion. Consequently,
it is common that once zinc levels become depleted, copper levels tend to rise.

Copper Toxicity

Problems associated with copper toxicity include: pyrrole disorder,
estrogen dominance, schizophrenia, depression, anxiety disorder,
chronic fatigue, migraines, liver toxicity, thyroid conditions, chronic
candida yeast infections, PMS, to name a few. Some research has
even implicated copper toxicity with Alzheimer’s Disease and with
cardiovascular disease. Perhaps one of the primary mechanisms
through which copper toxicity can damage tissues is through its
initiation of oxidative stress and free radical formation. Free copper
ions that are not bound to copper proteins such as ceruloplasmin,
are pro-oxidants, and are highly reactive.

Empirical research from clinicians, indicates that there are different
types of copper imbalances. For example, if there is a lot of free,
unbound copper present, this may cause a situation of nutritive
copper deficiency. Another copper imbalance is when high pyrroles
depress zinc levels, and copper levels concomintantly rise. If high
pyrroles are present, B6 will also be lost in high amounts. In a general
but very real sense, all forms of copper excess will affect zinc status,
due to the dualistic nature of zinc and copper.

Copper & Estrogen

It has been known for many years that copper can cause a rise in
estrogen, and conversely estrogen may raise copper. Estrogen
dominance has been extensively studied in its role in breast
cancer development. One possible, critical role that can cause
estrogen to become carcinogenic, is through its oxidation induced by
copper. 
Once oxidized, estrogen forms volatile hydroxyl radicals and
the associated DNA damage and “mutagenesis”.

Zinc Deficiency

As mentioned previously, pyrrole disorder will directly depress
zinc status, causing high levels of its excretion. When zinc is
lost, copper rises. Because of their essential roles in neuro-
transmitter synthesis, low zinc and high copper levels can
directly effect cognition, behavior and thought processes.
Zinc has been studied in biochemical reactions involving
calcium-driven, synaptic neurotransmission, as well as in
glutamate/GABA balance and with limbic brain function.

Zinc & Reproduction

Zinc is essential for steroidal hormone synthesis, and is a
well known catalyst for testosterone synthesis, as well as
leutinizing hormone. Zinc has demonstrated its ability to
prevent miscarriage and toxicity during pregnancy. The male
prostate gland reportedly contains the highest concentration
of zinc in the body.

Zinc & Brain Function

Much attention has been given to excitotoxicity, such as the
effects induced by MSG (monosodium glutamtate). Excess
stimulation of the excitatory neurotransmitter glutamate,
may cause severe physical and psychological reactions in
certain individuals. Zinc has been studied for its ability to
enhance GABA 
(glutamate’s antagonistic neurotransmitter)
activity and to suppress excess glutamate.

Studies on mice demonstrated that when depleted of zinc
for two weeks, the mice developed seizures, most likely due
to GABA deficiencies and glutamate excess.

There is an emerging body of evidence that demonstrates
that Alzheimer’s disease may involve copper toxicity and
zinc deficiency. Not only can excess copper cause zinc
depletion, but so can excess lead.

The hippocampus, a major part of the limbic brain, records
memories and is responsible for processing meaningful
experiences. Numerous studies site that if hippocampal
cells are deprived of zinc, the hippocampal cells die. In
addition to hippocampus cell death induced by zinc
deprivation, the amygdala, the other major limbic gland
experiences cell death as well, when deprived of zinc.

Green Fluorescent Protein

Chalfie M, Tu Y, Euskirchen G, Ward WW, Prasher DC.
Science. 1994 Feb 11;263(5148):802-5.
http://www.ncbi.nlm.nih.gov/pubmed/8303295

A complementary DNA for the Aequorea victoria green fluorescent protein (GFP)
produces a fluorescent product when expressed in prokaryotic (Escherichia coli)
or eukaryotic (Caenorhabditis elegans) cells. Because exogenous substrates and
cofactors are not required for this fluorescence, GFP expression can be used
to monitor gene expression and protein localization in living organisms.

http://en.wikipedia.org/wiki/Green_fluorescent_protein

The green fluorescent protein (GFP) is a protein composed of 238 amino acid
residues (26.9 kDa) that exhibits bright green fluorescence when exposed
to light in the blue to ultraviolet range. Although many other marine organisms
have similar green fluorescent proteins, GFP traditionally refers to the protein
first isolated from the jellyfish Aequorea victoria. The GFP from A. victoria
has a major excitation peak at a wavelength of 395 nm and a minor one at
475 nm. Its emission peak is at 509 nm, which is in the lower green portion
of the visible spectrum. The fluorescence quantum yield (QY) of GFP is 0.79.
The GFP from the sea pansy (Renilla reniformis) has a single major excitation
peak at 498 nm.

In cell and molecular biology, the GFP gene is frequently used as a reporter of
expression. In modified forms it has been used to make biosensors, and many
animals have been created that express GFP as a proof-of-concept that a gene
can be expressed throughout a given organism. The GFP gene can be introduced
into organisms and maintained in their genome through breeding, injection with a
viral vector, or cell transformation. To date, the GFP gene has been introduced
and expressed in many Bacteria, Yeast and other Fungi, fish (such as zebrafish),
plant, fly, and mammalian cells, including human. Martin Chalfie, Osamu Shimomura,
and Roger Y. Tsien were awarded the 2008 Nobel Prize in Chemistry on 10 October
2008 for their discovery and development of the green fluorescent protein.

http://www.conncoll.edu/ccacad/zimmer/GFP-ww/GFP-1.htm

In Aequorea victoria a protein called aequorin releases blue light upon binding
with calcium. This blue light is then totally absorbed by the GFP, which in turn
gives off the green light as in the animation below.

In 1994 GFP was cloned. Now GFP is found in laboratories all over the world where
it is used in every conceivable plant and animal. Flatworms, algae, E. coli and pigs
have all been made to fluoresce with GFP.

The importance of GFP was recognized in 2008 when the Nobel Committee awarded
Osamu Shimomura, Marty Chalfie and Roger Tsien the Chemistry Nobel Prize ”
for the discovery and development of the green fluorescent protein, GFP.”

Why is it so popular? Well, I like to think of GFP as the microscope of the twenty-
first century. Using GFP we can see when proteins are made, and where they can go.
This is done by joining the GFP gene to the gene of the protein of interest so that
when the protein is made it will have GFP hanging off it. Since GFP fluoresces, one
can shine light at the cell and wait for the distinctive green fluorescence associated
with GFP to appear.

A variant of yellow fluorescent protein with fast and efficient maturation for
cell-biological applications
T Nagai, K Ibata, E Sun Park, M Kubota, K Mikoshiba & A Miyawaki
Nature Biotechnology 20, 87 – 90 (2002)  http://dx.doi.org:/10.1038/nbt0102-87

The green fluorescent protein (GFP) from the jellyfish Aequorea victoria
has provided a myriad of applications for biological systems. Over the last
several years, mutagenesis studies have improved folding properties of GFP.
However, slow maturation is still a big obstacle to the use of GFP variants for
visualization. These problems are exacerbated when GFP variants are expressed
at 37°C and/or targeted to certain organelles. Thus, obtaining GFP variants that
mature more efficiently is crucial for the development of expanded research
applications. Among Aequorea GFP variants, yellow fluorescent proteins (YFPs)
are relatively acid-sensitive,and uniquely quenched by chloride ion (Cl−)3. For
YFP to be fully and stably fluorescent, mutations that decrease the sensitivity
to both pH and Cl− are desired. Here we describe the development of an
improved version of YFP named “Venus”. Venus contains a novel mutation,
F46L, which at 37°C greatly accelerates oxidation of the chromophore, the rate-
limiting step of maturation. As a result of other mutations, F64L/M153T/
V163A/S175G, Venus folds well and is relatively tolerant of exposure
to acidosis and Cl−. We succeeded in efficiently targeting a neuropeptide
Y-Venus fusion protein to the dense-core granules of PC12 cells. Its secretion
was readily monitored by measuring release of fluorescence into the medium.
The use of Venus as an acceptor allowed early detection of reliable signals of
fluorescence resonance energy transfer (FRET) for Ca2+ measurements in brain
slices. With the improved speed and efficiency of maturation and the increased
resistance to environment, Venus will enable fluorescent labelings that were not
possible before.

Rhodopsin-like Protein from the Purple Membrane of Halobacterium halobium
DIETER OESTERHELT &  WALTHER STOECKENIUS
Nature New Biology 29 Sep 1971; 233, 149-152  | http://dx.doi.org:/10.1038/
newbio233149a0

HALOPHILIC bacteria require high concentrations of sodium chloride and lower
concentrations of KCl and MgCl2 for growth. The cell membrane dissociates into
fragments of varying size when the salt is removed1. One characteristic fragment—
termed the “purple membrane” because of its characteristic deep purple colour—
has been isolated in relatively pure form from Halobacterium halobium. We can
now show that the purple colour is due to retinal bound to an opsin-like protein,
the only protein present in this membrane fragment.

References

Stoeckenius, W. , and Rowen, R. , J. Cell Biol., 34, 365 (1967).

Stoeckenius, W. , and Kunau, W. H. , J. Cell Biol., 38, 337 (1968).

Blaurock, A. E. , and Stoeckenius, W. , Nature New Biology, 233, 152 (1971).

Sehgal, S. N. , and Gibbons, N. E. , Canad. J. Microbiol., 6, 165 (1960).

Kelly, M. , Norgård, S. , and Liaach-Jensen, S. , Acta Chem. Scand., 2A, 2169 (1970).

Shapiro, A. L. , Vinnela, E. , and Maizel, jun., J. V. , Biochem. Biophys. Res.
Commun., 28, 815 (1967).

The monomerization of the Purple protein, a member of the GFP-family
Corning, Brooke

Green fluorescent protein (GFP) has been used extensively since its discovery
in the 1960s to report and visualize gene expression. For years it has been the only
known naturally occurring fluorescent pigment that is encoded by a single gene,
making it extremely useful in various fields of biology, because the expression of
this gene directly leads to the appearance of the fluorescent green color. Recently,
however, many more proteins with similar properties to GFP, and available in a
variety of colors, have been isolated from the class of marine organisms called
Anthozoa, which includes the corals. This increase in the availability of colored
proteins in GFP family in turn has expanded the number of available biotech-
nology applications. However, some of these newly discovered GFP-like
proteins do not have wild-type forms that readily allow for the creation of
fusion proteins, particularly because of oligomerization. It is widely accepted
that almost all members of the GFP-family form dimers or tetramers in their
functional forms.

This study investigates a GFP-ike protein, Purple, isolated from two species,
Galaxea fascicularis and Montipora efflorescens. Purple protein forms oligomers
when expressed, which would then interfere with the normal expression of a  protein
to be tagged in gene fusion experiments. We selectively mutated 3 amino acids,
which we believed were responsible for oligomerization in Purple. These 3
residues were chosen based on sequence similarities to a very similar protein,
a mutant form of the Rtms5 chromoprotein from Montipora efflorescens. While
we had hoped that the resulting triple-mutant Purple protein would form
monomers in vivo while retaining its purple coloration, this turned out to
be incorrect. The resulting mutants had lost their ability to turn purple. However,
we also determined that we had successfully changed the oligomerization
state of Purple by examining the relative molecular mass of one our
mutant proteins, which turned out to be half the size of the original
purple protein. It is possible that by adding additional mutations in
the future, the original spectral properties could be recovered. If
successful, this would further expand the utility of the GFP family.

Rhodopsin, also known as visual purple, from Ancient Greek ῥόδον
(rhódon, “rose”), due to its pinkish color, and ὄψις (ópsis, “sight”), is
a light-sensitive receptor protein. It is a biological pigment in photo-
receptor cells of the retina. Rhodopsin is the primary pigment found
in rod photoreceptors. Rhodopsins belong to the G-protein-coupled
receptor (GPCR) family. They are extremely sensitive to light, enabling
vision in low-light conditions. Exposed to light, the pigment
immediately photobleaches, and it takes about 45 minutes to regenerate
fully in humans. Its discovery was reported by German physiologist
Franz Christian Boll in 1876.

Read Full Post »

The Life and Work of Allan Wilson

Curator: Larry H. Bernstein, MD, FCAP

 

Allan Charles Wilson (18 October 1934 – 21 July 1991) was a Professor of Biochemistry at the University of California, Berkeley, a pioneer in the use of molecular approaches to understand evolutionary change and reconstruct phylogenies, and a revolutionary contributor to the study of human evolution. He was one of the most controversial figures in post-war biology; his work attracted a great deal of attention both from within and outside the academic world. He is the only New Zealander to have won the MacArthur Fellowship.

He is best known for experimental demonstration of the concept of the molecular clock (with his doctoral student Vincent Sarich), which was theoretically postulated by Linus Pauling and Emile Zuckerkandl, revolutionary insights into the nature of the molecular anthropology of higher primates and human evolution, called Mitochondrial Eve hypothesis (with his doctoral students Rebecca L. Cann and Mark Stoneking).

Allan Wilson was born in Ngaruawahia, New Zealand, and raised on his family’s rural dairy farm at Helvetia, Pukekohe, about twenty miles south of Auckland. At his local Sunday School, the vicar’s wife was impressed by young Allan’s interest in evolution and encouraged Allan’s mother to enroll him at the elite King’s College secondary school in Auckland. There he excelled in mathematics, chemistry, and sports.

Wilson already had an interest in evolution and biochemistry, but intended to be the first in his family to attend university by pursuing studies in agriculture and animal science. Wilson met Professor Campbell Percy McMeekan, a New Zealand pioneer in animal science, who suggested that Wilson attend the University of Otago in southern New Zealand to further his study in biochemistry rather than veterinary science. Wilson gained a BSc from the University of Otago in 1955, majoring in both zoology and biochemistry.

The bird physiologist Donald S. Farner met Wilson as an undergraduate at Otago and invited him to Washington State University at Pullman as his graduate student. Wilson obliged and completed a master’s degree in zoology at WSU under Farner in 1957, where he worked on the effects of photoperiod on the physiology of birds.

Wilson then moved to the University of California, Berkeley, to pursue his doctoral research. At the time the family thought Allan would only be gone two years. Instead, Wilson remained in the United States, gaining his PhD at Berkeley in 1961 under the direction of biochemist Arthur Pardee for work on the regulation of flavin biosynthesis in bacteria. From 1961 to 1964, Wilson studied as a post-doc under biochemist Nathan O. Kaplan at Brandeis University in Waltham, Massachusetts. In Kaplan’s lab, working with lactate and malate dehydrogenases, Wilson was first introduced to the nascent field of molecular evolution. Nate Kaplan was one of the very earliest pioneers to address phylogenetic problems with evidence from protein molecules, an approach that Wilson later famously applied to human evolution and primate relationships. After Brandeis, Wilson returned to Berkeley where he set up his own lab in the Biochemistry department, remaining there for the rest of his life.

Wilson joined the UC Berkeley faculty of biochemistry in 1964, and was promoted to full professor in 1972. His first major scientific contribution was published as Immunological Time-Scale For Hominid Evolution in the journal Science in December 1967. With his student Vincent Sarich, he showed that evolutionary relationships of the human species with other primates, in particular the Great Apes (chimpanzees, gorillas, and orangutans), could be inferred from molecular evidence obtained from living species, rather than solely from fossils of extinct creatures.

Their microcomplement fixation method (see complement system) measured the strength of the immune reaction between an antigen (serum albumin) from one species and an antibody raised against the same antigen in another species. The strength of the antibody-antigen reaction was known to be stronger between more closely related species: their innovation was to measure it quantitatively among many species pairs as an “immunological distance”. When these distances were plotted against the divergence times of species pair with well-established evolutionary histories, the data showed that the molecular difference increased linearly with time, in what was termed a “molecular clock”. Given this calibration curve, the time of divergence between species pairs with unknown or uncertain fossil histories could be inferred. Most controversially, their data suggested that divergence times between humans, chimpanzees, and gorillas were on the order of 3~5 million years, far less than the estimates of 9~30 million years accepted by conventional paleoanthropologists from fossil hominids such as Ramapithecus. This ‘recent origin’ theory of human/ape divergence remained controversial until the discovery of the “Lucy” fossils in 1974.

Wilson and another PhD student Mary-Claire King subsequently compared several lines of genetic evidence (immunology, amino acid differences, and protein electrophoresis) on the divergence of humans and chimpanzees, and showed that all methods agreed that the two species were >99% similar.[4][19] Given the large organismal differences between the two species in the absence of large genetic differences, King and Wilson argued that it was not structural gene differences that were responsible for species differences, but gene regulation of those differences, that is, the timing and manner in which near-identical gene products are assembled during embryology and development. In combination with the “molecular clock” hypothesis, this contrasted sharply with the accepted view that larger or smaller organismal differences were due to large or smaller rates of genetic divergence.

In the early 1980s, Wilson further refined traditional anthropological thinking with his work with PhD students Rebecca Cann and Mark Stoneking on the so-called “Mitochondrial Eve” hypothesis.[20] In his efforts to identify informative genetic markers for tracking human evolutionary history, he focused on mitochondrial DNA (mtDNA) — genes that are found in mitochondria in the cytoplasm of the cell outside the nucleus. Because of its location in the cytoplasm, mtDNA is passed exclusively from mother to child, the father making no contribution, and in the absence of genetic recombination defines female lineages over evolutionary timescales. Because it also mutates rapidly, it is possible to measure the small genetic differences between individual within species by restriction endonuclease gene mapping. Wilson, Cann, and Stoneking measured differences among many individuals from different human continental groups, and found that humans from Africa showed the greatest inter-individual differences, consistent with an African origin of the human species (the so-called “Out of Africa” hypothesis). The data further indicated that all living humans shared a common maternal ancestor, who lived in Africa only a few hundreds of thousands of years ago.

This common ancestor became widely known in the media and popular culture as the Mitochondrial Eve. This had the unfortunate and erroneous implication that only a single female lived at that time, when in fact the occurrence of a coalescent ancestor is a necessary consequence of population genetic theory, and the Mitochondrial Eve would have been only one of many humans (male and female) alive at that time.[2][3] This finding was, like his earlier results, not readily accepted by anthropologists. Conventional hypothesis was that various human continental groups had evolved from diverse ancestors, over several million of years since divergence from chimpanzees. The mtDNA data, however, strongly suggested that all humans descended from a common, quite recent, African mother.

Wilson became ill with leukemia, and after a bone marrow transplant, died on Sunday, 21 July 1991, at the Fred Hutchinson Memorial Cancer Research Center in Seattle. He had been scheduled to give the keynote address at an international conference the same day. He was 56, at the height of his scientific recognition and powers.

Wilson’s success can be attributed to his strong interest and depth of knowledge in biochemistry and evolutionary biology, his insistence of quantification of evolutionary phenomena, and has early recognition of new molecular techniques that could shed light on questions of evolutionary biology. After development of quantitative immunological methods, his lab was the first to recognize restriction endonuclease mapping analysis as a quantitative evolutionary genetic method, which led to his early use of DNA sequencing, and the then-nascent technique of PCR to obtain large DNA sets for genetic analysis of populations. He trained scores of undergraduate, graduate (34 people, 17 each of men and women, received their doctoral degrees in his lab), and post-doctoral students in molecular evolutionary biology, including sabbatical visitors from six continents. His lab published more than 300 technical papers, and was recognized as a mecca for those wishing to enter the field of molecular evolution in the 1970s and 1980s.

The Allan Wilson Centre for Molecular Ecology and Evolution was established in 2002 in his honour to advance knowledge of the evolution and ecology of New Zealand and Pacific plant and animal life, and human history in the Pacific. The Centre is under the Massey University, at Palmerston North, New Zealand, and is a national collaboration involving the University of Auckland, Victoria University of Wellington, the University of Otago, University of Canterbury and the New Zealand Institute for Plant and Food Research.

A 41-minutes documentary film of his life entitled Allan Wilson, Evolutionary: Biochemist, Biologist, Giant of Molecular Biology was released by Films Media Group in 2008.

 

Allan Charles Wilson. 18 October 1934 — 21 July 1991

Rebecca L. Cann

Department of Cell and Molecular Biology, University of Hawaii at Manoa, Biomedical Sciences Building T514, 1960 East–West Rd, Honolulu, HI 96822, USA

Abstract

Allan Charles Wilson was born on 18 October 1934 at Ngaruawahia, New Zealand. He died in Seattle, Washington, on 21 July 1991 while undergoing treatment for leukemia.  Allan was known as a pioneering and highly innovative biochemist, helping to define the field of molecular evolution and establish the use of a molecular clock to measure evolutionary change between living species. The molecular clock, a method of measuring the timescale of evolutionary change between two organisms on the basis of the number of mutations that they have accumulated since last sharing a common genetic ancestor, was an idea initially championed by Émile Zuckerkandl and Linus Pauling (Zuckerkandl & Pauling 1962), on the basis of their observations that the number of changes in an amino acid sequence was roughly linear with time in the aligned hemoglobin proteins of animals. Although it is now not unusual to see the words ‘molecular evolution’ and ‘molecular phylogeny’ together, when Allan formed his own biochemistry laboratory in 1964 at the University of California, Berkeley, many scientists in the field of evolutionary biology considered these ideas complete heresy. Allan’s death at the relatively young age of 56 years left behind his wife, Leona (deceased in 2009), a daughter, Ruth (b. 1961), and a son, David (b. 1964), as well his as mother, Eunice (deceased in 2002), a younger brother, Gary Wilson, and a sister, Colleen Macmillan, along with numerous nieces, nephews and cousins in New Zealand, Australia and the USA. In this short span of time, he trained more than 55 doctoral students and helped launch the careers of numerous postdoctoral fellows.

Allan Charles Wilson, Biochemistry; Molecular Biology: Berkeley

1934-1991

Professor

The sudden death of Allan Wilson, of leukemia, on 21 July 1991, at the age of 56, and at the height of his powers, robbed the Berkeley campus and the international scientific community of one of its most active and respected leaders.

Read Full Post »

Evolution of Myoglobin and Hemoglobin

Author and Curator: Larry H. Bernstein, MD, FCAP 

Nitric oxide dioxygenase function and mechanism of flavohemoglobin, hemoglobin, myoglobin and their associated reductases

Paul R. Gardner
Journal of Inorganic Biochemistry Jan 2005;  99(1): 247–266
http://dx.doi.org:/10.1016/j.jinorgbio.2004.10.003

Microbial flavohemoglobins (flavoHbs) and hemoglobins (Hbs) show large radical dotNO dioxygenation rate constants ranging from 745 to 2900 μM−1 s−1 suggesting a primal radical dotNO dioxygenase (NOD) (EC 1.14.12.17) function for the ancient Hb superfamily. Indeed, modern O2-transporting and storing mammalian red blood cell Hb and related muscle myoglobin (Mb) show vestigial radical dotNO dioxygenation activity with rate constants of 34–89 μM−1 s−1. In support of a NOD function, microbial flavoHbs and Hbs catalyze O2-dependent cellular radical dotNO metabolism, protect cells from radical dotNO poisoning, and are induced by radical dotNO exposures. Red blood cell Hb, myocyte Mb, and flavoHb-like activities metabolize radical dotNO in the vascular lumen, muscle, and other mammalian cells, respectively, decreasing radical dotNO signalling and toxicity. HbFe(III)–OOradical dot, HbFe(III)–OONO and protein-caged [HbFe(III)–Oradical dotradical dotNO2] are proposed intermediates in a reaction mechanism that combines both O-atoms of O2 with radical dotNO to form nitrate and HbFe(III). A conserved Hb heme pocket structure facilitates the dioxygenation reaction and efficient turnover is achieved through the univalent reduction of HbFe(III) by associated reductases. High affinity flavoHb and Hb heme ligands, and other inhibitors, may find application as antibiotics and antitumor agents that enhance the toxicity of immune cell-derived radical dotNO or as vasorelaxants that increase radical dotNO signaling.

http://ars.els-cdn.com/content/image/1-s2.0-S016201340400296X-gr1.sml

NO-NOD-NOR map

http://ars.els-cdn.com/content/image/1-s2.0-S016201340400296X-gr2.sml

http://ars.els-cdn.com/content/image/1-s2.0-S016201340400296X-gr3.sml

The evolution of the globin family genes: Concordance of stochastic and augmented maximum parsimony genetic distances for α hemoglobin, β hemoglobin and myoglobin phylogenies
R Holmquist, TH Jukes, H Moise, M Goodman, GW Moore
Journal of Molecular Biology Jul 1976; 105(1): 39–74
http://dx.doi.org:/10.1016/0022-2836(76)90194-7

We compare the amino acid sequences of 70 globing, representing the following families: (a) α hemoglobin chains; (b) β hemoglobin chains; (c) myoglobins; (d) two lamprey, a mollusc, and two plant globins. The comparisons show a convergence of maximal and minimal estimates of genetic differences as calculated respectively by the stochastic and maximum parsimony procedures, thus demonstrating for the first time the logical consistency and complementarity of the two procedures. Evolutionary rates are non-constant, varying over a range of 1 to 75 nucleotide replacements per 100 codons per 108 years. These rate differentials are resolved into two components (a) due to change in the number of codon sites free to fix mutations during the period of divergence of the species involved; (b) due to change in fixation intensity at each site. These two components also show non-uniformity along different lineages. Positive Darwinian natural selection can bring about an increase in either component, and negative or stabilizing selection in protein evolution can lead to decreases. Accelerated rates of globin evolution were found in lineages of cold-blooded vertebrates, some marsupials, and early placental mammals, while slower rates were found in warm-blooded vertebrates, especially higher primates. One manifestation of negative selection in the globins is that minimal 3-base type amino acid replacements occur less frequently than would be expected if base replacements had occurred and were accepted at random. The selection against these replacements is not due to atypical behavior with respect to the change in electrical charge involved in the replacements. Interestingly, the globins from the lamprey, sea hare and the legumes are as distant from one another as are α-hemoglobin and β-hemoglobin from myoglobin.

Hemoglobin Orthologs
http://www.bio.davidson.edu/Courses/Molbio/MolStudents/spring2005/Heiner/ortholog.html

Orthologs are sequences of genes that evolved from a common ancestor and can be traced evolutionarily through different species. By comparing the ortholog sequences of a specific gene between many species, the amino acid sequences which are conserved can be determined. These highly conserved sequences are important, because they provide information on which amino acids are essential to the protein structure and function.

Evolution of Hemoglobin

Hemoglobin is derived from the myoglobin protein, and ancestral species just had myoglobin for oxygen transport. 500 million years ago the myoglobin gene duplicated and part of the gene became hemoglobin. Lampreys are the most ancestral animal to have hemoglobin, and the ancestral version was composed of dimers instead of tetramers and was only weakly cooperative. 100 million years later, the hemoglobin gene duplicated again forming alpha and beta subunits. This form of derived hemoglobin is found in bony fish, reptiles, and mammals, which all have both alpha and beta subunits to form a tetramer (Mathews et al., 2000).

Conserved Sequences

When the amino acid sequences of myoglobin, the hemoglobin alpha subunit, and the hemoglobin beta subunit are compared, there are several amino acids that remain conserved between all three globins (Mathews et al., 2000). These amino acid sequences are considered truly essential, because they have remained unchanged throughout evolution, and therefore are fundamental to the function of the protein. These essential amino acids can be seen in Figure 1, which compares myoglobin, and the alpha and beta subunits of hemoglobin. The histidines in helix F8 and helix E7 are highly conserved. These histidines are located proximally and distally to the heme molecule and keep the heme molecule in place within the hemoglobin protein as seen in Figure 2 (Mathews et al., 2000). This shows that the position of the heme molecule within the globin protein is essential to its function. Likewise, the amino acids in the FG region are also highly conserved. This region of the protein is essential to the conformational change between the T to R states (Mathews et al., 2000). Additionally, the amino acids at the alpha-beta subunit interfaces are highly conserved, because they also affect the conformational change between the subunits, which regulates oxygen affinity and cooperativity. In general, the most highly conserved sequences are located within the interior of the hemoglobin protein where the subunits contact each other (Gribaldo et al., 2003).

A cartoon drawing of the structure of hemoglobin around heme molecule. The histadines in helix F8 and E7 interact directly with the heme molecule.  figure2

http://www.bio.davidson.edu/Courses/Molbio/MolStudents/spring2005/Heiner/figure2.jpg

Figure 2: A cartoon drawing of the structure of hemoglobin around heme molecule. The histadines in helix F8 and E7 interact directly with the heme molecule. http://www.aw-bc.com/mathews/ch07/fi7p5.htm  (permission pending).

Figure 1: The amino acid sequences of myoglobin, alpha subunit of hemoglobin, and beta subunit of hemoglobin. The amino acid sequences highlighted in tan are conserved between all three globins and the amino acid sequences highlighted in gray are conserved between alpha and beta hemoglobin. http://www.aw-bc.com/mathews/ch07/fi7p11.htm (permission pending).

http://www.bio.davidson.edu/Courses/Molbio/MolStudents/spring2005/Heiner/figure1.jpg

Figure 2: A cartoon drawing of the structure of hemoglobin around heme molecule. The histidines in helix F8 and E7 interact directly with the heme molecule. http://www.aw-bc.com/mathews/ch07/fi7p11.htm (permission pending).

Alpha Subunit of Hemoglobin

The alpha subunit of hemoglobin has several amino acid sequences that are conserved across many species and are essential to its function. The alpha subunit of hemoglobin is encoded by the 2 genes HBA1 and HBA2 both located on chromosome 16 (GeneCard, 2005). Click here to see the gene card for HBA1. To determine which amino acid sequences are conserved, I compared the orthologs of HBA1 in Homo sapiens (humans) to 5 additional species including, Xenopus tropicalis (African clawed frog), Danio rerio (Zebra fish), Gallus gallus (Red jungle fowl), Mus musculus (mouse), and Rattus norvegicus (rat) using the Ensembl program. Figure 3 shows the 6 orthologs aligned and the important conserved regions highlighted. The stars indicate amino acids that are conserved between all of the species. As a general observation, the mouse ortholog of HBA is the most similar to human HBA, because it is the most evolutionarily related. The amino acid sequences that are conserved in all globin proteins (highlighted in blue) can be seen in Figure 3. There are also several conserved amino acids that are specifically important to HBA structure (highlighted in red) including: the phenylalanine (F) at position 44, which is in direct contact with the heme group; tyrosine (Y) at position 142, which stabilizes the hemoglobin molecule by forming hydrogen bonds between two of the helices; and glycine (G) at position 26, which is small and therefore allows two of the helices to approach each other, which is important to the structure of hemoglobin (Natzke, 1998). Additionally, there are several proteins found in the alpha subunit that are involved in the movement of the alpha and beta subunits (also highlighted in red) including: the tyrosine (Y) at position 43, which interacts with the beta subunit during the R state, and the arginine (N) at position 143, which interacts with the beta subunit during the T state (Gribaldo et al., 2003).

Mutations

Looking at the effects mutated portions of a gene is also a good way to determine the function of highly conserved sequences. In hemoglobin, deleterious mutations are most common in the heme pockets of the protein and in the alpha and beta subunit interfaces (Mathews et al., 2000). There are several key mutations in highly conserved portions of HBA (highlighted in yellow) including: the substitution of histidine (H) at position 88 to tyrosine (Y), which disrupts the heme molecule leading to decreased oxygen affinity; the substitution of arginine (N) at position 143 to histidine (H), which eliminates a bond in the T state and therefore favors the R state, resulting in increased oxygen affinity; the substitution of proline (P) at position 97 to arginine (N), which alters the alpha-beta contact region and results in the disassociation of the hemoglobin complex; and the substitution of leucine (L) at position 138 for proline (P), which interrupts the helix formation and also results in the disassociation of the hemoglobin complex (Mathews et al., 2000).

Bar-headed Goose Hemoglobin

As mentioned on the previous page, the bar-headed goose has hemoglobin that is specifically adapted to high altitudes. The bar-headed goose hemoglobin has an increased oxygen affinity which allows it to live in low oxygen pressure environments (Liang et al., 2001). This increased oxygen affinity is the result of a mutation at position 121 in the alpha subunit, which is highly conserved in other species, from proline to alanine, as seen in Figure 4 (Liang et al., 2001). This substitution leaves a two-carbon gap between the alpha-beta dimer, which relaxes the T structure and allows it to bind oxygen more readily under lower pressures (Jessen et al. 1991). Thus, comparing orthologs can also be used to explain differences in the oxygen binding capabilities of hemoglobin in different species.

References

Ensembl. Ensembl Genome Browser. http://www.ensembl.org/. Accessed March 2005.

GeneCard. 2005. GeneCard for HBA1. http://genome-www.stanford.edu/cgi-bin/genecards/carddisp?HBA1&search=HBA&suff=txt. Accessed March 2005.

Gribaldo, Simonetta, Didier Casane, Philippe Lopez and Herve Philippe. 2003. Functional Divergence Prediction from Evolutionary Analysis: A Case Study of Vertebrate Hemoglobin. Molecular Biology and Evolution 20 (11): 1754-1759.

Jessen, Timm H et al. 1991. Adaptation of bird hemoglobins to high altitudes: Demonstration of molecular mechanism by protein engineering. Evolution 88: 6519-6522.

Liang, Yuhe et al. 2001. The Crystal Structure of Bar-headed Goose Hemoglobin in Deoxy Form: The Alloseteric Mechanism of a Hemoglobin Species with High Oxygen Affinity. Journal of Molecular Biology 313: 123-137.

Mathews, Christopher, Kensal Van Holde and Kevin Ahern. 2000. Biochemistry 3 rd edition. http://www.aw-bc.com/mathews/ch07/c07emhp.htm .   Accessed March 2005.

Natzke, Lisa. 1998. Hemoglobin. http://biology.kenyon.edu/BMB/Chime/Lisa/FRAMES/hemetext.htm. Accessed March 2005.

Divergence pattern and selective mode in protein evolution: the example of vertebrate myoglobins and hemoglobin chains.
Otsuka J1, Miyazaki K, Horimoto K.
J Mol Evol. 1993 Feb; 36(2):153-81.

The evolutionary relation of vertebrate myoglobin and the hemoglobin chains including the agnathan hemoglobin chain is investigated on the basis of a new view of amino acid changes that is developed by canonical discriminant analysis of amino acid residues at individual sites. In contrast to the clear discrimination of amino acid residues between myoglobin, hemoglobin alpha chain, and hemoglobin beta chain in warm-blood vertebrates, the three types of globins in the lower class of vertebrates show so much variation that they are not well discriminated. This is seen particularly at the sites that are ascertained in mammals to carry the amino acid residues participating in stabilizing the monomeric structure in myoglobin and the residues forming the subunit contacts in hemoglobin. At these sites, agnathan hemoglobin chains are evaluated to be intermediate between the myoglobin and hemoglobin chains of gnathostomes. The variation in the phylogenetically lower class of globins is also seen in the internal region; there the amino acid residues of myoglobin and hemoglobin chains in the phylogenetically higher class exhibit an example of parallel evolution at the molecular level. New quantities, the distance of sequence property between discriminated groups and the variation within each group, are derived from the values of discriminant functions along the peptide chain, and this set of quantities simply describes an overall feature of globins such that the distinction between the three types of globins has been clearer as the vertebrates have evolved to become jawed, landed, and warm-blooded. This result strongly suggests that the functional constraint on the amino acid sequence of a protein is changed by living conditions and that severe conditions constitute a driving force that creates a distinctive protein from a less-constrained protein.

The globin gene repertoire of lampreys: Convergent evolution of hemoglobin and myoglobin in jawed and jawless vertebrates
K Schwarze, KL Campbell, T Hankeln, JF Storz, FG Hoffmann and T Burmester
Mol Biol Evol (2014).  http://dx.doi.org:/10.1093/molbev/msu216

Agnathans (jawless vertebrates) occupy a key phylogenetic position for illuminating the evolution of vertebrate anatomy and physiology. Evaluation of the agnathan globin gene repertoire can thus aid efforts to reconstruct the origin and evolution of the globin genes of vertebrates, a superfamily that includes the well-known model proteins hemoglobin and myoglobin. Here we report a comprehensive analysis of the genome of the sea lamprey (Petromyzon marinus) which revealed 23 intact globin genes and two hemoglobin pseudogenes. Analyses of the genome of the Arctic lamprey (Lethenteron camtschaticum) identified 18 full length and five partial globin gene sequences. The majority of the globin genes in both lamprey species correspond to the known agnathan hemoglobins. Both genomes harbor two copies of globin X, an ancient globin gene that has a broad phylogenetic distribution in the animal kingdom. Surprisingly, we found no evidence for an ortholog of neuroglobin in the lamprey genomes. Expression and phylogenetic analyses identified an ortholog of cytoglobin in the lampreys; in fact, our results indicate that cytoglobin is the only orthologous vertebrate-specific globin that has been retained in both gnathostomes and agnathans. Notably, we also found two globins that are highly expressed in the heart of P. marinus, thus representing functional myoglobins. Both genes have orthologs in L. camtschaticum. Phylogenetic analyses indicate that these heart-expressed globins are not orthologous to the myoglobins of jawed vertebrates (Gnathostomata), but originated independently within the agnathans. The agnathan myoglobin and hemoglobin proteins form a monophyletic group to the exclusion of functionally analogous myoglobins and hemoglobins of gnathostomes, indicating that specialized respiratory proteins for O2 transport in the blood and O2 storage in the striated muscles evolved independently in both lineages. This dual convergence of O2-transport and O2-storage proteins in agnathans and gnathostomes involved the convergent co-option of different precursor proteins in the ancestral globin repertoire of vertebrates.

Globin evolution
Kent Holsinger
http://darwin.eeb.uconn.edu/eeb348/lecturenotes/molevol-multigene/node2.html

I’ve just pointed out the distinction between myoglobin and hemoglobin. You may also remember that hemoglobin is a multimeric protein consisting of four subunits, 2 α\alpha subunits and 2 β\beta subunits. What you may not know is that in humans there are actually two types of α\alpha hemoglobin and four types of β\beta hemoglobin, each coded by a different genetic locus (see Table 1). The five α\alpha -globin loci (α\alpha_1, α\alpha_2, ς\zeta, and two non-functional pseudogenes) are found in a cluster on chromosome 16. The six β\beta-globin loci (ε\epsilon, ϒ\gamma_G, ϒ\gamma_A, δ\delta, β\beta, and a pseudogene) are found in a cluster on chromosome 11. The myoglobin locus is on chromosome 22.

Table 1: Human hemoglobins arranged in developmental sequence. Adult hemoglobins composed of 2 and 2 subunits typically account for less than 3% of hemoglobins in adults (http://sickle.bwh.harvard.edu/hbsynthesis.html).

Not only do we have all of these different types of globin genes in our bodies, they’re all related to one another. Comparative sequence analysis has shown that vertebrate myoglobin and hemoglobins diverged from one another about 450 million years ago. Figure 1 shows a phylogenetic analysis of globin genes from humans, mice, and a variety of Archaea. Focus your attention on the part of the tree that has human and mouse sequences. You’ll notice two interesting things:

Human and mouse neuroglobins (Ngb) are more closely related to one another than they are to other globins, even those from the same species. The same holds true for cytoglobins (Cyg) and myoglobins (Mb).

Within the hemoglobins, only mouse β\beta-globin (Mouse HbB) is misplaced. All other α\alpha- and β\beta-globins group with the corresponding mouse and human loci.

This pattern is exactly what we expect as a result of duplication and divergence. Up to the time that a gene becomes duplicated, its evolutionary history matches the evolutionary history of the organisms containing it. Once there are duplicate copies, each follows an independent evolutionary history. Each traces the history of speciation and divergence. And over long periods duplicate copies of the same gene share more recent common ancestry with copies of the same gene in a different species than they do with duplicate genes in the same genome.

Figure 1: Evolution of globin genes in Archaea and mammals (from [2]).

http://darwin.eeb.uconn.edu/eeb348/lecturenotes/molevol-multigene/img11.png

Evolution of globin genes in Archaea and mammals

Evolution of globin genes in Archaea and mammals

A history of duplication and divergence in multigene families makes it important to distinguish between two classes of related loci: those that represent the same locus in different species and between which divergence is a result of species divergence are orthologs. Those that represent different loci and between which divergence occurred after duplication of an ancestral gene are paralogs. The β\beta-globin loci of humans and chickens are orthologous. The α\alpha $- and $\beta $-globin loci of any pair of taxa are paralogous.

As multigene families go, the globin family is relatively simple and easy to understand. There are only about a dozen loci involved, one isolated locus (myoglobin) and two clusters of loci ($\alpha- and β\beta-globins). You’ll find a diagram of the β\beta-globin cluster in Figure 2. As you can see the β\beta-globins are not only evolutionarily related to one another they occur relatively close to one another on chromosome 11 in humans.

Figure 2: Structure of the human β\beta-globin gene cluster. % identity refers to similarity to the mouse β\beta-globin sequence. From http://globin.cse.psu.edu/html/pip/betaglobin/iplot.ps  (retrieved 28 Nov 2006).

Other families are far more complex. Class I and class II MHC loci, for example are part of the same multigene family. Moreover, immunoglobulins, T-cell receptors, and, and MHC loci are part of a larger superfamily of genes, i.e., all are ultimately derived from a common ancestral gene by duplication and divergence. Table 2 lists a few examples of multigene families and superfamilies in the human genome and the number of proteins produced.

Table 2: A few gene families from the human genome (adapted from [5,6]).
Distribution and conservation of sequence

Distribution and conservation of sequence

https://encrypted-tbn3.gstatic.com/images?q=tbn:ANd9GcRHcfUpQ09ufj8cleSgnhDfQVUEHsTvnYGNxKaPa5wxMqNFzFU6

Distribution and conservation of sequence motifs throughout mammalian beta-globin gene clusters.A detailed map of the gene cluster is shown on the numbered line

evolutionary history of three hypothetical living species (C, D, and E)

evolutionary history of three hypothetical living species (C, D, and E)

the evolutionary history of three hypothetical living species (C, D, and E), inferred by comparing amino-acid differences in their myoglobin molecules.

http://media-1.web.britannica.com/eb-media/98/52998-004-A8682A5B.jpg

oxyhemoglobin dissociation curve

oxyhemoglobin dissociation curve

much higher affinity for oxygen than haemoglobin.

much higher affinity for oxygen than haemoglobin.

http://i.stack.imgur.com/WQJe9.jpg

myoglobin hs much higher oxygen affinity than Hb

Evolution of Myoglobin / Hemoglobin Proteins

Primitive Globin – Very primitive animals had only a myoglobin-like, single-chain ancestral globin for oxygen storage and were so small that they did not require a transport protein. Roughly 500 million years ago the ancestral myoglobin gene was duplicated. One copy became the ancestor of the myoglobin genes of all higher organisms. The other copy evolved into the gene for an oxygen transport protein and gave rise to the hemoglobins.

Most Primitive Hemoglobin – The most primitive animals to possess hemoglobin are the lampreys. Lamprey hemoglobin can form dimers but not tetramers and is only weakly cooperative. It represents a first step toward allosteric binding. Subsequently a second gene duplication must have occurred, giving rise to the ancestors of the present-day  and  hemoglobin chain families. This must have happened about 400 million years ago, at about the time of divergence of the sharks and bony fish. The evolutionary line of the bony fish led to the reptiles and eventually to the mammals, all carrying genes for both  and  globins and capable of forming tetrameric 22 hemoglobins. Further gene duplications have occurred in the hemoglobin line, leading to the embryonic forms  and , the fetal form, , and the infant form  (Figure 7.22).

Conserved Amino Acid Sequences – During the long evolution of the myoglobin/hemoglobin family of proteins, only a few amino acid residues have remained invariant (Figure 7.11). They include the histidines proximal and distal to the heme iron (F8 and E7- see Figure 7.5b) and Val FG5, which has been implicated in the hemoglobin deoxy/oxy conformation change. These may mark the truly essential positions in the molecule. Other regions highly conserved in hemoglobins are those near the 1 – 2 and 2 – 1 contacts. These parts of the molecule are most directly involved in the allosteric conformational change.

 http://web.squ.edu.om/med-lib/med_cd/e_cds/Electronic%20Study%20Guide%20of%20Biochemistry/ch07/c07emhp.htm

Read Full Post »

More Complexity in Protein Evolution

Author and Curator: Larry H. Bernstein, MD, FCAP 

Lactate dehydrogenase like crystallin: a potentially protective shield for Indian spiny-tailed lizard (Uromastix ltardwickit) lens against environmental stress?
A Atta, A Ilyas, Z Hashim, A Ahmed and S Zarina
The Protein Journal 2014; 33(2), p. 128-34.
http://dx.doi.org/10.1007/s10930-014-9543-4

Taxon specific lens crystallins in ve1iebrates are either similar or identical with various metabolic enzymes. These bifunctional crystallins serve as structural protein in lens along with their catalytic role. In the present study, we have partially purified and characterized lens crystallin from Indian spiny-tailed lizard (Uroma stix hardwickii). We have found lactate dehydrogenase (LDH) activity in lens indicating presence of an enzyme crystallin with dual functions. Taxon specific lens crystallins are product of gene sharing or gene duplication phenomenon where a pre-existing enzyme is recruited as lens crystallin in addition to structural role. In lens, same gene adopts refractive role in lens without modification or loss of pre-existing function during gene sharing phenomenon. Apart from conventional role of structural protein, LDH activity containing crystallin in Uromastix hardwickii lens is likely to have adaptive characteristics to offer protection against toxic effects of oxidative stress and ultraviolet light, hence justifying its recruitment. Taxon specific crystallins may serve as good models to understand structure-function relationship of these proteins.

αB-Crystallin and 27-kd Heat Shock Protein Are Regulated by Stress Conditions in the Central Nervous System and Accumulate in Rosenthal Fibers
T Iwaki, A Iwaki, J Tateishi, Y Sakaki, and JE Goldmant
Ameri J Pathol  1993; 143(2):487-495.

To understand the significance of the accumulation of αB-crystallin in Rosenthal fibers within astrocytes, the expression and metabolism of αB-crystallin in glioma cell lines were examined under the conditions of heat and oxidative stress. αB-crystallin mRNA was increased after both stresses, and αB-crystallin protein moved from a detergent-soluble to a detergent-insoluble form. In addition, Western blotting of Alexander’s  disease brain homogenates revealed that the 27-kd heat shock protein (HSP27), which is related to αB-crystallin, accumulates along with αB-crystallin. The presence of HSP27 in Rosenthal fibers was directly demonstrated by immunohistochemistry. Our results suggest that astrocytes in Alexander’s disease may be involved in an as yet unknown kind of stress reaction that causes the accumulation of αB-ccystallin and HSP27 and results in Rosenthal fiber formation.

α-Crystallin can function as a molecular chaperone
Joseph Horwitz
Proc. Nadl. Acad. Sci. USA Nov 1992; 89: 10449-10453. Biochemistry

The α-crystallins (αA and αB) are major lens structural proteins of the vertebrate eye that are related to the small heat shock protein family. In addition, crystallins (especially αB) are found in many cells organs outside the lens, and aα is overexpressed in several neurological disorders and in cell lines under stress conditions. Here I show that α-crystallin can function as a molecular chaperone. Stoichiometric amounts of αA and αB suppress thermally induced aggregation of various enzymes. In particular, α-crystalln is very efficient in suppressing the thermally induced aggregation of β- and y-crystallins, the two other major mammalian stuctural lens proteins. α-Crystallin was also effective in preventing aggregation and in refolding guanidine hydrochloride-denatured y-crystallin, as judged by circular dichroism spectroscopy. My results thus indicate that α-crystallin refracts light and protects proteins from aggregation in the transparent eye lens and that in nonlens cells α-crystallin may have other functions in addition to its capacity to suppress aggregation of proteins.

Gene sharing by δ-crystallin and argininosuccinate Iyase
J Piatigorsky, WE O’Brient, BL Norman, K Kalumuckt, GJ Wistow, T Borras, et al.
Proc. Natl. Acad. Sci. USA  May 1988; 85: 3479-3483. Evolution.

The lens structural protein δ-crystallin and the metabolic enzyme argininosuccinate lyase (ASL; Largininosuccinate argine-lyase, EC 4.3.2.1) have striking sequence similarity. We have demonstrated that duck δ-crystallin has enormously high ASL activity, while chicken δ-crystallin has lower but significant activity. The lenses of these birds had much greater ASL activity than other tissues, suggesting that ASL is being expressed at unusually high levels as a structural component. In Southern blots of human genomic DNA, chicken δ1-crystallin cDNA hybridized only to the human ASL gene; moreover, the two chicken δ-crystallin genes accounted for all the sequences in the chicken genome able to cross-hybridize with a human ASL cDNA, with preferential hybridization to the δ2 gene. Correlations of enzymatic activity and recent data on mRNA levels in the chicken lens suggest that ASL activity depends on expression of the δ2-crystallin gene. The data indicate that the same gene, at least in ducks, encodes two different functions, an enzyme (ASL) and a structural protein (δ-crystallin), although in chickens specialization and separation of functions may have occurred.

Gecko i-crystallin: How cellular retinol-binding protein became an eye lens ultraviolet filter
PJ L Werten, Beate Roll, DMF van Aalten, and WW de Jong
PNAS Mar 2000; 97(7): 3282–3287 http://pnas.org/cgi/doi/10.1073ypnas.050500597

Eye lenses of various diurnal geckos contain up to 12% i-crystallin. This protein is related to cellular retinol-binding protein type I (CRBP I) but has 3,4-didehydroretinol, rather than retinol, as a ligand. The 3,4-didehydroretinol gives the lens a yellow color, thus protecting the retina by absorbing short-wave radiation. i-Crystallin could be either the gecko’s housekeeping CRBP I, recruited for an additional function in the lens, or the specialized product of a duplicated CRBP I gene. The finding of the same CRBP I-like sequence in lens and liver cDNA of the gecko Lygodactylus picturatus now supports the former option. Comparison with i-crystallin of a distantly related gecko, Gonatodes vittatus, and with mammalian CRBP I, suggests that acquiring the additional lens function is associated with increased amino acid changes. Compared with the rat CRBP I structure, the i-crystallin model shows reduced negative surface charge, which might facilitate the required tight protein packing in the lens. Other changes may provide increased stability, advantageous for a long-living lens protein, without frustrating its role as retinol transporter outside the lens. Despite a number of replacements in the ligand pocket, recombinant i-crystallin binds 3,4-didehydroretinol and retinol with similar and high affinity (1.6 nM). Availability of ligand thus determines whether it binds 3,4-didehydroretinol, as in the lens, or retinol, in other tissues. i-Crystallin presents a striking example of exploiting the potential of an existing gene without prior duplication.

Expression of βA3/A1-crystallin in the developing and adult rat eye
G Parthasarathy, Bo Ma, C Zhang, C Gongora, JS Zigler, MK Duncan, D Sinha
J Molec Histol 2011; 42(1): 59-69. http://dx.doi.org:/10.1007/s10735-010-9307-1

Crystallins are very abundant structural proteins of the lens and are also expressed in other tissues. We have previously reported a spontaneous mutation in the rat βA3/A1-crystallin gene, termed Nuc1, which has a novel, complex, ocular phenotype. The current study was undertaken to compare the expression pattern of this gene during eye development in wild type and Nuc1 rats by in situ hybridization (ISH) and immunohistochemistry (IHC).
βA3/A1-crystallin expression was first detected in the eyes of both wild type and Nuc1 rats at embryonic (E) day 12.5 in the posterior portion of the lens vesicle, and remained limited to the lens fibers throughout fetal life.
After birth, βA3/A1-crystallin expression was also detected in the neural retina (specifically in the astrocytes and ganglion cells) and in the retinal pigmented epithelium (RPE).
This suggested that βA3/A1-crystallin is not only a structural protein of the lens, but has cellular function(s) in other ocular tissues.
In summary, expression of βA3/A1-crystallin is controlled differentially in various eye tissues with lens being the site of greatest expression.
Similar staining patterns, detected by ISH and IHC, in wild type and Nuc1 animals suggest that functional differences in the protein, rather than changes in mRNA/protein level of expression likely account for developmental abnormalities in Nuc1.

βA3/A1Crystallin controls anoikis-mediated cell death in astrocytes by modulating PI3K/AKT/mTOR and ERK survival pathways through the PKD/Bit1-signaling axis
B Ma, T Sen, L Asnaghi, M Valapala, F Yang, S Hose, D S McLeod, Y Lu, et la.
Cell Death and Disease 2011; 2(10). http://dx.doi.org:/10.1038/cddis.2011.100

During eye development, apoptosis is vital to the maturation of highly specialized structures such as the lens and retina. Several forms of apoptosis have been described, including anoikis, a form of apoptosis triggered by inadequate or inappropriate cell–matrix contacts. The anoikis regulators, Bit1 (Bcl-2 inhibitor of transcription-1) and protein kinase-D (PKD), are expressed in developing lens when the organelles are present in lens fibers, but are downregulated as active denucleation is initiated.
We have previously shown that in rats with a spontaneous mutation in the Cryba1 gene, coding for βA3/A1-crystallin, normal denucleation of lens fibers is inhibited. In rats with this mutation (Nuc1), both Bit1 and PKD remain abnormally high in lens fiber cells. To determine whether βA3/A1-crystallin has a role in anoikis, we induced anoikis in vitro and conducted mechanistic studies on astrocytes, cells known to express βA3/A1-crystallin.
The expression pattern of Bit1 in retina correlates temporally with the development of astrocytes. Our data also indicate that loss of βA3/A1-crystallin in astrocytes results in a failure of Bit1 to be trafficked to the Golgi, thereby suppressing anoikis. This loss of βA3/A1-crystallin also induces insulin-like growth factor-II, which increases cell survival and growth by modulating the phosphatidylinositol-3-kinase (PI3K)/AKT/mTOR and extracellular signal-regulated kinase pathways. We propose that βA3/A1-crystallin is a novel regulator of both life and death decisions in ocular astrocytes.

βA3/A1-crystallin in astroglial cells regulates retinal vascular remodeling during development
D Sinha, A Klise, Y Sergeev, S Hose, IA Bhutto, L Hackler Jr., T Malpic-llanos, et al.
Molec Cell Neurosci 2008; 37(1): 85-95.

http://dx.doi.org:/10.1016/j.mcn.2007.08.016

Vascular remodeling is a complex process critical to development of the mature vascular system. Astrocytes are known to be indispensable for initial formation of the retinal vasculature; our studies with the Nuc1 rat provide novel evidence that these cells are also essential in the retinal vascular remodeling process.
Nuc1 is a spontaneous mutation in the Sprague–Dawley rat originally characterized by nuclear cataracts in the heterozygote and microphthalmia in the homozygote. We report here that the Nuc1 allele results from mutation of the βA3/A1-crystallin gene, which in the neural retina is expressed only in astrocytes. We demonstrate striking structural abnormalities in Nuc1 astrocytes with profound effects on the organization of intermediate filaments. While vessels form in the Nuc1 retina, the subsequent remodeling process required to provide a mature vascular network is deficient. Our data implicate βA3/A1-crystallin as an important regulatory factor mediating vascular patterning and remodeling in the retina.

A developmental defect in astrocytes inhibits programmed regression of the hyaloid vasculature in the mammalian eye
C Zhang, L Asnaghi, C Gongora, B Patek, S Hose, Bo Ma, MA Fard, L Brako, et al.
Eur J Cell Biol 2011; 90(5): 440-448.
http://dx.doi.org:/10.1016/j.ejcb.2011.01.003

Previously we reported the novel observation that astrocytes ensheath the persistent hyaloid artery, both in the Nuc1 spontaneous mutant rat, and in human PFV (persistent fetal vasculature) disease (Developmental Dynamics 234:36–47, 2005). We now show that astrocytes isolated from both the optic nerve and retina of Nuc1 rats migrate faster than wild type astrocytes. Aquaporin 4 (AQP4), the major water channel in astrocytes, has been shown to be important in astrocyte migration. We demonstrate that AQP4 expression is elevated in the astrocytes in PFV conditions, and we hypothesize that this causes the cells to migrate abnormally into the vitreous where they ensheath the hyaloid artery. This abnormal association of astrocytes with the hyaloid artery may impede the normal macrophage-mediated remodeling and regression of the hyaloid system.

βA3/A1-crystallin is required for proper astrocyte template formation and vascular remodeling in the retina.
D Sinha; WJ Stark; M Valapala; IA Bhutto; M Cano; S Hose; GA Lutty; et al.  Transgenic research 2012; 21(5):1033-42.

Nuc1 is a spontaneous rat mutant resulting from a mutation in the Cryba1 gene, coding for βA3/A1-crystallin. Our earlier studies with Nuc1 provided novel evidence that astrocytes, which express βA3/A1-crystallin, have a pivotal role in retinal remodeling. The role of astrocytes in the retina is only beginning to be explored. One of the limitations in the field is the lack of appropriate animal models to better investigate the function of astrocytes in retinal health and disease. We have now established transgenic mice that overexpress the Nuc1 mutant form of Cryba1, specifically in astrocytes. Astrocytes in wild type mice show normal compact stellate structure, producing a honeycomb-like network. In contrast, in transgenics over-expressing the mutant (Nuc1) Cryba1 in astrocytes, bundle-like structures with abnormal patterns and morphology were observed. In the nerve fiber layer of the transgenic mice, an additional layer of astrocytes adjacent to the vitreous is evident. This abnormal organization of astrocytes affects both the superficial and deep retinal vascular density and remodeling. Fluorescein angiography showed increased venous dilation and tortuosity of branches in the transgenic retina, as compared to wild type. Moreover, there appear to be fewer interactions between astrocytes and endothelial cells in the transgenic retina than in normal mouse retina. Further, astrocytes overexpressing the mutant βA3/A1-crystallin migrate into the vitreous, and ensheath the hyaloid artery, in a manner similar to that seen in the Nuc1 rat. Together, these data demonstrate that developmental abnormalities of astrocytes can affect the normal remodeling process of both fetal and retinal vessels of the eye and that βA3/A1-crystallin is essential for normal astrocyte function in the retina.

Ontogeny of oxytocin and vasopressin receptor binding in the lateral septum in prairie and montane voles
Z. Wang, L.J. Young
Developmental Brain Research 1997; 104:191–195.

Adult prairie (Microtus ochrogaster). and montane voles (M. montanus). differ in the distribution of oxytocin OT. and vasopressin AVP receptor binding in the brain. The present study examined the ontogenetic pattern of these receptor bindings in the lateral septum in both species to determine whether adult differences in the receptor binding are derived from a common pattern in development. In both species, OT and AVP receptor binding in the lateral septum were detected neonatally, increased during development, and reached the adult level at weaning third week. The progression of OT and AVP receptor differed, as OT receptor binding increased continually until weaning while AVP receptor binding did not change in the first week, increased rapidly in the second week, and was sustained thereafter. For both receptors, the binding increased more rapidly in montane than in prairie voles, resulting in species differences in receptor binding at weaning and in adulthood. Together, these data indicate that OT and AVP could affect the brain during development in a peptide- and species-specific manner in voles.

Evolution of the vasopressin/oxytocin superfamily: Characterization of a cDNA encoding a vasopressin-related precursor, preproconopressin, from the mollusc Lymnaea stagnalis
RE Van Kesteren, AB Smit, RW Dirksi, ND De With, WPM Geraerts, and J Joosse
Proc. Nadl. Acad. Sci. USA May 1992; 89: 4593-4597. Neurobiology

Although the nonapeptide hormones vasopressin, oxytocin, and related peptides from vertebrates and some nonapeptides from invertebrates share similarities in amino acid sequence, their evolutionary relationships are not dear. To investigate this issue, we doned a cDNA encoding a vasopressin-related peptide, Lys-conopressin, produced in the central nervous system of the gastropod mollusc Lymnaea stagnalis. The predicted preproconopressin has the overall architecture of vertebrate preprovasopressins, with a signal peptide, Lys-conopressin, that is flanked at the C terminus by an amidation signal and a pair of basic residues, followed by a neurophysin domain. The Lymnaea neurophysin and the vertebrate neurophysins share high sequence identity, which includes the conservation of all 14 cysteine residues. In addition, the Lymnaea neurophysin possesses unique structural characteristics. It contains a putative N-linked glycosylation site at a position in the vertebrate neurophysins where a strictly conserved tyrosine residue, which plays an essential role in binding of the nonapptide hormones, is found. The C-terminal copeptin homologous extension of the Lymnaea neurophysin has low sequence identity with the vertebrate counterparts and is probably not cleaved from the prohormone, as are the mammalin copeptins. The conopressin gene is expressed in only a few neurons in both pedal ganglia of the central nervous system. The conopressin transcript is present in two sizes, due to alternative use of polyadenylylation signals. The data presented here demonstrate that the typical organization of the prohormones of the vasopressin/oxytocin superfamily must have been present in the common ancestors of vertebrates and invertebrates.

A common allele in the oxytocin receptor gene (OXTR) impacts prosocial temperament and human hypothalamic-limbic structure and function
H Tosta, B Kolachanaa, S Hakimia, H Lemaitrea, BA Verchinskia, et al.
PNAS Aug 3, 2010; 107(31): 13936–13941
http://pnas.org/cgi/doi/10.1073/pnas.1003296107

The evolutionarily highly conserved neuropeptide oxytocin is a key mediator of social and emotional behavior in mammals, including humans. A common variant (rs53576) in the oxytocin receptor gene (OXTR) has been implicated in social-behavioral phenotypes, such as maternal sensitivity and empathy, and with neuropsychiatric disorders associated with social impairment, but the intermediate neural mechanisms are unknown. Here, we used multimodal neuroimaging in a large sample of healthy human subjects to identify structural and functional alterations in OXTR risk allele carriers and their link to temperament. Activation and interregional coupling of the amygdala during the processing of emotionally salient social cues was significantly affected by genotype. In addition, evidence for structural alterations in key oxytocinergic regions emerged, particularly in the hypothalamus. These neural characteristics predicted lower levels of reward dependence, specifically in male risk allele carriers. Our findings identify sex-dependent mechanisms impacting the structure and function of hypothalamic-limbic circuits that are of potential clinical and translational significance.
Test of Association Between 10 SNPs in the Oxytocin Receptor Gene and Conduct Disorder
JT Sakai, TJ Crowley, MC Stallings, M McQueen, JK Hewitt, C Hopfer, et al.
Psychiatr Genet. 2012 Apr; 22(2): 99–102. http://dx.doi.org:/10.1097/YPG.0b013e32834c0cb2

Animal and human studies have implicated oxytocin (OXT) in affiliative and prosocial behaviors. We tested whether genetic variation in the OXT receptor (OXTR) gene is associated with conduct disorder (CD).
Utilizing a family-based sample of adolescent probands recruited from an adolescent substance abuse treatment program, control probands and their families (total sample n=1,750), we conducted three tests of association with CD and 10 SNPs (single nucleotide polymorphisms) in the OXTR gene: (1) family-based comparison utilizing the entire sample; (2) within-Whites, case control comparison of adolescent patients with CD and controls without CD; and (3) within-Whites case-control comparison of parents of patients and parents of controls.
Family-based association tests failed to show significant results (no results p<0.05). While strictly correcting for the number of tests (α=0.002), adolescent patients with CD did not differ significantly from adolescent controls in genotype frequency for the OXTR SNPs tested; similarly, comparison of OXTR genotype frequencies for parents failed to differentiate patient and control family type, except a trend association for rs237889 (p=0.004). In this sample, 10 SNPs in the OXTR gene were not significantly associated with CD.

Leu55Pro transthyretin accelerates subunit exchange and leads to rapid formation of hybrid tetramers
CA Keetch, EHC Bromley, MG McCammon, N Wang, J Christodoulou, CV Robinson
JBC  Oct 11, 2005 M508753200. http://jbc.org/cgi/doi/10.1074/jbc.M508753200

Transthyretin is a tetrameric protein associated with the commonest form of

systemic amyloid disease. Using isotopically labeled proteins and mass spectrometry we compared subunit exchange in wild-type transthyretin with that of the variant associated with the most aggressive form of the disease, Leu55Pro. Wild-type subunit exchange occurs via both monomers and dimers , while exchange via dimers is the dominant mechanism for the Leu55Pro variant. Since patients with the Leu55Pro mutation are heterozygous, expressing both proteins simultaneously, we also analyzed the subunit exchange reaction between wild-type and Leu55Pro tetramers . We found that hybrid tetramers containing two or three Leu55Pro subunits dominate in the early stages of the reaction. Surprisingly we also found that in the presence of Leu55Pro transthyretin, the rate of dissociation of wild-type transthyretin is increased. This implies interactions between the two proteins that accelerate the formation of hybrid tetramers, a result with important implications for transthyretin amyloidos is.

Beyond Genetic Factors in Familial Amyloidotic Polyneuropathy: Protein Glycation and the Loss of Fibrinogen’s Chaperone Activity
G da Costa, RA Gomes, A Guerreiro, E Mateus, E Monteiro, et al.
PLoS ONE 2011; 6(10): e24850. http://dx.doi.org:/10.1371/journal.pone.0024850

Familial amyloidotic polyneuropathy (FAP) is a systemic conformational disease characterized by extracellular amyloid fibril formation from plasma transthyretin (TTR). This is a crippling, fatal disease for which liver transplantation is the only effective therapy. More than 80 TTR point mutations are associated with amyloidotic diseases and the most widely accepted disease model relates TTR tetramer instability with TTR point mutations. However, this model fails to explain two observations. First, native TTR also forms amyloid in systemic senile amyloidosis, a geriatric disease. Second, age at disease onset varies by decades for patients bearing the same mutation and some mutation carrier individuals are asymptomatic throughout their lives. Hence, mutations only accelerate the process and non-genetic factors must play a key role in the molecular mechanisms of disease. One of these factors is protein glycation, previously associated with conformational diseases like Alzheimer’s and Parkinson’s. The glycation hypothesis in FAP is supported by our previous discovery of methylglyoxal-derived glycation of amyloid fibrils in FAP patients. Here we show that plasma proteins are differentially glycated by methylglyoxal in FAP patients and that fibrinogen is the main glycation target. Moreover, we also found that fibrinogen interacts with TTR in plasma. Fibrinogen has chaperone activity which is compromised upon glycation by methylglyoxal. Hence, we propose that methylglyoxal glycation hampers the chaperone activity of fibrinogen, rendering TTR more prone to aggregation, amyloid formation and ultimately, disease.

Aromatic Sulfonyl Fluorides Covalently Kinetically Stabilize Transthyretin to Prevent Amyloidogenesis while Affording a Fluorescent Conjugate
NP Grimster, S Connelly, A Baranczak, J Dong, …, JW Kelly
J Am Chem Soc. 2013 Apr 17; 135(15): 5656–5668. http://dx.doi.org:/10.1021/ja311729d

Molecules that bind selectively to a given protein and then undergo a rapid chemoselective reaction to form a covalent conjugate have utility in drug development. Herein a library of 1,3,4-oxadiazoles substituted at the 2 position with an aryl sulfonyl fluoride and at the 5 position with a substituted aryl known to have high affinity for the inner thyroxine binding subsite of transthyretin (TTR) were conceived of by structure-based design principles and were chemically synthesized. When bound in the thyroxine binding site, most of the aryl sulfonyl fluorides react rapidly and chemoselectively with the pKa-perturbed K15 residue, kinetically stabilizing TTR and thus preventing amyloid fibril formation, known to cause polyneuropathy. Conjugation t50s range from 1 to 4 min, ~ 1400 times faster than the hydrolysis reaction outside the thyroxine binding site. Xray crystallography confirms the anticipated binding orientation and sheds light on the sulfonyl fluoride activation leading to the sulfonamide linkage to TTR. A few of the aryl sulfonyl fluorides efficiently form conjugates with TTR in plasma. A few of the TTR covalent kinetic stabilizers synthesized exhibit fluorescence upon conjugation and therefore could have imaging applications as a consequence of the environment sensitive fluorescence of the chromophore.

Identification of S-sulfonation and S-thiolation of a novel transthyretin Phe33Cys variant from a patient diagnosed with familial transthyretin amyloidosis
A Lim, T Prokaeva, ME Mccomb, LH Connors, M Skinner, and CE Costello
Protein Science 2003; 12:1775–1786.
http://proteinscience.org/cgi/doi/10.1110/ps.0349703.

Familial transthyretin amyloidosis (ATTR) is an autosomal dominant disorder associated with a variant form of the plasma carrier protein transthyretin (TTR). Amyloid fibrils consisting of variant TTR, wild-type TTR, and TTR fragments deposit in tissues and organs. The diagnosis of ATTR relies on the identification of pathologic TTR variants in plasma of symptomatic individuals who have biopsy proven amyloid disease. Previously, we have developed a mass spectrometry-based approach, in combination with direct DNA sequence analysis, to fully identify TTR variants. Our methodology uses immunoprecipitation to isolate TTR from serum, and electrospray ionization and matrix-assisted laser desorption/ionization mass spectrometry (MS) peptide mapping to identify TTR variants and posttranslational modifications. Unambiguous identification of the amino acid substitution is performed using tandem MS (MS/MS) analysis and confirmed by direct DNA sequence analysis. The MS and MS/MS analyses also yield information about posttranslational modifications. Using this approach, we have recently identified a novel pathologic TTR variant. This variant has an amino acid substitution (Phe — Cys) at position 33. In addition, like the Cys10 present in the wild type and in this variant, the Cys33 residue was both S-sulfonated and S-thiolated (conjugated to cysteine, cysteinylglycine, and glutathione). These adducts may play a role in the TTR fibrillogenesis.

Evolutionary relationships of lactate dehydrogenases (LDHs) from mammals, birds, an amphibian, fish, barley, and bacteria: LDH cDNA sequences from Xenopus, pig, and rat
S Tsuji, MA Qureshi, EW Hou, WM Fitch, and S S.-L. Li
Proc. Natl. Acad. Sci. USA Sep 1994; 91: 9392-9396. Evolution

The nucleotide sequences of the cDNAs encoding LDH (EC 1.1.1.27) subunits LDH-A (muscle), LDH-B (liver), and LDH-C (oocyte) from Xenopus laevis, LDH-A (muscle) and LDH-B (heart) from pig, and LDH-B (heart) and LDH-C (testis) from rat were determined. These seven newly deduced amino acid sequences and 22 other published LDH sequences, and three unpublished fish LDH-A sequences kindly provided by G. N. Somero and D. A. Powers, were used to construct the most parsimonious phylogenetic tree of these 32 LDH subunits from mammals, birds, an amphibian, fish, barley, and bacteria. There have been at least six LDH gene duplications among the vertebrates. The Xenopus LDH-A, LDH-B, and LDH-C subunits are most closely related to each other and then are more closely related to vertebrate LDH-B than LDH-A. Three fish LDH-As, as well as a single LDH of lamprey, also seem to be more related to vertebrate LDH-B than to land vertebrate LDH-A. The mammalian LDH-C (testis) subunit appears to have diverged very early, prior to the divergence of vertebrate LDH-A and LDH-B subunits, as reported previously.

Evidence for neutral and selective processes in the recruitment of enzyme-crystallins in avian lenses
Graeme Wistow, Andrea Anderson, and Joram Piatigorsky
Proc. Natl. Acad. Sci. USA Aug 1990; 87: 6277-6280, Evolution

In apparent contrast to most other tissues, the ocular lenses in vertebrates show striking differences in protein composition between taxa, most notably in the recruitment of different enzymes as major structural proteins. This variability appears to be the result of at least partially neutral evolutionary processes, although there is also evidence for selective modification in molecular structure. Here we describe a bird, the chimney swift (Chaetura pelagica), that lacks δ-crystallin/ argininosuccinate lyase, usually the major crystallin of avian lenses. Clearly, δ-crystallin is not specifically required for a functionally effective avian lens. Furthermore the lens composition of the swift is more similar to that of the related hummingbirds than to that of the barn swallow (Hirundo rustica), suggesting that phylogeny is more important than environmental selection in the recruitment of crystallins. However differences in ε-crystallin/lactate dehydrogenase-B sequence between swift and hummingbird and other avian and reptilian species suggest that selective pressures may also be working at the molecular level. These differences also confirm the close relationship between swifts and hummingbirds.

Enzyme/crystallins and extremely high pyridine nucleotide levels in the eye lens.
Zigler, J. S., Jr.; Rao, P. V.
FASEB J. 1991; 3: 223-225.

Taxon-specific crystallins are proteins present in high abundance in the lens of phylogenetically restricted groups of animals. Recently it has been found that these proteins are actually enzymes which the lens has apparently adopted to serve as structural proteins. Most of these proteins have been shown to be identical to, or related to, oxidoreductases. In guinea pig lens, which contains zeta-crystallin, a protein with an NADPH dependent oxidoreductase activity, the levels of both NADPH and NADP* are extremely high and correlate with the concentration of zeta-crystallin. We report here nucleotide assays on lenses from vertebrates containing other enzyme/crystallins. In each case where the enzyme/crystallin is a pyridine nucleotide-binding protein the level of that particular nucleotide is extremely high in the lens. The presence of an enzyme/crystallin does not affect the lenticular concentrations of those nucleotides which are not specifically bound. The possibility that nucleotide binding may be a factor in the selection of some enzymes to serve as enzyme/crystallins is considered.

Comparison of stability properties of lactate dehydrogenase B4/ε-crystallin from different species
CEM Voorter, LTM Wintjes, PWH Heinstra, H Bloemendal and WW De Jong
Eur. J. Biochem. 1993; 211: 643-648

ε-Crystallin occurs as an abundant lens protein in many birds and in crocodiles and has been identified as heart-type lactate dehydrogenase (LDH-B4). Lens proteins have, due to their longevity and environmental conditions, extraordinary requirements for structural stability. To study lens protein stability, we compared various parameters of LDH-B4/ε-crystallin from lens and/or heart of duck, which has abundant amounts of this enzyme in its lenses, and of chicken and pig, which have no λ-crystallin. Measuring the thermostability of LDH-B4 from the different sources, the t50 values (temperature at which 50% of the enzyme activity remains after a 20-min period) for LDH-B4 from duck heart, duck lens and chicken heart were all found to be around 76°C whereas pig heart LDHB4 was less thermostable, having a t50 value of 625°C. A similar tendency was found with urea inactivation studies. Plotting the first-order rate constants obtained from inactivation kinetic plots against urea concentration, it was clear that LDH-B4 from pig heart was less stable in urea than the homologous enzymes from duck heart, chicken heart and duck lens. The duck and chicken enzymes were also much more resistant against proteolysis than the porcine enzyme. Therefore, it is concluded that avian LDH-B4 is structurally more stable than the homologous enzyme in mammals. This greater stability might make it suitable to function as a ε-crystallin, as in duck, but is not necessarily associated with high lens expression, as in chicken.

Duck lens ε-crystallin and lactate dehydrogenase B4 are identical: A single-copy gene product with two distinct functions
W Hendriks, JWM Mulders, MA Bibby, C Slingsby, H Bloemendal, and WW De Jong
Proc. Natl. Acad. Sci. USA Oct 1988; 85: 7114-7118. Biochemistry

To investigate whether or not duck lens ε-crystaliin and duck heart lactate dehydrogenase (LDH) B4 are the product of the same gene, we have isolated and sequenced cDNA clones of duck ε-crystallin. By using these clones we demonstrate that there is a single-copy Ldh-B gene in duck and in chicken. In the duck lens this gene is overexpressed, and its product is subject to posttranslational modification. Reconstruction of the evolutionary history of the LDH protein family reveals that the mammalian Ldh-C gene most probably originated from an ancestral Ldh-A gene and that the amino acid replacement rate in LDH-C is approximately 4 times the rate in LDH-A. Molecular modeling of LDH-B sequences shows that the increased thermostability of the avian tetramer might be explained by mutations that increase the number of ion pairs. Furthermore, the replacement of bulky side chains by glycines on the corners of the duck protein suggests an adaptation to facilitate close packing in the lens.

Lactate Dehydrogenase A as a Highly Abundant Eye Lens Protein in Platypus (Ornithorhynchus anatinus): Upsilon (υ)-Crystallin
T van Rheede,  R Amons, N Stewart, and WW de Jong
Mol. Biol. Evol. 2003; 20(06):994–998. http://dx.doi.org:/10.1093/molbev/msg116

Vertebrate eye lenses mostly contain two abundant types of proteins, the α-crystallins and the β/λ-crystallins. In addition, certain housekeeping enzymes are highly expressed as crystallins in various taxa. We now observed an unusual approximately 41-kd protein that makes up 16% to 18% of the total protein in the platypus eye lens. Its cDNA sequence was determined, which identified the protein as muscle-type lactate dehydrogenase A (LDH-A). It is the first observation of LDH-A as a crystallin, and we designate it upsilon (υ)-crystallin. Interestingly, the related heart-type LDH-B occurs as an abundant lens protein, known as ε-crystallin, in many birds and crocodiles. Thus, two members of the ldh gene family have independently been recruited as crystallins in different higher vertebrate lineages, suggesting that they are particularly suited for this purpose in terms of gene regulatory or protein structural properties. To establish whether platypus LDH-A/υ-crystallin has been under different selective constraints as compared with other vertebrate LDH-A sequences, we reconstructed the vertebrate Ldh-A gene phylogeny. No conspicuous rate deviations or amino acid replacements were observed.

Isozymes, moonlighting proteins and promiscous enzymes
M Nath Gupta, M Kapoor, AB Majumder and V Singh
Current Science Apr 2011; 100(8): 1152-1162.

The structures of isoenzymes differ and yet these catalyse the same type of reaction. These structures evolved to suit the physiological needs and are located in different parts of cells or tissues. Moonlighting proteins represent the same structure performing very different biological functions. Biological promiscuity reveals that the same active sites can catalyse different types of reactions. These three different phenomena, all illustrate similar evolutionary strategies. Viewed together, it emerges that biologists need to take a hard look at the ‘structure–function’ paradigm as well as the notions of biological specificity. Meanwhile, biotechnologists  continue to exploit the opportunities which ‘nonspecificity’ offers.

Read Full Post »

« Newer Posts - Older Posts »