ForestQC: Quality control on genetic variants from next-generation sequencing data using random forest.

Next-generation sequencing technology (NGS) enables the discovery of nearly all genetic variants present in a genome. A subset of these variants, however, may have poor sequencing quality due to limitations in NGS or variant callers. In genetic studies that analyze a large number of sequenced individuals, it is critical to detect and remove those variants with poor quality as they may cause spurious findings. In this paper, we present ForestQC, a statistical tool for performing quality control on variants identified from NGS data by combining a traditional filtering approach and a machine learning approach. Our software uses the information on sequencing quality, such as sequencing depth, genotyping quality, and GC contents, to predict whether a particular variant is likely to be false-positive. To evaluate ForestQC, we applied it to two whole-genome sequencing datasets where one dataset consists of related individuals from families while the other consists of unrelated individuals. Results indicate that ForestQC outperforms widely used methods for performing quality control on variants such as VQSR of GATK by considerably improving the quality of variants to be included in the analysis. ForestQC is also very efficient, and hence can be applied to large sequencing datasets. We conclude that combining a machine learning algorithm trained with sequencing quality information and the filtering approach is a practical approach to perform quality control on genetic variants from sequencing data

Structural analysis and corrosion studies on an ISO 5832-9 biomedical alloy with TiO2 sol–gel layers

The aim of this study was to demonstrate the relationship between the structural and corrosion properties of an ISO 5832-9 biomedical alloy modified with titanium dioxide (TiO2) layers. These layers were obtained via the sol–gel method by acid-catalyzed hydrolysis of titanium isopropoxide in isopropanol solution. To obtain TiO2 layers with different structural properties, the coated samples were annealed at temperatures of 200, 300, 400, 450, 500, 600 and 800 C for 2 h. For all the prepared samples, accelerated corrosion measurements were performed in Tyrode’s physiological solution using electrochemical methods. The most important corrosion parameters were determined: corrosion potential, polarization resistance, corrosion rate, breakdown and repassivation potentials. Corrosion damage was analyzed using scanning electron microscopy. Structural analysis was carried out for selected TiO2 coatings annealed at 200, 400, 600 and 800 C. In addition, the morphology, chemical composition, crystallinity, thickness and density of the deposited TiO2 layers were determined using suitable electron and X-ray measurement methods. It was shown that the structure and character of interactions between substrate and deposited TiO2 layers depended on annealing temperature. All the obtained TiO2 coatings exhibit anticorrosion properties, but these properties are related to the crystalline structure and character of substrate–layer interaction. From the point of view of corrosion, the best TiO2 sol–gel coatings for stainless steel intended for biomedical applications seem to be those obtained at 400 C.This study was supported by Grant No. N N507 501339 of the National Science Centre. The authors wish to express their thanks to J. Borowski (MEDGAL, Poland) for the Rex 734 alloy

Neighborhood Homogeneous Labelings of Graphs

Given a labeling of the vertices and edges of a graph, we define a type of homogeneity that requires that the neighborhood of every vertex contains the same number of each of the labels. This homogeneity constraint is a generalization of regularity – all such graphs are regular. We consider a specific condition in which both the edge and vertex label sets have two elements and every neighborhood contains two of each label. We show that vertex homogeneity implies edge homogeneity (so long as the number of edges in any neighborhood is four), and give two theorems describing how to build new homogeneous graphs (or multigraphs) from others. Keywords: vertex labeling; edge labeling; homogenous graph; regular graph 1

Contribution of common and rare variants to bipolar disorder susceptibility in extended pedigrees from population isolates.

Current evidence from case/control studies indicates that genetic risk for psychiatric disorders derives primarily from numerous common variants, each with a small phenotypic impact. The literature describing apparent segregation of bipolar disorder (BP) in numerous multigenerational pedigrees suggests that, in such families, large-effect inherited variants might play a greater role. To identify roles of rare and common variants on BP, we conducted genetic analyses in 26 Colombia and Costa Rica pedigrees ascertained for bipolar disorder 1 (BP1), the most severe and heritable form of BP. In these pedigrees, we performed microarray SNP genotyping of 838 individuals and high-coverage whole-genome sequencing of 449 individuals. We compared polygenic risk scores (PRS), estimated using the latest BP1 genome-wide association study (GWAS) summary statistics, between BP1 individuals and related controls. We also evaluated whether BP1 individuals had a higher burden of rare deleterious single-nucleotide variants (SNVs) and rare copy number variants (CNVs) in a set of genes related to BP1. We found that compared with unaffected relatives, BP1 individuals had higher PRS estimated from BP1 GWAS statistics (P = 0.001 ~ 0.007) and displayed modest increase in burdens of rare deleterious SNVs (P = 0.047) and rare CNVs (P = 0.002 ~ 0.033) in genes related to BP1. We did not observe rare variants segregating in the pedigrees. These results suggest that small-to-moderate effect rare and common variants are more likely to contribute to BP1 risk in these extended pedigrees than a few large-effect rare variants

Dlk1 in normal and abnormal hematopoiesis

Udgivelsesdato: AugDlk1 (Pref-1) is a transmembrane and secreted protein, which is a member of the epidermal growth factor-like family, homologous to Notch/Delta/Serrate. We have found by real-time RT-PCR that Dlk1 mRNA levels were high in CD34(+) cells in 10 of 12 MDS samples compared with CD34(+) cells from 11 normals. Also, Dlk1 mRNA was elevated in mononuclear, low density bone marrow cells from 11/38 MDS patients, 5/11 AML M6 and 2/4 AML M7 samples. Furthermore, 5/6 erythroleukemia and 2/2 megakaryocytic leukemia cell lines highly expressed Dlk1 mRNA. Levels of Dlk1 mRNA markedly increased during megakaryocytic differentiation of both CMK megakaryoblasts as well as normal CD34(+) hematopoietic stem cells. High serum levels of Dlk1 occurred in RA (4/10) and essential thrombocythemia (2/10) patients. Functional studies showed that forced expression of Dlk1 enhanced proliferation of K562 cells growing in 1% fetal bovine serum. Analysis of hematopoiesis of Dlk1 knockout mice suggested that Dlk1 contributed to granulocyte, megakaryocyte and B-cell clonogenic growth and was needed for generation of splenic B-cells. In summary, Dlk1 is overexpressed in selected samples of MDS (especially RA and RAEB) and AML (particularly M6, M7), and it appears to be associated with normal development of megakaryocytes and B cells

The biocompatibility of titanium in a buffer solution: compared effects of a thin film of TiO2 deposited by MOCVD and of collagen deposited from a gel

This study aims at evaluating the biocompatibility of titanium surfaces modified according two different ways: (i) deposition of a bio-inert, thin film of rutile TiO2 by chemical vapour deposition (MOCVD), and (ii) biochemical treatment with collagen gel, in order to obtain a bio-interactive coating. Behind the comparison is the idea that either the bio-inert or the bio-active coating has specific advantages when applied to implant treatment, such as the low price of the collagen treatment for instance. The stability in buffer solution was evaluated by open circuit potential (OCP) for medium time and cyclic voltametry. The OCP stabilized after 5104 min for all the specimens except the collagen treated sample which presented a stable OCP from the first minutes. MOCVD treated samples stabilized to more electropositive values. Numeric results were statistically analysed to obtain the regression equations for long time predictable evolution. The corrosion parameters determined from cyclic curves revealed that the MOCVD treatment is an efficient way to improve corrosion resistance. Human dermal fibroblasts were selected for cell culture tests, taking into account that these cells are present in all bio-interfaces, being the main cellular type of connective tissue. The cells grew on either type of surface without phenotype modification. From the reduction of yellow, water-soluble 3-(4,5-dimethyldiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT cytotoxicity test), MOCVD treated samples offer better viability than mechanically polished Ti and collagen treated samples as well. Cell spreading, as evaluated from microscope images processed by the program Sigma Scan, showed also enhancement upon surface modification. Depending on the experimental conditions, MOCVD deposited TiO2 exhibits different nanostructures that may influence biological behaviour. The results demonstrate the capacity of integration in simulated physiologic liquids for an implant pretreated by either method

Rare Copy Number Variants in \u3cem\u3eNRXN1\u3c/em\u3e and \u3cem\u3eCNTN6\u3c/em\u3e Increase Risk for Tourette Syndrome

Tourette syndrome (TS) is a model neuropsychiatric disorder thought to arise from abnormal development and/or maintenance of cortico-striato-thalamo-cortical circuits. TS is highly heritable, but its underlying genetic causes are still elusive, and no genome-wide significant loci have been discovered to date. We analyzed a European ancestry sample of 2,434 TS cases and 4,093 ancestry-matched controls for rare (\u3c 1% frequency) copy-number variants (CNVs) using SNP microarray data. We observed an enrichment of global CNV burden that was prominent for large (\u3e 1 Mb), singleton events (OR = 2.28, 95% CI [1.39–3.79], p = 1.2 × 10−3) and known, pathogenic CNVs (OR = 3.03 [1.85–5.07], p = 1.5 × 10−5). We also identified two individual, genome-wide significant loci, each conferring a substantial increase in TS risk (NRXN1 deletions, OR = 20.3, 95% CI [2.6–156.2]; CNTN6 duplications, OR = 10.1, 95% CI [2.3–45.4]). Approximately 1% of TS cases carry one of these CNVs, indicating that rare structural variation contributes significantly to the genetic architecture of TS

Policy Evaluation with Nonlinear Trended Outcomes: COVID-19 Vaccination Rates in the US

This paper points out some pitfalls in the use of two-way fixed effects (TWFE) regressions when outcome variables contain nonlinear or stochastic trend components. When a policy change shifts trend paths of outcome variables conventional TWFE estimation can distort results and invalidate inference. A robust solution is proposed by identifying determinants of dynamic club membership based on the idea of relative convergence, which can be assessed empirically by the so-called ‘logt’ test (Phillips & Sul, 2007a). Club membership in each time period is estimated by recursive regression, transforming outcome variables to statistically stable, stationary status. Time varying club membership can then be used to identify the determinants of club memberships by running a panel logit or ordered logit regression. This approach is applied to study COVID-19 vaccination data across 50 states and the District of Columbia (DC). A new weekly database is created to track individual state and DC vaccination policies and mandates over the period from March 2021 to February 2022. Initially two convergent clubs are identified. Later evidence of the vaccination rates across states reveals a single convergent club. The primary determinant of this merger of sub-clubs is found to be federal-level vaccine mandates

A note on comonotonicity and positivity of the control components of decoupled quadratic FBSDE

In this small note we are concerned with the solution of Forward-Backward Stochastic Differential Equations (FBSDE) with drivers that grow quadratically in the control component (quadratic growth FBSDE or qgFBSDE). The main theorem is a comparison result that allows comparing componentwise the signs of the control processes of two different qgFBSDE. As a byproduct one obtains conditions that allow establishing the positivity of the control process.Comment: accepted for publicatio