Point-estimating observer models for latent cause detection

The spatial distribution of visual items allows us to infer the presence of latent causes in the world. For instance, a spatial cluster of ants allows us to infer the presence of a common food source. However, optimal inference requires the integration of a computationally intractable number of world states in real world situations. For example, optimal inference about whether a common cause exists based on N spatially distributed visual items requires marginalizing over both the location of the latent cause and 2^N possible affiliation patterns (where each item may be affiliated or non-affiliated with the latent cause). How might the brain approximate this inference? We show that subject behaviour deviates qualitatively from Bayes-optimal, in particular showing an unexpected positive effect of N (the number of visual items) on the false-alarm rate. We propose several “point-estimating” observer models that fit subject behaviour better than the Bayesian model. They each avoid a costly computational marginalization over at least one of the variables of the generative model by “committing” to a point estimate of at least one of the two generative model variables. These findings suggest that the brain may implement partially committal variants of Bayesian models when detecting latent causes based on complex real world data.     

Convergence and constraint in eukaryotic release factor 1 (eRF1) domain 1: the evolution of stop codon specificity

Class 1 release factor in eukaryotes (eRF1) recognizes stop codons and promotes peptide release from the ribosome. The ‘molecular mimicry’ hypothesis suggests that domain 1 of eRF1 is analogous to the tRNA anticodon stem–loop. Recent studies strongly support this hypothesis and several models for specific interactions between stop codons and residues in domain 1 have been proposed. In this study we have sequenced and identified novel eRF1 sequences across a wide diversity of eukaryotes and re-evaluated the codon-binding site by bioinformatic analyses of a large eRF1 dataset. Analyses of the eRF1 structure combined with estimates of evolutionary rates at amino acid sites allow us to define the residues that are under structural (i.e. those involved in intramolecular interactions) versus non-structural selective constraints. Furthermore, we have re-assessed convergent substitutions in the ciliate variant code eRF1s using maximum likelihood-based phylogenetic approaches. Our results favor the model proposed by Bertram et al. that stop codons bind to three ‘cavities’ on the protein surface, although we suggest that the stop codon may bind in the opposite orientation to the original model. We assess the feasibility of this alternative binding orientation with a triplet stop codon and the eRF1 domain 1 structures using molecular modeling techniques.     

Moving in and out of Poverty: The Within-Individual Association between Socioeconomic Status and Juvenile Delinquency

A family’s SES can be changeable over time. This study was the first to investigate if such within-individual changes in family SES are associated with parallel fluctuations in boys’ delinquent behavior from childhood to adolescence. Participants were a community sample of boys and their caregivers (N = 503) who were assessed annually for ten consecutive years spanning ages 7–18. Fixed effects models revealed that changes in familial SES were related to changes in delinquency: Youths were more likely to offend during years in which their parents’ SES was lower than during years in which their parents’ SES was higher. Contrary to expectations, we found no evidence that this association was accounted for by families moving to different neighborhoods or by changes in parenting. Since within-individual models provide a stricter test of causality than between-individual models, these findings support claims that impacting familial SES may have a direct effect on youths’ delinquency.     

A Bayesian phase 2 model based adaptive design to optimise antivenom dosing: Application to a dose-finding trial for a novel Russell’s viper antivenom in Myanmar

For most antivenoms there is little information from clinical studies to infer the relationship between dose and efficacy or dose and toxicity. Antivenom dose-finding studies usually recruit too few patients (e.g. fewer than 20) relative to clinically significant event rates (e.g. 5%). Model based adaptive dose-finding studies make efficient use of accrued patient data by using information across dosing levels, and converge rapidly to the contextually defined ‘optimal dose’. Adequate sample sizes for adaptive dose-finding trials can be determined by simulation. We propose a model based, Bayesian phase 2 type, adaptive clinical trial design for the characterisation of optimal initial antivenom doses in contexts where both efficacy and toxicity are measured as binary endpoints. This design is illustrated in the context of dose-finding for Daboia siamensis (Eastern Russell’s viper) envenoming in Myanmar. The design formalises the optimal initial dose of antivenom as the dose closest to that giving a pre-specified desired efficacy, but resulting in less than a pre-specified maximum toxicity. For Daboia siamensis envenoming, efficacy is defined as the restoration of blood coagulability within six hours, and toxicity is defined as anaphylaxis. Comprehensive simulation studies compared the expected behaviour of the model based design to a simpler rule based design (a modified ‘3+3’ design). The model based design can identify an optimal dose after fewer patients relative to the rule based design. Open source code for the simulations is made available in order to determine adequate sample sizes for future adaptive snakebite trials. Antivenom dose-finding trials would benefit from using standard model based adaptive designs. Dose-finding trials where rare events (e.g. 5% occurrence) are of clinical importance necessitate larger sample sizes than current practice. We will apply the model based design to determine a safe and efficacious dose for a novel lyophilised antivenom to treat Daboia siamensis envenoming in Myanmar.     

Life-History Traits of the Miocene Hipparion concudense (Spain) Inferred from Bone Histological Structure

Histological analyses of fossil bones have provided clues on the growth patterns and life history traits of several extinct vertebrates that would be unavailable for classical morphological studies. We analyzed the bone histology of Hipparion to infer features of its life history traits and growth pattern. Microscope analysis of thin sections of a large sample of humeri, femora, tibiae and metapodials of Hipparion concudense from the upper Miocene site of Los Valles de Fuentidueña (Segovia, Spain) has shown that the number of growth marks is similar among the different limb bones, suggesting that equivalent skeletochronological inferences for this Hipparion population might be achieved by means of any of the elements studied. Considering their abundance, we conducted a skeletechronological study based on the large sample of third metapodials from Los Valles de Fuentidueña together with another large sample from the Upper Miocene locality of Concud (Teruel, Spain). The data obtained enabled us to distinguish four age groups in both samples and to determine that Hipparion concudense tended to reach skeletal maturity during its third year of life. Integration of bone microstructure and skeletochronological data allowed us to identify ontogenetic changes in bone structure and growth rate and to distinguish three histologic ontogenetic stages corresponding to immature, subadult and adult individuals. Data on secondary osteon density revealed an increase in bone remodeling throughout the ontogenetic stages and a lesser degree thereof in the Concud population, which indicates different biomechanical stresses in the two populations, likely due to environmental differences. Several individuals showed atypical growth patterns in the Concud sample, which may also reflect environmental differences between the two localities. Finally, classification of the specimens’ age within groups enabled us to characterize the age structure of both samples, which is typical of attritional assemblages.     

Predicting species distributions from herbarium collections: does climate bias in collection sampling influence model outcomes?

Aim Species distribution models and geographical information system (GIS) technologies are becoming increasingly important tools in conservation planning and decision‐making. Often the rich data bases of museums and herbaria serve as the primary data for predicting species distributions. Yet key assumptions about the primary data often are untested, and violation of such assumptions may have consequences for model predictions. For example, users of primary data assume that sampling has been random with respect to geography and environmental gradients. Here we evaluate the assumption that plant voucher specimens adequately sample the climatic gradient and test whether violation of this assumption influences model predictions. Location Bolivia and Ecuador. Methods Using 323,711 georeferenced herbarium collections and nine climatic variables, we predicted the distribution of 76 plant species using maximum entropy models (MAXENT) with training points that sampled the climate environments randomly and training points that reflected the climate bias in the herbarium collections. To estimate the distribution of species, MAXENT finds the distribution of maximum entropy (i.e. closest to uniform) subject to the constraint that the expected value for each environmental variable under the estimated distribution matches its empirical average. The experimental design included species that differed in geographical range and elevation; all species were modelled with 20 and 100 training points. We examined the influence of the number of training points and climate bias in training points, elevation and range size on model performance using analysis of variance models. Results We found that significant parts of the climatic gradient were poorly represented in herbarium collections for both countries. For the most part, existing climatic bias in collections did not greatly affect distribution predictions when compared with an unbiased data set. Although the effects of climate bias on prediction accuracy were found to be greater where geographical ranges were characterized by high spatial variation in the degree of climate bias (i.e. ranges where the bias of the various climates sampled by collections deviated considerably from the mean bias), the greatest influence on model performance was the number of presence points used to train the model. Main conclusions These results demonstrate that predictions of species distributions can be quite good despite existing climatic biases in primary data found in natural history collections, if a sufficiently large number of training points is available. Because of consistent overprediction of models, these results also confirm the importance of validating models with independent data or expert opinion. Failure to include independent model validation, especially in cases where training points are limited, may potentially lead to grave errors in conservation decision‐making and planning.     

Marginalized models for moderate to long series of longitudinal binary response data

Marginalized models (Heagerty, 1999, Biometrics 55, 688-698) permit likelihood-based inference when interest lies in marginal regression models for longitudinal binary response data. Two such models are the marginalized transition and marginalized latent variable models. The former captures within-subject serial dependence among repeated measurements with transition model terms while the latter assumes exchangeable or nondiminishing response dependence using random intercepts. In this article, we extend the class of marginalized models by proposing a single unifying model that describes both serial and long-range dependence. This model will be particularly useful in longitudinal analyses with a moderate to large number of repeated measurements per subject, where both serial and exchangeable forms of response correlation can be identified. We describe maximum likelihood and Bayesian approaches toward parameter estimation and inference, and we study the large sample operating characteristics under two types of dependence model misspecification. Data from the Madras Longitudinal Schizophrenia Study (Thara et al., 1994, Acta Psychiatrica Scandinavica 90, 329-336) are analyzed.     

Do Baseline P-Values Follow a Uniform Distribution in Randomised Trials?

Background

The theory has been put forward that if a null hypothesis is true, P-values should follow a Uniform distribution. This can be used to check the validity of randomisation.

Method

The theory was tested by simulation for two sample t tests for data from a Normal distribution and a Lognormal distribution, for two sample t tests which are not independent, and for chi-squared and Fisher’s exact test using small and using large samples.

Results

For the two sample t test with Normal data the distribution of P-values was very close to the Uniform. When using Lognormal data this was no longer true, and the distribution had a pronounced mode. For correlated tests, even using data from a Normal distribution, the distribution of P-values varied from simulation run to simulation run, but did not look close to Uniform in any realisation. For binary data in a small sample, only a few probabilities were possible and distribution was very uneven. With a sample of two groups of 1,000 observations, there was great unevenness in the histogram and a poor fit to the Uniform.

Conclusions

The notion that P-values for comparisons of groups using baseline data in randomised clinical trials should follow a Uniform distribution if the randomisation is valid has been found to be true only in the context of independent variables which follow a Normal distribution, not for Lognormal data, correlated variables, or binary data using either chi-squared or Fisher’s exact tests. This should not be used as a check for valid randomisation.     

Inferring Epidemiological Dynamics with Bayesian Coalescent Inference: The Merits of Deterministic and Stochastic Models

Estimation of epidemiological and population parameters from molecular sequence data has become central to the understanding of infectious disease dynamics. Various models have been proposed to infer details of the dynamics that describe epidemic progression. These include inference approaches derived from Kingman’s coalescent theory. Here, we use recently described coalescent theory for epidemic dynamics to develop stochastic and deterministic coalescent susceptible–infected–removed (SIR) tree priors. We implement these in a Bayesian phylogenetic inference framework to permit joint estimation of SIR epidemic parameters and the sample genealogy. We assess the performance of the two coalescent models and also juxtapose results obtained with a recently published birth–death-sampling model for epidemic inference. Comparisons are made by analyzing sets of genealogies simulated under precisely known epidemiological parameters. Additionally, we analyze influenza A (H1N1) sequence data sampled in the Canterbury region of New Zealand and HIV-1 sequence data obtained from known United Kingdom infection clusters. We show that both coalescent SIR models are effective at estimating epidemiological parameters from data with large fundamental reproductive number R₀ and large population size S₀. Furthermore, we find that the stochastic variant generally outperforms its deterministic counterpart in terms of error, bias, and highest posterior density coverage, particularly for smaller R₀ and S₀. However, each of these inference models is shown to have undesirable properties in certain circumstances, especially for epidemic outbreaks with R₀ close to one or with small effective susceptible populations.     

Ecological Niche Modelling Predicts Southward Expansion of Lutzomyia (Nyssomyia) flaviscutellata (Diptera: Psychodidae: Phlebotominae), Vector of Leishmania (Leishmania) amazonensis in South America,under Climate Change

Vector borne diseases are susceptible to climate change because distributions and densities of many vectors are climate driven. The Amazon region is endemic for cutaneous leishmaniasis and is predicted to be severely impacted by climate change. Recent records suggest that the distributions of Lutzomyia (Nyssomyia) flaviscutellata and the parasite it transmits, Leishmania (Leishmania) amazonensis, are expanding southward, possibly due to climate change, and sometimes associated with new human infection cases. We define the vector’s climatic niche and explore future projections under climate change scenarios. Vector occurrence records were compiled from the literature, museum collections and Brazilian Health Departments. Six bioclimatic variables were used as predictors in six ecological niche model algorithms (BIOCLIM, DOMAIN, MaxEnt, GARP, logistic regression and Random Forest). Projections for 2050 used 17 general circulation models in two greenhouse gas representative concentration pathways: “stabilization” and “high increase”. Ensemble models and consensus maps were produced by overlapping binary predictions. Final model outputs showed good performance and significance. The use of species absence data substantially improved model performance. Currently, L. flaviscutellata is widely distributed in the Amazon region, with records in the Atlantic Forest and savannah regions of Central Brazil. Future projections indicate expansion of the climatically suitable area for the vector in both scenarios, towards higher latitudes and elevations. L. flaviscutellata is likely to find increasingly suitable conditions for its expansion into areas where human population size and density are much larger than they are in its current locations. If environmental conditions change as predicted, the range of the vector is likely to expand to southeastern and central-southern Brazil, eastern Paraguay and further into the Amazonian areas of Bolivia, Peru, Ecuador, Colombia and Venezuela. These areas will only become endemic for L. amazonensis, however, if they have competent reservoir hosts and transmission dynamics matching those in the Amazon region.     

A Novel Tool Improves Existing Estimates of Recent Tuberculosis Transmission in Settings of Sparse Data Collection

In any setting, a proportion of incident active tuberculosis (TB) reflects recent transmission (“recent transmission proportion”), whereas the remainder represents reactivation. Appropriately estimating the recent transmission proportion has important implications for local TB control, but existing approaches have known biases, especially where data are incomplete. We constructed a stochastic individual-based model of a TB epidemic and designed a set of simulations (derivation set) to develop two regression-based tools for estimating the recent transmission proportion from five inputs: underlying TB incidence, sampling coverage, study duration, clustered proportion of observed cases, and proportion of observed clusters in the sample. We tested these tools on a set of unrelated simulations (validation set), and compared their performance against that of the traditional ‘n-1’ approach. In the validation set, the regression tools reduced the absolute estimation bias (difference between estimated and true recent transmission proportion) in the ‘n-1’ technique by a median [interquartile range] of 60% [9%, 82%] and 69% [30%, 87%]. The bias in the ‘n-1’ model was highly sensitive to underlying levels of study coverage and duration, and substantially underestimated the recent transmission proportion in settings of incomplete data coverage. By contrast, the regression models’ performance was more consistent across different epidemiological settings and study characteristics. We provide one of these regression models as a user-friendly, web-based tool. Novel tools can improve our ability to estimate the recent TB transmission proportion from data that are observable (or estimable) by public health practitioners with limited available molecular data.     

Advances on genetic rat models of epilepsy

Considering the suitability of laboratory rats in epilepsy research, we and other groups have been developing genetic models of epilepsy in this species. After epileptic rats or seizure-susceptible rats were sporadically found in outbred stocks, the epileptic traits were usually genetically-fixed by selective breeding. So far, the absence seizure models GAERS and WAG/Rij, audiogenic seizure models GEPR-3 and GEPR-9, generalized tonic-clonic seizure models IER, NER and WER, and Canavan-disease related epileptic models TRM and SER have been established. Dissection of the genetic bases including causative genes in these epileptic rat models would be a significant step toward understanding epileptogenesis. N-ethyl-N-nitrosourea (ENU) mutagenesis provides a systematic approach which allowed us to develop two novel epileptic rat models: heat-induced seizure susceptible (Hiss) rats with an Scn1a missense mutation and autosomal dominant lateral temporal epilepsy (ADLTE) model rats with an Lgi1 missense mutation. In addition, we have established episodic ataxia type 1 (EA1) model rats with a Kcna1 missense mutation derived from the ENU-induced rat mutant stock, and identified a Cacna1a missense mutation in a N-Methyl-N-nitrosourea (MNU)-induced mutant rat strain GRY, resulting in the discovery of episodic ataxia type 2 (EA2) model rats. Thus, epileptic rat models have been established on the two paths: ‘phenotype to gene’ and ‘gene to phenotype’. In the near future, development of novel epileptic rat models will be extensively promoted by the use of sophisticated genome editing technologies.     

Climate and pH Predict the Potential Range of the Invasive Apple Snail (Pomacea insularum) in the Southeastern United States

Predicting the potential range of invasive species is essential for risk assessment, monitoring, and management, and it can also inform us about a species’ overall potential invasiveness. However, modeling the distribution of invasive species that have not reached their equilibrium distribution can be problematic for many predictive approaches. We apply the modeling approach of maximum entropy (MaxEnt) that is effective with incomplete, presence-only datasets to predict the distribution of the invasive island apple snail, Pomacea insularum. This freshwater snail is native to South America and has been spreading in the USA over the last decade from its initial introductions in Texas and Florida. It has now been documented throughout eight southeastern states. The snail’s extensive consumption of aquatic vegetation and ability to accumulate and transmit algal toxins through the food web heighten concerns about its spread. Our model shows that under current climate conditions the snail should remain mostly confined to the coastal plain of the southeastern USA where it is limited by minimum temperature in the coldest month and precipitation in the warmest quarter. Furthermore, low pH waters (pH <5.5) are detrimental to the snail’s survival and persistence. Of particular note are low-pH blackwater swamps, especially Okefenokee Swamp in southern Georgia (with a pH below 4 in many areas), which are predicted to preclude the snail’s establishment even though many of these areas are well matched climatically. Our results elucidate the factors that affect the regional distribution of P. insularum, while simultaneously presenting a spatial basis for the prediction of its future spread. Furthermore, the model for this species exemplifies that combining climatic and habitat variables is a powerful way to model distributions of invasive species.     

A comparison of metrics for estimating phylogenetic signal under alternative evolutionary models

Several metrics have been developed for estimating phylogenetic signal in comparative data. These may be important both in guiding future studies on correlated evolution and for inferring broad-scale evolutionary and ecological processes (e.g., phylogenetic niche conservatism). Notwithstanding, the validity of some of these metrics is under debate, especially after the development of more sophisticated model-based approaches that estimate departure from particular evolutionary models (i.e., Brownian motion). Here, two of these model-based metrics (Blomberg’s K-statistics and Pagel’s λ) are compared with three statistical approaches [Moran’s I autocorrelation coefficient, coefficients of determination from the autoregressive method (ARM), and phylogenetic eigenvector regression (PVR)]. Based on simulations of a trait evolving under Brownian motion for a phylogeny with 209 species, we showed that all metrics are strongly, although non-linearly, correlated to each other. Our analyses revealed that statistical approaches provide valid results and may be still particularly useful when detailed phylogenies are unavailable or when trait variation among species is difficult to describe by more standard Brownian or O-U evolutionary models.     

Prediction of Hematopoietic Stem Cell Transplantation Related Mortality- Lessons Learned from the In-Silico Approach: A European Society for Blood and Marrow Transplantation Acute Leukemia Working Party Data Mining Study

Models for prediction of allogeneic hematopoietic stem transplantation (HSCT) related mortality partially account for transplant risk. Improving predictive accuracy requires understating of prediction limiting factors, such as the statistical methodology used, number and quality of features collected, or simply the population size. Using an in-silico approach (i.e., iterative computerized simulations), based on machine learning (ML) algorithms, we set out to analyze these factors. A cohort of 25,923 adult acute leukemia patients from the European Society for Blood and Marrow Transplantation (EBMT) registry was analyzed. Predictive objective was non-relapse mortality (NRM) 100 days following HSCT. Thousands of prediction models were developed under varying conditions: increasing sample size, specific subpopulations and an increasing number of variables, which were selected and ranked by separate feature selection algorithms. Depending on the algorithm, predictive performance plateaued on a population size of 6,611–8,814 patients, reaching a maximal area under the receiver operator characteristic curve (AUC) of 0.67. AUCs’ of models developed on specific subpopulation ranged from 0.59 to 0.67 for patients in second complete remission and receiving reduced intensity conditioning, respectively. Only 3–5 variables were necessary to achieve near maximal AUCs. The top 3 ranking variables, shared by all algorithms were disease stage, donor type, and conditioning regimen. Our findings empirically demonstrate that with regards to NRM prediction, few variables “carry the weight” and that traditional HSCT data has been “worn out”. “Breaking through” the predictive boundaries will likely require additional types of inputs.     

Multi-scale Inference of Interaction Rules in Animal Groups Using Bayesian Model Selection

Inference of interaction rules of animals moving in groups usually relies on an analysis of large scale system behaviour. Models are tuned through repeated simulation until they match the observed behaviour. More recent work has used the fine scale motions of animals to validate and fit the rules of interaction of animals in groups. Here, we use a Bayesian methodology to compare a variety of models to the collective motion of glass prawns (Paratya australiensis). We show that these exhibit a stereotypical ‘phase transition’, whereby an increase in density leads to the onset of collective motion in one direction. We fit models to this data, which range from: a mean-field model where all prawns interact globally; to a spatial Markovian model where prawns are self-propelled particles influenced only by the current positions and directions of their neighbours; up to non-Markovian models where prawns have ‘memory’ of previous interactions, integrating their experiences over time when deciding to change behaviour. We show that the mean-field model fits the large scale behaviour of the system, but does not capture the observed locality of interactions. Traditional self-propelled particle models fail to capture the fine scale dynamics of the system. The most sophisticated model, the non-Markovian model, provides a good match to the data at both the fine scale and in terms of reproducing global dynamics, while maintaining a biologically plausible perceptual range. We conclude that prawns’ movements are influenced by not just the current direction of nearby conspecifics, but also those encountered in the recent past. Given the simplicity of prawns as a study system our research suggests that self-propelled particle models of collective motion should, if they are to be realistic at multiple biological scales, include memory of previous interactions and other non-Markovian effects.     

Crowdsourcing Language Change with Smartphone Applications

Crowdsourcing linguistic phenomena with smartphone applications is relatively new. In linguistics, apps have predominantly been developed to create pronunciation dictionaries, to train acoustic models, and to archive endangered languages. This paper presents the first account of how apps can be used to collect data suitable for documenting language change: we created an app, Dialäkt Äpp (DÄ), which predicts users’ dialects. For 16 linguistic variables, users select a dialectal variant from a drop-down menu. DÄ then geographically locates the user’s dialect by suggesting a list of communes where dialect variants most similar to their choices are used. Underlying this prediction are 16 maps from the historical Linguistic Atlas of German-speaking Switzerland, which documents the linguistic situation around 1950. Where users disagree with the prediction, they can indicate what they consider to be their dialect’s location. With this information, the 16 variables can be assessed for language change. Thanks to the playfulness of its functionality, DÄ has reached many users; our linguistic analyses are based on data from nearly 60,000 speakers. Results reveal a relative stability for phonetic variables, while lexical and morphological variables seem more prone to change. Crowdsourcing large amounts of dialect data with smartphone apps has the potential to complement existing data collection techniques and to provide evidence that traditional methods cannot, with normal resources, hope to gather. Nonetheless, it is important to emphasize a range of methodological caveats, including sparse knowledge of users’ linguistic backgrounds (users only indicate age, sex) and users’ self-declaration of their dialect. These are discussed and evaluated in detail here. Findings remain intriguing nevertheless: as a means of quality control, we report that traditional dialectological methods have revealed trends similar to those found by the app. This underlines the validity of the crowdsourcing method. We are presently extending DÄ architecture to other languages.     

Maximum likelihood parentage assignment using quantitative genotypes

The cost of parentage assignment precludes its application in many selective breeding programmes and molecular ecology studies, and/or limits the circumstances or number of individuals to which it is applied. Pooling samples from more than one individual, and using appropriate genetic markers and algorithms to determine parental contributions to pools, is one means of reducing the cost of parentage assignment. This paper describes and validates a novel maximum likelihood (ML) parentage-assignment method, that can be used to accurately assign parentage to pooled samples of multiple individuals—previously published ML methods are applicable to samples of single individuals only—using low-density single nucleotide polymorphism (SNP) ‘quantitative’ (also referred to as ‘continuous’) genotype data. It is demonstrated with simulated data that, when applied to pools, this ‘quantitative maximum likelihood’ method assigns parentage with greater accuracy than established maximum likelihood parentage-assignment approaches, which rely on accurate discrete genotype calls; exclusion methods; and estimating parental contributions to pools by solving the weighted least squares problem. Quantitative maximum likelihood can be applied to pools generated using either a ‘pooling-for-individual-parentage-assignment’ approach, whereby each individual in a pool is tagged or traceable and from a known and mutually exclusive set of possible parents; or a ‘pooling-by-phenotype’ approach, whereby individuals of the same, or similar, phenotype/s are pooled. Although computationally intensive when applied to large pools, quantitative maximum likelihood has the potential to substantially reduce the cost of parentage assignment, even if applied to pools comprised of few individuals.Subject terms: Plant breeding, Agricultural genetics, Animal breeding, Genetic markers     

Neurokernel: An Open Source Platform for Emulating the Fruit Fly Brain

We have developed an open software platform called Neurokernel for collaborative development of comprehensive models of the brain of the fruit fly Drosophila melanogaster and their execution and testing on multiple Graphics Processing Units (GPUs). Neurokernel provides a programming model that capitalizes upon the structural organization of the fly brain into a fixed number of functional modules to distinguish between these modules’ local information processing capabilities and the connectivity patterns that link them. By defining mandatory communication interfaces that specify how data is transmitted between models of each of these modules regardless of their internal design, Neurokernel explicitly enables multiple researchers to collaboratively model the fruit fly’s entire brain by integration of their independently developed models of its constituent processing units. We demonstrate the power of Neurokernel’s model integration by combining independently developed models of the retina and lamina neuropils in the fly’s visual system and by demonstrating their neuroinformation processing capability. We also illustrate Neurokernel’s ability to take advantage of direct GPU-to-GPU data transfers with benchmarks that demonstrate scaling of Neurokernel’s communication performance both over the number of interface ports exposed by an emulation’s constituent modules and the total number of modules comprised by an emulation.