Bayesian Gaussian Copula Factor Models for Mixed Data.

Gaussian factor models have proven widely useful for parsimoniously characterizing dependence in multivariate data. There is a rich literature on their extension to mixed categorical and continuous variables, using latent Gaussian variables or through generalized latent trait models acommodating measurements in the exponential family. However, when generalizing to non-Gaussian measured variables the latent variables typically influence both the dependence structure and the form of the marginal distributions, complicating interpretation and introducing artifacts. To address this problem we propose a novel class of Bayesian Gaussian copula factor models which decouple the latent factors from the marginal distributions. A semiparametric specification for the marginals based on the extended rank likelihood yields straightforward implementation and substantial computational gains. We provide new theoretical and empirical justifications for using this likelihood in Bayesian inference. We propose new default priors for the factor loadings and develop efficient parameter-expanded Gibbs sampling for posterior computation. The methods are evaluated through simulations and applied to a dataset in political science. The models in this paper are implemented in the R package bfa.

Recovery from an acute infection in C. elegans requires the GATA transcription factor ELT-2.

The mechanisms involved in the recognition of microbial pathogens and activation of the immune system have been extensively studied. However, the mechanisms involved in the recovery phase of an infection are incompletely characterized at both the cellular and physiological levels. Here, we establish a Caenorhabditis elegans-Salmonella enterica model of acute infection and antibiotic treatment for studying biological changes during the resolution phase of an infection. Using whole genome expression profiles of acutely infected animals, we found that genes that are markers of innate immunity are down-regulated upon recovery, while genes involved in xenobiotic detoxification, redox regulation, and cellular homeostasis are up-regulated. In silico analyses demonstrated that genes altered during recovery from infection were transcriptionally regulated by conserved transcription factors, including GATA/ELT-2, FOXO/DAF-16, and Nrf/SKN-1. Finally, we found that recovery from an acute bacterial infection is dependent on ELT-2 activity.

B cells in rheumatoid synovitis.

In rheumatoid arthritis, T cells, B cells, macrophages, and dendritic cells invade the synovial membranes, establishing complex microstructures that promote inflammatory/tissue destructive lesions. B cell involvement has been considered to be limited to autoantibody production. However, recent studies suggest that B cells support rheumatoid disease through other mechanisms. A critical element of rheumatoid synovitis is the process of ectopic lymphoid neogenesis, with highly efficient lymphoid architectures established in a nonlymphoid tissue site. Rheumatoid synovitis recapitulates the pathways of lymph node formation, and B cells play a key role in this process. Furthermore, studies of rheumatoid lesions implanted in immunodeficient mice suggest that T cell activation in synovitis is B cell dependent, indicating the role played by B cells in presenting antigens and providing survival signals.

Polymorphisms in the kinesin-like factor 1 B gene and risk of epithelial ovarian cancer in Eastern Chinese women.

The kinesin-like factor 1 B (KIF1B) gene plays an important role in the process of apoptosis and the transformation and progression of malignant cells. Genetic variations in KIF1B may contribute to risk of epithelial ovarian cancer (EOC). In this study of 1,324 EOC patients and 1,386 cancer-free female controls, we investigated associations between two potentially functional single nucleotide polymorphisms in KIF1B and EOC risk by the conditional logistic regression analysis. General linear regression model was used to evaluate the correlation between the number of variant alleles and KIF1B mRNA expression levels. We found that the rs17401966 variant AG/GG genotypes were significantly associated with a decreased risk of EOC (adjusted odds ratio (OR) = 0.81, 95 % confidence interval (CI) = 0.68-0.97), compared with the AA genotype, but no associations were observed for rs1002076. Women who carried both rs17401966 AG/GG and rs1002076 AG/AA genotypes of KIF1B had a 0.82-fold decreased risk (adjusted 95 % CI = 0.69-0.97), compared with others. Additionally, there was no evidence of possible interactions between about-mentioned co-variants. Further genotype-phenotype correlation analysis indicated that the number of rs17401966 variant G allele was significantly associated with KIF1B mRNA expression levels (P for GLM = 0.003 and 0.001 in all and Chinese subjects, respectively), with GG carriers having the lowest level of KIF1B mRNA expression. Taken together, the rs17401966 polymorphism likely regulates KIF1B mRNA expression and thus may be associated with EOC risk in Eastern Chinese women. Larger, independent studies are warranted to validate our findings.

Explicit DNase sequence bias modeling enables high-resolution transcription factor footprint detection.

DNaseI footprinting is an established assay for identifying transcription factor (TF)-DNA interactions with single base pair resolution. High-throughput DNase-seq assays have recently been used to detect in vivo DNase footprints across the genome. Multiple computational approaches have been developed to identify DNase-seq footprints as predictors of TF binding. However, recent studies have pointed to a substantial cleavage bias of DNase and its negative impact on predictive performance of footprinting. To assess the potential for using DNase-seq to identify individual binding sites, we performed DNase-seq on deproteinized genomic DNA and determined sequence cleavage bias. This allowed us to build bias corrected and TF-specific footprint models. The predictive performance of these models demonstrated that predicted footprints corresponded to high-confidence TF-DNA interactions. DNase-seq footprints were absent under a fraction of ChIP-seq peaks, which we show to be indicative of weaker binding, indirect TF-DNA interactions or possible ChIP artifacts. The modeling approach was also able to detect variation in the consensus motifs that TFs bind to. Finally, cell type specific footprints were detected within DNase hypersensitive sites that are present in multiple cell types, further supporting that footprints can identify changes in TF binding that are not detectable using other strategies.

Drug-drug interaction studies of cardiovascular drugs involving P-glycoprotein, an efflux transporter, on the pharmacokinetics of edoxaban, an oral factor Xa inhibitor.

BACKGROUND: Edoxaban, an oral direct factor Xa inhibitor, is in development for thromboprophylaxis, including prevention of stroke and systemic embolism in patients with atrial fibrillation (AF). P-glycoprotein (P-gp), an efflux transporter, modulates absorption and excretion of xenobiotics. Edoxaban is a P-gp substrate, and several cardiovascular (CV) drugs have the potential to inhibit P-gp and increase drug exposure. OBJECTIVE: To assess the potential pharmacokinetic interactions of edoxaban and 6 cardiovascular drugs used in the management of AF and known P-gp substrates/inhibitors. METHODS: Drug-drug interaction studies with edoxaban and CV drugs with known P-gp substrate/inhibitor potential were conducted in healthy subjects. In 4 crossover, 2-period, 2-treatment studies, subjects received edoxaban 60 mg alone and coadministered with quinidine 300 mg (n = 42), verapamil 240 mg (n = 34), atorvastatin 80 mg (n = 32), or dronedarone 400 mg (n = 34). Additionally, edoxaban 60 mg alone and coadministered with amiodarone 400 mg (n = 30) or digoxin 0.25 mg (n = 48) was evaluated in a single-sequence study and 2-cohort study, respectively. RESULTS: Edoxaban exposure measured as area under the curve increased for concomitant administration of edoxaban with quinidine (76.7 %), verapamil (52.7 %), amiodarone (39.8 %), and dronedarone (84.5 %), and exposure measured as 24-h concentrations for quinidine (11.8 %), verapamil (29.1 %), and dronedarone (157.6 %) also increased. Administration of edoxaban with amiodarone decreased the 24-h concentration for edoxaban by 25.7 %. Concomitant administration with digoxin or atorvastatin had minimal effects on edoxaban exposure. CONCLUSION: Coadministration of the P-gp inhibitors quinidine, verapamil, and dronedarone increased edoxaban exposure. Modest/minimal effects were observed for amiodarone, atorvastatin, and digoxin.

A novel inflammatory biomarker, GlycA, associates with disease activity in rheumatoid arthritis and cardio-metabolic risk in BMI-matched controls

BACKGROUND: RA and CVD both have inflammation as part of the underlying biology. Our objective was to explore the relationships of GlycA, a measure of glycosylated acute phase proteins, with inflammation and cardiometabolic risk in RA, and explore whether these relationships were similar to those for persons without RA. METHODS: Plasma GlycA was determined for 50 individuals with mild-moderate RA disease activity and 39 controls matched for age, gender, and body mass index (BMI). Regression analyses were performed to assess relationships between GlycA and important markers of traditional inflammation and cardio-metabolic health: inflammatory cytokines, disease activity, measures of adiposity and insulin resistance. RESULTS: On average, RA activity was low (DAS-28 = 3.0 ± 1.4). Traditional inflammatory markers, ESR, hsCRP, IL-1β, IL-6, IL-18 and TNF-α were greater in RA versus controls (P < 0.05 for all). GlycA concentrations were significantly elevated in RA versus controls (P = 0.036). In RA, greater GlycA associated with disease activity (DAS-28; RDAS-28 = 0.5) and inflammation (RESR = 0.7, RhsCRP = 0.7, RIL-6 = 0.3: P < 0.05 for all); in BMI-matched controls, these inflammatory associations were absent or weaker (hsCRP), but GlycA was related to IL-18 (RhsCRP = 0.3, RIL-18 = 0.4: P < 0.05). In RA, greater GlycA associated with more total abdominal adiposity and less muscle density (Rabdominal-adiposity = 0.3, Rmuscle-density = -0.3, P < 0.05 for both). In BMI-matched controls, GlycA associated with more cardio-metabolic markers: BMI, waist circumference, adiposity measures and insulin resistance (R = 0.3-0.6, P < 0.05 for all). CONCLUSIONS: GlycA provides an integrated measure of inflammation with contributions from traditional inflammatory markers and cardio-metabolic sources, dominated by inflammatory markers in persons with RA and cardio-metabolic factors in those without.