Funding, Deals & Partnerships: BIOLOGICS & MEDICAL DEVICES; BioMed e-Series; Medicine and Life Sciences Scientific Journal – http://PharmaceuticalIntelligence.com
Genetic Testing in CVD and Precision Medicine, Volume 2 (Volume Two: Latest in Genomics Methodologies for Therapeutics: Gene Editing, NGS and BioInformatics, Simulations and the Genome Ontology), Part 1: Next Generation Sequencing (NGS)
In 2010, we introduced an approach to the evaluation of a personal genome in a clinical context (7). A patient with a family history of coronary artery disease (CAD) and sudden death was evaluated by a cardiac clinical team in conjunction with whole genome sequencing and interpretation. The genomic analysis revealed an increased genetic risk for myocardial infarction and type 2 diabetes. In addition, a pharmacogenomics analysis was performed to assess how the genetics of the patient might influence response to certain drugs, including lipid-lowering therapies and warfarin (7). This clinical assessment, which focused heavily on cardiovascular risk, suggested that whole genome sequencing might provide clinically relevant information for patients.
A 2011 joint statement from the Heart Rhythm Society and the European Heart Rhythm association recommended genetic testing as a class I indication for patients with a number of channelopathies and cardiomyopathies, including long QT syndrome (LQTS), arrhythmogenic right ventricular cardiomyopathy, familial dilated cardiomyopathy (DCM), and hypertrophic cardiomyopathy (HCM) (8). Similarly, a statement from the American Heart Association and the American College of Cardiology recommended genetic testing for HCM, DCM, and thoracic aortic aneurysms to facilitate familial cascade screening and deduce causative mutations 9, 10.
The diagnostic power of genetic testing is significant across the spectrum of CVDs, ranging from cardiomyopathies to life-threatening arrhythmias 10, 11, 12. In the clinic, genetic testing can:
1.
clarify disease diagnoses: genetic testing can help to clarify the diagnosis of diseases that cause similar clinical presentation (e.g., cardiac hypertrophy could be TTR amyloidosis, Fabry disease, or sarcomeric HCM);
2.
facilitate cascade screening: genetic testing can help to identify relatives at risk for CVD before disease symptoms manifest if a disease-associated variant is found in a proband and then screened for in relatives;
3.
direct more precise therapy: genetic testing can help physicians choose appropriate treatments and plan appropriate timing of those treatments. For example, inherited connective tissue disease due to variants in ACTA2, MYH11, or TGFBR2 might prompt consideration of surgical intervention at a smaller aortic aneurysm diameter (13); and
4.
identify patients for targeted therapies: targeted medical therapies, including antibody-based therapeutics, gene editing, and silencing technologies, are available or under development for several genetic diseases, including LQTS, Duchenne muscular dystrophy (DMD), TTR cardiac amyloidosis (14), and Fabry disease 13, 15.
8. Ackerman M.J., Priori S.G., Willems S. HRS/EHRA expert consensus statement on the state of genetic testing for the channelopathies and cardiomyopathies: this document was developed as a partnership between the Heart Rhythm Society (HRS) and the European Heart Rhythm Association (EHRA) Europace. 2011;13:1077–1109. [PubMed] [Google Scholar]
9. Gersh B.J., Maron B.J., Bonow R.O. 2011 ACCF/AHA guideline for the diagnosis and treatment of hypertrophic cardiomyopathy: executive summary: a report of the American College of Cardiology Foundation/American Heart Association Task Force on Practice Guidelines. J Am Coll Cardiol. 2011;58:2703–2738. [PubMed] [Google Scholar]
10. Harper A.R., Parikh V.N., Goldfeder R.L., Caleshu C., Ashley E.A. Delivering clinical grade sequencing and genetic test interpretation for cardiovascular medicine. Circ Cardiovasc Genet. 2017;10(2) [PubMed] [Google Scholar]
11. Walsh R., Thomson K.L., Ware J.S. Reassessment of Mendelian gene pathogenicity using 7,855 cardiomyopathy cases and 60,706 reference samples. Genet Med. 2017;19:192–203. [PMC free article] [PubMed] [Google Scholar]
12. Sturm A.C., Hershberger R.E. Genetic testing in cardiovascular medicine: current landscape and future horizons. Curr Opin Cardiol. 2013;28:317–325. [PubMed] [Google Scholar]
13. Caleshu C., Ashley E. Genetic testing for cardiovascular conditions predisposing to sudden death. In: Wilson M.G., Drezner J., editors. IOC Manual of Sports Cardiology. Wiley & Sons, Ltd; Hoboken, NJ: 2016. pp. 175–186. [Google Scholar]
14. Benson M.D., Dasgupta N.R., Rissing S.M., Smith J., Feigenbaum H. Safety and efficacy of a TTR specific antisense oligonucleotide in patients with transthyretin amyloid cardiomyopathy. Amyloid. 2017;24:217–223. [PubMed] [Google Scholar]
15. Parikh V.N., Ashley E.A. Next-generation sequencing in cardiovascular disease: present clinical applications and the horizon of precision medicine. Circulation. 2017;135:406–409. [PMC free article] [PubMed] [Google Scholar]
Individuals without angiographic CAD but with hiPRS remain at significantly elevated risk of mortality after cardiac catheterization
Reporter: Aviva Lev-Ari, PhD, RN
A genome-wide Polygenic risk scores (PRS) improves risk stratification when added to traditional risk factors and coronary angiography. Individuals without angiographic CAD but with hiPRS remain at significantly elevated risk of mortality.
Originally published12 Nov 2018 Circulation: Genomic and Precision Medicine. 2018;11:e002352
Abstract
Background:
Coronary artery disease (CAD) is influenced by genetic variation and traditional risk factors. Polygenic risk scores (PRS), which can be ascertained before the development of traditional risk factors, have been shown to identify individuals at elevated risk of CAD. Here, we demonstrate that a genome-wide PRS for CAD predicts all-cause mortality after accounting for not only traditional cardiovascular risk factors but also angiographic CAD itself.
Methods:
Individuals who underwent coronary angiography and were enrolled in an institutional biobank were included; those with prior myocardial infarction or heart transplant were excluded. Using a pruning-and-thresholding approach, a genome-wide PRS comprised of 139 239 variants was calculated for 1503 participants who underwent coronary angiography and genotyping. Individuals were categorized into high PRS (hiPRS) and low-PRS control groups using the maximally selected rank statistic. Stratified analysis based on angiographic findings was also performed. The primary outcome was all-cause mortality following the index coronary angiogram.
Results:
Individuals with hiPRS were younger than controls (66 years versus 69 years; P=2.1×10-5) but did not differ by sex, body mass index, or traditional risk-factor profiles. Individuals with hiPRS were at significantly increased risk of all-cause mortality after cardiac catheterization, adjusting for traditional risk factors and angiographic extent of CAD (hazard ratio, 1.6; 95% CI, 1.2–2.2; P=0.004). The strongest increase in risk of all-cause mortality conferred by hiPRS was seen among individuals without angiographic CAD (hazard ratio, 2.4; 95% CI, 1.1–5.5; P=0.04). In the overall cohort, adding hiPRS to traditional risk assessment improved prediction of 5-year all-cause mortality (area under the receiver-operating curve 0.70; 95% CI, 0.66–0.75 versus 0.66; 95% CI, 0.61–0.70; P=0.001).
Conclusions:
A genome-wide PRS improves risk stratification when added to traditional risk factors and coronary angiography. Individuals without angiographic CAD but with hiPRS remain at significantly elevated risk of mortality.
Scott M. Damrauer, MD, Department of Surgery, Hospital of the University of Pennsylvania, 3400 Spruce St, Silverstein 4, Philadelphia, PA 19104. Email Scott.Damrauer@uphs.upenn.edu
The HFE H63D variant confers an increased risk for hypertension, no increased risk for adverse cardiovascular events or substantial left ventricular remodeling
Reporter: Aviva Lev-Ari, PhD, RN
Conclusion:
The HFE H63D variant confers an increased risk for hypertension per allele and, given its frequency, accounts for a significant number of cases of hypertension. However, there was no increased risk for adverse cardiovascular events or substantial left ventricular remodeling.
HFE H63D Polymorphism and the Risk for Systemic Hypertension, Myocardial Remodeling, and Adverse Cardiovascular Events in the ARIC Study
Originally published12 Nov 2018Hypertension. 2018;0:HYPERTENSIONAHA.118.11730
Abstract
H63D has been identified as a novel locus associated with the development of hypertension. The quantitative risks for hypertension, cardiac remodeling, and adverse events are not well studied. We analyzed white participants from the ARIC study (Atherosclerosis Risk in Communities) with H63D genotyping (N=10 902). We related genotype status to prevalence of hypertension at each of 5 study visits and risk for adverse cardiovascular events. Among visit 5 participants (N=4507), we related genotype status to echocardiographic features. Frequencies of wild type (WT)/WT, H63D/WT, and H63D/H63D were 73%, 24.6%, and 2.4%. The average age at baseline was 54.9±5.7 years and 47% were men. Participants carrying the H63D variant had higher systolic blood pressure (P=0.004), diastolic blood pressure (0.012), and more frequently had hypertension (P<0.001). Compared with WT/WT, H63D/WT and H63D/H63D participants had a 2% to 4% and 4% to 7% absolute increase in hypertension risk at each visit, respectively. The population attributable risk of H63D for hypertension among individuals aged 45 to 64 was 3.2% (95% CI, 1.3–5.1%) and 1.3% (95% CI, 0.0–2.4%) among individuals >65 years. After 25 years of follow-up, there was no relationship between genotype status and any outcome (P>0.05). H63D/WT and H63D/H63D genotypes were associated with small differences in cardiac remodeling. In conclusion, the HFE H63D variant confers an increased risk for hypertension per allele and, given its frequency, accounts for a significant number of cases of hypertension. However, there was no increased risk for adverse cardiovascular events or substantial left ventricular remodeling.
Correspondence to Scott D. Solomon, Department of Medicine, Brigham and Women’s Hospital, Harvard Medical School, 75 Francis St, Boston, MA 02115. Email ssolomon@rics.bwh.harvard.edu
Aortic Stenosis (AS): Managed Surgically by Transcatheter Aortic Valve Replacement (TAVR) – Search Results for “TAVR” on NIH.GOV website, Top 16 pages
Reporter: Aviva Lev-Ari, PhD, RN
UPDATED on 9/24/2018
Sapien 3, CoreValve Evolut R on Par for Aortic Stenosis
Head-to-head trial also shows local, general anesthesia outcomes similar
by Ashley Lyles, Staff Writer, MedPage Today
This article is a collaboration between MedPage Today® and:
SAN DIEGO — Transfemoral transcatheter aortic valve replacement (TAVR) with the balloon-expandable Edwards Sapien 3 valve yields the same early outcomes as the self-expanding CoreValve Evolut R, regardless of anesthesia strategy, a two-by-two randomized trial showed.
In the valve comparison, the primary endpoint of all-cause mortality, stroke, moderate or severe prosthetic valve regurgitation, and permanent pacemaker implantation at 30 days met criteria for equivalence, with a composite rate of 27.2% with Evolut R and 26.1% with Sapien 3, Holger Thiele, MD, of University Hospital in Leipzig, Germany, reported here at the Transcatheter Cardiovascular Therapeutics meeting.
The researchers also evaluated the effects of anesthesia used during these procedures and found no significant difference. The composite endpoint at 30 days came out 27.0% for local anesthesia and 25.5% for general anesthesia.
“The SOLVE-TAVI trial is the first adequately powered randomized trial comparing local versus general anesthesia in patients with symptomatic aortic valve stenosis undergoing TAVR,” said Thiele in a press release. “Results indicate that local anesthesia is both safe and effective and may be a good option for those patients undergoing TAVR with an intermediate or high surgical risk.”
In the majority of aortic stenosis cases, it doesn’t matter which valve you choose, although there are still some cases, like heavy calcification, when it may be better to choose one valve over the other, noted panel discussant Molly Szerlip, MD, of Baylor Scott & White The Heart Group in McKinney, Texas.
The researchers evaluated 447 patients who were receiving care at German medical centers for severe symptomatic aortic stenosis and were at an intermediate- to high-surgical risk. The patients were randomized to have the Sapien 3 valve or CoreValve Evolut R and to either receive general or local anesthesia with conscious sedation.
The individual valve strategy findings again showed equivalence without superiority between Evolut R and Sapien 3 for mortality (2.8% vs 2.3%) and moderate or severe valve regurgitation (1.9% vs 1.4%). But for stroke Evolut R came out superior (0.5% vs 4.7%), and the two didn’t meet criteria for equivalence on pacemaker implantation (22.9% vs 19.0%, P=0.06 for equivalence).
“The rate of relevant valve regurgitation was low whereas permanent pacemaker rates are still relatively high,” the researchers wrote.
The anesthesia comparison endpoints all met the criteria for equivalence without superiority of general anesthesia over local anesthesia:
The concept of transcatheter balloon expandable valves was first introduced in the 1980s by a Danish researcher by the name of H. R. Anderson who began testing this idea on pigs. In 2002, Dr. Alain Cribier performed the first successful percutaneous aortic valve replacement on an inoperable patient. The first approval of TAVR for the indication of severe AS in prohibitive risk patients came in 2011. In 2012, the FDA approved TAVR in patients at high surgical risk. In 2015 the indication was expanded to include “valve-in-valve” procedure for failed surgical bioprosthetic valves. Most recently, in 2016 the FDA approved the SAPIEN valve for use in patients with severe AS at intermediate risk.
4 months ago – Scientists developed a novel technique that prevents a rare but often fatal complication that can arise during a heart valve procedure called …
Transcaval TAVR was developed at the NHLBI Cardiovascular Intervention Program and applied to patient care in collaboration with Dr. Adam Greenbaum at …
Mandatory Reporting of Clinical Trial Identifier Numbers on Claims . Q: Do organizations bill Medicare for all services related to the clinical trial …
de Biasi AR, Worku B, Skubas NJ, Salemi A. Transcatheter aortic valve replacement (TAVR) continues to garner considerable attention, especially as the …
Díez JG, Schechter M, Dougherty KG, Preventza O, Coselli JS. Transcatheter aortic valve replacement (TAVR) is a well-established method for replacing …
PubMed comprises more than 28 million citations for biomedical literature from MEDLINE, life science journals, and online books. Citations may include …
Cardiac rehabilitation after transcatheter aortic valve implantation compared to patients after valve replacement. Tarro Genta F(1), Tidu M, Bouslenko …
Heart surgery is done to correct problems with the heart. Many heart surgeries are done each year in the United States for various heart problems. The …
NIH study also shows that overall stroke risk is down from 10 years ago. Aspirin combined with the antiplatelet drug clopidogrel is no better than asp …
NIH study also shows that overall stroke risk is down from 10 years ago. Aspirin combined with the antiplatelet drug clopidogrel is no better than asp …
Cardiac catheterization (KATH-eh-ter-ih-ZA-shun) is a medical procedure used to diagnose and treat some heart conditions. A long, thin, flexible tube …
TCT-712 “Cusp Overlap” View Facilitates Accurate Fluoro-Guided Implantation of Self-Expanding Valve in TAVR. Zaid S, Raza A, Michev I, Ahmad H, Kaple …
Balloon expandable sheath for transfemoral aortic valve implantation: a viable option for patients with challenging access. Dimitriadis Z(1), Scholtz …
Author information: (1)Hôpital du Sacré-Coeur de Montréal, Montréal, Québec, Canada; Morristown Medical Center, Morristown, New Jersey; Cardiovascular …
Intraprocedural TAVR Annulus Sizing Using 3D TEE and the “Turnaround Rule”. Wiley BM, Kovacic JC, Basnet S, Makoto A, Chaudhry FA, Kini AS, Sharma SK …
1. Curr Pharm Des. 2016;22(13):1911-7. Bivalirudin or Heparin Treatment During Transcatheter Valve Interventions: Where are we Now? Mastoris I, Dangas …
Author information: (1)Hôpital du Sacré-Coeur de Montréal, Montréal, Québec, Canada; Morristown Medical Center, Morristown, New Jersey; Cardiovascular …
Intraprocedural TAVR Annulus Sizing Using 3D TEE and the “Turnaround Rule”. Wiley BM, Kovacic JC, Basnet S, Makoto A, Chaudhry FA, Kini AS, Sharma SK …
Clinical Trial Registry Numbers in MEDLINE ® /PubMed ® Records. A clinical trial is a research study using human volunteers to answer specific health …
Zhao A(1), Minhui H(2), Li X(1), Zhiyun X(1). Author information: (1)Department of Cardiothoracic Surgery, Changhai Hospital, Second Military Medical …
Zhao A(1), Minhui H(2), Li X(1), Zhiyun X(1). Author information: (1)Department of Cardiothoracic Surgery, Changhai Hospital, Second Military Medical …
PubMed comprises more than 26 million citations for biomedical literature from MEDLINE, life science journals, and online books. Citations may include …
Etiologies of Cardiovascular Diseases: Epigenetics, Genetics and Genomics: Request for Book Review Writing on Amazon.com, Volume 2 (Volume Two: Latest in Genomics Methodologies for Therapeutics: Gene Editing, NGS and BioInformatics, Simulations and the Genome Ontology), Part 1: Next Generation Sequencing (NGS)
1.2 Ribozymes and RNA Machines – Work of Jennifer A. Doudna
1.3 Genomics and Medicine: Contributions of Genetics and Genomics to Cardiovascular Disease Diagnoses
1.4 Genomics Orientations for Individualized Medicine, Volume One
1.4.1 CVD Epidemiology, Ethnic subtypes Classification, and Medication Response Variability: Cardiology, Genomics and Individualized Heart Care: Framingham Heart Study (65 y-o study) & Jackson Heart Study (15 y-o study)
1.4.2 What comes after finishing the Euchromatic Sequence of the Human Genome?
1.5 Genomics in Medicine – Establishing a Patient-Centric View of Genomic Data
PART 2
Epigenetics – Modifiable Factors Causing Cardiovascular Diseases
2.1 Diseases Etiology
2.1.1 Environmental Contributors Implicated as Causing Cardiovascular Diseases
2.1.2 Diet: Solids, Fluid Intake and Nutraceuticals
2.1.3 Physical Activity and Prevention of Cardiovascular Diseases
2.1.4 Psychological Stress and Mental Health: Risk for Cardiovascular Diseases
2.1.5 Correlation between Cancer and Cardiovascular Diseases
2.1.6 Medical Etiologies for Cardiovascular Diseases: Evidence-based Medicine – Leading DIAGNOSES of Cardiovascular Diseases, Risk Biomarkers and Therapies
2.1.7 Signaling Pathways
2.1.8 Proteomics and Metabolomics
2.1.9 Sleep and Cardiovascular Diseases
2.2 Assessing Cardiovascular Disease with Biomarkers
2.2.1 Issues in Genomics of Cardiovascular Diseases
2.2.2 Endothelium, Angiogenesis, and Disordered Coagulation
2.2.3 Hypertension BioMarkers
2.2.4 Inflammatory, Atherosclerotic and Heart Failure Markers
2.2.5 Myocardial Markers
2.3 Therapeutic Implications: Focus on Ca(2+) signaling, platelets, endothelium
2.3.1 The Centrality of Ca(2+) Signaling and Cytoskeleton Involving Calmodulin Kinases and Ryanodine Receptors in Cardiac Failure, Arterial Smooth Muscle, Post-ischemic Arrhythmia, Similarities and Differences, and Pharmaceutical Targets
2.3.2 EMRE in the Mitochondrial Calcium Uniporter Complex
2.3.3 Platelets in Translational Research 2: Discovery of Potential Anti-platelet Targets
2.3.4 The Final Considerations of the Role of Platelets and Platelet Endothelial Reactions in Atherosclerosis and Novel Treatments
2.3.5 Nitric Oxide Synthase Inhibitors (NOS-I)
2.3.6 Resistance to Receptor of Tyrosine Kinase
2.3.7 Oxidized Calcium Calmodulin Kinase and Atrial Fibrillation
2.3.8 Advanced Topics in Sepsis and the Cardiovascular System at its End Stage
2.4 Comorbidity of Diabetes and Aging
2.4.1 Heart and Aging Research in Genomic Epidemiology: 1700 MIs and 2300 coronary heart disease events among about 29 000 eligible patients
2.4.2 Pathophysiological Effects of Diabetes on Ischemic-Cardiovascular Disease and on Chronic Obstructive Pulmonary Disease (COPD)
2.4.3 Risks of Hypoglycemia in Diabetics with Chronic Kidney Disease (CKD)
2.4.4 Mitochondrial Mechanisms of Disease in Diabetes Mellitus
2.4.5 Mitochondria: More than just the “powerhouse of the cell”
2.4.6 Pathophysiology of GLP-1 in Type 2 Diabetes
2.4.7 Developments in the Genomics and Proteomics of Type 2 Diabetes Mellitus and Treatment Targets
2.4.8 CaKMII Inhibition in Obese, Diabetic Mice leads to Lower Blood Glucose Levels
2.4.9 Protein Target for Controlling Diabetes, Fractalkine: Mediator cell-to-cell Adhesion though CX3CR1 Receptor, Released from cells Stimulate Insulin Secretion
2.4.10 Peroxisome proliferator-activated receptor (PPAR-gamma) Receptors Activation: PPARγ transrepression for Angiogenesis in Cardiovascular Disease and PPARγ transactivation for Treatment of Diabetes
2.4.11 CABG or PCI: Patients with Diabetes – CABG Rein Supreme
2.4.12 Reversal of Cardiac Mitochondrial Dysfunction
2.4.13 BARI 2D Trial Outcomes
2.4.14 Overview of new strategy for treatment of T2DM: SGLT2 inhibiting oral antidiabetic agents
2.5 Drug Toxicity and Cardiovascular Diseases
2.5.1 Predicting Drug Toxicity for Acute Cardiac Events
2.5.2 Cardiotoxicity and Cardiomyopathy Related to Drugs Adverse Effects
2.5.3 Decoding myocardial Ca2+ signals across multiple spatial scales: A role for sensitivity analysis
2.5.4. Leveraging Mathematical Models to Understand Population Variability in Response to Cardiac Drugs: Eric Sobie, PhD
2.5.5 Exploiting mathematical models to illuminate electrophysiological variability between individuals.
2.5.6 Clinical Effects and Cardiac Complications of Recreational Drug Use: Blood pressure changes, Myocardial ischemia and infarction, Aortic dissection, Valvular damage, and Endocarditis, Cardiomyopathy, Pulmonary edema and Pulmonary hypertension, Arrhythmias, Pneumothorax and Pneumopericardium
2.6 Male and Female Hormonal Replacement Therapy: The Benefits and the Deleterious Effects on Cardiovascular Diseases
2.6.1 Testosterone Therapy for Idiopathic Hypogonadotrophic Hypogonadism has Beneficial and Deleterious Effects on Cardiovascular Risk Factors
2.6.2 Heart Risks and Hormones (HRT) in Menopause: Contradiction or Clarification?
2.6.3 Calcium Dependent NOS Induction by Sex Hormones: Estrogen
2.6.4 Role of Progesterone in Breast Cancer Progression
PART 3
Determinants of Cardiovascular Diseases Genetics, Heredity and Genomics Discoveries
Introduction
3.1 Why cancer cells contain abnormal numbers of chromosomes (Aneuploidy)
3.1.1 Aneuploidy and Carcinogenesis
3.2 Functional Characterization of Cardiovascular Genomics: Disease Case Studies @ 2013 ASHG
3.3 Leading DIAGNOSES of Cardiovascular Diseases covered in Circulation: Cardiovascular Genetics, 3/2010 – 3/2013
3.3.9: Genetics and Vascular Pathologies and Platelet Aggregation, Cardiac Troponin T in Serum
3.3.10: Genomics and Valvular Disease
3.4 Commentary on Biomarkers for Genetics and Genomics of Cardiovascular Disease
PART 4
Individualized Medicine Guided by Genetics and Genomics Discoveries
4.1 Preventive Medicine: Cardiovascular Diseases
4.1.1 Personal Genomics for Preventive Cardiology Randomized Trial Design and Challenges
4.2 Gene-Therapy for Cardiovascular Diseases
4.2.1 Genetic Basis of Cardiomyopathy
4.3 Congenital Heart Disease/Defects
4.4 Cardiac Repair: Regenerative Medicine
4.4.1 A Powerful Tool For Repairing Damaged Hearts
4.4.2 Modified RNA Induces Vascular Regeneration After a Heart
4.5 Pharmacogenomics for Cardiovascular Diseases
4.5.1 Blood Pressure Response to Antihypertensives: Hypertension Susceptibility Loci Study
4.5.2 Statin-Induced Low-Density Lipoprotein Cholesterol Reduction: Genetic Determinants in the Response to Rosuvastatin
4.5.3 SNPs in apoE are found to influence statin response significantly. Less frequent variants in PCSK9 and smaller effect sizes in SNPs in HMGCR
4.5.4 Voltage-Gated Calcium Channel and Pharmacogenetic Association with Adverse Cardiovascular Outcomes: Hypertension Treatment with Verapamil SR (CCB) vs Atenolol (BB) or Trandolapril (ACE)
4.5.5 Response to Rosuvastatin in Patients With Acute Myocardial Infarction: Hepatic Metabolism and Transporter Gene Variants Effect
4.5.6 Helping Physicians identify Gene-Drug Interactions for Treatment Decisions: New ‘CLIPMERGE’ program – Personalized Medicine @ The Mount Sinai Medical Center
4.5.7 Is Pharmacogenetic-based Dosing of Warfarin Superior for Anticoagulation Control?
Lysyl Oxidase (LOX) gene missense mutation causes Thoracic Aortic Aneurysm and Dissection (TAAD) in Humans because of inadequate cross-linking of collagen and elastin in the aortic wall
Mutation carriers may be predisposed to vascular diseases because of weakened vessel walls under stress conditions.
Reporter: Aviva Lev-Ari, PhD, RN
2.1.3.7 Lysyl Oxidase (LOX) gene missense mutation causes Thoracic Aortic Aneurysm and Dissection (TAAD) in Humans because of inadequate cross-linking of collagen and elastin in the aortic wall – Mutation carriers may be predisposed to vascular diseases because of weakened vessel walls under stress conditions, Volume 2 (Volume Two: Latest in Genomics Methodologies for Therapeutics: Gene Editing, NGS and BioInformatics, Simulations and the Genome Ontology), Part 2: CRISPR for Gene Editing and DNA Repair
Loss of function mutation in LOX causes thoracic aortic aneurysm and dissection in humans
Author contributions: V.S.L., R.P.M., N.Y.F., and N.O.S. designed research; V.S.L., C.M.H., and N.O.S. performed research; E.P.H., N.C., C.G.L., D.V., B.G.M.P., R.P.M., and N.Y.F. contributed new reagents/analytic tools; V.S.L., C.M.
Significance
The mechanical integrity of the arterial wall is dependent on a properly structured ECM. Elastin and collagen are key structural components of the ECM, contributing to the stability and elasticity of normal arteries. Lysyl oxidase (LOX) normally cross-links collagen and elastin molecules in the process of forming proper collagen fibers and elastic lamellae. Here, using whole-genome sequencing in humans and genome engineering in mice, we show that a missense mutation in LOX causes aortic aneurysm and dissection because of insufficient elastin and collagen cross-linking in the aortic wall. These findings confirm mutations in LOX as a cause of aortic disease in humans and identify LOX as a diagnostic and potentially therapeutic target.
Abstract
Thoracic aortic aneurysms and dissections (TAAD) represent a substantial cause of morbidity and mortality worldwide. Many individuals presenting with an inherited form of TAAD do not have causal mutations in the set of genes known to underlie disease. Using whole-genome sequencing in two first cousins with TAAD, we identified a missense mutation in the lysyl oxidase (LOX) gene (c.893T > G encoding p.Met298Arg) that cosegregated with disease in the family. Using clustered regularly interspaced short palindromic repeats (CRISPR)/clustered regularly interspaced short palindromic repeats-associated protein-9 nuclease (Cas9) genome engineering tools, we introduced the human mutation into the homologous position in the mouse genome, creating mice that were heterozygous and homozygous for the human allele. Mutant mice that were heterozygous for the human allele displayed disorganized ultrastructural properties of the aortic wall characterized by fragmented elastic lamellae, whereas mice homozygous for the human allele died shortly after parturition from ascending aortic aneurysm and spontaneous hemorrhage. These data suggest that a missense mutation in LOX is associated with aortic disease in humans, likely through insufficient cross-linking of elastin and collagen in the aortic wall. Mutation carriers may be predisposed to vascular diseases because of weakened vessel walls under stress conditions. LOX sequencing for clinical TAAD may identify additional mutation carriers in the future. Additional studies using our mouse model of LOX-associated TAAD have the potential to clarify the mechanism of disease and identify novel therapeutics specific to this genetic cause.
NEW YORK (GenomeWeb) – Researchers from Washington University School of Medicine have linked a LOX gene variant with aortic rupture and aneurysm.
As they reported in the online early edition of the Proceedings of the National Academy of Sciences yesterday, the researchers sequenced two first cousins from a family with a history of aortic ruptures and aneurysms to uncover a missense mutation in the lysyl oxidase (LOX) gene, which encodes a protein that cross-links elastin and collagen. When they used CRISPR/Cas9 genome engineering to introduce the mutation into a mouse model, mice heterogeneous for the mutation had disorganized aortic walls, while mice homozygous for the mutation died shortly after birth of ascending aneurysm and spontaneous hemorrhage, suggesting that the LOX variant might be causal.
Human Genetics and Childhood Diseases, Volume 2 (Volume Two: Latest in Genomics Methodologies for Therapeutics: Gene Editing, NGS and BioInformatics, Simulations and the Genome Ontology), Part 1: Next Generation Sequencing (NGS)
Human Genetics and Childhood Diseases
Curator: Larry H. Bernstein, MD, FCAP
Publication Roundup: HGMD
HGMD®, the Human Gene Mutation Database is used by scientists around the world to find information on reported genetic mutations. The papers below use the database to advance our understanding of disease, DNA dynamics, and more.
Local DNA dynamics shape mutational patterns of mononucleotide repeats in human genomes First author: Albino Bacolla
Scientists in the US and UK published results in Nucleic Acids Research of a detailed analysis of single-base substitutions and indels in the human genome. Their findings show that certain base positions are more susceptible to mutagenesis than others. They used HGMD Professional to find mutations in specific genomic regions for analysis; the paper includes charts showing mutation patterns, germline SNPs, and more from HGMD data.
High prevalence of CDH23 mutations in patients with congenital high-frequency sporadic or recessively inherited hearing loss First author: Kunio Mizutari
This Orphanet Journal of Rare Diseases paper from scientists in Japan sequenced 72 patients with unexplained hearing loss, finding several CDH23 mutations, some of which were novel. Mutations in the gene have been linked to Usher syndrome and other forms of hereditary hearing loss. The scientists used HGMD to find all known CDH23 mutations within nearly 70 coding regions.
Mutation analyses and prenatal diagnosis in families of X-linked severe combined immunodeficiency caused by IL2Rγ gene novel mutation First author: Q.L. Bai
In Genetics and Molecular Research, scientists report the utility of mutation analysis of the interleukin-2 receptor gamma gene to assess carrier status and perform prenatal diagnosis for X-linked severe combined immunodeficiency. They studied two high-risk families, along with 100 controls, to evaluate the approach. Sequence variation was determined using HGMD Professional and an X-SCID database, and a new mutation was discovered in the project.
Impact of glucocerebrosidase mutations on motor and nonmotor complications in Parkinson’s disease First author: Tomoko Oeda
Researchers from three hospitals in Japan published this Neurobiology of Aging report that may help stratify Parkinson’s disease patients by prognosis. They sequenced mutations in the GBA gene in 215 patients, finding that those who had mutations associated with Gaucher disease suffered dementia and psychosis much earlier than those who didn’t. The team found previously reported GBA mutations using HGMD Professional.
Comprehensive Genetic Characterization of a Spanish Brugada Syndrome Cohort First author: Elisabet Selga
In this PLoS One publication, scientists from a number of institutions in Spain examined genetic variation among patients with Brugada syndrome, a rare genetic cardiac arrhythmia. They sequenced 14 genes in 55 patients, identifying 61 variants and finding the subset that appear pathogenic. Variants were filtered against a number of databases, including HGMD.
Local DNA dynamics shape mutational patterns of mononucleotide repeats in human genomes
Single base substitutions (SBSs) and insertions/deletions are critical for generating population diversity and can lead both to inherited disease and cancer. Whereas on a genome-wide scale SBSs are influenced by cellular factors, on a fine scale SBSs are influenced by the local DNA sequence-context, although the role of flanking sequence is often unclear. Herein, we used bioinformatics, molecular dynamics and hybrid quantum mechanics/molecular mechanics to analyze sequence context-dependent mutagenesis at mononucleotide repeats (A-tracts and G-tracts) in human population variation and in cancer genomes. SBSs and insertions/deletions occur predominantly at the first and last base-pairs of A-tracts, whereas they are concentrated at the second and third base-pairs in G-tracts. These positions correspond to the most flexible sites along A-tracts, and to sites where a ‘hole’, generated by the loss of an electron through oxidation, is most likely to be localized in G-tracts. For A-tracts, most SBSs occur in the direction of the base-pair flanking the tracts. We conclude that intrinsic features of local DNA structure, i.e. base-pair flexibility and charge transfer, render specific nucleotides along mononucleotide runs susceptible to base modification, which then yields mutations. Thus, local DNA dynamics contributes to phenotypic variation and disease in the human population.
INTRODUCTION
Changes in human genomic DNA in the form of base substitutions and insertions/deletions (indels) are essential to ensure population diversity, adaptation to the environment, defense from pathogens and self-recognition; they are also a critical source of human inherited disease and cancer. On a genome-wide scale, base substitutions result from the combined action of several factors, including replication fidelity, lagging versus leading strand DNA synthesis, repair, recombination, replication timing, transcription, nucleosome occupancy, etc., both in the germline and in cancer (1–4). On a much finer scale [(over a few base pairs (bp)], rates of base substitutions may be strongly influenced by interrelationships between base–protein and base–base interactions. For example, the mutator role of activation-induced deaminase (AID) in B-cells during class-switch recombination and somatic hypermutation (5) targets preferentially cytosines within WRC (W: A|T; R: A|G) sequences (6), whereas apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) overexpression displays a preference for base substitutions at cytosines in TCW contexts (7). Other examples, such as the induction of C→T transitions at CG:CG dinucleotides by cytosine-5-methylation and the role of UV light in promoting base substitutions at pyrimidine dimers have been well documented (reviewed in (4,8)). More recently, complex patterns of base substitution at guanosines in cancer genomes have been found to correlate with changes in guanosine ionization potentials as a result of electronic interactions with flanking bases (9), suggesting a role for electron transfer and oxidation reactions in sequence-dependent mutagenesis. However, despite these advances, the increasing number of sequence-dependent patterns of mutation noted in genome-wide sequencing studies has met with a lack of understanding of most of the underlying mechanisms (10). Thus, a picture is emerging in which mutations are often heavily dependent on sequence-context, but for which our comprehension is limited.
Mononucleotide repeats comprise blocks of identical base pairs (A|T or C|G; hereafter referred to as A-tracts and G-tracts) and display distinct features: they are abundant in vertebrate genomes; mutations within the tracts occur more frequently than the genome-wide average; mutations generally increase with increasing tract length; length instability is a hallmark of mismatch repair-deficiency in cancers; and sequence polymorphism within the general population has been linked to phenotypic diversity (11–15). Thus, mononucleotide repeats appear ideal for addressing the question of sequence-dependent mutagenesis since base pairs within the tracts are flanked by identical neighbors. Both historic and recent investigations concur with the conclusion that a major source of mononucleotide repeat polymorphism is the occurrence of slippage (i.e. repeat misalignment) during semiconservative DNA replication, which gives rise to the addition or deletion of repeat units (11,12). An additional and equally important source of mutation has recently been suggested to arise from errors in DNA replication by translesion synthesis DNA polymerases, such as pol η and pol κ (13), also on slipped intermediates, leading to single base substitutions.
A key question that remains unanswered in these studies and which is relevant to the issue of sequence context-dependent mutagenesis is whether all base pairs within mononucleotide repeats display identical susceptibility to single base changes and whether indels (which are consequent to DNA breakage) occur randomly within the tracts.
Herein, we combine bioinformatics analyses on mononucleotide repeat variants from the 1000 Genomes Project and cancer genomes with molecular dynamics simulations and hybrid quantum mechanics/molecular mechanics calculations to address the question of sequence-dependent mutagenesis within these tracts. We show that mutations along both A-tracts and G-tracts are highly non-uniform. Specifically, both base substitutions and indels occur preferentially at the first and last bp of A-tracts, whereas they are concentrated between the second and third G:C base pairs in G-tracts. These positions coincide with the most flexible base pairs for A-tracts and with the preferential localization of a ‘hole’ that results when one electron is lost due to an oxidation reaction anywhere along G-tracts. Thus, despite the uniformity of sequence composition, mutations occur in a sequence-dependent context at homopolymeric runs according to a hierarchy that is imposed by both local DNA structural features and long-range base–base interactions. We also show that the repair processes leading to base substitution must differ between A- and G-tracts, since in the former, but not in the latter, base substitutions occur predominantly in the direction of the base immediately flanking the tracts. Additional sequence-dependent patterns of mutation are likely to arise from studies of more heterogeneous sequence combinations, possibly involving other aspects intrinsic to the structure of DNA.
RESULTS
Mononucleotide repeat variation is defined by tract length and flanking base composition
We define mononucleotide repeats in the GRCh37/hg19 (hg19) human genome assembly as uninterrupted runs of A:T and G:C base pairs (hereafter referred to as A-tracts and G-tracts, respectively) from 4 to 13 base pairs in length (Figure 1A). We retrieved a total of 48,767,945 A-tracts and 13,633,781 G-tracts, both of which displayed a biphasic distribution with an inflection point between tract lengths of 8 and 9 (bp) and with the number of runs declining with length more dramatically for G-tracts than for A-tracts (Figure 1B), as noted previously (29). Both the number of short tracts and the extent of decline varied with flanking base composition, TA[n]T runs being two- to three-fold more abundant than CA[n]Cs (Supplementary Figure S1A) and AG[n]As declining the most rapidly (Supplementary Figure S1B). Thus, mononucleotide runs exist as a collection of separate pools of sequences in extant human genomes, each maintained at distinctive rates of sequence stability, as determined by factors such as bp composition (A:T versus G:C), tract length and flanking sequence composition.
Mononucleotide repeat variation, evolutionary conservation and association with transcription. (A) The search algorithm was designed to retrieve runs of As or Ts (A-tracts) and Gs or Cs (G-tracts) length n (n = 4 to 13), along with their 5′ (n = 0) and 3′ (n = n + 1) nearest neighbors from hg19. Tract bases were numbered 5′ to 3′ with respect to the purine-rich sequence. The panel exemplifies the nomenclature for A- and G-tracts of length 4. (B) Logarithmic plot of the number of A-tracts (closed circles) and G-tracts (open circles) in hg19 as a function of length. (C) Normalized fractions of polymorphic tracts (F SNV) (number of SNVs divided by both hg19 number of tracts and n) from the 1KGP for A-tracts (closed circles) and G-tracts (open circles). (D) Radial plot of SNVs in the 1KGP at the 5′ and 3′ nearest neighbors of A-tracts. Periphery, tract length; horizontal axis, scale for the fraction of SNVs (F SNV). (E) Radial plot of SNVs in the 1KGP at the 5′ and 3′ nearest neighbors of G-tracts. (F) Percent difference in the numbers of A-tracts (closed circles) and G-tracts (open circles) between syntenic regions of hg19 and HN genomes. (G) The exponents of Benjamini-corrected P-values for A-tract-containing genes enriched in transcription-factor binding sites plotted as a function of A-tract length (triangles); each value represents the median of the top 11 USCS_TFBS terms. The percent A-tracts (closed circles) and G-tracts (open circles) intersecting genomic regions pulled-down by chromatin immunoprecipitation using antibodies against transcription factors are plotted as a function of tract length. (H) List of gene enrichment terms with a Benjamini-corrected P-value of <0.05 in common between genes containing A- and G-tracts of lengths 4–13, excluding the UCSC_TFBS terms.
We examined the extent of sequence variation in the human population by mapping 38,878,546 single nucleotide variants (SNVs) from 1092 haplotype-resolved genomes (the 1000 Genomes Project, 1KGP) (30) to the hg19 A- and G-tracts. The normalized fractions of polymorphic tracts (F SNV) were greater for G-tracts than A-tracts and both displayed Gaussian-type distributions, with maxima of 0.067 for G-tracts of length 8 and 0.017 for A-tracts of length 9 (Figure 1C). CA[n]C and AG[n]A runs displayed the highest F SNV values for A- and G-tracts, respectively (Supplementary Figure S1C and D), with F SNV values for AG[n]As attaining ∼0.10 at length 8. We conclude that flanking base composition influences the rates of SNV within mononucleotide runs and, as a consequence, their representation in the reference human genome.
F SNV values at the flanking 5′ and 3′ bp were similar between A- and G-tracts, except for minor differences for the least represented (i.e. longest) tracts and did not exceed 0.02 (Supplementary Figure S1E). These fractions are expected to be greater than at more distant positions from the tracts, based on previous data (29). SNVs at G-tracts, but not at A-tracts, were more frequent than at flanking base pairs. F SNVs for base pairs flanking short (≤8 bp) tracts were at least twice as high as those flanking long tracts; F SNVs also displayed distinct sequence preference with most (∼0.1) variants occurring at Ts 3′ of G-tracts (Figure 1D and E). In summary, SNVs at mononucleotide runs do not increase monotonically with length but peak at 8–9 bp. This behavior mirrors the genomic distributions, both with respect to the total number of tracts (Figure 1B) and the subsets flanked by specific-sequence combinations (Supplementary Figure S1A–D). Variation at flanking base pairs also displayed a biphasic pattern centered at a length of 8–9 bp, with a greater chance of variation adjacent to G- than A-tracts and with characteristic sequence preferences.
Long tracts are evolutionarily conserved and associated with high transcription
To assess whether more variable monosatellite runs (Figure 1C) might have undergone a greater reduction in number in extant humans relative to extinct hominids, we compared the number of A- and G-tracts between syntenic regions of five individuals comprising hg19 and three Neanderthal (HN) specimens (31). The difference between hg19 and HN was very small (<±2%) for the short tracts, but it displayed more negative values in hg19 with increasing tract length, which reached a maximum of −11.8 and −32.7% for A- and G-tracts, respectively, of length 9. Beyond this threshold, the numbers of tracts converged for A-tracts, whereas they were more abundant in hg19 for G-tracts >11 bp (Figure 1F). In summary, the largest difference in the number of mononucleotide runs between hg19 and HN sequences was centered at 9 bp for both A- and G-tracts, suggesting that the length distributions (Figure 1A and Supplementary Figure S1A and B) reflect distinct rates of evolutionary gains and losses due to differential sequence mutability (Figure 1C) as a function of length and flanking sequence composition (12).
The fact that long (>9 bp) mononucleotide runs display low variability in the human population (Figure 1C) and sequence conservation during evolutionary divergence (Figure 1F) raises the possibility that they might serve functional roles. Through gene enrichment analyses, we found that genes containing A- and G-tracts were enriched for genes associated with the term ‘UCSC_TFBS’, which pertains to transcripts harboring frequent transcription factor binding sites (32,33). For A-tract-containing genes, the median P-values for the top 11 UCSC_TFBS terms decreased from 2.95E-26 for tracts of length 4 to 5.22E-241 for tracts of length 13 (Figure 1G). The percent of A-tracts intersecting genomic fragments amplified from chromatin immunoprecipitation using transcription-factor binding antibodies (32,33) also increased from 8.7 to 9.9 from length 6 to 13, whereas it was constant (mean ± SD, 22.4 ± 1.1) for G-tracts (Figure1G). For gene classes excluding ‘UCSC_TFBS’, a search for categories enriched at P < 0.05 and common to all A- and G-tract-containing genes returned a set of 25 terms, 22 of which were associated with high levels of tissue-specific gene expression (Figure 1H). In summary, these analyses extend prior work (14) supporting a role for mononucleotide tracts in enhancing gene expression, a function that for A-tracts appears to increase with increasing tract length.
Repeat variability is highly skewed
Next we addressed whether bp along A- and G-tracts display equal probability and type of variation. In the 1KGP dataset, the number of SNVs at each position along both A- and G-tracts of length 4 was within a two-fold difference (144,000–240,000); for both types of sequence, transitions (i.e. A→G and G→A) were the predominant (51–78%) type of base substitution (Supplementary Figure S2A and B). However, with increasing length, the number of SNVs decreased up to 30-fold more drastically for G-tracts than for A-tracts, with increasing numbers of transversions (A→T and G→C|T) being predominant. Normalizing the data for the number of tracts genome-wide revealed that the extent of SNV varied by up to 10-fold, depending upon tract length and bp position. Specifically, the highest degree of variation was observed at the first and last A within the A-tracts (i.e. A1 and An), which underwent up to 61% A→T and 43% A→C transversions, respectively, at length 9 (Figure 2A). Likewise, for G-tracts, the most polymorphic sites were G3, followed by G2, for mid-size tracts of 8–10 bp, with 44% G→C transversions at G3 for tracts of length 8 (Figure2B). Thus, the extent of SNV at mononucleotide runs is grossly skewed in human genomes, both along the sequence itself and across tract length, which must account for the bell-shape behavior in F SNV for the tracts as a whole (Figure 1C).
Population variation spectra. (A) Variation spectra of A-tracts. Percent (number of SNVs at each position divided by the number of tracts in hg19 × 100) of A→T (black), A→C (red) and A→G (green) SNVs in the 1KGP dataset (left). Percent SNVs at A1 as a function of tract length (right). (B) Variation spectra of G-tracts. As in panel A with G→T (black), G→C (red) and G→A (cyan) (left). Percent SNVs at G3 as a function of tract length (right). (C) Percent A→T, A→C and A→G transitions at each position along A-tracts (stars) preceded and followed by a T (TA[n]T, left), C (CA[n]C), center) and G (GA[n]G, right) as a function of tract length. (D) Percent G→T, G→C and G→A transitions at each position along G-tracts (stars) preceded and followed by a T (TG[n]T, left), C (CG[n]C), center) and A (AG[n]A, right) as a function of tract length. (E) Percent transitions at base pairs (stars) preceding or following A-tracts (left) and G-tracts (right) as a function of tract length (n). *, mutated position.
We assessed whether SNV hypervariability was associated with specific combinations of nearest neighbors. For A-tracts flanked 5′ by a T, C or G, the highest percentage of SNVs was observed at A1 when preceded by a T, which reached 7.9% for TA[n] tracts of length 9 (Supplementary Figure S2C). By contrast, for 3′ T, C or G, the greatest effect was elicited by a C, with the highest percentage (7.1%) of SNVs at An for A[n]C tracts of length 9 (Supplementary Figure S2D). Therefore, flanking base pairs play a critical role both in the spectra and frequencies of SNVs at A-tracts. More detailed plots along A-tracts either preceded (Supplementary Figure S2E), followed (Supplementary Figure S2F) or preceded and followed (Figure 2C) by a T, C or G revealed the dramatic and long-range (up to 9–10 bp for the longest tracts, higher than the value of 4 bp predicted by mathematical models of slippage (11)) influence of flanking base pairs on variation spectra, in which up to 95% of the changes were in the direction of the base flanking the tract. Because the number of A-tracts preceded or followed by a specific base varies by up to three-fold (Supplementary Figure S2G), we conclude that for A-tracts, the overall mutation fractions and spectra are the result of at least three variables; length, position along the tract, and base composition of the 5′ and 3′ nearest-neighbors.
For G-tracts flanked 5′ by a T, C or A, high percentages (10–12%) of SNVs were observed at G1 for tracts preceded by a C, an effect that decreased with increasing tract length (Supplementary Figure S3A). This result, together with an exceedingly low number of G→A transitions at G1 for tracts not preceded by a C (Supplementary Figure S3C) relative to all tracts (Supplementary Figure S2B), is consistent with the known high mutability of CG:CG dinucleotides as a result of cytosine-5 methylation (9). The hypermutability at G2 was observed preferentially for tracts preceded by an A, and to a lesser extent T, whereas that at G3 was insensitive to flanking sequence composition. Likewise, G-tracts flanked 3′ by a T, C or A did not display marked sequence-dependent effects (Supplementary Figure S3B). Detailed plots of the SNV spectra along G-tracts either preceded (Supplementary Figure S3D), followed (Supplementary Figure S3E), or preceded and followed (Figure 2D) by a T, C or A revealed a noticeable effect only for 5′ T in association with G→T substitutions at G1for tracts of length ≥8. Thus, despite a consistent over-representation of G-tracts flanked 5′ by a T (Supplementary Figures S3F and S1B), which must account for the high absolute number of SNVs at G1 for TG[n] relative to AG[n] and CG[n] (Supplementary Figure S3G), nearest-neighbor base composition seems to play a lesser role in SNV spectra at G-tracts than at A-tracts.
With respect to SNVs at the flanking 5′ and 3′ nearest positions, no B→A or H→G substitutions (Figure 1A) were found above a length threshold of 9 for A-tracts and 8 for G-tracts (Figure 2E, gray shading) out of 5969 SNVs, implying that tract expansion by recruiting flanking base pairs is disfavored at these lengths. In summary, base substitution along mononucleotide repeats is strongly skewed towards the edges of A-tracts and within the 5′ half of G-tracts, with frequencies that peak at midsize lengths (8–9 bp). For A-tracts ≥7 bp, base substitution occurred almost exclusively in the direction of the flanking nearest-neighbors. Finally, base substitution at flanking bases did not contribute to tract expansion for mononucleotide runs longer than 8–9 bp.
Insertions and deletions display length and positional preference
In addition to SNVs, mononucleotide runs are polymorphic in length as a result of indels. Herein, we consider separately two types of indels: one in which tract length changes by ±1 and flanking bp composition is not altered (slippage); the other comprising all other cases involving the addition or removal of 1–200 bp (indels). Slippage is a widely accepted mutational mechanism (11–12,34), whereby DNA replication errors at reiterated DNA motifs cause changes in the number of motifs (most often +/−1). The normalized fractions of slippage in the 1KGP dataset peaked at lengths of 8 bp for A-tracts and 9 bp for G-tracts (Figure 3A), generating bell-shaped curves similar to those observed for SNVs (Figure1C) and with no differences in the highest fraction of ‘slipped’ tracts, which peaked at ∼0.02. By contrast, +1 slippage occurred more frequently than −1 slippage at A-tracts (Figure 3B). These results support recent studies on microsatellite repeats (12) and contrast with previous conclusions that slippage increases monotonically with tract length, and that the extent of slippage differs between A- and G-tracts (35,36).
Population insertions and deletions. (A) Normalized fractions of A-tracts (closed circles) and G-tracts (open circles) displaying +/−1 bp slippage in the 1KGP dataset as a function of tract length. Data were obtained by dividing the number of events by both the number of hg19 tracts and tract length (n). (B) Ratio of the number of +1 to −1 slippage for A-tracts (closed circles) and G-tracts (open circles). (C) Indels at A-tracts. For positions along the tracts (‘Tract’), ‘F Indel’ is the ratio between the number of indels and the number of tracts in hg19 multiplied by tract length. For the positions immediately flanking the tracts genomic coordinates (‘Before tract’ and ‘After tract’), ‘F Indel’ is the ratio between the number of indels and the number of tracts in hg19. (D) Indels at G-tracts, calculated as described in panel C. (E) Heatmap representation of insertions along A-tracts. The percent insertions (i.e. the number of insertions at each position divided by the number of tracts in hg19) (y-axis) plotted as a function of location (x-axis) from position 0 (insertion between the bp 5′ to the tract and the first bp of the tract) to position n + 1 (insertion between the bp 3′ to the last bp of the tract and the following bp) (see Figure 1A) and as a function of tract length (z-axis). (F) Heatmap representation of insertions along G-tracts.
With respect to indels, the normalized fractions were low (<1 × 10−3) along short (4–6 bp) A- and G-tracts, but rose to a plateau for longer tracts as reported earlier (11); this plateau was 10-fold higher for G-tracts (∼0.03) than for A-tracts (∼0.003) (Figure 3C and D). Indels also occurred more frequently (up to six-fold for A-tracts of length 11) at nearest-neighboring base pairs (‘Before tract’ and ‘After tract’ in Figure 3C and D) than along the tracts. Thus, contrary to SNVs and slippage, indels increased to a plateau with mononucleotide tract length.
We analyzed in detail the locations of insertions along the tracts and the flanking positions with respect to the 5′ to 3′ orientation of the tracts (Figure 1A). The normalized fractions demonstrated that insertions peaked at the 3′, and to a lesser extent 5′, ends of the longest A-tracts (Figure 3E), but remained low. For G-tracts, insertions occurred most efficiently at two locations (G2–3 and G5) (Figure 3F), they increased with tract length (up to ∼0.04), and attained ∼10-fold higher values than for A-tracts. In conclusion, insertion sites at A- and G-tracts followed the patterns observed for SNVs (Figure 2A and B), suggesting that factors associated with local DNA dynamics sensitize specific bases along the tracts to genetic alteration, inducing both SBS and indels.
Base pair flexibility and charge localization map to sites of sequence changes
To elucidate elements of intrinsic DNA dynamics that may be responsible for the biases in SNV and insertion sites, we performed molecular dynamics (MD) and hybrid quantum mechanics/molecular mechanics (QM/MM) simulations on model A[6], A[9], G[6] and G[9] duplex DNA fragments. We focused on water bridge coordination (Figure 4A), bp step flexibility, and for the G[6] and G[9], charge localization, as these properties are known to impact the susceptibility of DNA to base damage, repair and mutation. The fractions of one water coordination increased along the A[9] and A[6] structures in a 5′ to 3′ direction, irrespective of flanking sequence composition, in concert with a decrease in minor groove width (Figure 4B and Supplementary Figure S4A) as predicted (37). Vstep, a measure of bp structural fluctuation, displayed a prominent peak of ∼40 Å3deg3 at the 5′-TA-3′ step for both structures (Figure 4C and Supplementary Figure S4B), which together with low water occupancy points to 5′-TA-3′ being a preferred location for base modification and mutation. In the G[9] and G[6] structures water coordination involved mostly two-water bridges due to wide (∼14 Å) minor grooves (Figure 4Dand Supplementary Figure S4C), whereas flexibility was modest (∼20–22 Å3deg3, Figure 4E and Supplementary Figure S4D). Thus, bp dynamics are likely to impact mutations at A-tracts to a greater extent than at G-tracts. Guanine has the lowest ionization potential (IP) of all four bases and IP further decreases at guanine runs, rendering them targets for electron loss, charge localization, oxidation and eventually mutation (4,38). Because after electron loss the ensuing charge (hole) can migrate along the DNA double-helix and relocalize at specific guanines, we addressed whether the preferred sites of mutation along G-tracts, i.e. G2–3 and G5, would also be preferred sites for charge localization. The QM/MM determinations indicated that whereas for the short G[6] fragment the difference in the density-derived atomic partial charges (DDAPC) (i.e. the hole) localized most often (∼50%) to the first position (Figure 4F), for the long G[9] fragment charge localization shifted downstream (mostly to the second, but also to positions 6–7, Figure 4G). Importantly, the charge was found exclusively around the guanine rings (Figure 4H). Thus, the two main sites of sequence change along G-tracts, i.e. G2–3 and G5, coincide with positions where charge localization and hence one-electron oxidation reactions is predicted to occur most frequently. In summary, bp flexibility at A-tracts and charge transfer at G-tracts likely represent intrinsic DNA features underlying the bias in SNV and insertions at mononucleotide runs in human genomes.
MD and QM/MM simulations. (A) Molecular modeling of one (left) and two (right) minor groove water bridge coordination. (B) Fraction of one-water bridge occupancy (left axis) at A[9] DNA sequences flanked 5′ and 3′ by a T (black circles), C (red circles) or G (green circles). Minor groove widths (right axis), as determined from intrastrand phosphate-to-phosphate distances. (C) Vstep for A[9] DNA sequences, determined as the product of the square root of the eigenvalues (λi) described by the six bp step parameters shift, slide, rise, tilt, roll and twist; i.e. Vstep=∏6i=1λi−−√. (D) Fraction of one- (black circles) and two-water (red circles) bridge occupancy (left axis) at G[9] DNA sequences. Minor groove widths (right axis), as assessed from intrastrand phosphate-to-phosphate distances. (E) Vstep for G9 DNA sequences. (F) Average charge redistribution (open circles and right axis) for G[6] DNA structures upon vertical ionization, examined by calculating the difference on the density-derived atomic partial charges (DDAPC) for the neutral and negatively charged states. Histogram of the number of instances (left axis) in which the largest charge redistribution occurred at a specific position along the G[6] structures. (G) DDAPC for G[9] DNA structures (open circles and right axis) and histogram of the number of instances (left axis) in which the largest charge redistribution occurred at a specific position. (H) VMD rendering of a G[9] DNA structure displaying hole localization at G2. Capped base pairs were removed for clarity.
Position and orientation along nucleosome core particles modulate sequence variation
DNA wrapped around histones in nucleosomes is subject to local deformation (39), which may impact mutation. Thus, we analyzed the 1KGP SNVs at A- and G-tracts predicted to overlap with well-positioned nucleosome core particles (NCPs) (16). In hg19, the percentage of tracts that overlap with NCPs decreased moderately from ∼90% at length of 4 to 81% and 71% for A- and G-tracts of length 13, respectively (Figure 5A), suggesting that mononucleotide runs are not depleted in NCPs in human genomes as previously proposed (40). A-tracts of lengths 4–8 base pairs displayed distinctive peaks along the NCP surface in phase with the helical repeat of DNA (10.5 bp) and with minor grooves facing toward the inner protein core (lengths 4–5) (16) (Figure 5B and Supplementary Figure S5A). A-tracts of length of 9–13 bp exhibited only half (six) the peaks evident for the shorter tracts. For the G-tracts, only small peaks with no clear minor groove-inward-facing regions were detected (Supplementary Figure S5B).
Positioning along nucleosome core particles. (A) Percent of A-tract (open circles) and G-tract (closed circles) base pairs in hg19 overlapping with well-positioned NCP genomic coordinates as a function of tract length. (B) Counts of base pairs in hg19 A-tracts of length 5 overlapping with NCPs genomic regions as a function of distance from the histone octamer dyad axis. Minor groove-inward-facing regions (gray) were derived from the X-ray crystal structure of NCP147 (41). (C) Percent SNVs in the 1KGP dataset (left axis) at every bp along A-tracts of length 5 for tracts centered at maxima (black) and minima (gray) along NCPs (Figure 5B). Percent increase (right axis) of SNVs at minima relative to maxima (green). P-values for paired t-tests: 0.013 (*), 0.002 (**) and 4.7 × 10−6 (***). (D) Whisker plots of%SNVs (left axis) at A1 for A-tracts of length 5 centered at maxima and minima (black) along NCPs (Figure 5B). Percent difference (right axis) in the number of A-tracts of length 5 in hg19 preceded by C, T or G (red) between those centered at minima and those centered at maxima (Figure5B). (E) C-containing/G-containing ratios (see text) for G-tracts of length 5 in hg19 as a function of distance from the NCP dyad axis (black) and location of core histones (maroon and green). Peaks correspond to negative iSAT (i.e. tilt parameters multiplied by the corresponding sin θ) values (gray) (39). Ratios of%SNV at G1 (upshifted by 0.5 for clarity) between C-containing (5′-CCCCCG-3′ sequences on the hg19 forward strand) and G-containing (5′-CGGGGG-3′ sequences on the hg19 forward strand) (Figure 1A) CG[5] tracts mapping NCP Chip-seq genomic intervals (red) fitted by a non-parametric local regression (loess; sampling proportion, 0.100; polynomial degree, 3). (F) VMD rendering (top) of TATTT residues 34–38 (yellow) and the complementary AAATA residues 672–753 (pink) from the 1EQZ pdb nucleosomal crystal structure, corresponding to peak area from −40 to −36 in Figure 5E. The switch in G-tract (lengths of 5 and 7) orientation along NCPs (bottom) serves to position the C-containing strand on the outside (yellow) and, correspondingly, the G-containing strand on the inside (pink).
To assess if tract-positioning along NCPs influences SNVs, we selected A-tracts of lengths 5, 7 and 9 bp and G-tracts of lengths 5 and 7 bp whose central positions coincided with either the maxima or minima (41) (Figure 5B and Supplementary Figure S5A and B) and conducted pair-wiset-tests (330 total) between permutations of ‘categories’, including ‘tracts centered at maxima versus minima’, ‘position along the tracts’, ‘flanking sequence composition’, ‘specific NCP locations’ and ‘tract orientation’. For A-tracts, 79/207 (38%) significant pairs were found, 68 (86%) of which were related to differences between tracts centered at maxima versus minima, with a preponderance (63%) of tests displaying increased %SNVs at minima (Supplementary Figure S5C and E). For example, %SNVs at length 5 bp were greater at minima than at maxima at each position along the A-tracts (Figure 5C). A→C substitutions at A1 were more abundant at maxima than at minima (mean ± SD, 18.7 ± 0.7% at max and 17.6 ± 0.8% at min; P-value 0.001), whereas A→T substitutions at the same position displayed the opposite trend (mean ± SD, 18.4 ± 0.5% at max and 19.8 ± 1.1% at min; P-value 0.0005) (Figure 5D). A-tracts of length 7 also exhibited a similar pattern at A7 (Supplementary Figure S5H). The percentages of CA[5] and A[7]C tracts in hg19 centered at maxima were greater than at minima and the reverse was observed for the TA[5] and A[7T] tracts (Figure 5D and Supplementary Figure S5H). Thus, we conclude that positioning along the NCP surface of both the double-helical grooves and junctions with flanking base pairs influence SNVs along A-tracts. However, this influence is complex and for the most part, difficult to predict.
For G-tracts, most pairwise comparisons (18/34, 53%) indicated SNV variation according to sequence orientation (Supplementary Figure S5F and G). In hg19, the ratio of the numbers of G-tracts of lengths 5 and 7 for which the C-containing strand coincided with the forward sequence (downstream example sequence in Figure 1A) to the numbers of G-tracts for which the G-containing strand coincided with the forward sequence (upstream example sequence in Figure 1A) (C-containing/G-containing ratios) displayed a prominent 10.5-bp oscillation in phase with iSAT (Figure 5E), a measure of ‘inside’ and ‘outside’ bases, according to the bp step tilt parameter (39). Analysis of the helical path of a 146-bp DNA fragment wrapped around histones showed that the oscillation in the C-containing/G-containing ratios corresponds to a preference for guanine bases to face the protein core (Figure 5F). We analyzed the subset of G-tracts preceded by a 5′ C (i.e. CG[5]) to assess whether SNVs at G1, the position known to be mutable due to CpG methylation also oscillated with the C-containing/G-containing ratios. Oscillation in SNV-C-containing/SNV-G-containing values was evident, with peaks aligning to the hg19 troughs (Figure 5E) implying that the cytosines facing the protein surface harbor more variants than those facing away. We conclude that A- and G-tracts display preferential positioning (the former) and orientation (the latter) along NCPs, which in turn modulate the rate of sequence variation.
Mutations associated with human disease
Knowing that the first and last As of long A-tracts and G2–3 in G-tracts are the major sites of SNV in the human population, we addressed whether these features are also discernible in mutated mononucleotide tracts associated with human genetic disease. We collected 9,450,456 unique SBSs (both SBSs and SNVs refer to single base changes) from sequenced cancer genomes and normalized the percent mutations along A- and G-tracts to enable a direct comparison with the 1KGP dataset. For A-tracts (Figure 6A and Supplementary Figure S6A), SBSs displayed the same trend as the 1KGP data (Figure 2A) with respect to the bell-shape increase in mutations at A1 and An and the mutation spectra, although the susceptibility to mutation as a function of tract length attained greater values (6.36% for length 11 in cancer versus 4.15% for length 9 in the 1KGP datasets at A1). The first and last 3 bp also harbored more SBSs than in the 1KGP dataset for tracts >7 bp, a feature that we found to be due exclusively to a large cancer dataset (42) containing high-level microsatellite instability (MSI) samples (Supplementary Figure S6B and C), which are known to result from mismatch-repair deficiency (15). Thus, A-tracts display similar patterns of base substitution between the germline and somatic cancer tissues. For G-tracts, mutation spectra were characterized by G→T transversions at tract lengths >7, particularly at G1, the most frequently mutated position for tracts lengths up to 11 bp (Figure 6B and Supplementary Figure S6D). This trend persisted even when the high rates of methylation-mediated deamination mutations at the CG dinucleotide were removed (Supplementary Figure S6E). Thus, mutation patterns in cancer genomes contrast with those observed in the germline, both with respect to the most mutable position (G1 versus G2–3) and the types of base substitution (G→T in cancer genomes versus G→T and G→C in the germline).
Mutation patterns in cancer genomes. (A) Mutation spectra for SBSs at A-tracts. Percent values were obtained by dividing the total number of SBSs at each position by the number of tracts in hg19 and then multiplying by 3.2516 to equalize the percentage of A-tracts of length 4 between the cancer genomes and the 1KGP datasets. (B) Mutation spectra for SBSs at G-tracts in cancer genomes. Percent values were obtained as in (A) using a multiplication factor of 3.7419. (C) Normalized fractions of A-tracts (closed circles) and G-tracts (open circles) displaying +/−1 bp slippage, obtained by dividing the number of events by both the number of tracts in hg19 and tract length. (D) Indels at A-tracts, calculated as described in Figure 3C. (E) Indels at G-tracts, calculated as described in Figure3C. (F) Heatmap representation of insertions along G-tracts, as described in Figure 3E.
With respect to slippage, the fractions for A-tracts elicited an excess at lengths 9 and 10 bp relative to the 1KGP dataset, which was also due to the MSI-containing dataset. For G-tracts, the fractions peaked at length 8, as for the 1KGP dataset (Figures 3A and 6C), implying that the propensity to undergo slippage is indistinguishable between the germline and soma. Indels were also more abundant at flanking base pairs than along the tracts (Figure 6D and E), particularly for G-tracts of length >7, similar to the 1KGP dataset (Figure 3C and D). Detailed analyses of insertions revealed that both G1 and the preceding position were the most significant sites of mutation (F-values up to 0.08 at G1 for tracts of length 8) (Figure 6F). Thus, the 5′ end of long G-tracts is the most susceptible site for both SBSs and insertions in cancer genomes, in contrast to the germline where these occur within the runs, typically at G2–3.
We also extracted the mutated A- and G-tracts from the Human Gene Mutation Database (HGMD), a collection of >150,000 germline gene mutations associated with human inherited disease. A total of 1519 genes were mutated at A- or G-tracts out of a total of 3972 (38%); 3480 SBSs and 2866 slippage events were noted within these tracts, 85 and 46% of which were predicted to be disease-causing, respectively (Figure 7A and Supplementary Table S1). Ranking genes by the number of literature reports indicated that among the top 10 entries three were associated with cancer (BRCA1, BRCA2 and APC), two with hemophilia (F8 and F9), four with debilitating lesions of the skin (COL71A), muscle (DMD), lung (CFTR) and kidney (PKD1), with one causing hypercholesterolemia (LDLR) (Figure 7B). Thus, mutations within A- and G-tracts carry a high social burden by contributing to some of the most common human pathological conditions.
Mutation patterns in HGMD and model for sequence context-dependent changes. (A) Number of germline SBSs and slippage events (Slip.) at A- and G-tracts in HGMD. Gene alterations were classified as disease-causing mutation (DM), likely disease-causing mutation (DM?), disease-associated and putatively functional polymorphism (DFP), disease-associated polymorphism with additional supporting functional evidence (DP) and invitro/laboratory orinvivo functional polymorphism (FP). Codon changes (SIFT predictor) were classified as damaging (d), null (n), tolerated (t) and low-confidence prediction (l). (B) The 10 most commonly reported genes in HGMD with mutations at A- and G-tracts. Various mutated tracts were generally reported for the same gene in different reports. (C) Mutation spectra for SBSs at A- (left) and G-tracts (right) in HGMD. Percent values were obtained by dividing the total number of SBSs at each position by the number of tracts in hg19 exons. A|G→T (black), A|G→C (red), A→G (green), G→A (cyan). (D) Normalized fractions of A-tracts (closed circles) and G-tracts (open circles) displaying +/−1 bp slippage, obtained by dividing the total number of events by the number of tracts in hg19 exons and by tract length. (E) Model for sequence context-dependent changes at A-tracts (left) and G-tracts (right). *, site of base modification.
For both A- and G-tracts, SBSs occurred mostly at tract lengths of 4–7, with patterns more similar to those in the 1KGP than in the cancer datasets, both with respect to the location of the most mutable positions (first and last As and first/second Gs) and the types of base substitution (A→T and G→H) (Figure 7C and Supplementary Figure S6F). Likewise, slippage events peaked at tract lengths of 7–9 as observed in the 1KGP dataset (Figure 7D). In summary, the patterns of both SBSs and slippage in the HGMD dataset followed the trend observed in the 1KGP dataset, suggesting that germline variants at mononucleotide repeats leading to either population variation or human inherited disease may have arisen through similar mechanisms.
DISCUSSION
Why are specific A:T and G:C base pairs within A- and G-tracts more susceptible to sequence changes than their identical neighbors? For A-tracts, bp flexibility may play a role. Chemical damage to DNA, such as by hydroxyl radicals has been shown to be proportional to the geometrical solvent-accessible surface of the atomic groups, which increases with DNA flexibility (43). Along A-tracts flexibility is restricted, but it is high at both the 5′ and 3′ junctions. Thus, the fact that the highest rates of mutation coincide with the highest degree of flexibility at the 5′-TA-3′ bp step is consistent with the view that this position may be susceptible to DNA damage as a result of flexibility. Other sources of DNA dynamics are also likely to be relevant, such as sugar flexibility at the junctions, which increases with tract length (44). Chemical modification at these junctions may then lead to base substitution and indels, the latter as a result of strand breaks.
With respect to SNV mutation spectra, these were found mostly in the direction of flanking base composition above a length of 7–8 bp. We interpret this behavior in terms of DNA slippage along A-tracts when attempts are made during translesion synthesis (TLS) to bypass a damaged site (Figure 7Ei). Two scenarios may be considered to account for A→T transitions at A1. In the first, the last tract-template base would loop out into the polymerase active site permitting base-pairing and strand elongation (Figure 7Eii) using the tract-flanking base as a template (34,45–46). In the second (Figure 7Eiii), slippage would occur behind the polymerase, prompting extension past the newly created A*:T mispair generated by primer/template misalignment. Either pathway would yield a common intermediate (Figure 7Eiv) that contains the base complementary to the junction across from the damaged site upon slippage resolution (34). Following DNA synthesis (S) and/or repair (R) (Figure 7Ev and vi), this mispair will generate a base change that is always identical to the tract-flanking base.
For G-tracts, the high rates of G→T transversions at G1 in cancer genomes are also consistent with preferred chemical attack at this site due to high flexibility (Figure 7F top). Direct chemical attack at a guanine is known to result in stable products, such as 8-oxo-G and Fapy-G, both of which are known to yield G→T transversions (47–50). Thus, G1 may be the most susceptible site for such reactions for G-tracts of lengths ≥7 (Figure 7Fright), which in cancer genomes would become a mutation hotspot. In the germline, SNVs peaked inside G-tract base pairs, while mutational spectra were insensitive to flanking base composition; these events are inconsistent with a role for template misalignment and slippage as noted for A-tracts. Rather, the correspondence between hotspot mutations at G2–3 and G5 and the QM/MM simulations suggest a role for charge transfer. A large body of work during the past 20 years using computational, theoretical chemistry and biophysical techniques on short oligonucleotides, has shown that guanine is the most easily oxidizable base in DNA and that indeed a guanine radical cation can be generated through long-range hole transfer from an oxidant via one-electron oxidation mechanisms (51–55). GGG triplets were found to act as the most effective traps in hole transfer by both experimental and theoretical work (56–59), demonstrating that the resulting guanine radical cation (or its neutral deprotonated form) became rather delocalized, but it preferentially centered at the first and second G. These well-established patterns of chemical reactivity are consistent with our experimental observation of high mutation frequencies at G1 for short G-tracts and the results from QM/MM simulations on G6. For longer tracts, the downstream shift in mutation hotspots, i.e., G2–3 and G5, also correlate well with the charge localization predicted from QM/MM simulations, which explicitly included solvent effects and structural fluctuations. Thus, in conjunction with the constrained density functional theory (60), both the neutral and oxidized forms of a guanine nucleobase can be reliably constructed to infer the accurate determination of mutational patterns of mononucleotide repeats in human genomic DNA.
The compact organization of the sperm genome (61), and presumably low levels of oxidative stress in the germline, may enable guanine oxidization through one-electron oxidation reactions rather than by direct chemical attack, thereby favoring the formation of radical cations. A charge injected at G1 by electron loss would then migrate to neighboring guanines and localize at sites of low IP, such as G2 (Figure 7F left). Guanine radical cations are known to readily undergo further chemical modification leading to products such as 8-oxo-G, oxazolone, imidazolone, guanidinohydantoin, and spiroiminodyhydantoin (62) (M in Figure 7F), to yield G→T, G→C and G→A substitutions (4,63). Our model is in line with recent observations in which mutations at guanines within short G-runs (1–4 bp) correlate with sequence-dependent IPs at the target guanine in cancer genomes (9). Interestingly, these correlations were not observed in the germline (9). We interpret these composite observations as follows. The IP values for G-runs have been shown to decrease asymptotically with tract length, although the absolute values vary according to the methods and assumptions used (we obtained a value of 5.43 eV for both G[6] and G[9]) (64,65). We suggest that short G-runs with high IPs undergo one-electron oxidation reactions in the oxidative environment of cancer cells but would be refractory to such a mechanism in the germline (Figure 7Fright yellow and left white sectors). As length increases and IP values fall, G-runs would be attacked directly by oxidants abundant in tumor cells (Figure 7F orange sector), whereas oxidation will be limited to electron loss in the germline environment (Figure 7F left yellow sector).
These models (template misalignment for A-tracts and charge transfer for G-tracts) suggest a more complex scenario for mechanisms underlying mononucleotide repeat polymorphism in the human population than recently proposed (13), in which nucleotide misincorporation by error-prone polymerases is proposed as a primary source of mutations at both A- and G-tracts. As already stated, the directionality of SNVs toward tract-flanking bases in A-tracts and the hotspot mutations at G2–3, supports multiple and distinct mechanisms of base substitution at mononucleotide repeats.
Our analyses highlight additional information, including the lack of mutations in the direction of tract-base composition for base pairs flanking long tracts, the association with gene expression and the preference of guanines for the inner NCP surface, and extend prior observations (12) such as the bell-shape character of base substitution and slippage, whose mechanisms remain to be fully clarified. Finally, we document the contribution of mononucleotide mutagenesis to key aspects of human pathology beyond the well-established MSI instability in cancer (15), including hemophilia and tissue degeneration. Our collective work supports the conclusion that as the human genome undergoes evolutionary diversification and along the way suffers disease-associated mutations, oxidation reactions including charge transfer may play a prominent role.
Severe combined immunodeficiency diseases (SCIDs) are a group of primary immunodeficiency diseases characterized by a severe lack of T cells (or T cell dysfunction) caused by various gene abnormalities and accompanied by B cell dysfunction (WHO, 1992; Buckley et al., 1997). The incidence rates in infants were 1/75,000-1/10,0000 (WHO, 1992), but no morbidity statistics are available in China. The 2 genetic modes of SCID include X-linked recessive and autosomal recessive genetic inheritance. X-linked severe combined immunodeficiency (X-SCID) is the most common form, accounting for 50-60% of SCID cases (Noguchi et al., 1993). Immune system abnormalities in patients with X-SCID include T-B+NK-, in which T cells (CD3+) and natural killer (NK) cells (CD16+/CD56+) are absent or significantly reduced, and the number of B cells (CD19+) is normal or increased, causing reduced immunoglobulin production and class switching disorder (Buckley, 2004; Fischer et al., 2005). The IL- 2Rg gene mutation has been confirmed to be a major cause of X-SCID (Noguchi et al., 1993). In recent years, great progress has been made in understanding the pathogenesis of primary immunodeficiency disease and its application in clinical treatment, particularly regarding the development of critical care medicine and immune reconstruction technology. With timely control of infection and early bone marrow or stem cell transplantation, X-SCID patients can be treated, prolonging survival time. Therefore, early diagnosis of X-SCID is very important for patient treatment. Gene diagnosis has become a better early diagnosis or differential diagnosis method. In addition, familial X-SCID brings a great psychological burden to the relatives of patients. Ordinary chromosome analysis and immunological evaluation cannot be used for female carrier identification and fetal diagnosis, and gene diagnosis is the most effective method of carrier detection and prenatal diagnosis. In this study, we detected mutations in 2 families with X-SCID and identified 2 novel mutations, confirming the X-SCID pedigrees. Prenatal diagnosis was performed for the pregnant fetus in the mother of one of the probands based on gene diagnosis. Female individuals in this family were subjected to carrier detection.
IL2Rg gene mutation test Direct sequencing of 1-8 exons and the flanking region of the IL2Rg gene by PCR in family 1 showed that the 3rd exon of the proband contained the c.361-363delGAG heterozygous deletion mutation, which led to deletion of the 121st amino acid glutamate (p.E121del) in its coding product. There were no sequence variations in other coding regions or in the shear zone. The proband’s mother carried the same heterozygous mutation, while his father did not carry the mutation site (Figure 2a, b, c). This mutation was not observed in any cases of the control group, and this family was identified as an X-SCID family. The c.510-511insGAACT insertion heterozygous mutation was present in the 4th exon of the proband’s mother in family 2. This mutation was a 5-base repeat of GAACT, resulting in a change in amino acid 173 from tryptophan into a stop codon (p.W173X). While there were no sequence variations in other coding regions or in the shear zone, the patient’s father did not carry the mutation (see Figure 2d, e). We did not find this mutation in the healthy control group. We presumed that the 4th exon of the deceased child in family 2 contained the c.510-511insGAACT insertion mutation, leading to X-SCID symptoms, and thus we speculated that this family was an X-SCID pedigree. Prenatal diagnosis We verified the chorionic villus status of the fetus in family 1 using the PowerPlex 16 HS System kit. The results of prenatal diagnosis showed that the fetal tissue contained no maternal contamination and that this fetus was female. The results of prenatal diagnosis showed that there was no c.361-363delGAG (p.E121del) heterozygous mutation in the female fetus of family 1.
Figure 2. Sequencing graph of IL2Rg gene in 2 pedigrees with X-chain severe combined immunodeficiency. a.-c. Family 1. a. Normal control (rectangle indicates 3 edentulous bases of this patient). b. Proband carrying the c.361- 363delGAG (p.E121del) mutation (arrow indicates deletion of fragment connection sites). c. The proband’s mother contained a c.361-363delGAG (p.E121del) heterozygous mutation (arrow). d.-e. Family 2. d. The proband’s mother carried the c.510-511insGAACT (p.W173X) heterozygous mutation (arrow indicates that the reverse sequencing graph was positive). e. Normal control (rectangular box indicates 2 normal copies of GAACT (the mutation fragment was 3 copies). Carrier detection results For the c.361-363delGAG (p.E121del) site, the gene analysis results of the female individual in family 1 showed that I2 (proband’s grandmother) was a heterozygous carrier and that II3 (proband’s aunt) was a non-carrier and had no mutations.
IL-2 can combine with the IL-2 receptor (IL-2R) of the immune cell membrane. IL-2R is composed of 3 subunits, including the IL-2Ra chain (CD25), IL-2Rb chain (CD122), and IL- 2Rg chain (CD132). IL-2Rg functional units in common with IL-4, IL-7, IL-9, IL-15, IL-21, and other cytokine receptors, and these regions are referred to as the total chain (Li et al., 2000). The IL-2Rg chain can maintain the integrity of the IL-2R complex and is required for the internalization of the IL-2/IL-2R complex; it is also the link that contacts the cell membrane surface factor region and downstream cell signal transduction molecules. Therefore, the integrity of the IL-2Rg chain is vital for the immune function of an organism (Malka et al., 2008; Shi et al., 2009).
Mutations in the IL2Rg gene, which encodes IL-2Rg, were identified to be a major cause of X-SCID in 1993 (Noguchi et al., 1993). The IL2Rg gene is located on chromosome X q21.3-22, is 37.5 kb length, and contains 8 exons, which encode 369 IL-2Rg amino acids. The IL2Rg chain exhibits varying structural regions, such as the signal peptide [amino acids (AA) 1-22], extracellular domain (AA 23-262), transmembrane region (AA 263-283), and intracellular region (AA 284-369). The WSXWS motif is located in the extracellular region (AA 237-241), while Box 1 is located in the intracellular region (AA 286-294).
By the end of 2013, the Human Gene Mutation Database contained a total of 200 mutations in the IL2Rg gene (HGMD Professional 2013.4). The most common mutation types in the IL2Rg gene were the missense or nonsense mutations, which result from single base changes. A total of 100 missense or nonsense mutations have been identified, followed by insertion or deletion mutations in a total of 50 species. The 3rd most common type of mutations includes shear mutations in approximately 30 species. Eight exons contained mutations, and mutations in 3rd or 4th exons were the highest, accounting for a total mutation rate of 43% (86/200). According to the X-SCID gene database (IL2RGbase) (http://research.nhgri. nih.gov/scid/), the gene mutations in IL2Rg mainly occurred in the extracellular region of the IL2Rg chain (Fugmann et al., 1998). Zhang et al. (2013) reported that the IL2Rg gene mutations in 10 patients with X-SCID in China were located in the extracellular region. Two mutations reported in our study were also located in the extracellular region. The mutation of IL2Rg gene in family 1 was a codon mutation in the 3rd exon, resulting in a 3-base deletion. The c.361-363delGAG (p.E121del) mutation was located in the extracellular area of the IL- 2Rg subunit, and we inferred that the 121 glutamate deletion caused by the mutation would lead to changes in the structure of the peptide chain, affecting signal transmission and resulting in serious symptoms. The mutation of family 2 was a GAACT repeat of ILR2g gene; this repeat of 5 bases resulted in 173 codon changes from tryptophan into a stop codon. Generation of the peptide chain with the mutation lacked 196 amino acids compared to the normal chain, including the intracellular, transmembrane, and some extracellular regions, directly affecting the structure and function of receptors and causing disease. No studies have been reported regarding these 2 mutations. We combined with the mutation characteristics and clinical manifestations and diagnosed family 1 as X-SCID pedigrees. Although the patient in family 2 was deceased, it can be speculated that the 2 deceased patients in family 2 were X-SCID pedigrees caused by c.510-511insGAACT (W173X).
Prenatal diagnosis can accurately identify fetal situations and be used to avoid birth defects, which can also ease the anxiety of the pregnant mother. Gene diagnosis for pedigrees of patients based on DNA samples has advanced recently, particularly with the application of high-throughput sequencing technology (Alsina et al., 2013). We can now perform gene analysis for varied clinical infectious diseases for differential diagnosis. However, the effectiveness of prenatal diagnosis for pedigrees in which the proband is dead remains unclear. Because the gene mutations in the proband is unknown in these cases, the patient’s situation was only inferred by his mother’s genotypes. However, we considered that for the deceased, if we can define the mother was a pathogenic gene carrier, even if the proband is not X-SCID, the woman also has a risk of having X-SCID children and this pedigree may be X-linked recessive inheritance. Prenatal diagnosis may provide a choice for preventing the birth of patients in these families in the premise of informed consent.
Gene diagnosis of IL2Rg can also be used for carrier detection of suspected females in the family.
In the present study, we performed carrier detection of the patient’s grandmother and aunt in family 1 and determined that the patient’s pathogenic mutations were from his grandmother. His aunt did not inherit the pathogenic gene, and thus she was a non-carrier and her fertility will not be affected. In this study, we used direct sequencing of PCR products and identified IL2Rg gene mutations in 2 pedigrees with X-SCID. We found 2 unreported mutations in the IL2Rg gene, and prenatal diagnosis and carrier detection were conducted in 1 X-SCID family. Because the incidence rate of X-SCID is extremely low, it is difficult to promote the widespread use and application of genetic diagnosis. However, this study may provide some implications for the diagnosis of infants with immunodeficiency, and gene diagnosis techniques such as conventional or high-throughput sequencing should be used as soon as possible during pregnancy, which can be used to guide treatment. This method can also provide reliable prenatal diagnosis and carrier detection service for these families.
MEF2A gene mutations and susceptibility to coronary artery disease in the Chinese population
Coronary artery disease (CAD) has high morbidity and mortality rates worldwide. Thus, the pathogenesis of CAD has long been the focus of medical studies. Myocyte enhancer factor 2A (MEF2A) was first discovered as a CAD-related gene by Wang (2005) and Wang et al. (2003, 2005). Three mutation points in exon 7 of MEF2A were subsequently identified by Bhagavatula et al. (2004); however, Altshuler and Hirschhorn (2005) and Weng et al. (2005) predicted that the MEF2A gene lacked mutations. Zhou et al. (2006a,b) analyzed the mutations and polymorphisms in exons 7 and 11 of the MEF2A gene in the Han population in Beijing, and various rare mutations were found in exon 11 rather than in exon 7. The clinical significance of specific 21-bp deletions in MEF2A was also explored, and previous studies have shown mixed results. In this study, polymerase chain reaction-singlestrand conformation polymorphism (PCR-SSCP) and DNA sequencing were used to detect exon 11 of the MEF2A gene in samples collected from 210 CAD patients and 190 healthy controls and to investigate the function of the MEF2A gene in CAD pathogenesis and their correlation.
CAD, a common disease in China, is induced by multiple factors, such as genetics, the environment, and lifestyle. Thus, a multi-faceted approach is necessary in the study of CAD pathogenesis, particularly in molecular biology research, which is important for developing comprehensive treatment of CAD based on gene therapy. The MEF2A gene was first identified as a CAD-related gene through linkage analysis of a large family with CAD (9 of 13 patients developed MI) in 2003.
In this study, we found the following mutations: 1) codon 451G/T (147191) heterozygous or homozygous mutation; 2) loss of 1 (Q), 2 (QQ), 3 (QQP), 6 (425QQQQQQ430), and 7 (424QQQQQQQ430) amino acids (147108-147131); and 3) codon 435G/A (147143) heterozygous mutation. Among these mutations, the synonymous mutation at locus 147191 was confirmed by reference to the National Center for Biotechnology Information (NCBI) database to be a single nucleotide polymorphism, which was also demonstrated in our study by the extensive presence of this polymorphism in healthy controls. However, the heterozygous mutation at locus 147143 was only found in the genomes of CAD patients, and was therefore identified as a mutation.
Given that MEF2A is a CAD-related gene, the results of various studies are controversial among several countries. Weng et al. (2005) screened gene mutations in exon 11 of the MEF2A gene from 300 CAD patients and 1500 healthy controls. They hypothesized that the changes in 5-12 CAG repeats are genetic polymorphisms and that the 21-base deletion in exon 11 of the MEF2A gene did not induce autosomal dominant genetic CAD. Gonzalez et al. (2006) suggested that the CAG repeat polymorphism was independent of MI susceptibility in Spanish patients. Kajimoto et al. (2005) reported that the CAG repeat sequence was not correlated with MI susceptibility in Japanese patients. Horan et al. (2006) also found that the CAG repeat sequence was not associated with the susceptibility to early-onset familial CAD in an Irish population. Hsu et al. (2010) identified no correlation between the CAG repeat sequence and CAD susceptibility in the Taiwanese population. Dai et al. (2010) found that the structural change in exon 11 was not related to CAD in the Chinese Han population. Lieb et al. (2008) and Guella et al. (2009) hypothesized that MEF2A was independent of CAD. However, Yuan et al. (2006) and Han et al. (2007) suggested that the CAG repeat sequence was correlated with CAD because 9 CAG repeats was an independent predictor of CAD. Elhawari et al. (2010) and Maiolino et al. (2011) suggested that MEF2A is a susceptibility gene for CAD. Dai et al. (2013) showed that mutations in exon 12 are associated with the early onset of CAD in the Chinese population. Liu et al. (2012) failed to demonstrate a correlation between the CAG repeat sequence and CAD through case-control analysis, systematic review, and meta-analysis, but found that the 21- base deletion in exon 11 was strongly associated with CAD, and that genetic variations in MEF2A may be a relatively rare, but specific, pathogenic gene for CAD/MI. Kajimoto et al. (2005) reported 4-15 CAG repeats. However, only 4-11 CAG repeats were observed in our study, possibly because of genetic differences in patients in this study. Eleven CAG repeats were observed in most samples from the control group, and the proportion of 10, 9, and 8 repeats exceeded 1%. The heterozygous mutation at 147143, as well as the 4 and 5 CAG repeats, was only observed in CAD patients. Thus, we speculated that the CAG repeat sequence is correlated with CAD susceptibility, and the presence of 4 or 5 repeats may be a risk factor for CAD, which was inconsistent with the results obtained by Han et al. (2007). The inconsistency in these results may be explained by the differences in subjects and sample sizes among studies.
Impact of glucocerebrosidase mutations on motor and nonmotor complications in Parkinson’s disease
Here, we conducted a multicenter retrospective cohort analysis, and the data were investigated by survival time analysis to show the impact of GBA mutations on PD clinical course. We also investigated regional cerebral blood flow (rCBF) and cardiac sympathetic nerve degeneration of subjects with GBA mutations, compared with matched PD controls.
3.1. Subjects
Among the 224 eligible PD patients (the subjects were not related to each other), 9 subjects were excluded from the analysis (4 due to multiple system atrophy findings on subsequent brain MRI and 5 because of insufficient clinical information). Therefore, 215 PD patients [female, 52.1%; age, 66.7 ± 10.8 (mean ± standard deviation)] were analyzed. For non-PD healthy controls, 126 patients’ spouses (female, 58.7%; age, 67.3 ± 10.3) without a family history of PD or GD were enrolled.
3.2. GBA mutations and risk ratios for PD
In the PD subjects, we identified 10 nonsynonymous and 2 synonymous GBA variants. Within the nonsynonymous variants, 7 mutations were previously reported in GD [R120W, L444P-A456P-V460 (RecNciI), L444P, D409H, A384D, D380N, and444L(1447-1466 del 20, insTG)] as GD-associated mutations. Three nonsynonymous mutations have never been reported in GD patients [I(-20)V, I489V, and there was one novel mutation (Y11H)].
GD-associated GBA mutations were found in 19 of the 215 (8.8%) PD patients but none in the healthy controls. The risk of PD development relative to these GD-associated mutations was estimated as an OR of 25.1 [95% confidence interval (CI), 1.50–420,p = 0.0001] with 0-cell correction. The nonsynonymous mutations that were not reported in GD patients had no association with PD development (p = 0.506; OR, 1.3; 95% CI, 0.7–2.6) ( Table 1). Four subjects had double mutations. For subsequent analyses, 2 subjects with double mutations of I (-20)V and K466K were adopted to the group of mutations unreported in GD, and 2 subjects with double mutations of R120W and I(-20)V, and of R120W and L336L were adopted to the group of GD-associated mutations.
Table 1.Frequency of glucocerebrosidase gene allele in Parkinson’s disease patients and controls
3.3. Clinical features of PD patients by GBA mutation groups
The clinical features of PD patients with GD-associated mutations, those with mutations unreported in GD, and those without mutations are shown in Table 2. In the GD-associated mutation group, females, those with a family history and those with dementia (DSM IV) were significantly more frequent than those in the no-mutation group (p = 0.047, 0.012, and 0.020, respectively). The age of PD onset was lower in patients with GD-associated mutations (55.2 ± 9.9 years ± standard deviation), compared with those without mutations (59.3 ± 11.5), although the statistical difference was not significant. There were no differences in clinical manifestations between subjects with mutations unreported in GD and those without mutations, except for dopamine agonist dosage (p = 0.026) ( Table 2).
Table 2.Epidemiological and clinical features of PD patients with Gaucher disease–associated GBA mutations, those with mutations previously unreported in GD and those without mutations
3.4. Survival time analyses to develop dementia, psychosis, dyskinesia, and wearing-off
Time to develop clinical outcomes (dementia, psychosis, dyskinesia, and wearing-off) was compared in 19 subjects with GD-associated mutations, 29 with mutations unreported in GD, and 167 without mutation. The median observation time was 6.0 years. The subjects with GD-associated mutations showed a significantly earlier development of dementia and psychosis, compared with subjects without mutation (p < 0.001 and p = 0.017) ( Supplementary Table e-1, Fig. 1A and B). We rereviewed the clinical record of the subject who showed early dementia (defined by DSM IV) ( Fig. 1A) and made sure it did not satisfy the criteria of DLB ( McKeith et al., 2005).
Fig. 1.
Kaplan–Meier curves of dementia and psychosis in Parkinson’s disease (PD) patients with Gaucher disease (GD)-associated glucocerebrosidase gene (GBA) mutations and those without mutations. PD patients with GD-associated GBA mutations and those without GBA mutations were compared to investigate the time taken to develop dementia (A) and psychosis (B). Because of insufficient information in several patients, the numbers in each analysis were different. The patients with and without mutations were 17 and 165 (A), 18 and 165 (B) against a total of 19 and 167. DSM IV, Diagnostic and Statistical Manual of Mental Disorders, revised fourth edition. p-Values were calculated by log-rank tests.
The associations of GBA mutations and these symptoms were estimated as HRs, adjusting for sex and age at PD onset. HRs were 8.3 for dementia (95% CI, 3.3–20.9; p < 0.001) and 3.1 for psychosis (95% CI, 1.5–6.4; p = 0.002). The time until development of wearing-off and dyskinesia complications was not statistically significant, with HRs of 1.5 (95% CI, 0.8–3.1; p = 0.219) and 1.9 (95% CI, 0.9–4.1; p = 0.086) ( Table 3).
Table 3.Hazard ratios of GBA pathogenic mutations for clinical symptoms
Model
Clinical feature
Hazard ratio
95% CI
p
1
Dementia (DSM-IV)
8.3
3.3–20.9
<0.001
2
Psychosis
3.1
1.5–6.4
0.002
3
Wearing-off
1.5
0.8–3.1
0.219
4
Dyskinesia
1.9
0.9–4.1
0.086
Each model was adjusted for sex and age at onset.
Key: CI, confidence interval; DSM-IV; The Diagnostic and Statistical Manual of Mental Disorders part 1IV; GBA, glucocerebrosidase.
Subjects with mutations unreported in GD did not show significant differences in time to develop all 4 outcomes, compared with no mutation subjects. Therefore, subjects with GD-unreported mutations were regarded as subjects without GBA mutations in further analyses.
3.5. rCBF on SPECT in patients with GD-associated GBA mutations
We conducted pixel-by-pixel comparisons of rCBF on SPECT between PD subjects with mutations (cases) and sex-, age-, and disease duration-matched PD subjects without any mutations in GBA (controls). Four controls were adopted for each case (except for a 34-year-old female case who was matched to a control), and in total 12 cases (female 50%, age at SPECT mean ± standard error (SE); 58.9 ± 3.3 years, disease duration at SPECT 7.3 ± 1.5 years) and 45 controls (female 64.4%, age at SPECT mean ± SE; 61.0 ± 1.3 years, disease duration at SPECT 7.1 ± 0.7 years) were analyzed. As a result, a significantly lower rCBF was seen in the cases compared to the controls in the bilateral parietal cortex, including the precuneus ( Fig. 2).
Fig. 2.
Regional cerebral blood flow in the group with GD-associated mutations compared with the matched Parkinson’s disease group without mutations. Regions with lower regional cerebral blood flow in the group with GD-associated mutations displayed on an anatomic reference map. Abbreviation: GD, Gaucher disease.
3.6. H/M ratios on MIBG scintigraphy in patients with GD-associated GBA mutations
Cardiac MIBG scintigraphy visualizes catecholaminergic terminals in vivo that are reduced as well as brain dopaminergic neurons in PD patients. We also investigated MIBG scintigraphy between 16 cases (female 68.8%, age at examination mean ± SE; 60.2 ± 2.6 years, disease duration at examination 6.2 ± 1.2 years) and sex-, age- and disease duration-matched 61 controls [(63.8 %, age 62.0 ± 1.1 years, disease duration 5.5 ± 0.6 years) (1:4 except for 1 young 34-year-old female case who was matched to a control)]. In the results, both early and late H/M ratios declined in both groups and did not show any significant differences (p = 0.309 and 0.244) ( Supplementary Table e-2).
4. Discussion
4.1. Contributions of GD-associated GBA mutations to the development of PD
In the analysis of 215 PD patients and 126 non-PD controls, we identified 10 nonsynonymous heterozygous GBA mutations, including 1 novel mutation. Among these mutations, 7 were GD-associated, and the patients carrying these mutations represented 8.8% of the PD cohort. No significant association was found between the GD-unreported mutations and PD development, which suggests that only the GD-associated mutations are a genetic risk for PD. According to a worldwide multicenter analysis of 1883 fully sequenced PD patients, 7% of the GD-associated mutations are found in non-Ashkenazi Jewish PD patients ( Sidransky et al., 2009). Although the mutation frequency in the present study was similar to previous results, the OR of GD-associated heterozygous mutations (25.1) was significantly greater than the OR (5.43) of other ethnic cohorts (Sidransky et al., 2009) and was consistent with an OR of 28.0 from a previous Japanese report ( Mitsui et al., 2009). These results, taken together, suggest the possibility thatGBA mutations are at a distinct risk for PD in the Japanese population. However, a larger Japanese cohort study is required to confirm this.
4.2. Cross-sectional clinical figures of PD with GBA mutations
4.3. Impact of GBA mutations on the clinical course of PD
To investigate the impact of GBA mutations on the clinical course of PD, a prospective-designed study over a long period is preferred. Although there has been a few longitudinally designed study to date, follow-up clinical data for a median of 6 years of 121 PD cases from a community-based incident cohort was recently reanalyzed; results demonstrate that progression to dementia defined by DSM IV (HR 5.7) and Hoehn and Yahr stage 3 (HR 3.2) are significantly earlier in 4 GBA mutation-carrier patients compared with 117 patients with wild-type GBA ( Winder-Rhodes et al., 2013). A 2-year follow-up clinical report of 28 heterozygous GBA carriers who were recruited from relatives of GD-patients shows slight but significant deterioration of cognition and smelling, compared to healthy controls ( Beavan et al., 2015). Brockmann et al. (2015)assessed motor and nonmotor symptoms including cognitive and mood disturbances for 3 years in 20 PD patients with GBA mutations and showed a more rapid disease progression of motor impairment and cognitive decline in GBA mutation cases comparing to sporadic PD controls. The current long-term retrospective cohort study up to 12 years reinforced these results. It revealed that dementia and psychosis developed significantly earlier in subjects with GD-associated mutations compared with those without mutation, and the HRs of GBA mutations were estimated at 8.3 for dementia and 23.1 for psychosis, with adjustments for sex and PD onset age. In contrast, the results showed no significant difference in developing wearing-off and dyskinesia.
In this study, we also investigated whether GD-unreported mutations affected the clinical course of PD. In both cross-sectional and survival time analyses, the mutations unreported in GD carried no increased burden on clinical symptoms such as dementia, psychosis, wearing-off, and dyskinesia.
4.4. Reduced rCBF in PD with GBA mutations compared with matched PD controls
We found a significantly decreased rCBF, reflecting decreased synaptic activity, in the bilateral parietal cortex including the precuneus, in subjects with GD-associated mutations compared with matched subjects without mutations. The pattern of reduced rCBF was very similar to the pattern of H215O positron-emission tomography that Goker-Alpan et.al. (2012) reported, showing decreased resting rCBF in the lateral parietal association cortex and the precuneus bilaterally in GD subjects with parkinsonism (7 subjects with homozygous or compound heterozygous GBA mutations), compared with 11 PD without GBA mutations. Results suggest that PD with heterozygous GBAmutations and GD patients presenting parkinsonism had a common reduced pattern of rCBF. Interestingly, in their study, rCBF in the precuneus—but not in the lateral parietal cortex—correlated with IQ, suggesting that the involvement of the precuneus is critical for defining GBA-associated patterns.
4.5. Reduced cardiac MIBG H/M ratios as well as matched PD controls
We also showed that cardiac MIBG H/M ratios in subjects with GD-associated mutations were lower than the cutoff point for PD discrimination (Sawada et al., 2009), suggesting that postganglionic sympathetic nerve terminals to the epicardium were denervated, as well as in PD without mutations.
4.6. Mechanisms of impact on PD clinical course by GD-associated GBA mutations
Experimental studies suggesting a bidirectional pathogenic loop between α-synuclein and glucocerebrosidase have been accumulated (Fishbein et al., 2014, Gegg et al., 2012, Mazzulli et al., 2011, Noelker et al., 2015, Schondorf et al., 2014 and Uemura et al., 2015). Loss of glucocerebrosidase function compromises α-synuclein degradation in lysosome, whereas aggregated α-synuclein inhibits normal lysosomal function of glucocerebrosidase. The pathogenic loop may facilitate neurodegeneration in GD-associated PD brain, resulting in early development of dementia or psychosis as shown in the present study. Several recent researches propose the possibility that the similar mechanism as in PD with GBA mutations exists even in idiopathic PD brain ( Alcalay et al., 2015, Chiasserini et al., 2015, Gegg et al., 2012 and Murphy et al., 2014). On the other hand, the impacts of GD-associated GBA mutations for the development of motor complications such as wearing-off and dyskinesia were not statistically significant, suggesting other pathophysiological mechanisms in the striatal circuit brought out after long-term therapy especially by l-dopa.
4.7. Limitations
Our study has several limitations. In the design of the study, we assumed that the sample size was 215 (PD patients) for survival time analyses and investigated 224 PD patients. We assumed that the mutation prevalence would be 9.4%, and in fact, we found 19 patients with mutations (8.5%) of the 224 patients. Based on these figures, we estimated the risk ratios of heterozygous GBA mutations for the risk of PD development and PD clinical symptoms as ORs in the cross-sectional multivariate analyses, although the 95% CIs were broad. More of subject numbers will be needed to determine robust risk ratios.
Comprehensive Genetic Characterization of a Spanish Brugada Syndrome Cohort
Brugada syndrome (BrS) was identified as a new clinical entity in 1992 [1]. Six years later, the first genetic basis for the disease was identified, with the discovery of genetic variations inSCN5A [2]. Nowadays, more than 300 pathogenic variations in this first gene are known to be associated with BrS [3]. SCN5A encodes for the α subunit of the cardiac voltage-dependent sodium channel (Nav1.5), which is responsible for inward sodium current (INa), and thus plays an essential role in phase 0 of the cardiac action potential (AP). Genetic variations in this gene can explain around 20–25% of BrS cases [3].
Since BrS was classified as a genetic disease, several other genes have been described to confer BrS-susceptibility [4–7]. Pathogenic variations have been mainly described in: 1) genes encoding proteins that modulate Nav1.5 function, and 2) other calcium and potassium channels and their regulatory subunits. All these proteins participate, either directly or indirectly, in the development of the cardiac AP. Although the incidence of pathogenic variations in these BrS-associated genes is low [6], it is considered that, among all of them, they could provide a genetic diagnosis for up to an extra 5–10% of BrS cases. Hence, altogether, a genetic diagnosis can be achieved approximately in 35% of clinically diagnosed BrS patients.
Other types of genetic abnormalities have been suggested to explain the remaining percentage of undiagnosed patients. Indeed, multiplex ligation-dependent probe amplification (MLPA) has allowed the detection of large-scale gene rearrangements involving one or several exons ofSCN5A in BrS cases. However, the low proportion of BrS patients carrying large genetic imbalances identified to date suggests that this type of rearrangements will provide a genetic diagnosis for a modest percentage of BrS cases [8–10].
BrS has been associated with an increased risk of sudden cardiac death (SCD), despite the reported variability in disease penetrance and expressivity [11]. The prevalence of BrS is estimated at about 1.34 cases per 100 000 individuals per year, with a higher incidence in Asia than in the United States and Europe [12]. However, the dynamic nature of the typical electrocardiogram (ECG) and the fact that it is often concealed, hinder the diagnosis of BrS. Therefore, an exhaustive genetic testing and subsequent family screening may prove to be crucial in identifying silent carriers. A large percentage of these pathogenic variation carriers are clinically asymptomatic, and may be at risk of SCD, which is, sometimes, the first manifestation of the disease [13].
In the present work, we aimed to determine the spectrum and prevalence of genetic variations in BrS-susceptibility genes in a Spanish cohort diagnosed with BrS, and to identify variation carriers among relatives, which would enable the adoption of preventive measures to avoid SCD in their families.
Table 1. Demographics of the 55 Spanish BrS patients included in the study.
The table shows the demographic characteristics of all the patients included in the study. Numbers in parentheses represent the relative percentages for each condition. T1 ECG refers to Type 1 BrS diagnostic electrocardiogram (ECG), obtained either spontaneously, or after drug challenge. The information regarding both the electrophysiological studies (EPS) and the treatment was not available for all the patients. Two of the patients that didn’t receive any treatment died, and were not taken into account for the calculations of percentages (+2 dead). ICD, intracardiac cardioverter defibrillator.
Table 2. Characteristics of the Spanish BrS patients carrying rare genetic variations.
The table shows the clinical characteristics of the probands who carried rare genetic variations in SCN5A, SCN2B, or RANGRF. All of them are potentially pathogenic except that found in RANGRF, which is of unknown significance (see discussion). All the potentially pathogenic variations (PPVs) that had been previously reported, except p.P1725L and p.R1898C, had been identified in BrS patients. p.P1725L had been associated with Long QT Syndrome and p.R1898C was found in Exome Variant Server with a MAF of 0.0079%. No rare variations were identified in the control population. Patient’s age is expressed in years. Bold identifies the patients carrying variations that had not been described previously. M, male; F, female; S, syncope; ICD, intracardiac cardioverter defibrillator; UK, unknown; EPS, electrophysiological studies (+, positive response;-, negative response; N/P, not performed). The two patients who carried two PPVs each are identified by a and b, respectively.
We performed a genetic screening of 14 genes (SCN5A, CACNA1C, CACNB2, GPD1L,SCN1B, SCN2B, SCN3B, SCN4B, KCNE3, RANGRF, HCN4, KCNJ8, KCND3, and KCNE1L), which allowed the identification of 61 genetic variations in our cohort. Of these, 20 were classified as potentially pathogenic variations (PPVs), one variation of unknown significance, and 40 common or synonymous variants considered benign.
The 20 PPVs were found in 18 of the 55 patients (32.7% of the patients, 83.3% males; Table 2). Sixteen patients (88.9%) carried one PPV, and two patients (11.1%) carried two different PPVs each. Nineteen out of the 20 PPVs identified were localized in SCN5A and one in SCN2B.
The vast majority of the PPVs identified were missense (70%). We also detected 2 nonsense variations (10%), 3 insertions or deletions causing frameshifts (15%), and one splicing variation (5%). The three frameshifts (p.R569Pfs*151, p.E625Rfs*95 and p.R1623Efs*7) were identified in SCN5A. These were not found in any of the databases consulted (see Methods), and were thus considered potentially pathogenic (see below). The other 16 rare variations identified inSCN5A had been previously described, and hence were also considered potentially pathogenic. Fourteen of them had been identified in BrS patients. Of these, 6 had also been identified in individuals diagnosed with other cardiac electric diseases (i.e. Sick Sinus Syndrome, Long QT Syndrome, Sudden Unexplained Nocturnal Death Syndrome or Idiopathic Ventricular Fibrillation [2,15,16,20,21,25]). The other 2, p.P1725L and p.R1898C, had only been associated with Long QT Syndrome or found in Exome Variant Server with a MAF of 0.0079%, respectively. Furthermore, we identified a variation in SCN2B (c.632A>G in exon 4 of the gene, resulting in p.D211G) which was considered pathogenic. This patient was included within our cohort, but the functional characterization of channels expressing SCN2B p.D211G was object of a previous study from our group [7]. We also identified a nonsense variation in RANGRFwhich has been formerly reported as rare genetic variation of unknown significance [29].
Additionally, we screened the relatives of those probands carrying a PPV. We analysed a total of 129 relatives, 69 of which (53.5%) were variation carriers. Genotype-phenotype correlations evidenced that 8 of the families displayed complete penetrance (S3 Table). Additionally, no relatives were available for one of the probands carrying a PPV, thus hampering genotype-phenotype correlation assessment. The other 12 families showed incomplete penetrance.
MLPA analysis
The 37 patients with negative results after the genetic screening of the 14 BrS-associated genes underwent MLPA analyses of SCN5A. This technique did not reveal any large exon deletion or duplication in this gene for any of the patients.
SCN5A p.R569Pfs*151 (c.1705dupC), a novel PPV
A 41-year-old asymptomatic male presented a type 3 BrS ECG which was suggestive of BrS. Flecainide challenge unmasked a type 1 BrS ECG (Fig 1A, left), which was also spontaneously observed sometimes during medical follow up. Sequencing of SCN5A revealed a duplication of a cytosine at position 1705 (c.1705dupC; Fig 1A, right), which originated a frameshift that lead to a truncated Nav1.5 channel (p.R569Pfs*151). The proband’s sister also carried this duplication, but had never presented signs of arrhythmogenesis. The proband’s twin daughters were also variation carriers, displayed normal ECGs and, to date, are asymptomatic (Fig 1A, middle). Thus, p.R569Pfs*151 represents a novel genetic alteration in the Nav1.5 channel that could potentially lead to BrS, but with incomplete penetrance.
Fig 1. Characteristics of the probands carrying non-reported potentially pathogenic variations (PPVs) in SCN5A and their families.
Left: Electrocardiograms of the probands: (A) patient carrying the p.R569Pfs*151 variation, showing the ST elevation characteristic of BrS in V1 at the time of the flecainide test; (B) patient carrying the p.E625Rfs*95 variation, showing the spontaneous ST elevation characteristic of BrS in V1 and V2; and (C) patient carrying the p.R1623Efs*7 variation, showing the spontaneous ST elevation characteristic of BrS in V1 and V2. Middle: Family pedigrees. Open symbols designate clinically normal subjects, filled symbols mark clinically affected individuals and question marks identify subjects without an available clinical diagnosis. Plus signs indicate the carriers of the PPVs and minus signs, non-carriers. The crosses mark deceased individuals and arrows identify the proband. Right: Detail of the electropherograms obtained after SCN5Asequence analysis of a control subject (left panels) and of the probands (right panels).
A 51-year-old asymptomatic male was diagnosed with BrS since he presented a spontaneous ST segment elevation in leads V1 and V2 characteristic of type 1 BrS ECG (Fig 1B, left). The sequencing of SCN5A evidenced an adenine duplication at position 1872 (c.1872dupA, Fig 1B, right). This genetic variation results in a truncated Nav1.5 channel (p.E625Rfs*95). The genetic analysis of the proband’s relatives proved that only her mother carried the variation (Fig 1B, middle). She was asymptomatic, but a BrS ECG was unmasked upon ajmaline challenge. The proband’s sister was found dead in her crib at 6 months of age, which suggests that her death might be compatible with BrS. Therefore, the p.E625Rfs*95 variation in the Nav1.5 channel represents a novel genetic alteration potentially causing BrS.
SCN5A p.R1623Efs*7 (c.4867delC), a novel PPV
The proband, a 31-year-old male, was admitted to hospital after suffering a syncope. His baseline 12-lead ECG showed a ST segment elevation in leads V1 and V2 that strongly suggested BrS type 1 (Fig 1C, left). A deletion of the cytosine at position 4867 (c.4867delC) was observed upon SCN5A sequencing (Fig 1C, right). This base deletion leads to a frameshift that originates a truncated Nav1.5 channel (p.R1623Efs*7). Genetic screening of his parents and sisters evidenced that none of them carried this novel variation (Fig 1C, middle). None of them had presented any signs of arrhythmogenicity, nor had a BrS ECG. Nevertheless, in uterogenetic analysis of one of his daughters proved that she had inherited the variation. She died when she was 1 year of age of non-arrhythmogenic causes. Hence, the p.R1623Efs*7 variation in the Nav1.5 channel is a novel genetic alteration originated de novo in the proband that could potentially lead to BrS.
Synonymous and common genetic variations portrayal
In our cohort, we identified 40 single nucleotide variations which were common genetic variants and/or synonymous variants (S2 Table). Twenty-nine had a minor allele frequency (MAF) over 1%, and were thus considered common genetic variants.
We also identified 11 variants with MAF less than 1%. Of them, 9 were synonymous variants, what made us assume that they were not disease-causing. Four of these synonymous variants were not found in any of the databases consulted, and thus their MAF was considered to be less than 1%. Each of these synonymous variations was identified in 1 patient of the cohort. A similar proportion of individuals carrying these novel variations was detected upon sequencing of 300 healthy Spanish individuals (600 alleles). The remaining 2 variants were missense, and although they had either a MAF of less than 1% or an unknown MAF according to the Exome Variant Server and dbSNP websites, they were common in our cohort (29.2 and 50%, respectively; S2 Table), and a similar MAF was detected in a Spanish cohort of healthy individuals (26.7% and 48.8%, respectively).
Influence of phenotype and age on PPV discovery
To assess if a connection existed between the probands’ phenotype and the PPV detection yield, we classified the patients in our cohort according to their ECG (spontaneous or induced type 1), the presence of BrS cases within their families, and the presence/absence of symptoms. Even though the overall PPV detection yield was 32.7%, it was even higher for symptomatic patients (Fig 2). Indeed, in this group of patients, having a family history of BrS was identified as a factor for increased PPV discovery yield. In the case of absence of BrS in the family, the variation discovery yield was almost double for those patients having a spontaneous type 1 BrS ECG than for patients with drug-induced type 1 ECG (45.5% vs 25%, respectively). In addition, we identified a PPV in 44.4% of the asymptomatic patients who presented family history of BrS and a spontaneous type 1 BrS ECG. When the patient presented drug-induced type 1 ECG or in the absence of family history of BrS, the PPV discovery yield was of around 15%.
Fig 2. Influence of the phenotype on PPV discovery yield.
Bar graph comparing the PPV detection yield in 8 different clinical categories (stated below the graph). Each bar shows the total number of patients for each clinical category divided in those with a PPV (black) and those without an identified PPV (white). The number of patients (in brackets) and percentages are given. Pos, positive; Neg, negative; Spont, spontaneous type 1 BrS ECG; Drug, drug-induced type 1 BrS ECG; n, number of patients.
We also investigated the role of age on the PPV occurrence. No significant age differences were observed between variation carriers and non-carriers (38.6±10.3 and 43.5±14.4, respectively, p = 0.16). However, the PPV discovery yield was higher for patients with ages between 30 and 50 years: out of the total of patients carrying a PPV, 83.3% of the patients were in this age range, while 11.1% were younger and 5.6% were older patients (Fig 3A, upper panel). The PPV discovery yield was significantly higher for symptomatic than for asymptomatic patients (42.3% vs 24.1%, respectively; Fig 3A, lower panels).
Fig 3. Influence of the age on PPVs discovery yield.
(A) Pie charts showing the distribution of patients in the overall population as well as in the categories of symptomatic and asymptomatic patients regarding PPV discovery. The percentage and the number of patients (in brackets) are given for each group. The small pie charts correspond to the age distribution of patients with an identified PPV. (B) Bar graphs of the PPV detection yields obtained for each of the age groups (< 30 years, 30–50 years and > 50 years). Numbers inside each bar correspond to the number of patients carrying a PPV for each category and the percentages represent the variation detection yield.
Noteworthy, in the 30–50 age range, 52.9% (9/17) of the symptomatic patients and 35.3% (6/17) of asymptomatic patients carried one PPV (Fig 3B, middle). Additionally, 40% (2/5) of the symptomatic young patients (< 30 years) were variation carriers, while no PPVs were identified in asymptomatic patients within this age range.
Overall, 55 unrelated Spanish patients clinically diagnosed with BrS were included in our study.Table 1 shows the demographics of this cohort, and Table 2 and S1 Table show the clinical and genetic characteristics of all the patients included in the study. The mean age at clinical diagnosis was of 41.9±13.3 years. Although the majority of patients were males (74.5%), their age at diagnosis was not different than that of females (41.8±12.1 years and 42.3±16.3 years, respectively; p = 0.92). A type 1 BrS ECG was present spontaneously in 37 patients (67.3%), and drug challenge revealed a type 1 BrS ECG for the remaining 18 patients (32.7%). Almost half of the patients had experienced symptoms, including 2 SCD and 4 aborted SCD. Patients who had not previously experienced any signs of arrhythmogenicity despite having a BrS ECG were considered asymptomatic. Comparison of symptomatic vs asymptomatic patients evidenced a similar percentage of males (73.1% and 75.9%, respectively). However, the mean age at diagnosis was different between the two groups of patients (37.7±14.3 and 45.7±11.4, respectively; p<0.05).
Discussion
To the best of our knowledge, this is the first comprehensive genetic evaluation of 14 BrS-susceptibility genes and MLPA of SCN5A in a Spanish cohort. Well delimited BrS cohorts from Japan, China, Greece and even Spain have been genetically studied [24,30–32]. Additionally, an international compendium of BrS genetic variations identified in more than 2100 unrelated patients from different countries was published in 2010 [3]. However, all these studies screenedSCN5A exclusively. In 2012, Crotti et al. reported the spectrum and prevalence of genetic variations in 12 BrS-susceptibility genes in a BrS cohort [5]. However, this study included patients of different ethnicity. Here, we report the analysis of 14 genes which has been conducted on a well-defined BrS cohort of the same ethnicity.
Our results confirm that SCN5A is still the most prevalent gene associated with BrS. Indeed,SCN5A-mediated BrS in our cohort (30.9%) is higher than the proportion described in other European reports [3,23], where a potentially causative variation is identified in only 20–25% of BrS patients. The reason for this discrepancy is unclear but could point towards a higher prevalence of SCN5A PPVs in the Spanish population or to selection bias. Additionally, we identified a genetic variation in SCN2B (c.632A>G, which results in p.D211G). We have formerly published the comprehensive electrophysiological characterization of this variation, and showed that indeed this variation could be responsible of the phenotype of the patient, thus linking SCN2B with BrS for the first time [7]. Also, we identified a variation in RANGRF. This variation (c.181G>T leading to p.E61X) had been previously reported in a Danish atrial fibrillation cohort [33]. Surprisingly, the authors reported an incidence of 0.4% for this variation in the healthy Danish population, which brought into question its pathogenicity. Our finding of this variation in an asymptomatic patient displaying a type 2 BrS ECG also points toward considering it as a rare genetic variation with a potential modifier effect on the phenotype but not clearly responsible for the disease [29].
No PPVs were identified in the other genes tested. Certainly, it is well accepted that the contribution of these genes to the disease is minor, and thus should only be considered under special circumstances [13,34]. In addition, recent studies have questioned the causality of variations identified in some of these minority genes [35].
We also used the MLPA technique for the detection of large exon duplications and/or deletions in SCN5A in patients without PPVs, and no large rearrangements were identified. This is in accordance with previous reports, which revealed that such imbalances are uncommon [8–10].
Kapplinger et al. [3] reported a predominance of PPVs in transmembrane regions of Nav1.5. Indeed, it has been proposed that most rare genetic variations in interdomain linkers may be considered as non-pathogenic [36]. In contrast, PPVs identified in this study are mainly located in extracellular loops and cytosolic linker regions of Nav1.5 (Fig 4). Additionally, 2 of our non-previously reported frameshifts are located in the DI-DII linker. These 2 genetic variations lead to truncated proteins, which would lack around 75% of the protein sequence, and thus are presupposed to be pathogenic.
Fig 4. Nav1.5 channel scheme showing the relative position of the SCN5A PPVs identified in our cohort.
Open symbols indicate already described variations and closed symbols locate novel variations reported in this study. DI to DIV designate the 4 domains of the protein, and numbers 1–6 identify the different segments within each domain. Crosses mark the voltage sensor.
In our cohort, we have identified 40 synonymous or common genetic variations, 4 of which have not been previously reported. These variations are gradually becoming more and more important in the explanation of certain phenotypes of genetic diseases. Only a few common variations identified here are already published as phenotypic modifiers [37,38]. The effect of these and other common variants identified in our cohort on BrS phenotype should be further studied.
Unexpectedly, almost 40% (7/18) of the PPV carriers did not present signs of arrhythmogenicity. We also performed genotype-phenotype correlations of the PPVs identified in the families (S3 Table). These studies uncovered relatives, most of whom were young individuals, who carried a familial variation but had never exhibited any clinical manifestations of the disease. This is in agreement with Crotti et al. and Priori et al. [5,23], who postulated that a positive genetic testing result is not always associated with the presence of symptoms. Indeed, the existence of asymptomatic patients carrying genetic variations described to cause a severe Nav1.5 channel dysfunction has been reported [39]. The identification of silent carriers is of paramount importance since it allows the adoption of preventive measures before any lethal episode takes place. Unknown environmental factors, medication and modifier genes have been suggested to influence and/or predispose to arrhythmogenesis [11]. Hence, this group of patients has to be cautiously followed in order to avoid fatal events.
Our studies on the connection between patients’ phenotype and the PPV detection yield highlighted the presence of symptoms as a factor for an increased variation discovery yield. Within the group of symptomatic individuals, a PPV was identified in a higher proportion of patients displaying a spontaneous type 1 BrS ECG than for patients showing a drug-induced ECG. Likewise, within the asymptomatic patients with family history of BrS, those who presented spontaneous type 1 BrS ECG carried a PPV more often than those with a drug-induced ECG (Fig 2). Referring to age, the vast majority (17/20, 85%) of the PPVs were identified in patients around their fourth decade of age (30–50 years). This is in accordance with the accepted mean age of disease manifestation. Moreover, in this age range, more than 50% of the patients who presented symptoms carried a variation that could be pathogenic (Fig 3). Importantly, 35.3% of asymptomatic patients of around 40 years of age also carried one of such variations. These data highlight the importance of performing a genetic test even in the absence of clinical manifestations of the disease, and particularly when in the 30–50 years range, which is in accordance with consensus recommendations [13,34].
In conclusion, we have analysed for the first time 14 BrS-susceptibility genes and performed MLPA of SCN5A in a Spanish BrS cohort. Our cohort showed male prevalence with a mean age of disease manifestation around 40 years. BrS in this cohort was almost exclusivelySCN5A-mediated. The mean PPV discovery yield in our Spanish BrS patients is higher than that described for other BrS cohorts (32.7% vs 20–25%, respectively), and is even higher for patients in the 30–50 years age range (up to 53% for symptomatic patients). All these evidences support the genetic testing, at least of SCN5A, in all clinically well diagnosed BrS patients.
Study Limitations
First of all, drug challenge tests were not performed for all the relatives who were asymptomatic variation carriers. This fact hampered their clinical diagnosis and represents an impediment to definitely assess the link between PPVs and BrS. These patients are nowadays under follow-up.
New PPVs have been identified in our cohort. The clinical information available for the families suggests that these new variations could be pathogenic. Still, in vitro studies of these variations are required in order to evaluate their functional effects and verify their pathogenic role. Additionally, genotyping in an independent cohort would help reduce the likelihood of type I (false positive) error in genetic variant discovery.
We have to acknowledge that the study set is relatively small. Consequently, the classification of patients according to the different clinical categories rendered rather small sub-groups, which may lead to over-interpretation of the results. Future studies will be directed to the genetic screening of additional Spanish BrS patients, which will probably reinforce the significance of the tendencies observed here.
Sets of co-expressed Genes influence Blood Pressure Regulation: Genome-wide Association and mRNA expression @US National Heart, Lung, and Blood Institute
Reporter: Aviva Lev-Ari, PhD, RN
NHLBI-led Team Untangles Gene Networks Involved in Blood Pressure Regulation
NEW YORK (GenomeWeb) – Using network approaches, researchers from the US National Heart, Lung, and Blood Institute and their colleagues combined genome-wide association and mRNA expression data to home in on sets of co-expressed genes that appear to influence blood pressure regulation.
As they reported in Molecular Systems Biology today, the NHLBI-led team drew on data from more than 3,600 people participating in the Framingham Heart Study to identify four potentially causal gene modules and key driver genes contained within them.
“Our work was able to pinpoint several gene networks closely linked to the regulation of blood pressure,” first author Tianxiao Huan from NHLBI said in a statement.
In addition, Huan and her colleagues traced the function of one key driver gene — SH2B3 — to response to angiotensin II infusion in a mouse model, indicating that this approach may help identify new treatment targets.
For this study, Huan and her colleagues examined the gene expression profiles of 3,679 Framingham Heart Study participants of European descent who were not taking an antihypertensive drug. They correlated gene expression changes they observed in this cohort with systolic blood pressure, diastolic blood pressure, and hypertension and, after accounting for age, BMI, gender, and other factors, came up with 83 associated genes.
At the same time, the researchers constructed gene co-expression networks from that gene expression data to develop gene co-expression network modules that they then also correlated to blood pressure phenotypes. Of these 27 gene co-expression network modules, seven were significantly associated with either systolic or diastolic blood pressure, the researchers said.
While that set of 83 blood pressure-related genes wasn’t significantly enriched for any gene ontology terms, the seven gene co-expression network modules were linked to a variety of functions, including chromatin modification, immune cell-mediated cytotoxicity, inflammatory response, and more. This suggested to the researchers that genes involved in a range of biological processes are tightly co-regulated with respect to blood pressure.
Using a SNP set enrichment analysis approach, the researchers found that four of the gene co-expression network modules appeared to be potentially causal and that more than a dozen genes in those modules appeared to contribute to their association with blood pressure regulation.
For instance, one SNP, dubbed rs3184504, had been linked with blood pressure through a genome-wide association study, and it is linked with the expression of four genes in the set of genetically inferred causal blood pressure genes.
Using blood Bayesian networks and protein-protein interaction networks other groups had developed, Huan and her colleagues further zoomed in on key driver genes by testing whether the surrounding region of each gene in those four gene co-expression network modules was enriched for other potentially causal blood pressure genes.
These top key driver genes, they noted, were involved in subnetworks that appeared to regulate blood pressure-related genes.
For example, a missense SNP in an exon of SH2B3 has been associated with blood pressure and hypertension in a GWAS and is linked to expression changes in 10 other genes the researchers identified. These genes, Huan and her colleagues said, were enriched for activity in the intracellular signaling cascade, T-cell activation, and T-cell differentiation. This SH2B3-subnework was also enriched for genes known to be linked to blood pressure.
Previous work had linked SH2B3 to blood pressure regulation, Huan and her colleagues said, but how it had its effect wasn’t clear.
Mice lacking the SH2B3 gene, they noted, had normal baseline blood pressure, though it became elevated in response to a low dose of angiotensin II, an effect not seen in wild-type mice.
In addition, RNA sequencing of the whole-blood transcriptomes from wild-type and Sh2b3-/- mice indicated that more than 2,240 genes were differentially expressed between the two, especially ones involved in immune and inflammatory response. These genes significantly overlapped with the SH2B3 genetic subnetwork, and those overlapping genes were enriched for ones involved in the intracellular signaling cascade and T-cell activation, the researchers reported.
“Moving forward, it should be possible to study additional key driver genes in this way, which should help in our efforts to identify novel targets for the prevention and treatment of hypertension,” Huan added.
This curation deals with a large and specialized branch of medicine that grew since the mid 20th century in concert with the developments in genetics and as a result of a growing population, with large urban populations, increasing problems of premature deliveries. The problems of prematurity grew very preterm to very low birth weight babies with special problems. While there were nurseries, the need for intensive care nurseries became evident in the 1960s, and the need for perinatal care of pregnant mothers also grew as a result of metabolic problems of the mother, intrauterine positioning of the fetus, and increasing numbers of teen age pregnancies as well as nutritional problems of the mother. There was also a period when the manufacturers of nutritional products displaced the customary use of breast feeding, which was consequential. This discussion is quite comprehensive, as it involves a consideration of the heart, the lungs, the brain, and the liver, to a large extent, and also the kidneys and skeletal development.
It is possible to outline, with a proportionate emphasis based on frequency and severity, this as follows:
Genetic and metabolic diseases
Nervous system
Cardiovascular
Pulmonary
Skeletal – bone and muscle
Hematological
Liver
Esophagus, stomach, and intestines
Kidneys
Immune system
Fetal Development
Gestation is the period of time between conception and birth when a baby grows and develops inside the mother’s womb. Because it’s impossible to know exactly when conception occurs, gestational age is measured from the first day of the mother’s last menstrual cycle to the current date. It is measured in weeks. A normal gestation lasts anywhere from 37 to 41 weeks.
Week 5 is the start of the “embryonic period.” This is when all the baby’s major systems and structures develop. The embryo’s cells multiply and start to take on specific functions. This is called differentiation. Blood cells, kidney cells, and nerve cells all develop. The embryo grows rapidly, and the baby’s external features begin to form.
Week 6-9: Brain forms into five different areas. Some cranial nerves are visible. Eyes and ears begin to form. Tissue grows that will the baby’s spine and other bones. Baby’s heart continues to grow and now beats at a regular rhythm. Blood pumps through the main vessels. Your baby’s brain continues to grow. The lungs start to form. Limbs look like paddles. Essential organs begin to grow.
Weeks 11-18: Limbs extended. Baby makes sucking motion. Movement of limbs. Liver and pancreas produce secretions. Muscle and bones developing.
Week 19-21: Baby can hear. Mom feels baby – and quickening.
The original cause of most genetic metabolic disorders is a gene mutation that occurred many, many generations ago. The gene mutation is passed along through the generations, ensuring its preservation.
Each inherited metabolic disorder is quite rare in the general population. Considered all together, inherited metabolic disorders may affect about 1 in 1,000 to 2,500 newborns. In certain ethnic populations, such as Ashkenazi Jews (Jews of central and eastern European ancestry), the rate of inherited metabolic disorders is higher.
Hundreds of inherited metabolic disorders have been identified, and new ones continue to be discovered. Some of the more common and important genetic metabolic disorders include:
Lysosomal storage disorders : Lysosomes are spaces inside cells that break down waste products of metabolism. Various enzyme deficiencies inside lysosomes can result in buildup of toxic substances, causing metabolic disorders including:
Hurler syndrome (abnormal bone structure and developmental delay)
Tay-Sachs disease (progressive weakness in a months-old child, progressing to severe nerve damage; the child usually lives only until age 4 or 5)
Gauchers disease and others
Galactosemia: Impaired breakdown of the sugar galactose leads to jaundice, vomiting, and liver enlargement after breast or formula feeding by a newborn.
Maple syrup urine disease: Deficiency of an enzyme called BCKD causes buildup of amino acids in the body. Nerve damage results, and the urine smells like syrup.
Phenylketonuria (PKU): Deficiency of the enzyme PAH results in high levels of phenylalanine in the blood. Mental retardation results if the condition is not recognized.
Glycogen storage diseases: Problems with sugar storage lead to low blood sugar levels, muscle pain, and weakness.
Metal metabolism disorders: Levels of trace metals in the blood are controlled by special proteins. Inherited metabolic disorders can result in protein malfunction and toxic accumulation of metal in the body:
Wilson disease (toxic copper levels accumulate in the liver, brain, and other organs)
Hemochromatosis (the intestines absorb excessive iron, which builds up in the liver, pancreas, joints, and heart, causing damage)
Organic acidemias: methylmalonic acidemia and propionic acidemia.
Urea cycle disorders: ornithine transcarbamylase deficiency and citrullinemia
New variations in the galactose-1-phosphate uridyltransferase (GALT) gene
Clinical and molecular spectra in galactosemic patients from neonatal screening in northeastern Italy: Structural and functional characterization of new variations in the galactose-1-phosphate uridyltransferase (GALT) gene
E Viggiano, A Marabotti, AP Burlina, C Cazzorla, MR D’Apice, et al.
Gene 559 (2015) 112–118 http://dx.doi.org/10.1016/j.gene.2015.01.013
Galactosemia (OMIM 230400) is a rare autosomal recessive inherited disorder caused by deficiency of galactose-1-phosphate uridyltransferase (GALT; OMIM 606999) activity. The incidence of galactosemia is 1 in 30,000–60,000, with a prevalence of 1 in 47,000 in the white population. Neonates with galactosemia can present acute symptoms, such as severe hepatic and renal failure, cataract and sepsis after milk introduction. Dietary restriction of galactose determines the clinical improvement in these patients. However, despite early diagnosis by neonatal screening and dietary treatment, a high percentage of patients develop long-term complications such as cognitive disability, speech problems, neurological and/or movement disorders and, in females, ovarian dysfunction.
With the benefit of early diagnosis by neonatal screening and early therapy, the acute presentation of classical galactosemia can be prevented. The objectives of the current study were to report our experience with a group of galactosemic patients identified through the neonatal screening programs in northeastern Italy during the last 30 years.
No neonatal deaths due to galactosemia complications occurred after the introduction of the neonatal screening program. However, despite the early diagnosis and dietary treatment, the patients with classical galactosemia showed one or more long-term complications.
A total of 18 different variations in the GALT gene were found in the patient cohort: 12 missense, 2 frameshift, 1 nonsense, 1 deletion, 1 silent variation, and 1 intronic. Six (p.R33P, p.G83V, p.P244S, p.L267R, p.L267V, p.E271D) were new variations. The most common variation was p.Q188R (12 alleles, 31.5%), followed by p.K285N (6 alleles, 15.7%) and p.N314D (6 alleles, 15.7%). The other variations comprised 1 or 2 alleles. In the patients carrying a new mutation, the biochemical analysis of GALT activity in erythrocytes showed an activity of < 1%. In silico analysis (SIFT, PolyPhen-2 and the computational analysis on the static protein structure) showed potentially damaging effects of the six new variations on the GALT protein, thus expanding the genetic spectrum of GALT variations in Italy. The study emphasizes the difficulty in establishing a genotype–phenotype correlation in classical galactosemia and underlines the importance of molecular diagnostic testing prior to making any treatment.
Diagnosis and Management of Hereditary Hemochromatosis
Hereditary hemochromatosis (HH) is a diagnosis most commonly made in patients with elevated iron indices (transferrin saturation and ferritin), and HFE genetic mutation testing showing C282Y homozygosity.
The HFE mutation is believed to result in clinical iron overload through altering hepcidin levels resulting in increased iron absorption.
The most common clinical complications of HH include cirrhosis, diabetes, nonischemic cardiomyopathy, and hepatocellular carcinoma.
Liver biopsy should be performed in patients with HH if the liver enzymes are elevated or serum ferritin is greater than 1000 mg/L. This is useful to determine the degree of iron overload and stage the fibrosis.
Treatment of HH with clinical iron overload involves a combination of phlebotomy and/or chelation therapy. Liver transplantation should be considered for patients with HH-related decompensated cirrhosis.
Health economic evaluation of plasma oxysterol screening in the diagnosis of Niemann–Pick Type C disease among intellectually disabled using discrete event simulation
Background: Recently a less invasive method of screening and diagnosing Niemann–Pick C (NP-C) disease has emerged. This approach involves the use of a metabolic screening test (oxysterol assay) instead of the current practice of clinical assessment of patients suspected of NP-C (review of medical history, family history and clinical examination for the signs and symptoms). Our objective is to compare costs and outcomes of plasma oxysterol screening versus current practice in diagnosis of NP-C disease among intellectually disabled (ID) patients using decision-analytic methods. Methods: A discrete event simulation model was conducted to follow ID patients through the diagnosis and treatment of NP-C, forecast the costs and effectiveness for a cohort of ID patients and compare the outcomes and costs in two different arms of the model: plasma oxysterol screening and routine diagnosis procedure (anno 2013) over 5 years of follow up. Data from published sources and clinical trials were used in simulation model. Unit costs and quality-adjusted life-years (QALYs) were discounted at a 3% annual rate in the base case analysis. Deterministic and probabilistic sensitivity analyses were conducted. Results: The outcomes of the base case model showed that using plasma oxysterol screening for diagnosis of NP-C disease among ID patients is a dominant strategy. It would result in lower total cost and would slightly improve patients’ quality of life. The average amount of cost saving was $3642 CAD and the incremental QALYs per each individual ID patient in oxysterol screening arm versus current practice of diagnosis NP-C was 0.0022 QALYs. Results of sensitivity analysis demonstrated robustness of the outcomes over the wide range of changes in model inputs. Conclusion: Whilst acknowledging the limitations of this study, we conclude that screening ID children and adolescents with oxysterol tests compared to current practice for the diagnosis of NP-C is a dominant strategy with clinical and economic benefits. The less costly, more sensitive and specific oxysterol test has potential to save costs to the healthcare system while improving patients’ quality of life and may be considered as a routine tool in the NP-C diagnosis armamentarium for ID. Further research is needed to elucidate its effectiveness in patients presenting characteristics other than ID in childhood and adolescence.
Neurological and Behavioral Disorders
Estrogen receptor signaling during vertebrate development
Estrogen receptors are expressed and their cognate ligands produced in all vertebrates, indicative of important and conserved functions. Through evolution estrogen has been involved in controlling reproduction, affectingboth the development of reproductive organs and reproductive behavior. This review broadly describes the synthesis of estrogens and the expression patterns of aromatase and the estrogen receptors, in relation to estrogen functions in the developing fetus and child. We focus on the role of estrogens for the development of reproductive tissues, as well as non-reproductive effects on the developing brain. We collate data from human, rodent, bird and fish studies and highlight common and species-specific effects of estrogen signaling on fetal development. Morphological malformations originating from perturbed estrogen signaling in estrogen receptor and aromatase knockout mice are discussed, as well as the clinical manifestations of rare estrogen receptor alpha and aromatase gene mutations in humans. This article is part of a Special Issue entitled: Nuclear receptors in animal development.
Memory function and hippocampal volumes in preterm born very-low-birth-weight (VLBW) young adults
The hippocampi are regarded as core structures for learning and memory functions, which is important for daily functioning and educational achievements. Previous studies have linked reduction in hippocampal volume to working memory problems in very low birth weight (VLBW; ≤1500 g) children and reduced general cognitive ability in VLBW adolescents. However, the relationship between memory function and hippocampal volume has not been described in VLBW subjects reaching adulthood. The aim of the study was to investigate memory function and hippocampal volume in VLBW young adults, both in relation to perinatal risk factors and compared to term born controls, and to look for structure–function relationships. Using Wechsler Memory Scale-III and MRI, we included 42 non-disabled VLBW and 61 control individuals at age 19–20 years, and related our findings to perinatal risk factors in the VLBW-group. The VLBW young adults achieved lower scores on several subtests of the Wechsler Memory Scale-III, resulting in lower results in the immediate memory indices (visual and auditory), the working memory index, and in the visual delayed and general memory delayed indices, but not in the auditory delayed and auditory recognition delayed indices. The VLBW group had smaller absolute and relative hippocampal volumes than the controls. In the VLBW group inferior memory function, especially for the working memory index, was related to smaller hippocampal volume, and both correlated with lower birth weight and more days in the neonatal intensive care unit (NICU). Our results may indicate a structural–functional relationship in the VLBW group due to aberrant hippocampal development and functioning after preterm birth.
The relation of infant attachment to attachment and cognitive and behavioural outcomes in early childhood
Background: In China, research on the relation of mother–infant attachment to children’s development is scarce. Aims: This study sought to investigate the relation of mother–infant attachment to attachment, cognitive and behavioral development in young children. Study design: This study used a longitudinal study design. Subjects: The subjects included healthy infants (n=160) aged 12 to 18 months. Outcome measures: Ainsworth’s “Strange Situation Procedure” was used to evaluate mother–infant attachment types. The attachment Q-set (AQS) was used to evaluate the attachment between young children and their mothers. The Bayley scale of infant development-second edition (BSID-II) was used to evaluate cognitive developmental level in early childhood. Achenbach’s child behavior checklist (CBCL) for 2- to 3-year-oldswas used to investigate behavioral problems. Results: In total, 118 young children (73.8%) completed the follow-up; 89.7% of infants with secure attachment and 85.0% of infants with insecure attachment still demonstrated this type of attachment in early childhood (κ = 0.738, p b 0.05). Infants with insecure attachment collectively exhibited a significantly lower mental development index (MDI) in early childhood than did infants with secure attachment, especially the resistant type. In addition, resistant infants were reported to have greater social withdrawal, sleep problems and aggressive behavior in early childhood. Conclusion: There is a high consistency in attachment development from infancy to early childhood. Secure mother–infant attachment predicts a better cognitive and behavioral outcome; whereas insecure attachment, especially the resistant attachment, may lead to a lower cognitive level and greater behavioral problems in early childhood.
representations of the HPA axis
representations of limbic stress-integrative pathways from the prefrontal cortex, amygdala and hippocampus
Fetal programming of schizophrenia: Select mechanisms
Mounting evidence indicates that schizophrenia is associated with adverse intrauterine experiences. An adverse or suboptimal fetal environment can cause irreversible changes in brain that can subsequently exert long-lasting effects through resetting a diverse array of biological systems including endocrine, immune and nervous. It is evident from animal and imaging studies that subtle variations in the intrauterine environment can cause recognizable differences in brain structure and cognitive functions in the offspring. A wide variety of environmental factors may play a role in precipitating the emergent developmental dysregulation and the consequent evolution of psychiatric traits in early adulthood by inducing inflammatory, oxidative and nitrosative stress (IO&NS) pathways, mitochondrial dysfunction, apoptosis, and epigenetic dysregulation. However, the precise mechanisms behind such relationships and the specificity of the risk factors for schizophrenia remain exploratory. Considering the paucity of knowledge on fetal programming of schizophrenia, it is timely to consolidate the recent advances in the field and put forward an integrated overview of the mechanisms associated with fetal origin of schizophrenia.
NMDA receptor dysfunction in autism spectrum disorders
Autism spectrum disorders (ASDs) represent neurodevelopmental disorders characterized by two core symptoms;
(1) impaired social interaction and communication, and
(2) restricted and repetitive behaviors, interests, and activities.
ASDs affect ~ 1% of the population, and are considered to be highly genetic in nature. A large number (~600) of ASD-related genetic variations have been identified (sfari.org), and target gene functions are apparently quite diverse. However, some fall onto common pathways, including synaptic function and chromosome remodeling, suggesting that core mechanisms may exist.
Abnormalities and imbalances in neuronal excitatory and inhibitory synapses have been implicated in diverse neuropsychiatric disorders including autism spectrum disorders (ASDs). Increasing evidence indicates that dysfunction of NMDA receptors (NMDARs) at excitatory synapses is associated with ASDs. In support of this, human ASD-associated genetic variations are found in genes encoding NMDAR subunits. Pharmacological enhancement or suppression of NMDAR function ameliorates ASD symptoms in humans. Animal models of ASD display bidirectional NMDAR dysfunction, and correcting this deficit rescues ASD-like behaviors. These findings suggest that deviation of NMDAR function in either direction contributes to the development of ASDs, and that correcting NMDAR dysfunction has therapeutic potential for ASDs.
Among known synaptic proteins implicated in ASD are metabotropic glutamate receptors (mGluRs). Functional enhancement and suppression of mGluR5 are associated with fragile X syndrome and tuberous sclerosis, respectively, which share autism as a common phenotype. More recently, ionotropic glutamate receptors, namely NMDA receptors (NMDARs) and AMPA receptors (AMPARs), have also been implicated in ASDs. In this review, we will focus on NMDA receptors and summarize evidence supporting the hypothesis that NMDAR dysfunction contributes to ASDs, and, by extension, that correcting NMDAR dysfunction has therapeutic potential for ASDs. ASD-related human NMDAR genetic variants.
Chemokines roles within the hippocampus
IL-1 mediates stress-induced activation of the HPA axis
A systemic model of the beneficial role of immune processes in behavioral and neural plasticity
Three Classes of Glutamate Receptors
Clinical studies on ASDs have identified genetic variants of NMDAR subunit genes. Specifically, de novo mutations have been identified in the GRIN2B gene, encoding the GluN2B subunit. In addition, SNP analyses have linked both GRIN2A (GluN2A subunit) and GRIN2B with ASDs. Because assembled NMDARs contain four subunits, each with distinct properties, ASD-related GRIN2A/ GRIN2B variants likely alter the functional properties of NMDARs and/or NMDAR-dependent plasticity.
Pharmacological modulation of NMDAR function can improve ASD symptoms. D-cycloserine (DCS), an NMDAR agonist, significantly ameliorates social withdrawal and repetitive behavior in individuals with ASD. These results suggest that reduced NMDAR function may contribute to the development of ASDs in humans.
We can divide animal studies into two groups. The first group consists of animals in which NMDAR modulators were shown to normalize both NMDAR dysfunction and ASD-like behaviors, establishing strong association between NMDARs and ASD phenotypes (Fig.). In the second group, NMDAR modulators were shown to rescue ASD-like behaviors, but NMDAR dysfunction and its correction have not been demonstrated.
ASD models with data showing rescue of both NMDAR dysfunction and ASD like behaviors Mice lacking neuroligin-1, an excitatory postsynaptic adhesion molecule, show reduced NMDAR function in the hippocampus and striatum, as evidenced by a decrease in NMDA/AMPA ratio and long-term potentiation (LTP). Neuroligin-1 is thought to enhance synaptic NMDAR function, by directly interacting with and promoting synaptic localization of NMDARs.
Fig not shown.
Bidirectional NMDAR dysfunction in animal models of ASD. Animal models of ASD with bidirectional NMDAR dysfunction can be positioned on either side of an NMDAR function curve. Model animals were divided into two groups.
Group 1: NMDAR modulators normalize both NMDAR dysfunction and ASD-like behaviors (green).
Group 2: NMDAR modulators rescue ASD-like behaviors, but NMDAR dysfunction and its rescue have not been demonstrated (orange). Note that Group 2 animals are tentatively placed on the left-hand side of the slope based on the observed DCS rescue of their ASD-like phenotypes, but the directions of their NMDAR dysfunctions remain to be experimentally determined.
ASD models with data showing rescue of ASD-like behaviors but no demonstrated NMDAR dysfunction
Tbr1 is a transcriptional regulator, one of whose targets is the gene encoding the GluN2B subunit of NMDARs. Mice haploinsufficient for Tbr1 (Tbr1+/-) show structural abnormalities in the amygdala and limited GluN2B induction upon behavioral stimulation. Both systemic injection and local amygdalar infusion of DCS rescue social deficits and impaired associative memory in Tbr1+/- mice. However, reduced NMDAR function and its DCS-dependent correction have not been demonstrated.
Spatial working memory and attention skills are predicted by maternal stress during pregnancy
Introduction: Experimental evidence in rodents shows that maternal stress during pregnancy (MSDP) negatively impacts spatial learning and memory in the offspring. We aim to investigate the association between MSDP (i.e., life events) and spatial working memory, as well as attention skills (attention shifting and attention focusing), in humans. The moderating roles of child sex, maternal anxiety during pregnancy and postnatal care are also investigated. Methods: Participants were 236mother–child dyads that were followed from the second trimester of pregnancy until 4 years postpartum. Measurements included questionnaires and independent observations. Results: MSDP was negatively associated with attention shifting at 18monthswhen concurrent maternal anxiety was low. MSDP was associated with poorer spatial working memory at 4 years of age, but only for boys who experienced poorer postnatal care. Conclusion: Consistent with results observed in rodents, MSDP was found to be associated with spatial working memory and attention skills. These results point to postnatal care and maternal anxiety during pregnancy as potential targets for interventions that aim to buffer children from the detrimental effects of MSDP.
Acute and massive bleeding from placenta previa and infants’ brain damage
Background: Among the causes of third trimester bleeding, the impact of placenta previa on cerebral palsy is not well known. Aims: To clarify the effect ofmaternal bleeding fromplacenta previa on cerebral palsy, and in particular when and how it occurs. Study design: A descriptive study. Subjects: Sixty infants born to mothers with placenta previa in our regional population-based study of 160,000 deliveries from 1998 to 2012. Premature deliveries occurring atb26 weeks of gestation and placenta accrete were excluded. Outcome measures: Prevalence of cystic periventricular leukomalacia (PVL) and cerebral palsy (CP). Results: Five infants had PVL and 4 of these infants developed CP (1/40,000 deliveries). Acute and massive bleeding (>500 g) within 8 h) occurred at around 30–31 weeks of gestation, and was severe enough to deliver the fetus. None of the 5 infants with PVL underwent antenatal corticosteroid treatment, and 1 infant had mild neonatal hypocapnia with a PaCO2 < 25 mm Hg. However, none of the 5 PVL infants showed umbilical arterial academia with pH < 7.2, an abnormal fetal heart rate monitoring pattern, or neonatal hypotension. Conclusions: Our descriptive study showed that acute and massive bleeding from placenta previa at around 30 weeks of gestation may be a risk factor for CP, and requires careful neonatal follow-up. The underlying process connecting massive placental bleeding and PVL requires further investigation.
Impact of bilirubin-induced neurologic dysfunction on neurodevelopmental outcomes
Extreme neonatal hyperbilirubinemia has long been known to cause the clinical syndrome of kernicterus, or chronic bilirubin encephalopathy (CBE). Kernicterus most usually is characterized by choreoathetoid cerebral palsy (CP), impaired upward gaze, and sensorineural hearing loss, whereas cognition is relatively spared. The chronic condition of kernicterus may be, but is not always, preceded in the acute stage by acute bilirubin encephalopathy (ABE). This acute neonatal condition is also due to hyperbilirubinemia, and is characterized by lethargy and abnormal behavior, evolving to frank neonatal encephalopathy, opisthotonus, and seizures. Less completely defined is the syndrome of bilirubin-induced neurologic dysfunction (BIND).
Bilirubin-induced neurologic dysfunction (BIND) is the constellation of neurologic sequelae following milder degrees of neonatal hyperbilirubinemia than are associated with kernicterus. Clinically, BIND may manifest after the neonatal period as developmental delay, cognitive impairment, disordered executive function, and behavioral and psychiatric disorders. However, there is controversy regarding the relative contribution of neonatal hyperbilirubinemia versus other risk factors to the development of later neurodevelopmental disorders in children with BIND. In this review, we focus on the empiric data from the past 25 years regarding neurodevelopmental outcomes and BIND, including specific effects on developmental delay, cognition, speech and language development, executive function, and the neurobehavioral disorders, such as attention deficit/hyperactivity disorder and autism.
As noted in a technical report by the American Academy of Pediatrics Subcommittee on Hyperbilirubinemia, “it is apparent that the use of a single total serum bilirubin level to predict long-term outcomes is inadequate and will lead to conflicting results”. As described above, this has certainly been the case in research to date. To clarify how hyperbilirubinemia influences neurodevelopmental outcome, more sophisticated consideration is needed both of how to assess bilirubin exposure leading to neurotoxicity, and of those comorbid conditions which may lower the threshold for brain injury.
For example, premature infants are known to be especially susceptible to bilirubin neurotoxicity, with kernicterus reported following TB levels far lower than the threshold expected in term neonates. Similarly, among extremely preterm neonates, BBC is proportional to gestational age, meaning that the most premature infants have the highest UB, even for similar TB levels. Thus, future studies must be adequately powered to examine preterm infants separately from term infants, and should consider not just peak TB, but also BBC, as independent variables in neonates with hyperbilirubinemia. Similarly, an analysis by the NICHD NRN found that, among ELBW infants, higher UB levels were associated with a higher risk of death or NDI. However, increased TB levels were only associated with death or NDI in unstable infants. Again, UB or BBC appeared to be more useful than TB.
Are the neuromotor disabilities of bilirubin-induced neurologic dysfunction disorders related to the cerebellum and its connections?
Investigators have hypothesized a range of subcortical neuropathology in the genesis of bilirubin induced neurologic dysfunction (BIND). The current review builds on this speculation with a specific focus on the cerebellum and its connections in the development of the subtle neuromotor disabilities of BIND. The focus on the cerebellum derives from the following observations:
(i) the cerebellum is vulnerable to bilirubin-induced injury; perhaps the most vulnerable region within the central nervous system;
(ii) infants with cerebellar injury exhibit a neuromotor phenotype similar to BIND; and (iii) the cerebellum has extensive bidirectional circuitry projections to motor and non-motor regions of the brain-stem and cerebral cortex that impact a variety of neurobehaviors.
Future study using advanced magnetic resonance neuroimaging techniques have the potential to shed new insights into bilirubin’s effect on neural network topology via both structural and functional brain connectivity measurements.
Bilirubin-induced neurologic damage is most often thought of in terms of severe adverse neuromotor (dystonia with or without athetosis) and auditory (hearing impairment or deafness) sequelae. Observed together, they comprise the classic neurodevelopmental phenotype of chronic bilirubin encephalopathy or kernicterus, and may also be seen individually as motor or auditory predominant subtypes. These injuries reflect both a predilection of bilirubin toxicity for neurons (relative to glial cells) and the regional topography of bilirubin-induced neuronal damage characterized by prominent involvement of the globus pallidus, subthalamic nucleus, VIII cranial nerve, and cochlear nucleus.
It is also asserted that bilirubin neurotoxicity may be associated with other less severe neurodevelopmental disabilities, a condition termed “subtle kernicterus” or “bilirubin-induced neurologic dysfunction” (BIND). BIND is defined by a constellation of “subtle neurodevelopmental disabilities without the classical findings of kernicterus that, after careful evaluation and exclusion of other possible etiologies, appear to be due to bilirubin neurotoxicity”. These purportedly include:
(i) mild-to-moderate disorders of movement (e.g., incoordination, clumsiness, gait abnormalities, disturbances in static and dynamic balance, impaired fine motor skills, and ataxia); (ii) disturbances in muscle tone; and
(iii) altered sensorimotor integration. Isolated disturbances of central auditory processing are also included in the spectrum of BIND.
Cerebellar vulnerability to bilirubin-induced injury
Cerebellar injury phenotypes and BIND
Cerebellar projections
Transverse section of cerebellum and brainstem
Transverse section of cerebellum and brain-stem from a 34 gestational-week premature kernicteric infant formalin-fixed for two weeks. Yellow staining is evident in the cerebellar dentate nuclei (upper arrow) and vestibular nuclei at the pontomedullary junction (lower arrowhead). Photo is courtesy of Mahmdouha Ahdab-Barmada and reprinted with permission from Taylor-Francis Group (Ahdab Barmada M. The neuropathology of kernicterus: definitions and debate. In: Maisel MJ, Watchko JF editors. Neonatal jaundice. Amsterdam: Harwood Academic Publishers; 2000. p. 75e88
Whether cerebellar injury is primal or an integral part of disturbed neural circuitry in bilirubin-induced CNS damage is unclear. Movement disorders, however, are increasingly recognized to arise from abnormalities of neuronal circuitry rather than localized, circumscribed lesions. The cerebellum has extensive bidirectional circuitry projections to an array of brainstem nuclei and the cerebral cortex that modulate and refine motor activities. In this regard, the cerebellum is characteristically subdivided into three lobes based on neuroanatomic and phylogenetic criteria as well as by their primary afferent and efferent connections. They include:
(i) flocculonodular lobe (archicerebellum);
(ii) anterior lobe (paleocerebellum); and
(iii) posterior lobe (neocerebellum).
The archicerebellum, the oldest division phylogenically, receives extensive input from the vestibular system and is therefore also known as the vestibulocerebellum and is important for equilibrium control. The paleocerebellum, also a primitive region, receives extensive somatosensory input from the spinal cord, including the anterior and posterior spinocerebellar pathways that convey unconscious proprioception, and is therefore also known as the spinocerebellum. The neocerebellum is the most recently evolved region, receives most of the input from the cerebral cortex, and is thus termed the cerebrocerebellum. This area has greatly expanded in association with the extensive development of the cerebral cortex in mammals and especially primates. To cause serious longstanding dysfunction, cerebellar injury must typically involve the deep cerebellar nuclei and their projections.
Schematic of the bidirectional connectivity between the cerebellum and other
Schematic of the bidirectional connectivity between the cerebellum and other brain regions including the cerebral cortex. Most cerebro-cerebellar afferent projections pass through the basal (anterior or ventral) pontine nuclei and intermediate cerebellar peduncle, whereas most cerebello-cerebral efferent projections pass through the dentate and ventrolateral thalamic nuclei. DCN, deep cerebellar nuclei; RN, red nucleus; ATN, anterior thalamic nucleus; PFC, prefrontal cortex; MC, motor cortex; PC, parietal cortex; TC, temporal cortex; STN, subthalamic nucleus; APN, anterior pontine nuclei. Reprinted under the terms of the Creative Commons Attribution License from D’Angelo E, Casali S. Seeking a unified framework for cerebellar function and dysfunction: from circuit to cognition. Front Neural Circuits 2013; 6:116.
Given the vulnerability of the cerebellum to bilirubin-induced injury, cerebellar involvement should also be evident in classic kernicterus, contributing to neuromotor deficits observed therein. It is of interest, therefore, that cerebellar damage may play a role in the genesis of bilirubin-induced dystonia, a prominent neuromotor feature of chronic bilirubin encephalopathy in preterm and term neonates alike. This complex movement disorder is characterized by involuntary sustained muscle contractions that result in abnormal position and posture. Moreover, dystonia that is brief in duration results in chorea, and, if brief and repetitive, leads to athetosis ‒ conditions also classically observed in kernicterus. Recent evidence suggests that dystonic movements may depend on disruption of both basal ganglia and cerebellar neuronal networks, rather than isolated dysfunction of only one motor system.
Dystonia is also a prominent feature in Gunn rat pups and neonatal Ugt1‒/‒-deficient mice both robust models of kernicterus. The former is used as an experimental model of dystonia. Although these models show basal ganglia injury, the sine qua non of bilirubin-induced murine neuropathology is cerebellar damage and resultant cerebellar hypoplasia.
Studies are needed to define more precisely the motor network abnormalities in kernicterus and BIND. Magnetic resonance imaging (MRI) has been widely used in evaluating infants at risk for bilirubin-induced brain injury using conventional structural T1-and T2-weighted imaging. Infants with chronic bilirubin encephalopathy often demonstrate abnormal bilateral, symmetric, high-signal intensity on T2-weighted MRI of the globus pallidus and subthalamic nucleus, consistent with the neuropathology of kernicterus. Early postnatal MRI of at-risk infants, although frequently showing increased T1-signal in these regions, may give false-positive findings due to the presence of myelin in these structures.
Diffusion tensor imaging and tractography could be used to delineate long-term changes involving specific white matter pathways, further elucidating the neural basis of long-term disability in infants and children with chronic bilirubin encephalopathy and BIND. It will be equally valuable to use blood oxygen level-dependent (BOLD) “resting state” functional MRI to study intrinsic connectivity in order to identify vulnerable brain networks in neonates with kernicterus and BIND. Structural networks of the CNS (connectome) and functional network topology can be characterized in infants with kernicterus and BIND to determine disease-related pattern(s) with respect to both long- and short-range connectivity. These findings have the potential to shed novel insights into the pathogenesis of these disorders and their impact on complex anatomical connections and resultant functional deficits.
Audiologic impairment associated with bilirubin-induced neurologic damage
Hyperbilirubinemia affects up to 84% of term and late preterm infants in the first week of life. The elevation of total serum/plasma bilirubin (TB) levels is generally mild, transitory, and, for most children, inconsequential. However, a subset of infants experiences lifelong neurological sequelae. Although the prevalence of classic kernicterus has fallen steadily in the USA in recent years, the incidence of jaundice in term and premature infants has increased, and kernicterus remains a significant problem in the global arena. Bilirubin-induced neurologic dysfunction (BIND) is a spectrum of neurological injury due to acute or sustained exposure of the central nervous system(CNS) to bilirubin. The BIND spectrum includes kernicterus, acute bilirubin encephalopathy, and isolated neural pathway dysfunction.
Animal studies have shown that unconjugated bilirubin passively diffuses across cell membranes and the blood‒brain barrier (BBB), and bilirubin not removed by organic anion efflux pumps accumulates within the cytoplasm and becomes toxic. Exposure of neurons to bilirubin results in increased oxidative stress and decreased neuronal proliferation and presynaptic neuro-degeneration at central glutaminergic synapses. Furthermore, bilirubin administration results in smaller spiral ganglion cell bodies, with decreased cellular density and selective loss of large cranial nerve VIII myelinated fibers. When exposed to bilirubin, neuronal supporting cells have been found to secrete inflammatory markers, which contribute to increased BBB permeability and bilirubin loading.
The jaundiced Gunn rat is the classic animal model of bilirubin toxicity. It is homozygous for a premature stop codon within the gene for UDP-glucuronosyltransferase family 1 (UGT1). The resultant gene product has reduced bilirubin-conjugating activity, leading to a state of hyperbilirubinemia. Studies with this rat model have led to the concept that impaired calcium homeostasis is an important mechanism of neuronal toxicity, with reduced expression of calcium-binding proteins in affected cells being a sensitive index of bilirubin-induced neurotoxicity. Similarly, application of bilirubin to cultured auditory neurons from brainstem cochlear nuclei results in hyperexcitability and excitotoxicity.
The auditory pathway and normal auditory brainstem response (ABR).
The auditory pathway and normal auditory brain-stem response (ABR). The ipsilateral (green) and contralateral (blue) auditory pathways are shown, with structures that are known to be affected by hyperbilirubinemia highlighted in red. Roman numerals in parentheses indicate corresponding waves in the normal human ABR (inset). Illustration adapted from the “Ear Anatomy” series by Robert Jackler and Christine Gralapp, with permission.
Beyond the traditional recognized areas of fulminant injury to the globus pallidus as seen in infants with kernicterus, other vulnerable areas include the cerebellum, hippocampus, and subthalamic nuclear bodies as well as certain cranial nerves. The hippocampus is a brain region that is particularly affected by age related morphological changes. It is generally assumed that a loss in hippocampal volume results in functional deficits that contribute to age-related cognitive deficits. Lower grey matter volumes within the limbic-striato-thalamic circuitry are common to other etiological mechanisms of subtle neurologic injury. Lower grey matter volumes in the amygdala, caudate, frontal and medial gyrus are found in schizophrenia and in the putamen in autism. Thus, in terms of brain volumetrics, schizophrenia and autism spectrum disorders have a clear degree of overlap that may reflect shared etiological mechanisms. Overlap with injuries observed in infants with BIND raises the question about how these lesions are arrived at in the context of the impact of common etiologies.
Stress-induced perinatal and transgenerational epigenetic programming of brain development and mental health
Research efforts during the past decades have provided intriguing evidence suggesting that stressful experiences during pregnancy exert long-term consequences on the future mental wellbeing of both the mother and her baby. Recent human epidemiological and animal studies indicate that stressful experiences in utero or during early life may increase the risk of neurological and psychiatric disorders, arguably via altered epigenetic regulation. Epigenetic mechanisms, such as miRNA expression, DNA methylation, and histone modifications are prone to changes in response to stressful experiences and hostile environmental factors. Altered epigenetic regulation may potentially influence fetal endocrine programming and brain development across several generations. Only recently, however, more attention has been paid to possible transgenerational effects of stress. In this review we discuss the evidence of transgenerational epigenetic inheritance of stress exposure in human studies and animal models. We highlight the complex interplay between prenatal stress exposure, associated changes in miRNA expression and DNA methylation in placenta and brain and possible links to greater risks of schizophrenia, attention deficit hyperactivity disorder, autism, anxiety- or depression-related disorders later in life. Based on existing evidence, we propose that prenatal stress, through the generation of epigenetic alterations, becomes one of the most powerful influences on mental health in later life. The consideration of ancestral and prenatal stress effects on lifetime health trajectories is critical for improving strategies that support healthy development and successful aging.
Sensitive time-windows for susceptibility in neurodevelopmental disorders
Many neurodevelopmental disorders (NDDs) are characterized by age-dependent symptom onset and regression, particularly during early postnatal periods of life. The neurobiological mechanisms preceding and underlying these developmental cognitive and behavioral impairments are, however, not clearly understood. Recent evidence using animal models for monogenic NDDs demonstrates the existence of time-regulated windows of neuronal and synaptic impairments. We propose that these developmentally-dependent impairments can be unified into a key concept: namely, time-restricted windows for impaired synaptic phenotypes exist in NDDs, akin to critical periods during normal sensory development in the brain. Existence of sensitive time-windows has significant implications for our understanding of early brain development underlying NDDs and may indicate vulnerable periods when the brain is more susceptible to current therapeutic treatments.
Fig (not shown)
Misregulated mechanisms underlying spine morphology in NDDs. Several proteins implicated in monogenic NDDs (highlighted in red) are linked to the regulation of the synaptic cytoskeleton via F-actin through different Rho-mediated signaling pathways (highlighted in green). Mutations in OPHN1, TSC1/2, FMRP, p21-activated kinase (PAK) are directly linked to human NDDs of intellectual disability. For instance, point mutations in OPHN1 and a PAK isoform are linked to non-syndromic mental retardation, whereas mutations or altered expression of TSC1/2 and FMRP are linked to TSC and FXS, respectively. Cytoplasmic interacting protein (CYFIP) and LIM-domain kinase 1 (LIMK1) are known to interact with FMRP and PAK, respectively [105]. LIMK1 is one of many dysregulated proteins contributing to the NDD Williams syndrome. Mouse models are available for all highlighted (red) proteins and reveal specific synaptic and behavioral deficits. Local protein synthesis in synapses, dendrites and glia is also regulated by proteins such as TSC1/2 and the FMRP/CYFIP complex. Abbreviations: 4EBP, 4E binding protein; eIF4E, eukaryotic translation initiation factor 4E.
Fig (not shown)
Sensitive time-windows, synaptic phenotypes and NDD gene targets. Sensitive time-windows exist in neural circuits, during which gene targets implicated in NDDs are normally expressed. Misregulation of these genes can affect multiple synaptic phenotypes during a restricted developmental period. The effect upon synaptic phenotypes is dependent upon the temporal expression of these NDD genes and their targets. (a) Expression outside a critical period of development will have no effect upon synaptic phenotypes. (b,c) A temporal expression pattern that overlaps with the onset (b) or closure (c) of a known critical period can alter the synaptic phenotype during that developmental time-window.
Outstanding questions
(1) Can treatment at early presymptomatic stages in animal models for NDDs prevent or ease the later synaptic, neuronal, and behavioral impairments?
(2) Are all sensory critical periods equally misregulated in mouse models for a specific NDD? Are there different susceptibilities for auditory, visual and somatosensory neurocircuits that reflect the degree of impairments observed in patients?
(3) If one critical period is missed or delayed during formation of a layer-specific connection in a network, does the network overcome this misregulated connectivity or plasticity window?
(4) In monogenic NDDs, does the severity of misregulating one particular time-window for synaptic establishment during development correlate with the importance of that gene for that synaptic circuit?
(5) Why do critical periods close in brain development?
(6) What underlies the regression of some altered synaptic phenotypes in Fmr1-KO mice?
(7) Can the concept of susceptible time-windows be applied to other NDDs, including schizophrenia and Tourette’s syndrome?
There is an increased interest in methods of objective cardiac output measurement in critically ill patients. Several techniques are available for measurement of cardiac output in children, although this remains very complex in newborns. Cardiac output monitoring could provide essential information to guide hemodynamic management. An overview is given of various methods of cardiac output monitoring with advantages and major limitations of each technology together with a short explanation of the basic principles.
Fick principle
According to the Fick principle the volume of blood flow in a given period equals the amount of substance entering the blood stream in the same period divided by the difference in concentrations of the substrate upstream respectively downstream to the point of entry in the circulation. This substance can be oxygen (O2-Fick) or carbon dioxide (CO2-FICK), so cardiac output can be calculated by dividing measured pulmonary oxygen uptake by the arteriovenous oxygen concentration difference. The direct O2-Fick method is regarded as gold standard in cardiac output monitoring in a research setting, despite its limitations. When the Fick principle is applied for carbon dioxide (CO2 Fick), the pulmonary carbon dioxide exchange is divided by the venoarterial CO2 concentration difference to calculate cardiac output.
In the modified CO2 Fick method pulmonary CO2 exchange is measured at the endotracheal tube. Measurement of total CO2 concentration in blood is more complex and simultaneous sampling of arterial and central venous blood is required. However, frequent blood sampling will result in an unacceptable blood loss in the neonatal population.
Blood flow can be calculated if the change in concentration of a known quantity of injected indicator is measured in time distal to the point of injection, so an indicator dilution curve can be obtained. Cardiac output can then be calculated with the use of the Stewart–Hamilton equation. Several indicators are used, such as indocyanine green, Evans blue and brilliant red in dye dilution, cold solutions in thermodilution, lithium in lithium dilution, and isotonic saline in ultrasound dilution.
Cardiovascular adaptation to extra uterine life
Alice Lawford, Robert MR Tulloh
Paediatrics And Child Health 2014; 25(1): 1-6.
The adaptation to extra uterine life is of interest because of its complexity and the ability to cause significant health concerns. In this article we describe the normal changes that occur and the commoner abnormalities that are due to failure of normal development and the effect of congenital cardiac disease. Abnormal development may occur as a result of problems with the mother, or with the fetus before birth. After birth it is essential to determine whether there is an underlying abnormality of the fetal pulmonary or cardiac development and to determine the best course of management of pulmonary hypertension or congenital cardiac disease. Causes of underdevelopment, maldevelopment and maladaptation are described as are the causes of critical congenital heart disease. The methods of diagnosis and management are described to allow the neonatologist to successfully manage such newborns.
Fetal vascular structures that exist to direct blood flow
Fetal structure
Function
Arterial duct
Connects pulmonary artery to the aorta and shunts blood right to left; diverting flow away from fetal lungs
Foramen ovale
Opening between the two atria thatdirects blood flow returning to right
atrium through the septal wall into the left atrium bypassing lungs
Ductus venosus
Receives oxygenated blood fromumbilical vein and directs it to the
inferior vena cava and right atrium
Umbilical arteries
Carrying deoxygenated blood fromthe fetus to the placenta
Umbilical vein
Carrying oxygenated blood from theplacenta to the fetus
Maternal causes of congenital heart disease
Maternal disorders
rubella, SLE, diabetes mellitus
Maternal drug use
Warfarin, alcohol
Chromosomal abnormality
Down, Edward, Patau, Turner, William, Noonan
Fetal and Neonatal Circulation The fetal circulation is specifically adapted to efficiently exchange gases, nutrients, and wastes through placental circulation. Upon birth, the shunts (foramen ovale, ductus arteriosus, and ductus venosus) close and the placental circulation is disrupted, producing the series circulation of blood through the lungs, left atrium, left ventricle, systemic circulation, right heart, and back to the lungs.
Clinical monitoring of systemic hemodynamics in critically ill newborns
Circulatory failure is a major cause of mortality and morbidity in critically ill newborn infants. Since objective measurement of systemic blood flow remains very challenging, neonatal hemodynamics is usually assessed by the interpretation of various clinical and biochemical parameters. An overview is given about the predictive value of the most used indicators of circulatory failure, which are blood pressure, heart rate, urine output, capillary refill time, serum lactate concentration, central–peripheral temperature difference, pH, standard base excess, central venous oxygen saturation and color.
Key guidelines
➢ The clinical assessment of cardiac output by the interpretation of indirect parameters of systemic blood flow is inaccurate, irrespective of the level of experience of the clinician
➢ Using blood pressure to diagnose low systemic blood flow will consequently mean that too many patients will potentially be undertreated or overtreated, both with substantial risk of adverse effects and iatrogenic damage.
➢ Combining different clinical hemodynamic parameters enhances the predictive value in the detection of circulatory failure, although accuracy is still limited.
➢ Variation in time (trend monitoring) might possibly be more informative than individual, static values of clinical and biochemical parameters to evaluate the adequacy of neonatal circulation.
Monitoring oxygen saturation and heart rate in the early neonatal period
Pulse oximetry is commonly used to assist clinicians in assessment and management of newly born infants in the delivery room (DR). In many DRs, pulse oximetry is now the standard of care for managing high risk infants, enabling immediate and dynamic assessment of oxygenation and heart rate. However, there is little evidence that using pulse oximetry in the DR improves short and long term outcomes. We review the current literature on using pulse oximetry to measure oxygen saturation and heart rate and how to apply current evidence to management in the DR.
Practice points
Understand how SpO2 changes in the first minutes after birth.
Apply a sensor to an infant’s right wrist as soon as possible after birth.
Attach sensor to infant then to oximeter cable.
Use two second averaging and maximum sensitivity.
Using pulse oximetry assists clinicians:
Assess changes in HR in real time during transition.
Assess oxygenation and titrate the administration of oxygen to maintain oxygenation within the appropriate range for SpO2 during the first minutes after birth.
Research directions
What are the appropriate centiles to target during the minutes after birth to prevent hypoxia and hyperoxia: 25th to 75th, or 10th to 90th, or just the 50th (median)?
Can the inspired oxygen be titrated against the SpO2 to keep the SpO2 in the ‘normal range’?
Does the use of centile charts in the DR for HR and oxygen saturation reduce the rate of hyperoxia when infants are treated with oxygen.
Does the use of pulse oximetry immediately after birth improve short term outcomes, e.g. efficacy of immediate respiratory support, intubation rates in the DR, percentage of inspired oxygen, rate of use of adrenalin or chest compressions, duration of hypoxia/hyperoxia and bradycardia.
Does the use of pulse oximetry in the DR improve short term respiratory and long term neurodevelopmental outcomes for preterm infants, e.g. rate of intubation, use of surfactant, and duration of ventilation, continuous positive airway pressure, or supplemental oxygen?
Can all modern pulse oximeters be used effectively in the DR or do some have a longer delay before giving an accurate signal and more movement artefact?
Would a longer averaging time result in more stable data?
Peripheral haemodynamics in newborns: Best practice guidelines
Peripheral hemodynamics refers to blood flow, which determines oxygen and nutrient delivery to the tissues. Peripheral blood flow is affected by vascular resistance and blood pressure, which in turn varies with cardiac function. Arterial oxygen content depends on the blood hemoglobin concentration (Hb) and arterial pO2; tissue oxygen delivery depends on the position of the oxygen-dissociation curve, which is determined by temperature and the amount of adult or fetal hemoglobin. Methods available to study tissue perfusion include near-infrared spectroscopy, Doppler flowmetry, orthogonal polarization spectral imaging and the peripheral perfusion index. Cardiac function, blood gases, Hb, and peripheral temperature all affect blood flow and oxygen extraction. Blood pressure appears to be less important. Other factors likely to play a role are the administration of vasoactive medications and ventilation strategies, which affect blood gases and cardiac output by changing the intrathoracic pressure.
graphic
NIRS with partial venous occlusion to measure venous oxygen saturation
NIRS with partial venous occlusion to measure venous oxygen saturation. Taken from Yoxall and Weindling
Schematic representation of the biphasic relationship between oxygen delivery and oxygen consumption in tissue
graphic
Schematic representation of the biphasic relationship between oxygen delivery and oxygen consumption in tissue. (a) oxygen delivery (DO2). (b) As DO2 decreases, VO2 is dependent on DO2. The slope of the line indicates the FOE, which in this case is about 0.50. (c) The slope of the line indicates the FOE in the normal situation where oxygenation is DO2 independent, usually < 0.35
The oxygen-dissociation curve
graphic
The oxygen-dissociation curve
Considerable information about the response of the peripheral circulation has been obtained using NIRS with venous occlusion. Although these measurements were validated against blood co-oximetry in human adults and infants, they can only be made intermittently by a trained operator and are thus not appropriate for general clinical use. Further research is needed to find other better measures of peripheral perfusion and oxygenation which may be easily and continuously monitored, and which could be useful in a clinical setting.
Peripheral oxygenation and management in the perinatal period
The mechanisms for the adequate provision of oxygen to the peripheral tissues are complex. They involve control of the microcirculation and peripheral blood flow, the position of the oxygen dissociation curve including the proportion of fetal and adult hemoglobin, blood gases and viscosity. Systemic blood pressure appears to have little effect, at least in the non-shocked state. The adequate delivery of oxygen (DO2) depends on consumption (VO2), which is variable. The balance between VO2 and DO2 is given by fractional oxygen extraction (FOE ¼ VO2/DO2). FOE varies from organ to organ and with levels of activity. Measurements of FOE for the whole body produce a range of about 0.15-0.33, i.e. the body consumes 15-33% of oxygen transported.
Fig (not shown)
Biphasic relationship between oxygen delivery (DO2) and oxygen consumption (VO2) in tissue. Dotted lines show fractional oxygen extraction (FOE). ‘A’ indicates the normal situation when VO2 is independent ofDO2 and FOE is about 0.30. AsDO2 decreases in the direction of the arrow, VO2 remains independent of DO2 until the critical point is reached at ‘B’; in this illustration, FOE is about 0.50. The slope of the dotted line indicates the FOE (¼ VO2/DO2), which increases progressively as DO2 decreases.
Relationship between haemoglobin F fraction (HbF) and peripheral fractional oxygen extraction
Graphic
(A)Relationship between haemoglobin F fraction (HbF) and peripheral fractional oxygen extraction in anaemic and control infants. (From Wardle et al.) (B) HbF synthesis and concentration. (From Bard and Widness.) (C) Oxygen dissociation curve.
Peripheral fractional oxygen extraction in babies
graphic
Peripheral fractional oxygen extraction in babies with asymptomatic or symptomatic anemia compared to controls. Bars represent the median for each group. (From Wardle et al.)
Practice points
Peripheral tissue DO2 is complex: cardiac function, blood gases, Hb concentration and the proportion of HbF, and peripheral temperature all play a part in determining blood flow and oxygen extraction in the sick, preterm infant. Blood pressure appears to be less important.
Other factors likely to play a role are the administration of vasoactive medications and ventilation strategies, which affect blood gases and cardiac output by changing intrathoracic pressure.
Central blood pressure is a poor surrogate measurement for the adequacy of DO2 to the periphery. Direct measurement, using NIRS, laser Doppler flowmetry or other means, may give more useful information.
Reasons for total hemoglobin concentration (Hb) being a relatively poor indicator of the adequacy of the provision of oxygen to the tissues:
Hb is only indirectly related to red blood cell volume, which may be a better indicator of the body’s oxygen delivering capacity.
Hb-dependent oxygen availability depends on the position of the oxygen-hemoglobin dissociation curve.
An individual’s oxygen requirements vary with time and from organ to organ. This means that DO2 also needs to vary.
It is possible to compensate for a low Hb by increasing cardiac output and ventilation, and so the ability to compensate for anemia depends on an individual’s cardio-respiratory reserve as well as Hb.
The normal decrease of Hb during the first few weeks of life in both full-term and preterm babies usually occurs without symptoms or signs of anemia or clinical consequences.
The relationship between VO2 and DO2 is complex and various factors need to be taken into account, including the position of the oxygen dissociation curve, determined by the proportion of HbA and HbF, temperature and pH. Furthermore, diffusion of oxygen from capillaries to the cell depends on the oxygen tension gradient between erythrocytes and the mitochondria, which depends on microcirculatory conditions, e.g. capillary PO2, distance of the cell from the capillary (characterized by intercapillary distances) and the surface area of open capillaries. The latter can change rapidly, for example, in septic shock where arteriovenous shunting occurs associated with tissue hypoxia in spite of high DO2 and a low FOE.
Changes in local temperature deserve particular consideration. When the blood pressure is low, there may be peripheral vasoconstriction with decreased local perfusion and DO2. However, the fall in local tissue temperature would also be expected to be associated with a decreased metabolic rate and a consequent decrease in VO2. Thus a decreased DO2 may still be appropriate for tissue needs.
Pulmonary
Accurate Measurements of Oxygen Saturation in Neonates: Paired Arterial and Venous Blood Analyses
Oxygen saturation (So2) measurements (functional measurement, So2; and fractional measurement, oxyhemoglobin [Hbo2]) and monitoring are commonly investigated as a method of assessing oxygenation in neonates. Differences exist between the So2 and Hbo2 when blood tests are performed, and clinical monitors indicate So2 values. Oxyhemoglobin will decrease with the increased levels of carbon monoxide hemoglobin (Hbco) and methemo-globin (MetHb), and it is the most accurate measurements of oxygen (O2) association of hemoglobin (Hb). Pulse oximeter (for pulse oximetry saturation [Spo2] measurement) is commonly used in neonates. However, it will not detect the changes of Hb variations in the blood for accurate So2 measurements. Thus, the measurements from clinical oximeters should be used with caution. In neonates, fetal hemoglobin (HbF) accounts for most of the circulating Hb in their blood. Fetal hemoglobin has a high O2 affinity, thus releases less O2 to the body tissues, presenting a left-shifted Hbo2 dissociation curve.5,6 To date, however, limited data are available with HbF correction, for accurate arterial and venous (AV) So2 measurements (arterial oxygen saturation [Sao2] and venous oxygen saturation [Svo2]) in neonates, using paired AV blood samples.
In a study of critically ill adult patients, increased pulmonary CO production and elevation in arterial Hbco but not venous Hbco were documented by inflammatory stimuli inducing pulmonary heme oxygenase–1. In normal adults, venous Hbco level might be slightly higher than or equal to arterial Hbco because of production of CO by enzyme heme oxygenase–2, which is predominantly produced in the liver and spleen. However, hypoxia or pulmonary inflammation could induce heme oxygenase–1 to increase endogenous CO, thus elevating pulmonary arterial and systemic arterial Hbco levels in adults. Both endogenous and exogenous CO can suppress proliferation of pulmonary smooth muscles, a significant consideration for the prevention of chronic lung diseases in newborns. Despite these considerations, a later study in healthy adults indicated that the AV differences in Hbco were from technical artifacts and perhaps from inadequate control of different instruments. Thus, further studies are needed to provide more definitive answers for the AV differences of Hbco for adults and neonates with acute and chronic lung diseases.
Methemoglobin is an indicator of Hb oxidation and is essential for accurate measurement of Hbo2, So2, and oxygenation status. No evidence exists to show the AV MetHb difference, although this difference was elucidated with the potential changes of MetHb with different O2 levels. Methemoglobin can be increased with nitric oxide (NO) therapy, used in respiratory distress syndrome (RDS) to reduce pulmonary hypertension and during heart surgery. Nitric oxide, in vitro, is an oxidant of Hb, with increased O2 during ischemia reperfusion. In hypoxemic conditions in vivo, nitrohemoglobin is a product generated by vessel responsiveness to nitrovasodilators. Nitro-hemoglobin can be spontaneously reversible in vivo, requiring no chemical agents or reductase. However, when O2 levels were increased experimentally in vitro following acidic conditions (pH 6.5) to simulate reperfusion conditions, MetHb levels were increased for the hemolysates (broken red cells). Nitrite-induced oxidation of Hb was associated with an increase in red blood cell membrane rigidity, thus contributing to Hb breakdown. A newer in vitro study of whole blood cells, however, concluded that MetHb formation is not dependent on increased O2 levels. Additional studies are needed to examine in vivo reperfusion of O2 and MetHb effects.
Purpose: The aim of this study was to examine the accuracy of arterial oxygen saturation (Sao2) and venous oxygen saturation (Svo2) with paired arterial and venous (AV) blood in relation to pulse oximetry saturation (Spo2) and oxyhemoglobin (Hbo2) with fetal hemoglobin determination, and their Hbo2 dissociation curves. Method: Twelve preterm neonates with gestational ages ranging from 27 to 34 weeks at birth, who had umbilical AV lines inserted, were investigated. Analyses were performed with 37 pairs of AV blood samples by using a blood volume safety protocol. Results: The mean differences between Sao2 and Svo2, and AV Hbo2 were both 6 percent (F6.9 and F6.7 percent, respectively), with higher Svo2 than those reported for adults. Biases were 2.1 – 0.49 for Sao2, 2.0 – 0.44 for Svo2, and 3.1 – 0.45 for Spo2, compared against Hbo2. With left-shifted Hbo2 dissociation curves in neonates, for the critical values of oxygen tension values between 50 and 75 millimeters of mercury, Hbo2 ranged from 92 to 93.4 percent; Sao2 ranged from 94.5 to 95.7 percent; and Spo2 ranged from 93.7 to 96.3 percent (compared to 85–94 percent in healthy adults). Conclusions: In neonates, both left-shifted Hbo2 dissociation curve and lower AV differences of oxygen saturation measurements indicated low flow of oxygen to the body tissues. These findings demonstrate the importance of accurate assessment of oxygenation statues in neonates.
In these neonates, the mean AV blood differences for both So2 and Hbo2 were about 6 percent, which was much lower than those reported for healthy adults (23 percent) for O2 supply and demand. In addition, with very high levels of HbF releasing less O2 to the body tissue, the results of blood analyses are worrisome for these critically ill neonates for low systemic oxygen states. O’Connor and Hall determined AV So2 in neonates without HbF determination. Much of the AV So2 difference is dependent on Svo2 measurement. The ranges of Svo2 spanned for 35 percent, and the ranges of Sao2 spanned 6 percent in these neonates. The greater intervals for Svo2 measurements contribute to greater sensitivity for the measurements (than Sao2 measurements) in responding to nursing care and changes of O2 demand. Thus, Svo2 measurement is essential for better assessment of oxygenation status in neonates.
The findings of this study on AV differences of So2 were limited with very small number of paired AV blood samples. However, critically ill neonates need accurate assessment of oxygenation status because of HbF, which releases less O2 to the tissues. Decreased differences of AV So2 measurements added further possibilities of lower flow of O2 to the body tissues and demonstrated the greater need to accurately assess the proper oxygenation in the neonates. The findings of this study continued to clarify the accuracy of So2 measurements for neonates. Additional studies are needed to examine So2 levels in neonates to further validate these findings by using larger sample sizes.
Neonatal ventilation strategies and long-term respiratory outcomes
Long-term respiratory morbidity is common, particularly in those born very prematurely and who have developed bronchopulmonary dysplasia (BPD), but it does occur in those without BPD and in infants born at term. A variety of neonatal strategies have been developed, all with short-term advantages, but meta-analyses of randomized controlled trials (RCTs) have demonstrated that only volume-targeted ventilation and prophylactic high-frequency oscillatory ventilation (HFOV) may reduce BPD. Few RCTs have incorporated long-term follow-up, but one has demonstrated that prophylactic HFOV improves respiratory and functional outcomes at school age, despite not reducing BPD. Results from other neonatal interventions have demonstrated that any impact on BPD may not translate into changes in long-term outcomes. All future neonatal ventilation RCTs should have long-term outcomes rather than BPD as their primary outcome if they are to impact on clinical practice.
A Model Analysis of Arterial Oxygen Desaturation during Apnea in Preterm Infants
Rapid arterial O2 desaturation during apnea in the preterm infant has obvious clinical implications but to date no adequate explanation for why it exists. Understanding the factors influencing the rate of arterial O2 desaturation during apnea (_SSaO2 ) is complicated by the non-linear O2 dissociation curve, falling pulmonary O2 uptake, and by the fact that O2 desaturation is biphasic, exhibiting a rapid phase (stage 1) followed by a slower phase when severe desaturation develops (stage 2). Using a mathematical model incorporating pulmonary uptake dynamics, we found that elevated metabolic O2 consumption accelerates _SSaO2 throughout the entire desaturation process. By contrast, the remaining factors have a restricted temporal influence: low pre-apneic alveolar PO2 causes an early onset of desaturation, but thereafter has little impact; reduced lung volume, hemoglobin content or cardiac output, accelerates _SSaO2 during stage 1, and finally, total blood O2 capacity (blood volume and hemoglobin content) alone determines _SSaO2 during stage 2. Preterm infants with elevated metabolic rate, respiratory depression, low lung volume, impaired cardiac reserve, anemia, or hypovolemia, are at risk for rapid and profound apneic hypoxemia. Our insights provide a basic physiological framework that may guide clinical interpretation and design of interventions for preventing sudden apneic hypoxemia.
A novel approach to study oxidative stress in neonatal respiratory distress syndrome
Oxidative stress is an imbalance between the systemic manifestation of reactive oxygen species and a biological system’s ability to readily detoxify the reactive intermediates or to repair the resulting damage. It is a physiological event in the fetal-to-neonatal transition, which is actually a great stress to the fetus. These physiological changes and processes greatly increase the production of free radicals, which must be controlled by the antioxidant defense system, the maturation of which follows the course of the gestation. This could lead to several functional alterations with important repercussions for the infants. Adequately mature and healthy infants are able to tolerate this drastic change in the oxygen concentration. A problem occurs when the intrauterine development is incomplete or abnormal. Preterm or intrauterine growth retarded (IUGR) and low birth weight neonates are typically of this kind. An oxidant/antioxidant imbalance in infants is implicated in the pathogenesis of the major complications of prematurity including respiratory distress syndrome (RDS), necrotizing enterocolitis (NEC), chronic lung disease, retinopathy of prematurity and intraventricular hemorrhage (IVH).
Background: Respiratory distress syndrome of the neonate (neonatal RDS) is still an important problem in treatment of preterm infants. It is accompanied by inflammatory processes with free radical generation and oxidative stress. The aim of study was to determine the role of oxidative stress in the development of neonatal RDS. Methods: Markers of oxidative stress and antioxidant activity in umbilical cord blood were studied in infants with neonatal respiratory distress syndrome with reference to healthy newborns. Results: Status of markers of oxidative stress (malondialdehyde, protein carbonyl and 8-hydroxy-2-deoxy guanosine) showed a significant increase with depleted levels of total antioxidant capacity in neonatal RDS when compared to healthy newborns. Conclusion: The study provides convincing evidence of oxidative damage and diminished antioxidant defenses in newborns with RDS. Neonatal RDS is characterized by damage of lipid, protein and DNA, which indicates the augmentation of oxidative stress. General significance: The identification of the potential biomarker of oxidative stress consists of a promising strategy to study the pathophysiology of neonatal RDS.
Neonatal respiratory distress syndrome represents the major lung complications of newborn babies. Preterm neonates suffer from respiratory distress syndrome (RDS) due to immature lungs and require assisted ventilation with high concentrations of oxygen. The pathogenesis of this disorder is based on the rapid formation of the oxygen reactive species, which surpasses the detoxification capacity of antioxidative defense system. The high chemical reactivity of free radical leads to damage to a variety of cellular macro molecules including proteins, lipids and nucleic acid. This results in cell injury and may induce respiratory cell death.
Malondialdehyde (MDA) is one of the final products of polyunsaturated fatty acids peroxidation. The present study showed increased concentration of MDA in neonates with respiratory disorders than that of control in consonance with the reported study.
Anemia, Apnea of Prematurity, and Blood Transfusions
The etiology of apnea of prematurity is multifactorial; however, decreased oxygen carrying capacity may play a role. The respiratory neuronal network in neonates is immature, particularly in those born preterm, as demonstrated by their paradoxical response to hypoxemia. Although adults increase the minute ventilation in response to hypoxemia, newborns have a brief increase in ventilation followed by periodic breathing, respiratory depression, and occasionally cessation of respiratory effort. This phenomenon may be exacerbated by anemia in preterm newborns, where a decreased oxygen carrying capacity may result in decreased oxygen delivery to the central nervous system, a decreased efferent output of the respiratory neuronal network, and an increase in apnea.
Objective Compare the frequency and severity of apneic events in very low birth weight (VLBW) infants before and after blood transfusions using continuous electronic waveform analysis. Study design We continuously collected waveform, heart rate, and oxygen saturation data from patients in all 45 neonatal intensive care unit beds at the University of Virginia for 120 weeks. Central apneas were detected using continuous computer processing of chest impedance, electrocardiographic, and oximetry signals. Apnea was defined as respiratory pauses of >10, >20, and >30 seconds when accompanied by bradycardia (<100 beats per minute) and hypoxemia (<80% oxyhemoglobin saturation as detected by pulse oximetry). Times of packed red blood cell transfusions were determined from bedside charts. Two cohorts were analyzed. In the transfusion cohort, waveforms were analyzed for 3 days before and after the transfusion for all VLBW infants who received a blood transfusion while also breathing spontaneously. Mean apnea rates for the previous 12 hours were quantified and differences for 12 hours before and after transfusion were compared. In the hematocrit cohort, 1453 hematocrit values from all VLBW infants admitted and breathing spontaneously during the time period were retrieved, and the association of hematocrit and apnea in the next 12 hours was tested using logistic regression. Results Sixty-seven infants had 110 blood transfusions during times when complete monitoring data were available. Transfusion was associated with fewer computer-detected apneic events (P < .01). Probability of future apnea occurring within 12 hours increased with decreasing hematocrit values (P < .001). Conclusions Blood transfusions are associated with decreased apnea in VLBW infants, and apneas are less frequent at higher hematocrits.
Bronchopulmonary dysplasia: The earliest and perhaps the longest lasting obstructive lung disease in humans
Bronchopulmonary dysplasia (BPD) is one of the most important sequelae of premature birth and the most common form of chronic lung disease of infancy, an umbrella term for a number of different diseases that evolve as a consequence of a neonatal respiratory disorder. BPD is defined as the need for supplemental oxygen for at least 28 days after birth, and its severity is graded according to the respiratory support required at 36 post-menstrual weeks.
BPD was initially described as a chronic respiratory disease occurring in premature infants exposed to mechanical ventilation and oxygen supplementation. This respiratory disease (later named “old BPD”) occurred in relatively large premature newborn and, from a pathological standpoint, it was characterized by intense airway inflammation, disruption of normal pulmonary structures and lung fibrosis.
Bronchopulmonary dysplasia (BPD) is one of the most important sequelae of premature birth and the most common form of chronic lung disease of infancy. From a clinical standpoint BPD subjects are characterized by recurrent respiratory symptoms, which are very frequent during the first years of life and, although becoming less severe as children grow up, they remain more common than in term-born controls throughout childhood, adolescence and into adulthood. From a functional point of view BPD subjects show a significant airflow limitation that persists during adolescence and adulthood and they may experience an earlier and steeper decline in lung function during adulthood. Interestingly, patients born prematurely but not developing BPD usually fare better, but they too have airflow limitations during childhood and later on, suggesting that also prematurity per se has life-long detrimental effects on pulmonary function. For the time being, little is known about the presence and nature of pathological mechanisms underlying the clinical and functional picture presented by BPD survivors. Nonetheless, recent data suggest the presence of persistent neutrophilic airway inflammation and oxidative stress and it has been suggested that BPD may be sustained in the long term by inflammatory pathogenic mechanisms similar to those underlying COPD. This hypothesis is intriguing but more pathological data are needed. A better understanding of these pathogenetic mechanisms, in fact, may be able to orient the development of novel targeted therapies or prevention strategies to improve the overall respiratory health of BPD patients.
We have a limited understanding of the presence and nature of pathological mechanisms in the lung of BPD survivors. The possible role of asthma-like inflammation has been investigated because BPD subjects often present with recurrent wheezing and other symptoms resembling asthma during their childhood and adolescence. But BPD subjects have normal or lower than normal exhaled nitric oxide levels and exhaled air temperatures, whereas they are higher than normal in asthmatic patients.
Of all obstructive lung diseases in humans, BPD has the earliest onset and is possibly the longest lasting. Given its frequent association with other conditions related to preterm birth (e.g. growth retardation, pulmonary hypertension, neurodevelopmental delay, hearing defects, and retinopathy of prematurity), it often warrants a multidisciplinary management.
Effects of Sustained Lung Inflation, a lung recruitment maneuver in primary acute respiratory distress syndrome, in respiratory and cerebral outcomes in preterm infants
Background: Sustained Lung Inflation (SLI) is a maneuver of lung recruitment in preterm newborns at birth that can facilitate the achieving of larger inflation volumes, leading to the clearance of lung fluid and formation of functional residual capacity (FRC). Aim: To investigate if Sustained Lung Inflation (SLI) reduces the need of invasive procedures and iatrogenic risks. Study design: 78 newborns (gestational age ≤ 34 weeks, weighing ≤ 2000 g) who didn’t breathe adequately at birth and needed to receive SLI in addition to other resuscitation maneuvers (2010 guidelines). Subjects: 78 preterm infants born one after the other in our department of Neonatology of Catania University from 2010 to 2012. Outcome measures: The need of intubation and surfactant, the ventilation required, radiological signs, the incidence of intraventricular hemorrhage (IVH), periventricular leukomalacia, retinopathy in prematurity from III to IV plus grades, bronchopulmonary dysplasia, patent ductus arteriosus, pneumothorax and necrotizing enterocolitis. Results: In the SLI group infants needed less intubation in the delivery room (6% vs 21%; p b 0.01), less invasive mechanical ventilation (14% vs 55%; p ≤ 0.001) and shorter duration of ventilation (9.1 days vs 13.8 days; p ≤ 0.001). There wasn’t any difference for nasal continuous positive airway pressure (82% vs 77%; p = 0.43); but there was less surfactant administration (54% vs 85%; p ≤ 0.001) and more infants received INSURE (40% vs 29%; p=0.17). We didn’t found any differences in the outcomes, except for more mild intraventricular hemorrhage in the SLI group (23% vs 14%; p = 0.15; OR= 1.83). Conclusion: SLI is easier to perform even with a single operator, it reduces the necessity of more complicated maneuvers and surfactant without statistically evident adverse effects.
Long-term respiratory consequences of premature birth at less than 32 weeks of gestation
Chronic respiratory morbidity is a common adverse outcome of very premature birth, particularly in infants who had developed bronchopulmonary dysplasia (BPD). Prematurely born infants who had BPD may require supplementary oxygen at home for many months and affected infants have increased healthcare utilization until school age. Chest radiograph abnormalities are common; computed tomography of the chest gives predictive information in children with ongoing respiratory problems. Readmission to hospital is common, particularly for those who have BPD and suffer respiratory syncytial virus lower respiratory infections (RSV LRTIs). Recurrent respiratory symptoms requiring treatment are common and are associated with evidence of airways obstruction and gas trapping. Pulmonary function improves with increasing age, but children with BPD may have ongoing airflow limitation. Lung function abnormalities may be more severe in those who had RSV LRTIs, although this may partly be explained by worse premorbid lung function. Worryingly, lung function may deteriorate during the first year. Longitudinal studies are required to determine if there is catch up growth.
Long-term pulmonary outcomes of patients with bronchopulmonary dysplasia
Bronchopulmonary dysplasia (BPD) is the commonest cause of chronic lung disease in infancy. The incidence of BPD has remained unchanged despite many advances in neonatal care. BPD starts in the neonatal period but its effects can persist long term. Premature infants with BPD have a greater incidence of hospitalization, and continue to have a greater respiratory morbidity and need for respiratory medications, compared to those without BPD. Lung function abnormalities, especially small airway abnormalities, often persist. Even in the absence of clinical symptoms, BPD survivors have persistent radiological abnormalities and presence of emphysema has been reported on chest computed tomography scans. Concern regarding their exercise tolerance remains. Long-term effects of BPD are still unknown, but given reports of a more rapid decline in lung function and their susceptibility to develop chronic obstructive pulmonary disease phenotype with aging, it is imperative that lung function of survivors of BPD be closely monitored.
Neonatal ventilation strategies and long-term respiratory outcomes
Long-term respiratory morbidity is common, particularly in those born very prematurely and who have developed bronchopulmonary dysplasia (BPD), but it does occur in those without BPD and in infants born at term. A variety of neonatal strategies have been developed, all with short-term advantages, but meta-analyses of randomized controlled trials (RCTs) have demonstrated that only volume-targeted ventilation and prophylactic high-frequency oscillatory ventilation (HFOV) may reduce BPD. Few RCTs have incorporated long-term follow-up, but one has demonstrated that prophylactic HFOV improves respiratory and functional outcomes at school age, despite not reducing BPD. Results from other neonatal interventions have demonstrated that any impact on BPD may not translate into changes in long-term outcomes. All future neonatal ventilation RCTs should have long-term outcomes rather than BPD as their primary outcome if they are to impact on clinical practice.
Prediction of neonatal respiratory distress syndrome in term pregnancies by assessment of fetal lung volume and pulmonary artery resistance index
Mohamed Laban, GM Mansour, MSE Elsafty, AS Hassanin, SS EzzElarab
International Journal of Gynecology and Obstetrics 128 (2015) 246–250 http://dx.doi.org/10.1016/j.ijgo.2014.09.018
Objective: To develop reference cutoff values for mean fetal lung volume (FLV) and pulmonary artery resistance index (PA-RI) for prediction of neonatal respiratory distress syndrome (RDS) in low-risk term pregnancies. Methods: As part of a cross-sectional study, women aged 20–35 years were enrolled and admitted to a tertiary hospital in Cairo, Egypt, for elective repeat cesarean at 37–40 weeks of pregnancy between January 1, 2012, and July 31, 2013. FLV was calculated by virtual organ computer-aided analysis, and PA-RI was measured by Doppler ultrasonography before delivery. Results: A total of 80 women were enrolled. Neonatal RDS developed in 11 (13.8%) of the 80 newborns. Compared with neonates with RDS, healthy neonates had significantly higher FLVs (P b 0.001) and lower PA-RIs (P b 0.001). Neonatal RDS is less likely with FLV of at least 32 cm3 or PA-RI less than or equal to 0.74. Combining these two measures improved the accuracy of prediction. Conclusion: The use of either FLV or PA-RI predicted neonatal RDS. The predictive value increased when these two measures were combined
Pulmonary surfactant – a front line of lung host defense, 2003 JCI0318650.f2
Pulmonary hypertension in bronchopulmonary dysplasia
Pulmonary hypertension (PH) is a common complication of neonatal respiratory diseases, including bronchopulmonary dysplasia (BPD), and recent studies have increased aware- ness that PH worsens the clinical course, morbidity and mortality of BPD. Recent evidence indicates that up to 18% of all extremely low-birth-weight infants will develop some degree of PH during their hospitalization, and the incidence rises to 25–40% of the infants with established BPD. Risk factors are not yet well understood, but new evidence shows that fetal growth restriction is a significant predictor of PH. Echocardiography remains the primary method for evaluation of BPD-associated PH, and the development of standardized screening timelines and techniques for identification of infants with BPD-associated PH remains an important ongoing topic of investigation. The use of pulmonary vasodilator medications, such as nitric oxide, sildenafil, and others, in the BPD population is steadily growing, but additional studies are needed regarding their long-term safety and efficacy. An update on pharmacologic approaches to bronchopulmonary dysplasia
Bronchopulmonary dysplasia (BPD) is the most prevalent long-term morbidity in surviving extremely preterm infants and is linked to increased risk of reactive airways disease, pulmonary hypertension, post-neonatal mortality, and adverse neurodevelopmental outcomes. BPD affects approximately 20% of premature newborns, and up to 60% of premature infants born before completing 26 weeks of gestation. It is characterized by the need for assisted ventilation and/or supplemental oxygen at 36 weeks postmenstrual age. Approaches to prevention and treatment of BPD have evolved with improved understanding of its pathogenesis. This review will focus on recent advancements and detail current research in pharmacotherapy for BPD. The evidence for both current and potential future experimental therapies will be reviewed in detail. As our understanding of the complex and multifactorial pathophysiology of BPD changes, research into these current and future approaches must continue to evolve.
Methylxanthines
Diuretics and bronchodilators
Corticosteroids
Macrolide antibiotics
Recombinant human Clara cell 10-kilodalton protein(rhCC10)
The nature, biological characteristics, and contribution to organ physiology of skeletal stem cells are not completely determined. Chan et al. and Worthley et al. demonstrate that a stem cell for skeletal tissues, and a system of more restricted, downstream progenitors, can be identified in mice and demonstrate its role in skeletal tissue maintenance and regeneration.
The groundbreaking concept that bone, cartilage, marrow adipocytes, and hematopoiesis-supporting stroma could originate from a common progenitor and putative stem cell was surprising at the time when it was formulated (Owen and Friedenstein, 1988). The putative stem cell, nonhematopoietic in nature, would be found in the postnatal bone marrow stroma, generate tissues previously thought of as foreign to each other, and support the turnover of tissues and organs that self-renew at a much slower rate compared to other tissues associated with stem cells (blood, epithelia). This concept also connected bone and bone marrow as parts of a single-organ system, implying their functional interplay. For many years, the evidence underpinning the concept has been incomplete.
While multipotency of stromal progenitors has been demonstrated by in vivo transplantation experiments, self-renewal, the defining property of a stem cell, has not been easily demonstrated until recently in humans (Sacchetti et al., 2007) and mice (Mendez-Ferrer et al., 2010). Meanwhile, a confusing and plethoric terminology has been introduced into the literature, which diverted and confounded the search for a skeletal stem cell and its physiological significance (Bianco et al., 2013).
Two studies in this issue of Cell (Chan et al., 2015; Worthley et al., 2015), using a combination of rigorous single-cell analyses and lineage tracing technologies, mark significant steps toward rectifying the course of skeletal stem cell discovery by making several important points, within and beyond skeletal physiology.
First, a stem cell for skeletal tissues, and a system of more restricted, downstream progenitors can in fact be identified and linked to defined phenotype(s) in the mouse. The system is framed conceptually, and approached experimentally, similar to the hematopoietic system.
Second, based on its assayable functions and potential, the stem cell at the top of the hierarchy is defined as a skeletal stem cell (SSC). As noted earlier (Sacchetti et al., 2007) (Bianco et al., 2013), this term clarifies, well beyond semantics, that the range of tissues that the self-renewing stromal progenitor (originally referred to as an ‘‘osteogenic’’ or ‘‘stromal’’ stem cell) (Owen and Friedenstein, 1988) can actually generate in vivo, overlaps with the range of tissues that make up the skeleton.
Third, these cells are spatially restricted, local residents of the bone/bone marrow organ. The systemic circulation is not a sizable contributor to their recruitment to locally deployed functions.
Fourth, a native skeletogenic potential is inherent to the system of progenitor/ stem cells found in the skeleton, and internally regulated by bone morphogenetic protein (BMP) signaling. This is reflected in the expression of regulators and antagonists of BMP signaling within the system, highlighting potential feedback mechanisms modulating expansion or quiescence of specific cell compartments.
Fifth, in cells isolated from other tissues, an assayable skeletogenic potential is not inherent: it can only be induced de novo by BMP reprogramming. These two studies (Chan et al., 2015, Worthley et al., 2015) corroborate the classical concept of ‘‘determined’’ and ‘‘inducible’’ skeletal progenitors (Owen and Friedenstein, 1988): the former residing in the skeleton, the latter found in nonskeletal tissues; the former capable of generating skeletal tissues, in vivo and spontaneously, the latter requiring reprogramming signals in order to acquire a skeletogenic capacity; the former operating in physiological bone formation, the latter in unwanted, ectopic bone formation in diseases such as fibrodysplasia ossificans progressiva.
To optimize our ability to obtain specific skeletal tissues for medical application, the study by Chan et al. offers a glimpse of another facet of the biology of SSC lineages and progenitors. Chan et al. show that a homogeneous cell population inherently committed to chondrogenesis can alter its output to generate bone if cotransplanted with multipotent progenitors. Conversely, osteogenic cells can be shifted to a chondrogenic fate by blockade of vascular endothelial growth factor receptor, consistent with the avascular and hypoxic milieu of cartilage. This has two important implications:
commitment is flexible in the system;
the choir is as important as the soloist and can modulate the solo tune.
Reversibility and population behavior thus emerge as two features that may be characteristic, albeit not unique, of the stromal system, resonating with conceptually comparable evidence in the human system.
The two studies by Chan et al. and Worthely et al. emphasize the relevance not only of their new data, but also of a proper concept of a skeletal stem cell per se, for proper clinical use. Confusion arising from improper conceptualization of skeletal stem cells has markedly limited clinical development of skeletal stem cell biology.
Gremlin 1 Identifies a Skeletal Stem Cell with Bone, Cartilage, and Reticular Stromal Potential
The stem cells that maintain and repair the postnatal skeleton remain undefined. One model suggests that perisinusoidal mesenchymal stem cells (MSCs) give rise to osteoblasts, chondrocytes, marrow stromal cells, and adipocytes, although the existence of these cells has not been proven through fate-mapping experiments. We demonstrate here that expression of the bone morphogenetic protein (BMP) antagonist gremlin 1 defines a population of osteochondroreticular (OCR) stem cells in the bone marrow. OCR stem cells self-renew and generate osteoblasts, chondrocytes, and reticular marrow stromal cells, but not adipocytes. OCR stem cells are concentrated within the metaphysis of long bones not in the perisinusoidal space and are needed for bone development, bone remodeling, and fracture repair. Grem1 expression also identifies intestinal reticular stem cells (iRSCs) that are cells of origin for the periepithelial intestinal mesenchymal sheath. Grem1 expression identifies distinct connective tissue stem cells in both the bone (OCR stem cells) and the intestine (iRSCs).
Identification and Specification of the Mouse Skeletal Stem Cell
How are skeletal tissues derived from skeletal stem cells? Here, we map bone, cartilage, and stromal development from a population of highly pure, postnatal skeletal stem cells (mouse skeletal stem cells, mSSCs) to their downstream progenitors of bone, cartilage, and stromal tissue. We then investigated the transcriptome of the stem/progenitor cells for unique gene-expression patterns that would indicate potential regulators of mSSC lineage commitment. We demonstrate that mSSC niche factors can be potent inducers of osteogenesis, and several specific combinations of recombinant mSSC niche factors can activate mSSC genetic programs in situ, even in nonskeletal tissues, resulting in de novo formation of cartilage or bone and bone marrow stroma. Inducing mSSC formation with soluble factors and subsequently regulating the mSSC niche to specify its differentiation toward bone, cartilage, or stromal cells could represent a paradigm shift in the therapeutic regeneration of skeletal tissues.
Bone mesenchymal development
Bone mesenchymal development
The bone-remodeling cycle
Nuclear receptor modulation – Role of coregulators in selective estrogen receptor modulator (SERM) actions
Selective estrogen receptor modulators (SERMs) are a class of small-molecule chemical compounds that bind to estrogen receptor (ER) ligand binding domain (LBD) with high affinity and selectively modulate ER transcriptional activity in a cell- and tissue-dependent manner. The prototype of SERMs is tamoxifen, which has agonist activity in bone, but has antagonist activity in breast. Tamoxifen can reduce the risk of breast cancer and, at same time, prevent osteoporosis in postmenopausal women. Tamoxifen is widely prescribed for treatment and prevention of breast cancer. Mechanistically the activity of SERMs is determined by the selective recruitment of coactivators and corepressors in different cell types and tissues. Therefore, understanding the coregulator function is the key to understanding the tissue selective activity of SERMs.
Hematopoietic
Hematopoietic Stem Cell Arrival Triggers Dynamic Remodeling of the Perivascular Niche
Hematopoietic stem and progenitor cells (HSPCs) can reconstitute and sustain the entire blood system. We generated a highly specific transgenic reporter of HSPCs in zebrafish. This allowed us to perform high resolution live imaging on endogenous HSPCs not currently possible in mammalian bone marrow. Using this system, we have uncovered distinct interactions between single HSPCs and their niche. When an HSPC arrives in the perivascular niche, a group of endothelial cells remodel to form a surrounding pocket. This structure appears conserved in mouse fetal liver. Correlative light and electron microscopy revealed that endothelial cells surround a single HSPC attached to a single mesenchymal stromal cell. Live imaging showed that mesenchymal stromal cells anchor HSPCs and orient their divisions. A chemical genetic screen found that the compound lycorine promotes HSPC-niche interactions during development and ultimately expands the stem cell pool into adulthood. Our studies provide evidence for dynamic niche interactions upon stem cell colonization.
Neonatal anemia and the need for red blood cell (RBC) transfusions are very common in neonatal intensive care units. Neonatal anemia can be due to blood loss, decreased RBC production, or increased destruction of erythrocytes. Physiologic anemia of the newborn and anemia of prematurity are the two most common causes of anemia in neonates. Phlebotomy losses result in much of the anemia seen in extremely low birthweight infants (ELBW). Accepting a lower threshold level for transfusion in ELBW infants can prevent these infants being exposed to multiple donors.
Red blood cell (RBC) transfusions are administered to neonates and premature infants using poorly defined indications that may result in unintentional adverse consequences. Blood products are often manipulated to limit potential adverse events, and meet the unique needs of neonates with specific diagnoses. Selection of RBCs for small volume (5–20 mL/kg) transfusions and for massive transfusion, defined as extracorporeal bypass and exchange transfusions, are of particular concern to neonatologists. Mechanisms and therapeutic treatments to avoid transfusion are another area of significant investigation. RBCs collected in anticoagulant additive solutions and administered in small aliquots to neonates over the shelf life of the product can decrease donor exposure and has supplanted the use of fresh RBCs where each transfusion resulted in a donor exposure. The safety of this practice has been documented and procedures established to aid transfusion services in ensuring that these products are available. Less well established are the indications for transfusion in this population; hemoglobin or hematocrit alone are insufficient indications unless clinical criteria (e.g. oxygen desaturation, apnea and bradycardia, poor weight gain) also augment the justification to transfuse. Comorbidities increase oxygen consumption demands in these infants and include bronchopulmonary dysplasia, rapid growth and cardiac dysfunction. Noninvasive methods or assays have been developed to measure tissue oxygenation; however, a true measure of peripheral oxygen offloading is needed to improve transfusion practice and determine the value of recombinant products that stimulate erythropoiesis. The development of such noninvasive methods is especially important since randomized, controlled clinical trials to support specific practices are often lacking, due at least in part, to the difficulty of performing such studies in tiny infants. The Effect of Blood Transfusion on the Hemoglobin Oxygen Dissociation Curve of Very Early Preterm Infants During the First Week of Life
This study was conducted during the first week of life to determine the changes in Ps0 (PO2 required to achieve a saturation of 50% at pH 7.4 and 37~ and the proportions of fetal hemoglobin (I-IbF) and adult hemoglobin (HbA) prior to and after transfusion in very early preterm infants. Eleven infants with a gestational age <–27 weeks have been included in study. The hemoglobin dissociation curve and the Ps0 was determined by Hemox-analyser. Liquid chromatography was also performed to determine the proportions of HbF and HbA. The mean gestational age of the 11 infants was 25.1 weeks (-+1 weeks) and their mean birth weight was 736 g (-+125 g). They received 26.9 mL/kg of packed red cells. The mean Ps0 prior and after transfusion was 18.5 +- 0.8 and 21.0 + 1 mm Hg (P = .0003) while the mean percentage of HbF was 92.9 -+ 1.1 and 42.6 -+ 5.7%, respectively. The data of this study show a decrease of hemoglobin oxygen affinity as a result of blood transfusion in very early preterm infants prone to O 2 toxicity. The shift in HbO 2 curve after transfusion should be taken into consideration when oxygen therapy is being regulated for these infants.
Effect of neonatal hemoglobin concentration on long-term outcome of infants affected by fetomaternal hemorrhage
Background: Fetomaternal hemorrhage (FMH) can cause severe morbidity. However, perinatal risk factors for long-term poor outcome due to FMH have not been extensively studied. Aims: To determine which FMH infants are likely to have neurological sequelae. Study design: A single-center retrospective observational study. Perinatal factors, including demographic characteristics, Kleihauer–Betke test, blood gas analysis, and neonatal blood hemoglobin concentration ([Hb]), were analyzed in association with long-term outcomes. Subjects: All 18 neonates referred to a Neonatal Intensive Care Unit of Kagoshima City Hospital and diagnosed with FMH during a 15-year study period. All had a neonatal [Hb] b7.5 g/dL and 15 of 17 neonates tested had Kleihauer–Betke test result N4.0%. Outcome measures: Poor long-term outcome was defined as any of the following determined at 12 month old or more: cerebral palsy, mental retardation, attention deficit/hyperactivity disorder, and epilepsy. Results: Nine of the 18 neonates exhibited poor outcomes. Among demographic characteristics and blood variables compared between two groups with poor and favorable outcomes, significant differences were observed in [Hb] (3.6 ± 1.4 vs. 5.4 ± 1.1 g/dL, P = 0.01), pH (7.09 ± 0.11 vs. 7.25 ± 0.13, P = 0.02) and base deficits (17.5 ± 5.4 vs. 10.4 ± 6.0 mmol/L, P = 0.02) in neonatal blood, and a number of infants with [Hb] ≤ 4.5 g/dL (78%[7/9] vs. 22%[2/9], P= 0.03), respectively. The base deficit in neonatal arterial blood increased significantly with decreasing neonatal [Hb]. Conclusions: Severe anemia causing severe base deficit is associated with neurological sequelae in FMH infants
Clinical and hematological presentation among Indian patients with common hemoglobin variants
Background: Co-inheritance of structural hemoglobin variants like HbS, HbD Punjab and HbE can lead to a variable clinical presentation and only few cases have been described so far in the Indian population. Methods: We present the varied clinical and hematological presentation of 22 cases (HbSD Punjab disease-15, HbSE disease-4, HbD Punjab E disease-3) referred to us for diagnosis. Results: Two of the 15 HbSDPunjab disease patients had moderate crisis, one presented with mild hemolytic anemia; however, the other 12 patients had a severe clinical presentation with frequent blood transfusion requirements, vaso occlusive crisis, avascular necrosis of the femur and febrile illness. The 4 HbSE disease patients had a mild to moderate presentation. Two of the 3 HbD Punjab E patients were asymptomatic with one patient’s sibling having a mild presentation. The hemoglobin levels of the HbSD Punjab disease patients ranged from 2.3 to 8.5 g/dl and MCV from 76.3 to 111.6 fl. The hemoglobin levels of the HbD Punjab E and HbSE patients ranged from 10.8 to 11.9 and 9.8 to 10.0 g/dl whereas MCV ranged from 67.1 to 78.2 and 74.5 to 76.0 fl respectively. Conclusions: HbSD Punjab disease patients should be identified during newborn screening programs and managed in a way similar to sickle cell disease. Couple at risk of having HbSD Punjab disease children may be given the option of prenatal diagnosis in subsequent pregnancies.
Sickle cell anemia is the most common hemoglobinopathy seen across the world. It is caused by a point mutation in the 6th codon of the beta (β) globin gene leading to the substitution of the amino acid glutamic acid to valine. The sickle gene is frequently seen in Africa, some Mediterranean countries, India, Middle East—Saudi Arabia and North America. In India the prevalence of hemoglobin S (HbS) carriers varies from 2 to 40% among different population groups and HbS is mainly seen among the scheduled tribe, scheduled caste and other backward class populations in the western, central and parts of eastern and southern India. Sickle cell anemia has a variable clinical presentation in India with the most severe clinical presentation seen in central India whereas patients in the western region show a mild to moderate clinical presentation.
Hemoglobin D Punjab (HbD Punjab) (also known as HbD Los-Angeles, HbD Portugal, HbD North Carolina, D Oak Ridge and D Chicago) is another hemoglobin variant due to a point mutation in codon 121 of the β globin gene resulting in the substitution of the amino acid glutamic acid to glycine. It is a widely distributed hemoglobin with a relatively low prevalence of 0.86% in the Indo-Pak subcontinent, 1–3% in north-western India, 1–3% in the Black population in the Caribbean and North America and has also been reported among the English. It accounts for 55.6% of all the Hb variants seen in the Xenjiang province of China.
Hemoglobin E (HbE) is the most common abnormal hemoglobin in Southeast Asia. In India, the frequency ranges from 4% to 51% in the north eastern region and 3% to 4% in West Bengal in the east. The HbE mutation (β26 GAG→AAG) creates an alternative splice site and the βE chain is insufficiently synthesized, hence the phenotype of this disorder is that of a mild form of β thalassemia.
Though these 3 structural variants are prevalent in different regions of India, their interaction is increasingly seen in all states of the country due to migration of people to different regions for a better livelihood. There are very few reports on interaction of these commonly seen Hb variants and the phenotypic–genotypic presentation of these cases is important for genetic counseling and management.
HbF of patients with HbSD Punjab disease with variable clinical severity. The HbF values of 4 patients are not included as they were post blood transfusion
The genotypes of the patients were confirmed by restriction enzyme digestion and ARMS (Fig). Patients 1 to 15 were characterized as compound heterozygous for HbS and HbD Punjab whereas patients 16 to 19 were characterized as compound heterozygous for HbS and HbE. Patient nos. 20 to 22 were characterized as compound heterozygous for HbE and HbD Punjab.
Molecular characterization of HbS and HbDPunjab by restriction enzyme digestion and of HbE by ARMS.
Molecular characterization of HbS and HbDPunjab by restriction enzyme digestion and of HbE by ARMS.
The 3 common β globin gene variants of hemoglobin, HbS, HbE and HbD Punjab are commonly seen in India, with HbS having a high prevalence in the central belt and some parts of western, eastern and southern India, HbE in the eastern and north eastern region whereas HbD is mostly seen in the north western part of India. These hemoglobin variants have been reported in different population groups. However, with migration and intermixing of the different populations from different geographic regions, occasional cases of HbSD Punjab and HbSE are being reported. There are several HbD variants like HbD Punjab, HbD Iran, HbD Ibadan. However, of these only HbD Punjab interacts with HbS to form a clinically significant condition as the glutamine residue facilitates polymerization of HbS. HbD Iran and HbD Ibadan are non-interacting and produce benign conditions like the sickle cell trait. The first case of HbSD Punjab disease was a brother and sister considered to have atypical sickle cell disease in 1934. This family was further reinvestigated and reported as the first case of HbD Los Angeles which has the same mutation as the HbD Punjab. Serjeant et al. reported HbD Punjab in an English parent in 6 out of 11 HbSD-Punjab disease cases. This has been suggested to be due to the stationing of nearly 50,000 British troops on the Indian continent for a period of 200 y and the introduction into Britain of their Anglo-Indian children.
HbSD Punjab disease shows a similar pattern to HbS homozygous on alkaline hemoglobin electrophoresis but can be differentiated on acid agar gel electrophoresis and on HPLC. In HbSD Punjab disease cases, the peripheral blood films show anisocytosis, poikilocytosis, target cells and irreversibly sickled cells. Values of HbF and HbA2 are similar to those in sickle homozygous cases. HbSD Punjab disease is characterized by a moderately severe hemolytic anemia.
Twenty-one cases of HbSDPunjab were reported by Serjeant of which 16 were reported by different workers among patients originating from Caucasian, Spanish, Australian, Irish, English, Portuguese, Black, American, Venezuelan, Caribbean, Mexican, Turkish and Jamaican backgrounds. Yavarian et al. 2009 reported a multi centric origin of HbD Punjab which in combination with HbS results in sickle cell disease. Patel et al. 2010 have also reported 12 cases of HbSD Punjab from the Orissa state of eastern India. Majority of these cases were symptomatic, presenting with chronic hemolytic anemia and frequent painful crises.
HbF levels >20% were seen in 4 out of our 11 clinically severe patients of HbSD-Punjab disease with the mean HbF levels of 16.8% in 8 clinically severe patients, while 3 clinically severe patients were post transfused. However, the 3 patients with a mild to moderate clinical presentation showed a mean HbF level of 8.6%. This is in contrast to the relatively milder clinical presentation associated with high HbF seen in patients with sickle cell anemia. This was also reported by Adekile et al. 2010 in 5 cases of HbS-DLos Angeles where high HbF did not ameliorate the severe clinical presentation seen in these patients.
These 15 cases of HbSDPunjab disease give us an overall idea of the severe clinical presentation of the disease in different regions of India. However the HbDPunjabE cases were milder or asymptomatic and the HbSE cases were moderately symptomatic. Since most of the cases of HbSDPunjab disease were clinically severe, it is important to pick up these cases during newborn screening and enroll them into a comprehensive care program with the other sickle cell disease patients with introduction of therapeutic interventions such as penicillin prophylaxis if required and pneumococcal immunization. In fact, 2 of our cases (No. 6 and 7) were identified during newborn screening for sickle cell disorders. The parents can be given information on home care and educated to detect symptoms that may lead to serious medical emergencies. The parents of these patients as well as the couples who are at risk of having a child with HbSDPunjab disease could also be counseled about the option of prenatal diagnosis in subsequent pregnancies. It is thus important to document the clinical and hematological presentation of compound heterozygotes with these common β globin chain variants.
Common Hematologic Problems in the Newborn Nursery
Common RBC disorders include hemolytic disease of the newborn, anemia, and polycythemia. Another clinically relevant hematologic issue in neonates to be covered herein is thrombocytopenia. Disorders of white blood cells will not be reviewed.
KEY POINTS
(1) Early clinical jaundice or rapidly developing hyperbilirubinemia are often signs of hemolysis, the differential diagnosis of which commonly includes immune-mediated disorders, red-cell enzyme deficiencies, and red-cell membrane defects.
(2) Knowledge of the maternal blood type and antibody screen is critical in identifying non-ABO alloantibodies in the maternal serum that may pose a risk for severe hemolytic disease in the newborn.
(3) Moderate to severe thrombocytopenia in an otherwise well-appearing newborn strongly suggests immune-mediated (alloimmune or autoimmune) thrombocytopenia.
Hemolytic conditions in the neonate
1. Immune-mediated (positive direct Coombs test) a. Rhesus blood group: Anti-D, -c, -C, -e, -E, CW, and several others
b. Non-Rhesus blood groups: Kell, Duffy, Kidd, Xg, Lewis, MNS, and others
c. ABO blood group: Anti-A, -B
2. Red blood cell (RBC) enzyme defects
a. Glucose-6-phosphate dehydrogenase (G6PD) deficiency
b. Pyruvate kinase deficiency
c. Others
3. RBC membrane defects
a. Hereditary spherocytosis
b. Elliptocytosis
c. Stomatocytosis
d. Pyknocytosis
e. Others
4. Hemoglobinopathies
a. alpha-thalassemia
b. gamma-thalassemia
Standard maternal antibody screeningAlloantibody Blood Group
D, C, c, E, e, f, CW, V Rhesus
K, k, Kpa, Jsa Kell
Fya, Fyb Duffy
Jka, Jkb Kidd
Xga Xg
Lea, Leb Lewis
S, s, M, N MNS
P1 P
Lub Lutheran
Non-ABO alloantibodies reported to cause moderate to severe hemolytic disease of the newbornWithin Rh system: Anti-D, -c, -C, -Cw, -Cx, -e, -E, -Ew, -ce, -Ces, -Rh29, -Rh32, -Rh42, -f, -G, -Goa, -Bea, -Evans, -Rh17, -Hro, -Hr, -Tar, -Sec, -JAL, -STEM
G6PD9 and pyruvate kinase (PK) deficiency are the 2 most common red-cell enzyme disorders associated with marked neonatal hyperbilirubinemia. Of these, G6PD deficiency is the more frequently encountered and it remains an important cause of kernicterus worldwide, including the United States, Canada, and the United Kingdom, the prevalence in Western countries a reflection in part of immigration patterns and intermarriage. The risk of kernicterus in G6PD deficiency also relates to the potential for unexpected rapidly developing extreme hyperbilirubinemia in this disorder associated with acute severe hemolysis.
Red Blood Cell Membrane Defects
Establishing a diagnosis of RBC membrane defects is classically based on the development of Coombs-negative hyperbilirubinemia, a positive family history, and abnormal RBC smear, albeit it is often difficult because newborns normally exhibit a marked variation in red-cell membrane size and shape. Spherocytes, however, are not often seen on RBC smears of hematologically normal newborns and this morphologic abnormality, when prominent, may yield a diagnosis of hereditary spherocytosis (HS) in the immediate neonatal period. Given that approximately 75% of families affected with hereditary spherocytosis manifest an autosomal dominant phenotype, a positive family history can often be elicited and provide further support for this diagnosis. More recently, Christensen and Henry highlighted the use of an elevated mean corpuscular hemoglobin concentration (MCHC) (>36.0 g/dL) and/or elevated ratio of MCHC to mean corpuscular volume, the latter they term the “neonatal HS index” (>0.36, likely >0.40) as screening tools for HS. An index of greater than 0.36 had 97% sensitivity, greater than 99% specificity, and greater than 99% negative predictive value for identifying HS in neonates. Christensen and colleagues also provided a concise update of morphologic RBC features that may be helpful in diagnosing this and other underlying hemolytic conditions in newborns.
The diagnosis of HS can be confirmed using the incubated osmotic fragility test when coupled with fetal red-cell controls or eosin-5-maleimide flow cytometry. One must rule out symptomatic ABO hemolytic disease by performing a direct Coombs test, as infants so affected also may manifest prominent micro-spherocytosis. Moreover, HS and symptomatic ABO hemolytic disease can occur in the same infant and result in severe hyperbilirubinemia and anemia. Of other red-cell membrane defects, only hereditary elliptocytosis, stomato-cytosis, and infantile pyknocytosis have been reported to exhibit significant hemolysis in the newborn period. Hereditary elliptocytosis and stomatocytosis are both rare. Infantile pyknocytosis, a transient red-cell membrane abnormality manifesting itself during the first few months of life, is more common.
Risk factors for bilirubinneurotoxicityIsoimmune hemolytic disease
G6PD deficiency
Asphyxia
Sepsis
Acidosis
Albumin less than 3.0 g/dL
Data from Maisels MJ, Bhutani VK, Bogen D, et al. Hyperbilirubinemia in the newborn infant > or 535 weeks’ gestation: an update with clarifications. Pediatrics 2009; 124:1193–8.
Polycythemia
Polycythemia (venous hematocrit 65%) in seen in infants across a range of conditions associated with active erythropoiesis or passive transfusion.76,77 They include, among others, placental insufficiency, the infant of a diabetic mother, recipient in twin-twin transfusion syndrome, and several aneuploidies, including trisomy. The clinical concern related to polycythemia is the risk for microcirculatory complications of hyperviscosity. However, determining which polycythemic infants are hyperviscous and when to intervene is a challenge.
Liver
Metabolic disorders presenting as liver disease
Germaine Pierre, Efstathia Chronopoulou
Paediatrics and Child Health 2013; 23(12): 509-514
The liver is a highly metabolically active organ and many inherited metabolic disorders have hepatic manifestations. The clinical presentation in these patients cannot usually be distinguished from liver disease due to acquired causes like infection, drugs or hematological disorders. Manifestations include acute and chronic liver failure, cholestasis and hepatomegaly. Metabolic causes of acute liver failure in childhood can be as high as 35%. Certain disorders like citrin deficiency and Niemann-Pick C disease may present in infancy with self-limiting cholestasis before presenting in later childhood or adulthood with irreversible disease. This article reviews important details from the history and clinical examination when evaluating the pediatric patient with suspected metabolic disease, the specialist and genetic tests when investigating, and also discusses specific disorders, their clinical course and treatment. The role of liver transplantation is also briefly discussed. Increased awareness of this group of disorders is important as in many cases, early diagnosis leads to early intervention with improved outcome. Diagnosis also allows genetic counselling and future family planning.
Adult liver disorders caused by inborn errors of metabolism: Review and update
Inborn errors of metabolism (IEMs) are a group of genetic diseases that have protean clinical manifestations and can involve several organ systems. The age of onset is highly variable but IEMs afflict mostly the pediatric population. However, in the past decades, the advancement in management and new therapeutic approaches have led to the improvement in IEM patient care. As a result, many patients with IEMs are surviving into adulthood and developing their own set of complications. In addition, some IEMs will present in adulthood. It is important for internists to have the knowledge and be familiar with these conditions because it is predicted that more and more adult patients with IEMs will need continuity of care in the near future. The review will focus on Wilson disease, alpha-1 antitrypsin deficiency, citrin deficiency, and HFE-associated hemochromatosis which are typically found in the adult population. Clinical manifestations and pathophysiology, particularly those that relate to hepatic disease as well as diagnosis and management will be discussed in detail.
Inborn errors of metabolism (IEMs) are a group of genetic diseases characterized by abnormal processing of biochemical reactions, resulting in accumulation of toxic substances that could interfere with normal organ functions, and failure to synthesize essential compounds. IEMs are individually rare, but collectively numerous. The clinical presentations cover a broad spectrum and can involve almost any organ system. The age of onset is highly variable but IEMs afflict mostly the pediatric population.
Wilson disease is an autosomal recessive genetic disorder of copper metabolism. It is characterized by an abnormal accumulation of inorganic copper in various tissues, most notably in the liver and the brain, especially in the basal ganglia. The disease was first described in 1912 by Kinnier Wilson, and affects between 1 in 30,000 and 1 in 100,000 individuals. Clinical features are variable and depend on the extent and the severity of copper deposition. Typically, patients tend to develop hepatic disease at a younger age than the neuropsychiatric manifestations. Individuals withWilson disease eventually succumb to complications of end stage liver disease or become debilitated from neurological problems, if they are left untreated.
The clinical presentations of Wilson disease are varied affecting many organ systems. However, the overwhelming majority of cases display hepatic and neurologic symptoms. In general, patients with hepatic disease present between the first and second decades of life although patients as young as 3 years old or over 50 years old have also been reported. The most common modes of presentations are acute self-limited hepatitis and chronic active hepatitis that are indistinguishable from other hepatic disorders although liver aminotransferases are generally much lower than in autoimmune or viral hepatitis. Acute fulminant hepatic failure is less common but is observed in approximately 3% of all cases of acute liver failure. Symptoms of acute liver failure include jaundice, coagulopathy, and hepatic encephalopathy. Cirrhosis can develop over time and may be clinically silent. Hepatocellular carcinoma (HCC) is rarely associated with Wilson disease, but may occur in the setting of cirrhosis and chronic inflammation.
Copper is an essential element, and is required for the proper functioning of various proteins and enzymes. The total body content of copper in a healthy adult individual is approximately 70–100 mg, while the daily requirements are estimated to be between 1 and 5 mg. Absorption occurs in the small intestine. Copper is taken up to the hepatocytes via the copper transporter hTR1. Once inside the cell, copper is bound to various proteins including metallothionein and glutathione, however, it is the metal chaperone, ATOX1 that helps direct copper to the ATP7B protein for intracellular transport and excretion. At the steady state, copper will be bound to ATP7B and is then incorporated to ceruloplasmin and secreted into the systemic circulation. When the cellular copper concentration arises, ATP7B protein will be redistributed from the trans-Golgi network to the prelysosomal vesicles facilitating copper excretion into the bile. The molecular defects in ATP7B lead to a reduction of copper excretion. Excess copper is accumulated in the liver causing tissue injury. The rate of accumulation of copper varies among individuals, and it may depend on other factors such as alcohol consumption, or viral hepatitis infections. If the liver damage is not severe, patients will accumulate copper in various tissues including the brain, the kidney, the eyes, and the musculoskeletal system leading to clinical disease. A failure of copper to incorporate into ceruloplasmin leads to secretion of the unsteady protein that has a shorter half-life, resulting in the reduced concentrations of ceruloplasmin seen in most patients with Wilson disease.
Wilson disease used to be a progressive fatal condition during the first half of the 20th century because there was no effective treatment available at that time. Penicillamine was the first pharmacologic agent introduced in 1956 for treating this condition. Penicillamine is a sulfhydryl-bearing amino acid cysteine doubly substituted with methyl groups. This drug acts as a chelating agent that promotes the urinary excretion of copper. It is rapidly absorbed in the gastrointestinal track, and over 80% of circulating penicillamine is excreted via the kidneys. Although it is very effective, approximately 10%–50% of Wilson disease patients with neuropsychiatric presentations may experience worsening of their symptoms, and often times the worsening symptoms may not be reversible.
Alpha1-antitrypsin deficiency
Alpha1-antitrypsin deficiency (AATD) is one of the most common genetic liver diseases in children and adults, affecting 1 in 2000 to 1 in 3000 live births worldwide. It is transmitted in an autosomal co-dominant fashion with variable expressivity. Alpha1 antitrypsin (A1AT) is a member of the serine protease inhibitor (SERPIN) family. Its function is to counteract the proteolytic effect of neutrophil elastase and other neutrophil proteases. Mutations in the SERPINA1, the gene encoding A1AT, result in changes in the protein structure with the PiZZ phenotype being the most common cause of liver and lung disease-associated AATDs. Although, it classically causes early onset chronic obstructive pulmonary disease (COPD) in adults, liver disease characterized by chronic inflammation, hepatic fibrosis, and cirrhosis is not uncommon in the adult population. Decreased plasma concentration of A1AT predisposes lung tissue to be more susceptible to injury from protease enzymes. However, the underlying mechanism of liver injury is different, and is believed to be caused by accumulation of polymerized mutant A1AT in the hepatocyte endoplasmic reticulum (ER). Currently, there is no specific treatment for liver disease-associated AATD, but A1AT augmentation therapy is available for patients affected with pulmonary involvement.
A1AT is a single-chain, 52-kDa polypeptide of approximately 394 amino acids [56]. It is synthesized in the liver, circulates in the plasma, and functions as an inhibitor of neutrophil elastase and other proteases such as cathepsin G, and proteinase 3. A1AT has a globular shape composed of two central β sheets surrounded by a small β sheet and nine α helices. The pathophysiology underlying liver disease is thought to be a toxic gain-of-function mutation associated with the PiZZ phenotypes. This hypothesis has been supported by the fact that null alleles which produce no detectable plasma A1AT, are not associated with liver disease. In addition, the transgenic mouse model of AATD PiZZ developed periodic acid-Schiff-positive diastase-resistant intrahepatic globule early in life similar to AATD patients. The PiZZ phenotype results in the blockade of the final processing of A1AT in the liver, as only 15% of the A1AT reaches the circulation whereas 85% of non-secreted protein is accumulated in the hepatocytes.
Citrin deficiency
Citrin deficiency is a relatively newly-defined autosomal recessive disease. It encompasses two different sub-groups of patients, neonatal intrahepatic cholestasis caused by citrin deficiency (NICCD), and adult onset citrullinemia type 2 (CTLN 2).
AGC2 exports aspartate out of the mitochondrial matrix in exchange for glutamate and a proton. Thus, this protein has an important role in ureagenesis and gluconeogenesis. In CTLN2, a defect in this protein is believed to limit the supply of aspartate for the formation of argininosuccinate in the cytosol resulting in impairment of ureagenesis. Interestingly, the mouse model of citrin deficiency (Ctrn−/−) fails to develop symptoms of CTLN2 suggesting that the mitochondrial aspartate is not the only source of ureagenesis. However, it should be noted that the rodent liver expresses higher glycerol-phosphate shuttle activity than the human counterpart. With the intact glycerol-phosphate dehydrogenase, it can compensate for the deficiency of AGC2, as demonstrated by the AGC2 and glycerol-phosphate dehydrogenase double knock-out mice that exhibit similar features to those observed in human CTLN2.
HFE-associated hemochromatosis
HFE-associated hemochromatosis is an inborn error of iron metabolism characterized by excessive iron storage resulting in tissue and organ damage. It is the most common autosomal recessive disorder in the Caucasian population, affecting 0.3%–0.5% of individuals of Northern European descent. The term “hemochromatosis” was coined in 1889 by the German pathologist Friedrich Daniel Von Recklinghausen, who described it as bronze stain of organs caused by a blood borne pigment.
The classic clinical triad of cirrhosis, diabetes, and bronze skin pigmentation is rarely observed nowadays given the early recognition, diagnosis, and treatment of this condition. The most common presenting symptoms are nonspecific including weakness, lethargy, and arthralgia.
The liver is a major site of iron storage in healthy individuals and as such it is the organ that is universally affected in HFE-associated hemochromatosis. Elevation of liver aminotransferases indicative of hepatocyte injury is the most common mode of presentation and it can be indistinguishable from other causes of hepatitis. Approximately 15%–40% of patients with HFE-associated hemochromatosis have other liver conditions, including chronic viral hepatitis B or C infection, nonalcoholic fatty liver disease, and alcoholic liver disease.
The review deals with genetic, regulatory and clinical aspects of iron homeostasis and hereditary hemochromatosis. Hemochromatosis was first described in the second half of the 19th century as a clinical entity characterized by excessive iron overload in the liver. Later, increased absorption of iron from the diet was identified as the pathophysiological hallmark. In the 1970s genetic evidence emerged supporting the apparent inheritable feature of the disease. And finally in 1996 a new “hemochromato-sis gene” called HFE was described which was mutated in about 85% of the patients. From the year2000 onward remarkable progress was made in revealing the complex molecular regulation of iron trafficking in the human body and its disturbance in hemochromatosis. The discovery of hepcidin and ferroportin and their interaction in regulating the release of iron from enterocytes and macrophages to plasma were important milestones. The discovery of new, rare variants of non-HFE-hemochromatosis was explained by mutations in the multicomponent signal transduction pathway controlling hepcidin transcription. Inhibited transcription induced by the altered function of mutated gene products, results in low plasma levels of hepcidin which facilitate entry of iron from enterocytes into plasma. In time this leads to progressive accumulation of iron and subsequently development of disease in the liver and other parenchymatous organs. Being the major site of excess iron storage and hepcidin synthesis the liver is a cornerstone in maintaining normal systemic iron homeostasis. Its central pathophysiological role in HFE-hemochromatosis with downgraded hepcidin synthesis, was recently shown by the finding that liver transplantation normalized the hepcidin levels in plasma and there was no sign of iron accumulation in the new liver.
Gastrointestinal
Decoding the enigma of necrotizing enterocolitis in premature infants
Necrotizing enterocolitis (NEC) is an enigmatic disease that affects primarily premature infants. It often occurs suddenly and when it occurs, treatment attempts at treatment often fail and results in death. If the infant survives, there is a significant risk of long term sequelae including neurodevelopmental delays. The pathophysiology of NEC is poorly understood and thus prevention has been difficult. In this review, we will provide an overview of why progress may be slow in our understanding of this disease, provide a brief review diagnosis, treatment and some of the current concepts about the pathophysiology of this disease.
Necrotizing enterocolitis (NEC) has been reported since special care units began to house preterm infants .With the advent of modern neonatal intensive care approximately 40 years ago, the occurrence and recognition of the disease markedly increased. It is currently the most common and deadly gastro-intestinal illness seen in preterm infants. Despite major efforts to better understand, treat and prevent this devastating disease, little if any progress has been made during these 4 decades. Underlying this lack of progress is the fact that what is termed “NEC” is likely more than one disease, or mimicked by other diseases, each with a different etiopathogenesis.
Human gut microbiome
Term or near term infants with “NEC” when compared to matched controls usually have occurrence of their disease in the first week after birth, have a significantly higher frequency of prolonged rupture of membranes, chorio-amnionitis, Apgar score <7 at 1 and 5 min, respiratory problems, congenital heart disease, hypoglycemia, and exchange transfusions. When a “NEC” like illness presents in term or near term infants, it should be noted that these are likely to be distinct in pathogenesis than the most common form of NEC and should be differentiated as such.
The infants who suffer primary ischemic necrosis are term or near term infants (although this can occur in preterms) who have concomitant congenital heart disease, often related to poor left ventricular output or obstruction. Other factors that have been associated with primary ischemia are maternal cocaine use, hyperviscosity caused by polycythemia or a severe antecedent hypoxic–ischemic event. Whether the dis-ease entity that results from this should be termed NEC can be debated on historical grounds, but the etiology is clearly different from the NEC seen in most preterm infants.
The pathogenesis of NEC is uncertain, and the etiology seems to be multifactorial. The “classic” form of NEC is highly associated with prematurity; intestinal barrier immaturity, immature immune response, and an immature regulation of intestinal blood flow (Fig.). Although genetics appears to play a role, the environment, especially a dysbiotic intestinal microbiota acting in concert with host immaturities predisposes the preterm infant to disruption of the intestinal epithelia, increased permeability of tight junctions, and release of inflammatory mediators that leads to intestinal mucosa injury and therefore development of necrotizing enterocolitis.
NEC is a multifactorial disease
What causes NEC? NEC is a multifactorial disease with an interaction of several etiophathologies
It is clear from this review that there are several entities that have been described as NEC. What is also clear is that despite having some overlap in the final parts of the pathophysiologic cascade that lead to necrosis, the disease that is most commonly seen in the preterm infant is likely to have an origin that differs markedly from that seen in term infants with congenital heart disease or severe hypoxic–ischemic injury. Thus, epidemiologic studies will need to differentiate these entities, if the aim is to dissect common features that are most highly associated with development of the disease. At this juncture, we areleft with more of a population based preventative approach, where the use of human milk, evidence based feeding guide-lines, considerations for microbial therapy once these are proved safe and effective and approved as such by regulatory authorities, and perhaps even measures that prevent prematurity will have a major impact on this devastating disease.
Influenced by the microbiota, intestinal epithelial cells (IECs) elaborate cytokines
Influenced by the microbiota, intestinal epithelial cells (IECs) elaborate cytokines, including thymic stromal lymphoprotein (TSLP), transforming growthfactor (TGF), and interleukin-10 (IL-10), that can influence pro-inflammatory cytokine production by dendritic cells (DC) and macrophages present in the laminapropria (GALT) and Peyer’s patches. Signals from commensal organisms may influence tissue-specific functions, resulting in T-cell expansion and regulation of the numbers of Th-1,
Th-2, and Th-3 cells. Also modulated by the microbiota, other IEC derived factors, including APRIL (a proliferation-inducing ligand),B-cell activating factor (BAFF), secretory leukocyte peptidase inhibitor (SLPI), prostaglandin E2(PGE2), and other metabolites, directly regulate functions ofboth antigen presenting cells and lymphocytes in the intestinal ecosystem. NK: natural killer cell; LN: lymph node; DC: dendritic cells.Modified from R. Sharma, C. Young, M. Mshvildadze, J. Neu, Intestinal microbiota does it play a role in diseases of the neonate? NeoReviews 10 (4) (2009)e166, with permission
Cross-talk between monocyte.macrophage cells and T.NK lymphocytes
Current Issues in the Management of Necrotizing Enterocolitis
Necrotizing enterocolitis is almost exclusively a disease of prematurity, with 90% of all cases occurring in premature infants and 90% of those infants weighing less than 2000 g. Prematurity is the only risk factor for necrotizing enterocolitis consistently identified in case control studies and the disease is rare in countries where prematurity is uncommon such as Japan and Sweden. When necrotizing enterocolitis does occur in full-term infants, it appears to by a somewhat different disease, typically associated with some predisposing condition.
NEC occurs in one to three in 1,000 live births and most commonly affects babies born between 30-32 weeks. It is most often diagnosed during the second week of life and occurs more often in previously fed infants. The mortality from NEC has been cited as 10% to 50% of all NEC cases. Surgical mortality has decreased over the last several decades from 70% to between 20 and 50%. The incremental cost per case of acute hospital care is estimated at $74 to 186 thousand compared to age matched controls, not including additional costs of long term care for the infants’ with lifelong morbidity. Survivors may develop short bowel syndrome, recurrent bouts of catheter-related sepsis, malabsorption, malnutrition, and TPN induced liver failure.
Although extensive research concerning the pathophysiology of necrotizing enterocolitis has occurred, a complete understanding has not been fully elucidated. The classic histologic finding is coagulation necrosis; present in over 90% of specimens. This finding suggests the importance of ischemia in the pathogenesis of NEC. Inflammation and bacterial overgrowth also are present. These findings support the assumptions by Kosloske that NEC occurs by the interaction of 3 events:
intestinal ischemia,
colonization by pathogenic bacteria and
excess protein substrate in the intestinal lumen.
Additionally, the immunologic immaturity of the neonatal gut has been implicated in the development of NEC. Reparative tissue changes including epithelial regeneration, formation of granulation tissue and fibrosis, and mixed areas of acute and chronic inflammatory changes suggest that the pathogenesis of NEC may involve a chronic process of injury and repair.
Premature newborns born prior to the 32nd week of gestational age may have compromised intestinal peristalsis and decreased motility. These motility problems may lead to poor clearance of bacteria, and subsequent bacterial overgrowth. Premature infants also have an immature intestinal tract in terms of immunologic immunity.
There are fewer functional B lymphocytes present and the ability to produce sufficient secretory IgA is reduced. Pepsin, gastric acid and mucus are also not produced as well in prematurity. All of these factors may contribute to the limited proliferation of intestinal flora and the decreased binding of these flora to mucosal cells (Fig).
Role of nitric oxide in the pathogenesis of NEC
Role of nitric oxide in the pathogenesis of NEC.
Characteristics of the immature gut leading to increased risk of necrotizing enterocolitis
Characteristics of the immature gut leading to increased risk of necrotizing enterocolitis.
As understanding of the pathophysiology of necrotizing enterocolitis continues to evolve, a unifying concept is emerging. Initially, there is likely a subclinical insult leading to NEC. This may arise from a brief episode of hypoxia or infection. With colonization of the intestines, bacteria bind to the injured mucosa eliciting an inflammatory response which leads to further inflammation.
Intestinal Microbiota Development in Preterm Neonates and Effect of Perinatal Antibiotics
Objectives Assess the establishment of the intestinal microbiota in very low birth-weight preterm infants and to evaluate the impact of perinatal factors, such as delivery mode and perinatal antibiotics. Study design We used 16S ribosomal RNA gene sequence-based microbiota analysis and quantitative polymerase chain reaction to evaluate the establishment of the intestinal microbiota. We also evaluated factors affecting the microbiota, during the first 3 months of life in preterm infants (n = 27) compared with full-term babies (n = 13). Results Immaturity affects the microbiota as indicated by a reduced percentage of the family Bacteroidaceae during the first months of life and by a higher initial percentage of Lactobacillaceae in preterm infants compared with full term infants. Perinatal antibiotics, including intrapartum antimicrobial prophylaxis, affects the gut microbiota, as indicated by increased Enterobacteriaceae family organisms in the infants.
Human gut microbiome
Conclusions Prematurity and perinatal antibiotic administration strongly affect the initial establishment of microbiota with potential consequences for later health.
Ischemia and necrotizing enterocolitis: where, when, and how
While it is accepted that ischemia contributes to the pathogenesis of necrotizing enterocolitis (NEC), three important questions regarding this role subsist. First, where within the intestinal circulation does the vascular pathophysiology occur? It is most likely that this event begins within the intramural microcirculation, particularly the small arteries that pierce the gut wall and the submucosal arteriolar plexus insofar as these represent the principal sites of resistance regulation in the gut. Mucosal damage might also disrupt the integrity or function of downstream villous arterioles leading to damage thereto; thereafter, noxious stimuli might ascend into the submucosal vessels via downstream venules and lymphatics. Second, when during the course of pathogenesis does ischemia occur? Ischemia is unlikely to the sole initiating factor of NEC; instead, it is more likely that ischemia is triggered by other events, such as inflammation at the mucosal surface. In this context, it is likely that ischemia plays a secondary, albeit critical role in disease extension. Third, how does the ischemia occur? Regulation of vascular resistance within newborn intestine is principally determined by a balance between the endothelial production of the vasoconstrictor peptide endothelin-1 (ET-1) and endothelial production of the vasodilator free radical nitric oxide (NO). Under normal conditions, the balance heavily favors NO-induced vasodilation, leading to a low resting resistance and high rate of flow. However, factors that disrupt endothelial cell function, eg, ischemia-reperfusion, sustained low-flow perfusion, or proinflammatory mediators, alter the ET-1:NO balance in favor of constriction. The unique ET-1–NO interaction thereafter might facilitate rapid extension of this constriction, generating a viscous cascade wherein ischemia rapidly extends into larger portions of the intestine.
Schematic representation of the intestinal microcirculation
Schematic representation of the intestinal microcirculation. Small mesenteric arteries pierce the muscularis layers and terminate in the submucosa where they give rise to 1A (1st order) arterioles. 2A (2nd order) arterioles arise from the 1A. Although not shown here, these 2A arterioles connect merge with several 1A arterioles, thus generating an arteriolar plexus, or manifold that serves to pressurize the terminal downstream microvasculature. 3A (3rd order) arterioles arise from the 2A and proceed to the mucosa, giving off a 4A branch just before descent into the mucosa. This 4A vessel travels to the muscularis layers. Each 3A vessel becomes the single arteriole perfusing each villus.
Collectively, these studies indicate that disruption of endothelial cell function has the potential to disrupt the normal balance between NO and ET-1 within the newborn intestinal circulation, and that such an event can generate significant ischemia. In this context, it is important to note that NO and ET-1 each regulate the expression and activity of the other. An increased [NO] within the microvascular environment reduces ET-1 expression and compromises ligand binding to the ETA receptor (thus decreasing its contractile efficacy), while ET-1 compromises eNOS expression. Thus, factors that upset the balance between NO and ET-1 will have an immediate and direct effect on vascular tone, but also exert an additional indirect effect by extenuating the disruption of balance between these two factors.
It is not difficult to construct a hypothesis that links the perturbations of I/R and sustained low-flow perfusion with an initial inflammatory insult. Initiation of an inflammatory process at the mucosal–luminal interface could have a direct impact on villus and mucosal 3A arterioles, damaging arteriolar integrity and disrupting villus hemodynamics. Ascent of proinflammatory mediators to the submucosal 1A–2A arteriolar plexus could occur via draining venules and lymphatics, generating damage to vascular effector systems therein; these mediators might include cytokines and platelet activating factor, as these elements have been recovered from human infants with NEC. This event, coupled with a generalized loss of 3A flow throughout a large portion of the mucosal surface, could compromise flow rate within the submucosal arteriolar plexus.
Necrotizing enterocolitis (NEC) is a leading cause of death among patients in the neonatal intensive care unit, carrying a mortality rate of 15e30%. Its pathogenesis is multifactorial and involves an over reactive response of the immune system to an insult. This leads to increased intestinal permeability, bacterial translocation, and sepsis. There are many inflammatory mediators involved in this process, but thus far none has been shown to be a suitable target for preventive or therapeutic measures. NEC usually occurs in the second week of life after the initiation of enteral feeds, and the diagnosis is made based on physical examination findings, laboratory studies, and abdominal radiographs. Neonates with NEC are followed with serial abdominal examinations and radiographs, and may require surgery or primary peritoneal drainage for perforation or necrosis. Many survivors are plagued with long term complications including short bowel syndrome, abnormal growth, and neurodevelopmental delay. Several evidence-based strategies exist that may decrease the incidence of NEC including promotion of human breast milk feeding, careful feeding advancement, and prophylactic probiotic administration in at-risk patients. Prevention is likely to have the greatest impact on decreasing mortality and morbidity related to NEC, as little progress has been made with regard to improving outcomes for neonates once the disease process is underway.
Immune Deficiencies
Primary immunodeficiencies: A rapidly evolving story
The characterization of primary immunodeficiencies (PIDs) in human subjects is crucial for a better understanding of the biology of the immune response. New achievements in this field have been possible in light of collaborative studies; attention paid to new phenotypes, infectious and otherwise; improved immunologic techniques; and use of exome sequencing technology. The International Union of Immunological Societies Expert Committee on PIDs recently reported on the updated classification of PIDs. However, new PIDs are being discovered at an ever-increasing rate. A series of 19 novel primary defects of immunity that have been discovered after release of the International Union of Immunological Societies report are discussed here. These new findings highlight the molecular pathways that are associated with clinical phenotypes and suggest potential therapies for affected patients.
Combined Immunodeficiencies
T-cell receptor a gene mutation: T-cell receptor ab1 T-cell depletion
T cells comprise 2 distinct lineages that express either ab or gd T-cell receptor (TCR) complexes that perform different tasks in immune responses. During T-cell maturation, the precise order and efficacy of TCR gene rearrangements determine the fate of the cells. Productive β-chain gene rearrangement produces a pre-TCR on the cell surface in association with pre-Tα invariant peptide (β-selection). Pre-TCR signals promote α-chain recombination and transition to a double-positive stage (CD41CD81). This is the prerequisite for central tolerance achieved through positive and negative selection of thymocytes.
Ras homolog gene family member H deficiency: Loss of naive T cells and persistent human papilloma virus infections
MST1 deficiency: Loss of naive T cells
New insight into the role of MST1 as a critical regulator of T-cell homing and function was provided by the characterization of 8 patients from 4 unrelated families who had homozygous nonsense mutations in STK4, the gene encoding MST1. MST1 was originally identified as an ubiquitously expressed kinase with structural homology to yeast Ste. MST1 is the mammalian homolog of the Drosophila Hippo protein, controlling cell growth, apoptosis, and tumorigenesis. It has both proapoptotic and antiapoptotic functions.
Lymphocyte-specific protein tyrosine kinase deficiency: T-cell deficiency with CD41 lymphopenia
Defects in pre-TCR– and TCR-mediated signaling lead to aberrant T-cell development and function (Fig). One of the earliest biochemical events occurring after engagement of the (pre)-TCR is the activation of lymphocyte-specific protein tyrosine kinase (LCK), a member of the SRC family of protein tyrosine kinases. This kinase then phosphorylates immunoreceptor tyrosine-based activation motifs of intracellular domains of CD3 subunits. Phosphorylated immunoreceptor tyrosine-based activation motifs recruit z-chain associated protein kinase of 70 kDa, which, after being phosphorylated by LCK, is responsible for activation of critical downstream events. Major consequences include activation of the membrane-associated enzyme phospholipase Cg1, activation of the mitogen-activated protein kinase, nuclear translocation of nuclear factor kB (NFkB), and Ca21/Mg21 mobilization. Through these pathways, LCK controls T-cell development and activation. In mice lacking LCK, T-cell development in the thymus is profoundly blocked at an early double-negative stage.
TCR signaling
TCR signaling. Multiple signal transduction pathways are stimulated through the TCR. These pathways collectively activate transcription factors that organize T-cell survival, proliferation, differentiation, homeostasis, and migration. Mutant molecules in patients with TCR-related defects are indicated in red.
Idiopathic CD41 lymphopenia (ICL) is a very heterogeneous clinical entity that is defined, by default, by persistent CD41 T-cell lymphopenia (<300 cells/mL or <20% of total T cells) in the absence of HIV infection or any other known cause of immunodeficiency.
Well-Defined Syndromes with Immunodeficiency
Wiskott-Aldrich syndrome protein–interacting protein deficiency: Wiskott-Aldrich syndrome-like phenotype
In hematopoietic cells Wiskott-Aldrich syndrome protein (WASP) is stabilized through forming a complex with WASP interacting protein (WIP).
Phospholipase Cg2 gain-of-function mutations: Cold urticaria, immunodeficiency, and autoimmunity/autoinflammatory
This is a unique phenotype, sharing features of antibody deficiency, autoinflammatory diseases, and immune dysregulatory disorders, making its classification difficult. Two recent studies validated the pleiotropy of genetic alterations in the same gene.
Predominantly Antibody Defects
Defect in the p85a subunit of phosphoinositide 3-kinase: Agammaglobulinemia and absent B cells
CD21 deficiency: Hypogammaglobulinemia
LPS-responsive beige-like anchor deficiency:
Hypogammaglobulinemia with autoimmunity and
early colitis
Defects Of Immune Dysregulation
Pallidin deficiency: Hermansky-Pudlak syndrome type 9
CD27 deficiency: Immune dysregulation and
persistent EBV infection
Congenital Defects Of Phagocyte Number, Function, Or Both
Interferon-stimulated gene 15 deficiency: Mendelian susceptibility to mycobacterial diseases
Minichromosome maintenance complex component 4 deficiency: NK cell deficiency associated with growth retardation and adrenal insufficiency
Autoinflammatory Disorders
A disintegrin and metalloproteinase 17 deficiency: Inflammatory skin and bowel disease
Cross-talk between monocyte.macrophage cells and T.NK lymphocytes
Cross-talk between monocyte/macrophage cells and T/NK lymphocytes. Genes in the IL-12/IFN-g pathway are particularly important for protection against mycobacterial disease. IRF8 is an IFN-g–inducible transcription factor required for the induction of various target genes, including IL-12. The NF-kB essential modulator (NEMO) mutations in the LZ domain impair CD40-NEMO–dependent pathways. Some gp91phox mutations specifically abolish the respiratory burst in monocyte-derived macrophages. ISG15 is secreted by neutrophils and potentiates IFN-g production by NK/T cells. Genetic defects that preclude monocyte development (eg, GATA2) can also predispose to mycobacterial infections (not shown). Mutant molecules in patients with unusual susceptibility to infection are indicated in red.
The field of PIDs is advancing at full speed in 2 directions. New genetic causes of known PIDs are being discovered (eg, CD21 and TRIF). Moreover, new phenotypes qualify as PIDs with the identification of a first genetic cause (eg, generalized pustular psoriasis). Recent findings contribute fundamental knowledge about immune system biology and its perturbation in disease. They are also of considerable clinical benefit for the patients and their families. A priority is to further translate these new discoveries into improved diagnostic methods and more effective therapeutic strategies, promoting the well-being of patients with PIDs.
In the last years, advances in molecular genetics and immunology have resulted in the identification of a growing number of genes causing primary immunodeficiencies (PIDs) in human subjects and a better understanding of the pathophysiology of these disorders. Characterization of the molecular mechanisms of PIDs has also facilitated the development of novel diagnostic assays based on analysis of the expression of the protein encoded by the PID-specific gene. Pilot newborn screening programs for the identification of infants with severe combined immunodeficiency have been initiated. Finally, significant advances have been made in the treatment of PIDs based on the use of subcutaneous immunoglobulins, hematopoietic cell transplantation from unrelated donors and cord blood, and gene therapy. In this review we will discuss the pathogenesis, diagnosis, and treatment of PIDs, with special attention to recent advances in the field.
Heartbeat Protein Deficits Linked to Sudden Cardiac Deaths
Reporter: Aviva Lev-Ari, PhD, RN
If only the molecular mechanisms underpinning arrhythmias were better understood, it would be possible to develop more targeted drugs. To study these mechanisms, and hopefully direct research to more effective anti-arrhythmia drugs, researchers at NYU Langone Medical Center decided to take a close look at a protein known as Pcp4, a known regulator of the heart’s rhythm. The researchers found that when Pcp4 is relatively scarce, as it is when expression of the Pcp4 gene is disrupted, ventricular arrhythmias may result.
The results of this investigation into Pcp4’s functions appeared October 8 in The Journal of Clinical Investigation, in an article entitled, “PCP4 regulates Purkinje cell excitability and cardiac rhythmicity.” In this article, the researchers described how they were able to isolate cardiac Purkinje cells in a mouse model of cardiomyopathy and show that Pcp4 expression was down-regulated in the diseased hearts.