Towards pan-genome read alignment to improve variation calling

Background

Typical human genome differs from the reference genome at 4-5 million sites. This diversity is increasingly catalogued in repositories such as ExAC/gnomAD, consisting of >15,000 whole-genomes and >126,000 exome sequences from different individuals. Despite this enormous diversity, resequencing data workflows are still based on a single human reference genome. Identification and genotyping of genetic variants is typically carried out on short-read data aligned to a single reference, disregarding the underlying variation.

Results

We propose a new unified framework for variant calling with short-read data utilizing a representation of human genetic variation – a pan-genomic reference. We provide a modular pipeline that can be seamlessly incorporated into existing sequencing data analysis workflows. Our tool is open source and available online: https://gitlab.com/dvalenzu/PanVC.

Conclusions

Our experiments show that by replacing a standard human reference with a pan-genomic one we achieve an improvement in single-nucleotide variant calling accuracy and in short indel calling accuracy over the widely adopted Genome Analysis Toolkit (GATK) in difficult genomic regions.

     

Estimation of Canine Leishmania Infection Prevalence in Six Cities of the Algerian Littoral Zone Using a Bayesian Approach

A large-scale study on canine Leishmania infection (CanL) was conducted in six localities along a west-east transect in the Algerian littoral zone (Tlemcen, Mostaganem, Tipaza, Boumerdes, Bejaia, Jijel) and covering two sampling periods. In total 2,184 dogs were tested with an indirect fluorescent antibody test (IFAT) and a direct agglutination test (DAT). Combined multiple-testing and several statistical methods were compared to estimate the CanL true prevalence and tests characteristics (sensitivity and specificity). The Bayesian full model showed the best fit and yielded prevalence estimates between 11% (Mostaganem, first period) and 38% (Bejaia, second period). Sensitivity of IFAT varied (in function of locality) between 86% and 88% while its specificity varied between 65% and 87%. DAT was less sensitive than IFAT but showed a higher specificity (between 80% and 95% in function of locality or/and season). A general increasing trend of the CanL prevalence was noted from west to east. A concordance between the present results and the incidence of human cases of visceral leishmaniasis was observed, where also a maximum was recorded for Bejaia. The results of the present study highlight the dangers when using IFAT as a gold standard.     

Volatile secondary metabolites of Micromeria dalmatica Benth. (Lamiaceae): biosynthetical and chemotaxonomical aspects

Analysis by GC and GC/MS of the essential oil obtained from above-ground parts of Micromeria dalmatica Benth. allowed the identification of 116 components, comprising 93.6% of the total oil composition. The major compounds are 3-oxygenated p-menthane monoterpenes and were identified as pulegone (29.6%), menthone (11.7%), and piperitenone (10.8%). The chemical composition of this and additional 30 oils obtained from selected Micromeria Benth. taxa were compared by using multivariate statistical analysis (agglomerative hierarchical cluster analysis and principal component analysis (PCA)). The results of statistical analyses, as well as the domination of different concurrent p-menthane-skeleton-type monoterpene biosynthetical sub-branches in the compared M. dalmatica samples, implied the occurrence of at least two different chemotypes of the mentioned species.     

Chemical Composition of Hypericum rumeliacum Boiss. Essential Oil. A New Chemotype of This Pharmacologically Valuable Species?

Analysis by GC and GC/MS of the essential‐oil samples obtained from dry above‐ground parts of Hypericum rumeliacum Boiss . (collected in the flowering and fruit‐forming vegetative stages) allowed the identification of 212 components in total, comprising ≥97.8% of the total oil composition. In the flowering phase, the major identified volatile compounds were undecane (6.6%), dodecanal (10.8%), and germacrene D (14.1%), whereas α‐pinene (7.3%), β‐pinene (26.1%), (Z)‐β‐ocimene (8.5%), (E)‐β‐ocimene (10.2%), bicyclogermacrene (7.7%), and germacrene D (15.1%) were dominant in the fruit‐forming phase. Some of the minor constituents found in the studied oil samples (e.g., a homologous series of four 6‐alkyl‐5,6‐dihydro‐2H‐pyran‐2‐ones, i.e., massoia dodeca‐, trideca‐, tetradeca‐, and hexadecalactones) have a restricted occurrence in the Plant Kingdom, and their presence in Hypericum L. spp. has not been previously reported. The chemical compositions of the herein studied additional 34 oils obtained from selected Hypericum taxa were compared using multivariate statistical analysis (agglomerative hierarchical cluster analysis and principal component analysis). The results of these statistical analyses could not be used to either confirm or discard the existence of different H. rumeliacum chemotypes. However, they have implied that the volatile profile of this plant species is determined by the stage of its phenological development.     

The future in and of criticality assessments

In recent literature, the concept of criticality aspires to provide a multifaceted risk assessment of resource supply shortage. However, most existing methodologies for the criticality assessment of raw materials are restricted to a fixed temporal and spatial reference system. They provide a snapshot in time of the equilibrium between supply and demand/economic importance and do not account for temporal changes of their indicators. The static character of criticality assessments limits the use of criticality methodologies to short‐term policy making of raw materials. In the current paper, we argue for an enhancement of the criticality framework to account for three key dynamic characteristics, namely changes of social, technical, and economic features; consideration of the spatial dimension in site‐specific assessments; and impact of changing governance frameworks. We illustrate how these issues were addressed in studies outside of the field of criticality and identify the dynamic parameters that influence resource supply and demand based on a review of studies that belong to the general field of resource supply and demand. The parameters are grouped in seven categories: extraction, social, economic, technical, policy, market dynamics, and environmental. We explore how these parameters were considered in the reviewed studies and propose ways and specific examples of addressing the dynamic effects in the criticality indicators. Furthermore, we discuss the current work on future scenarios to provide reference points for indicator benchmarks. The insights and guidelines derived from the review and our recommendations for future research set the foundations for an enhanced dynamic and site‐specific criticality assessment framework.     

The high road and the low road: trafficking choices in plants

Polar transport of the signaling molecule auxin is critical for plant development and depends on both the polar distribution of auxin efflux carriers, which pump auxin out of the cell and the alignment of these polarized cells. Two papers in this issue of Cell (Michniewicz et al., 2007; Jaillais et al., 2007) address how polar transport of these carriers occurs and describe the endosomal pathways involved.     

Changes in antioxidant status of heart muscle tissue in experimental diabetes in rabbits

The present study was designed to evaluate the oxidative stress-related parameters in alloxan-induced diabetes in rabbits. After 3, 6, 12 and 24 weeks of hyperglycaemia the enzymatic and non-enzymatic factors were measured in heart tissue of diabetic and control groups. Superoxide dismutase and glutathione peroxidase activities and the contents of total sulfhydryl compounds significantly increased at all time intervals. Catalase activity increased initially (after 3 and 6 weeks), decreased after 12 weeks and increased again at the 24th week of the experiment. Glutathione reductase activity increased initially (at 3rd week), decreased below control level after 6 and 12 weeks, then increased again. Ascorbic acid concentration decreased after 3 and 6 weeks, and increased at the 12th and 24th weeks. The level of lipid peroxidation products was reduced after 3, 6 and 12 weeks of the experiment. After 24 weeks it was significantly elevated. These data suggest that hyperglycaemia induces oxidative stress in the heart but the defense mechanisms in the heart tissue are fairly efficacious against oxidative injury.     

GREEN-DB: a framework for the annotation and prioritization of non-coding regulatory variants from whole-genome sequencing data

Non-coding variants have long been recognized as important contributors to common disease risks, but with the expansion of clinical whole genome sequencing, examples of rare, high-impact non-coding variants are also accumulating. Despite recent advances in the study of regulatory elements and the availability of specialized data collections, the systematic annotation of non-coding variants from genome sequencing remains challenging. Here, we propose a new framework for the prioritization of non-coding regulatory variants that integrates information about regulatory regions with prediction scores and HPO-based prioritization. Firstly, we created a comprehensive collection of annotations for regulatory regions including a database of 2.4 million regulatory elements (GREEN-DB) annotated with controlled gene(s), tissue(s) and associated phenotype(s) where available. Secondly, we calculated a variation constraint metric and showed that constrained regulatory regions associate with disease-associated genes and essential genes from mouse knock-outs. Thirdly, we compared 19 non-coding impact prediction scores providing suggestions for variant prioritization. Finally, we developed a VCF annotation tool (GREEN-VARAN) that can integrate all these elements to annotate variants for their potential regulatory impact. In our evaluation, we show that GREEN-DB can capture previously published disease-associated non-coding variants as well as identify additional candidate disease genes in trio analyses.     

uPA-silica-Particles (SP-uPA): a novel analytical system to investigate uPA-uPAR interaction and to test synthetic uPAR antagonists as potential cancer therapeutics

The urokinase-type plasminogen activation system, including the serine protease uPA (urokinase-type plasminogen activator) and its cell surface receptor (uPAR, CD87), are important key molecules in tumor invasion and metastasis. Besides its proteolytic function, binding of uPA to uPAR on tumor cells exerts various cell responses such as migration, adhesion, proliferation, and differentiation. Hence, the uPA/uPAR system is a potential target for tumor therapy. We have designed a new generation of uPA-derived synthetic cyclic peptides suited to interfere with the binding of uPA to uPAR and present a new technology involving micro silica particles coated with uPA (SP-uPA) and reacting with recombinant soluble uPAR (suPAR), to rapidly assess the antagonistic potential of uPA-peptides by flow cytofluorometry (FACS). For this, we used silica particles of 10 microm in diameter to which HMW-uPA is coupled using the EDC/NHS method. Soluble, recombinant suPAR was added and the interaction of SP-uPA with suPAR verified by reaction with monoclonal antibody HD13.1 directed to uPAR, followed by a cyan dye (cy5)-labeled antibody directed against mouse IgG. Thereby it was possible to test naturally occurring ligands of uPAR (HMW-uPA, ATF) as well as highly effective, synthetic cyclic uPA-derived peptides (cyclo21,29[D-Cys21Cys29]-UPA21-30, cyclo21,29[D-Cys21Nle28Cys29]-uPA21-30, cyclo21,29[D-Cys(21)2-Nal24Cys29]-uPA21-30, and cyclo21,29[D-Cys21Orn23Thi24Thi25Cys29]-uPA21-30. The results obtained with the noncellular SP-uPA/uPAR system are highly comparable to those obtained with a cellular system involving FITC-uPA and the promyeloid cell line U937 as the source of uPAR.