•Biparental inheritance pattern described by Luo et al. enters in sharp conflict with available evidence.•Recent claims on mtDNA lack methodological support.•An extraordinary unknown molecular ...mechanism would be necessary to avoid mtDNA multi-haploid cells.•Artificial scenarios can be conceived that mimic biparental inheritance patterns.
A breakthrough article published in PNAS by Luo et al. challenges a central dogma in biology which states that the mitochondrial DNA (mtDNA) in humans is inherited exclusively from the mother. We re-analyzed original FASTQ files and results reported by Luo et al. to investigate methodological issues (e.g. nuclear mitochondrial DNA or NUMTs, DNA rearrangements) that could lead to biological misinterpretations. A comprehensive analysis of their data reveals several methodological and analytical issues that must be carefully addressed before challenging the current paradigm. We first show that the probability of the findings described by the authors is extremely small (most likely below 10−37). The sequencing replicates from the same donors show aberrations in the variants detected that need further investigation to exclude contributions from other sources or methodological artifacts. Applying the principle of reductio ad absurdum, we demonstrate that the nuclear factor invoked by the authors to explain the phenomenon would need to be extraordinarily complex and precise to preclude linear accumulation of mtDNA lineages across generations, which would make the appearance of mixed haplotypes a much more frequent event in the population. We discuss alternate scenarios that explain findings of the same nature as reported by Luo et al., in the context of in-vitro fertilization and therapeutic mtDNA replacement ooplasmic transplantation.
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.
This study sought to investigate the new, very frequent splicing variant KIV-2 4733G>A on Lp(a) and CAD.
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.
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.
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.
The human leukocyte antigen (HLA) locus is associated with more complex diseases than any other locus in the human genome. In many diseases, HLA explains more heritability than all other known loci ...combined. In silico HLA imputation methods enable rapid and accurate estimation of HLA alleles in the millions of individuals that are already genotyped on microarrays. HLA imputation has been used to define causal variation in autoimmune diseases, such as type I diabetes, and in human immunodeficiency virus infection control. However, there are few guidelines on performing HLA imputation, association testing, and fine mapping. Here, we present a comprehensive tutorial to impute HLA alleles from genotype data. We provide detailed guidance on performing standard quality control measures for input genotyping data and describe options to impute HLA alleles and amino acids either locally or using the web-based Michigan Imputation Server, which hosts a multi-ancestry HLA imputation reference panel. We also offer best practice recommendations to conduct association tests to define the alleles, amino acids, and haplotypes that affect human traits. Along with the pipeline, we provide a step-by-step online guide with scripts and available software ( https://github.com/immunogenomics/HLA_analyses_tutorial ). This tutorial will be broadly applicable to large-scale genotyping data and will contribute to defining the role of HLA in human diseases across global populations.
To enable direct comparison of ancestry background in different studies, we developed LASER to estimate individual ancestry by placing either sezquenced or genotyped samples in a common ancestry ...space, regardless of the sequencing strategy or genotyping array used to characterize each sample. Here we describe the LASER server to facilitate application of the method to a wide range of genetic studies. The server provides genetic ancestry estimation for different geographic regions and user-friendly interactive visualization of the results.
The LASER server is freely accessible at http://laser.sph.umich.edu/.
dtaliun@umich.edu or wangcl@gis.a-star.edu.sg.
Supplementary data are available at Bioinformatics online.
Understanding the genetic basis of human diseases and traits is dependent on the identification and accurate genotyping of genetic variants. Deep whole-genome sequencing (WGS), the gold standard ...technology for SNP and indel identification and genotyping, remains very expensive for most large studies. Here, we quantify the extent to which array genotyping followed by genotype imputation can approximate WGS in studies of individuals of African, Hispanic/Latino, and European ancestry in the US and of Finnish ancestry in Finland (a population isolate). For each study, we performed genotype imputation by using the genetic variants present on the Illumina Core, OmniExpress, MEGA, and Omni 2.5M arrays with the 1000G, HRC, and TOPMed imputation reference panels. Using the Omni 2.5M array and the TOPMed panel, ≥90% of bi-allelic single-nucleotide variants (SNVs) are well imputed (r2 > 0.8) down to minor-allele frequencies (MAFs) of 0.14% in African, 0.11% in Hispanic/Latino, 0.35% in European, and 0.85% in Finnish ancestries. There was little difference in TOPMed-based imputation quality among the arrays with >700k variants. Individual-level imputation quality varied widely between and within the three US studies. Imputation quality also varied across genomic regions, producing regions where even common (MAF > 5%) variants were consistently not well imputed across ancestries. The extent to which array genotyping and imputation can approximate WGS therefore depends on reference panel, genotype array, sample ancestry, and genomic location. Imputation quality by variant or genomic region can be queried with our new tool, RsqBrowser, now deployed on the Michigan Imputation Server.
This manuscript presents—across multiple ancestries, genotype arrays, and imputation reference panels—the minor-allele frequency thresholds above which array genotyping and imputation accurately capture genetic variants genotyped with deep whole-genome sequencing.
HDL-mediated cholesterol efflux capacity (CEC) may protect from cardiovascular disease. Thus, we aimed to identify its genetic and non-genetic determinants.
We measured CEC to 2% apolipoprotein ...B-depleted serum using BODIPY-cholesterol and cAMP-stimulated J774A.1 macrophages using serum samples from 4,981 participants in the German Chronic Kidney Disease (GCKD) study. Variance of CEC explained by clinical and biochemical parameters in a multivariable linear regression model was calculated by proportional marginal variance decomposition. A genome-wide association study with 7,746,917 variants was performed based on an additive genetic model. The main model was adjusted for age, sex and principal components 1-10. Further models were selected for sensitivity analysis and to reduce residual variance by known CEC pathways.
Variables that explained 1% and more of the variance of CEC were concentrations of triglycerides (12.9%), HDL-cholesterol (11.8%), LDL-cholesterol (3.0%), apolipoprotein A-IV (2.8%), PCSK9 (1.0%), and eGFR (1.0%). The KLKB1 (chr4) and APOE/C1 (chr19) loci were genome-wide significantly (p < 5x10−8) associated with CEC in our main model (p = 8.8x10−10 and p = 3.3x10−10, respectively). KLKB1 remained significantly associated after additional adjustment for either kidney parameters, HDL-cholesterol, triglycerides or apolipoprotein A-IV concentrations, while the APOE/C1 locus was not significantly associated anymore after adjustment for triglycerides. Adjustment for triglycerides also revealed an association with the CLSTN2 locus (chr3; p = 6.0x10−9).
We identified HDL-cholesterol and triglycerides as the main determinants of CEC. Furthermore, we newly found a significant association of CEC with the KLKB1 and the CLSTN2 locus and confirmed the association with the APOE/C1 locus, likely mediated by triglycerides.
Display omitted
•HDL-C and triglycerides are the main biochemical determinants of cholesterol efflux capacity each explaining >11% of variance.•In a GWAS, 2 novel genetic loci (KLKB1 &CLSTN2) and the known APOE/C1 locus are associated with cholesterol efflux capacity.•The KLKB1 locus is significant irrespective of adjustment for kidney function, HDL-C, triglycerides and apolipoprotein A-IV.•The association with the CLSTN2 locus is suggestive and reaches genome-wide significance with adjustment for triglycerides.
Lp(a) concentrations represent a major cardiovascular risk factor and are almost entirely controlled by one single locus (LPA). However, many genetic factors in LPA governing the enormous variance of ...Lp(a) levels are still unknown. Since up to 70% of the LPA coding sequence are located in a difficult to access hypervariable copy number variation named KIV-2, we hypothesized that it may contain novel functional variants with pronounced effects on Lp(a) concentrations. We performed a large scale mutation analysis in the KIV-2 using an extreme phenotype approach.
We compiled an discovery set of 123 samples showing discordance between LPA isoform phenotype and Lp(a) concentrations and controls. Using ultra-deep sequencing, we identified a splice site variant (G4925A) in preferential association with the smaller LPA isoforms. Follow-up in a European general population (n = 2892) revealed an exceptionally high carrier frequency of 22.1% in the general population. The variant explains 20.6% of the Lp(a) variance in carriers of low molecular weight (LMW) apo(a) isoforms (P = 5.75e-38) and reduces Lp(a) concentrations by 31.3 mg/dL. Accordingly the odds ratio for cardiovascular disease was reduced from 1.39 95% confidence interval (CI): 1.17-1.66, P = 1.89e-04 for wildtype LMW individuals to 1.19 95%CI: 0.92; 1.56, P = 0.19 in LMW individuals who were additionally positive for G4925A. Functional studies point towards a reduction of splicing efficiency by this novel variant.
A highly frequent but until now undetected variant in the LPA KIV-2 region is strongly associated with reduced Lp(a) concentrations and reduced cardiovascular risk in LMW individuals.
High-throughput technologies, such as next-generation sequencing, have turned molecular biology into a data-intensive discipline, requiring bioinformaticians to use high-performance computing ...resources and carry out data management and analysis tasks on large scale. Workflow systems can be useful to simplify construction of analysis pipelines that automate tasks, support reproducibility and provide measures for fault-tolerance. However, workflow systems can incur significant development and administration overhead so bioinformatics pipelines are often still built without them. We present the experiences with workflows and workflow systems within the bioinformatics community participating in a series of hackathons and workshops of the EU COST action SeqAhead. The organizations are working on similar problems, but we have addressed them with different strategies and solutions. This fragmentation of efforts is inefficient and leads to redundant and incompatible solutions. Based on our experiences we define a set of recommendations for future systems to enable efficient yet simple bioinformatics workflow construction and execution.
Background
Chronic kidney disease (CKD) represents a chronic proinflammatory state and is associated with very high cardiovascular risk. Apolipoprotein A‐IV (apoA‐IV) has antiatherogenic, ...antioxidative, anti‐inflammatory and antithrombotic properties and levels increase significantly during the course of CKD.
Objectives
We aimed to investigate the association between apoA‐IV and all‐cause mortality and cardiovascular outcomes in the German Chronic Kidney Disease study.
Methods
This was a prospective cohort study including 5141 Caucasian patients with available apoA‐IV measurements and CKD. The majority of the patients had an estimated glomerular filtration rate (eGFR) of 30–60 ml/min/1.73m2 or an eGFR >60 ml/min/1.73m2 in the presence of overt proteinuria. Median follow‐up was 6.5 years. The association of apoA‐IV with comorbidities at baseline and endpoints during follow‐up was modelled adjusting for major confounders.
Results
Mean apoA‐IV concentrations of the entire cohort were 28.9 ± 9.8 mg/dl. Patients in the highest apoA‐IV quartile had the lowest high‐sensitivity C‐reactive protein values despite the highest prevalence of diabetes, albuminuria and the lowest eGFR. Each 10 mg/dl higher apoA‐IV translated into lower odds of prevalent cardiovascular disease (1289 cases, odds ratio = 0.80, 95% confidence interval CI 0.72–0.86, p = 0.0000003). During follow‐up, each 10 mg/dl higher apoA‐IV was significantly associated with a lower risk for all‐cause mortality (600 cases, hazard ratio HR = 0.81, 95% CI 0.73–0.89, p = 0.00004), incident major adverse cardiovascular events (506 cases, HR = 0.88, 95% CI 0.79–0.99, p = 0.03) and death or hospitalizations due to heart failure (346 cases, HR = 0.84, 95% CI 0.73–0.96, p = 0.01).
Conclusions
These data support a link between elevated apoA‐IV concentrations and reduced inflammation in moderate CKD. ApoA‐IV appears to be an independent risk marker for reduced all‐cause mortality, cardiovascular events and heart failure in a large cohort of patients with CKD.
Abstract
Background
There is little information in haemodialysis (HD) patients on whether temporal changes in serum calcium, phosphate or intact parathyroid hormone (iPTH) are associated with ...mortality.
Methods
We analysed associations of phosphate, total calcium and iPTH with all-cause and cardiovascular mortality in 8817 incident HD patients from the European second Analyzing Data, Recognizing Excellence and Optimizing Outcomes (AROii) cohort enrolled in 2007–09, which were prospectively followed for a median of 3 years, using time-dependent Cox proportional hazards models. We evaluated changes in risk over time depending on changes in phosphate, calcium or iPTH.
Results
The association of phosphate and iPTH with all-cause mortality was U-shaped, with the lowest risk ranges between 1.20 and 1.89 mmol/L for phosphate and between 239 and 710 ng/L for iPTH. For total calcium, the associations were J-shaped, with an increased risk for all-cause mortality at levels >2.36 mmol/L. Lowest risk ranges for cardiovascular mortality did not change markedly for all three parameters. If iPTH was below the lowest risk range at baseline (iPTH <239 ng/L), a subsequent increase in levels was associated with improved survival. For phosphate, an increase or decrease out of the lowest risk range was associated with increased mortality risk. For calcium, this was only the case when the values increased above the lowest risk range.
Conclusion
In the AROii cohort, the ranges of bone mineral biomarkers associated with the lowest mortality ranges were largely consistent with the current Kidney Disease: Improving Global Outcomes chronic kidney disease–mineral and bone disorder guideline recommendations. Allowing a suppressed iPTH to increase was associated with a lower mortality, whereas shifts of phosphate or calcium outside the lowest risk range increased mortality.