Strong scaling of general-purpose molecular dynamics simulations on GPUs

We describe a highly optimized implementation of MPI domain decomposition in a GPU-enabled, general-purpose molecular dynamics code, HOOMD-blue (Anderson and Glotzer, arXiv:1308.5587). Our approach is inspired by a traditional CPU-based code, LAMMPS (Plimpton, J. Comp. Phys. 117, 1995), but is implemented within a code that was designed for execution on GPUs from the start (Anderson et al., J. Comp. Phys. 227, 2008). The software supports short-ranged pair force and bond force fields and achieves optimal GPU performance using an autotuning algorithm. We are able to demonstrate equivalent or superior scaling on up to 3,375 GPUs in Lennard-Jones and dissipative particle dynamics (DPD) simulations of up to 108 million particles. GPUDirect RDMA capabilities in recent GPU generations provide better performance in full double precision calculations. For a representative polymer physics application, HOOMD-blue 1.0 provides an effective GPU vs. CPU node speed-up of 12.5x.Comment: 30 pages, 14 figure

A kinematically exact finite element formulation of planar elastic-plastic frames

A finite element formulation of finite deformation static analysis of plane elastic-plastic frames subjected to static loads is presented, in which the only function to be interpolated is the rotation of the centroid axis of the beam. One of the advantages of such a formulation is that the problem of the field-consistency does not arise. Exact non-linear kinematic relationships of the finite-strain beam theory are used, which assume the Bernoulli hypothesis of plane cross-sections. Finite displacements and rotations as well as finite extensional and bending strains are accounted for. The effects of shear strains and non-conservative loads are at present neglected, yet they can simply be incorporated in the formulation. Because the potential energy of internal forces does not exist with elastic-plastic material, the principle of virtual work is introduced as the basis of the finite element formulation. A generalized principle of virtual work is proposed in which the displacements, rotation, extensional and bending strains, and the Lagrangian multipliers are independent variables. By exploiting the special structure of the equations of the problem, the displacements, the strains and the multipliers are eliminated from the generalized principle of virtual work. A novel principle is obtained in which the rotation becomes the only function to be approximated in its finite element implementation. It is shown that (N-1)-point numerical integration must be employed in conjunction with N-node interpolation polynomials for the rotation, and the Lobatto rule is recommended. Regarding the integration over the cross-section, it is demonstrated by numerical examples that, due to discontinuous integrands, no integration order defined as `computationally efficient yet accurate enough' could be suggested. The theoretical findings and a nice performance of the derived finite elements are illustrated by numerical examples

The genetic basis of onset age in schizophrenia: evidence and models

Schizophrenia is a heritable neurocognitive disorder affecting about 1% of the population, and usually has an onset age at around 21–25 in males and 25–30 in females. Recent advances in genetics have helped to identify many common and rare variants for the liability to schizophrenia. Earlier evidence appeared to suggest that younger onset age is associated with higher genetic liability to schizophrenia. Clinical longitudinal research also found that early and very-early onset schizophrenia are associated with poor clinical, neurocognitive, and functional profiles. A recent study reported a heritability of 0.33 for schizophrenia onset age, but the genetic basis of this trait in schizophrenia remains elusive. In the pre-Genome-Wide Association Study (GWAS) era, genetic loci found to be associated with onset age were seldom replicated. In the post-Genome-Wide Association Study era, new conceptual frameworks are needed to clarify the role of onset age in genetic research in schizophrenia, and to identify its genetic basis. In this review, we first discussed the potential of onset age as a characterizing/subtyping feature for psychosis, and as an important phenotypic dimension of schizophrenia. Second, we reviewed the methods, samples, findings and limitations of previous genetic research on onset age in schizophrenia. Third, we discussed a potential conceptual framework for studying the genetic basis of onset age, as well as the concepts of susceptibility, modifier, and “mixed” genes. Fourth, we discussed the limitations of this review. Lastly, we discussed the potential clinical implications for genetic research of onset age of schizophrenia, and how future research can unveil the potential mechanisms for this trait

Predicting first-episode psychosis patients who will never relapse over 10 years.

BACKGROUND: Although relapse in psychosis is common, a small proportion of patients will not relapse in the long term. We examined the proportion and predictors of patients who never relapsed in the 10 years following complete resolution of positive symptoms from their first psychotic episode. METHOD: Patients who previously enrolled in a 12-month randomized controlled trial on medication discontinuation and relapse following first-episode psychosis (FEP) were followed up after 10 years. Relapse of positive symptoms was operationalized as a change from a Clinical Global Impression scale positive score of <3 for at least 3 consecutive months to a score of ⩾3 (mild or more severe). Baseline predictors included basic demographics, premorbid functioning, symptoms, functioning, and neurocognitive functioning. RESULTS: Out of 178 first-episode patients, 37 (21%) never relapsed during the 10-year period. Univariate predictors (p ⩽ 0.1) of patients who never relapsed included a duration of untreated psychosis (DUP) ⩽30 days, diagnosed with non-schizophrenia spectrum disorders, having less severe negative symptoms, and performing better in logical memory immediate recall and verbal fluency tests. A multivariate logistic regression analysis further suggested that the absence of any relapsing episodes was significantly related to better short-term verbal memory, shorter DUP, and non-schizophrenia spectrum disorders. CONCLUSIONS: Treatment delay and neurocognitive function are potentially modifiable predictors of good long-term prognosis in FEP. These predictors are informative as they can be incorporated into an optimum risk prediction model in the future, which would help with clinical decision making regarding maintenance treatment in FEP

The Burden of Invasive Bacterial Infections in Pemba, Zanzibar

BACKGROUND: We conducted a surveillance study to determine the leading causes of bloodstream infection in febrile patients seeking treatment at three district hospitals in Pemba Island, Zanzibar, Tanzania, an area with low malaria transmission. METHODS: All patients above two months of age presenting to hospital with fever were screened, and blood was collected for microbiologic culture and malaria testing. Bacterial sepsis and malaria crude incidence rates were calculated for a one-year period and were adjusted for study participation and diagnostic sensitivity of blood culture. RESULTS: Blood culture was performed on 2,209 patients. Among them, 166 (8%) samples yielded bacterial growth; 87 (4%) were considered as likely contaminants; and 79 (4%) as pathogenic bacteria. The most frequent pathogenic bacteria isolated were Salmonella Typhi (n = 46; 58%), followed by Streptococcus pneumoniae (n = 12; 15%). The crude bacteremia rate was 6/100,000 but when adjusted for potentially missed cases the rate may be as high as 163/100,000. Crude and adjusted rates for S. Typhi infections and malaria were 4 and 110/100,000 and 4 and 47/100,000, respectively. Twenty three (51%), 22 (49%) and 22 (49%) of the S. Typhi isolates were found to be resistant toward ampicillin, chloramphenicol and cotrimoxazole, respectively. Multidrug resistance (MDR) against the three antimicrobials was detected in 42% of the isolates. CONCLUSIONS: In the presence of very low malaria incidence we found high rates of S. Typhi and S. pneumoniae infections on Pemba Island, Zanzibar. Preventive measures such as vaccination could reduce the febrile disease burden

Cost of Illness Due to Typhoid Fever in Pemba, Zanzibar, East Africa

The aim of this study was to estimate the economic burden of typhoid fever in Pemba, Zanzibar, East Africa. This study was an incidence-based cost-of-illness analysis from a societal perspective. It covered new episodes of blood culture-confirmed typhoid fever in patients presenting at the outpatient or inpatient departments of three district hospitals between May 2010 and December 2010. Cost of illness was the sum of direct costs and costs for productivity loss. Direct costs covered treatment, travel, and meals. Productivity costs were loss of income by patients and caregivers. The analysis included 17 episodes. The mean age of the patients, was 23 years (range=5-65, median=22). Thirty-five percent were inpatients, with a mean of 4.75 days of hospital stay (range=3-7, median=4.50). The mean cost for treatment alone during hospital care was US$21.97 at 2010 prices (US$ 1=1,430.50 Tanzanian Shilling\u2500TSH). The average societal cost was US$154.47 per typhoid episode. The major expenditure was productivity cost due to lost wages of US$ 128.02 (83%). Our results contribute to the further economic evaluation of typhoid fever vaccination in Zanzibar and other sub-Saharan African countries

Whole exome sequencing coupled with unbiased functional analysis reveals new Hirschsprung disease genes

Background: Hirschsprung disease (HSCR), which is congenital obstruction of the bowel, results from a failure of enteric nervous system (ENS) progenitors to migrate, proliferate, differentiate, or survive within the distal intestine. Previous studies that have searched for genes underlying HSCR have focused on ENS-related pathways and genes not fitting the current knowledge have thus often been ignored. We identify and validate novel HSCR genes using whole exome sequencing (WES), burden tests, in silico prediction, unbiased in vivo analyses of the mutated genes in zebrafish, and expression analyses in zebrafish, mouse, and human. Results: We performed de novo mutation (DNM) screening on 24 HSCR trios. We identify 28 DNMs in 21 different genes. Eight of the DNMs we identified occur in RET, the main HSCR gene, and the remaining 20 DNMs reside in genes not reported in the ENS. Knockdown of all 12 genes with missense or loss-of-function DNMs showed that the orthologs of four genes (DENND3, NCLN, NUP98, and TBATA) are indispensable for ENS development in zebrafish, and these results were confirmed by CRISPR knockout. These genes are also expressed in human and mouse gut and/or ENS progenitors. Importantly, the encoded proteins are linked to neuronal processes shared by the central nervous system and the ENS. Conclusions: Our data open new fields of investigation into HSCR pathology and provide novel insights into the development of the ENS. Moreover, the study demonstrates that functional analyses of genes carrying DNMs are warranted to delineate the full genetic architecture of rare complex diseases

Fine Mapping of the NRG1 Hirschsprung's Disease Locus

The primary pathology of Hirschsprung's disease (HSCR, colon aganglionosis) is the absence of ganglia in variable lengths of the hindgut, resulting in functional obstruction. HSCR is attributed to a failure of migration of the enteric ganglion precursors along the developing gut. RET is a key regulator of the development of the enteric nervous system (ENS) and the major HSCR-causing gene. Yet the reduced penetrance of RET DNA HSCR-associated variants together with the phenotypic variability suggest the involvement of additional genes in the disease. Through a genome-wide association study, we uncovered a ∼350 kb HSCR-associated region encompassing part of the neuregulin-1 gene (NRG1). To identify the causal NRG1 variants contributing to HSCR, we genotyped 243 SNPs variants on 343 ethnic Chinese HSCR patients and 359 controls. Genotype analysis coupled with imputation narrowed down the HSCR-associated region to 21 kb, with four of the most associated SNPs (rs10088313, rs10094655, rs4624987, and rs3884552) mapping to the NRG1 promoter. We investigated whether there was correlation between the genotype at the rs10088313 locus and the amount of NRG1 expressed in human gut tissues (40 patients and 21 controls) and found differences in expression as a function of genotype. We also found significant differences in NRG1 expression levels between diseased and control individuals bearing the same rs10088313 risk genotype. This indicates that the effects of NRG1 common variants are likely to depend on other alleles or epigenetic factors present in the patients and would account for the variability in the genetic predisposition to HSCR