Means and covariance functions for geostatistical compositional data: an axiomatic approach

This work focuses on the characterization of the central tendency of a sample of compositional data. It provides new results about theoretical properties of means and covariance functions for compositional data, with an axiomatic perspective. Original results that shed new light on the geostatistical modeling of compositional data are presented. As a first result, it is shown that the weighted arithmetic mean is the only central tendency characteristic satisfying a small set of axioms, namely continuity, reflexivity and marginal stability. Moreover, this set of axioms also implies that the weights must be identical for all parts of the composition. This result has deep consequences on the spatial multivariate covariance modeling of compositional data. In a geostatistical setting, it is shown as a second result that the proportional model of covariance functions (i.e., the product of a covariance matrix and a single correlation function) is the only model that provides identical kriging weights for all components of the compositional data. As a consequence of these two results, the proportional model of covariance function is the only covariance model compatible with reflexivity and marginal stability

Radiative Transfer for Exoplanet Atmospheres

Remote sensing of the atmospheres of distant worlds motivates a firm understanding of radiative transfer. In this review, we provide a pedagogical cookbook that describes the principal ingredients needed to perform a radiative transfer calculation and predict the spectrum of an exoplanet atmosphere, including solving the radiative transfer equation, calculating opacities (and chemistry), iterating for radiative equilibrium (or not), and adapting the output of the calculations to the astronomical observations. A review of the state of the art is performed, focusing on selected milestone papers. Outstanding issues, including the need to understand aerosols or clouds and elucidating the assumptions and caveats behind inversion methods, are discussed. A checklist is provided to assist referees/reviewers in their scrutiny of works involving radiative transfer. A table summarizing the methodology employed by past studies is provided.Comment: 7 pages, no figures, 1 table. Filled in missing information in references, main text unchange

A theoretical account of cue averaging in the rodent head direction system

Head direction (HD) cell responses are thought to be derived from a combination of internal (or idiothetic) and external (or allothetic) sources of information. Recent work from the Jeffery laboratory shows that the relative influence of visual versus vestibular inputs upon the HD cell response depends on the disparity between these sources. In this paper, we present simulation results from a model designed to explain these observations. The model accurately replicates the Knight et al. data. We suggest that cue conflict resolution is critically dependent on plastic remapping of visual information onto the HD cell layer. This remap results in a shift in preferred directions of a subset of HD cells, which is then inherited by the rest of the cells during path integration. Thus, we demonstrate how, over a period of several minutes, a visual landmark may gain cue control. Furthermore, simulation results show that weaker visual landmarks fail to gain cue control as readily. We therefore suggest a second longer term plasticity in visual projections onto HD cell areas, through which landmarks with an inconsistent relationship to idiothetic information are made less salient, significantly hindering their ability to gain cue control. Our results provide a mechanism for reliability-weighted cue averaging that may pertain to other neural systems in addition to the HD system

Serological Markers Suggest Heterogeneity of Effectiveness of Malaria Control Interventions on Bioko Island, Equatorial Guinea

BACKGROUND: In order to control and eliminate malaria, areas of on-going transmission need to be identified and targeted for malaria control interventions. Immediately following intense interventions, malaria transmission can become more heterogeneous if interventions are more successful in some areas than others. Bioko Island, Equatorial Guinea, has been subject to comprehensive malaria control interventions since 2004. This has resulted in substantial reductions in the parasite burden, although this drop has not been uniform across the island. METHODS/PRINCIPAL FINDINGS: In 2008, filter paper blood samples were collected from 7387 people in a cross-sectional study incorporating 18 sentinel sites across Bioko, Equatorial Guinea. Antibodies were measured to P. falciparum Apical Membrane Antigen-1 (AMA-1) by Enzyme Linked Immunosorbent Assay (ELISA). Age-specific seropositivity rates were used to estimate seroconversion rates (SCR). Analysis indicated there had been at least a 60% decline in SCR in four out of five regions on the island. Changes in SCR showed a high degree of congruence with changes in parasite rate (PR) and with regional reductions in all cause child mortality. The mean age adjusted concentration of anti-AMA-1 antibodies was mapped to identify areas where individual antibody responses were higher than expected. This approach confirmed the North West of the island as a major focus of continuing infection and an area where control interventions need to be concentrated or re-evaluated. CONCLUSION/INTERPRETATION: Both SCR and PR revealed heterogeneity in malaria transmission and demonstrated the variable effectiveness of malaria control measures. This work confirms the utility of serological analysis as an adjunct measure for monitoring transmission. Age-specific seroprevalence based evidence of changes in transmission over time will be of particular value when no baseline data are available. Importantly, SCR data provide additional evidence to link malaria control activities to contemporaneous reductions in all-cause child mortality

Optimizing topological cascade resilience based on the structure of terrorist networks

Complex socioeconomic networks such as information, finance and even terrorist networks need resilience to cascades - to prevent the failure of a single node from causing a far-reaching domino effect. We show that terrorist and guerrilla networks are uniquely cascade-resilient while maintaining high efficiency, but they become more vulnerable beyond a certain threshold. We also introduce an optimization method for constructing networks with high passive cascade resilience. The optimal networks are found to be based on cells, where each cell has a star topology. Counterintuitively, we find that there are conditions where networks should not be modified to stop cascades because doing so would come at a disproportionate loss of efficiency. Implementation of these findings can lead to more cascade-resilient networks in many diverse areas.Comment: 26 pages. v2: In review at Public Library of Science ON

Accurate path integration in continuous attractor network models of grid cells

Grid cells in the rat entorhinal cortex display strikingly regular firing responses to the animal's position in 2-D space and have been hypothesized to form the neural substrate for dead-reckoning. However, errors accumulate rapidly when velocity inputs are integrated in existing models of grid cell activity. To produce grid-cell-like responses, these models would require frequent resets triggered by external sensory cues. Such inadequacies, shared by various models, cast doubt on the dead-reckoning potential of the grid cell system. Here we focus on the question of accurate path integration, specifically in continuous attractor models of grid cell activity. We show, in contrast to previous models, that continuous attractor models can generate regular triangular grid responses, based on inputs that encode only the rat's velocity and heading direction. We consider the role of the network boundary in the integration performance of the network and show that both periodic and aperiodic networks are capable of accurate path integration, despite important differences in their attractor manifolds. We quantify the rate at which errors in the velocity integration accumulate as a function of network size and intrinsic noise within the network. With a plausible range of parameters and the inclusion of spike variability, our model networks can accurately integrate velocity inputs over a maximum of ~10–100 meters and ~1–10 minutes. These findings form a proof-of-concept that continuous attractor dynamics may underlie velocity integration in the dorsolateral medial entorhinal cortex. The simulations also generate pertinent upper bounds on the accuracy of integration that may be achieved by continuous attractor dynamics in the grid cell network. We suggest experiments to test the continuous attractor model and differentiate it from models in which single cells establish their responses independently of each other

Structural variation in the chicken genome identified by paired-end next-generation DNA sequencing of reduced representation libraries

<p>Abstract</p> <p>Background</p> <p>Variation within individual genomes ranges from single nucleotide polymorphisms (SNPs) to kilobase, and even megabase, sized structural variants (SVs), such as deletions, insertions, inversions, and more complex rearrangements. Although much is known about the extent of SVs in humans and mice, species in which they exert significant effects on phenotypes, very little is known about the extent of SVs in the 2.5-times smaller and less repetitive genome of the chicken.</p> <p>Results</p> <p>We identified hundreds of shared and divergent SVs in four commercial chicken lines relative to the reference chicken genome. The majority of SVs were found in intronic and intergenic regions, and we also found SVs in the coding regions. To identify the SVs, we combined high-throughput short read paired-end sequencing of genomic reduced representation libraries (RRLs) of pooled samples from 25 individuals and computational mapping of DNA sequences from a reference genome.</p> <p>Conclusion</p> <p>We provide a first glimpse of the high abundance of small structural genomic variations in the chicken. Extrapolating our results, we estimate that there are thousands of rearrangements in the chicken genome, the majority of which are located in non-coding regions. We observed that structural variation contributes to genetic differentiation among current domesticated chicken breeds and the Red Jungle Fowl. We expect that, because of their high abundance, SVs might explain phenotypic differences and play a role in the evolution of the chicken genome. Finally, our study exemplifies an efficient and cost-effective approach for identifying structural variation in sequenced genomes.</p

Obesity resistant mechanisms in the Lean polygenic mouse model as indicated by liver transcriptome and expression of selected genes in skeletal muscle

<p>Abstract</p> <p>Background</p> <p>Divergently selected Lean and Fat mouse lines represent unique models for a polygenic form of resistance and susceptibility to obesity development. Previous research on these lines focused mainly on obesity-susceptible factors in the Fat line. This study aimed to examine the molecular basis of obesity-resistant mechanisms in the Lean line by analyzing various fat depots and organs, the liver transcriptome of selected metabolic pathways, plasma and lipid homeostasis and expression of selected skeletal muscle genes.</p> <p>Results</p> <p>Expression profiling using our custom Steroltalk v2 microarray demonstrated that Lean mice exhibit a higher hepatic expression of cholesterol biosynthesis genes compared to the Fat line, although this was not reflected in elevation of total plasma or liver cholesterol. However, FPLC analysis showed that protective HDL cholesterol was elevated in Lean mice. A significant difference between the strains was also found in bile acid metabolism. Lean mice had a higher expression of <it>Cyp8b1</it>, a regulatory enzyme of bile acid synthesis, and the <it>Abcb11 </it>bile acid transporter gene responsible for export of acids to the bile. Additionally, a higher content of blood circulating bile acids was observed in Lean mice. Elevated HDL and upregulation of some bile acids synthesis and transport genes suggests enhanced reverse cholesterol transport in the Lean line - the flux of cholesterol out of the body is higher which is compensated by upregulation of endogenous cholesterol biosynthesis. Increased skeletal muscle <it>Il6 </it>and <it>Dio2 </it>mRNA levels as well as increased activity of muscle succinic acid dehydrogenase (SDH) in the Lean mice demonstrates for the first time that changes in muscle energy metabolism play important role in the Lean line phenotype determination and corroborate our previous findings of increased physical activity and thermogenesis in this line. Finally, differential expression of <it>Abcb11 </it>and <it>Dio2 </it>identifies novel strong positional candidate genes as they map within the quantitative trait loci (QTL) regions detected previously in crosses between the Lean and Fat mice.</p> <p>Conclusion</p> <p>We identified novel candidate molecular targets and metabolic changes which can at least in part explain resistance to obesity development in the Lean line. The major difference between the Lean and Fat mice was in increased liver cholesterol biosynthesis gene mRNA expression, bile acid metabolism and changes in selected muscle genes' expression in the Lean line. The liver <it>Abcb11 </it>and muscle <it>Dio2 </it>were identified as novel positional candidate genes to explain part of the phenotypic difference between the Lean and Fat lines.</p

Genetic Analysis of HIV-1 Subtypes in Nairobi, Kenya

Background: Genetic analysis of a viral infection helps in following its spread in a given population, in tracking the routes of infection and, where applicable, in vaccine design. Additionally, sequence analysis of the viral genome provides information about patterns of genetic divergence that may have occurred during viral evolution. Objective: In this study we have analyzed the subtypes of Human Immunodeficiency Virus -1 (HIV-1) circulating in a diverse sample population of Nairobi, Kenya. Methodology: 69 blood samples were collected from a diverse subject population attending the Aga Khan University Hospital in Nairobi, Kenya. Total DNA was extracted from peripheral blood mononuclear cells (PBMCs), and used in a Polymerase Chain Reaction (PCR) to amplify the HIV gag gene. The PCR amplimers were partially sequenced, and alignment and phylogenetic analysis of these sequences was performed using the Los Alamos HIV Database. Results: Blood samples from 69 HIV-1 infected subjects from varying ethnic backgrounds were analyzed. Sequence alignment and phylogenetic analysis showed 39 isolates to be subtype A, 13 subtype D, 7 subtype C, 3 subtype AD and CRF01_AE, 2 subtype G and 1 subtype AC and 1 AG. Deeper phylogenetic analysis revealed HIV subtype A sequences to be highly divergent as compared to subtypes D and C. Conclusion: Our analysis indicates that HIV-1 subtypes in the Nairobi province of Kenya are dominated by a genetically diverse clade A. Additionally, the prevalence of highly divergent, complex subtypes, intersubtypes, and the recombinant forms indicates viral mixing in Kenyan population, possibly as a result of dual infections