Cardiac Genetic Predisposition in Sudden Infant Death Syndrome
Cardiac Genetic Predisposition in Sudden Infant Death Syndrome
Cardiac Genetic Predisposition in Sudden Infant Death Syndrome
11, 2018
ABSTRACT
BACKGROUND Sudden infant death syndrome (SIDS) is a leading cause of postneonatal mortality. Genetic heart dis-
eases (GHDs) underlie some cases of SIDS.
OBJECTIVES This study aimed to determine the spectrum and prevalence of GHD-associated mutations as a potential
monogenic basis for SIDS.
METHODS A cohort of 419 unrelated SIDS cases (257 male; average age 2.7 1.9 months) underwent whole exome
sequencing and a targeted analysis of 90 GHD-susceptibility genes. The yield of “potentially informative,” ultra-rare
variants (minor allele frequency <0.00005) in GHD-associated genes was assessed.
RESULTS Overall, 53 of 419 (12.6%) SIDS cases had $1 “potentially informative,” GHD-associated variant. The yield was
14.9% (21 of 141) for mixed-European ancestry cases and 11.5% (32 of 278) for European ancestry SIDS cases.
Infants older than 4 months were more likely to host a “potentially informative” GHD-associated variant. There was
significant overrepresentation of ultra-rare nonsynonymous variants in European SIDS cases (18 of 278 [6.5%]) versus
European control subjects (30 of 973 [3.1%]; p ¼ 0.013) when combining all 4 major cardiac channelopathy genes
(KCNQ1, KCNH2, SCN5A, and RYR2). According to the American College of Medical Genetics guidelines, only 18 of
419 (4.3%) SIDS cases hosted a “pathogenic” or “likely pathogenic” variant.
CONCLUSIONS Less than 15% of more than 400 SIDS cases had a “potentially informative” variant in a
GHD-susceptibility gene, predominantly in the 4- to 12-month age group. Only 4.3% of cases possessed immediately
clinically actionable variants. Consistent with previous studies, ultra-rare, nonsynonymous variants within the major
cardiac channelopathy-associated genes were overrepresented in SIDS cases in infants of European ethnicity. These
findings have major implications for the investigation of SIDS cases and families. (J Am Coll Cardiol 2018;71:1217–27)
© 2018 the American College of Cardiology Foundation. Published by Elsevier. All rights reserved.
From the aDepartments of Cardiovascular Medicine (Division of Heart Rhythm Services), Pediatrics (Division of Pediatric Cardi-
ology), and Molecular Pharmacology & Experimental Therapeutics (Windland Smith Rice Sudden Death Genomics Laboratory),
Mayo Clinic, Rochester, Minnesota; bMolecular and Clinical Sciences Research Institute, St. George’s, University of London,
London, United Kingdom; cCardiology Clinical Academic Group, St. George’s University Hospitals’ NHS Foundation Trust, Lon-
don, United Kingdom; dMedical and Molecular Genetics, Guy’s Hospital, King’s College London, London, United Kingdom; eMRC
Human Genetics Unit, University of Edinburgh, Edinburgh, United Kingdom; fRoyal Infirmary of Edinburgh, Edinburgh, United
Kingdom; gCentre for Child and Adolescent Health, Bristol Medical School, University of Bristol, Bristol, United Kingdom;
h
Department of Cellular Pathology, St George’s, University of London, London, United Kingdom; iDepartment of Cellular
Pathology’, St. George’s University Hospitals’ NHS Foundation Trust, London, United Kingdom; jHistopathology Department,
Sheffield Children’s Hospital, Sheffield, United Kingdom; kHonorary Senior Lecturer, University of Sheffield, Sheffield, United
Kingdom; lDepartment of Cardiology, The Heart Centre, Rigshospitalet, Copenhagen University Hospital, Copenhagen, Denmark;
and the mDepartment of Forensic Medicine, Faculty of Medical Sciences, University of Copenhagen, Copenhagen, Denmark. This
work was supported by the Eunice Kennedy Shriver National Institute of Child Health & Human Development of the National In-
Listen to this manuscript’s
stitutes of Health (grant no. R01HD042569 to Dr. Ackerman) and by the British Heart Foundation (BHF Clinical Research Training
audio summary by
Fellowship FS/13/78/30520 to Drs. Wong and Behr). The content is solely the responsibility of the authors and does not necessarily
JACC Editor-in-Chief
represent the official views of the National Institutes of Health. Mr. Tester is supported by the Mayo Clinic Windland Smith Rice
Dr. Valentin Fuster.
Comprehensive Sudden Cardiac Death Program. Dr. Wong is supported by additional funds from Biotronik and Cardiac Risk in the
Young. Dr. Behr is supported by the Higher Education Funding Council for England; is a consultant for Medtronic; has received
research funding from Biotronik; and has received funds from The Robert Lancaster Memorial Fund sponsored by McColl’s RG Ltd.
Dr. Ackerman is a consultant for Audentes Therapeutics, Boston Scientific, Gilead Sciences, Invitae, Medtronic, MyoKardia, and
St. Jude Medical; is supported by the Mayo Clinic Windland Smith Rice Comprehensive Sudden Cardiac Death Program; and the Mayo
Clinic have an equity- or royalty-based relationship with AliveCor, Blue Ox Health, and StemoniX, although none of these entities
S METHODS
ABBREVIATIONS udden infant death syndrome (SIDS) is
AND ACRONYMS the sudden unexpected death of an in-
fant <1 year of age that remains unex- STUDY GROUP. The SIDS cohort (N ¼ 427) consisted
GHD = genetic heart disease
plained despite comprehensive clinical and of 95 coroners’ cases from the United Kingdom
gnomAD = Genome
pathological investigations (1). SIDS repre- (London, Sheffield, Edinburgh, and Bristol) and 332
Aggregation Database
sents 70% to 80% of all sudden unexpected coroner-, medical examiner–, or forensic pathologist–
LQTS = long QT syndrome
infant deaths with an incidence of 0.4 in referred cases collected from 6 ethnically and
MAF = minor allele frequency
1,000 live births in the United Kingdom and geographically diverse U.S. population groups.
NSV = nonsynonymous variant
0.5 in 1,000 live births in the United States Because of the lack of uniformity in procedures and
PCA = principal component
(2,3). The peak incidence occurs between 2 reporting among medical examiner offices in the
analysis
and 4 months of age and is more common United States, minor differences in protocol may
SIDS = sudden infant death
syndrome
in boys. Such infant deaths are commonly exist. Nonetheless, both gross and histological ex-
associated with environmental risk factors aminations of all major organs were performed, and
WES = whole exome
sequencing such as co-sleeping or prone sleeping posi- all cases satisfied our enrollment criteria, which
tion (4). Despite successful targeted risk included the following: 1) sudden unexplained death
reduction campaigns, the number of SIDS cases has of an infant <1 year of age; 2) European descent; and
plateaued, and SIDS remains the leading cause of 3) a comprehensive negative medicolegal autopsy
postneonatal death (4). including a negative toxicology screen result and
death scene investigation. Infants with asphyxia or a
SEE PAGE 1228
specific disease causing death were excluded.
A triple-risk model for SIDS suggests the conver- Ethnicity was self-reported by the referring coroner
gence of the vulnerable infant in the setting of or medical examiner. This anonymous autopsy study
exogenous stressors during a critical development had only limited medical information available such
period (5) (Central Illustration). Although many path- as the sex, ethnicity, age at the time of death, and
ophysiological theories have been proposed, decisive sleep position. This study complies with the Decla-
pathogenic substrates or mechanisms triggering an ration of Helsinki; locally appointed ethics commit-
infant’s sudden demise remain unclear (6–9). Several tees including Mayo Clinic’s Institutional Review
studies have implicated both common and rare ge- Board approved the research protocol. Some of the
netic variants involved in autonomic function, 332 samples from the United States were included in
infection, and cardiac repolarization (10–14). In driven, specific candidate gene mutational analysis
addition, potentially lethal genetic heart diseases (10,18,19,23–28). Of the 332 U.S. cases, 58 had been
(GHDs) including long QT syndrome (LQTS), Brugada analyzed previously for variants in SCN5A (10),
syndrome, catecholaminergic polymorphic ventricu- KCNQ1 (18), KCNH2 (18), RYR2 (19), SNTA1 (23),
lar tachycardia, and hypertrophic cardiomyopathy KCNJ8 (24), Cx43 (25), GPD1L (26), CAV3 (27), SCN1B
have been implicated as monogenic causes for a small (28), SCN2B (28), SCN3B (28), and SCN4B (28). An
proportion of SIDS cases (10,13,15–27). additional 25 of the 332 cases were also analyzed for
However, fewer than 100 investigations of genetic RYR2 (19), and an additional 145 of the 332 cases were
variations in population-based SIDS cohorts have also analyzed for SNTA1 (23), KCNJ8 (24), Cx43 (25),
been published to date, largely on the basis of GPD1L (26), CAV3 (27), SCN1B (28), SCN2B (28), SCN3B
hypothesis-driven, candidate gene– or pathway- (28), and SCN4B (28). None of the 95 cases from the
based approaches that recognize established patho- United Kingdom have been published previously.
biological risk factors for SIDS, with an average cohort CONTROL GROUP. A total of 973 control exomes (509
size of just 125 SIDS cases (13). In the present study, female, 464 male) from the ICR1000 U.K. exome se-
using whole exome sequencing (WES), we conducted ries and the 1958 Birth Cohort study were included for
a GHD-associated gene-specific analysis on a cohort case-control analysis (29). As previously reported,
of more than 400 unrelated SIDS cases. exome sequencing was performed using the Illumina
were involved in this study in any way. Dr. Simpson has a part-time contract of service with Genomics. All other authors have
reported that they have no relationships relevant to the contents of this paper to disclose. Mr. Tester and Dr. Wong contributed
equally to this work and are joint first authors. Drs. Simpson, Behr, and Ackerman contributed equally to this work and are joint
senior authors.
Manuscript received November 8, 2017; revised manuscript received December 15, 2017, accepted January 8, 2018.
JACC VOL. 71, NO. 11, 2018 Tester et al. 1219
MARCH 20, 2018:1217–27 Cardiac Genetic Analyses in SIDS
C ENTR AL I LL U STRA T I O N Whole Exome Sequencing and a Targeted Analysis of 90 GHD-Susceptibility Genes
The triple-risk hypothesis for SIDS highlighting genetic heart disease (GHD) as a potential explanation for infant vulnerability, and our whole exome sequencing
strategy to detect American College of Medical Genetics and Genomics (ACMG) guideline–predicated “pathogenic” or “likely pathogenic” variants in SIDS cases.
MAF ¼ minor allele frequency; SIDS ¼ sudden infant death syndrome.
TruSeq and Illumina instruments (Illumina, San index DNA adaptors (Agilent) with a single “T” base
Diego, California) (29). overhang at the 3 0 end were ligated, and the resulting
constructs were purified using AMPure SPRI (solid-
WHOLE EXOME SEQUENCING. Genomic DNA iso- phase reversible immobilization) beads (Agencourt
lated from each SIDS case underwent WES at the KCL- Bioscience Corporation, Beverly, Massachusetts). The
GSTT Biomedical Research Centre Genomics Platform adapter-modified DNA fragments were enriched by 4
in London, United Kingdom or Mayo Clinic’s Medical cycles of polymerase chain reaction using SureSelect
Genome Facility in Rochester, Minnesota. forward and SureSelect ILM Pre-Capture Indexing
Paired-end libraries were prepared following the reverse (Agilent) primers. The concentration and size
manufacturer’s protocol (Agilent Technologies, Inc., distribution of the libraries was determined on an
Santa Clara, California) using the Bravo Automated Agilent Bioanalyzer DNA 1000 chip.
Liquid Handling Platform from Agilent. Briefly, 1 to Whole exon capture was performed using the
3 m g of genomic DNA was fragmented to 150 to 200 bp protocol for Agilent’s Sure SelectXT Human All Exon
using the Covaris E210 sonicator (Covaris, Inc., V5þUTR kit. Briefly, 750 ng of the prepared library
Woburn, Massachusetts). The ends were repaired, was incubated with whole exon biotinylated RNA
and an “A” base was added to the 3 0 ends. Paired-end capture baits supplied in the kit for 24 h at 65 C.
1220 Tester et al. JACC VOL. 71, NO. 11, 2018
The captured DNA:RNA hybrids were recovered using common variants located outside of regions of the
Dynabeads MyOne Streptavidin T1 (Invitrogen Dynal, genome where there is extensive linkage disequilib-
Thermo Fisher Scientific, Waltham, Massachusetts). rium was used to estimate relatedness within the
The DNA was eluted from the beads and purified us- study cohort and ethnic ancestry alongside the con-
ing Ampure XP beads (Agencourt). The purified cap- trol group (32). Estimation was undertaken using the
ture products were then amplified using the first 2 dimensions of a multidimensional scaling using
SureSelect Post-Capture Indexing forward and index Euclidean distance undertaken with the King soft-
polymerase chain reaction reverse primers (Agilent) ware package.
for 12 cycles. ANCESTRY CONFIRMATION. To avoid potential con-
Libraries were pooled at equimolar concentrations founding secondary to population stratification
and loaded onto paired-end flow cells at concentra- resulting from genetic admixture, a principal compo-
tions of 7 to 8 pM to generate cluster densities of nent analysis (PCA) was performed. The PCA served
600,000 to 800,000/mm 2 following Illumina’s stan- only for the rare variant analysis between European
dard protocol using the Illumina cBot and HiSeq SIDS cases and European control subjects. The PCA
paired-end cluster kit version 3. Each lane of a HiSeq data were not used for attributing causality to iden-
flow cell produced 21 to 39 Gb of sequence. The level tified variants where ethnically matched control
of sample pooling was controlled by the size of the subjects would not be necessary for variant adjudi-
capture region and the desired depth of coverage. cation. SIDS cases and control subjects forming a ho-
The flow cells were sequenced as 101 2 paired- mogeneous cluster on the first 2 components were
end reads on an Illumina HiSeq 2000 using TruSeq included in the case-control rare variant analysis.
SBS sequencing kit version 3 and HiSeq data collec-
tion version 2.0.12.0 software. Base-calling was per- GENETIC HEART DISEASE GENE-SPECIFIC VARIANT
formed using Illumina’s RTA version 1.17.21.3. ANALYSIS. Genes known to be associated with car-
The FASTQ files underwent quality control checks diac channelopathy (LQTS, catecholaminergic poly-
using FASTQC (Babraham Bioinformatics, Babraham morphic ventricular tachycardia, Brugada syndrome)
Institute, Cambridge, United Kingdom). The Illumina susceptibility and cardiomyopathy (hypertrophic car-
paired-end reads were aligned to the GRCh37 (hg19) diomyopathy, dilated cardiomyopathy, arrhythmo-
human reference genome using Novoalign (Novocraft genic cardiomyopathy) susceptibility (N ¼ 90) (Online
Technologies, Selandor, Malaysia). Single-sample Table 1) were evaluated for the presence of “ultra-rare”
variant calling with the Genome Analysis Toolkit nonsynonymous variants (NSVs) with a minor allele
(GATK version 3.2-2, Broad Institute, Cambridge, frequency (MAF) <0.00005 (1 in 20,000 alleles)
Massachusetts) (30), and the resulting gVCFs (genomic derived from the Genome Aggregation Database (33). A
variant call format) subsequently underwent multi- comparison of yield was undertaken for ultra-rare
sample genotyping and variant quality score recali- NSVs in SIDS cases of PCA-determined European
bration. Genotypes were excluded if the quality ancestry versus European control subjects across all 90
control was <15 or there were fewer than 4 reads GHD-susceptibility genes and the 4 “major” channel-
supporting the call. Further filtering of variant sites opathy genes (KCNQ1, KCNH2, SCN5A, and RYR2).
was performed to exclude sites with missingness >0.1 All putative loss of function variants (i.e., a “radical”
in cases or control subjects. Variants were annotated variant: frameshift, nonsense, and essential splice-site
with respect to the genes in which they reside with variants) or missense variants with a previously
Annovar, allele frequencies were obtained from established abnormal in vitro function characteriza-
the Exome Aggregation Consortium database, and tion that resided within any of the 90 GHD-associated
Combined Annotation Dependent Depletion (CADD) genes and all ultra-rare, missense variants residing in
scores were derived from the CADD server (31). any of the 4 “major” channelopathy genes were
considered to be “potentially informative” variants
QUALITY CONTROL COVERAGE ANALYSIS AND that would be appropriate for investigation of their
PRINCIPAL COMPONENT ANALYSIS FOR RELATEDNESS significance in a family. Such variants were confirmed
AND ETHNICITY. Coverage across the exome was using standard Sanger sequencing techniques. The
assessed using the bedtools software package (Quin- American College of Medical Genetics and Genomics
lan Laboratory, University of Utah, Salt Lake City, standards and guidelines for the interpretation of
Utah), and cases were excluded from further analysis sequence variants were used to assist in the classifi-
if <75% of the GENCODE (GENCODE Project, an in- cation of our genetic findings further among all ultra-
ternational collaboration)–defined protein coding rare (MAF <0.00005) NSVs identified across the 90
exome was covered by <20 reads. A set of 3,847 GHD-associated genes (34).
JACC VOL. 71, NO. 11, 2018 Tester et al. 1221
MARCH 20, 2018:1217–27 Cardiac Genetic Analyses in SIDS
T A B L E 1 Summary of the Sudden Infant Death Syndrome F I G U R E 1 Yield of Ultra-Rare and “Potentially Informative” GHD-Associated
Cohort Demographics Gene Variants
compared with the Fisher exact or chi-square tests. Overall, a total of 285 unique, ultra-rare NSVs (256
Probability values were determined on the basis of missense, 23 putative loss of function [12 frameshift,
2-sided tests considered significant at p < 0.05. 8 splice errors, 3 nonsense], and 6 in-frame indels)
Analysis was conducted with SPSS version 18.0 were identified in 194 of 419 (46.3%) SIDS cases
software (SPSS, Chicago, Illinois). overall (Figure 1). Further, 45 of 278 (16.2%) European
ancestry cases and 25 of 141 (17.7%) mixed-European
RESULTS ancestry cases hosted >1, ultra-rare NSVs.
These ultra-rare NSVs resided in 68 of the 90 GHD-
DEMOGRAPHICS. WES was performed in 427 SIDS associated genes (21 of 31 channelopathy-associated
cases. However, quality control metrics excluded 7 genes and 47 of 59 cardiomyopathy-associated
cases because of insufficient exome coverage and 1 genes). The gene-specific yields for the overall
individual from a half-sibling pair (Online Figure 1). cohort and the European and mixed European sub-
The cohort therefore consisted of 419 cases (257 male, sets are shown in Online Table 2.
162 female; average age 2.7 1.9 months) with a Of the 285 unique, ultra-rare NSVs identified, 57
skewed bell-shaped distribution of age (Online (20%) were considered “potentially informative.” A
Figure 2). The epidemiologically higher-risk age total of 25 of 57 (43.9%) were missense variants in the
group of 2 to 4 months (58.9%) and male sex (61%) 4 major channelopathy genes, 23 of 57 (40.3%) were
accounted for the majority of the cases. The PCA putative loss of function variants, and 10 of 57 (17.5%)
demonstrated a wide distribution of ancestral origins were variants previously reported in published pa-
with 278 cases (173 male, 105 female) considered to be pers as having an abnormal in vitro functional
of European ancestry and 141 cases (84 male, 57 phenotype (Table 2). Overall, 53 of 419 (12.6%) SIDS
female) considered to be of mixed-European ancestry cases hosted at least 1 “potentially informative”
(Online Figure 3, Table 1). Sleep characteristics were variant (Figure 1). Four of 419 cases (0.95%) had 2
known in 54% of the cohort (Table 1). There were no “potentially informative” variants. The yield was
significant differences in demographics and sleep 14.9% (21 of 141) for mixed-European ancestry cases
characteristics between the European and mixed- and 11.5% (32 of 278) for European ancestry SIDS cases
European ancestry cases. (Figure 1).
1222 Tester et al. JACC VOL. 71, NO. 11, 2018
T A B L E 2 Continued
ACMG ¼ American College of Medical Genetics and Genomics; MAF ¼ minor allele frequency; NSV ¼ nonsynonymous variant; SIDS ¼ sudden infant death syndrome; VUS ¼ variant of uncertain significance.
There were no significant differences in the yield of European (12 of 278 [4.3%]) and the mixed-European
either ultra-rare NSVs among all 90 GHD genes or cohorts (6 of 141 [4.3%]).
“potentially informative” variants when comparing CASE-CONTROL ANALYSIS. Consistent with previ-
sex, sleep position (supine vs. prone), or co-sleeping ous studies, there was significant overrepresentation
(yes vs. no) in either the overall or stratified study of ultra-rare NSVs in European SIDS cases (18 of 278
groups (Table 3). However, there was a significantly [6.5%]) versus European control subjects (30 of 973
higher yield of “potentially informative” GHD- [3.1%]; p ¼ 0.013) when combining all 4 major cardiac
associated genetic variants in those infants who channelopathy genes (KCNQ1, KCNH2, SCN5A, and
died at >4 months of age (15 of 65 [23.1%]) compared RYR2). (Figure 3). However, there was no significant
with those younger than 4 months of age (37 of 354 difference in yield between cases and control subjects
[10.4%]; p ¼ 0.0075) (Figure 2). for any specific gene.
Following further vetting using the strict American
College of Medical Genetics and Genomics guidelines, DISCUSSION
only 18 of the 285 ultra-rare NSVs achieved a “path-
ogenic” or “likely pathogenic” designation and were This paper reports results derived from a whole
identified in 18 (4.3%) of the 419 SIDS cases (Table 2, exome molecular autopsy with GHD gene-specific
Figure 1). There was no difference in yield of “path- analysis of a large cohort of unrelated SIDS cases.
ogenic” or “likely pathogenic” variants between the Previous postmortem genetic studies implicated
T A B L E 3 The Effect of Various Demographics on the Yield of Genetic Heart Disease–Associated Gene Variants
5 nelopathy genes.
Recently, Hertz et al. (22) reported a 34% yield of
4
3.1 3.2 “variants with likely functional effects” following a
3 genetic analysis of GHD-associated genes in only 47
2 1.7 sudden unexpected deaths in infancy cases. Howev-
1.1 1.4
er, given the rarity of GHDs in the general population,
1 0.7 0.5 0.7
0.2 we believe that their definition of rarity (MAF <1%)
0 was unacceptably and erroneously high, thus causing
4 Major KCNQ1 KCNH2 SCN5A RYR2
Genes an overestimated burden of potentially pathogenic
European SIDS (n = 278) European Controls (n = 973) variants in their SIDS cohort. In fact, of their 16
“pathogenic” variants, only 1 novel RYR2 variant
would have been deemed “potentially informative”
A bar graph depicting the percentage of yield of ultra-rare (minor allele
frequency <0.00005) variants in major cardiac channelopathy genes (KCNQ1, KCNH2, by our robust criteria.
SCN5A, and RYR2). SIDS ¼ sudden infant death syndrome. In 2017, Neubauer et al. (38) reported a yield
of “potentially causative” variants in 20% of their 155
JACC VOL. 71, NO. 11, 2018 Tester et al. 1225
MARCH 20, 2018:1217–27 Cardiac Genetic Analyses in SIDS
European SIDS cases following WES and genetic bases that are largely different genetically and
interrogation of their 192-gene focused panel that mechanistically from sudden death occurring after
comprised both cardiovascular-associated and meta- the age of 1 year.
bolic disorder–associated genes. The majority of their Several risk factors for SIDS have been established.
seemingly genotype-positive infants had a variant One could hypothesize that vulnerable infants dying
with “likely functional effects” in genes associated of SIDS without the presence of additional risk factors
with a cardiac channelopathy (9%) or cardiomyopathy are more likely to host a highly penetrant monogenic
(7%). However, most of these variants represent cause for their death compared with infants exposed
missense variants within “minor” genes (38). In fact, to additional environmental risk factors. Yet no sig-
only 2.6% of their cases hosted what we would nificant differences in the yield of “potentially
consider a “potentially informative” variant by our informative” GHD gene variants associated with sex,
strict definition. sleep position, or bed sharing were observed. How-
Although our study supports the utility of WES to ever, a significant age effect on the yield, where 23%
identify potential sudden death-causing variants of those infants older than 4 months of age hosted a
within established or potential sudden death– “potentially informative” GHD-associated variant
susceptibility genes, the challenge of the WES-based compared with only 10% of the infants younger than
molecular autopsy does not lie in the identification 4 months of age, was observed. These data support
of variants, but rather in the adjudication of their the potential stratification of those SIDS cases that
potential pathogenicity. Accurate variant classifica- may benefit most from postmortem genetic testing of
tion is crucial to enable proper counseling of surviv- the major channelopathy- or cardiomyopathy-
ing family members. Erroneously or prematurely associated genes.
adjudicating ambiguous variants as pathogenic has The significant overrepresentation of ultra-rare
the potential to harm patients and their families. NSVs within the 4 major channelopathy genes asso-
Tragically, this became a reality for 1 family described ciated with either inheritable LQTS (KCNQ1, KCNH2,
by Ackerman et al. (36) recently, as they dealt with SCN5A) or catecholaminergic polymorphic ventricular
the disastrous consequences of unnecessary treat- tachycardia (RYR2) observed in our European Cauca-
ment on the basis of an erroneously interpreted sian SIDS cases compared with ethnically matched
variant in KCNQ1. Thus, overattribution of SIDS control subjects (6.5% vs. 3.1%; p ¼ 0.013) supports
deaths to GHDs has significant implications for the that cardiac channelopathies may represent the un-
immediate family, and we urge extreme caution in derlying pathogenic basis for some SIDS cases and
variant interpretation. When such cautionary advice that post-mortem genetic testing of the 4 major
was heeded, <5% of more than 400 SIDS cases had channelopathy-associated genes may be warranted in
either a “pathogenic” or “likely pathogenic” variant cases of SIDS.
in 1 of 90 GHD-susceptibility genes, a percentage that
is substantially lower than in previous extrapolations STUDY LIMITATIONS. In our study, we used a strict
of the prevalence of either channelopathic or MAF cutoff of 0.005% (i.e., 1 in 20,000 alleles or 1
cardiomyopathic SIDS. This parallels our experience in 10,000 individuals). Although using a stringent
of the “molecular autopsy” in unexplained threshold could reduce the possibility of identifying
sudden death where stringent variant evaluation variants that would be deemed too common in the
results in a significant reduction of numbers of population to cause a rare disease such as LQTS, it
“likely pathogenic” and “pathogenic” variants of could also prevent the identification of potentially
clinical utility (39). important functionally significant variants that
Using a similarly stringent variant analysis, we could play a role in SIDS pathogenesis. For example,
observed previously a 13% (40 of 302) yield of ultra- the p.R176W-KCNH2 variant was identified in a
rare “pathogenic” or “likely pathogenic” variants single European SIDS case in our cohort. This
within sudden death–susceptibility genes among variant has been associated with LQTS previously, it
302 autopsy-negative cases of sudden arrhythmic has been demonstrated to have a functional effect
death syndrome in persons who died at an age >1 by in vitro assays, and it has been considered a
year (median age 24 years), compared with a founder LQT2 mutation in the Finnish population
significantly (p ¼ 0.00002) lower yield of 4.3% (18 (40,41). Although this variant meets the current
of 419) in our SIDS cohort (39). These data suggest American College of Medical Genetics and Genomics
that most SIDS cases stem from pathobiological guideline classification of “pathogenic” variant, its
1226 Tester et al. JACC VOL. 71, NO. 11, 2018
Finnish European individuals) exceeds our stringent Behr, St. George’s, University of London, Cranmer
cutoff and was therefore not included in our Terrace, London SW17 0RE, United Kingdom. E-mail:
analysis. [email protected]. OR Dr. Michael J. Ackerman, Mayo
Clinic, 200 First Street Southwest, Rochester,
CONCLUSIONS Minnesota 55905. E-mail: [email protected].
REFERENCES
1. Krous HF, Beckwith JB, Byard RW, et al. Sudden 9. Yun AJ, Lee PY. Sudden death among infants 16. Glengarry JM, Crawford J, Morrow PL,
infant death syndrome and unclassified sudden and adults: companion disorders of maladaptive Stables SR, Love DR, Skinner JR. Long QT molec-
infant deaths: a definitional and diagnostic sympathetic bias. Med Hypotheses 2004;62: ular autopsy in sudden infant death syndrome.
approach. Pediatrics 2004;114:234–8. 857–60. Arch Dis Child 2014;99:635–40.
2. Matthews T, MacDorman MF. Infant mortality 10. Ackerman MJ, Siu BL, Sturner WQ, et al. 17. Arnestad M, Crotti L, Rognum TO, et al. Prev-
statistics from the 2010 period linked birth/infant Postmortem molecular analysis of SCN5A defects alence of long-QT syndrome gene variants in
death data set. Natl Vital Stat Rep 2013;18:1–26. in sudden infant death syndrome. JAMA 2001; sudden infant death syndrome. Circulation 2007;
3. Unexplained deaths in infancy: England and 286:2264–9. 115:361–7.
Wales: 2009. Office for National Statistics 2011. 18. Tester DJ, Ackerman MJ. Sudden infant
11. Kinney HC, Richerson GB, Dymecki SM,
Available at: http://www.ons.gov.uk/ons/ death syndrome: how significant are the cardiac
Darnall RA, Nattie EE. The brainstem and seroto-
dcp171778_227450.pdf. Accessed January 24, channelopathies? Cardiovasc Res 2005;67:
nin in the sudden infant death syndrome. Annu
2018. 388–96.
Rev Pathol 2009;4:517–50.
4. Moon RY, Horne RSC, Hauck FR. Sudden infant 19. Tester DJ, Dura M, Carturan E, et al.
death syndrome. Lancet 2007;370:1578–87. 12. Weese-Mayer DE, Ackerman MJ, Marazita ML,
A mechanism for sudden infant death syndrome
Berry-Kravis EM. Sudden infant death syndrome:
5. Filiano JJ, Kinney HC. A perspective on neuro- (SIDS): Stress-induced leak via ryanodine re-
review of implicated genetic factors. Am J Med
pathologic findings in victims of the sudden infant ceptors. Heart Rhythm 2007;4:733–9.
Genet 2007;143:771–88.
death syndrome: the triple-risk model. Biol 20. Van Norstrand DW, Ackerman MJ. Sudden in-
Neonate 1994;65:194–7. 13. Van Norstrand DW, Ackerman MJ. Genomic risk fant death syndrome: do ion channels play a role?
factors in sudden infant death syndrome. Genome Heart Rhythm 2009;6:272–8.
6. Valdes-Dapena MA. Editorial: Sudden, unex-
pected and unexplained death in infancy—a status Med 2010;2:86.
21. Brion M, Allegue C, Santori M, et al. Sarcomeric
report—1973. N Engl J Med 1973;289:1195–7. 14. Salomonis N. Systems-level perspective of gene mutations in sudden infant death syndrome
7. Schwartz PJ. Cardiac sympathetic innervation sudden infant death syndrome. Pediatr Res 2014; (SIDS). Forensic Sci Int 2012;219:278–81.
and the sudden infant death syndrome: a possible 76:220–9.
22. Hertz CL, Christiansen SL, Larsen MK, et al.
pathogenic link. Am J Med 1976;60:167–72.
15. Schwartz PJ, Stramba-Badiale M, Segantini A, Genetic investigations of sudden unexpected
8. Guntheroth WG. Theories of cardiovascular et al. Prolongation of the QT interval and the deaths in infancy using next-generation
causes in sudden infant death syndrome. J Am Coll sudden infant death syndrome. N Engl J Med sequencing of 100 genes associated with cardiac
Cardiol 1989;14:443–7. 1998;338:1709–14. diseases. Eur J Hum Genet 2016;24:817–22.
JACC VOL. 71, NO. 11, 2018 Tester et al. 1227
MARCH 20, 2018:1217–27 Cardiac Genetic Analyses in SIDS
23. Cheng J, Van Norstrand DW, Medeiros- admixed populations. Am J Hum Genet 2008; the CASQ2 gene in patients with CPVT: implication
Domingo A, et al. Alpha1-syntrophin mutations 83:132–5. for genetic counselling and clinical management.
identified in sudden infant death syndrome cause Hum Mut 2011;32:995–9.
33. Exome Aggregation Consortium (ExAC), Cam-
an increase in late cardiac sodium current. Circ
bridge, MA. Available at: http://exac. 43. Fressart V, Duthoit G, Donal E, et al. Desmo-
Arrhythm Electrophysiol 2009;2:667–76.
broadinstitute.org. Accessed December 30, 2017. somal gene analysis in arrhythmogenic right ven-
24. Tester DJ, Tan B-H, Medeiros-Domingo A, tricular dysplasia/cardiomyopathy: spectrum of
34. Richards S, Aziz N, Bale S, et al. Standards and
Song C, Makielski JC, Ackerman MJ. Loss-of- mutations and clinical impact in practice. Europace
guidelines for the interpretation of sequence var-
function mutations in the KCNJ8-encoded Kir6.1 2010;12:861–8.
iants: a joint consensus recommendation of the
KATP channel and sudden infant death syndrome. 44. Kapplinger JD, Tester DJ, Salisbury BA, et al.
American College of Medical Genetics and Geno-
Circ Cardiovasc Genet 2011;4:510–5. Spectrum and prevalence of mutations from the
mics and the Association for Molecular Pathology.
25. Van Norstrand DW, Asimaki A, Rubinos C, et al. Genet Med 2015;17:405–24. first 2,500 consecutive unrelated patients referred
Connexin43 mutation causes heterogeneous gap for the FAMILION long QT syndrome genetic test.
35. Ackerman MJ. Genetic purgatory and the car-
junction loss and sudden infant death. Circulation Heart Rhythm 2009;6:1297–303.
diac channelopathies: exposing the variants of
2012;125:474–81. 45. Napolitano C, Priori SG, Schwartz PJ, et al.
uncertain/unknown significance (VUS) issue. Heart
Rhythm 2015;12:2325–31. Genetic testing in the long QT syndrome: devel-
26. Van Norstrand DW, Valdivia CR, Tester DJ,
opment and validation of an efficient approach to
et al. Molecular and functional characteriza- 36. Ackerman JP, Bartos DC, Kapplinger JD,
tion of novel glycerol-3-phosphate dehydro- genotyping in clinical practice. JAMA 2005;294:
Tester DJ, Delisle BP, Ackerman MJ. The promise
2975–80.
genase 1 like gene (GPD1-L) mutations in and peril of precision medicine: phenotyping still
sudden infant death syndrome. Circulation matters most. Mayo Clin Proc 2016;91:1606–16. 46. Crotti L, Lundquist AL, Insolia R, et al. KCNH2-
2007;116:2253–9. K897T is a genetic modifier of latent congenital
37. Lek M, Karczewski KJ, Minikel EV, et al. Anal- long-QT syndrome. Circulation 2005;112:1251–8.
27. Cronk LB, Ye B, Kaku T, et al. Novel mechanism ysis of protein-coding genetic variation in 60,706
for sudden infant death syndrome: Persistent late humans. Nature 2016;536:285–91. 47. Haghighi K, Kolokathis F, Gramolini AO, et al.
sodium current secondary to mutations in A mutation in the human phospholamban gene,
caveolin-3. Heart Rhythm 2007;4:161–6. 38. Neubauer J, Lecca MR, Russo G, et al. Post- deleting arginine 14, results in lethal, hereditary
mortem whole-exome analysis in a large sudden cardiomyopathy. Proc Natl Acad Sci U S A 2006;
28. Tan B-H, Pundi KN, Van Norstrand DW, et al. infant death syndrome cohort with a focus on 103:1388–93.
Sudden infant death syndrome-associated muta- cardiovascular and metabolic genetic diseases. Eur
tions in the sodium channel beta subunits. Heart 48. Medeiros-Domingo A, Bhuiyan ZA, Tester DJ,
J Hum Genet 2017;25:404–9.
Rhythm 2010;7:771–8. et al. The RYR2-Encoded Ryanodine Receptor/
39. Lahrouchi N, Raju H, Lodder EM, et al. Utility Calcium Release Channel in patients diagnosed
29. Ruark E, Münz M, Renwick A, et al. The of post-mortem genetic testing in cases of sudden previously with either catecholaminergic poly-
ICR1000 UK exome series: a resource of gene arrhythmic death syndrome. J Am Coll Cardiol morphic ventricular tachycardia or genotype
variation in an outbred population. F1000Res 2017;69:2134–45. negative, exercise-induced long QT Syndrome: a
2015;4:883. comprehensive open reading frame mutational
40. Marjamaa A, Salomaa V, Newton-Cheh C,
30. McKenna A, Hanna M, Banks E, et al. The et al. High prevalence of four long QT syndrome analysis. J Am Coll Cardiol 2009;54:2065–74.
Genome Analysis Toolkit: a MapReduce framework founder mutations in the Finnish population. Ann
for analyzing next-generation DNA sequencing Med 2009;41:234–40. KEY WORDS genetic heart diseases,
data. Genome Res 2010;20:1297–303.
41. Lahti AL, Kujala VJ, Chapman H, et al. Model molecular autopsy, sudden infant death
31. Combined Annotation Dependent Depletion for long QT syndrome type 2 using human iPS cells syndrome, whole exome sequencing
(CADD) scores. Available at: http://cadd.gs. demonstrates arrhythmogenic characteristics in
washington.edu. Accessed January 24, 2018. cell culture. Dis Model Mech 2012;5:220–30.
A PP END IX For supplemental tables and
32. Price AL, Weale ME, Patterson N, et al. 42. Roux-Buisson N, Rendu J, Denjoy I, et al. figures, please see the online version of this
Long-range LD can confound genome scans in Functional analysis reveals splicing mutations of paper.