The Secreted Ly6/uPAR-Related Protein 1 (Slurp1) Modulates Corneal Angiogenic Inflammation Via NF-κB Signaling.

Purpose: Previously we demonstrated that the secreted Ly-6/uPAR related protein 1 (SLURP1), abundantly expressed in the corneal epithelium (CE) and secreted into the tear fluid, serves as an antiangiogenic molecule. Here we describe the Slurp1-null (Slurp1X-/-) mouse corneal response to silver nitrate (AgNO3) cautery. Methods: Five days after AgNO3 cautery, we compared the wild-type (WT) and Slurp1X-/- mouse (1) corneal neovascularization (CNV) and immune cell influx by whole-mount immunofluorescent staining for CD31 and CD45, (2) macrophage and neutrophil infiltration by flow cytometry, and (3) gene expression by quantitative RT-PCR. Quantitative RT-PCR, immunofluorescent staining, and immunoblots were employed to evaluate the expression, phosphorylation status, and subcellular localization of NF-κB pathway components. Results: Unlike the WT, the Slurp1X-/- corneas displayed denser CNV in response to AgNO3 cautery, with more infiltrating macrophages and neutrophils and greater upregulation of the transcripts encoding VEGFA, MMP2, IL-1b, and vimentin. At 2, 7, and 10 days after AgNO3 cautery, Slurp1 expression was significantly downregulated in the WT corneas. Compared with the WT, naive Slurp1X-/- CE displayed increased phosphorylation of IKK(a/b), elevated phosphorylation of IκB with decreased amounts of total IκB, and higher phosphorylation of NF-κB, suggesting that NF-κB signaling is constitutively active in naive Slurp1X-/- corneas. Conclusions: Enhanced angiogenic inflammation in AgNO3 cauterized Slurp1X-/- corneas and constitutively active status of NF-κB signaling in the absence of Slurp1 suggest that Slurp1 modulates corneal angiogenic inflammation via NF-κB signaling.

Bayesian Statistics for Medical Devices: Progress Since 2010.

The use of Bayesian statistics to support regulatory evaluation of medical devices began in the late 1990s. We review the literature, focusing on recent developments of Bayesian methods, including hierarchical modeling of studies and subgroups, borrowing strength from prior data, effective sample size, Bayesian adaptive designs, pediatric extrapolation, benefit-risk decision analysis, use of real-world evidence, and diagnostic device evaluation. We illustrate how these developments were utilized in recent medical device evaluations. In Supplementary Material, we provide a list of medical devices for which Bayesian statistics were used to support approval by the US Food and Drug Administration (FDA), including those since 2010, the year the FDA published their guidance on Bayesian statistics for medical devices. We conclude with a discussion of current and future challenges and opportunities for Bayesian statistics, including artificial intelligence/machine learning (AI/ML) Bayesian modeling, uncertainty quantification, Bayesian approaches using propensity scores, and computational challenges for high dimensional data and models.

Comparing a Bayesian Approach (BEST) with the Two One-Sided t-Tests (TOSTs) for Bioequivalence Studies.

The two one-sided t-tests (TOST) procedure has been used to evaluate average bioequivalence (BE). As a regulatory standard, it is crucial that TOST distinguish BE from not-BE (NBE) when BE data are not lognormal. TOST was compared with a Bayesian procedure (BEST by Kruschke) in simulated datasets of test/reference ratios (T/R) which were BE and NBE, wherein (1) log(T/R) or T-R were normally distributed, (2) sample sizes ranged 10-50, and (3) extreme log(T/R) or T-R values were randomly included in datasets. The 90% "credible interval" (CrI) from BEST is a Bayesian alternative of the 90% confidence interval (CI) of TOST and it can be derived from a posterior distribution that is more reflective of the observed mean log(T/R) distribution that often deviates from normality. In the absence of extreme T/R values, both methods agreed BE when observed T/R were lognormal. BEST more accurately concluded BE or NBE, while requiring fewer subjects, when observed log(T/R) or T-R were normal in the presence of extreme values. Overall, TOST and BEST perform comparably on lognormal T/R, while BEST is more accurate, requiring fewer subjects when datasets are normal for T-R or contain extreme values. Of note, the normally distributed datasets only rarely contain extreme values. Our results imply that when BEST and TOST yield different BE assessment results from bioequivalent products, TOST may disadvantage applicants when T/R are not lognormal and/or include extreme T/R values. Application of BEST can address the situation when T/R are not lognormal or include extreme data values. Application of BEST to BE data can be considered a useful alternative to TOST for evaluation of BE and for efficient development of BE formulations.

Ectopic expression of T in the paraxial mesoderm disrupts somite maturation in the mouse.

T is the founding member of the T-box family of transcription factors; family members are critical for cell fate decisions and tissue morphogenesis throughout the animal kingdom. T is expressed in the primitive streak and notochord with mouse mutant studies revealing its critical role in mesoderm formation in the primitive streak and notochord integrity. We previously demonstrated that misexpression of Tbx6 in the paraxial and lateral plate mesoderm results in embryos resembling Tbx15 and Tbx18 nulls. This, together with results from in vitro transcriptional assays, suggested that ectopically expressed Tbx6 can compete with endogenously expressed Tbx15 and Tbx18 at the binding sites of target genes. Since T-box proteins share a similar DNA binding domain, we hypothesized that misexpressing T in the paraxial and lateral plate mesoderm would also interfere with the endogenous Tbx15 and Tbx18, causing embryonic phenotypes resembling those seen upon Tbx6 expression in the somites and limbs. Interestingly, ectopic T expression led to distinct embryonic phenotypes, specifically, reduced-sized somites in embryos expressing the highest levels of T, which ultimately affects axis length and neural tube morphogenesis. We further demonstrate that ectopic T leads to ectopic expression of Tbx6 and Mesogenin 1, known targets of T. These results suggests that ectopic T expression contributes to the phenotype by activating its own targets rather than via a straight competition with endogenous T-box factors.

Secreted Ly-6/uPAR-related protein-1 (SLURP1) is a pro-differentiation factor that stalls G1-S transition during corneal epithelial cell cycle progression.

PURPOSE: Previously we demonstrated that the secreted Ly-6/uPAR related protein-1 (SLURP1), abundantly expressed in the corneal epithelium (CE) and secreted into the tear fluid, serves as an anti-inflammatory and anti-angiogenic molecule. Here we describe the Slurp1-null (Slurp1X-/-) mouse corneal phenotype for the first time. METHODS: We compared the 10-week-old wild type (WT) and Slurp1X-/- mouse corneal (i) histology by hematoxylin-eosin and periodic acid-Schiff's reagent staining, (ii) cell proliferation by immunostaining for Ki67, (iii) cell adhesion molecules by immunostaining for desmosomal and tight junction proteins, (iv) barrier function by fluorescein staining and (v) wound-healing by epithelial debridement. Effect of SLURP1 on cell cycle was quantified in human corneal limbal epithelial (HCLE) cells engineered to express SLURP1 (HCLE-SLURP1). RESULTS: WT and Slurp1X-/- corneal histology was largely comparable, other than a few loosely attached superficial cells in Slurp1X-/- corneas. Compared with the WT, Slurp1X-/- corneas displayed (i) increase in Ki67+ cells, (ii) altered expression and/or localization of tight junction proteins Tjp1 and Pard3, and desmosomal Dsp, (iii) increased superficial fragility and (iv) slower CE wound healing. HCLE-SLURP1 cells displayed (i) decrease in Ki67+ cells, (ii) increased cell number doubling time, (iii) stalling in G1-S phase transition during cell cycle, and (iv) downregulation of cyclins CCNE and CCND1/D2, cyclin-dependent kinases CDK4 and CDK6, and upregulation of CDK inhibitor p15/CDKN2B. CONCLUSIONS: Collectively, these results elucidate that Slurp1X-/- CE cell homeostasis is altered and suggest that SLURP1 is a pro-differentiation factor that stalls G1-S transition during cell cycle progression by downregulating cyclins and upregulating p15/CDKN2B.

BMP6 Regulates Corneal Epithelial Cell Stratification by Coordinating Their Proliferation and Differentiation and Is Upregulated in Pterygium.

Purpose: Proper balance between cell proliferation and differentiation is essential for corneal epithelial (CE) stratification and homeostasis. Although bone morphogenetic protein-6 (BMP6) is known to be expressed in the CE for over 25 years, its function in this tissue remains unknown. Here, we test the hypothesis that BMP6 promotes CE cell stratification and homeostasis by regulating their proliferation and differentiation. Methods: We employed postnatal day-12 (PN-12), PN-14, PN-20, and PN-90 mouse eyes; human corneal limbal epithelial (HCLE) cells; and ocular surface fibrovascular disease pterygium tissues to evaluate the role of BMP6 in CE proliferation, differentiation, and pathology by RT-qPCR, immunoblots, and/or immunofluorescent staining. Cell proliferation was quantified by immunostaining for Ki67. Results: Coincident with the mouse CE stratification between PN-12 and PN-20, BMP6 was significantly upregulated and the BMP6 antagonist Noggin downregulated. Mature CE retained high BMP6 and low Noggin expression at PN-90. BMP6 and its receptors BMPR1A and BMPR2 were upregulated during in vitro stratification of HCLE cells. Consistent with its anti-proliferative role, exogenous BMP6 suppressed HCLE cell proliferation, downregulated cyclin-D1 and cyclin-D2, and upregulated cell-cycle inhibitors Krüppel-like factor 4 (KLF4) and p21. BMP6 also upregulated the desmosomal cadherins desmoplakin and desmoglein in HCLE cells, consistent with its pro-differentiation role. Human pterygium displayed significant upregulation of BMP6 coupled with downregulation of Noggin and cell-cycle suppressors KLF4 and p21. Conclusions: BMP6 coordinates CE stratification and homeostasis by regulating their proliferation and differentiation. BMP6 is significantly upregulated in human pterygium concurrent with downregulation of Noggin, KLF4, and p21.

The secreted Ly-6/uPAR related protein-1 (SLURP1) stabilizes epithelial cell junctions and suppresses TNF-α-induced cytokine production.

The secreted Ly-6/uPAR related protein-1 (SLURP1) is an anti-angiogenic and anti-inflammatory peptide highly expressed by the mucosal epithelial cells. SLURP1 is abundantly expressed by the corneal epithelial cells and is significantly downregulated when these cells are transformed and adapted for culture in vitro. Here we studied the effect of overexpressing SLURP1 in Human Corneal Limbal Epithelial (HCLE) cells cultured in vitro. The expression of DSP1, DSG1, TJP1 and E-Cadherin was significantly upregulated in two different SLURP1-overexpressing HCLE cell (HCLE-SLURP1) clones. HCLE-SLURP1 cells also displayed a significant decrease in tumor necrosis factor-α (TNF-α)-induced upregulation of (i) IL-8 from 7.4- to 2.9- and 2.1-fold, (ii) IL-1ß from 4.9- to 3.9- and 2.9-fold, (iii) CXCL1 from 9- to 3.3- and 5.5-fold, and (iv) CXCL2 from 4.8- to 2.1- and 2.8-fold. ELISAs revealed a concomitant decrease in IL-8 levels in cell culture supernatants from 789 pg/ml in the control, to 503 and 352 pg/ml in HCLE-SLURP1 cells. Consistently, cytosolic IκB expression was elevated in HCLE-SLURP1 cells with a concurrent suppression of TNF-α-activated nuclear translocation of NF-κB. Collectively, these results elucidate the beneficial effects of SLURP1 in stabilizing the HCLE intercellular junctions and suppressing the TNF-α-induced upregulation of inflammatory cytokines by suppressing NF-κB nuclear translocation.

Intraflagellar transporter protein 140 (IFT140), a component of IFT-A complex, is essential for male fertility and spermiogenesis in mice.

Intraflagellar transport (IFT) is a conserved mechanism essential for the assembly and maintenance of most eukaryotic cilia and flagella. However, little is known about its role in sperm flagella formation and male fertility. IFT140 is a component of IFT-A complex. In mouse, it is highly expressed in the testis. Ift140 gene was inactivated specifically in mouse spermatocytes/spermatids. The mutant mice did not show any gross abnormalities, but all were infertile and associated with significantly reduced sperm number and motility. Multiple sperm morphological abnormalities were discovered, including amorphous heads, short/bent flagella and swollen tail tips, as well as vesicles along the flagella due to spermiogenesis defects. The epididymides contained round bodies of cytoplasm derived from the sloughing of the cytoplasmic lobes and residual bodies. Knockout of Ift140 did not significantly affect testicular expression levels of selective IFT components but localization of IFT27 and IFT88, two components of IFT-B complex, was changed. Our findings demonstrate that IFT140 is a key regulator for male fertility and normal spermiogenesis in mice. It not only plays a role in sperm flagella assembling, but is also involved in critical assembly of proteins that interface between the germ cell plasma and the Sertoli cell.

Analysis of in vitro evolution reveals the underlying distribution of catalytic activity among random sequences.

The emergence of catalytic RNA is believed to have been a key event during the origin of life. Understanding how catalytic activity is distributed across random sequences is fundamental to estimating the probability that catalytic sequences would emerge. Here, we analyze the in vitro evolution of triphosphorylating ribozymes and translate their fitnesses into absolute estimates of catalytic activity for hundreds of ribozyme families. The analysis efficiently identified highly active ribozymes and estimated catalytic activity with good accuracy. The evolutionary dynamics follow Fisher's Fundamental Theorem of Natural Selection and a corollary, permitting retrospective inference of the distribution of fitness and activity in the random sequence pool for the first time. The frequency distribution of rate constants appears to be log-normal, with a surprisingly steep dropoff at higher activity, consistent with a mechanism for the emergence of activity as the product of many independent contributions.

Clinical Correlates of Benefit From Radium-223 Therapy in Metastatic Castration Resistant Prostate Cancer.

BACKGROUND: We sought to identify potential clinical variables associated with outcomes after radium-223 therapy in routine practice. METHODS: Consecutive non-trial mCRPC patients who received ≥1 dose of radium dichloride-223 at four academic and one community urology-specific cancer centers from May 2013 to June 2014 were retrospectively identified. Association of baseline and on-therapy clinical variables with number of radium doses received and clinical outcomes including overall survival were analyzed using chi-square statistics, cox proportional hazards, and Kaplan-Meier methods. Bone Scan Index (BSI) was derived from available bone scans using EXINI software. RESULTS: One hundred and forty-five patients were included. Radium-223 was administered for six cycles in 74 patients (51%). One-year survival in this heavily pre-treated population was 64% (95%CI: 54-73%). In univariate and multivariate analysis, survival was highly associated with receiving all six doses of Radium-223. Receipt of six doses was associated with ECOG PS of 0-1, lower baseline PSA & pain level, no prior abiraterone/enzalutamide, <5 BSI value, and normal alkaline phosphatase. In patients who reported baseline pain (n = 72), pain declined in 51% after one dose and increased in 7%. PSA declined ≥50% in 16% (18/110). Alkaline phosphatase declined ≥25% in 48% (33/69) and ≥50% in 16/69 patients. BSI declined in 17 (68%) of the 25 patients who had bone scan available at treatment follow-up. Grade ≥3 neutropenia, anemia, and thrombocytopenia occurred in 4% (n = 114), 4% (n = 125), and 5% (n = 123), respectively. CONCLUSIONS: Patients earlier in their disease course with <5 BSI, low pain score, and good ECOG performance status are optimal candidates for radium-223. Radium-223 therapy is well tolerated with most patients reporting declines in pain scores and BSI. Prostate 77:479-488, 2017. © 2016 Wiley Periodicals, Inc.

Utilizing national and international registries to enhance pre-market medical device regulatory evaluation.

Regulatory decisions are made based on the assessment of risk and benefit of medical devices at the time of pre-market approval and subsequently, when post-market risk-benefit balance needs reevaluation. Such assessments depend on scientific evidence obtained from pre-market studies, post-approval studies, post-market surveillance studies, patient perspective information, as well as other real world data such as national and international registries. Such registries provide real world evidence and are playing a more and more important role in enhancing the safety and effectiveness evaluation of medical devices. While these registries provide large quantities of data reflecting real world practice and can potentially reduce the cost of clinical trials, challenges arise concerning (1) data quality adequate for regulatory decision-making, (2) bias introduced at every stage and aspect of study, (3) scientific validity of study designs, and (4) reliability and interpretability of study results. This article will discuss related statistical and regulatory challenges and opportunities with examples encountered in medical device regulatory reviews.

Computational analysis of fitness landscapes and evolutionary networks from in vitro evolution experiments.

In vitro selection experiments in biochemistry allow for the discovery of novel molecules capable of specific desired biochemical functions. However, this is not the only benefit we can obtain from such selection experiments. Since selection from a random library yields an unprecedented, and sometimes comprehensive, view of how a particular biochemical function is distributed across sequence space, selection experiments also provide data for creating and analyzing molecular fitness landscapes, which directly map function (phenotypes) to sequence information (genotypes). Given the importance of understanding the relationship between sequence and functional activity, reliable methods to build and analyze fitness landscapes are needed. Here, we present some statistical methods to extract this information from pools of RNA molecules. We also provide new computational tools to construct and study molecular fitness landscapes.

Statistical innovations in the medical device world sparked by the FDA.

The world of medical devices while highly diverse is extremely innovative, and this facilitates the adoption of innovative statistical techniques. Statisticians in the Center for Devices and Radiological Health (CDRH) at the Food and Drug Administration (FDA) have provided leadership in implementing statistical innovations. The innovations discussed include: the incorporation of Bayesian methods in clinical trials, adaptive designs, the use and development of propensity score methodology in the design and analysis of non-randomized observational studies, the use of tipping-point analysis for missing data, techniques for diagnostic test evaluation, bridging studies for companion diagnostic tests, quantitative benefit-risk decisions, and patient preference studies.

Tailoring and understanding the mechanical properties of nanoparticle-shelled bubbles.

One common approach to generate lightweight materials with high specific strength and stiffness is the incorporation of stiff hollow microparticles (also known as bubbles or microballoons) into a polymeric matrix. The mechanical properties of these composites, also known as syntactic foams, greatly depend on those of the hollow microparticles. It is critical to precisely control the properties of these bubbles to fabricate lightweight materials that are suitable for specific applications. In this paper, we present a method to tailor the mechanical properties and response of highly monodisperse nanoparticle-shelled bubbles using thermal treatment. We characterize the mechanical properties of individual as-assembled bubbles as well as those of thermally treated ones using nanoindentation and quantitative in situ compression tests. As-assembled bubbles display inelastic response, whereas thermally treated bubbles behave elastically. We also show that the stiffness and strength of bubbles are enhanced significantly, as much as 12 and 14 times that of the as-assembled bubbles, respectively, via thermal treatment. We complement the experimental results with finite element analysis (FEA) to understand the effect of shell thickness nonuniformity as well as the inelasticity on the mechanical response and fracture behavior of these bubbles. We demonstrate that the failure mechanism of bubbles incorporated into a polymer composite depends on the structure of the bubbles.

The paradox of dual roles in the RNA world: resolving the conflict between stable folding and templating ability.

The hypothesized dual roles of RNA as both information carrier and biocatalyst during the earliest stages of life require a combination of features: good templating ability (for replication) and stable folding (for ribozymes). However, this poses the following paradox: well-folded sequences are poor templates for copying, but poorly folded sequences are unlikely to be good ribozymes. Here, we describe a strategy to overcome this dilemma through G:U wobble pairing in RNA. Unlike Watson-Crick base pairs, wobble pairs contribute highly to the energetic stability of the folded structure of their sequence, but only slightly, if at all, to the stability of the folded reverse complement. Sequences in the RNA World might thereby combine stable folding of the ribozyme with an unstructured, reverse-complementary genome, resulting in a "division of labor" between the strands. We demonstrate this strategy using computational simulations of RNA folding and an experimental model of early replication, nonenzymatic template-directed RNA primer extension. Additional study is needed to solve other problems associated with a complete replication cycle, including separation of strands after copying. Interestingly, viroid RNA sequences, which have been suggested to be relics of an RNA World (Diener, Proc Natl Acad Sci USA 86:9370-9374, 1989), also show significant asymmetry in folding energy between the infectious (+) and template (-) strands due to G:U pairing, suggesting that this strategy may even be used by replicators in the present day.

Merogonic stages of Theileria cervi in mule deer (Odocoileus hemionus).

In February 2012, 12 farmed mule deer (Odocoileus hemionus) were moved from a facility in southwestern Oklahoma to a facility in southeastern Oklahoma that housed 100 farmed white-tailed deer (Odocoileus virginianus). Between the third and fifth weeks, 9 of the 12 mule deer had died, 4 of which were submitted for necropsy. The deer were heavily infested with Amblyomma americanum (lone star ticks). Hematologic data from 1 deer revealed severe anemia, leukocytosis, and intraerythrocytic hemoparasites consistent with Theileria spp. Microscopically, the liver, lymph nodes, and spleen contained multifocally distributed, enlarged monocytic cells whose cytoplasm was replaced by developing meronts in various stages of merogony. It appears that, upon arrival, the Theileria cervi-naïve mule deer became infested with large numbers of Theileria-infected lone star ticks leading to massive exposure of the mule deer to sporozoites of the protozoan, resulting in an acute hemolytic crisis and fatalities. The merogonic stages of T. cervi are also described. The lack of earlier reports of merogony may be due to the fact that only a single, short-lived, merogonic cycle follows exposure to sporozoites and thus merogonic stages are demonstrable for only a short period. Polymerase chain reaction testing of paraffin-embedded tissue yielded a 507-bp amplicon sequence that was 100% identical with the sequence of T. cervi previously reported from white-tailed deer in Oklahoma and from elk in Wisconsin and Indiana.

Comprehensive experimental fitness landscape and evolutionary network for small RNA.

The origin of life is believed to have progressed through an RNA world, in which RNA acted as both genetic material and functional molecules. The structure of the evolutionary fitness landscape of RNA would determine natural selection for the first functional sequences. Fitness landscapes are the subject of much speculation, but their structure is essentially unknown. Here we describe a comprehensive map of a fitness landscape, exploring nearly all of sequence space, for short RNAs surviving selection in vitro. With the exception of a small evolutionary network, we find that fitness peaks are largely isolated from one another, highlighting the importance of historical contingency and indicating that natural selection would be constrained to local exploration in the RNA world.

Site specific probability of passive acoustic detection of humpback whale calls from single fixed hydrophones.

Passive acoustic monitoring of marine mammal calls is an increasingly important method for assessing population numbers, distribution, and behavior. A common mistake in the analysis of marine mammal acoustic data is formulating conclusions about these animals without first understanding how environmental properties such as bathymetry, sediment properties, water column sound speed, and ocean acoustic noise influence the detection and character of vocalizations in the acoustic data. The approach in this paper is to use Monte Carlo simulations with a full wave field acoustic propagation model to characterize the site specific probability of detection of six types of humpback whale calls at three passive acoustic monitoring locations off the California coast. Results show that the probability of detection can vary by factors greater than ten when comparing detections across locations, or comparing detections at the same location over time, due to environmental effects. Effects of uncertainties in the inputs to the propagation model are also quantified, and the model accuracy is assessed by comparing calling statistics amassed from 24,690 humpback units recorded in the month of October 2008. Under certain conditions, the probability of detection can be estimated with uncertainties sufficiently small to allow for accurate density estimates.