Elucidating Emergence and Transmission of Multidrug-Resistant Tuberculosis in Treatment Experienced Patients by Whole Genome Sequencing

Background

Understanding the emergence and spread of multidrug-resistant tuberculosis (MDR-TB) is crucial for its control. MDR-TB in previously treated patients is generally attributed to the selection of drug resistant mutants during inadequate therapy rather than transmission of a resistant strain. Traditional genotyping methods are not sufficient to distinguish strains in populations with a high burden of tuberculosis and it has previously been difficult to assess the degree of transmission in these settings. We have used whole genome analysis to investigate M. tuberculosis strains isolated from treatment experienced patients with MDR-TB in Uganda over a period of four years.

Methods and Findings

We used high throughput genome sequencing technology to investigate small polymorphisms and large deletions in 51 Mycobacterium tuberculosis samples from 41 treatment-experienced TB patients attending a TB referral and treatment clinic in Kampala. This was a convenience sample representing 69% of MDR-TB cases identified over the four year period. Low polymorphism was observed in longitudinal samples from individual patients (2-15 SNPs). Clusters of samples with less than 50 SNPs variation were examined. Three clusters comprising a total of 8 patients were found with almost identical genetic profiles, including mutations predictive for resistance to rifampicin and isoniazid, suggesting transmission of MDR-TB. Two patients with previous drug susceptible disease were found to have acquired MDR strains, one of which shared its genotype with an isolate from another patient in the cohort.

Conclusions

Whole genome sequence analysis identified MDR-TB strains that were shared by more than one patient. The transmission of multidrug-resistant disease in this cohort of retreatment patients emphasises the importance of early detection and need for infection control. Consideration should be given to rapid testing for drug resistance in patients undergoing treatment to monitor the emergence of resistance and permit early intervention to avoid onward transmission.     

Unsupervised Clustering of Subcellular Protein Expression Patterns in High-Throughput Microscopy Images Reveals Protein Complexes and Functional Relationships between Proteins

Protein subcellular localization has been systematically characterized in budding yeast using fluorescently tagged proteins. Based on the fluorescence microscopy images, subcellular localization of many proteins can be classified automatically using supervised machine learning approaches that have been trained to recognize predefined image classes based on statistical features. Here, we present an unsupervised analysis of protein expression patterns in a set of high-resolution, high-throughput microscope images. Our analysis is based on 7 biologically interpretable features which are evaluated on automatically identified cells, and whose cell-stage dependency is captured by a continuous model for cell growth. We show that it is possible to identify most previously identified localization patterns in a cluster analysis based on these features and that similarities between the inferred expression patterns contain more information about protein function than can be explained by a previous manual categorization of subcellular localization. Furthermore, the inferred cell-stage associated to each fluorescence measurement allows us to visualize large groups of proteins entering the bud at specific stages of bud growth. These correspond to proteins localized to organelles, revealing that the organelles must be entering the bud in a stereotypical order. We also identify and organize a smaller group of proteins that show subtle differences in the way they move around the bud during growth. Our results suggest that biologically interpretable features based on explicit models of cell morphology will yield unprecedented power for pattern discovery in high-resolution, high-throughput microscopy images.     

Weight-Bearing Locomotion in the Developing Opossum,Monodelphis domestica following Spinal Transection: Remodeling of Neuronal Circuits Caudal to Lesion

Complete spinal transection in the mature nervous system is typically followed by minimal axonal repair, extensive motor paralysis and loss of sensory functions caudal to the injury. In contrast, the immature nervous system has greater capacity for repair, a phenomenon sometimes called the infant lesion effect. This study investigates spinal injuries early in development using the marsupial opossum Monodelphis domestica whose young are born very immature, allowing access to developmental stages only accessible in utero in eutherian mammals. Spinal cords of Monodelphis pups were completely transected in the lower thoracic region, T10, on postnatal-day (P)7 or P28 and the animals grew to adulthood. In P7-injured animals regrown supraspinal and propriospinal axons through the injury site were demonstrated using retrograde axonal labelling. These animals recovered near-normal coordinated overground locomotion, but with altered gait characteristics including foot placement phase lags. In P28-injured animals no axonal regrowth through the injury site could be demonstrated yet they were able to perform weight-supporting hindlimb stepping overground and on the treadmill. When placed in an environment of reduced sensory feedback (swimming) P7-injured animals swam using their hindlimbs, suggesting that the axons that grew across the lesion made functional connections; P28-injured animals swam using their forelimbs only, suggesting that their overground hindlimb movements were reflex-dependent and thus likely to be generated locally in the lumbar spinal cord. Modifications to propriospinal circuitry in P7- and P28-injured opossums were demonstrated by changes in the number of fluorescently labelled neurons detected in the lumbar cord following tracer studies and changes in the balance of excitatory, inhibitory and neuromodulatory neurotransmitter receptors’ gene expression shown by qRT-PCR. These results are discussed in the context of studies indicating that although following injury the isolated segment of the spinal cord retains some capability of rhythmic movement the mechanisms involved in weight-bearing locomotion are distinct.     

Spatio-Temporal Distribution of Injured Elephants in Masai Mara and the Putative Negative and Positive Roles of the Local Community

Background

Very few studies have ever focused on the elephants that are wounded or killed as local communities attempt to scare these animals away from their settlements and farms, or on the cases in which local people take revenge after elephants have killed or injured humans. On the other hand, local communities live in close proximity to elephants and hence can play a positive role in elephant conservation by informing the authorities of the presence of injured elephants.

Methodology/Principal Findings

Between 2007 and 2011, 129 elephants were monitored in Masai Mara (Kenya), of which 54 had various types of active (intentionally caused) or passive (non-intentionally caused) injuries. Also studied were 75 random control samples of apparently unaffected animals. The observed active injuries were as expected biased by age, with adults suffering more harm; on the other hand, no such bias was observed in the case of passive injuries. Bias was also observed in elephant sex since more males than females were passively and actively injured. Cases of passive and active injuries in elephants were negatively related to the proximity to roads and farms; the distribution of injured elephants was not affected by the presence of either human settlements or water sources. Overall more elephants were actively injured during the dry season than the wet season as expected. Local communities play a positive role by informing KWS authorities of the presence of injured elephants and reported 43% of all cases of injured elephants.

Conclusions

Our results suggest that the negative effect of local communities on elephants could be predicted by elephant proximity to farms and roads. In addition, local communities may be able to play a more positive role in elephant conservation given that they are key informants in the early detection of injured elephants.     

Effects of Photophase and Altitude on Oviposition Rhythm of the Himalayan Strains of Drosophila Ananassae

The effects of varying photophase and altitude of origin on the phase angle difference (Ψ) of the circadian rhythm of oviposition during entrainment to light‐dark (LD) cycles and the aftereffects of such photophases on the period of the free‐running rhythm (τ) in constant darkness (DD) were evaluated in two Himalayan strains of Drosophila ananassae, the high‐altitude (HA) strain from Badrinath (5,123 m above sea level=ASL) and the low‐altitude (LA) strain from Firozpur (179 m ASL). The Ψ (i.e., the hours from lights‐on of the LD cycle to oviposition median) of both strains was determined in LD cycles in which the photophase at 100 lux varied from 6 to 18 h/24 h. The HA strain was entrained by all LD cycles except the one with 6 h photophase in which it was weakly rhythmic, but the LA strain was entrained by only three LD cycles with photophases of 10, 12, and 14 h, but photophases of 6, 8, 16, and 18 h rendered it arrhythmic. Lights‐off transition of LD cycles was the phase‐determining signal for both strains as oviposition medians of the HA strain occurred～6 h prior to lights‐off, while those of the LA strain occurred～1 h after lights‐off. The Ψ of the HA strain increased from～2 h in 8 h photophase to～11 h in 18 h photophase, while that of the LA strain increased from～11 h in 10 h photophase to～15 h in 14 h photophase. The aftereffects of photophase of the prior entraining LD cycles on τ in DD were determined by transferring flies from LD cycles to DD. The τ of the HA strain increased from～19 to～25 h when transferred to DD from LD 8:16 and LD 18:6 cycles, respectively, whereas the τ of the LA strain increased from～26 to～28 h when transferred to DD from LD 10:14 and LD 14:10 cycles, respectively. Thus, these results demonstrate that the photophases of entraining LD cycles and the altitude of origin affected several parameters of entrainment and the period of the free‐running rhythm of these strains.     

Molecular Characterization of Staphylococcus aureus from Patients with Surgical Site Infections at Mulago Hospital in Kampala,Uganda

Background

The prevalence of Methicillin resistant Staphylococcus aureus (MRSA) is progressively increasing globally with significant regional variation. Understanding the Staphylococcus aureus lineages is crucial in controlling nosocomial infections. Recent studies on S. aureus in Uganda have revealed an escalating burden of MRSA. However, the S. aureus genotypes circulating among patients are not known. Here, we report S. aureus lineages circulating in patients with surgical site infections (SSI) at Mulago National hospital, Kampala, Uganda.

Methods

A cross-sectional study involving 314 patients with SSI at Mulago National Hospital was conducted from September 2011 to April 2012. Pus swabs from the patients’ SSI were processed using standard microbiological procedures. Methicillin sensitive Staphylococcus aureus (MSSA) and MRSA were identified using phenotypic tests and confirmed by PCR-detection of the nuc and mecA genes, respectively. SCCmec genotypes were determined among MRSA isolates using multiplex PCR. Furthermore, to determine lineages, spa sequence based-genotyping was performed on all S. aureus isolates.

Results

Of the 314 patients with SSI, S. aureus accounted for 20.4% (64/314), of which 37.5% (24/64) were MRSA. The predominant SCCmec types were type V (33.3%, 8/24) and type I (16.7%, 4/24). The predominant spa lineages were t645 (17.2%, 11/64) and t4353 (15.6%, 10/64), and these were found to be clonally circulating in all the surgical wards. On the other hand, lineages t064, t355, and t4609 were confined to the obstetrics and gynecology wards. A new spa type (t10277) was identified from MSSA isolate. On multivariate logistic regression analysis, cancer and inducible clindamycin resistance remained as independent predictors of MRSA-SSI.

Conclusion

SCCmec types I and V are the most prevalent MRSA mecA types from the patients’ SSI. The predominant spa lineages (t645 and t4353) are clonally circulating in all the surgical wards, calling for strengthening of infection control practices at Mulago National Hospital.     

Mortality,Morbidity and Health-Seeking Behaviour during the Ebola Epidemic 2014–2015 in Monrovia Results from a Mobile Phone Survey

Between March 2014 and July 2015 at least 10,500 Ebola cases including more than 4,800 deaths occurred in Liberia, the majority in Monrovia. However, official numbers may have underestimated the size of the outbreak. Closure of health facilities and mistrust in existing structures may have additionally impacted on all-cause morbidity and mortality. To quantify mortality and morbidity and describe health-seeking behaviour in Monrovia, Médecins sans Frontières (MSF) conducted a mobile phone survey from December 2014 to March 2015. We drew a random sample of households in Monrovia and conducted structured mobile phone interviews, covering morbidity, mortality and health-seeking behaviour from 14 May 2014 until the day of the survey. We defined an Ebola-related death as any death meeting the Liberian Ebola case definition. We calculated all-cause and Ebola-specific mortality rates. The sample consisted of 6,813 household members in 905 households. We estimated a crude mortality rate (CMR) of 0.33/10,000 persons/day (95%CI:0.25–0.43) and an Ebola-specific mortality rate of 0.06/10,000 persons/day (95%-CI:0.03–0.11). During the recall period, 17 Ebola cases were reported including those who died. In the 30 days prior to the survey 277 household members were reported sick; malaria accounted for 54% (150/277). Of the sick household members, 43% (122/276) did not visit any health care facility. The mobile phone-based survey was found to be a feasible and acceptable alternative method when data collection in the community is impossible. CMR was estimated well below the emergency threshold of 1/10,000 persons/day. Non-Ebola-related mortality in Monrovia was not higher than previous national estimates of mortality for Liberia. However, excess mortality directly resulting from Ebola did occur in the population. Importantly, the small proportion of sick household members presenting to official health facilities when sick might pose a challenge for future outbreak detection and mitigation. Substantial reported health-seeking behaviour outside of health facilities may also suggest the need for adapted health messaging and improved access to health care.     

Development and Two-Year Follow-Up Evaluation of a Training Workshop for the Large Preventive Positive Psychology Happy Family Kitchen Project in Hong Kong

Evidence-based practice and capacity-building approaches are essential for large-scale health promotion interventions. However, there are few models in the literature to guide and evaluate training of social service workers in community settings. This paper presents the development and evaluation of the “train-the-trainer” workshop (TTT) for the first large scale, community-based, family intervention projects, entitled “Happy Family Kitchen Project” (HFK) under the FAMILY project, a Hong Kong Jockey Club Initiative for a Harmonious Society. The workshop aimed to enhance social workers’ competence and performance in applying positive psychology constructs in their family interventions under HFK to improve family well-being of the community they served. The two-day TTT was developed and implemented by a multidisciplinary team in partnership with community agencies to 50 social workers (64% women). It focused on the enhancement of knowledge, attitude, and practice of five specific positive psychology themes, which were the basis for the subsequent development of the 23 family interventions for 1419 participants. Acceptability and applicability were enhanced by completing a needs assessment prior to the training. The TTT was evaluated by trainees’ reactions to the training content and design, changes in learners (trainees) and benefits to the service organizations. Focus group interviews to evaluate the workshop at three months after the training, and questionnaire survey at pre-training, immediately after, six months, one year and two years after training were conducted. There were statistically significant increases with large to moderate effect size in perceived knowledge, self-efficacy and practice after training, which sustained to 2-year follow-up. Furthermore, there were statistically significant improvements in family communication and well-being of the participants in the HFK interventions they implemented after training. This paper offers a practical example of development, implementation and model-based evaluation of training programs, which may be helpful to others seeking to develop such programs in diverse communities.     

Midpoint attractors and species richness: Modelling the interaction between environmental drivers and geometric constraints

We introduce a novel framework for conceptualising, quantifying and unifying discordant patterns of species richness along geographical gradients. While not itself explicitly mechanistic, this approach offers a path towards understanding mechanisms. In this study, we focused on the diverse patterns of species richness on mountainsides. We conjectured that elevational range midpoints of species may be drawn towards a single midpoint attractor – a unimodal gradient of environmental favourability. The midpoint attractor interacts with geometric constraints imposed by sea level and the mountaintop to produce taxon‐specific patterns of species richness. We developed a Bayesian simulation model to estimate the location and strength of the midpoint attractor from species occurrence data sampled along mountainsides. We also constructed midpoint predictor models to test whether environmental variables could directly account for the observed patterns of species range midpoints. We challenged these models with 16 elevational data sets, comprising 4500 species of insects, vertebrates and plants. The midpoint predictor models generally failed to predict the pattern of species midpoints. In contrast, the midpoint attractor model closely reproduced empirical spatial patterns of species richness and range midpoints. Gradients of environmental favourability, subject to geometric constraints, may parsimoniously account for elevational and other patterns of species richness.