Search

Under stress, bacteria display selective growth propensities on solid and liquid media. For Mycobacterium tuberculosis (Mtb), differentially culturable bacteria have been found in tuberculosis (TB) patient sputa. We hypothesized that antibiotic treatment can induce selective culturability in Mtb. We investigated the effect of exposure to TB drugs on Mtb culturability using clinical samples from an ongoing TB patient cohort and by conducting several in vitro experiments with a diverse set of Mtb strains. In patients, serial sputa were more likely to generate Mtb-positive cultures in liquid as opposed to solid medium, with this liquid culture bias extending up to 5 months post diagnosis. Experimentally, there was a disparity between bacterial time to positivity (TTP) and colony forming units (CFUs) when Mtb was exposed to isoniazid (INH) and rifampicin (RIF) alone or in combination. Cultures recovered from RIF treatment yielded more CFUs on agar plates, but INH-treated cultures had a faster TTP in liquid. Follow up experiments using a fluorescently labelled laboratory strain of Mtb revealed that CFUs of INH-treated bacteria were too low, thus over-estimating the killing effect of the treatment. Here we provide evidence that drug exposure affects culturability on solid medium in diverse Mtb strains, which has implications for treatment monitoring and drug-pathogen interaction studies., Competing Interests: Declarations. Competing interests: The authors declare no competing interests., (© 2025. The Author(s).)

Background: Mycobacterium tuberculosis complex (MTBC) species evolve slowly, so isolates from individuals linked in transmission often have identical or nearly identical genomes, making it difficult to reconstruct transmission chains. Finding additional sources of shared MTBC variation could help overcome this problem. Previous studies have reported MTBC diversity within infected individuals; however, whether within-host variation improves transmission inferences remains unclear. Here, we aimed to quantify within-host MTBC variation and assess whether such information improves transmission inferences., Methods: We conducted a retrospective genomic epidemiology study in which we reanalysed publicly available sequence data from household transmission studies published in PubMed from database inception until Jan 31, 2024, for which both genomic and epidemiological contact data were available, using household membership as a proxy for transmission linkage. We quantified minority variants (ie, positions with two or more alleles each supported by at least five-fold coverage and with a minor allele frequency of 1% or more) outside of PE and PPE genes, within individual samples and shared across samples. We used receiver operator characteristic (ROC) curves to compare the performance of a general linear model for household membership that included shared minority variants and one that included only fixed genetic differences., Findings: We identified three MTBC household transmission studies with publicly available whole-genome sequencing data and epidemiological linkages: a household transmission study in Vitória, Brazil (Colangeli et al), a retrospective population-based study of paediatric tuberculosis in British Columbia, Canada (Guthrie et al), and a retrospective population-based study in Oxfordshire, England (Walker et al). We found moderate levels of minority variation present in MTBC sequence data from cultured isolates that varied significantly across studies: mean 168·6 minority variants (95% CI 151·4-185·9) for the Colangeli et al dataset, 5·8 (1·5-10·2) for Guthrie et al (p<0·0001, Wilcoxon rank sum test, vs Colangeli et al), and 7·1 (2·4-11·9) for Walker et al (p<0·0001, Wilcoxon rank sum test, vs Colangeli et al). Isolates from household pairs shared more minority variants than did randomly selected pairs of isolates: mean 97·7 shared minority variants (79·1-116·3) versus 9·8 (8·6-11·0) in Colangeli et al, 0·8 (0·1-1·5) versus 0·2 (0·1-0·2) in Guthrie et al, and 0·7 (0·1-1·3) versus 0·2 (0·2-0·2) in Walker et al (all p<0·0001, Wilcoxon rank sum test). Shared within-host variation was significantly associated with household membership (odds ratio 1·51 [95% CI 1·30-1·71], p<0·0001), for one standard deviation increase in shared minority variants. Models that included shared within-host variation versus models without within-host variation improved the accuracy of predicting household membership in all three studies: area under the ROC curve 0·95 versus 0·92 for the Colangeli et al study, 0·99 versus 0·95 for the Guthrie et al study, and 0·93 versus 0·91 for the Walker et al study., Interpretation: Within-host MTBC variation persists through culture of sputum and could enhance the resolution of transmission inferences. The substantial differences in minority variation recovered across studies highlight the need to optimise approaches to recover and incorporate within-host variation into automated phylogenetic and transmission inference., Funding: National Institutes of Health., Competing Interests: Declaration of interests JRA reports support from the US National Institutes of Health (NIH); receiving donated Cepheid disposable open cartridges from Cepheid to Stanford University for NIH-funded tuberculosis research; participation in two Data Safety Monitoring Boards for NIH-funded tuberculosis trials; and participation in the Science Advisory Board for an NIH-funded tuberculosis network. JC reports grants from the NIH, Valneva/Butantan Institution, Coalition for Epidemic Preparedness Innovations–Sabin Institute, MSD, and Sanofi Pasteur for clinical and epidemiological studies; honoraria from Pfizer; participation in the Brazil Advisory Board for the mRNA-1273 vaccine Moderna–Zodiac, the Latin America and Brazilian Advisory Board for paxlovid (nirmatrelvir–ritonavir; Pfizer), the Brazil Advisory Board for the Qdenga vaccine (Takeda), and the Takeda Global Dengue Steering Committee. All other authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Culturing and genomic sequencing of Mycobacterium tuberculosis (MTB) from tuberculosis (TB) cases is the basis for many research and clinical applications. The alternative, culture-free sequencing from diagnostic samples, is promising but poses challenges to obtain and analyse the MTB genome. Paradoxically, culture is assumed to impose a diversity bottleneck, which, if true, would entail unexplored consequences. To unravel this paradox we generate high-quality genomes of sputum-culture pairs from two different settings after developing a workflow for sequencing from sputum and a tailored bioinformatics analysis. Careful downstream comparisons reveal sources of sputum-culture incongruences due to false positive/negative variation associated with factors like low input MTB DNA or variable genomic depths. After accounting for these factors, contrary to the bottleneck dogma, we identify a 97% variant agreement within sputum-culture pairs, with a high correlation also in the variants' frequency (0.98). The combined analysis from five different settings and more than 100 available samples shows that our results can be extrapolated to different TB epidemic scenarios, demonstrating that for the cases tested culture accurately mirrors clinical samples., (© 2024. The Author(s).)

Persons living with HIV are known to be at increased risk of developing tuberculosis (TB) disease upon infection with Mycobacterium tuberculosis (Mtb). However, it has remained unclear how HIV co-infection affects subsequent Mtb transmission from these patients. Here, we customized a Bayesian phylodynamic framework to estimate the effects of HIV co-infection on the Mtb transmission dynamics from sequence data. We applied our model to four Mtb genomic datasets collected in sub-Saharan African countries with a generalized HIV epidemic. Our results confirm that HIV co-infection is a strong risk factor for developing active TB. Additionally, we demonstrate that HIV co-infection is associated with a reduced effective reproductive number for TB. Stratifying the population by CD4+ T-cell count yielded similar results, suggesting that, in this context, CD4+ T-cell count is not a better predictor of Mtb transmissibility than HIV infection status alone. Together, our genome-based analyses complement observational household contact studies, and more firmly establish the negative association between HIV co-infection and Mtb transmissibility., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Windels et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Background: Because M. tuberculosis evolves slowly, transmission clusters often contain multiple individuals with identical consensus genomes, making it difficult to reconstruct transmission chains. Finding additional sources of shared M. tuberculosis variation could help overcome this problem. Previous studies have reported M. tuberculosis diversity within infected individuals; however, whether within-host variation improves transmission inferences remains unclear., Methods: To evaluate the transmission information present in within-host M. tuberculosis variation, we re-analyzed publicly available sequence data from three household transmission studies, using household membership as a proxy for transmission linkage between donor-recipient pairs., Findings: We found moderate levels of minority variation present in M. tuberculosis sequence data from cultured isolates that varied significantly across studies (mean: 6, 7, and 170 minority variants above a 1% minor allele frequency threshold, outside of PE/PPE genes). Isolates from household members shared more minority variants than did isolates from unlinked individuals in the three studies (mean 98 shared minority variants vs. 10; 0.8 vs. 0.2, and 0.7 vs. 0.2, respectively). Shared within-host variation was significantly associated with household membership (OR: 1.51 [1.30,1.71], for one standard deviation increase in shared minority variants). Models that included shared within-host variation improved the accuracy of predicting household membership in all three studies as compared to models without within-host variation (AUC: 0.95 versus 0.92, 0.99 versus 0.95, and 0.93 versus 0.91)., Interpretation: Within-host M. tuberculosis variation persists through culture and could enhance the resolution of transmission inferences. The substantial differences in minority variation recovered across studies highlights the need to optimize approaches to recover and incorporate within-host variation into automated phylogenetic and transmission inference., Funding: NIAID: 5K01AI173385., Competing Interests: Declaration of interest The authors report no conflict of interest.

Background: Experimental data show that drug-resistance-conferring mutations are often associated with a decrease in the replicative fitness of bacteria in vitro, and that this fitness cost can be mitigated by compensatory mutations; however, the role of compensatory evolution in clinical settings is less clear. We assessed whether compensatory evolution was associated with increased transmission of rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa., Methods: We did a genomic epidemiological study by analysing available M tuberculosis isolates and their associated clinical data from individuals routinely diagnosed with rifampicin-resistant tuberculosis in primary care and hospitals in Khayelitsha, Cape Town, South Africa. Isolates were collected as part of a previous study. All individuals diagnosed with rifampicin-resistant tuberculosis and with linked biobanked specimens were included in this study. We applied whole-genome sequencing, Bayesian reconstruction of transmission trees, and phylogenetic multivariable regression analysis to identify individual and bacterial factors associated with the transmission of rifampicin-resistant M tuberculosis strains., Findings: Between Jan 1, 2008, and Dec 31, 2017, 2161 individuals were diagnosed with multidrug-resistant or rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa. Whole-genome sequences were available for 1168 (54%) unique individual M tuberculosis isolates. Compensatory evolution was associated with smear-positive pulmonary disease (adjusted odds ratio 1·49, 95% CI 1·08-2·06) and a higher number of drug-resistance-conferring mutations (incidence rate ratio 1·38, 95% CI 1·28-1·48). Compensatory evolution was also associated with increased transmission of rifampicin-resistant disease between individuals (adjusted odds ratio 1·55; 95% CI 1·13-2·12), independent of other patient and bacterial factors., Interpretation: Our findings suggest that compensatory evolution enhances the in vivo fitness of drug-resistant M tuberculosis genotypes, both within and between patients, and that the in vitro replicative fitness of rifampicin-resistant M tuberculosis measured in the laboratory correlates with the bacterial fitness measured in clinical settings. These results emphasise the importance of enhancing surveillance and monitoring efforts to prevent the emergence of highly transmissible clones capable of rapidly accumulating new drug resistance mutations. This concern becomes especially crucial at present, because treatment regimens incorporating novel drugs are being implemented., Funding: Funding for this study was provided by a Swiss and South Africa joint research award (grant numbers 310030&#95;188888, CRSII5&#95;177163, and IZLSZ3&#95;170834), the European Research Council (grant number 883582), and a Wellcome Trust fellowship (to HC; reference number 099818/Z/12/Z). ZS-D was funded through a PhD scholarship from the South African National Research Foundation and RMW was funded through the South African Medical Research Council., Competing Interests: Declaration of interest RMW acknowledges baseline support from the South African Medical Research Council. All other authors declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY 4.0 license. Published by Elsevier Ltd.. All rights reserved.)

Treatment of multidrug-resistant or rifampicin-resistant tuberculosis (MDR/RR-TB), although improved in recent years with shorter, more tolerable regimens, remains largely standardized and based on limited drug susceptibility testing (DST). More individualized treatment with expanded DST access is likely to improve patient outcomes. To assess the potential of TB drug resistance prediction based on whole-genome sequencing (WGS) to provide more effective treatment regimens, we applied current South African treatment recommendations to a retrospective cohort of MDR/RR-TB patients from Khayelitsha, Cape Town. Routine DST and clinical data were used to retrospectively categorize patients into a recommended regimen, either a standardized short regimen or a longer individualized regimen. Potential regimen changes were then described with the addition of WGS-derived DST. WGS data were available for 1274 MDR/RR-TB patient treatment episodes across 2008 to 2017. Among 834 patients initially eligible for the shorter regimen, 385 (46%) may have benefited from reduced drug dosage or removing ineffective drugs when WGS data were considered. A further 187 (22%) patients may have benefited from more effective adjusted regimens. Among 440 patients initially eligible for a longer individualized regimen, 153 (35%) could have been switched to the short regimen. Overall, 305 (24%) patients had MDR/RR-TB with second-line TB drug resistance, where the availability of WGS-derived DST would have allowed more effective treatment individualization. These data suggest considerable benefits could accrue from routine access to WGS-derived resistance prediction. Advances in culture-free sequencing and expansion of the reference resistance mutation catalogue will increase the utility of WGS resistance prediction.

A rapid and accurate diagnostic assay represents an important means to detect Mycobacterium tuberculosis , identify drug-resistant strains and ensure treatment success. Currently employed techniques to diagnose drug-resistant tuberculosis include slow phenotypic tests or more rapid molecular assays that evaluate a limited range of drugs. Whole-genome-sequencing-based approaches can detect known drug-resistance-conferring mutations and novel variations; however, the dependence on growing samples in culture, and the associated delays in achieving results, represents a significant limitation. As an alternative, targeted sequencing strategies can be directly performed on clinical samples at high throughput. This study proposes a targeted sequencing assay to rapidly detect drug-resistant strains of M. tuberculosis using the Nanopore MinION sequencing platform. We designed a single-tube assay that targets nine genes associated with drug resistance to seven drugs and two phylogenetic-determining regions to determine strain lineage and tested it in nine clinical isolates and six sputa. The study's main aim is to calibrate MinNION variant calling to detect drug-resistance-associated mutations with different frequencies to match the accuracy of Illumina (the current gold-standard sequencing technology) from both culture and sputum samples. After calibrating Nanopore MinION variant calling, we demonstrated 100% agreement between Illumina WGS and our MinION set up to detect known drug resistance and phylogenetic variants in our dataset. Importantly, other variants in the amplicons are also detected, decreasing the recall. We identify minority variants and insertions/deletions as crucial bioinformatics challenges to fully reproduce Illumina WGS results.

Background: South Africa has a high burden of rifampicin-resistant tuberculosis (including multidrug-resistant [MDR] tuberculosis), with increasing rifampicin-monoresistant (RMR) tuberculosis over time. Resistance acquisition during first-line tuberculosis treatment could be a key contributor to this burden, and HIV might increase the risk of acquiring rifampicin resistance. We assessed whether HIV during previous treatment was associated with RMR tuberculosis and resistance acquisition among a retrospective cohort of patients with MDR or rifampicin-resistant tuberculosis., Methods: In this retrospective cohort study, we included all patients routinely diagnosed with MDR or rifampicin-resistant tuberculosis in Khayelitsha, Cape Town, South Africa, between Jan 1, 2008, and Dec 31, 2017. Patient-level data were obtained from a prospective database, complemented by data on previous tuberculosis treatment and HIV from a provincial health data exchange. Stored MDR or rifampicin-resistant tuberculosis isolates from patients underwent whole-genome sequencing (WGS). WGS data were used to infer resistance acquisition versus transmission, by identifying genomically unique isolates (single nucleotide polymorphism threshold of five). Logistic regression analyses were used to assess factors associated with RMR tuberculosis and genomic uniqueness., Findings: The cohort included 2041 patients diagnosed with MDR or rifampicin-resistant tuberculosis between Jan 1, 2008, and Dec 31, 2017; of those, 463 (22·7%) with RMR tuberculosis and 1354 (66·3%) with previous tuberculosis treatment. In previously treated patients, HIV positivity during previous tuberculosis treatment versus HIV negativity (adjusted odds ratio [OR] 2·07, 95% CI 1·35-3·18), and three or more previous tuberculosis treatment episodes versus one (1·96, 1·21-3·17) were associated with RMR tuberculosis. WGS data showing MDR or rifampicin-resistant tuberculosis were available for 1169 patients; 360 (30·8%) isolates were identified as unique. In previously treated patients, RMR tuberculosis versus MDR tuberculosis (adjusted OR 4·96, 3·40-7·23), HIV positivity during previous tuberculosis treatment (1·71, 1·03-2·84), and diagnosis in 2013-17 (1·42, 1·02-1·99) versus 2008-12, were associated with uniqueness. In previously treated patients with RMR tuberculosis, HIV positivity during previous treatment (adjusted OR 5·13, 1·61-16·32) was associated with uniqueness as was female sex (2·50 [1·18-5·26])., Interpretation: These data suggest that HIV contributes to rifampicin-resistance acquisition during first-line tuberculosis treatment and that this might be driving increasing RMR tuberculosis over time. Large-scale prospective cohort studies are required to further quantify this risk., Funding: Swiss National Science Foundation, South African National Research Foundation, and Wellcome Trust., Competing Interests: HC reports research grants from the US National Institutes of Health and the European and Developing Countries Clinical Trials Partnership outside the submitted work. All other authors declare no competing interests., (© 2021 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY 4.0 license.)

The coronavirus disease 2019 (COVID-19) pandemic has affected the world radically since 2020. Spain was one of the European countries with the highest incidence during the first wave. As a part of a consortium to monitor and study the evolution of the epidemic, we sequenced 2,170 samples, diagnosed mostly before lockdown measures. Here, we identified at least 500 introductions from multiple international sources and documented the early rise of two dominant Spanish epidemic clades (SECs), probably amplified by superspreading events. Both SECs were related closely to the initial Asian variants of SARS-CoV-2 and spread widely across Spain. We inferred a substantial reduction in the effective reproductive number of both SECs due to public-health interventions (R e  < 1), also reflected in the replacement of SECs by a new variant over the summer of 2020. In summary, we reveal a notable difference in the initial genetic makeup of SARS-CoV-2 in Spain compared with other European countries and show evidence to support the effectiveness of lockdown measures in controlling virus spread, even for the most successful genetic variants., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Background: Direct whole-genome sequencing of Mycobacterium tuberculosis from clinical specimens will be a major breakthrough in tuberculosis diagnosis and control. To date, direct whole-genome sequencing has never been used in genomic epidemiology, and its accuracy in transmission inference remains unknown. We investigated the technical challenges imposed by direct whole-genome sequencing, and used it to infer transmission clusters and predict drug resistance., Methods: Using an optimised workflow, we did direct whole-genome sequencing for 37 clinical specimens from 23 tuberculosis patients. Nine sputum samples from nine patients who were infected with different non-tuberculous mycobacteria and culture-negative for tuberculosis were used as controls in the qPCR assays and pre-sequencing runs. Additionally, 780 clinical isolates in the region of Comunidad Valenciana (Spain) were whole-genome sequenced between Jan 1, 2014, and Dec 31, 2016. We analysed the genomic variants to build a tuberculosis transmission network for the region, including the clinical specimens, and to predict drug susceptibility profiles., Findings: After sequencing 37 clinical specimens, 28 specimens (22 [85%] of 26 smear-positive and six [55%] of 11 smear-negative) met the quality criteria for downstream analysis. All 28 clinical specimens clustered with their matching culture isolates, with a median distance of 0 single nucleotide polymorphisms. Of the 28 clinical specimens, 16 (57%) were accurately assigned to ten transmission clusters in the region, and 12 (43%) were unique cases. Transmission inferences and drug-susceptibility predictions from direct whole-genome sequencing data were concordant with sequences from corresponding cultures and phenotypic drug-susceptibility testing. Complete genomic analysis, within a week of specimen receipt, cost €217 per sample (excluding personnel costs)., Interpretation: Direct whole-genome sequencing could be used to accurately delineate transmission clusters of tuberculosis and conduct culture-independent surveillance. Compared with conventional approaches, direct whole-genome sequencing allows researchers to do real-time genomic epidemiology and drug resistance surveillance in settings where culture and drug susceptibility testing are not available., Funding: European Research Council; Ministerio de Ciencia, Innovación y Universidades (Spanish Government)., (Copyright © 2020 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY 4.0 license. Published by Elsevier Ltd.. All rights reserved.)

The administration of a high fat content diet is an accelerating factor for metabolic syndrome, impaired glucose tolerance, and early type 2 diabetes. The present study aims to assess the impact of a high fat diet on tuberculosis progression and microbiota composition in an experimental animal model using a C3HeB/FeJ mouse strain submitted to single or multiple consecutive aerosol infections. These models allowed us to study the protection induced by Bacillus Calmette-Guérin vaccination as well as by the natural immunity induced by chemotherapy after a low dose Mycobacterium tuberculosis infection. Our results show that a high fat diet is able to trigger a pro-inflammatory response, which results in a faster progression toward active tuberculosis and an impaired protective effect of BCG vaccination, which is not the case for natural immunity. This may be related to dysbiosis and a reduction in the Firmicutes/Bacteroidetes ratio in the gut microbiota caused by a decrease in the abundance of the Porphyromonadaceae family and, in particular, the Barnesiella genus. It should also be noted that a high fat diet is also related to an increase in the genera Alistipes, Parasuterella, Mucispirillum , and Akkermansia , which have previously been related to dysbiotic processes. As diabetes mellitus type 2 is a risk factor for developing tuberculosis, these findings may prove useful in the search for new prophylactic strategies for this population subset., (Copyright © 2019 Arias, Goig, Cardona, Torres-Puente, Díaz, Rosales, Garcia, Tapia, Comas, Vilaplana and Cardona.)

Bacillus thuringiensis parasporal crystal proteins (Cry proteins) are insecticidal pore-forming toxins that bind to specific receptor molecules on the brush border membrane of susceptible insect midgut cells to exert their toxic action. In the Colorado potato beetle (CPB), a coleopteran pest, we previously proposed that interaction of Cry3Aa toxin with a CPB ADAM10 metalloprotease is an essential part of the mode of action of this toxin. Here, we annotated the gene sequence encoding an ADAM10 metalloprotease protein (CPB-ADAM10) in the CPB genome sequencing project, and using RNA interference gene silencing we demonstrated that CPB-ADAM10 is a Cry3Aa toxin functional receptor in CPB. Cry3Aa toxicity was significantly lower in CPB-ADAM10 silenced larvae and in vitro toxin pore-forming ability was greatly diminished in lipid planar bilayers fused with CPB brush border membrane vesicles (BBMVs) prepared from CPB-ADAM10 silenced larvae. In accordance with our previous data that indicated this toxin was a substrate of ADAM10 in CPB, Cry3Aa toxin membrane-associated proteolysis was altered when CPB BBMVs lacked ADAM10. The functional validation of CPB-ADAM10 as a Cry3Aa toxin receptor in CPB expands the already recognized role of ADAM10 as a pathogenicity determinant of pore-forming toxins in humans to an invertebrate species., (© 2016 The Royal Entomological Society.)

Colonial medical reports claimed that tuberculosis (TB) was largely unknown in Africa prior to European contact, providing a "virgin soil" for spread of TB in highly susceptible populations previously unexposed to the disease [1, 2]. This is in direct contrast to recent phylogenetic models which support an African origin for TB [3-6]. To address this apparent contradiction, we performed a broad genomic sampling of Mycobacterium tuberculosis in Ethiopia. All members of the M. tuberculosis complex (MTBC) arose from clonal expansion of a single common ancestor [7] with a proposed origin in East Africa [3, 4, 8]. Consistent with this proposal, MTBC lineage 7 is almost exclusively found in that region [9-11]. Although a detailed medical history of Ethiopia supports the view that TB was rare until the 20(th) century [12], over the last century Ethiopia has become a high-burden TB country [13]. Our results provide further support for an African origin for TB, with some genotypes already present on the continent well before European contact. Phylogenetic analyses reveal a pattern of serial introductions of multiple genotypes into Ethiopia in association with human migration and trade. In place of a "virgin soil" fostering the spread of TB in a previously naive population, we propose that increased TB mortality in Africa was driven by the introduction of European strains of M. tuberculosis alongside expansion of selected indigenous strains having biological characteristics that carry a fitness benefit in the urbanized settings of post-colonial Africa., (Copyright © 2015 The Authors. Published by Elsevier Ltd.. All rights reserved.)