Similar Risk of Kidney Failure among Patients with Blinding Diseases Who Receive Ranibizumab, Aflibercept, and Bevacizumab: An Observational Health Data Sciences and Informatics Network Study.

PURPOSE: To characterize the incidence of kidney failure associated with intravitreal anti-VEGF exposure; and compare the risk of kidney failure in patients treated with ranibizumab, aflibercept, or bevacizumab. DESIGN: Retrospective cohort study across 12 databases in the Observational Health Data Sciences and Informatics (OHDSI) network. SUBJECTS: Subjects aged ≥ 18 years with ≥ 3 monthly intravitreal anti-VEGF medications for a blinding disease (diabetic retinopathy, diabetic macular edema, exudative age-related macular degeneration, or retinal vein occlusion). METHODS: The standardized incidence proportions and rates of kidney failure while on treatment with anti-VEGF were calculated. For each comparison (e.g., aflibercept versus ranibizumab), patients from each group were matched 1:1 using propensity scores. Cox proportional hazards models were used to estimate the risk of kidney failure while on treatment. A random effects meta-analysis was performed to combine each database's hazard ratio (HR) estimate into a single network-wide estimate. MAIN OUTCOME MEASURES: Incidence of kidney failure while on anti-VEGF treatment, and time from cohort entry to kidney failure. RESULTS: Of the 6.1 million patients with blinding diseases, 37 189 who received ranibizumab, 39 447 aflibercept, and 163 611 bevacizumab were included; the total treatment exposure time was 161 724 person-years. The average standardized incidence proportion of kidney failure was 678 per 100 000 persons (range, 0-2389), and incidence rate 742 per 100 000 person-years (range, 0-2661). The meta-analysis HR of kidney failure comparing aflibercept with ranibizumab was 1.01 (95% confidence interval [CI], 0.70-1.47; P = 0.45), ranibizumab with bevacizumab 0.95 (95% CI, 0.68-1.32; P = 0.62), and aflibercept with bevacizumab 0.95 (95% CI, 0.65-1.39; P = 0.60). CONCLUSIONS: There was no substantially different relative risk of kidney failure between those who received ranibizumab, bevacizumab, or aflibercept. Practicing ophthalmologists and nephrologists should be aware of the risk of kidney failure among patients receiving intravitreal anti-VEGF medications and that there is little empirical evidence to preferentially choose among the specific intravitreal anti-VEGF agents. FINANCIAL DISCLOSURES: Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Ranitidine Use and Incident Cancer in a Multinational Cohort.

Importance: Ranitidine, the most widely used histamine-2 receptor antagonist (H2RA), was withdrawn because of N-nitrosodimethylamine impurity in 2020. Given the worldwide exposure to this drug, the potential risk of cancer development associated with the intake of known carcinogens is an important epidemiological concern. Objective: To examine the comparative risk of cancer associated with the use of ranitidine vs other H2RAs. Design, Setting, and Participants: This new-user active comparator international network cohort study was conducted using 3 health claims and 9 electronic health record databases from the US, the United Kingdom, Germany, Spain, France, South Korea, and Taiwan. Large-scale propensity score (PS) matching was used to minimize confounding of the observed covariates with negative control outcomes. Empirical calibration was performed to account for unobserved confounding. All databases were mapped to a common data model. Database-specific estimates were combined using random-effects meta-analysis. Participants included individuals aged at least 20 years with no history of cancer who used H2RAs for more than 30 days from January 1986 to December 2020, with a 1-year washout period. Data were analyzed from April to September 2021. Exposure: The main exposure was use of ranitidine vs other H2RAs (famotidine, lafutidine, nizatidine, and roxatidine). Main Outcomes and Measures: The primary outcome was incidence of any cancer, except nonmelanoma skin cancer. Secondary outcomes included all cancer except thyroid cancer, 16 cancer subtypes, and all-cause mortality. Results: Among 1 183 999 individuals in 11 databases, 909 168 individuals (mean age, 56.1 years; 507 316 [55.8%] women) were identified as new users of ranitidine, and 274 831 individuals (mean age, 58.0 years; 145 935 [53.1%] women) were identified as new users of other H2RAs. Crude incidence rates of cancer were 14.30 events per 1000 person-years (PYs) in ranitidine users and 15.03 events per 1000 PYs among other H2RA users. After PS matching, cancer risk was similar in ranitidine compared with other H2RA users (incidence, 15.92 events per 1000 PYs vs 15.65 events per 1000 PYs; calibrated meta-analytic hazard ratio, 1.04; 95% CI, 0.97-1.12). No significant associations were found between ranitidine use and any secondary outcomes after calibration. Conclusions and Relevance: In this cohort study, ranitidine use was not associated with an increased risk of cancer compared with the use of other H2RAs. Further research is needed on the long-term association of ranitidine with cancer development.

Newly identified lineages of porcine hemagglutinating encephalomyelitis virus exhibit respiratory phenotype.

Swine pathogens have a long history of zoonotic transmission to humans, occasionally leading to sustained outbreaks or pandemics. Through a retrospective epidemiological study of swine populations in China, we describe novel lineages of porcine hemagglutinating encephalomyelitis virus (PHEV) complex coronaviruses (CoVs) that cause exclusively respiratory symptoms with no signs of the neurological symptoms typically associated with classical PHEV infection. Through large-scale epidemiological surveillance, we show that these novel lineages have circulated in at least eight provinces in southeastern China. Phylogenetic and recombination analyses of twenty-four genomes identified two major viral lineages causing respiratory symptoms with extensive recombination within them, between them, and between classical PHEV and the novel respiratory variant PHEV (rvPHEV) lineages. Divergence times among the sampled lineages in the PHEV virus complex date back to 1886-1958 (mean estimate 1928), with the two major rvPHEV lineages separating approximately 20 years later. Many rvPHEV viruses show amino acid substitutions at the carbohydrate-binding site of hemagglutinin esterase (HE) and/or have lost the cysteine required for HE dimerization. This resembles the early adaptation of human CoVs, where HE lost its hemagglutination ability to adapt to growth in the human respiratory tract. Our study represents the first report of the evolutionary history of rvPHEV circulating in swine and highlights the importance of characterizing CoV diversity and recombination in swine to identify pathogens with outbreak potential that could threaten swine farming.

Individualising intensive systolic blood pressure reduction in hypertension using computational trial phenomaps and machine learning: a post-hoc analysis of randomised clinical trials.

BACKGROUND: The cardiovascular benefits of intensive systolic blood pressure control vary across clinical populations tested in large randomised clinical trials. We aimed to evaluate the application of machine learning to clinical trials of patients without and with type 2 diabetes to define the personalised cardiovascular benefit of intensive control of systolic blood pressure. METHODS: In SPRINT, a trial of intensive (systolic blood pressure <120 mm Hg) versus standard (systolic blood pressure <140 mm Hg) systolic blood pressure control in patients without type 2 diabetes, we defined a phenotypic representation of the study population using 59 baseline variables. We extracted personalised treatment effect estimates for the primary outcome, time-to-first major adverse cardiovascular event (MACE; cardiovascular death, myocardial infarction or acute coronary syndrome, stroke, and acute decompensated heart failure), through iterative Cox regression analyses providing average hazard ratio (HR) estimates weighted for the phenotypic distance of each participant from the index patient of each iteration. Next, we trained an extreme gradient boosting algorithm (known as XGBoost) to predict the personalised effect of intensive systolic blood pressure control using features most consistently linked to increased personalised benefit, before evaluating its performance in the ACCORD BP trial of patients with type 2 diabetes randomly assigned to receive intensive versus standard systolic blood pressure control. We stratified patients based on their predicted treatment effect, and key demographic groups (age, sex, cardiovascular disease, and smoking). We assessed the presence of heterogeneity with an interaction test, and assessed the performance of the algorithm in a simulation analysis of SPRINT in the presence or absence of an artificially introduced heterogeneous treatment effect. FINDINGS: From SPRINT, we included all 9361 study participants (mean age 67·9 years [SD 9·4], 3332 [35·6%] female) who underwent randomisation to either intensive (n=4678) or standard (n=4683) treatment. The median individualised HR for MACE was 0·63 (IQR 0·53-0·78). An eight-feature tool built for this analysis to predict personalised benefit in SPRINT was externally tested in ACCORD BP (4733 participants (mean age 62·7 years [SD 6·7], 2258 [47·7%] female), wherein it successfully identified individuals with differential benefit from intensive versus standard systolic blood pressure control (adjusted HR for MACE of 0·70 [95% CI 0·55-0·90] in individuals with above-median MACE benefit versus 1·05 [95% CI 0·84-1·32] for below-median predicted benefit; pinteraction=0·0184). Subgroup analysis based on age (<65 years: HR 0·89 [95% CI 0·71-1·12]; ≥65 years: 0·85 [0·67-1·09]), sex (male: 0·89 [0·72-1·10]; female: 0·85 [0·65-1·10]), established cardiovascular disease (no: 0·89 [0·70-1·14]; yes: 0·84 [0·67-1·06]), or active smoking (no: 0·85 [0·71-1·02]; yes: 1·01 [0·64-1·60]) did not identify groups with heterogeneity of treatment effect. In a simulation analysis of SPRINT, the proposed algorithm detected groups with heterogeneous treatment effects in the presence, but not absence, of simulated subgroup differences. INTERPRETATION: By use of machine learning to define an individual's personalised benefit through phenotypic representations of clinical trials, we created a practical tool for individualising the selection of intensive versus standard systolic blood pressure control in patients without and with type 2 diabetes. FUNDING: National Heart, Lung, and Blood Institute of the US National Institutes of Health.

Comparative risk of thrombosis with thrombocytopenia syndrome or thromboembolic events associated with different covid-19 vaccines: international network cohort study from five European countries and the US.

OBJECTIVE: To quantify the comparative risk of thrombosis with thrombocytopenia syndrome or thromboembolic events associated with use of adenovirus based covid-19 vaccines versus mRNA based covid-19 vaccines. DESIGN: International network cohort study. SETTING: Routinely collected health data from contributing datasets in France, Germany, the Netherlands, Spain, the UK, and the US. PARTICIPANTS: Adults (age ≥18 years) registered at any contributing database and who received at least one dose of a covid-19 vaccine (ChAdOx1-S (Oxford-AstraZeneca), BNT162b2 (Pfizer-BioNTech), mRNA-1273 (Moderna), or Ad26.COV2.S (Janssen/Johnson & Johnson)), from December 2020 to mid-2021. MAIN OUTCOME MEASURES: Thrombosis with thrombocytopenia syndrome or venous or arterial thromboembolic events within the 28 days after covid-19 vaccination. Incidence rate ratios were estimated after propensity scores matching and were calibrated using negative control outcomes. Estimates specific to the database were pooled by use of random effects meta-analyses. RESULTS: Overall, 1 332 719 of 3 829 822 first dose ChAdOx1-S recipients were matched to 2 124 339 of 2 149 679 BNT162b2 recipients from Germany and the UK. Additionally, 762 517 of 772 678 people receiving Ad26.COV2.S were matched to 2 851 976 of 7 606 693 receiving BNT162b2 in Germany, Spain, and the US. All 628 164 Ad26.COV2.S recipients from the US were matched to 2 230 157 of 3 923 371 mRNA-1273 recipients. A total of 862 thrombocytopenia events were observed in the matched first dose ChAdOx1-S recipients from Germany and the UK, and 520 events after a first dose of BNT162b2. Comparing ChAdOx1-S with a first dose of BNT162b2 revealed an increased risk of thrombocytopenia (pooled calibrated incidence rate ratio 1.33 (95% confidence interval 1.18 to 1.50) and calibrated incidence rate difference of 1.18 (0.57 to 1.8) per 1000 person years). Additionally, a pooled calibrated incidence rate ratio of 2.26 (0.93 to 5.52) for venous thrombosis with thrombocytopenia syndrome was seen with Ad26.COV2.S compared with BNT162b2. CONCLUSIONS: In this multinational study, a pooled 30% increased risk of thrombocytopenia after a first dose of the ChAdOx1-S vaccine was observed, as was a trend towards an increased risk of venous thrombosis with thrombocytopenia syndrome after Ad26.COV2.S compared with BNT162b2. Although rare, the observed risks after adenovirus based vaccines should be considered when planning further immunisation campaigns and future vaccine development.

Seek COVER: using a disease proxy to rapidly develop and validate a personalized risk calculator for COVID-19 outcomes in an international network.

BACKGROUND: We investigated whether we could use influenza data to develop prediction models for COVID-19 to increase the speed at which prediction models can reliably be developed and validated early in a pandemic. We developed COVID-19 Estimated Risk (COVER) scores that quantify a patient's risk of hospital admission with pneumonia (COVER-H), hospitalization with pneumonia requiring intensive services or death (COVER-I), or fatality (COVER-F) in the 30-days following COVID-19 diagnosis using historical data from patients with influenza or flu-like symptoms and tested this in COVID-19 patients. METHODS: We analyzed a federated network of electronic medical records and administrative claims data from 14 data sources and 6 countries containing data collected on or before 4/27/2020. We used a 2-step process to develop 3 scores using historical data from patients with influenza or flu-like symptoms any time prior to 2020. The first step was to create a data-driven model using LASSO regularized logistic regression, the covariates of which were used to develop aggregate covariates for the second step where the COVER scores were developed using a smaller set of features. These 3 COVER scores were then externally validated on patients with 1) influenza or flu-like symptoms and 2) confirmed or suspected COVID-19 diagnosis across 5 databases from South Korea, Spain, and the United States. Outcomes included i) hospitalization with pneumonia, ii) hospitalization with pneumonia requiring intensive services or death, and iii) death in the 30 days after index date. RESULTS: Overall, 44,507 COVID-19 patients were included for model validation. We identified 7 predictors (history of cancer, chronic obstructive pulmonary disease, diabetes, heart disease, hypertension, hyperlipidemia, kidney disease) which combined with age and sex discriminated which patients would experience any of our three outcomes. The models achieved good performance in influenza and COVID-19 cohorts. For COVID-19 the AUC ranges were, COVER-H: 0.69-0.81, COVER-I: 0.73-0.91, and COVER-F: 0.72-0.90. Calibration varied across the validations with some of the COVID-19 validations being less well calibrated than the influenza validations. CONCLUSIONS: This research demonstrated the utility of using a proxy disease to develop a prediction model. The 3 COVER models with 9-predictors that were developed using influenza data perform well for COVID-19 patients for predicting hospitalization, intensive services, and fatality. The scores showed good discriminatory performance which transferred well to the COVID-19 population. There was some miscalibration in the COVID-19 validations, which is potentially due to the difference in symptom severity between the two diseases. A possible solution for this is to recalibrate the models in each location before use.

Characteristics and Outcomes of Over 300,000 Patients with COVID-19 and History of Cancer in the United States and Spain.

BACKGROUND: We described the demographics, cancer subtypes, comorbidities, and outcomes of patients with a history of cancer and coronavirus disease 2019 (COVID-19). Second, we compared patients hospitalized with COVID-19 to patients diagnosed with COVID-19 and patients hospitalized with influenza. METHODS: We conducted a cohort study using eight routinely collected health care databases from Spain and the United States, standardized to the Observational Medical Outcome Partnership common data model. Three cohorts of patients with a history of cancer were included: (i) diagnosed with COVID-19, (ii) hospitalized with COVID-19, and (iii) hospitalized with influenza in 2017 to 2018. Patients were followed from index date to 30 days or death. We reported demographics, cancer subtypes, comorbidities, and 30-day outcomes. RESULTS: We included 366,050 and 119,597 patients diagnosed and hospitalized with COVID-19, respectively. Prostate and breast cancers were the most frequent cancers (range: 5%-18% and 1%-14% in the diagnosed cohort, respectively). Hematologic malignancies were also frequent, with non-Hodgkin's lymphoma being among the five most common cancer subtypes in the diagnosed cohort. Overall, patients were aged above 65 years and had multiple comorbidities. Occurrence of death ranged from 2% to 14% and from 6% to 26% in the diagnosed and hospitalized COVID-19 cohorts, respectively. Patients hospitalized with influenza (n = 67,743) had a similar distribution of cancer subtypes, sex, age, and comorbidities but lower occurrence of adverse events. CONCLUSIONS: Patients with a history of cancer and COVID-19 had multiple comorbidities and a high occurrence of COVID-19-related events. Hematologic malignancies were frequent. IMPACT: This study provides epidemiologic characteristics that can inform clinical care and etiologic studies.

Thirty-Day Outcomes of Children and Adolescents With COVID-19: An International Experience.

OBJECTIVES: To characterize the demographics, comorbidities, symptoms, in-hospital treatments, and health outcomes among children and adolescents diagnosed or hospitalized with coronavirus disease 2019 (COVID-19) and to compare them in secondary analyses with patients diagnosed with previous seasonal influenza in 2017-2018. METHODS: International network cohort using real-world data from European primary care records (France, Germany, and Spain), South Korean claims and US claims, and hospital databases. We included children and adolescents diagnosed and/or hospitalized with COVID-19 at age <18 between January and June 2020. We described baseline demographics, comorbidities, symptoms, 30-day in-hospital treatments, and outcomes including hospitalization, pneumonia, acute respiratory distress syndrome, multisystem inflammatory syndrome in children, and death. RESULTS: A total of 242 158 children and adolescents diagnosed and 9769 hospitalized with COVID-19 and 2 084 180 diagnosed with influenza were studied. Comorbidities including neurodevelopmental disorders, heart disease, and cancer were more common among those hospitalized with versus diagnosed with COVID-19. Dyspnea, bronchiolitis, anosmia, and gastrointestinal symptoms were more common in COVID-19 than influenza. In-hospital prevalent treatments for COVID-19 included repurposed medications (<10%) and adjunctive therapies: systemic corticosteroids (6.8%-7.6%), famotidine (9.0%-28.1%), and antithrombotics such as aspirin (2.0%-21.4%), heparin (2.2%-18.1%), and enoxaparin (2.8%-14.8%). Hospitalization was observed in 0.3% to 1.3% of the cohort diagnosed with COVID-19, with undetectable (n < 5 per database) 30-day fatality. Thirty-day outcomes including pneumonia and hypoxemia were more frequent in COVID-19 than influenza. CONCLUSIONS: Despite negligible fatality, complications including hospitalization, hypoxemia, and pneumonia were more frequent in children and adolescents with COVID-19 than with influenza. Dyspnea, anosmia, and gastrointestinal symptoms could help differentiate diagnoses. A wide range of medications was used for the inpatient management of pediatric COVID-19.

Baseline characteristics, management, and outcomes of 55,270 children and adolescents diagnosed with COVID-19 and 1,952,693 with influenza in France, Germany, Spain, South Korea and the United States: an international network cohort study.

Objectives To characterize the demographics, comorbidities, symptoms, in-hospital treatments, and health outcomes among children/adolescents diagnosed or hospitalized with COVID-19. Secondly, to describe health outcomes amongst children/adolescents diagnosed with previous seasonal influenza. Design International network cohort. Setting Real-world data from European primary care records (France/Germany/Spain), South Korean claims and US claims and hospital databases. Participants Diagnosed and/or hospitalized children/adolescents with COVID-19 at age <18 between January and June 2020; diagnosed with influenza in 2017-2018. Main outcome measures Baseline demographics and comorbidities, symptoms, 30-day in-hospital treatments and outcomes including hospitalization, pneumonia, acute respiratory distress syndrome (ARDS), multi-system inflammatory syndrome (MIS-C), and death. Results A total of 55,270 children/adolescents diagnosed and 3,693 hospitalized with COVID-19 and 1,952,693 diagnosed with influenza were studied. Comorbidities including neurodevelopmental disorders, heart disease, and cancer were all more common among those hospitalized vs diagnosed with COVID-19. The most common COVID-19 symptom was fever. Dyspnea, bronchiolitis, anosmia and gastrointestinal symptoms were more common in COVID-19 than influenza. In-hospital treatments for COVID-19 included repurposed medications (<10%), and adjunctive therapies: systemic corticosteroids (6.8% to 37.6%), famotidine (9.0% to 28.1%), and antithrombotics such as aspirin (2.0% to 21.4%), heparin (2.2% to 18.1%), and enoxaparin (2.8% to 14.8%). Hospitalization was observed in 0.3% to 1.3% of the COVID-19 diagnosed cohort, with undetectable (N<5 per database) 30-day fatality. Thirty-day outcomes including pneumonia, ARDS, and MIS-C were more frequent in COVID-19 than influenza. Conclusions Despite negligible fatality, complications including pneumonia, ARDS and MIS-C were more frequent in children/adolescents with COVID-19 than with influenza. Dyspnea, anosmia and gastrointestinal symptoms could help differential diagnosis. A wide range of medications were used for the inpatient management of pediatric COVID-19.

Association of Ticagrelor vs Clopidogrel With Net Adverse Clinical Events in Patients With Acute Coronary Syndrome Undergoing Percutaneous Coronary Intervention.

Importance: Current guidelines recommend ticagrelor as the preferred P2Y12 platelet inhibitor for patients with acute coronary syndrome (ACS), primarily based on a single large randomized clinical trial. The benefits and risks associated with ticagrelor vs clopidogrel in routine practice merits attention. Objective: To determine the association of ticagrelor vs clopidogrel with ischemic and hemorrhagic events in patients undergoing percutaneous coronary intervention (PCI) for ACS in clinical practice. Design, Setting, and Participants: A retrospective cohort study of patients with ACS who underwent PCI and received ticagrelor or clopidogrel was conducted using 2 United States electronic health record-based databases and 1 nationwide South Korean database from November 2011 to March 2019. Patients were matched using a large-scale propensity score algorithm, and the date of final follow-up was March 2019. Exposures: Ticagrelor vs clopidogrel. Main Outcomes and Measures: The primary end point was net adverse clinical events (NACE) at 12 months, composed of ischemic events (recurrent myocardial infarction, revascularization, or ischemic stroke) and hemorrhagic events (hemorrhagic stroke or gastrointestinal bleeding). Secondary outcomes included NACE or mortality, all-cause mortality, ischemic events, hemorrhagic events, individual components of the primary outcome, and dyspnea at 12 months. The database-level hazard ratios (HRs) were pooled to calculate summary HRs by random-effects meta-analysis. Results: After propensity score matching among 31 290 propensity-matched pairs (median age group, 60-64 years; 29.3% women), 95.5% of patients took aspirin together with ticagrelor or clopidogrel. The 1-year risk of NACE was not significantly different between ticagrelor and clopidogrel (15.1% [3484/23 116 person-years] vs 14.6% [3290/22 587 person-years]; summary HR, 1.05 [95% CI, 1.00-1.10]; P = .06). There was also no significant difference in the risk of all-cause mortality (2.0% for ticagrelor vs 2.1% for clopidogrel; summary HR, 0.97 [95% CI, 0.81-1.16]; P = .74) or ischemic events (13.5% for ticagrelor vs 13.4% for clopidogrel; summary HR, 1.03 [95% CI, 0.98-1.08]; P = .32). The risks of hemorrhagic events (2.1% for ticagrelor vs 1.6% for clopidogrel; summary HR, 1.35 [95% CI, 1.13-1.61]; P = .001) and dyspnea (27.3% for ticagrelor vs 22.6% for clopidogrel; summary HR, 1.21 [95% CI, 1.17-1.26]; P < .001) were significantly higher in the ticagrelor group. Conclusions and Relevance: Among patients with ACS who underwent PCI in routine clinical practice, ticagrelor, compared with clopidogrel, was not associated with significant difference in the risk of NACE at 12 months. Because the possibility of unmeasured confounders cannot be excluded, further research is needed to determine whether ticagrelor is more effective than clopidogrel in this setting.

The emergence of SARS-CoV-2 in Europe and North America.

Accurate understanding of the global spread of emerging viruses is critical for public health responses and for anticipating and preventing future outbreaks. Here we elucidate when, where, and how the earliest sustained severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission networks became established in Europe and North America. Our results suggest that rapid early interventions successfully prevented early introductions of the virus from taking hold in Germany and the United States. Other, later introductions of the virus from China to both Italy and Washington state, United States, founded the earliest sustained European and North America transmission networks. Our analyses demonstrate the effectiveness of public health measures in preventing onward transmission and show that intensive testing and contact tracing could have prevented SARS-CoV-2 outbreaks from becoming established in these regions.

Comparative safety and effectiveness of alendronate versus raloxifene in women with osteoporosis.

Alendronate and raloxifene are among the most popular anti-osteoporosis medications. However, there is a lack of head-to-head comparative effectiveness studies comparing the two treatments. We conducted a retrospective large-scale multicenter study encompassing over 300 million patients across nine databases encoded in the Observational Medical Outcomes Partnership (OMOP) Common Data Model (CDM). The primary outcome was the incidence of osteoporotic hip fracture, while secondary outcomes were vertebral fracture, atypical femoral fracture (AFF), osteonecrosis of the jaw (ONJ), and esophageal cancer. We used propensity score trimming and stratification based on an expansive propensity score model with all pre-treatment patient characteritistcs. We accounted for unmeasured confounding using negative control outcomes to estimate and adjust for residual systematic bias in each data source. We identified 283,586 alendronate patients and 40,463 raloxifene patients. There were 7.48 hip fracture, 8.18 vertebral fracture, 1.14 AFF, 0.21 esophageal cancer and 0.09 ONJ events per 1,000 person-years in the alendronate cohort and 6.62, 7.36, 0.69, 0.22 and 0.06 events per 1,000 person-years, respectively, in the raloxifene cohort. Alendronate and raloxifene have a similar hip fracture risk (hazard ratio [HR] 1.03, 95% confidence interval [CI] 0.94-1.13), but alendronate users are more likely to have vertebral fractures (HR 1.07, 95% CI 1.01-1.14). Alendronate has higher risk for AFF (HR 1.51, 95% CI 1.23-1.84) but similar risk for esophageal cancer (HR 0.95, 95% CI 0.53-1.70), and ONJ (HR 1.62, 95% CI 0.78-3.34). We demonstrated substantial control of measured confounding by propensity score adjustment, and minimal residual systematic bias through negative control experiments, lending credibility to our effect estimates. Raloxifene is as effective as alendronate and may remain an option in the prevention of osteoporotic fracture.

Positive Selection in CD8+ T-Cell Epitopes of Influenza Virus Nucleoprotein Revealed by a Comparative Analysis of Human and Swine Viral Lineages.

UNLABELLED: Numerous experimental studies have demonstrated that CD8(+) T cells contribute to immunity against influenza by limiting viral replication. It is therefore surprising that rigorous statistical tests have failed to find evidence of positive selection in the epitopes targeted by CD8(+) T cells. Here we use a novel computational approach to test for selection in CD8(+) T-cell epitopes. We define all epitopes in the nucleoprotein (NP) and matrix protein (M1) with experimentally identified human CD8(+) T-cell responses and then compare the evolution of these epitopes in parallel lineages of human and swine influenza viruses that have been diverging since roughly 1918. We find a significant enrichment of substitutions that alter human CD8(+) T-cell epitopes in NP of human versus swine influenza virus, consistent with the idea that these epitopes are under positive selection. Furthermore, we show that epitope-altering substitutions in human influenza virus NP are enriched on the trunk versus the branches of the phylogenetic tree, indicating that viruses that acquire these mutations have a selective advantage. However, even in human influenza virus NP, sites in T-cell epitopes evolve more slowly than do nonepitope sites, presumably because these epitopes are under stronger inherent functional constraint. Overall, our work demonstrates that there is clear selection from CD8(+) T cells in human influenza virus NP and illustrates how comparative analyses of viral lineages from different hosts can identify positive selection that is otherwise obscured by strong functional constraint. IMPORTANCE: There is a strong interest in correlates of anti-influenza immunity that are protective against diverse virus strains. CD8(+) T cells provide such broad immunity, since they target conserved viral proteins. An important question is whether T-cell immunity is sufficiently strong to drive influenza virus evolution. Although many studies have shown that T cells limit viral replication in animal models and are associated with decreased symptoms in humans, no studies have proven with statistical significance that influenza virus evolves under positive selection to escape T cells. Here we use comparisons of human and swine influenza viruses to rigorously demonstrate that human influenza virus evolves under pressure to fix mutations in the nucleoprotein that promote escape from T cells. We further show that viruses with these mutations have a selective advantage since they are preferentially located on the "trunk" of the phylogenetic tree. Overall, our results show that CD8(+) T cells targeting nucleoprotein play an important role in shaping influenza virus evolution.

Estimation for general birth-death processes.

Birth-death processes (BDPs) are continuous-time Markov chains that track the number of "particles" in a system over time. While widely used in population biology, genetics and ecology, statistical inference of the instantaneous particle birth and death rates remains largely limited to restrictive linear BDPs in which per-particle birth and death rates are constant. Researchers often observe the number of particles at discrete times, necessitating data augmentation procedures such as expectation-maximization (EM) to find maximum likelihood estimates. For BDPs on finite state-spaces, there are powerful matrix methods for computing the conditional expectations needed for the E-step of the EM algorithm. For BDPs on infinite state-spaces, closed-form solutions for the E-step are available for some linear models, but most previous work has resorted to time-consuming simulation. Remarkably, we show that the E-step conditional expectations can be expressed as convolutions of computable transition probabilities for any general BDP with arbitrary rates. This important observation, along with a convenient continued fraction representation of the Laplace transforms of the transition probabilities, allows for novel and efficient computation of the conditional expectations for all BDPs, eliminating the need for truncation of the state-space or costly simulation. We use this insight to derive EM algorithms that yield maximum likelihood estimation for general BDPs characterized by various rate models, including generalized linear models. We show that our Laplace convolution technique outperforms competing methods when they are available and demonstrate a technique to accelerate EM algorithm convergence. We validate our approach using synthetic data and then apply our methods to cancer cell growth and estimation of mutation parameters in microsatellite evolution.

High-dimensional, massive sample-size Cox proportional hazards regression for survival analysis.

Survival analysis endures as an old, yet active research field with applications that spread across many domains. Continuing improvements in data acquisition techniques pose constant challenges in applying existing survival analysis methods to these emerging data sets. In this paper, we present tools for fitting regularized Cox survival analysis models on high-dimensional, massive sample-size (HDMSS) data using a variant of the cyclic coordinate descent optimization technique tailored for the sparsity that HDMSS data often present. Experiments on two real data examples demonstrate that efficient analyses of HDMSS data using these tools result in improved predictive performance and calibration.

The spread of type 2 Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) in North America: a phylogeographic approach.

The emergence and spread of Type 2 Porcine Reproductive and Respiratory Syndrome virus (Type 2 PRRSV) in North America is heavily influenced by the multiple site production system used in the hog industry. However, it is unclear how anthropogenic factors such has this have shaped the current spatial distribution of PRRSV genotypes. We employed Bayesian phylogeographic analyses of 7040 ORF5 sequences to reveal the recent geographical spread of Type 2 PRRSV in North America. The directions and intensities in our inferred virus traffic network closely mirror the hog transportation. Most notably, we reveal multiple viral introductions from Canada into the United States causing a major shift in virus genetic composition in the Midwest USA that went unnoticed by the regular surveillance and field epidemiological studies. Overall, these findings provide important insights into the dynamics of Type 2 PRRSV evolution and spread that will facilitate programs for control and prevention.

Phylogeography of dengue virus serotype 4, Brazil, 2010-2011.

Dengue virus serotype 4 (DENV-4) reemerged in Roraima State, Brazil, 28 years after it was last detected in the country in 1982. To study the origin and evolution of this reemergence, full-length sequences were obtained for 16 DENV-4 isolates from northern (Roraima, Amazonas, Pará States) and northeastern (Bahia State) Brazil during the 2010 and 2011 dengue virus seasons and for an isolate from the 1982 epidemic in Roraima. Spatiotemporal dynamics of DENV-4 introductions in Brazil were applied to envelope genes and full genomes by using Bayesian phylogeographic analyses. An introduction of genotype I into Brazil from Southeast Asia was confirmed, and full genome phylogeographic analyses revealed multiple introductions of DENV-4 genotype II in Brazil, providing evidence for >3 introductions of this genotype within the last decade: 2 from Venezuela to Roraima and 1 from Colombia to Amazonas. The phylogeographic analysis of full genome data has demonstrated the origins of DENV-4 throughout Brazil.

Phylogeography and population dynamics of dengue viruses in the Americas.

Changes in Dengue virus (DENV) disease patterns in the Americas over recent decades have been attributed, at least in part, to repeated introduction of DENV strains from other regions, resulting in a shift from hypoendemicity to hyperendemicity. Using newly sequenced DENV-1 and DENV-3 envelope (E) gene isolates from 11 Caribbean countries, along with sequences available on GenBank, we sought to document the population genetic and spatiotemporal transmission histories of the four main invading DENV genotypes within the Americas and investigate factors that influence the rate and intensity of DENV transmission. For all genotypes, there was an initial invasion phase characterized by rapid increases in genetic diversity, which coincided with the first confirmed cases of each genotype in the region. Rapid geographic dispersal occurred upon each genotype's introduction, after which individual lineages were locally maintained, and gene flow was primarily observed among neighboring and nearby countries. There were, however, centers of viral diversity (Barbados, Puerto Rico, Colombia, Suriname, Venezuela, and Brazil) that were repeatedly involved in gene flow with more distant locations. For DENV-1 and DENV-2, we found that a "distance-informed" model, which posits that the intensity of virus movement between locations is inversely proportional to the distance between them, provided a better fit than a model assuming equal rates of movement between all pairs of countries. However, for DENV-3 and DENV-4, the more stochastic "equal rates" model was preferred.

A Bayesian model for the common effects of multiple predictors on mixed outcomes.

We propose a Bayesian multivariate model in which a single linear combination of the covariates predict multiple outcomes simultaneously. The single linear combination is a data-derived score along the lines of the Apache or Charlson index scores for critically ill patients, the Karnofsky or Eastern Cooperative Oncology Group score for cancer patients or Euro-score for cardiac patients that may be used to predict multiple outcomes. Outcomes may be discrete or continuous and we use a composition of generalized linear models for the marginal distribution for each outcome. We explain how to set the prior distribution and we use Markov chain Monte Carlo methods to calculate the posterior distribution. We propose two types of expanded models to diagnose whether each outcome indeed has predictor effects common with the other outcomes, and whether a particular predictor is commonly predictive for all outcomes. We determine a final model based on the diagnostic models. The method is applied to a study yielding multiple psychometric outcomes of mixed type measured in young people living with human immunodeficiency virus.