Deep sequencing of candidate genes identified 14 variants associated with smoking abstinence in an ethnically diverse sample.

Despite the large public health toll of smoking, genetic studies of smoking cessation have been limited with few discoveries of risk or protective loci. We investigated common and rare variant associations with success in quitting smoking using a cohort from 8 randomized controlled trials involving 2231 participants and a total of 10,020 common and 24,147 rare variants. We identified 14 novel markers including 6 mapping to genes previously related to psychiatric and substance use disorders, 4 of which were protective (CYP2B6 (rs1175607105), HTR3B (rs1413172952; rs1204720503), rs80210037 on chr15), and 2 of which were associated with reduced cessation (PARP15 (rs2173763), SCL18A2 (rs363222)). The others mapped to areas associated with cancer including FOXP1 (rs1288980) and ZEB1 (rs7349). Network analysis identified significant canonical pathways for the serotonin receptor signaling pathway, nicotine and bupropion metabolism, and several related to tumor suppression. Two novel markers (rs6749438; rs6718083) on chr2 are flanked by genes associated with regulation of bodyweight. The identification of novel loci in this study can provide new targets of pharmacotherapy and inform efforts to develop personalized treatments based on genetic profiles.

Variants in genes coding for collagen type IV α-chains are frequent causes of persistent, isolated hematuria during childhood.

BACKGROUND: Children with persistent, isolated microscopic hematuria typically undergo a limited diagnostic workup and are monitored for signs of kidney disease in long-term longitudinal follow-up, which can delay diagnosis and allow disease progression in some cases. METHODS: To determine the clinical utility of genetic screening in this population, we performed targeted genetic testing using a custom, 32-gene next-generation sequencing panel for progressive kidney disease on children referred to the Texas Children's Hospital Pediatric Nephrology clinic for persistent, microscopic hematuria (n = 30; cohort 1). Patients with microscopic hematuria identified by urinalysis on at least two separate occasions were eligible for enrollment, but those with other evidence of kidney disease were excluded. Results were analyzed for sequence variants using the American College of Medical Genetics and Genomics (ACMG) guideline for data interpretation and were validated using a secondary analysis of a dataset of children with hematuria and normal kidney function who had undergone genetic testing as part of an industry-sponsored program (cohort 2; n = 67). RESULTS: In cohort 1 33% of subjects (10/30) had pathogenic or likely pathogenic (P/LP) variants in the type IV collagen genes (COL4A3/A4/A5), and 10% (3/30) had variants of uncertain significance in these genes. The high diagnostic rate in type IV collagen genes was confirmed in cohort 2, where 27% (18/67) of subjects had P/LP variants in COL4A3/A4/A5 genes. CONCLUSIONS: Children with persistent, isolated microscopic hematuria have a high likelihood of having pathogenic variants in type IV collagen genes and genetic screening should be considered. A higher resolution version of the Graphical abstract is available as Supplementary information.

Recurring germline mosaicism in a family due to reversion of an inherited derivative chromosome 8 from an 8;21 translocation with interstitial telomeric sequences.

BACKGROUND: Mosaicism for chromosomal structural abnormalities, other than marker or ring chromosomes, is rarely inherited. METHODS: We performed cytogenetics studies and breakpoint analyses on a family with transmission of mosaicism for a derivative chromosome 8 (der(8)), resulting from an unbalanced translocation between the long arms of chromosomes 8 and 21 over three generations. RESULTS: The proband and his maternal half-sister had mosaicism for a der(8) cell line leading to trisomy of the distal 21q, and both had Down syndrome phenotypic features. Mosaicism for a cell line with the der(8) and a normal cell line was also detected in a maternal half-cousin. The der(8) was inherited from the maternal grandmother who had four abnormal cell lines containing the der(8), in addition to a normal cell line. One maternal half-aunt had the der(8) and an isodicentric chromosome 21 (idic(21)). Sequencing studies revealed microhomologies at the junctures of the der(8) and idic(21) in the half-aunt, suggesting a replicative mechanism in the rearrangement formation. Furthermore, interstitial telomeric sequences (ITS) were identified in the juncture between chromosomes 8 and 21 in the der(8). CONCLUSION: Mosaicism in the proband, his half-sister and half-cousin resulting from loss of chromosome 21 material from the der(8) appears to be a postzygotic event due to the genomic instability of ITS and associated with selective growth advantage of normal cells. The reversion of the inherited der(8) to a normal chromosome 8 in this family resembles revertant mosaicism of point mutations. We propose that ITS could mediate recurring revertant mosaicism for some constitutional chromosomal structural abnormalities.

Implementation of preemptive DNA sequence-based pharmacogenomics testing across a large academic medical center: The Mayo-Baylor RIGHT 10K Study.

PURPOSE: The Mayo-Baylor RIGHT 10K Study enabled preemptive, sequence-based pharmacogenomics (PGx)-driven drug prescribing practices in routine clinical care within a large cohort. We also generated the tools and resources necessary for clinical PGx implementation and identified challenges that need to be overcome. Furthermore, we measured the frequency of both common genetic variation for which clinical guidelines already exist and rare variation that could be detected by DNA sequencing, rather than genotyping. METHODS: Targeted oligonucleotide-capture sequencing of 77 pharmacogenes was performed using DNA from 10,077 consented Mayo Clinic Biobank volunteers. The resulting predicted drug response-related phenotypes for 13 genes, including CYP2D6 and HLA, affecting 21 drug-gene pairs, were deposited preemptively in the Mayo electronic health record. RESULTS: For the 13 pharmacogenes of interest, the genomes of 79% of participants carried clinically actionable variants in 3 or more genes, and DNA sequencing identified an average of 3.3 additional conservatively predicted deleterious variants that would not have been evident using genotyping. CONCLUSION: Implementation of preemptive rather than reactive and sequence-based rather than genotype-based PGx prescribing revealed nearly universal patient applicability and required integrated institution-wide resources to fully realize individualized drug therapy and to show more efficient use of health care resources.

CYP2C8, CYP2C9, and CYP2C19 Characterization Using Next-Generation Sequencing and Haplotype Analysis: A GeT-RM Collaborative Project.

Pharmacogenetic tests typically target selected sequence variants to identify haplotypes that are often defined by star (∗) allele nomenclature. Due to their design, these targeted genotyping assays are unable to detect novel variants that may change the function of the gene product and thereby affect phenotype prediction and patient care. In the current study, 137 DNA samples that were previously characterized by the Genetic Testing Reference Material (GeT-RM) program using a variety of targeted genotyping methods were recharacterized using targeted and whole genome sequencing analysis. Sequence data were analyzed using three genotype calling tools to identify star allele diplotypes for CYP2C8, CYP2C9, and CYP2C19. The genotype calls from next-generation sequencing (NGS) correlated well to those previously reported, except when novel alleles were present in a sample. Six novel alleles and 38 novel suballeles were identified in the three genes due to identification of variants not covered by targeted genotyping assays. In addition, several ambiguous genotype calls from a previous study were resolved using the NGS and/or long-read NGS data. Diplotype calls were mostly consistent between the calling algorithms, although several discrepancies were noted. This study highlights the utility of NGS for pharmacogenetic testing and demonstrates that there are many novel alleles that are yet to be discovered, even in highly characterized genes such as CYP2C9 and CYP2C19.

Accuracy and applications of sequencing and genotyping approaches for CYP2A6 and homologous genes.

OBJECTIVES: We evaluated multiple genotyping/sequencing approaches in a homologous region of chromosome 19, and investigated associations of two common 3'-UTR CYP2A6 variants with activity in vivo. METHODS: Individuals (n = 1704) of European and African ancestry were phenotyped for the nicotine metabolite ratio (NMR), an index of CYP2A6 activity, and genotyped/sequenced using deep amplicon exon sequencing, SNP array, genotype imputation and targeted capture sequencing. Amplicon exon sequencing was the gold standard to which other methods were compared within-individual for CYP2A6, CYP2A7, CYP2A13, and CYP2B6 exons to identify highly discordant positions. Linear regression models evaluated the association of CYP2A6*1B and rs8192733 genotypes (coded additively) with logNMR. RESULTS: All approaches were ≤2.6% discordant with the gold standard; discordant calls were concentrated at few positions. Fifteen positions were discordant in >10% of individuals, with 12 appearing in regions of high identity between homologous genes (e.g. CYP2A6 and CYP2A7). For six, allele frequencies in our study and online databases were discrepant, suggesting errors in online sources. In the European-ancestry group (n = 935), CYP2A6*1B and rs8192733 were associated with logNMR (P < 0.001). A combined model found main effects of both variants on increasing logNMR. Similar trends were found in those of African ancestry (n = 506). CONCLUSION: Multiple genotyping/sequencing approaches used in this chromosome 19 region contain genotyping/sequencing errors, as do online databases. Gene-specific primers and SNP array probes must consider gene homology; short-read sequencing of related genes in a single reaction should be avoided. Using improved sequencing approaches, we characterized two gain-of-function 3'-UTR variants, including the relatively understudied rs8192733.

Targeted Genotyping in Clinical Pharmacogenomics: What Is Missing?

Clinical pharmacogenomic testing typically uses targeted genotyping, which only detects variants included in the test design and may vary among laboratories. To evaluate the potential patient impact of genotyping compared with sequencing, which can detect common and rare variants, an in silico targeted genotyping panel was developed based on the variants most commonly included in clinical tests and applied to a cohort of 10,030 participants who underwent sequencing for CYP1A2, CYP2C19, CYP2C9, CYP2D6, CYP3A4, CYP3A5, DPYD, SLCO1B1, TPMT, UGT1A1, and VKORC1. The results of in silico targeted genotyping were compared with the clinically reported sequencing results. Of the 10,030 participants, 2780 (28%) had at least one potentially clinically relevant variant/allele identified by sequencing that would not have been detected in a standard targeted genotyping panel. The genes with the largest number of participants with variants only detected by sequencing were SLCO1B1, DPYD, and CYP2D6, which affected 13%, 6.3%, and 3.5% of participants, respectively. DPYD (112 variants) and CYP2D6 (103 variants) had the largest number of unique variants detected only by sequencing. Although targeted genotyping detects most clinically significant pharmacogenomic variants, sequencing-based approaches are necessary to detect rare variants that collectively affect many patients. However, efforts to establish pharmacogenomic variant classification systems and nomenclature to accommodate rare variants will be required to adopt sequencing-based pharmacogenomics.

Pharmacogenomics of Nicotine Metabolism: Novel CYP2A6 and CYP2B6 Genetic Variation Patterns in Alaska Native and American Indian Populations.

INTRODUCTION: Alaska Native and American Indian (AN/AI) populations have higher tobacco use prevalence than other ethnic/racial groups. Pharmacogenetic testing to tailor tobacco cessation treatment may improve cessation rates. This study characterized polymorphic variations among AN/AI people in genes associated with metabolism of nicotine and drugs used for tobacco cessation. METHODS: Recruitment of AN/AI individuals represented six subgroups, five geographic subgroups throughout Alaska and a subgroup comprised of AIs from the lower 48 states living in Alaska. We sequenced the CYP2A6 and CYP2B6 genes to identify known and novel gain, reduced, and loss-of-function alleles, including structural variation (eg, gene deletions, duplications, and hybridizations). RESULTS: Variant allele frequencies differed substantially between AN/AI subgroups. The gene deletion CYP2A6*4 and reduced function CYP2A6*9 alleles were found at high frequency in Northern/Western subgroups and in Lower 48/Interior subgroups, respectively. The reduced function CYP2B6*6 allele was observed in all subgroups and a novel, predicted reduced function CYP2B6 variant was found at relatively high frequency in the Southeastern subgroup. CONCLUSIONS: Diverse CYP2A6 and CYP2B6 variation among the subgroups highlight the need for comprehensive pharmacogenetic testing to guide tobacco cessation therapy for AN/AI populations. IMPLICATIONS: Nicotine metabolism is largely determined by CYP2A6 genotype, and variation in CYP2A6 activity has altered the treatment success in other populations. These findings suggest pharmacogenetic-guided smoking cessation drug treatment could provide benefit to this unique population seeking tobacco cessation therapy.

Whole-Organ Genomic Characterization of Mucosal Field Effects Initiating Bladder Carcinogenesis.

We used whole-organ mapping to study the locoregional molecular changes in a human bladder containing multifocal cancer. Widespread DNA methylation changes were identified in the entire mucosa, representing the initial field effect. The field effect was associated with subclonal low-allele frequency mutations and a small number of DNA copy alterations. A founder mutation in the RNA splicing gene, ACIN1, was identified in normal mucosa and expanded clonally with an additional 21 mutations in progression to carcinoma. The patterns of mutations and copy number changes in carcinoma in situ and foci of carcinoma were almost identical, confirming their clonal origins. The pathways affected by the DNA copy alterations and mutations, including the Kras pathway, were preceded by the field changes in DNA methylation, suggesting that they reinforced mechanisms that had already been initiated by methylation. The results demonstrate that DNA methylation can serve as the initiator of bladder carcinogenesis.

Pharmacogenomic Next-Generation DNA Sequencing: Lessons from the Identification and Functional Characterization of Variants of Unknown Significance in CYP2C9 and CYP2C19.

CYP2C9 and CYP2C19 are highly polymorphic pharmacogenes; however, clinically actionable genetic variability in drug metabolism due to these genes has been limited to a few common alleles. The identification and functional characterization of less-common open reading frame sequence variation might help to individualize therapy with drugs that are substrates for the enzymes encoded by these genes. The present study identified seven uncharacterized variants each in CYP2C9 and CYP2C19 using next-generation sequence data for 1013 subjects, and functionally characterized the encoded proteins. Constructs were created and transiently expressed in COS-1 cells for the assay of protein concentration and enzyme activities using fluorometric substrates and liquid chromatography- tandem mass spectrometry with tolbutamide (CYP2C9) and (S)-mephenytoin (CYP2C19) as prototypic substrates. The results were compared with the SIFT, Polyphen, and Provean functional prediction software programs. Cytochrome P450 oxidoreductase (CPR) activities were also determined. Positive correlations were observed between protein content and fluorometric enzyme activity for variants of CYP2C9 (P < 0.05) and CYP2C19 (P < 0.0005). However, CYP2C9 709G>C and CYP2C19 65A>G activities were much lower than predicted based on protein content. Substrate intrinsic clearance values for CYP2C9 218C>T, 343A>C, and CYP2C19 337G>A, 518C>T, 556C>T, and 557G>A were less than 25% of wild-type allozymes. CPR activity levels were similar for all variants. In summary, sequencing of CYP2C9 and CYP2C19 in 1013 subjects identified low-frequency variants that had not previously been functionally characterized. In silico predictions were not always consistent with functional assay results. These observations emphasize the need for high-throughput methods for pharmacogene variant mutagenesis and functional characterization.

Blood pressure signature genes and blood pressure response to thiazide diuretics: results from the PEAR and PEAR-2 studies.

BACKGROUND: Recently, 34 genes had been associated with differential expression relative to blood pressure (BP)/ hypertension (HTN). We hypothesize that some of the genes associated with BP/HTN are also associated with BP response to antihypertensive treatment with thiazide diuretics. METHODS: We assessed these 34 genes for association with differential expression to BP response to thiazide diuretics with RNA sequencing in whole blood samples from 150 hypertensive participants from the Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) and PEAR-2 studies. PEAR white and PEAR-2 white and black participants (n = 50 for each group) were selected based on the upper and lower quartile of BP response to hydrochlorothiazide (HCTZ) and to chlorthalidone. RESULTS: FOS, DUSP1 and PPP1R15A were differentially expressed across all cohorts (meta-analysis p-value < 2.0 × 10- 6), and responders to HCTZ or chlorthalidone presented up-regulated transcripts. Rs11065987 in chromosome 12, a trans-eQTL for expression of FOS, PPP1R15A and other genes, is also associated with BP response to HCTZ in PEAR whites (SBP: ß = - 2.1; p = 1.7 × 10- 3; DBP: ß = - 1.4; p = 2.9 × 10- 3). CONCLUSIONS: These findings suggest FOS, DUSP1 and PPP1R15A as potential molecular determinants of antihypertensive response to thiazide diuretics. TRIAL REGISTRATION: NCT00246519 , NCT01203852 www.clinicaltrials.gov.

Genetic Variants Influencing Plasma Renin Activity in Hypertensive Patients From the PEAR Study (Pharmacogenomic Evaluation of Antihypertensive Responses).

BACKGROUND: Plasma renin is an important regulator of blood pressure (BP). Plasma renin activity (PRA) has been shown to correlate with variability in BP response to antihypertensive agents. We conducted a genome-wide association study to identify single-nucleotide polymorphisms (SNPs) associated with baseline PRA using data from the PEAR study (Pharmacogenomic Evaluation of Antihypertensive Responses). METHODS: Multiple linear regression analysis was performed in 461 whites and 297 blacks using an additive model, adjusting for age, sex, and ancestry-specific principal components. Top SNPs were prioritized by testing the expected direction of association for BP response to atenolol and hydrochlorothiazide. Top regions from the BP response prioritization were tested for functional evidence through differences in gene expression by genotype using RNA sequencing data. Regions with functional evidence were assessed for replication with baseline PRA in an independent study (PEAR-2). RESULTS: Our top SNP rs3784921 was in the SNN-TXNDC11 gene region. The G allele of rs3784921 was associated with higher baseline PRA (ß=0.47; P=2.09×10-6) and smaller systolic BP reduction in response to hydrochlorothiazide (ß=2.97; 1-sided P=0.006). In addition, TXNDC11 expression differed by rs3784921 genotype (P=0.007), and rs1802409, a proxy SNP for rs3784921 (r2=0.98-1.00), replicated in PEAR-2 (ß=0.15; 1-sided P=0.038). Additional SNPs associated with baseline PRA that passed BP response prioritization were in/near the genes CHD9, XIRP2, and GHR. CONCLUSIONS: We identified multiple regions associated with baseline PRA that were prioritized through BP response signals to 2 mechanistically different antihypertensive drugs. CLINICAL TRIAL REGISTRATION: URL: https://www.clinicaltrials.gov. Unique identifier: NCT00246519.

Novel CYP2A6 diplotypes identified through next-generation sequencing are associated with in-vitro and in-vivo nicotine metabolism.

OBJECTIVES: Smoking patterns and cessation rates vary widely across smokers and can be influenced by variation in rates of nicotine metabolism [i.e. cytochrome P450 2A6 (CYP2A6), enzyme activity]. There is high heritability of CYP2A6-mediated nicotine metabolism (60-80%) owing to known and unidentified genetic variation in the CYP2A6 gene. We aimed to identify and characterize additional genetic variants at the CYP2A6 gene locus. METHODS: A new CYP2A6-specific sequencing method was used to investigate genetic variation in CYP2A6. Novel variants were characterized in a White human liver bank that has been extensively phenotyped for CYP2A6. Linkage and haplotype structure for the novel single nucleotide polymorphisms (SNPs) were assessed. The association between novel five-SNP diplotypes and nicotine metabolism rate was investigated. RESULTS: Seven high-frequency (minor allele frequencies ≥6%) noncoding SNPs were identified as important contributors to CYP2A6 phenotypes in a White human liver bank (rs57837628, rs7260629, rs7259706, rs150298687 (also denoted rs4803381), rs56113850, rs28399453, and rs8192733), accounting for two times more variation in in-vitro CYP2A6 activity relative to the four established functional CYP2A6 variants that are frequently tested in Whites (CYP2A6*2, *4, *9, and *12). Two pairs of novel SNPs were in high linkage disequilibrium, allowing us to establish five-SNP diplotypes that were associated with CYP2A6 enzyme activity (rate of nicotine metabolism) in-vitro in the liver bank and in-vivo among smokers. CONCLUSION: The novel five-SNP diplotype may be useful to incorporate into CYP2A6 genotype models for personalized prediction of nicotine metabolism rate, cessation success, and response to pharmacotherapies.

Whole Transcriptome Sequencing Analyses Reveal Molecular Markers of Blood Pressure Response to Thiazide Diuretics.

Thiazide diuretics (TD) are commonly prescribed anti-hypertensives worldwide. However, <40% of patients treated with thiazide monotherapy achieve BP control. This study uses whole transcriptome sequencing to identify novel molecular markers associated with BP response to TD. We assessed global RNA expression levels in whole blood samples from 150 participants, representing patients in the upper and lower quartile of BP response to TD from the Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) (50 whites) and from PEAR-2 (50 whites and 50 blacks). In each study cohort, we performed poly-A RNA-sequencing in baseline samples from 25 responders and 25 non-responders to hydrochlorothiazide (HCTZ) or chlorthalidone. At FDR adjusted p-value < 0.05, 29 genes were differentially expressed in relation to HCTZ or chlorthalidone BP response in whites. For each differentially expressed gene, replication was attempted in the alternate white group and PEAR-2 blacks. CEBPD (meta-analysis p = 1.8 × 10-11) and TSC22D3 (p = 1.9 × 10-9) were differentially expressed in all 3 cohorts, and explain, in aggregate, 21.9% of response variability to TD. This is the first report of the use of transcriptome-wide sequencing data to identify molecular markers of antihypertensive drug response. These findings support CEBPD and TSC22D3 as potential biomarkers of BP response to TD.

Breast cancer chemoprevention pharmacogenomics: Deep sequencing and functional genomics of the ZNF423 and CTSO genes.

Our previous GWAS using samples from the NSABP P-1 and P-2 selective estrogen receptor modulator (SERM) breast cancer prevention trials identified SNPs in ZNF423 and near CTSO that were associated with breast cancer risk during SERM chemoprevention. We have now performed Next Generation DNA sequencing to identify additional SNPs that might contribute to breast cancer risk and to extend our observation that SNPs located hundreds of bp from estrogen response elements (EREs) can alter estrogen receptor alpha (ERα) binding in a SERM-dependent fashion. Our study utilized a nested case-control cohort selected from patients enrolled in the original GWAS, with 199 cases who developed breast cancer during SERM therapy and 201 matched controls who did not. We resequenced approximately 500 kb across both ZNF423 and CTSO, followed by functional genomic studies. We identified 4079 SNPs across ZNF423 and 3876 across CTSO, with 9 SNPs in ZNF423 and 12 in CTSO with p < 1E-02 that were within 500 bp of an ERE motif. The rs746157 (p = 8.44E-04) and rs12918288 SNPs (p = 3.43E-03) in intron 5 of ZNF423, were in linkage equilibrium and were associated with alterations in ER-binding to an ERE motif distant from these SNPs. We also studied all nonsynonymous SNPs in both genes and observed that one nsSNP in ZNF423 displayed decreased protein expression. In conclusion, we identified additional functional SNPs in ZNF423 that were associated with SNP and SERM-dependent alternations in ER binding and transcriptional regulation for an ERE at a distance from the SNPs, thus providing novel insight into mechanisms of SERM effect.

Genome-Wide Prioritization and Transcriptomics Reveal Novel Signatures Associated With Thiazide Diuretics Blood Pressure Response.

BACKGROUND: Thiazide diuretics are among the most commonly prescribed antihypertensives. However, <50% of thiazide-treated patients achieve blood pressure (BP) control. Herein, we used different omics (genomics and transcriptomics) to identify novel biomarkers of thiazide diuretics BP response. METHODS AND RESULTS: Genome-wide analysis included 228 white hypertensives with BP determined at baseline and after 9 weeks of hydrochlorothiazide. Single-nucleotide polymorphisms with P <5×10-5 were prioritized according to their biological function, using RegulomeDB, haploreg, and Genome-Wide Annotation of Variants. The results from the prioritization approach revealed rs10995 as the most likely functional single-nucleotide polymorphism, among single-nucleotide polymorphisms tested, that has been associated with hydrochlorothiazide BP response. The rs10995 G-allele was associated with better BP response to hydrochlorothiazide versus noncarriers (Δ systolic BP/Δ diastolic BP: -12.3/-8.2 versus -6.8/-3.5 mm Hg, respectively, Δ systolic BP P=3×10-4, Δ diastolic BP P=5×10-5). This association was replicated in independent participants treated with chlorthalidone. In addition, rs10995 G-allele was associated with increased mRNA expression of VASP (vasodilator-stimulated phosphoprotein). Moreover, baseline expression of the VASP mRNA was significantly higher in 25 good responders to hydrochlorothiazide compared with 25 poor responders (P=0.01). This finding was replicated in independent participants treated with chlorthalidone (P=0.04). Last, allelic-specific expression analysis revealed a significant but modest imbalance with rs10995 and rs10156, a single-nucleotide polymorphism in high linkage disequilibrium (r2=0.7) with rs10995, which both could contribute to the observed genetic effects by affecting VASP mRNA expression. CONCLUSIONS: This study highlights the strength of using different omics to identify novel biomarkers of drug response and suggests VASP as a potential determinant of thiazide diuretics BP response. CLINICAL TRIAL REGISTRATION: URL: http://www.clinicaltrials.gov. Unique identifiers: NCT00246519 and NCT01203852.

Postmortem genetic screening for the identification, verification, and reporting of genetic variants contributing to the sudden death of the young.

Each year in the United States, thousands of cases of sudden and unexpected deaths of infants, children, and young adults are assigned an undetermined cause of death after postmortem investigation and autopsy. Heritable genetic variants have been suggested as the cause of up to a third of sudden death (SD) cases. Elucidation of the genetic variants involved in SD cases is important to not only help establish cause and manner of death of these individuals, but to also aid in determining whether familial genetic testing should be considered. Previously, these types of postmortem screenings have not been a feasible option for most county medical examiners' and coroners' offices. We sequenced full exons of 64 genes associated with SD in the largest known cohort (351) of infant and young SD decedents using massively parallel sequencing at <$600 per sample. Genetic variants were assessed through literature review and clinical evaluation by a multidisciplinary consortium of experts. Thirteen individuals (3.7%), eight infants (2.8% of those <1 yr of age) and five children/young adults (7.0% of those >1 yr of age), were found to have a reportable genetic variant contributing to SD. These percentages represent an estimate lower than those previously reported. Overall yields and results likely vary between studies due to differences in evaluation techniques and reporting. Additionally, we recommend ongoing assessment of data, including nonreported novel variants, as technology and literature continually advance. This study demonstrates a strategy to implement molecular autopsies in medicolegal investigations of young SD decedents.