Search

Introduction: Accurate tools to inform individual prognosis in patients with autosomal dominant polycystic kidney disease (ADPKD) are lacking. Here, we report an artificial intelligence (AI)-generated method for routinely measuring total kidney volume (TKV)., Methods: An ensemble U-net algorithm was created using the nnUNet approach. The training and internal cross-validation cohort consisted of all 1.5T magnetic resonance imaging (MRI) data acquired using 5 different MRI scanners (454 kidneys, 227 scans) in the CYSTic consortium, which was first manually segmented by a single human operator. As an independent validation cohort, we utilized 48 sequential clinical MRI scans with reference results of manual segmentation acquired by 6 individual analysts at a single center. The tool was then implemented for clinical use and its performance analyzed., Results: The training or internal validation cohort was younger (mean age 44.0 vs. 51.5 years) and the female-to-male ratio higher (1.2 vs. 0.94) compared to the clinical validation cohort. The majority of CYSTic patients had PKD1 mutations (79%) and typical disease (Mayo Imaging class 1, 86%). The median DICE score on the clinical validation data set between the algorithm and human analysts was 0.96 for left and right kidneys with a median TKV error of -1.8%. The time taken to manually segment kidneys in the CYSTic data set was 56 (±28) minutes, whereas manual corrections of the algorithm output took 8.5 (±9.2) minutes per scan., Conclusion: Our AI-based algorithm demonstrates performance comparable to manual segmentation. Its rapidity and precision in real-world clinical cases demonstrate its suitability for clinical application., (Crown Copyright © 2023 Published by Elsevier Inc. on behalf of the International Society of Nephrology.)

Objective: To test the performance of a novel machine learning-based breast density tool. The tool utilizes a convolutional neural network to predict the BI-RADS based density assessment of a study. The clinical density assessments of 33,000 mammographic examinations (164,000 images) from one academic medical center (Site A) were used for training., Materials and Methods: This was an IRB approved HIPAA compliant study performed at two academic medical centers. The validation data set was composed of 500 studies from one site (Site A) and 700 from another (Site B). At Site A, each study was assessed by three breast radiologists and the majority (consensus) assessment was used as truth. At Site B, if the tool agreed with the clinical reading, then it was considered to have correctly predicted the clinical reading. In cases where the tool and the clinical reading disagreed, then the study was evaluated by three radiologists and the consensus reading was used as the clinical reading., Results: For the classification into the four categories of the Breast Imaging Reporting and Data System (BI-RADS®), the AI classifier had an accuracy of 84.6% at Site A and 89.7% at Site B. For binary classification (dense vs. non-dense), the AI classifier had an accuracy of 94.4% at Site A and 97.4% at Site B. In no case did the classifier disagree with the consensus reading by more than one category., Conclusions: The automated breast density tool showed high agreement with radiologists' assessments of breast density., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Sven Schoenherr, Martin Seebass, MingDe Lin and Malte Westerhoff are employees of Visage Imaging, owner and potential distributor of the described software tool., (Copyright © 2023 Elsevier Inc. All rights reserved.)

The benefits of large-scale genetic studies for healthcare of the populations studied are well documented, but these genetic studies have traditionally ignored people from some parts of the world, such as South Asia. Here we describe whole genome sequence (WGS) data from 4806 individuals recruited from the healthcare delivery systems of Pakistan, India and Bangladesh, combined with WGS from 927 individuals from isolated South Asian populations. We characterize population structure in South Asia and describe a genotyping array (SARGAM) and imputation reference panel that are optimized for South Asian genomes. We find evidence for high rates of reproductive isolation, endogamy and consanguinity that vary across the subcontinent and that lead to levels of rare homozygotes that reach 100 times that seen in outbred populations. Founder effects increase the power to associate functional variants with disease processes and make South Asia a uniquely powerful place for population-scale genetic studies., (© 2023. The Author(s).)

Background: Lipoprotein(a) (Lp(a)) concentrations are a major independent risk factor for coronary artery disease (CAD) and are mainly determined by variation in LPA. Up to 70% of the LPA coding sequence is located in the hypervariable kringle IV type 2 (KIV-2) region. It is hardly accessible by conventional technologies, but may contain functional variants., Objectives: This study sought to investigate the new, very frequent splicing variant KIV-2 4733G>A on Lp(a) and CAD., Methods: We genotyped 4733G>A in the GCKD (German Chronic Kidney Disease) study (n = 4,673) by allele-specific polymerase chain reaction, performed minigene assays, identified proxy single nucleotide polymorphisms and used them to characterize its effect on CAD by survival analysis in UK Biobank (n = 440,234). Frequencies in ethnic groups were assessed in the 1000 Genomes Project., Results: The 4733G>A variant (38.2% carrier frequency) was found in most isoform sizes. It reduces allelic expression without abolishing protein production, lowers Lp(a) by 13.6 mg/dL (95% CI: 12.5-14.7; P < 0.0001) and is the strongest variance-explaining factor after the smaller isoform. Splicing of minigenes was modified. Compound heterozygosity (4.6% of the population) for 4733G>A and 4925G>A, another KIV-2 splicing mutation, reduces Lp(a) by 31.8 mg/dL and most importantly narrows the interquartile range by 9-fold (from 42.1 to 4.6 mg/dL) when compared to the wild type. In UK Biobank 4733G>A alone and compound heterozygosity with 4925G>A reduced HR for CAD by 9% (95% CI: 7%-11%) and 12% (95% CI: 7%-16%) (both P < 0.001). Frequencies in ethnicities differ notably., Conclusions: Functional variants in the previously inaccessible LPA KIV-2 region cooperate in determining Lp(a) variance and CAD risk. Even a moderate but lifelong genetic Lp(a) reduction translates to a noticeable CAD risk reduction., Competing Interests: Funding Support and Author Disclosures The study was supported by the Austrian Science Fund (FWF) project P31458-B34 and the D•A•CH Advancement Award Lipidology 2015 (supported by the Christine Katharine Schmitz Foundation) of the D•A•CH-Society for Prevention of Cardiovascular Diseases (to S.C.). The GCKD study is supported by the German Ministry of Education and Research (Bundesministerium für Bildung und Forschung, grants 01ER 0804, 01ER 0818, 01ER 0819, 01ER 0820, and 01ER 0821) and the KfH Foundation for Preventive Medicine (Kuratorium für Heimdialyse und Nierentransplantation e.V.–Stiftung Präventivmedizin) and corporate sponsors (see the GCKD website). Ms Schachtl-Riess has received support from the Dr Legerlotz Foundation. Dr Köttgen has received support from the German Research Foundation (grant KO 3598/5-1). Dr Kronenberg has received support from the Austrian Science Fund (project W-1253DK HOROS); has received lecture fees from Novartis, Amgen, and Kaneka; and has served on the advisory boards of Amgen and Kaneka. Drs Kronenberg and Coassin have received support from the Lipoprotein(a) Center And Research InstitutE [Lp(a)CARE] for their lipoprotein(a) research. All other authors have reported that they have no relationships relevant to the contents of this paper to disclose., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: SARS-CoV-2 entry in human cells depends on angiotensin-converting enzyme 2, which can be upregulated by inhibitors of the renin-angiotensin system (RAS). We aimed to test our hypothesis that discontinuation of chronic treatment with ACE-inhibitors (ACEIs) or angiotensin II receptor blockers (ARBs) mitigates the course o\f recent-onset COVID-19., Methods: ACEI-COVID was a parallel group, randomised, controlled, open-label trial done at 35 centres in Austria and Germany. Patients aged 18 years and older were enrolled if they presented with recent symptomatic SARS-CoV-2 infection and were chronically treated with ACEIs or ARBs. Patients were randomly assigned 1:1 to discontinuation or continuation of RAS inhibition for 30 days. Primary outcome was the maximum sequential organ failure assessment (SOFA) score within 30 days, where death was scored with the maximum achievable SOFA score. Secondary endpoints were area under the death-adjusted SOFA score (AUC SOFA ), mean SOFA score, admission to the intensive care unit, mechanical ventilation, and death. Analyses were done on a modified intention-to-treat basis. This trial is registered with ClinicalTrials.gov, NCT04353596., Findings: Between April 20, 2020, and Jan 20, 2021, 204 patients (median age 75 years [IQR 66-80], 37% females) were randomly assigned to discontinue (n=104) or continue (n=100) RAS inhibition. Within 30 days, eight (8%) of 104 died in the discontinuation group and 12 (12%) of 100 patients died in the continuation group (p=0·42). There was no significant difference in the primary endpoint between the discontinuation and continuation group (median [IQR] maximum SOFA score 0·00 (0·00-2·00) vs 1·00 (0·00-3·00); p=0·12). Discontinuation was associated with a significantly lower AUC SOFA (0·00 [0·00-9·25] vs 3·50 [0·00-23·50]; p=0·040), mean SOFA score (0·00 [0·00-0·31] vs 0·12 [0·00-0·78]; p=0·040), and 30-day SOFA score (0·00 [10-90th percentile, 0·00-1·20] vs 0·00 [0·00-24·00]; p=0·023). At 30 days, 11 (11%) in the discontinuation group and 23 (23%) in the continuation group had signs of organ dysfunction (SOFA score ≥1) or were dead (p=0·017). There were no significant differences for mechanical ventilation (10 (10%) vs 8 (8%), p=0·87) and admission to intensive care unit (20 [19%] vs 18 [18%], p=0·96) between the discontinuation and continuation group., Interpretation: Discontinuation of RAS-inhibition in COVID-19 had no significant effect on the maximum severity of COVID-19 but may lead to a faster and better recovery. The decision to continue or discontinue should be made on an individual basis, considering the risk profile, the indication for RAS inhibition, and the availability of alternative therapies and outpatient monitoring options., Funding: Austrian Science Fund and German Center for Cardiovascular Research., Competing Interests: Declaration of interests AB received research funding from Pfizer and Medtronic as well as speaker honoraria from Bayer, Boehringer Ingelheim, Edwards, Medtronic and Novartis. KDR received a research grant from Daiichi-Sankyo Europe GmbH. UM received fees for participation in a data and safety monitoring board from Thermosome. All other authors declare no competing interests., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Within-species contamination is a major issue in sequencing studies, especially for mitochondrial studies. Contamination can be detected by analyzing the nuclear genome or by inspecting polymorphic sites in the mitochondrial genome (mtDNA). Existing methods using the nuclear genome are computationally expensive, and no appropriate tool for detecting sample contamination in large-scale mtDNA data sets is available. Here we present haplocheck, a tool that requires only the mtDNA to detect contamination in both targeted mitochondrial and whole-genome sequencing studies. Our in silico simulations and amplicon mixture experiments indicate that haplocheck detects mtDNA contamination accurately and is independent of the phylogenetic distance within a sample mixture. By applying haplocheck to The 1000 Genomes Project Consortium data, we further evaluate the application of haplocheck as a fast proxy tool for nDNA-based contamination detection using the mtDNA and identify the mitochondrial copy number within a mixture as a critical component for the overall accuracy. The haplocheck tool is available both as a command-line tool and as a cloud web service producing interactive reports that facilitates the navigation through the phylogeny of contaminated samples., (© 2021 Weissensteiner et al.; Published by Cold Spring Harbor Laboratory Press.)

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes) 1 . In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Pharmaceutical drugs targeting dyslipidemia and cardiovascular disease (CVD) may increase the risk of fatty liver disease and other metabolic disorders. To identify potential novel CVD drug targets without these adverse effects, we perform genome-wide analyses of participants in the HUNT Study in Norway (n = 69,479) to search for protein-altering variants with beneficial impact on quantitative blood traits related to cardiovascular disease, but without detrimental impact on liver function. We identify 76 (11 previously unreported) presumed causal protein-altering variants associated with one or more CVD- or liver-related blood traits. Nine of the variants are predicted to result in loss-of-function of the protein. This includes ZNF529:p.K405X, which is associated with decreased low-density-lipoprotein (LDL) cholesterol (P = 1.3 × 10 -8 ) without being associated with liver enzymes or non-fasting blood glucose. Silencing of ZNF529 in human hepatoma cells results in upregulation of LDL receptor and increased LDL uptake in the cells. This suggests that inhibition of ZNF529 or its gene product should be prioritized as a novel candidate drug target for treating dyslipidemia and associated CVD.

Background: The concentrations of the highly atherogenic lipoprotein(a) [Lp(a)] are mainly genetically determined by the LPA gene locus. However, up to 70% of the coding sequence is located in the complex so-called kringle IV type 2 (KIV-2) copy number variation, a region hardly accessible by common genotyping and sequencing technologies. Despite its size, little is known about genetic variants in this complex region. The R21X variant is a functional variant located in this region, but it has never been analyzed in large cohorts., Methods: We typed R21X in 10,910 individuals from three European populations using a newly developed high-throughput allele-specific qPCR assay. R21X allelic location was determined by separating the LPA alleles using pulsed-field gel electrophoresis (PFGE) and typing them separately. Using GWAS data, we identified a proxy SNP located outside of the KIV-2. Linkage disequilibrium was determined both statistically and by long-range haplotyping using PFGE. Worldwide frequencies were determined by reanalyzing the sequencing data of the 1000 Genomes Project with a dedicated pipeline., Results: R21X carriers (frequency 0.016-0.021) showed significantly lower mean Lp(a) concentrations (- 11.7 mg/dL [- 15.5; - 7.82], p = 3.39e-32). The variant is located mostly on medium-sized LPA alleles. In the 1000 Genome data, R21X mostly occurs in Europeans and South Asians, is absent in Africans, and shows varying frequencies in South American populations (0 to 0.022). Of note, the best proxy SNP was another LPA null mutation (rs41272114, D' = 0.958, R 2  = 0.281). D' was very high in all 1000G populations (0.986-0.996), although rs41272114 frequency varies considerably (0-0.182). Co-localization of both null mutations on the same allele was confirmed by PFGE-based long-range haplotyping., Conclusions: We performed the largest epidemiological study on an LPA KIV-2 variant so far, showing that it is possible to assess LPA KIV-2 mutations on a large scale. Surprisingly, in all analyzed populations, R21X was located on the same haplotype as the splice mutation rs41272114, creating "double-null" LPA alleles. Despite being a nonsense variant, the R21X status does not provide additional information beyond the rs41272114 genotype. This has important implications for studies using LPA loss-of-function mutations as genetic instruments and emphasizes the complexity of LPA genetics.

Haplotype phasing is a fundamental problem in medical and population genetics. Phasing is generally performed via statistical phasing in a genotyped cohort, an approach that can yield high accuracy in very large cohorts but attains lower accuracy in smaller cohorts. Here we instead explore the paradigm of reference-based phasing. We introduce a new phasing algorithm, Eagle2, that attains high accuracy across a broad range of cohort sizes by efficiently leveraging information from large external reference panels (such as the Haplotype Reference Consortium; HRC) using a new data structure based on the positional Burrows-Wheeler transform. We demonstrate that Eagle2 attains a ∼20× speedup and ∼10% increase in accuracy compared to reference-based phasing using SHAPEIT2. On European-ancestry samples, Eagle2 with the HRC panel achieves >2× the accuracy of 1000 Genomes-based phasing. Eagle2 is open source and freely available for HRC-based phasing via the Sanger Imputation Service and the Michigan Imputation Server.

We describe a reference panel of 64,976 human haplotypes at 39,235,157 SNPs constructed using whole-genome sequence data from 20 studies of predominantly European ancestry. Using this resource leads to accurate genotype imputation at minor allele frequencies as low as 0.1% and a large increase in the number of SNPs tested in association studies, and it can help to discover and refine causal loci. We describe remote server resources that allow researchers to carry out imputation and phasing consistently and efficiently.