Blended Genome Exome (BGE) as a Cost Efficient Alternative to Deep Whole Genomes or Arrays.

Genomic scientists have long been promised cheaper DNA sequencing, but deep whole genomes are still costly, especially when considered for large cohorts in population-level studies. More affordable options include microarrays + imputation, whole exome sequencing (WES), or low-pass whole genome sequencing (WGS) + imputation. WES + array + imputation has recently been shown to yield 99% of association signals detected by WGS. However, a method free from ascertainment biases of arrays or the need for merging different data types that still benefits from deeper exome coverage to enhance novel coding variant detection does not exist. We developed a new, combined, "Blended Genome Exome" (BGE) in which a whole genome library is generated, an aliquot of that genome is amplified by PCR, the exome regions are selected and enriched, and the genome and exome libraries are combined back into a single tube for sequencing (33% exome, 67% genome). This creates a single CRAM with a low-coverage whole genome (2-3x) combined with a higher coverage exome (30-40x). This BGE can be used for imputing common variants throughout the genome as well as for calling rare coding variants. We tested this new method and observed >99% r 2 concordance between imputed BGE data and existing 30x WGS data for exome and genome variants. BGE can serve as a useful and cost-efficient alternative sequencing product for genomic researchers, requiring ten-fold less sequencing compared to 30x WGS without the need for complicated harmonization of array and sequencing data.

CR1 variants contribute to FSGS susceptibility across multiple populations.

Focal segmental glomerulosclerosis (FSGS) is a common cause of nephrotic syndrome with an annual incidence in the United States in African-Americans compared to European-Americans of 24 cases and 5 cases per million, respectively. Among glomerular diseases in Europe and Latin-America, FSGS was the second most frequent diagnosis, and in Asia the fifth. We expand previous efforts in understanding genetics of FSGS by performing a case-control study involving ethnically-diverse groups FSGS cases (726) and a pool of controls (13,994), using panel sequencing of approximately 2,500 podocyte-expressed genes. Through rare variant association tests, we replicated known risk genes - KANK1, COL4A4, and APOL1. A novel significant association was observed for the gene encoding complement receptor 1 (CR1). High-risk rare variants in CR1 in the European-American cohort were commonly observed in Latin- and African-Americans. Therefore, a combined rare and common variant analysis was used to replicate the CR1 association in non-European populations. The CR1 risk variant, rs17047661, gives rise to the Sl1/Sl2 (R1601G) allele that was previously associated with protection against cerebral malaria. Pleiotropic effects of rs17047661 may explain the difference in allele frequencies across continental ancestries and suggest a possible role for genetically-driven alterations of adaptive immunity in the pathogenesis of FSGS.

CHARR efficiently estimates contamination from DNA sequencing data.

DNA sample contamination is a major issue in clinical and research applications of whole-genome and -exome sequencing. Even modest levels of contamination can substantially affect the overall quality of variant calls and lead to widespread genotyping errors. Currently, popular tools for estimating the contamination level use short-read data (BAM/CRAM files), which are expensive to store and manipulate and often not retained or shared widely. We propose a metric to estimate DNA sample contamination from variant-level whole-genome and -exome sequence data called CHARR, contamination from homozygous alternate reference reads, which leverages the infiltration of reference reads within homozygous alternate variant calls. CHARR uses a small proportion of variant-level genotype information and thus can be computed from single-sample gVCFs or callsets in VCF or BCF formats, as well as efficiently stored variant calls in Hail VariantDataset format. Our results demonstrate that CHARR accurately recapitulates results from existing tools with substantially reduced costs, improving the accuracy and efficiency of downstream analyses of ultra-large whole-genome and exome sequencing datasets.

Systematic evaluation of genome sequencing for the diagnostic assessment of autism spectrum disorder and fetal structural anomalies.

Short-read genome sequencing (GS) holds the promise of becoming the primary diagnostic approach for the assessment of autism spectrum disorder (ASD) and fetal structural anomalies (FSAs). However, few studies have comprehensively evaluated its performance against current standard-of-care diagnostic tests: karyotype, chromosomal microarray (CMA), and exome sequencing (ES). To assess the clinical utility of GS, we compared its diagnostic yield against these three tests in 1,612 quartet families including an individual with ASD and in 295 prenatal families. Our GS analytic framework identified a diagnostic variant in 7.8% of ASD probands, almost 2-fold more than CMA (4.3%) and 3-fold more than ES (2.7%). However, when we systematically captured copy-number variants (CNVs) from the exome data, the diagnostic yield of ES (7.4%) was brought much closer to, but did not surpass, GS. Similarly, we estimated that GS could achieve an overall diagnostic yield of 46.1% in unselected FSAs, representing a 17.2% increased yield over karyotype, 14.1% over CMA, and 4.1% over ES with CNV calling or 36.1% increase without CNV discovery. Overall, GS provided an added diagnostic yield of 0.4% and 0.8% beyond the combination of all three standard-of-care tests in ASD and FSAs, respectively. This corresponded to nine GS unique diagnostic variants, including sequence variants in exons not captured by ES, structural variants (SVs) inaccessible to existing standard-of-care tests, and SVs where the resolution of GS changed variant classification. Overall, this large-scale evaluation demonstrated that GS significantly outperforms each individual standard-of-care test while also outperforming the combination of all three tests, thus warranting consideration as the first-tier diagnostic approach for the assessment of ASD and FSAs.

CHARR efficiently estimates contamination from DNA sequencing data.

DNA sample contamination is a major issue in clinical and research applications of whole genome and exome sequencing. Even modest levels of contamination can substantially affect the overall quality of variant calls and lead to widespread genotyping errors. Currently, popular tools for estimating the contamination level use short-read data (BAM/CRAM files), which are expensive to store and manipulate and often not retained or shared widely. We propose a new metric to estimate DNA sample contamination from variant-level whole genome and exome sequence data, CHARR, Contamination from Homozygous Alternate Reference Reads, which leverages the infiltration of reference reads within homozygous alternate variant calls. CHARR uses a small proportion of variant-level genotype information and thus can be computed from single-sample gVCFs or callsets in VCF or BCF formats, as well as efficiently stored variant calls in Hail VDS format. Our results demonstrate that CHARR accurately recapitulates results from existing tools with substantially reduced costs, improving the accuracy and efficiency of downstream analyses of ultra-large whole genome and exome sequencing datasets.

Rare coding variation provides insight into the genetic architecture and phenotypic context of autism.

Some individuals with autism spectrum disorder (ASD) carry functional mutations rarely observed in the general population. We explored the genes disrupted by these variants from joint analysis of protein-truncating variants (PTVs), missense variants and copy number variants (CNVs) in a cohort of 63,237 individuals. We discovered 72 genes associated with ASD at false discovery rate (FDR) ≤ 0.001 (185 at FDR ≤ 0.05). De novo PTVs, damaging missense variants and CNVs represented 57.5%, 21.1% and 8.44% of association evidence, while CNVs conferred greatest relative risk. Meta-analysis with cohorts ascertained for developmental delay (DD) (n = 91,605) yielded 373 genes associated with ASD/DD at FDR ≤ 0.001 (664 at FDR ≤ 0.05), some of which differed in relative frequency of mutation between ASD and DD cohorts. The DD-associated genes were enriched in transcriptomes of progenitor and immature neuronal cells, whereas genes showing stronger evidence in ASD were more enriched in maturing neurons and overlapped with schizophrenia-associated genes, emphasizing that these neuropsychiatric disorders may share common pathways to risk.

Large-scale sequencing identifies multiple genes and rare variants associated with Crohn's disease susceptibility.

Genome-wide association studies (GWASs) have identified hundreds of loci associated with Crohn's disease (CD). However, as with all complex diseases, robust identification of the genes dysregulated by noncoding variants typically driving GWAS discoveries has been challenging. Here, to complement GWASs and better define actionable biological targets, we analyzed sequence data from more than 30,000 patients with CD and 80,000 population controls. We directly implicate ten genes in general onset CD for the first time to our knowledge via association to coding variation, four of which lie within established CD GWAS loci. In nine instances, a single coding variant is significantly associated, and in the tenth, ATG4C, we see additionally a significantly increased burden of very rare coding variants in CD cases. In addition to reiterating the central role of innate and adaptive immune cells as well as autophagy in CD pathogenesis, these newly associated genes highlight the emerging role of mesenchymal cells in the development and maintenance of intestinal inflammation.

Case Report: Supernormal Vascular Aging in Leningrad Siege Survivors.

Age-related changes in the vascular system play an important role in the biological age and lifespan of a person and maybe affected from an early age onward. One of the indicators of changes in the vascular system is arterial wall stiffness and its main measure, i.e., carotid-femoral pulse wave velocity (cfPWV). We examined arterial wall stiffness in a sample of 305 Leningrad Siege survivors to assess how hunger and stressful conditions during fetal development and early childhood affected the state of the cardiovascular system at a later age and what factors may neutralize the negative impact sustained in early childhood. Here, we presented an evaluation of two unique patients with supernormal vascular aging (SUPERNOVA) phenotype from this cohort and described the details of congruence between hereditary resistance and practiced lifestyle yielding slower biological aging rate.

Genotype imputation and polygenic score estimation in northwestern Russian population.

Numerous studies demonstrated the lack of transferability of polygenic score (PGS) models across populations and the problem arising from unequal presentation of ancestries across genetic studies. However, even within European ancestry there are ethnic groups that are rarely presented in genetic studies. For instance, Russians, being one of the largest, diverse, and yet understudied group in Europe. In this study, we evaluated the reliability of genotype imputation for the Russian cohort by testing several commonly used imputation reference panels (e.g. HRC, 1000G, HGDP). HRC, in comparison with two other panels, showed the most accurate results based on both imputation accuracy and allele frequency concordance between masked and imputed genotypes. We built polygenic score models based on GWAS results from the UK biobank, measured the explained phenotypic variance in the Russian cohort attributed to polygenic scores for 11 phenotypes, collected in the clinic for each participant, and finally explored the role of allele frequency discordance between the UK biobank and the study cohort in the resulting PGS performance.

Thirteen Cases of Valeryl Fentanyl in Michigan: A Call for Expanding Opioid Testing.

ABSTRACT: In this report, we describe 13 cases of drug overdose in Michigan in which valeryl fentanyl was found in postmortem blood. Valeryl fentanyl is a schedule I opioid that is rarely found in drug overdoses in the United States. Although little data exist on the mortality and morbidity associated with valeryl fentanyl, its molecular structure indicates that it would be less potent than fentanyl.When analyzing blood samples for valeryl fentanyl, samples from peripheral sites were sometimes negative for quantitative levels; however, samples from central sites in the same decedent were positive. This could indicate unique pharmacokinetics for valeryl fentanyl, which could have implications for other fentanyl analogs. Given the paucity of pharmacodynamic information, the prohibition of its use, the potential to buttress law enforcement efforts in monitoring drug trafficking trends, and to determine the efficacy of current regulations, laboratories should test for valeryl fentanyl. When testing for valeryl fentanyl, and likely other fentanyl analogs, the site of sample collection is important: central sources of blood are preferred to peripheral sources.

Longitudinal Opioid Surveillance Project Involving Toxicologic Analysis of Postmortem Specimens from 9 Counties in Michigan Suggests the Discovery of New High-Intensity Drug Trafficking Areas.

ABSTRACT: Acetyl fentanyl (AF) is a Schedule I fentanyl analog that has been increasingly seen in heroin and fentanyl polydrug toxicity overdoses in Michigan (MI). Drug users are often unaware of the presence of AF in their drugs because it is often sold mixed into or disguised as heroin. High levels of AF in heroin drug products can cause increased incidence of overdose. This article describes data from a longitudinal opioid surveillance program and details 102 decedents in MI who were found to have evidence of heroin in their postmortem blood. A large portion of these decedents were also found to have evidence of fentanyl and AF. Our data further show significant overlap in incidence rates of AF and heroin-related overdose deaths in several MI counties, suggesting that AF is becoming enmeshed in heroin trafficking. Furthermore, we report unprecedented high incidence rates of AF and heroin-related overdose deaths in Calhoun county, and we propose that it is a high-intensity drug trafficking area. Highways US-131 and US-31 are likely used to transport these drugs. More study is needed into the drug trafficking trends in MI to ascertain drug sources and monitor the ever developing and dangerous polydrug heroin combinations.

Exome sequencing in patient-parent trios suggests new candidate genes for early-onset primary sclerosing cholangitis.

BACKGROUND & AIMS: Primary sclerosing cholangitis (PSC) is a rare bile duct disease strongly associated with inflammatory bowel disease (IBD). Whole-exome sequencing (WES) has contributed to understanding the molecular basis of very early-onset IBD, but rare protein-altering genetic variants have not been identified for early-onset PSC. We performed WES in patients diagnosed with PSC ≤ 12 years to investigate the contribution of rare genetic variants to early-onset PSC. METHODS: In this multicentre study, WES was performed on 87 DNA samples from 29 patient-parent trios with early-onset PSC. We selected rare (minor allele frequency < 2%) coding and splice-site variants that matched recessive (homozygous and compound heterozygous variants) and dominant (de novo) inheritance in the index patients. Variant pathogenicity was predicted by an in-house developed algorithm (GAVIN), and PSC-relevant variants were selected using gene expression data and gene function. RESULTS: In 22 of 29 trios we identified at least 1 possibly pathogenic variant. We prioritized 36 genes, harbouring a total of 54 variants with predicted pathogenic effects. In 18 genes, we identified 36 compound heterozygous variants, whereas in the other 18 genes we identified 18 de novo variants. Twelve of 36 candidate risk genes are known to play a role in transmembrane transport, adaptive and innate immunity, and epithelial barrier function. CONCLUSIONS: The 36 candidate genes for early-onset PSC need further verification in other patient cohorts and evaluation of gene function before a causal role can be attributed to its variants.

Whole exome sequencing analyses reveal gene-microbiota interactions in the context of IBD.

OBJECTIVE: Both the gut microbiome and host genetics are known to play significant roles in the pathogenesis of IBD. However, the interaction between these two factors and its implications in the aetiology of IBD remain underexplored. Here, we report on the influence of host genetics on the gut microbiome in IBD. DESIGN: To evaluate the impact of host genetics on the gut microbiota of patients with IBD, we combined whole exome sequencing of the host genome and whole genome shotgun sequencing of 1464 faecal samples from 525 patients with IBD and 939 population-based controls. We followed a four-step analysis: (1) exome-wide microbial quantitative trait loci (mbQTL) analyses, (2) a targeted approach focusing on IBD-associated genomic regions and protein truncating variants (PTVs, minor allele frequency (MAF) >5%), (3) gene-based burden tests on PTVs with MAF <5% and exome copy number variations (CNVs) with site frequency <1%, (4) joint analysis of both cohorts to identify the interactions between disease and host genetics. RESULTS: We identified 12 mbQTLs, including variants in the IBD-associated genes IL17REL, MYRF, SEC16A and WDR78. For example, the decrease of the pathway acetyl-coenzyme A biosynthesis, which is involved in short chain fatty acids production, was associated with variants in the gene MYRF (false discovery rate <0.05). Changes in functional pathways involved in the metabolic potential were also observed in participants carrying rare PTVs or CNVs in CYP2D6, GPR151 and CD160 genes. These genes are known for their function in the immune system. Moreover, interaction analyses confirmed previously known IBD disease-specific mbQTLs in TNFSF15. CONCLUSION: This study highlights that both common and rare genetic variants affecting the immune system are key factors in shaping the gut microbiota in the context of IBD and pinpoints towards potential mechanisms for disease treatment.

Prevalence and Effect of Genetic Risk of Thromboembolic Disease in Inflammatory Bowel Disease.

BACKGROUND AND AIMS: The largest cause of mortality in patients with inflammatory bowel disease (IBD) remains thromboembolic disease (TED). Recent reports have demonstrated that both monogenic and polygenic factors contribute to TED and 10% of healthy subjects are genetically at high risk for TED. Our aim was to utilize whole-exome sequencing and genome-wide genotyping to determine the proportion of IBD patients genetically at risk for TED and investigate the effect of genetic risk of TED in IBD. METHODS: The TED polygenic risk score was calculated from genome-wide genotyping. Thrombophilia pathogenic variants were extracted from whole-exome sequencing. In total, 792 IBD patients had both whole-exome sequencing and genotyping data. We defined patients at genetically high risk for TED if they had a high TED polygenic risk score or carried at least 1 thrombophilia pathogenic variant. RESULTS: We identified 122 of 792 IBD patients (15.4%) as genetically high risk for TED. Among 715 of 792 subjects whose documented TED status were available, 63 of the 715 patients (8.8%) had TED events. Genetic TED risk was significantly associated with increased TED event (odds ratio, 2.5; P = .0036). In addition, we confirmed an additive effect of monogenic and polygenic risk on TED (P = .0048). Patients with high TED genetic risk more frequently had thrombosis at multiple sites (78% vs 42%, odds ratio, 3.96; P = .048). CONCLUSIONS: Genetic risk (both poly- and monogenic) was significantly associated with TED history. Our results suggest that genetic traits identify approximately 1 in 7 patients with IBD who will experience 2.5-fold or greater risk for TED.

A structural variation reference for medical and population genetics.

Structural variants (SVs) rearrange large segments of DNA1 and can have profound consequences in evolution and human disease2,3. As national biobanks, disease-association studies, and clinical genetic testing have grown increasingly reliant on genome sequencing, population references such as the Genome Aggregation Database (gnomAD)4 have become integral in the interpretation of single-nucleotide variants (SNVs)5. However, there are no reference maps of SVs from high-coverage genome sequencing comparable to those for SNVs. Here we present a reference of sequence-resolved SVs constructed from 14,891 genomes across diverse global populations (54% non-European) in gnomAD. We discovered a rich and complex landscape of 433,371 SVs, from which we estimate that SVs are responsible for 25-29% of all rare protein-truncating events per genome. We found strong correlations between natural selection against damaging SNVs and rare SVs that disrupt or duplicate protein-coding sequence, which suggests that genes that are highly intolerant to loss-of-function are also sensitive to increased dosage6. We also uncovered modest selection against noncoding SVs in cis-regulatory elements, although selection against protein-truncating SVs was stronger than all noncoding effects. Finally, we identified very large (over one megabase), rare SVs in 3.9% of samples, and estimate that 0.13% of individuals may carry an SV that meets the existing criteria for clinically important incidental findings7. This SV resource is freely distributed via the gnomAD browser8 and will have broad utility in population genetics, disease-association studies, and diagnostic screening.

Large-Scale Exome Sequencing Study Implicates Both Developmental and Functional Changes in the Neurobiology of Autism.

We present the largest exome sequencing study of autism spectrum disorder (ASD) to date (n = 35,584 total samples, 11,986 with ASD). Using an enhanced analytical framework to integrate de novo and case-control rare variation, we identify 102 risk genes at a false discovery rate of 0.1 or less. Of these genes, 49 show higher frequencies of disruptive de novo variants in individuals ascertained to have severe neurodevelopmental delay, whereas 53 show higher frequencies in individuals ascertained to have ASD; comparing ASD cases with mutations in these groups reveals phenotypic differences. Expressed early in brain development, most risk genes have roles in regulation of gene expression or neuronal communication (i.e., mutations effect neurodevelopmental and neurophysiological changes), and 13 fall within loci recurrently hit by copy number variants. In cells from the human cortex, expression of risk genes is enriched in excitatory and inhibitory neuronal lineages, consistent with multiple paths to an excitatory-inhibitory imbalance underlying ASD.

Autism spectrum disorder and attention deficit hyperactivity disorder have a similar burden of rare protein-truncating variants.

The exome sequences of approximately 8,000 children with autism spectrum disorder (ASD) and/or attention deficit hyperactivity disorder (ADHD) and 5,000 controls were analyzed, finding that individuals with ASD and individuals with ADHD had a similar burden of rare protein-truncating variants in evolutionarily constrained genes, both significantly higher than controls. This motivated a combined analysis across ASD and ADHD, identifying microtubule-associated protein 1A (MAP1A) as a new exome-wide significant gene conferring risk for childhood psychiatric disorders.

A complicated case of bowel obstruction with sepsis and methamphetamine toxicity in a child with pica.

In this report, a pediatric case of bowel obstruction with sepsis complicated by methamphetamine toxicity is described. The decedent, an eleven-year-old female with a clinical history of pica, was found unresponsive in her home and pronounced dead following unsuccessful resuscitative efforts. Radiologic imaging showed multiple radio-opaque foreign objects in the stomach and bowel. Autopsy revealed a green leafy substance, coins and other metallic items, folded paper, and plastic in her stomach and bowels. Postmortem iliac blood and urine tested positive for amphetamine and methamphetamine. While the decedent's medical history and autopsy findings provided evidence consistent with bowel obstruction with sepsis due to the ingestion of foreign materials, the high methamphetamine concentration was suggestive of concurrent methamphetamine toxicity. Unique complications associated with this case include the phenomenon that methamphetamine toxicity and bowel obstruction can present similarly in children and the reported opinion that accidental drug ingestion is uncommon in children over the age of five. This case emphasizes that the age range for suspected accidental drug ingestion should be expanded for those with pica, as these patients, despite being older, may not be able to differentiate between what they should and should not ingest. Furthermore, when treating a pediatric patient with pica that appears to present with bowel obstruction, unintentional drug ingestion should also be considered, particularly if there is a suspicion that the child lives in a household where drugs are abused, given the prospect that drug toxicity can present similarly.