Discovering study-specific gene regulatory networks

This article has been made available through the Brunel Open Access Publishing Fund.Microarrays are commonly used in biology because of their ability to simultaneously measure thousands of genes under different conditions. Due to their structure, typically containing a high amount of variables but far fewer samples, scalable network analysis techniques are often employed. In particular, consensus approaches have been recently used that combine multiple microarray studies in order to find networks that are more robust. The purpose of this paper, however, is to combine multiple microarray studies to automatically identify subnetworks that are distinctive to specific experimental conditions rather than common to them all. To better understand key regulatory mechanisms and how they change under different conditions, we derive unique networks from multiple independent networks built using glasso which goes beyond standard correlations. This involves calculating cluster prediction accuracies to detect the most predictive genes for a specific set of conditions. We differentiate between accuracies calculated using cross-validation within a selected cluster of studies (the intra prediction accuracy) and those calculated on a set of independent studies belonging to different study clusters (inter prediction accuracy). Finally, we compare our method's results to related state-of-the art techniques. We explore how the proposed pipeline performs on both synthetic data and real data (wheat and Fusarium). Our results show that subnetworks can be identified reliably that are specific to subsets of studies and that these networks reflect key mechanisms that are fundamental to the experimental conditions in each of those subsets

Expression of the BCR-ABL1 Gene in Patients with Chronic Myeloproliferative Diseases with Signs of Progression

Background. The V617F mutation of JAK2 is known to manifest in Ph-negative chronic myeloproliferative diseases (cMPD), such as polycythemia vera, thrombocythemia, and myelofibrosis. These diseases not infrequently advance into more aggressive forms up to acute leukemia. As the progression mechanism is still unknown, its study retains a high priority. JAK2 carrying the V617F mutation is believed to cause constant activation of V(D)J recombinase in myeloid tumor cells in cMPD patients. Aberrant activation of V(D)J recombinase in tumor cells in cMPD patients can lead to t(9;22)(q34;q11) chromosomal rearrangement. Aim. To study the expression of BCR-ABL1 resulting from translocation t(9;22)(q34;q11) in cMPD patients at the progression stage in order to test the suggested hypothesis. Materials & Methods. The BCR–ABL1 expression was assessed in peripheral blood granulocytes in cMPD patients by real-time PCR. The JAK2 V617F mutation was identified by quantitative allele-specific PCR. The JAK2 exon 12 mutations were determined using Sanger direct sequencing of PCR products. Results. The BCR-ABL1 expression was discovered in 29 % of patients with cMPD progression. The BCR-ABL1 expression in these patients correlated with hepatosplenomegaly and hyperleukocytosis. Conclusion. In a significant proportion of cMPD patients the disease progression can be associated with activation of the BCR-ABL expression

Global genome diversity of the Leishmania donovani complex.

Protozoan parasites of the Leishmania donovani complex - L. donovani and L. infantum - cause the fatal disease visceral leishmaniasis. We present the first comprehensive genome-wide global study, with 151 cultured field isolates representing most of the geographical distribution. L. donovani isolates separated into five groups that largely coincide with geographical origin but vary greatly in diversity. In contrast, the majority of L. infantum samples fell into one globally-distributed group with little diversity. This picture is complicated by several hybrid lineages. Identified genetic groups vary in heterozygosity and levels of linkage, suggesting different recombination histories. We characterise chromosome-specific patterns of aneuploidy and identified extensive structural variation, including known and suspected drug resistance loci. This study reveals greater genetic diversity than suggested by geographically-focused studies, provides a resource of genomic variation for future work and sets the scene for a new understanding of the evolution and genetics of the Leishmania donovani complex

Nontypable Haemophilus influenzae Displays a Prevalent Surface Structure Molecular Pattern in Clinical Isolates

Non-typable Haemophilus influenzae (NTHi) is a Gram negative pathogen that causes acute respiratory infections and is associated with the progression of chronic respiratory diseases. Previous studies have established the existence of a remarkable genetic variability among NTHi strains. In this study we show that, in spite of a high level of genetic heterogeneity, NTHi clinical isolates display a prevalent molecular feature, which could confer fitness during infectious processes. A total of 111 non-isogenic NTHi strains from an identical number of patients, isolated in two distinct geographical locations in the same period of time, were used to analyse nine genes encoding bacterial surface molecules, and revealed the existence of one highly prevalent molecular pattern (lgtF+, lic2A+, lic1D+, lic3A+, lic3B+, siaA−, lic2C+, ompP5+, oapA+) displayed by 94.6% of isolates. Such a genetic profile was associated with a higher bacterial resistance to serum mediated killing and enhanced adherence to human respiratory epithelial cells

A new world malaria map: Plasmodium falciparum endemicity in 2010

Background: transmission intensity affects almost all aspects of malaria epidemiology and the impact of malaria on human populations. Maps of transmission intensity are necessary to identify populations at different levels of risk and to evaluate objectively options for disease control. To remain relevant operationally, such maps must be updated frequently. Following the first global effort to map Plasmodium falciparum malaria endemicity in 2007, this paper describes the generation of a new world map for the year 2010. This analysis is extended to provide the first global estimates of two other metrics of transmission intensity for P. falciparum that underpin contemporary questions in malaria control: the entomological inoculation rate (PfEIR) and the basic reproductive number (PfR). Methods: annual parasite incidence data for 13,449 administrative units in 43 endemic countries were sourced to define the spatial limits of P. falciparum transmission in 2010 and 22,212 P. falciparum parasite rate (PfPR) surveys were used in a model-based geostatistical (MBG) prediction to create a continuous contemporary surface of malaria endemicity within these limits. A suite of transmission models were developed that link PfPR to PfEIR and PfR and these were fitted to field data. These models were combined with the PfPR map to create new global predictions of PfEIR and PfR. All output maps included measured uncertainty. Results: an estimated 1.13 and 1.44 billion people worldwide were at risk of unstable and stable P. falciparum malaria, respectively. The majority of the endemic world was predicted with a median PfEIR of less than one and a median PfRc of less than two. Values of either metric exceeding 10 were almost exclusive to Africa. The uncertainty described in both PfEIR and PfR was substantial in regions of intense transmission. Conclusions: the year 2010 has a particular significance as an evaluation milestone for malaria global health policy. The maps presented here contribute to a rational basis for control and elimination decisions and can serve as a baseline assessment as the global health community looks ahead to the next series of milestones targeted at 20

Sarilumab in patients admitted to hospital with severe or critical COVID-19: a randomised, double-blind, placebo-controlled, phase 3 trial

Background: Elevated proinflammatory cytokines are associated with greater COVID-19 severity. We aimed to assess safety and efficacy of sarilumab, an interleukin-6 receptor inhibitor, in patients with severe (requiring supplemental oxygen by nasal cannula or face mask) or critical (requiring greater supplemental oxygen, mechanical ventilation, or extracorporeal support) COVID-19. Methods: We did a 60-day, randomised, double-blind, placebo-controlled, multinational phase 3 trial at 45 hospitals in Argentina, Brazil, Canada, Chile, France, Germany, Israel, Italy, Japan, Russia, and Spain. We included adults (≥18 years) admitted to hospital with laboratory-confirmed SARS-CoV-2 infection and pneumonia, who required oxygen supplementation or intensive care. Patients were randomly assigned (2:2:1 with permuted blocks of five) to receive intravenous sarilumab 400 mg, sarilumab 200 mg, or placebo. Patients, care providers, outcome assessors, and investigators remained masked to assigned intervention throughout the course of the study. The primary endpoint was time to clinical improvement of two or more points (seven point scale ranging from 1 [death] to 7 [discharged from hospital]) in the modified intention-to-treat population. The key secondary endpoint was proportion of patients alive at day 29. Safety outcomes included adverse events and laboratory assessments. This study is registered with ClinicalTrials.gov, NCT04327388; EudraCT, 2020-001162-12; and WHO, U1111-1249-6021. Findings: Between March 28 and July 3, 2020, of 431 patients who were screened, 420 patients were randomly assigned and 416 received placebo (n=84 [20%]), sarilumab 200 mg (n=159 [38%]), or sarilumab 400 mg (n=173 [42%]). At day 29, no significant differences were seen in median time to an improvement of two or more points between placebo (12·0 days [95% CI 9·0 to 15·0]) and sarilumab 200 mg (10·0 days [9·0 to 12·0]; hazard ratio [HR] 1·03 [95% CI 0·75 to 1·40]; log-rank p=0·96) or sarilumab 400 mg (10·0 days [9·0 to 13·0]; HR 1·14 [95% CI 0·84 to 1·54]; log-rank p=0·34), or in proportions of patients alive (77 [92%] of 84 patients in the placebo group; 143 [90%] of 159 patients in the sarilumab 200 mg group; difference −1·7 [−9·3 to 5·8]; p=0·63 vs placebo; and 159 [92%] of 173 patients in the sarilumab 400 mg group; difference 0·2 [−6·9 to 7·4]; p=0·85 vs placebo). At day 29, there were numerical, non-significant survival differences between sarilumab 400 mg (88%) and placebo (79%; difference +8·9% [95% CI −7·7 to 25·5]; p=0·25) for patients who had critical disease. No unexpected safety signals were seen. The rates of treatment-emergent adverse events were 65% (55 of 84) in the placebo group, 65% (103 of 159) in the sarilumab 200 mg group, and 70% (121 of 173) in the sarilumab 400 mg group, and of those leading to death 11% (nine of 84) were in the placebo group, 11% (17 of 159) were in the sarilumab 200 mg group, and 10% (18 of 173) were in the sarilumab 400 mg group. Interpretation: This trial did not show efficacy of sarilumab in patients admitted to hospital with COVID-19 and receiving supplemental oxygen. Adequately powered trials of targeted immunomodulatory therapies assessing survival as a primary endpoint are suggested in patients with critical COVID-19. Funding: Sanofi and Regeneron Pharmaceuticals

Multi-messenger observations of a binary neutron star merger

On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ~1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40+8-8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 Mo. An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ~40 Mpc) less than 11 hours after the merger by the One- Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ~10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ~9 and ~16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta

Multi-messenger Observations of a Binary Neutron Star Merger

On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ∼ 1.7 {{s}} with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of {40}-8+8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 {M}ȯ . An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ∼ 40 {{Mpc}}) less than 11 hours after the merger by the One-Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ∼10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ∼ 9 and ∼ 16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC 4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta.</p