Subsurface Microbial Diversity in Deep-Granitic-Fracture Water in Colorado

A microbial community analysis using 16S rRNA gene sequencing was performed on borehole water and a granite rock core from Henderson Mine, a >1,000-meter-deep molybdenum mine near Empire, CO. Chemical analysis of borehole water at two separate depths (1,044 m and 1,004 m below the mine entrance) suggests that a sharp chemical gradient exists, likely from the mixing of two distinct subsurface fluids, one metal rich and one relatively dilute; this has created unique niches for microorganisms. The microbial community analyzed from filtered, oxic borehole water indicated an abundance of sequences from iron-oxidizing bacteria (Gallionella spp.) and was compared to the community from the same borehole after 2 weeks of being plugged with an expandable packer. Statistical analyses with UniFrac revealed a significant shift in community structure following the addition of the packer. Phospholipid fatty acid (PLFA) analysis suggested that Nitrosomonadales dominated the oxic borehole, while PLFAs indicative of anaerobic bacteria were most abundant in the samples from the plugged borehole. Microbial sequences were represented primarily by Firmicutes, Proteobacteria, and a lineage of sequences which did not group with any identified bacterial division; phylogenetic analyses confirmed the presence of a novel candidate division. This “Henderson candidate division” dominated the clone libraries from the dilute anoxic fluids. Sequences obtained from the granitic rock core (1,740 m below the surface) were represented by the divisions Proteobacteria (primarily the family Ralstoniaceae) and Firmicutes. Sequences grouping within Ralstoniaceae were also found in the clone libraries from metal-rich fluids yet were absent in more dilute fluids. Lineage-specific comparisons, combined with phylogenetic statistical analyses, show that geochemical variance has an important effect on microbial community structure in deep, subsurface systems.     

Pheromone-encoding mRNA is transported to the yeast mating projection by specific RNP granules

Association of messenger RNAs with large complexes such as processing bodies (PBs) plays a pivotal role in regulating their translation and decay. Little is known about other possible functions of these assemblies. Exposure of haploid yeast cells, carrying mating type “a,” to “α pheromone” stimulates polarized growth resulting in a “shmoo” projection; it also induces synthesis of “a pheromone,” encoded by MFA2. In this paper, we show that, in response to α pheromone, MFA2 mRNA is assembled with two types of granules; both contain some canonical PB proteins, yet they differ in size, localization, motility, and sensitivity to cycloheximide. Remarkably, one type is involved in mRNA transport to the tip of the shmoo, whereas the other—in local translation in the shmoo. Normal assembly of these granules is critical for their movement, localization, and for mating. Thus, MFA2 mRNAs are transported to the shmoo tip, in complex with PB-like particles, where they are locally translated.     

Quantitative and qualitative beta diversity measures lead to different insights into factors that structure microbial communities

The assessment of microbial diversity and distribution is a major concern in environmental microbiology. There are two general approaches for measuring community diversity: quantitative measures, which use the abundance of each taxon, and qualitative measures, which use only the presence/absence of data. Quantitative measures are ideally suited to revealing community differences that are due to changes in relative taxon abundance (e.g., when a particular set of taxa flourish because a limiting nutrient source becomes abundant). Qualitative measures are most informative when communities differ primarily by what can live in them (e.g., at high temperatures), in part because abundance information can obscure significant patterns of variation in which taxa are present. We illustrate these principles using two 16S rRNA-based surveys of microbial populations and two phylogenetic measures of community beta diversity: unweighted UniFrac, a qualitative measure, and weighted UniFrac, a new quantitative measure, which we have added to the UniFrac website (http://bmf.colorado.edu/unifrac). These studies considered the relative influences of mineral chemistry, temperature, and geography on microbial community composition in acidic thermal springs in Yellowstone National Park and the influences of obesity and kinship on microbial community composition in the mouse gut. We show that applying qualitative and quantitative measures to the same data set can lead to dramatically different conclusions about the main factors that structure microbial diversity and can provide insight into the nature of community differences. We also demonstrate that both weighted and unweighted UniFrac measurements are robust to the methods used to build the underlying phylogeny.     

The meaning of significant mean group differences for biomarker discovery

Over the past decade, biomarker discovery has become a key goal in psychiatry to aid in the more reliable diagnosis and prognosis of heterogeneous psychiatric conditions and the development of tailored therapies. Nevertheless, the prevailing statistical approach is still the mean group comparison between “cases” and “controls,” which tends to ignore within-group variability. In this educational article, we used empirical data simulations to investigate how effect size, sample size, and the shape of distributions impact the interpretation of mean group differences for biomarker discovery. We then applied these statistical criteria to evaluate biomarker discovery in one area of psychiatric research—autism research. Across the most influential areas of autism research, effect size estimates ranged from small (d = 0.21, anatomical structure) to medium (d = 0.36 electrophysiology, d = 0.5, eye-tracking) to large (d = 1.1 theory of mind). We show that in normal distributions, this translates to approximately 45% to 63% of cases performing within 1 standard deviation (SD) of the typical range, i.e., they do not have a deficit/atypicality in a statistical sense. For a measure to have diagnostic utility as defined by 80% sensitivity and 80% specificity, Cohen’s d of 1.66 is required, with still 40% of cases falling within 1 SD. However, in both normal and nonnormal distributions, 1 (skewness) or 2 (platykurtic, bimodal) biologically plausible subgroups may exist despite small or even nonsignificant mean group differences. This conclusion drastically contrasts the way mean group differences are frequently reported. Over 95% of studies omitted the “on average” when summarising their findings in their abstracts (“autistic people have deficits in X”), which can be misleading as it implies that the group-level difference applies to all individuals in that group. We outline practical approaches and steps for researchers to explore mean group comparisons for the discovery of stratification biomarkers.     

An Antidote to the Imager's Fallacy,or How to Identify Brain Areas That Are in Limbo

Traditionally, fMRI data are analyzed using statistical parametric mapping approaches. Regardless of the precise thresholding procedure, these approaches ultimately divide the brain in regions that do or do not differ significantly across experimental conditions. This binary classification scheme fosters the so-called imager''s fallacy, where researchers prematurely conclude that region A is selectively involved in a certain cognitive task because activity in that region reaches statistical significance and activity in region B does not. For such a conclusion to be statistically valid, however, a test on the differences in activation across these two regions is required. Here we propose a simple GLM-based method that defines an “in-between” category of brain regions that are neither significantly active nor inactive, but rather “in limbo”. For regions that are in limbo, the activation pattern is inconclusive: it does not differ significantly from baseline, but neither does it differ significantly from regions that do show significant changes from baseline. This pattern indicates that measurement was insufficiently precise. By directly testing differences in activation, our procedure helps reduce the impact of the imager''s fallacy. The method is illustrated using concrete examples.     

Confidence Intervals for fMRI Activation Maps

Neuroimaging activation maps typically color voxels to indicate whether the blood oxygen level-dependent (BOLD) signals measured among two or more experimental conditions differ significantly at that location. This data presentation, however, omits information critical for interpretation of experimental results. First, no information is represented about trends at voxels that do not pass the statistical test. Second, no information is given about the range of probable effect sizes at voxels that do pass the statistical test. This leads to a fundamental error in interpreting activation maps by naïve viewers, where it is assumed that colored, “active” voxels are reliably different from uncolored “inactive” voxels. In other domains, confidence intervals have been added to data graphics to reduce such errors. Here, we first document the prevalence of the fundamental error of interpretation, and then present a method for solving it by depicting confidence intervals in fMRI activation maps. Presenting images where the bounds of confidence intervals at each voxel are coded as color allows readers to visually test for differences between “active” and “inactive” voxels, and permits for more proper interpretation of neuroimaging data. Our specific graphical methods are intended as initial proposals to spur broader discussion of how to present confidence intervals for fMRI data.     

Analysis of 567,758 randomized controlled trials published over 30 years reveals trends in phrases used to discuss results that do not reach statistical significance

The power of language to modify the reader’s perception of interpreting biomedical results cannot be underestimated. Misreporting and misinterpretation are pressing problems in randomized controlled trials (RCT) output. This may be partially related to the statistical significance paradigm used in clinical trials centered around a P value below 0.05 cutoff. Strict use of this P value may lead to strategies of clinical researchers to describe their clinical results with P values approaching but not reaching the threshold to be “almost significant.” The question is how phrases expressing nonsignificant results have been reported in RCTs over the past 30 years. To this end, we conducted a quantitative analysis of English full texts containing 567,758 RCTs recorded in PubMed between 1990 and 2020 (81.5% of all published RCTs in PubMed). We determined the exact presence of 505 predefined phrases denoting results that approach but do not cross the line of formal statistical significance (P < 0.05). We modeled temporal trends in phrase data with Bayesian linear regression. Evidence for temporal change was obtained through Bayes factor (BF) analysis. In a randomly sampled subset, the associated P values were manually extracted. We identified 61,741 phrases in 49,134 RCTs indicating almost significant results (8.65%; 95% confidence interval (CI): 8.58% to 8.73%). The overall prevalence of these phrases remained stable over time, with the most prevalent phrases being “marginally significant” (in 7,735 RCTs), “all but significant” (7,015), “a nonsignificant trend” (3,442), “failed to reach statistical significance” (2,578), and “a strong trend” (1,700). The strongest evidence for an increased temporal prevalence was found for “a numerical trend,” “a positive trend,” “an increasing trend,” and “nominally significant.” In contrast, the phrases “all but significant,” “approaches statistical significance,” “did not quite reach statistical significance,” “difference was apparent,” “failed to reach statistical significance,” and “not quite significant” decreased over time. In a random sampled subset of 29,000 phrases, the manually identified and corresponding 11,926 P values, 68,1% ranged between 0.05 and 0.15 (CI: 67. to 69.0; median 0.06). Our results show that RCT reports regularly contain specific phrases describing marginally nonsignificant results to report P values close to but above the dominant 0.05 cutoff. The fact that the prevalence of the phrases remained stable over time indicates that this practice of broadly interpreting P values close to a predefined threshold remains prevalent. To enhance responsible and transparent interpretation of RCT results, researchers, clinicians, reviewers, and editors may reduce the focus on formal statistical significance thresholds and stimulate reporting of P values with corresponding effect sizes and CIs and focus on the clinical relevance of the statistical difference found in RCTs.

The power of language to modify the reader’s perception of interpreting biomedical results cannot be underestimated. An analysis of more than half a million randomized controlled trials reveals that researchers are using appealing phrases to describe non-significant findings as if they were below the p=0.05 significance threshold.     

Bacterial Diversity in a Mine Water Treatment Plant

We investigated the microbial community in a pilot plant for treatment of acid mine water by biological ferrous iron oxidation using clone library analysis and calculated statistical parameters for further characterization. The microbial community in the plant was conspicuously dominated by a group of Betaproteobacteria affiliated with “Ferribacter polymyxa”.     

Computational Screening of Tip and Stalk Cell Behavior Proposes a Role for Apelin Signaling in Sprout Progression

Angiogenesis involves the formation of new blood vessels by sprouting or splitting of existing blood vessels. During sprouting, a highly motile type of endothelial cell, called the tip cell, migrates from the blood vessels followed by stalk cells, an endothelial cell type that forms the body of the sprout. To get more insight into how tip cells contribute to angiogenesis, we extended an existing computational model of vascular network formation based on the cellular Potts model with tip and stalk differentiation, without making a priori assumptions about the differences between tip cells and stalk cells. To predict potential differences, we looked for parameter values that make tip cells (a) move to the sprout tip, and (b) change the morphology of the angiogenic networks. The screening predicted that if tip cells respond less effectively to an endothelial chemoattractant than stalk cells, they move to the tips of the sprouts, which impacts the morphology of the networks. A comparison of this model prediction with genes expressed differentially in tip and stalk cells revealed that the endothelial chemoattractant Apelin and its receptor APJ may match the model prediction. To test the model prediction we inhibited Apelin signaling in our model and in an in vitro model of angiogenic sprouting, and found that in both cases inhibition of Apelin or of its receptor APJ reduces sprouting. Based on the prediction of the computational model, we propose that the differential expression of Apelin and APJ yields a “self-generated” gradient mechanisms that accelerates the extension of the sprout.     

UniFrac: a new phylogenetic method for comparing microbial communities

We introduce here a new method for computing differences between microbial communities based on phylogenetic information. This method, UniFrac, measures the phylogenetic distance between sets of taxa in a phylogenetic tree as the fraction of the branch length of the tree that leads to descendants from either one environment or the other, but not both. UniFrac can be used to determine whether communities are significantly different, to compare many communities simultaneously using clustering and ordination techniques, and to measure the relative contributions of different factors, such as chemistry and geography, to similarities between samples. We demonstrate the utility of UniFrac by applying it to published 16S rRNA gene libraries from cultured isolates and environmental clones of bacteria in marine sediment, water, and ice. Our results reveal that (i) cultured isolates from ice, water, and sediment resemble each other and environmental clone sequences from sea ice, but not environmental clone sequences from sediment and water; (ii) the geographical location does not correlate strongly with bacterial community differences in ice and sediment from the Arctic and Antarctic; and (iii) bacterial communities differ between terrestrially impacted seawater (whether polar or temperate) and warm oligotrophic seawater, whereas those in individual seawater samples are not more similar to each other than to those in sediment or ice samples. These results illustrate that UniFrac provides a new way of characterizing microbial communities, using the wealth of environmental rRNA sequences, and allows quantitative insight into the factors that underlie the distribution of lineages among environments.     

UniFrac: a New Phylogenetic Method for Comparing Microbial Communities

We introduce here a new method for computing differences between microbial communities based on phylogenetic information. This method, UniFrac, measures the phylogenetic distance between sets of taxa in a phylogenetic tree as the fraction of the branch length of the tree that leads to descendants from either one environment or the other, but not both. UniFrac can be used to determine whether communities are significantly different, to compare many communities simultaneously using clustering and ordination techniques, and to measure the relative contributions of different factors, such as chemistry and geography, to similarities between samples. We demonstrate the utility of UniFrac by applying it to published 16S rRNA gene libraries from cultured isolates and environmental clones of bacteria in marine sediment, water, and ice. Our results reveal that (i) cultured isolates from ice, water, and sediment resemble each other and environmental clone sequences from sea ice, but not environmental clone sequences from sediment and water; (ii) the geographical location does not correlate strongly with bacterial community differences in ice and sediment from the Arctic and Antarctic; and (iii) bacterial communities differ between terrestrially impacted seawater (whether polar or temperate) and warm oligotrophic seawater, whereas those in individual seawater samples are not more similar to each other than to those in sediment or ice samples. These results illustrate that UniFrac provides a new way of characterizing microbial communities, using the wealth of environmental rRNA sequences, and allows quantitative insight into the factors that underlie the distribution of lineages among environments.     

Non‐zero‐sum neutrality test for the tropical rain forest community using long‐term between‐census data

For community ecologists, “neutral or not?” is a fundamental question, and thus, rejecting neutrality is an important first step before investigating the deterministic processes underlying community dynamics. Hubbell''s neutral model is an important contribution to the exploration of community dynamics, both technically and philosophically. However, the neutrality tests for this model are limited by a lack of statistical power, partly because the zero‐sum assumption of the model is unrealistic. In this study, we developed a neutrality test for local communities that implements non‐zero‐sum community dynamics and determines the number of new species (N _sp) between observations. For the non‐zero‐sum neutrality test, the model distributed the expected N _sp, as calculated by extensive simulations, which allowed us to investigate the neutrality of the observed community by comparing the observed N _sp with distributions of the expected N _sp derived from the simulations. For this comparison, we developed a new “non‐zero‐sum N _sp test,” which we validated by running multiple neutral simulations using different parameter settings. We found that the non‐zero‐sum N _sp test rejected neutrality at a near‐significance level, which justified the validity of our approach. For an empirical test, the non‐zero‐sum N _sp test was applied to real tropical tree communities in Panama and Malaysia. The non‐zero‐sum N _sp test rejected neutrality in both communities when the observation interval was long and N _sp was large. Hence, the non‐zero‐sum N _sp test is an effective way to examine neutrality and has reasonable statistical power to reject the neutral model, especially when the observed N _sp is large. This unique and simple approach is statistically powerful, even though it only employs two temporal sequences of community data. Thus, this test can be easily applied to existing datasets. In addition, application of the test will provide significant benefits for detecting changing biodiversity under climate change and anthropogenic disturbance.     

Structural Insights into the Mechanism of Negative Regulation of Single-box High Mobility Group Proteins by the Acidic Tail Domain

The Drosophila and plant (maize) functional counterparts of the abundant vertebrate chromosomal protein HMGB1 (HMG-D and ZmHMGB1, respectively) differ from HMGB1 in having a single HMG box, as well as basic and acidic flanking regions that vary greatly in length and charge. We show that despite these variations, HMG-D and ZmHMGB1 exist in dynamic assemblies in which the basic HMG boxes and linkers associate with their intrinsically disordered, predominantly acidic, tails in a manner analogous to that observed previously for HMGB1. The DNA-binding surfaces of the boxes and linkers are occluded in “auto-inhibited” forms of the protein, which are in equilibrium with transient, more open structures that are “binding-competent.” This strongly suggests that the mechanism of auto-inhibition may be a general one. HMG-D and ZmHMGB1 differ from HMGB1 in having phosphorylation sites in their tail and linker regions. In both cases, in vitro phosphorylation of serine residues within the acidic tail stabilizes the assembled form, suggesting another level of regulation for interaction with DNA, chromatin, and other proteins that is not possible for the uniformly acidic (hence unphosphorylatable) tail of HMGB1.     

Community species diversity mediates the trade‐off between aboveground and belowground biomass for grasses and forbs in degraded alpine meadow,Tibetan Plateau

Although many empirical experiments have shown that increasing degradation results in lower aboveground biomass (AGB), our knowledge of the magnitude of belowground biomass (BGB) for individual plants is a prerequisite for accurately revealing the biomass trade‐off in degraded grasslands. Here, by linking the AGB and BGB of individual plants, species in the community, and soil properties, we explored the biomass partitioning patterns in different plant functional groups (grasses of Stipa capillacea and forbs of Anaphalis xylorhiza). Our results indicated that 81% and 60% of the biomass trade‐off variations could be explained by environmental factors affecting grasses and forbs, respectively. The change in community species diversity dominated the biomass trade‐off via either direct or indirect effects on soil properties and biomass. However, the community species diversity imparted divergent effects on the biomass trade‐off for grasses (scored at −0.72) and forbs (scored at 0.59). Our findings suggest that plant communities have evolved two contrasting strategies of biomass allocation patterns in degraded grasslands. These are the “conservative” strategy in grasses, in which plants with larger BGB trade‐off depends on gigantic roots for soil resources, and the “opportunistic” strategy in forbs, in which plants can adapt to degraded lands using high variation and optimal biomass allocation.     

Distinct Phyllosphere Bacterial Communities on Arabidopsis Wax Mutant Leaves

The phyllosphere of plants is inhabited by diverse microorganisms, however, the factors shaping their community composition are not fully elucidated. The plant cuticle represents the initial contact surface between microorganisms and the plant. We thus aimed to investigate whether mutations in the cuticular wax biosynthesis would affect the diversity of the phyllosphere microbiota. A set of four Arabidopsis thaliana eceriferum mutants (cer1, cer6, cer9, cer16) and their respective wild type (Landsberg erecta) were subjected to an outdoor growth period and analysed towards this purpose. The chemical distinctness of the mutant wax phenotypes was confirmed by gas chromatographic measurements. Next generation amplicon pyrosequencing of the bacterial communities showed distinct community patterns. This observation was supported by denaturing gradient gel electrophoresis experiments. Microbial community analyses revealed bacterial phylotypes that were ubiquitously present on all plant lines (termed “core” community) while others were positively or negatively affected by the wax mutant phenotype (termed “plant line-specific“ community). We conclude from this study that plant cuticular wax composition can affect the community composition of phyllosphere bacteria.     

Faecal Microbiota of Forage-Fed Horses in New Zealand and the Population Dynamics of Microbial Communities following Dietary Change

The effects of abrupt dietary transition on the faecal microbiota of forage-fed horses over a 3-week period were investigated. Yearling Thoroughbred fillies reared as a cohort were exclusively fed on either an ensiled conserved forage-grain diet (“Group A”; n = 6) or pasture (“Group B”; n = 6) for three weeks prior to the study. After the Day 0 faecal samples were collected, horses of Group A were abruptly transitioned to pasture. Both groups continued to graze similar pasture for three weeks, with faecal samples collected at 4-day intervals. DNA was isolated from the faeces and microbial 16S and 18S rRNA gene amplicons were generated and analysed by pyrosequencing. The faecal bacterial communities of both groups of horses were highly diverse (Simpson’s index of diversity >0.8), with differences between the two groups on Day 0 (P<0.017 adjusted for multiple comparisons). There were differences between Groups A and B in the relative abundances of four genera, BF311 (family Bacteroidaceae; P = 0.003), CF231 (family Paraprevotellaceae; P = 0.004), and currently unclassified members within the order Clostridiales (P = 0.003) and within the family Lachnospiraceae (P = 0.006). The bacterial community of Group A horses became similar to Group B within four days of feeding on pasture, whereas the structure of the archaeal community remained constant pre- and post-dietary change. The community structure of the faecal microbiota (bacteria, archaea and ciliate protozoa) of pasture-fed horses was also identified. The initial differences observed appeared to be linked to recent dietary history, with the bacterial community of the forage-fed horses responding rapidly to abrupt dietary change.     

Statistical Analysis of Individual Participant Data Meta-Analyses: A Comparison of Methods and Recommendations for Practice

Background

Individual participant data (IPD) meta-analyses that obtain “raw” data from studies rather than summary data typically adopt a “two-stage” approach to analysis whereby IPD within trials generate summary measures, which are combined using standard meta-analytical methods. Recently, a range of “one-stage” approaches which combine all individual participant data in a single meta-analysis have been suggested as providing a more powerful and flexible approach. However, they are more complex to implement and require statistical support. This study uses a dataset to compare “two-stage” and “one-stage” models of varying complexity, to ascertain whether results obtained from the approaches differ in a clinically meaningful way.

Methods and Findings

We included data from 24 randomised controlled trials, evaluating antiplatelet agents, for the prevention of pre-eclampsia in pregnancy. We performed two-stage and one-stage IPD meta-analyses to estimate overall treatment effect and to explore potential treatment interactions whereby particular types of women and their babies might benefit differentially from receiving antiplatelets. Two-stage and one-stage approaches gave similar results, showing a benefit of using anti-platelets (Relative risk 0.90, 95% CI 0.84 to 0.97). Neither approach suggested that any particular type of women benefited more or less from antiplatelets. There were no material differences in results between different types of one-stage model.

Conclusions

For these data, two-stage and one-stage approaches to analysis produce similar results. Although one-stage models offer a flexible environment for exploring model structure and are useful where across study patterns relating to types of participant, intervention and outcome mask similar relationships within trials, the additional insights provided by their usage may not outweigh the costs of statistical support for routine application in syntheses of randomised controlled trials. Researchers considering undertaking an IPD meta-analysis should not necessarily be deterred by a perceived need for sophisticated statistical methods when combining information from large randomised trials.     

Extensive Description and Comparison of Human Supra-Gingival Microbiome in Root Caries and Health

Knowledge of the polymicrobial etiology of root caries is limited. To conduct a comprehensive research study on root caries, we utilized 454-pyrosequencing of 16S rRNA gene libraries and quantitative PCR to compare supra-gingival bacterial communities from healthy sites and carious sites of 21 patients with root caries (Patient-controls and Patient-cases) and the sites of 21 healthy individuals (Healthy-controls) from two nursing homes. Healthy-controls and Patient-cases showed no significant differences in terms of biomass, species richness, and species diversity. However, as for beta diversity based on either community membership metric (unweighted UniFrac) or community structure metric (weighted UniFrac), Healthy-controls and Patient-cases were clearly distinguished from each other, appearing more variable in the community membership and structure in root caries microbiome but relatively conserved in the health microbiome. The Patient-controls group was at an intermediate stage between Healthy-controls and Patient-cases, but was more inclined to the former. Demonstrated in both relative abundance and prevalence of species in health and root caries, Propionibacterium acidifaciens, Streptococcus mutans, Olsenella profusa, Prevotella multisaccharivorax, and Lactobacillus crispatus were found to be most associated with root caries, whereas Delftia acidovorans, Bacteroidetes[G-2] sp., Lachnospiraceae[G-3] sp., and Prevotella intermedia are most associated with health. Our study provides a basis for further elucidating the microbial etiology of root caries in the elderly.     

Patient Knowledge on Malaria Symptoms Is a Key to Promoting Universal Access of Patients to Effective Malaria Treatment in Palawan,the Philippines

Introduction

Palawan, where health care facilities are still limited, is one of the most malaria endemic provinces in the Philippines. Since 1999, microscopists (community health workers) have been trained in malaria diagnosis and feasibility of early diagnosis and treatments have been enhanced throughout the province. To accelerate the universal access of malaria patients to diagnostic testing in Palawan, positive health seeking behavior should be encouraged when malaria infection is suspected.

Methods

In this cross-sectional study, structured interviews were carried out with residents (N = 218) of 20 remote malaria-endemic villages throughout Palawan with a history of suspected malaria from January to February in 2012. Structural equation modeling (SEM) was conducted to determine factors associated with appropriate treatment, which included: (1) socio-demographic characteristics; (2) proximity to a health facility; (3) health seeking behavior; (4) knowledge on malaria; (5) participation in community awareness-raising activities.

Results

Three factors independently associated with appropriate treatment were identified by SEM (CMIN = 10.5, df = 11, CFI = 1.000, RMSEA = .000): “living near microscopist” (p < 0.001), “not living near private pharmacy” (p < 0.01), and “having severe symptoms” (p < 0.01). “Severe symptoms” were positively correlated with more “knowledge on malaria symptoms” (p < 0.001). This knowledge was significantly increased by attending “community awareness-raising activities by microscopists” (p < 0.001).

Conclusions

In the resource-limited settings, microscopists played a significant role in providing appropriate treatment to all participants with severe malaria symptoms. However, it was considered that knowledge on malaria symptoms made participants more aware of their symptoms, and further progressed self-triage. Strengthening this recognition sensitivity and making residents aware of nearby microscopists may be the keys to accelerating universal access to effective malaria treatment in Palawan.     

Wild Sex in Zebrafish: Loss of the Natural Sex Determinant in Domesticated Strains

Sex determination can be robustly genetic, strongly environmental, or genetic subject to environmental perturbation. The genetic basis of sex determination is unknown for zebrafish (Danio rerio), a model for development and human health. We used RAD-tag population genomics to identify sex-linked polymorphisms. After verifying this “RAD-sex” method on medaka (Oryzias latipes), we studied two domesticated zebrafish strains (AB and TU), two natural laboratory strains (WIK and EKW), and two recent isolates from nature (NA and CB). All four natural strains had a single sex-linked region at the right tip of chromosome 4, enabling sex genotyping by PCR. Genotypes for the single nucleotide polymorphism (SNP) with the strongest statistical association to sex suggested that wild zebrafish have WZ/ZZ sex chromosomes. In natural strains, “male genotypes” became males and some “female genotypes” also became males, suggesting that the environment or genetic background can cause female-to-male sex reversal. Surprisingly, TU and AB lacked detectable sex-linked loci. Phylogenomics rooted on D. nigrofasciatus verified that all strains are monophyletic. Because AB and TU branched as a monophyletic clade, we could not rule out shared loss of the wild sex locus in a common ancestor despite their independent domestication. Mitochondrial DNA sequences showed that investigated strains represent only one of the three identified zebrafish haplogroups. Results suggest that zebrafish in nature possess a WZ/ZZ sex-determination mechanism with a major determinant lying near the right telomere of chromosome 4 that was modified during domestication. Strains providing the zebrafish reference genome lack key components of the natural sex-determination system but may have evolved variant sex-determining mechanisms during two decades in laboratory culture.