ALBERT

Description: Background: Algorithms designed to detect complex genetic disease associations are initially evaluated using simulated datasets. Typical evaluations vary constraints that influence the correct detection of underlying models (i.e. number of loci, heritability, and minor allele frequency). Such studies neglect to account for model architecture (i.e. the unique specification and arrangement of penetrance values comprising the genetic model), which alone can influence the detectability of a model. In order to design a simulation study which efficiently takes architecture into account, areliable metric is needed for model selection. Results: We evaluate three metrics as predictors of relative model detection difficulty derived from previous works: (1) Penetrance table variance (PTV), (2) customized odds ratio (COR), and (3) our own Ease of Detection Measure (EDM), calculated from the penetrance values and respective genotype frequencies of each simulated genetic model. We evaluate the reliability of these metrics across three very different data search algorithms, each with the capacity to detect epistatic interactions. We find that a model's EDM and COR are each stronger predictors of model detection success than heritability. Conclusions: This study formally identifies and evaluates metrics which quantify model detection difficulty. We utilize these metrics to intelligently select models from a population of potential architectures. This allows for an improved simulation study design which accounts for differences in detection difficulty attributed to model architecture. We implement the calculation and utilization of EDM and COR into GAMETES, an algorithm which rapidly and precisely generates pure, strict, n-locus epistatic models.

Description: Background: Geneticists who look beyond single locus disease associations require additional strategies for the detection of complex multi-locus effects. Epistasis, a multi-locus masking effect, presents a particular challenge, and has been the target of bioinformatic development. Thorough evaluation of new algorithms calls for simulation studies in which known disease models are sought. To date, the best methods for generating simulated multi-locus epistatic models rely on genetic algorithms. However, such methods are computationally expensive, difficult to adapt to multiple objectives, and unlikely to yield models with a precise form of epistasis which we refer to as pure and strict. Purely and strictly epistatic models constitute the worst-case in terms of detecting disease associations, since such associations may only be observed if all n-loci are included in the disease model. This makes them an attractive gold standard for simulation studies considering complex multi-locus effects. Results: We introduce GAMETES, a user-friendly software package and algorithm which generates complex biallelic single nucleotide polymorphism (SNP) disease models for simulation studies. GAMETES rapidly and precisely generates random, pure, strict n-locus models with specified genetic constraints. These constraints include heritability, minor allele frequencies of the SNPs, and population prevalence. GAMETES also includes a simple dataset simulation strategy which may be utilized to rapidly generate an archive of simulated datasets for given genetic models. We highlight the utility and limitations of GAMETES with an example simulation study using MDR, an algorithm designed to detect epistasis. Conclusions: GAMETES is a fast, flexible, and precise tool for generating complex n-locus models with random architectures. While GAMETES has a limited ability to generate models with higher heritabilities, it is proficient at generating the lower heritability models typically used in simulation studies evaluating new algorithms. In addition, the GAMETES modeling strategy may be flexibly combined with any dataset simulation strategy. Beyond dataset simulation, GAMETES could be employed to pursue theoretical characterization of genetic models and epistasis.

Description: Background: Glutamyl queuosine-tRNAAsp synthetase (GluQ-RS) is a paralog of the catalytic domain of glutamyl-tRNA synthetase and catalyzes the formation of glutamyl-queuosine on the wobble position of tRNAAsp. Here we analyze the transcription of its gene in Shigella flexneri, where it is found downstream of dksA, which encodes a transcriptional regulator involved in stress responses. Results: The genomic organization, dksA-gluQ-rs, is conserved in more than 40 bacterial species. RT-PCR assays show co-transcription of both genes without a significant change in transcript levels during growth of S. flexneri. However, mRNA levels of the intergenic region changed during growth, increasing at stationary phase, indicating an additional level of control over the expression of gluQ-rs gene. Transcriptional fusions with lacZ as a reporter gene only produced beta-galactosidase activity when the constructs included the dksA promoter, indicating that gluQ-rs do not have a separate promoter. Using bioinformatics, we identified a putative transcriptional terminator between dksA and gluQ-rs. Deletion or alteration of the predicted terminator resulted in increased expression of the lacZ reporter compared with cells containing the wild type terminator sequence. Analysis of the phenotype of a gluQ-rs mutant suggested that it may play a role in some stress responses, since growth of the mutant was impaired in the presence of osmolytes. Conclusions: The results presented here, show that the expression of gluQ-rs depends on the dksA promoter, and strongly suggest the presence and the functionality of a transcriptional terminator regulating its expression. Also, the results indicate a link between glutamyl-queuosine synthesis and stress response in Shigella flexneri.

Description: Background: Inflammation and increased ceramide concentrations characterise adipose tissue of obese women with high liver fat content compared to equally obese women with normal liver fat content. The present study characterises enzymes involved in ceramide metabolism in subcutaneous and intra-abdominal adipose tissue. Methods: Pathways leading to increased ceramide concentrations in inflamed versus non-inflamed adipose tissue were investigated by quantifying expression levels of key enzymes involved in ceramide metabolism. Sphingomyelinases (sphingomyelin phosphodiesterases SMPD1-3) were investigated further using immunohistochemistry to establish their location within adipose tissue, and their mRNA expression levels were determined in subcutaneous and intra-abdominal adipose tissue from both non-obese and obese subject. Results: Gene expression levels of sphingomyelinases, enzymes that hydrolyze sphingomyelin to ceramide, rather than enzymes involved in de novo ceramide synthesis, were higher in inflamed compared to non-inflamed adipose tissue of obese women (with high and normal liver fat contents respectively). Sphingomyelinases were localised to both macrophages and adipocytes, but also to blood vessels and to extracellular regions surrounding vessels within adipose tissue. Expression levels of SMPD3 mRNA correlated significantly with concentrations of different ceramides and sphingomyelins. In both non-obese and obese subjects SMPD3 mRNA levels were higher in the more inflamed intra-abdominal compared to the subcutaneous adipose tissue depot. Conclusions: Generation of ceramides within adipose tissue as a result of sphingomyelinase action may contribute to inflammation in human adipose tissue.

Description: Background: A sedentary lifestyle predisposes to cardiometabolic diseases. Lifestyle changes such as increased physical activity improve a range of cardiometabolic risk factors. The objective of this study was to examine whether functional changes in adipose tissue were related to these improvements. Methods: Seventy-three sedentary, overweight (mean BMI 29.9 +/- 3.2 kg/m2) and abdominally obese, but otherwise healthy men and women (67.6 +/- 0.5 years) from a randomised controlled trial of physical activity on prescription over a 6-month period were included (control n = 43, intervention n = 30). Detailed examinations were carried out at baseline and at follow-up, including fasting blood samples, a comprehensive questionnaire and subcutaneous adipose tissue biopsies for fatty acid composition analysis (n = 73) and quantification of mRNA expression levels of 13 candidate genes (n = 51), including adiponectin, leptin and inflammatory cytokines. Results: At follow-up, the intervention group had a greater increase in exercise time (+137 min/week) and a greater decrease in body fat mass (1.5 kg) compared to the control subjects (changes of 0 min/week and 0.5 kg respectively). Circulating concentrations of adiponectin were unchanged, but those of leptin decreased significantly more in the intervention group (1.8 vs 1.1 ng/mL for intervention vs control, P 〈 0.05). The w6-polyunsaturated fatty acid content, in particular linoleic acid (18:2w6), of adipose tissue increased significantly more in the intervention group, but the magnitude of the change was small (+0.17 vs +0.02 percentage points for intervention vs control, P 〈 0.05). Surprisingly leptin mRNA levels in adipose tissue increased in the intervention group (+107% intervention vs 20% control, P 〈 0.05), but changes in expression of the remaining genes did not differ between the groups. Conclusions: After a 6-month period of increased physical activity in overweight elderly individuals, circulating leptin concentrations decreased despite increased levels of leptin mRNA in adipose tissue. Otherwise, only minor changes occurred in adipose tissue, although several improvements in metabolic parameters accompanied the modest increase in physical activity.

Description: Background: Burkholderia pseudomallei and Burkholderia mallei are gram-negative pathogens responsible for the diseases melioidosis and glanders, respectively. Both species cause disease in humans and animals and have been designated as category B select agents by the Centers for Disease Control and Prevention (CDC). Burkholderia thailandensis is a closely related bacterium that is generally considered avirulent for humans. While it can cause disease in rodents, the B. thailandensis 50% lethal dose (LD50) is typically¿¿¿104-fold higher than the B. pseudomallei and B. mallei LD50in mammalian models of infection. Here we describe an alternative to mammalian hosts in the study of virulence and host-pathogen interactions of these Burkholderia species. Results: Madagascar hissing cockroaches (MH cockroaches) possess a number of qualities that make them desirable for use as a surrogate host, including ease of breeding, ease of handling, a competent innate immune system, and the ability to survive at 37°C. MH cockroaches were highly susceptible to infection with B. pseudomallei, B. mallei and B. thailandensis and the LD50 was 105 cfu. B. pseudomallei, B. mallei, and B. thailandensis cluster 1 type VI secretion system (T6SS-1) mutants were all attenuated in MH cockroaches, which is consistent with previous virulence studies conducted in rodents. B. pseudomallei mutants deficient in the other five T6SS gene clusters, T6SS-2 through T6SS-6, were virulent in both MH cockroaches and hamsters. Hemocytes obtained from MH cockroaches infected with B. pseudomallei harbored numerous intracellular bacteria, suggesting that this facultative intracellular pathogen can survive and replicate inside of MH cockroach phagocytic cells. The hemolymph extracted from these MH cockroaches also contained multinuclear giant cells (MNGCs) with intracellular B. pseudomallei, which indicates that infected hemocytes can fuse while flowing through the insect¿s open circulatory system in vivo. Conclusions: The results demonstrate that MH cockroaches are an attractive alternative to mammals to study host-pathogen interactions and may allow the identification of new Burkholderia virulence determinants. The importance of T6SS-1 as a virulence factor in MH cockroaches and rodents suggests that the primary role of this secretion system is to target evasion of the innate immune system.

Description: Current models of landscape response to Holocene climate change in midcontinent North America largely reconcile Earth orbital and atmospheric climate forcing with pollen-based forest histories on the east and eolian chronologies in Great Plains grasslands on the west. However, thousands of sand dunes spread across 12,000 km2 in eastern upper Michigan (EUM), more than 500 km east of the present forest-prairie ecotone, present a challenge to such models. We use 65 optically stimulated luminescence (OSL) ages on quartz sand deposited in silt caps (n = 8) and dunes (n = 57) to document eolian activity in EUM. Dune building was widespread ca. 10–8 ka, indicating a sharp, sustained decline in forest cover during that period. This decline was roughly coincident with hydrologic closure of the upper Great Lakes, but temporally inconsistent with most pollen-based models that imply canopy closure throughout the Holocene. Early Holocene forest openings are rarely recognized in pollen sums from EUM because faint signatures of non-arboreal pollen are largely obscured by abundant and highly mobile pine pollen. Early Holocene spikes in nonarboreal pollen are recorded in cores from small ponds, but suggest only a modest extent of forest openings. OSL dating of dune emplacement provides a direct, spatially explicit archive of greatly diminished forest cover during a very dry climate in eastern midcontinent North America ca. 10–8 ka.

Description: The abundant iron sulfide mineral pyrite has been shown to catalytically produce hydrogen peroxide (H2O2) and hydroxyl radical (.OH) in slurries of oxygenated water. Understanding the formation and fate of these reactive oxygen species is important to biological and ecological systems as exposure can lead to deleterious health effects, but also in environmental engineering during the optimization of remediation approaches for possible treatment of contaminated waste streams. This study presents the use of the amino acid phenylalanine (Phe) to monitor the kinetics of pyrite-induced .OH formation through rates of hydroxylation forming three isomers of tyrosine (Tyr) - ortho-, meta-, and para-Tyr. Results indicate that about 50% of the Phe loss results in Tyr formation, and that these products further react with .OH at rates comparable to Phe. The loss of Phe appeared to be pseudo first-order in [Phe] as a function of time, but for the first time showed that rates were much less than first-order as a function of initial substrate concentration, [Phe]. These results can be rationalized by considering that the effective concentration of .OH in solution is lower at higher level of reactant and that an increasing fraction of .OH is consumed by Phe-degradation products as a function of time. A simplified first-order model was created to describe Phe loss in pyrite slurries which incorporates the initial [Phe], a first-order dependence on pyrite surface area, the assumption that all Phe degradation products compete equally for the limited supply of highly reactive .OH, and a flux that is related to the release of H2O2 from the pyrite surface through the Fenton reaction. An empirically derived rate constant, Kpyr, was introduced to describe a variable .OH-reactivity for different batches of pyrite. Both the simplified first-order kinetic model, and a more detailed numerical simulation, yielded results that compare well to the observed kinetic data describing the effects of variations in concentrations of both initial [Phe] and pyrite. This work supports the use of Phe as a useful probe to assess the formation of .OH in the presence of pyrite, and its possible utility for similar applications with other minerals.

Description: Fault rocks can function as barriers to subsurface fluid flow and affect the storage of CO2 in geological structures. Even though flow across faults often involves more than one fluid phase, it is typically modeled using only single-phase functions due to a lack of fault rock relative permeability data and complexities in incorporating two-phase flow properties into flow simulations. Here we present two-phase fluid flow data for cataclastic fault rocks in porous sandstone from the 90-Fathom fault (northeast England). The study area represents a field analogue for North Sea saline aquifers of Permian–Triassic age that are currently being considered for CO2 storage. We use the experimental data to populate a synthetic model of a faulted saline aquifer to assess the impact of these fault rocks on CO2 injection. We show that even fault rocks with low clay contents and very limited quartz cementation can act as major baffles to the flow of a non-wetting phase if realistic two-phase properties are taken into account. Consequently, pressure may increase far more rapidly in the storage compartment during CO2 injection than anticipated based on models that only incorporate absolute fault rock permeabilities. To avoid high pressures, which may lead to hydrofracturing and CO2 leakage, either more complex injection strategies need to be adopted or seismic data acquired to ensure the absence of faults in aquifers selected for CO2 storage.

Description: Background: Army ants are dominant invertebrate predators in tropical and subtropical terrestrial ecosystems. Their close relatives within the dorylomorph group of ants are also highly specialized predators, although much less is known about their biology. We analyzed molecular data generated from 11 nuclear genes to infer a phylogeny for the major dorylomorph lineages, and incorporated fossil evidence to infer divergence times under a relaxed molecular clock. Results: Because our results indicate that one subfamily and several genera of dorylomorphs are non-monophyletic, we propose to subsume the six previous dorylomorph subfamilies into a single subfamily, Dorylinae. We find the monophyly of Dorylinae to be strongly supported and estimate the crown age of the group at 87 (74-101) million years. Our phylogenetic analyses provide only weak support for army ant monophyly and also call into question a previous hypothesis that army ants underwent a fundamental split into New World and Old World lineages. Outside the army ants, our phylogeny reveals for the first time many old, distinct lineages in the Dorylinae. The genus Cerapachys is shown to be non-monophyletic and comprised of multiple lineages scattered across the Dorylinae tree. We recover, with strong support, novel relationships among these Cerapachys-like clades and other doryline genera, but divergences in the deepest parts of the tree are not well resolved. We find the genus Sphinctomyrmex, characterized by distinctive abdominal constrictions, to consist of two separate lineages with convergent morphologies, one inhabiting the Old World and the other the New World tropics. Conclusions: While we obtain good resolution in many parts of the Dorylinae phylogeny, relationships deep in the tree remain unresolved, with major lineages joining each other in various ways depending upon the analytical method employed, but always with short internodes. This may be indicative of rapid radiation in the early history of the Dorylinae, but additional molecular data and more complete species sampling are needed for confirmation. Our phylogeny now provides a basic framework for comparative biological analyses, but much additional study on the behavior and morphology of doryline species is needed, especially investigations directed at the non-army ant taxa.