ALBERT

Description: Interactions between modular domains and short linear motifs (3–10 amino acids peptide stretches) are crucial for cell signaling. The motifs typically reside in the disordered regions of the proteome and the interactions are often transient, allowing for rapid changes in response to changing stimuli. The properties that make domain-motif interactions suitable for cell signaling also make them difficult to capture experimentally and they are therefore largely underrepresented in the known protein-protein interaction networks. Most of the knowledge on domain-motif interactions is derived from low-throughput studies, although there exist dedicated high-throughput methods for the identification of domain-motif interactions. The methods include arrays of peptides or proteins, display of peptides on phage or yeast, and yeast-two-hybrid experiments. We here provide a survey of scalable methods for domain-motif interaction profiling. These methods have frequently been applied to a limited number of ubiquitous domain families. It is now time to apply them to a broader set of peptide binding proteins, to provide a comprehensive picture of the linear motifs in the human proteome and to link them to their potential binding partners. Despite the plethora of methods, it is still a challenge for most approaches to identify interactions that rely on post-translational modification or context dependent or conditional interactions, suggesting directions for further method development.

Description: Background: Digoxin intoxication results in predominantly digestive, cardiac and neurological symptoms. This case is outstanding in that the intoxication occurred in a nonagenarian and induced severe, extensively documented visual symptoms as well as dysphagia and proprioceptive illusions. Moreover, it went undiagnosed for a whole month despite close medical follow-up, illustrating the difficulty in recognizing drug-induced effects in a polymorbid patient.Case presentationDigoxin 0.25 mg qd for atrial fibrillation was prescribed to a 91-year-old woman with an estimated creatinine clearance of 18 ml/min. Over the following 2–3 weeks she developed nausea, vomiting and dysphagia, snowy and blurry vision, photopsia, dyschromatopsia, aggravated pre-existing formed visual hallucinations and proprioceptive illusions. She saw her family doctor twice and visited the eye clinic once until, 1 month after starting digoxin, she was admitted to the emergency room. Intoxication was confirmed by a serum digoxin level of 5.7 ng/ml (reference range 0.8–2 ng/ml). After stopping digoxin, general symptoms resolved in a few days, but visual complaints persisted. Examination by the ophthalmologist revealed decreased visual acuity in both eyes, 4/10 in the right eye (OD) and 5/10 in the left eye (OS), decreased color vision as demonstrated by a score of 1/13 in both eyes (OU) on Ishihara pseudoisochromatic plates, OS cataract, and dry age-related macular degeneration (ARMD). Computerized static perimetry showed non-specific diffuse alterations suggestive of either bilateral retinopathy or optic neuropathy. Full-field electroretinography (ERG) disclosed moderate diffuse rod and cone dysfunction and multifocal ERG revealed central loss of function OU. Visual symptoms progressively improved over the next 2 months, but multifocal ERG did not. The patient was finally discharged home after a 5 week hospital stay. Conclusion: This case is a reminder of a complication of digoxin treatment to be considered by any treating physician. If digoxin is prescribed in a vulnerable patient, close monitoring is mandatory. In general, when facing a new health problem in a polymorbid patient, it is crucial to elicit a complete history, with all recent drug changes and detailed complaints, and to include a drug adverse reaction in the differential diagnosis.

Description: Background: Biological pathways are descriptive diagrams of biological processes widely used for functional analysis of differentially expressed genes or proteins. Primary data analysis, such as quality control, normalisation, and statistical analysis, is often performed in scripting languages like R, Perl, and Python. Subsequent pathway analysis is usually performed using dedicated external applications. Workflows involving manual use of multiple environments are time consuming and error prone. Therefore, tools are needed that enable pathway analysis directly within the same scripting languages used for primary data analyses. Existing tools have limited capability in terms of available pathway content, pathway editing and visualisation options, and export file formats. Consequently, making the full-fledged pathway analysis tool PathVisio available from various scripting languages will benefit researchers. Results: We developed PathVisioRPC, an XMLRPC interface for the pathway analysis software PathVisio. PathVisioRPC enables creating and editing biological pathways, visualising data on pathways, performing pathway statistics, and exporting results in several image formats in multiple programming environments.We demonstrate PathVisioRPC functionalities using examples in Python. Subsequently, we analyse a publicly available NCBI GEO gene expression dataset studying tumour bearing mice treated with cyclophosphamide in R. The R scripts demonstrate how calls to existing R packages for data processing and calls to PathVisioRPC can directly work together. To further support R users, we have created RPathVisio simplifying the use of PathVisioRPC in this environment. We have also created a pathway module for the microarray data analysis portal ArrayAnalysis.org that calls the PathVisioRPC interface to perform pathway analysis. This module allows users to use PathVisio functionality online without having to download and install the software and exemplifies how the PathVisioRPC interface can be used by data analysis pipelines for functional analysis of processed genomics data. Conclusions: PathVisioRPC enables data visualisation and pathway analysis directly from within various analytical environments used for preliminary analyses. It supports the use of existing pathways from WikiPathways or pathways created using the RPC itself. It also enables automation of tasks performed using PathVisio, making it useful to PathVisio users performing repeated visualisation and analysis tasks. PathVisioRPC is freely available for academic and commercial use at http://projects.bigcat.unimaas.nl/pathvisiorpc.

Description: Background: Mycoplasma pneumoniae is a common pathogen that causes upper and lower respiratory tract infections in people of all ages, responsible for up to 40 % of community-acquired pneumonias. It also causes a wide array of extrapulmonary infections and autoimmune phenomena. Phylogenetic studies of the organism have been generally restricted to specific genes or regions of the genome, because whole genome sequencing has been completed for only 4 strains. To better understand the physiology and pathogenicity of this important human pathogen, we performed comparative genomic analysis of 15 strains of M. pneumoniae that were isolated between the 1940s to 2009 from respiratory specimens and cerebrospinal fluid originating from the USA, China and England. Results: Illumina MiSeq whole genome sequencing was performed on the 15 strains and all genome sequences were completed. Results from the comparative genomic analysis indicate that although about 1500 SNP and indel variants exist between type1 and type 2 strains, there is an overall high degree of sequence similarity among the strains (〉99 % identical to each other). Within the two subtypes, conservation of most genes, including the CARDS toxin gene and arginine deiminase genes, was observed. The major variation occurs in the P1 and ORF6 genes associated with the adhesin complex. Multiple hsdS genes (encodes S subunit of type I restriction enzyme) with variable tandem repeat copy numbers were found in all 15 genomes. Conclusions: These data indicate that despite conclusions drawn from 16S rRNA sequences suggesting rapid evolution, the M. pneumoniae genome is extraordinarily stable over time and geographic distance across the globe with a striking lack of evidence of horizontal gene transfer.

Description: Background: Distinguishing between viable and dead bacteria in animal and urban effluents is a major challenge. Among existing methods, propidium monoazide (PMA)-qPCR is a promising way to quantify viable cells. However, its efficiency depends on the composition of the effluent, particularly on total suspended solids (TSS)) and on methodological parameters. The aim of this study was evaluate the influence of three methodological factors (concentration of PMA, incubation time and photoactivation time) on the efficiency of PMA-qPCR to quantify viable and dead cells of Listeria monocytogenes used as a microorganism model, in two piggery effluents (manure and lagoon effluent containing 20 and 0.4 TSS g.kg −1 , respectively). An experimental design strategy (Doehlert design and desirability function) was used to identify the experimental conditions to achieve optimal PMA-qPCR results. Results: The quantification of viable cells of L. monocytogenes was mainly influenced by the concentration of PMA in the manure and by the duration of photoactivation in the lagoon effluent. Optimal values differed with the matrix: 55 μM PMA, 5 min incubation and 56 min photoactivation for manure and 20 μM PMA, 20 min incubation and 30 min photoactivation for lagoon effluent. Applied to five manure and four lagoon samples, these conditions resulted in satisfactory quantification of viable and dead cells. Conclusion: PMA-qPCR can be used on undiluted turbid effluent with high levels of TSS, provided preliminary tests are performed to identify the optimal conditions.

Description: Background: Walnut (Juglans regia, 2n = 32, approximately 606 Mb per 1C genome) is an economically important tree crop. Resistance to anthracnose, caused by Colletotrichum gloeosporioides, is a major objective of walnut genetic improvement in China. The recently developed specific length amplified fragment sequencing (SLAF-seq) is an efficient strategy that can obtain large numbers of markers with sufficient sequence information to construct high-density genetic maps and permits detection of quantitative trait loci (QTLs) for molecular breeding. Results: SLAF-seq generated 161.64 M paired-end reads. 153,820 SLAF markers were obtained, of which 49,174 were polymorphic. 13,635 polymorphic markers were sorted into five segregation types and 2,577 markers of them were used to construct genetic linkage maps: 2,395 of these fell into 16 linkage groups (LGs) for the female map, 448 markers for the male map, and 2,577 markers for the integrated map. Taking into account the size of all LGs, the marker coverage was 2,664.36 cM for the female map, 1,305.58 cM for the male map, and 2,457.82 cM for the integrated map. The average intervals between two adjacent mapped markers were 1.11 cM, 2.91 cM and 0.95 cM for three maps, respectively. ‘SNP_only’ markers accounted for 89.25 % of the markers on the integrated map. Mapping markers contained 5,043 single nucleotide polymorphisms (SNPs) loci, which corresponded to two SNP loci per SLAF marker. According to the integrated map, we used interval mapping (Logarithm of odds, LOD 〉 3.0) to detect our quantitative trait. One QTL was detected for anthracnose resistance. The interval of this QTL ranged from 165.51 cM to 176.33 cM on LG14, and ten markers in this interval that were above the threshold value were considered to be linked markers to the anthracnose resistance trait. The phenotypic variance explained by each marker ranged from 16.2 to 19.9 %, and their LOD scores varied from 3.22 to 4.04. Conclusions: High-density genetic maps for walnut containing 16 LGs were constructed using the SLAF-seq method with an F1 population. One QTL for walnut anthracnose resistance was identified based on the map. The results will aid molecular marker-assisted breeding and walnut resistance genes identification.

Description: Background: Detecting and quantifying isoforms from RNA-seq data is an important but challenging task. The problem is often ill-posed, particularly at low coverage. One promising direction is to exploit several samples simultaneously. Results: We propose a new method for solving the isoform deconvolution problem jointly across several samples. We formulate a convex optimization problem that allows to share information between samples and that we solve efficiently. We demonstrate the benefits of combining several samples on simulated and real data, and show that our approach outperforms pooling strategies and methods based on integer programming. Conclusion: Our convex formulation to jointly detect and quantify isoforms from RNA-seq data of multiple related samples is a computationally efficient approach to leverage the hypotheses that some isoforms are likely to be present in several samples. The software and source code are available at http://cbio.ensmp.fr/flipflop.

Description: Background: Dual-specificity phosphatase-5 (DUSP5) plays a central role in vascular development and disease. We present a p-nitrophenol phosphate (pNPP) based enzymatic assay to screen for inhibitors of the phosphatase domain of DUSP5. Methods: pNPP is a mimic of the phosphorylated tyrosine on the ERK2 substrate (pERK2) and binds the DUSP5 phosphatase domain with a K m of 7.6 ± 0.4 mM. Docking followed by inhibitor verification using the pNPP assay identified a series of polysulfonated aromatic inhibitors that occupy the DUSP5 active site in the region that is likely occupied by the dual-phosphorylated ERK2 substrate tripeptide (pThr-Glu-pTyr). Secondary assays were performed with full length DUSP5 with ERK2 as substrate. Results: The most potent inhibitor has a naphthalene trisulfonate (NTS) core. A search for similar compounds in a drug database identified suramin, a dimerized form of NTS. While suramin appears to be a potent and competitive inhibitor (25 ± 5 μM), binding to the DUSP5 phosphatase domain more tightly than the monomeric ligands of which it is comprised, it also aggregates. Further ligand-based screening, based on a pharmacophore derived from the 7 Å separation of sulfonates on inhibitors and on sulfates present in the DUSP5 crystal structure, identified a disulfonated and phenolic naphthalene inhibitor (CSD 3 _2320) with IC 50 of 33 μM that is similar to NTS and does not aggregate. Conclusions: The new DUSP5 inhibitors we identify in this study typically have sulfonates 7 Å apart, likely positioning them where the two phosphates of the substrate peptide (pThr-Glu-pTyr) bind, with one inhibitor also positioning a phenolic hydroxyl where the water nucleophile may reside. Polysulfonated aromatic compounds do not commonly appear in drugs and have a tendency to aggregate. One FDA-approved polysulfonated drug, suramin, inhibits DUSP5 and also aggregates. Docking and modeling studies presented herein identify polysulfonated aromatic inhibitors that do not aggregate, and provide insights to guide future design of mimics of the dual-phosphate loops of the ERK substrates for DUSPs.

Description: Background: Although principal component analysis (PCA) is widely used for the dimensional reduction of biomedical data, interpretation of PCA results remains daunting. Most existing interpretation methods attempt to explain each principal component (PC) in terms of a small number of variables by generating approximate PCs with mainly zero loadings. Although useful when just a few variables dominate the population PCs, these methods can perform poorly on genomic data, where interesting biological features are frequently represented by the combined signal of functionally related sets of genes. While gene set testing methods have been widely used in supervised settings to quantify the association of groups of genes with clinical outcomes, these methods have seen only limited application for testing the enrichment of gene sets relative to sample PCs. Results: We describe a novel approach, principal component gene set enrichment (PCGSE), for unsupervised gene set testing relative to the sample PCs of genomic data. The PCGSE method computes the statistical association between gene sets and individual PCs using a two-stage competitive gene set test. To demonstrate the efficacy of the PCGSE method, we use simulated and real gene expression data to evaluate the performance of various gene set test statistics and significance tests. Conclusions: Gene set testing is an effective approach for interpreting the PCs of high-dimensional genomic data. As shown using both simulated and real datasets, the PCGSE method can generate biologically meaningful and computationally efficient results via a two-stage, competitive parametric test that correctly accounts for inter-gene correlation.

Description: Background: The cascade computer model (CCM) was designed as a machine-learning feature platform for prediction of drug diffusivity from the mucoadhesive formulations. Three basic models (the statistical regression model, the K nearest neighbor model and the modified version of the back propagation neural network) in CCM operate sequentially in close collaboration with each other, employing the estimated value obtained from the afore-positioned base model as an input value to the next-positioned base model in the cascade.The effects of various parameters on the pharmacological efficacy of a female controlled drug delivery system (FcDDS) intended for prevention of women from HIV-1 infection were evaluated using an in vitro apparatus “Simulant Vaginal System” (SVS). We used computer simulations to explicitly examine the changes in drug diffusivity from FcDDS and determine the prognostic potency of each variable for in vivo prediction of formulation efficacy. The results obtained using the CCM approach were compared with those from individual multiple regression model. Results: CCM significantly lowered the percentage mean error (PME) and enhanced r 2 values as compared with those from the multiple regression models. It was noted that CCM generated the PME value of 21.82 at 48169 epoch iterations, which is significantly improved from the PME value of 29.91 % at 118344 epochs by the back propagation network model. The results of this study indicated that the sequential ensemble of the classifiers allowed for an accurate prediction of the domain with significantly lowered variance and considerably reduces the time required for training phase. Conclusion: CCM is accurate, easy to operate, time and cost-effective, and thus, can serve as a valuable tool for prediction of drug diffusivity from mucoadhesive formulations. CCM may yield new insights into understanding how drugs are diffused from the carrier systems and exert their efficacies under various clinical conditions.

Description: Background: In family-based association analysis, each family is typically ascertained from a single proband, which renders the effects of ascertainment bias heterogeneous among family members. This is contrary to case–control studies, and may introduce sample or ascertainment bias. Statistical efficiency is affected by ascertainment bias, and careful adjustment can lead to substantial improvements in statistical power. However, genetic association analysis has often been conducted using family-based designs, without addressing the fact that each proband in a family has had a great influence on the probability for each family member to be affected.MethodWe propose a powerful and efficient statistic for genetic association analysis that considered the heterogeneity of ascertainment bias among family members, under the assumption that both prevalence and heritability of disease are available. With extensive simulation studies, we showed that the proposed method performed better than the existing methods, particularly for diseases with large heritability. Results: We applied the proposed method to the genome-wide association analysis of Alzheimer’s disease. Four significant associations with the proposed method were found. Conclusion: Our significant findings illustrated the practical importance of this new analysis method.

Description: Background: The full-fat flaxseed (Linum usitatissimum L.) meal has obtained relatively new flourished concept as food or feedstuff for the development of healthier products. It provides favorable balance of polyunsaturated, monounsaturated and saturated fatty acids. However, flaxseed meal may be susceptible to oxidation by exposure to various storage conditions which is extremely undesirable and produces toxic compounds to human health. Another consideration in the application of flaxseed meal relates to the presence of anti-nutritional compounds that need to be minimized using appropriate processing method. The present research work was conducted to evaluate the impact of extrusion processing conditions and storage of full-fat flaxseed meal on functional characteristics such as α-linolenic acid content, lipid peroxidation and sensory attributes. Methods: The raw flaxseed meal was analyzed for cyanogenic glycosides, tannin and mucilage anti-nutritional compounds. Fatty acids composition was quantified by gas chromatography. The meal was extruded at barrel exit temperature (100–140 °C), screw speed (50–150 rpm), feed rate (30–90 kg/h) and feed moisture (10–30 %) for reduction of anti-nutritional compounds. The raw and extruded meals were stored for a ninety-day period under room conditions (20–25 °C). Lipid peroxidation was analyzed by peroxide, free fatty acids, conjugated dienes, total volatiles and malondialdehyde assay. Color, aroma and overall acceptability attributes were evaluated by sensory multiple comparison tests. Results: The raw flaxseed meal possessed significant amount of anti-nutritional compounds, lipid and α-linolenic acid contents. The extrusion processing at high barrel exit temperature (140 °C) significantly reduced the cyanogenic compounds (84 %), tannin (73 %) and mucilage (27 %) in the flaxseed meal. The α-linolenic acid content and lipid peroxidation did not significantly change after extrusion processing or during storage at the end of 60 days. Fluctuations in sensory attributes occurred during storage, but at the end of 90 days, only the extruded samples presented negative effect and showed lowest consumer acceptability. Conclusions: The present study suggested that extrusion of flaxseed meal at optimum conditions and stored for 60 days did not change the stability of full-fat flaxseed meal and can be used as supplement or ingredient for the production of various healthier products.

Description: Background: Leukoencephalopathy with Vanishing White Matter (VWM) is an autosomal recessive disorder caused by germline mutations in the genes EIF2B1-5, which encode the 5 subunits of the eukaryotic translation initiation factor eIF2B. To date, analysis of the biochemical effects of mutations in the EIF2B2-5 genes has been carried out, but no study has been performed on mutations in the EIF2B1 gene. This gene encodes eIF2Bα, the smallest subunit in eIF2B which has an important role in both the structure and regulation of the eIF2B complex. Methods: eIF2B subunits were overexpressed in HEK293 cells and isolated from the resulting cell lysates by affinity chromatography. Formation of the eIF2B complex and binding of its substrate, eIF2, was assessed by western blot. Assays of the guanine nucleotide exchange (GEF) activity were also carried out. Results: Of the 5 eIF2Bα mutations studied, we found 3 that showed loss or reduction of binding of eIF2Bα to the rest of the complex, one with increased GEF activity, and one where no effects on activity or complex formation were observed. Conclusions: This is the first study on eIF2Bα VWM mutations. We show that some mutations cause expected decreases in GEF activity or complex formation, similar to a majority of observed VWM mutations. However, we also observe some unexpected changes which hint at other effects of these mutations on as yet undescribed functions of eIF2B.

Description: Background: Ataxia with oculomotor apraxia type 1 is an autosomal-recessive neurodegenerative disorder characterized by a childhood onset of slowly progressive cerebellar ataxia, followed by oculomotor apraxia and a severe primary motor peripheral axonal motor neuropathy. Ataxia with oculomotor apraxia type 1 is caused by bi-allelic mutations in APTX (chromosome 9p21.1).Case presentationOur patient has a clinical presentation that is typical for ataxia with oculomotor apraxia type 1 with no particularly severe phenotype. Multiplex Ligation-dependent Probe Amplification analysis resulted in the identification of a homozygous deletion of all coding APTX exons (3 to 9). SNP array analysis using the Illumina Infinium CytoSNP-850 K microarray indicated that the deletion was about 62 kb. Based on the SNP array results, the breakpoints were found using direct sequence analysis: c.-5 + 1225_*44991del67512, p.0?. Both parents were heterozygous for the deletion. Homozygous complete APTX deletions have been described in literature for two other patients. We obtained a sample from one of these two patients and characterized the deletion (156 kb) as c.-23729_*115366del155489, p.0?, including the non-coding exons 1A and 2 of APTX. The more severe phenotype reported for this patient is not observed in our patient. It remains unclear whether the larger size of the deletion (156 kb vs 62 kb) plays a role in the phenotype (no extra genes are deleted). Conclusion: Here we described an ataxia with oculomotor apraxia type 1 patient who has a homozygous deletion of the complete coding region of APTX. In contrast to the patient with the large deletion, our patient does not have a severe phenotype. More patients with deletions of APTX are required to investigate a genotype-phenotype effect.

Description: Background: Maize is a major staple food crop globally and contains various concentrations of vitamins. Folates are essential water-soluble B-vitamins that play an important role as one-carbon (C1) donors and acceptors in organisms. To gain an understanding of folate metabolism in maize, we performed an intensive in silico analysis to screen for genes involved in folate metabolism using publicly available databases, followed by examination of the transcript expression patterns and profiling of the folate derivatives in the kernels of two maize inbred lines. Results: A total of 36 candidate genes corresponding to 16 folate metabolism-related enzymes were identified. The maize genome contains all the enzymes required for folate and C1 metabolism, characterized by highly conserved functional domains across all the other species investigated. Phylogenetic analysis revealed that these enzymes in maize are conserved throughout evolution and have a high level of similarity with those in sorghum and millet. The LC-MS analyses of two maize inbred lines demonstrated that 5-methyltetrahydrofolate was the major form of folate derivative in young seeds, while 5-formyltetrahydrofolate in mature seeds. Most of the genes involved in folate and C1 metabolism exhibited similar transcriptional expression patterns between these two maize lines, with the highest transcript abundance detected on day after pollination (DAP) 6 and the decreased transcript abundance on DAP 12 and 18. Compared with the seeds on DAP 30, 5-methyltetrahydrofolate was decreased and 5-formyltetrahydrofolate was increased sharply in the mature dry seeds. Conclusions: The enzymes involved in folate and C1 metabolism are conserved between maize and other plant species. Folate and C1 metabolism is active in young developing maize seeds at transcriptional levels.

Description: Background: The nucleotide second messengers cAMP and c-di-GMP allow many bacteria, including the human intestinal pathogen Vibrio cholerae, to respond to environmental stimuli with appropriate physiological adaptations. In response to limitation of specific carbohydrates, cAMP and its receptor CRP control the transcription of genes important for nutrient acquisition and utilization; c-di-GMP controls the transition between motile and sessile lifestyles often, but not exclusively, through transcriptional mechanisms. In this study, we investigated the convergence of cAMP and c-di-GMP signaling pathways in regulating the expression of gbpA. GbpA is a colonization factor that participates in the attachment of V. cholerae to N-acetylglucosamine-containing surfaces in its native aquatic environment and the host intestinal tract. Results: We show that c-di-GMP inhibits gbpA activation in a fashion independent of the known transcription factors that directly sense c-di-GMP. Interestingly, inhibition of gbpA activation by c-di-GMP only occurs during growth on non-PTS dependent nutrient sources. Consistent with this result, we show that CRP binds to the gbpA promoter in a cAMP-dependent manner in vitro and drives transcription of gbpA in vivo. The interplay between cAMP and c-di-GMP does not broadly impact the CRP-cAMP regulon, but occurs more specifically at the gbpA promoter. Conclusions: These findings suggest that c-di-GMP directly interferes with the interaction of CRP-cAMP and the gbpA promoter via an unidentified regulator. The use of two distinct second messenger signaling mechanisms to regulate gbpA transcription may allow V. cholerae to finely modulate GbpA production, and therefore colonization of aquatic and host surfaces, in response to discrete environmental stimuli.

Description: Background: Appropriate utilization and compliance of Surgical Safety Checklist reduces occurrence of perioperative surgical complications and improve patient outcomes. However, data on compliance of surgical checklists are scarce in the study area. Therefore, the aim of this study was to evaluate compliance of checklist completion and its barrier for utilization at University of Gondar Hospital, Northwest Ethiopia. Methods: A prospective observational study was conducted among 282 patients undergoing elective and emergency surgery from January to March 2013. Compliance and completeness rate with implementation of Sign-in, Time-out, and Sign-out domains was computed with SPSS 20 package. Results: A total of 282 operations were performed and checklists were utilized in 39.7 % (112/282) of cases. Among these, most checklists were employed during emergency procedures (61.6 %) that need general anesthesia (75.9 %) in department of surgery (58.9 %). The overall compliance and completeness rate were 39.7 and 63.4 % respectively. The sign-in, time-out and sign-out were missed in 30.5 % (273/896), 35.4 % (436/1,232) and 45.7 % (307/672) respectively. The main reasons cited for non-user were lack of previous training (45.1 %) and lack of cooperation among surgical team members (21.6 %).Conclusions and recommendationsThe completeness rate was satisfactory but the overall compliance rate was suboptimal. An instrument that is used 40 % of the time has been a fairly basic introduction without significant reinforcement training. Moreover, frequent use of the checklist during emergency cases has been deemed to be of value by clinicians. Supplementary training and attention to actual checklist use would be indicated to ensure that this valuable tool could be used more routinely and improve communication. Conducting regular audit of checklist utilization is also recommended.

Description: Background: Butterflies of the subtribe Mycalesina have radiated successfully in almost all habitat types in Africa, Madagascar, the Indian subcontinent, Indo-China and Australasia. Studies aimed at understanding the reasons behind the evolutionary success of this spectacular Old World butterfly radiation have been hampered by the lack of a stable phylogeny for the group. Here, we have reconstructed a robust phylogenetic framework for the subtribe using 10 genes from 195 exemplar taxa. Results: We recovered seven well supported clades within the subtribe corresponding to the five traditional genera (Lohora, Heteropsis, Hallelesis, Bicyclus, Mycalesis), one as recently revised (Mydosama) and one newly revised genus (Culapa). The phylogenetic relationships of these mycalesine genera have been robustly established for the first time. Within the proposed phylogenetic framework, we estimated the crown age of the subtribe to be 40 Million years ago (Mya) and inferred its ultimate origin to be in Asia. Our results reveal both vicariance and dispersal as factors responsible for the current widespread distribution of the group in the Old World tropics. We inferred that the African continent has been colonized at least twice by Asian mycalesines within the last 26 and 23 Mya. In one possible scenario, an Asian ancestor gave rise to Heteropsis on continental Africa, which later dispersed into Madagascar and most likely back colonised Asia. The second colonization of Africa by Asian ancestors resulted in Hallelesis and Bicyclus on continental Africa, the descendants of which did not colonise other regions but rather diversified only in continental Africa. The genera Lohora and Mydosama are derivatives of ancestors from continental Asia. Conclusion: Our proposed time-calibrated phylogeny now provides a solid framework within which we can implement mechanistic studies aimed at unravelling the ecological and evolutionary processes that culminated in the spectacular radiation of mycalesines in the Old World tropics.

Description: Background: In many domains, scientists build complex simulators of natural phenomena that encode their hypotheses about the underlying processes. These simulators can be deterministic or stochastic, fast or slow, constrained or unconstrained, and so on. Optimizing the simulators with respect to a set of parameter values is common practice, resulting in a single parameter setting that minimizes an objective subject to constraints. Results: We propose algorithms for post optimization posterior evaluation (POPE) of simulators. The algorithms compute and visualize all simulations that can generate results of the same or better quality than the optimum, subject to constraints. These optimization posteriors are desirable for a number of reasons among which are easy interpretability, automatic parameter sensitivity and correlation analysis, and posterior predictive analysis. Our algorithms are simple extensions to an existing simulation-based inference framework called approximate Bayesian computation. POPE is applied two biological simulators: a fast and stochastic simulator of stem-cell cycling and a slow and deterministic simulator of tumor growth patterns. Conclusions: POPE allows the scientist to explore and understand the role that constraints, both on the input and the output, have on the optimization posterior. As a Bayesian inference procedure, POPE provides a rigorous framework for the analysis of the uncertainty of an optimal simulation parameter setting.

Description: Background: Invasive species can be a major threat to native biodiversity and the number of invasive plant species is increasing across the globe. Population genetic studies of invasive species can provide key insights into their invasion history and ensuing evolution, but also for their control. Here we genetically characterise populations of Impatiens glandulifera, an invasive plant in Europe that can have a major impact on native plant communities. We compared populations from the species’ native range in Kashmir, India, to those in its invaded range, along a latitudinal gradient in Europe. For comparison, the results from 39 other studies of genetic diversity in invasive species were collated. Results: Our results suggest that I. glandulifera was established in the wild in Europe at least twice, from an area outside of our Kashmir study area. Our results further revealed that the genetic diversity in invasive populations of I. glandulifera is unusually low compared to native populations, in particular when compared to other invasive species. Genetic drift rather than mutation seems to have played a role in differentiating populations in Europe. We find evidence of limitations to local gene flow after introduction to Europe, but somewhat less restrictions in the native range. I. glandulifera populations with significant inbreeding were only found in the species’ native range and invasive species in general showed no increase in inbreeding upon leaving their native ranges. In Europe we detect cases of migration between distantly located populations. Human activities therefore seem to, at least partially, have facilitated not only introductions, but also further spread of I. glandulifera across Europe. Conclusions: Although multiple introductions will facilitate the retention of genetic diversity in invasive ranges, widespread invasive species can remain genetically relatively invariant also after multiple introductions. Phenotypic plasticity may therefore be an important component of the successful spread of Impatiens glandulifera across Europe.

Description: Background: In the male germline, neonatal prospermatogonia give rise to spermatogonia, which include stem cell population (undifferentiated spermatogonia) that supports continuous spermatogenesis in adults. Although the levels of DNA methyltransferases change dynamically in the neonatal and early postnatal male germ cells, detailed genome-wide DNA methylation profiles of these cells during the stem cell formation and differentiation have not been reported. Results: To understand the regulation of spermatogonial stem cell formation and differentiation, we examined the DNA methylation and gene expression dynamics of male mouse germ cells at the critical stages: neonatal prospermatogonia, and early postntal (day 7) undifferentiated and differentiating spermatogonia. We found large partially methylated domains similar to those found in cancer cells and placenta in all these germ cells, and high levels of non-CG methylation and 5-hydroxymethylcytosines in neonatal prospermatogonia. Although the global CG methylation levels were stable in early postnatal male germ cells, and despite the reported scarcity of differential methylation in the adult spermatogonial stem cells, we identified many regions showing stage-specific differential methylation in and around genes important for stem cell function and spermatogenesis. These regions contained binding sites for specific transcription factors including the SOX family members. Conclusions: Our findings show a distinctive and dynamic regulation of DNA methylation during spermatogonial stem cell formation and differentiation in the neonatal and early postnatal testes. Furthermore, we revealed a unique accumulation and distribution of non-CG methylation and 5hmC marks in neonatal prospermatogonia. These findings contrast with the reported scarcity of differential methylation in adult spermatogonial stem cell differentiation and represent a unique phase of male germ cell development.

Description: Background: Natural and artificial selection following domestication has led to the existence of more than a hundred pig breeds, as well as incredible variation in phenotypic traits. Berkshire pigs are regarded as having superior meat quality compared to other breeds. As the meat production industry seeks selective breeding approaches to improve profitable traits such as meat quality, information about genetic determinants of these traits is in high demand. However, most of the studies have been performed using trained sensory panel analysis without investigating the underlying genetic factors. Here we investigate the relationship between genomic composition and this phenotypic trait by scanning for signatures of positive selection in whole-genome sequencing data. Results: We generated genomes of 10 Berkshire pigs at a total of 100.6 coverage depth, using the Illumina Hiseq2000 platform. Along with the genomes of 11 Landrace and 13 Yorkshire pigs, we identified genomic variants of 18.9 million SNVs and 3.4 million Indels in the mapped regions. We identified several associated genes related to lipid metabolism, intramuscular fatty acid deposition, and muscle fiber type which attribute to pork quality (TG, FABP1, AKIRIN2, GLP2R, TGFBR3, JPH3, ICAM2, and ERN1) by applying between population statistical tests (XP-EHH and XP-CLR). A statistical enrichment test was also conducted to detect breed specific genetic variation. In addition, de novo short sequence read assembly strategy identified several candidate genes (SLC25A14, IGF1, PI4KA, CACNA1A) as also contributing to lipid metabolism. Conclusions: Results revealed several candidate genes involved in Berkshire meat quality; most of these genes are involved in lipid metabolism and intramuscular fat deposition. These results can provide a basis for future research on the genomic characteristics of Berkshire pigs.

Description: Background: The transposable element (TE) content of the genomes of plant species varies from near zero in the genome of Utricularia gibba to more than 80 % in many species. It is not well understood whether this variation in genome composition results from common mechanisms or stochastic variation. The major obstacles to investigating mechanisms of TE evolution have been a lack of comparative genomic data sets and efficient computational methods for measuring differences in TE composition between species. In this study, we describe patterns of TE evolution in 14 species in the flowering plant family Asteraceae and 1 outgroup species in the Calyceraceae to investigate phylogenetic patterns of TE dynamics in this important group of plants. Results: Our findings indicate that TE families in the Asteraceae exhibit distinct patterns of non-neutral evolution, and that there has been a directional increase in copy number of Gypsy retrotransposons since the origin of the Asteraceae. Specifically, there is marked increase in Gypsy abundance at the origin of the Asteraceae and at the base of the tribe Heliantheae. This latter shift in genome composition has had a significant impact on the diversity and abundance distribution of TEs in a lineage-specific manner. Conclusions: We show that the TE-driven expansion of plant genomes can be facilitated by just a few TE families, and is likely accompanied by the modification and/or replacement of the TE community. Importantly, large shifts in TE composition may be correlated with major of phylogenetic transitions.

Description: Background: Acute kidney injury (AKI) is a severe disease with high morbidity and mortality. Methods that promote repair of the injured kidney have been extensively investigated. Cell-based therapy with mesenchymal stem cells or renal progenitor cells (RPCs) resident in the kidney has appeared to be an effective strategy for the treatment of AKI. Embryonic stem cells or induced pluripotent stem cells (iPSCs) are also utilized for AKI recovery. However, the therapeutic effect of iPSC-derived RPCs for AKI has yet to be determined. Methods: In this study, we induced iPSCs differentiation into RPCs using a nephrogenic cocktail of factors combined with the renal epithelial cell growth medium. We then established the rat ischemia–reperfusion injury (IR) model and transplanted the iPSC-derived RPCs into the injured rats in combination with the hydrogel. Next, we examined the renal function-related markers and renal histology to assess the therapeutic effect of the injected cells. Moreover, we investigated the mechanism by which iPSC-derived RPCs affect AKI caused by IR. Results: We showed that the differentiation efficiency of iPSCs to RPCs increased when cultured with renal epithelial cell growth medium after stimulation with a nephrogenic cocktail of factors. The transplantation of iPSC-derived RPCs decreased the levels of biomarkers indicative of renal injury and attenuated the necrosis and apoptosis of renal tissues, but resulted in the up-regulation of renal tubules formation, cell proliferation, and the expression of pro-renal factors. Conclusion: Our results revealed that iPSC-derived RPCs can protect AKI rat from renal function impairment and severe tubular injury by up-regulating the renal tubules formation, promoting cell proliferation, reducing apoptosis, and regulating the microenvironment in the injured kidney.

Description: Background: Hospital accreditation programs are internationally widespread and consume increasingly scarce health resources. However, we lack tools to consistently identify suitable indicators to assess and monitor accreditation outcomes. We describe the development and validation of such a tool. Results: Using Australian accreditation standards as our reference point we: reviewed the research evidence for potential indicators; looked for links with existing external indicators; and assessed relevant state and federal policies. We allocated provisional scores, on a five point Likert scale, to the five accountability criteria in the tool: research; accuracy; proximity; no adverse effects; and specificity. An expert panel validated the use of the purpose designed indicator assessment tool. The panel identified hand hygiene compliance rates as a suitable process indicator, and hospital acquired Staphylococcus aureus infection (SAB) rates as an outcome indicator, with the hypothesis that improved hand hygiene compliance rates and lower SAB rates would correlate with accreditation performance. Conclusions: This new tool can be used to identify, analyse, and compare accreditation indicators. Using infection control indicators such as hand hygiene compliance and SAB rates to measure accreditation effectiveness has merit, and their efficacy can be determined by comparing accreditation scores with indicator outcomes. To verify the tool as a robust instrument, testing is needed in other health service domains, both in Australia and internationally. This tool provides health policy makers with an important means for assessing the accreditation programs which form a critical part of the national patient safety and quality framework.

Description: Background: Recently, a lot of energy devices in the surgical field, especially in the liver surgery, have been developed, and a fine tip LigaSure™, Dolphin Tip Sealer/Divider (DT-SD) also has been used frequently to dissect liver parenchyma as well as ultrasonically activated device (USAD). However, the utility of this instrument for liver dissection (LD) is still unknown. Moreover, to reduce bleeding during LD, a half-grip technique (HGT) was contrived. We herein report an experimental study in swine model to evaluate the feasibility and effectiveness of HGT using DT-SD for LD. Methods: The swine model experiment was carried out under general anesthesia by veterinarians. LD was performed repeatedly by DT-SD with the HGT (Group A, n = 6), or the conventional clamp-crush technique (CCT) (Group B, n = 6), and by variable mode USAD (Group C, n = 6). The dissection length and depth (cm) as well as bleeding volume (g) were measured carefully, and the dissection area (cm 2 ) and speed (cm 2 /min) were calculated precisely. Histological examinations of the dissection surfaces were also executed. Mann–Whitney’s U test was used for Statistical analyses with variance at a significance level of 0.05. Results: Among the three groups, the three averages of dissection lengths were unexpectedly equalized to 8.3 cm. The dissection area (cm 2 ) was 9.9 ± 5.1 in Group A, 9.8 ± 4.7 in Group B, and 9.9 ± 4.5 in Group C. The mean blood loss during LD was 10.6 ± 14.8 g in Group A, 41.4 ± 39.2 g in Group B, and 34.3 ± 39.2 g in Group C. For Group A, the bleeding rate was the least, 0.9 ± 1.0 g/cm 2 , and the average depth of coagulation was the thickest, 1.47 ± 0.29 mm, among the three groups (p 

Description: Background: Streptococcus pneumoniae Neuraminidase A (NanA) is a multi-domain protein anchored to the bacterial surface. Upstream of the catalytic domain of NanA is a domain that conforms to the sialic acid-recognising CBM40 family of the CAZY (carbohydrate-active enzymes) database. This domain has been identified to play a critical role in allowing the bacterium to promote adhesion and invasion of human brain microvascular endothelial cells, and hence may play a key role in promoting bacterial meningitis. In addition, the CBM40 domain has also been reported to activate host chemokines and neutrophil recruitment during infection. Results: Crystal structures of both apo- and holo- forms of the NanA CBM40 domain (residues 121 to 305), have been determined to 1.8 Å resolution. The domain shares the fold of other CBM40 domains that are associated with sialidases. When in complex with α2,3- or α2,6-sialyllactose, the domain is shown to interact only with the terminal sialic acid. Significantly, a deep acidic pocket adjacent to the sialic acid-binding site is identified, which is occupied by a lysine from a symmetry-related molecule in the crystal. This pocket is adjacent to a region that is predicted to be involved in protein-protein interactions. Conclusions: The structural data provide the details of linkage-independent sialyllactose binding by NanA CBM40 and reveal striking surface features that may hold the key to recognition of binding partners on the host cell surface. The structure also suggests that small molecules or sialic acid analogues could be developed to fill the acidic pocket and hence provide a new therapeutic avenue against meningitis caused by S. pneumoniae.

Description: Background: Body fluids such as saliva and tears from patients with hepatitis B virus (HBV) infection are known as infectious agents. The infectivity of feces from patients with HBV infection has not been established. The aim of this study was to determine whether feces from HBV carriers can be a source of HBV infection. Methods: Thirty-three children and 17 adults (ages 0–49 years, median age 13 years) who were chronically infected with HBV were enrolled. The levels of HBV DNA in the feces from these patients were quantified by real-time PCR, and the levels of fecal HBsAg were measured. Isolated human hepatocytes from chimeric mice with humanized livers were co-cultured with serum, tears and feces from the HBV carriers. Four chimeric mice were inoculated intravenously with sterilized feces from HBV carriers. Results: HBV DNA was detected in the feces of 37 (74 %) of the 50 patients. The fecal HBV DNA levels ranged from 2.8 to 8.4 log copies/mL (mean ± SD  =  5.6 ± 1.2 log copies/mL). A significant correlation was observed in the levels of HBV DNA between serum and feces (r  =  0.54, p 

Description: Background: Cetylpyridinium chloride (CPC) and sodium fluoride augment oral hygiene by inactivating bacteria and inhibiting enamel demineralisation, respectively. However, there are few reports in the literature documenting the antibacterial efficacy of their combined use in mouthrinses. We have used six experimental systems to compare the antibacterial effects of mouthrinses containing 0.075 % CPC (test rinse, TR) or 0.075 % CPC with sodium fluoride (test fluoride rinse, TFR). Results: Effects against planktonic bacteria were determined using viable counting (for Streptococcus mutans and salivary bacteria), a redox dye (for Actinomyces viscosus and salivary bacteria) and viable counting (for ex vivo oral rinses). Effects against saliva-derived biofilms were quantified using confocal microscopy and differential viable counting. Inhibition of biofilm formation was evaluated by pre-treating hydroxyapatite coupons with mouthrinses prior to inoculation. Otherwise-identical controls without CPC (control rinse and control fluoride rinse, CR and CFR, respectively), were included throughout. Compared to the controls, TFR and TR demonstrated significant antimicrobial effects in the redox assays, by viable counts (〉3 log reductions) and in oral rinse samples (〉1.25 log reductions, p  3 log difference, p 

Description: Background: Root rot caused by Aphanomyces euteiches is one of the most destructive pea diseases while a distantly related species P. pisi has been recently described as the agent of pea and faba bean root rot. These two oomycete pathogens with different pathogenicity factor repertories have both evolved specific mechanisms to infect pea. However, little is known about the genes and mechanisms of defence against these pathogens in pea. In the present study, the transcriptomic response of pea to these two pathogens was investigated at two time points during early phase of infection using a Medicago truncatula microarray. Results: Of the 37,976 genes analysed, 574 and 817 were differentially expressed in response to A. euteiches at 6 hpi and 20 hpi, respectively, while 544 and 611 genes were differentially regulated against P. pisi at 6 hpi and 20 hpi, respectively. Differentially expressed genes associated with plant immunity responses were involved in cell wall reinforcement, hormonal signalling and phenylpropanoid metabolism. Activation of cell wall modification, regulation of jasmonic acid biosynthesis and induction of ethylene signalling pathway were among the common transcriptional responses to both of these oomycetes. However, induction of chalcone synthesis and the auxin pathway were specific transcriptional changes against A. euteiches. Conclusions: Our results demonstrate a global view of differentially expressed pea genes during compatible interactions with P. pisi and A. euteiches at an early phase of infection. The results suggest that distinct signalling pathways are triggered in pea by these two pathogens that lead to common and specific immune mechanisms in response to these two oomycetes. The generated knowledge may eventually be used in breeding pea varieties with resistance against root rot disease.

Description: Background: There is general consensus that appropriate development and use of information and communication technologies (ICT) are crucial in the delivery of effective primary care (PC). Several countries are defining policies to support and promote a structural change of the health care system through the introduction of ICT. This study analyses the state of development of basic ICT in PC systems of 31 European countries with the aim to describe the extent of, and main purposes for, computer use by General Practitioners (GPs) across Europe. Additionally, trends over time have been analysed. Methods: Descriptive statistical analysis was performed on data from the QUALICOPC (Quality and Costs of Primary Care in Europe) survey, to describe the geographic differences in the general use of computer, and in specific computerized clinical functions for different health-related purposes such as prescribing, medication checking, generating health records and research for medical information on the Internet. Results: While all the countries have achieved a near-universal adoption of a computer in their primary care practices, with only a few countries near or under the boundary of 90 %, the computerisation of primary care clinical functions presents a wide variability of adoption within and among countries and, in several cases (such as in the southern and central-eastern Europe), a large room for improvement. Conclusions: At European level, more efforts could be done to support southern and central-eastern Europe in closing the gap in adoption and use of ICT in PC. In particular, more attention seems to be need on the current usages of the computer in PC, by focusing policies and actions on the improvement of the appropriate usages that can impact on quality and costs of PC and can facilitate an interconnected health care system. However, policies and investments seem necessary but not sufficient to achieve these goals. Organizational, behavioural and also networking aspects should be taken in consideration.

Description: Background: In humans it is unknown if the composition of the gut microbiota alters the risk of Plasmodium falciparum infection or the risk of developing febrile malaria once P. falciparum infection is established. Here we collected stool samples from a cohort composed of 195 Malian children and adults just prior to an intense P. falciparum transmission season. We assayed these samples using massively parallel sequencing of the 16S ribosomal RNA gene to identify the composition of the gut bacterial communities in these individuals. During the ensuing 6-month P. falciparum transmission season we examined the relationship between the stool microbiota composition of individuals in this cohort and their prospective risk of both P. falciparum infection and febrile malaria. Results: Consistent with prior studies, stool microbial diversity in the present cohort increased with age, although the overall microbiota profile was distinct from cohorts in other regions of Africa, Asia and North America. Age-adjusted Cox regression analysis revealed a significant association between microbiota composition and the prospective risk of P. falciparum infection; however, no relationship was observed between microbiota composition and the risk of developing febrile malaria once P. falciparum infection was established. Conclusions: These findings underscore the diversity of gut microbiota across geographic regions, and suggest that strategic modulation of gut microbiota composition could decrease the risk of P. falciparum infection in malaria-endemic areas, potentially as an adjunct to partially effective malaria vaccines.

Description: Background: It remains unclear to what extent midgut rotation determines human intestinal topography and pathology. We reinvestigated the midgut during its looping and herniation phases of development, using novel 3D visualization techniques. Results: We distinguished 3 generations of midgut loops. The topography of primary and secondary loops was constant, but that of tertiary loops not. The orientation of the primary loop changed from sagittal to transverse due to the descent of ventral structures in a body with a still helical body axis. The 1 st secondary loop (duodenum, proximal jejunum) developed intraabdominally towards a left-sided position. The 2 nd secondary loop (distal jejunum) assumed a left-sided position inside the hernia before returning, while the 3 rd and 4 th secondary loops retained near-midline positions. Intestinal return into the abdomen resembled a backward sliding movement. Only after return, the 4 th secondary loop (distal ileum, cecum) rapidly “slid” into the right lower abdomen. The seemingly random position of the tertiary small-intestinal loops may have a biomechanical origin. Conclusions: The interpretation of “intestinal rotation” as a mechanistic rather than a descriptive concept underlies much of the confusion accompanying the physiological herniation. We argue, instead, that the concept of “en-bloc rotation” of the developing midgut is a fallacy of schematic drawings. Primary, secondary and tertiary loops arise in a hierarchical fashion. The predictable position and growth of secondary loops is pre-patterned and determines adult intestinal topography. We hypothesize based on published accounts that malrotations result from stunted development of secondary loops.

Description: Background: The hemibiotroph Colletotrichum lentis, causative agent of anthracnose on Lens culinaris (lentil) was recently described as a new species. During its interaction with the host plant, C. lentis likely secretes numerous effector proteins, including toxins to alter the plant’s innate immunity, thereby gaining access to the host tissues for nutrition and reproduction. Results: In silico analysis of 2000 ESTs generated from C. lentis-infected lentil leaf tissues identified 15 candidate effectors. In planta infection stage-specific gene expression waves among candidate effectors were revealed for the appressorial penetration phase, biotrophic phase and necrotrophic phase. No sign of positive selection pressure [ω (dN/dS) 

Description: Background: Improved understanding of bacterial-fungal interactions in the rhizosphere should assist in the successful application of bacteria as biological control agents against fungal pathogens of plants, providing alternatives to chemicals in sustainable agriculture. Rhizoctonia solani is an important soil-associated fungal pathogen and its chemical treatment is not feasible or economic. The genomes of the plant-associated bacteria Serratia proteamaculans S4 and Serratia plymuthica AS13 have been sequenced, revealing genetic traits that may explain their diverse plant growth promoting activities and antagonistic interactions with R. solani. To understand the functional response of this pathogen to different bacteria and to elucidate whether the molecular mechanisms that the fungus exploits involve general stress or more specific responses, we performed a global transcriptome profiling of R. solani Rhs1AP anastomosis group 3 (AG-3) during interaction with the S4 and AS13 species of Serratia using RNA-seq. Results: Approximately 104,504 million clean 75-100 bp paired-end reads were obtained from three libraries, each in triplicate (AG3-Control, AG3-S4 and AG3-AS13). Transcriptome analysis revealed that approximately 10 % of the fungal transcriptome was differentially expressed during challenge with Serratia. The numbers of S4- and AS13-specific differentially expressed genes (DEG) were 866 and 292 respectively, while there were 1035 common DEGs in the two treatment groups. Four hundred and sixty and 242 genes respectively had values of log 2 fold-change 〉 3 and for further analyses this cut-off value was used. Functional classification of DEGs based on Gene Ontology enrichment analysis and on KEGG pathway annotations revealed a general shift in fungal gene expression in which genes related to xenobiotic degradation, toxin and antioxidant production, energy, carbohydrate and lipid metabolism and hyphal rearrangements were subjected to transcriptional regulation. Conclusions: This RNA-seq profiling generated a novel dataset describing the functional response of the phytopathogen R. solani AG3 to the plant-associated Serratia bacteria S4 and AS13. Most genes were regulated in the same way in the presence of both bacterial isolates, but there were also some strain-specific responses. The findings in this study will be beneficial for further research on biological control and in depth exploration of bacterial-fungal interactions in the rhizosphere.

Description: Background: Hemolysins are crucial virulence factors which help pathogens to survive and persist in the host. This study investigated whether common electrolytes will affect the hemolysins of Candida species. The hemolysins from 25 Candida isolates were investigated using a plate specially designed for Candida species in the presence of three electrolytes, CaCl 2 , NaCl and KCl, at different concentrations. The hemolytic activity was determined after 48 h and the hemolytic index was calculated. Results: All three electrolytes caused a decrease in the hemolytic activity. Significant differences existed between varying concentrations of NaCl, while no significant differences existed for the CaCl 2 and KCl groups. Additionally, the peripheral hemolytic index was highly correlated with the hemolytic index (r = 0.656, p 

Description: Background: The radioactive isotope Strontium-90 ( 90 Sr) may be released as a component of fallout from nuclear accidents, or in the event of a radiological incident such as detonation of an improvised nuclear device, and if ingested poses a significant health risk to exposed individuals. In order to better understand the response to 90 Sr, using an easily attainable and standard biodosimetry sample fluid, we analyzed the global transcriptomic response of blood cells in an in vivo model system. Results: We injected C57BL/6 mice with a solution of 90SrCl2 and followed them over a 30-day period. At days 4, 7, 9, 25 and 30, we collected blood and isolated RNA for microarray analyses. These days corresponded to target doses in a range from 1–5 Gy. We investigated changes in mRNA levels using microarrays, and changes in specific microRNA (miRNA) predicted to be involved in the response using qRT-PCR. We identified 8082 differentially expressed genes in the blood of mice exposed to 90 Sr compared with controls. Common biological functions were affected throughout the study, including apoptosis of B and T lymphocytes, and atrophy of lymphoid organs. Cellular functions such as RNA degradation and lipid metabolism were also affected during the study. The broad down regulation of genes observed in our study suggested a potential role for miRNA in gene regulation. We tested candidate miRNAs, mmu-miR-16, mmu-miR-124, mmu-miR-125 and mmu-mir-21; and found that all were induced at the earliest time point, day 4. Conclusions: Our study is the first to report the transcriptomic response of blood cells to the internal emitter 90 Sr in mouse and a possible role for microRNA in gene regulation after 90 Sr exposure. The most dramatic effect was observed on gene expression related to B-cell development and RNA maintenance. These functions were affected by genes that were down regulated throughout the study, suggesting severely compromised antigen response, which may be a result of the deposition of the radioisotope proximal to the hematopoietic compartment in bone.

Description: Background: A complete genome sequence is an essential tool for the genetic improvement of wheat. Because the wheat genome is large, highly repetitive and complex due to its allohexaploid nature, the International Wheat Genome Sequencing Consortium (IWGSC) chose a strategy that involves constructing bacterial artificial chromosome (BAC)-based physical maps of individual chromosomes and performing BAC-by-BAC sequencing. Here, we report the construction of a physical map of chromosome 6B with the goal of revealing the structural features of the third largest chromosome in wheat. Results: We assembled 689 informative BAC contigs (hereafter reffered to as contigs) representing 91 % of the entire physical length of wheat chromosome 6B. The contigs were integrated into a radiation hybrid (RH) map of chromosome 6B, with one linkage group consisting of 448 loci with 653 markers. The order and direction of 480 contigs, corresponding to 87 % of the total length of 6B, were determined. We also characterized the contigs that contained a part of the nucleolus organizer region or centromere based on their positions on the RH map and the assembled BAC clone sequences. Analysis of the virtual gene order along 6B using the information collected for the integrated map revealed the presence of several chromosomal rearrangements, indicating evolutionary events that occurred on chromosome 6B. Conclusions: We constructed a reliable physical map of chromosome 6B, enabling us to analyze its genomic structure and evolutionary progression. More importantly, the physical map should provide a high-quality and map-based reference sequence that will serve as a resource for wheat chromosome 6B.

Description: Background: Base-resolution methylome data generated by whole-genome bisulfite sequencing (WGBS) is often used to segment the genome into domains with distinct methylation levels. However, most segmentation methods include many parameters to be carefully tuned and/or fail to exploit the unsurpassed resolution of the data. Furthermore, there is no simple method that displays the composition of the domains to grasp global trends in each methylome. Results: We propose to use changepoint detection for domain demarcation based on base-resolution methylome data. While the proposed method segments the methylome in a largely comparable manner to conventional approaches, it has only a single parameter to be tuned. Furthermore, it fully exploits the base-resolution of the data to enable simultaneous detection of methylation changes in even contrasting size ranges, such as focal hypermethylation and global hypomethylation in cancer methylomes. We also propose a simple plot termed methylated domain landscape (MDL) that globally displays the size, the methylation level and the number of the domains thus defined, thereby enabling one to intuitively grasp trends in each methylome. Since the pattern of MDL often reflects cell lineages and is largely unaffected by data size, it can serve as a novel signature of methylome. Conclusions: Changepoint detection in base-resolution methylome data followed by MDL plotting provides a novel method for methylome characterization and will facilitate global comparison among various WGBS data differing in size and even species origin.

Description: Background: Bi-parental mapping populations have been commonly utilized to identify and characterize quantitative trait loci (QTL) controlling resistance to soybean cyst nematode (SCN, Heterodera glycines Ichinohe). Although this approach successfully mapped a large number of SCN resistance QTL, it captures only limited allelic diversity that exists in parental lines, and it also has limitations for genomic resolution. In this study, a genome-wide association study (GWAS) was performed using a diverse set of 553 soybean plant introductions (PIs) belonging to maturity groups from III to V to detect QTL/genes associated with SCN resistance to HG Type 0. Results: Over 45,000 single nucleotide polymorphism (SNP) markers generated by the SoySNP50K iSelect BeadChip (http// www.soybase.org ) were utilized for analysis. GWAS identified 14 loci distributed over different chromosomes comprising 60 SNPs significantly associated with SCN resistance. Results also confirmed six QTL that were previously mapped using bi-parental populations, including the rhg1 and Rhg4 loci. GWAS identified eight novel QTL, including QTL on chromosome 10, which we have previously mapped by using a bi-parental population. In addition to the known loci for four simple traits, such as seed coat color, flower color, pubescence color, and stem growth habit, two traits, like lodging and pod shattering, having moderately complex inheritance have been confirmed with great precision by GWAS. Conclusions: The study showed that GWAS can be employed as an effective strategy for identifying complex traits in soybean and for narrowing GWAS-defined genomic regions, which facilitates positional cloning of the causal gene(s).

Description: Background: Unlike mammals, zebrafish have the ability to regenerate damaged parts of their central nervous system (CNS) and regain functionality of the affected area. A better understanding of the molecular mechanisms involved in zebrafish regeneration may therefore provide insight into how CNS repair might be induced in mammals. Although many studies have described differences in gene expression in zebrafish during CNS regeneration, the regulatory mechanisms underpinning the differential expression of these genes have not been examined. Results: We used microarrays to analyse and integrate the mRNA and microRNA (miRNA) expression profiles of zebrafish retina after optic nerve crush to identify potential regulatory mechanisms that underpin central nerve regeneration. Bioinformatic analysis identified 3 miRNAs and 657 mRNAs that were differentially expressed after injury. We then combined inverse correlations between our miRNA expression and mRNA expression, and integrated these findings with target predictions from TargetScan Fish to identify putative miRNA-gene target pairs. We focused on two over-expressed miRNAs (miR-29b and miR-223), and functionally validated seven of their predicted gene targets using RT-qPCR and luciferase assays to confirm miRNA-mRNA binding. Gene ontology analysis placed the miRNA-regulated genes (eva1a, layna, nefmb, ina, si:ch211-51a6.2, smoc1, sb:cb252) in key biological processes that included cell survival/apoptosis, ECM-cytoskeleton signaling, and heparan sulfate proteoglycan binding, Conclusion: Our results suggest a key role for miR-29b and miR-223 in zebrafish regeneration. The identification of miRNA regulation in a zebrafish injury model provides a framework for future studies in which to investigate not only the cellular processes required for CNS regeneration, but also how these mechanisms might be regulated to promote successful repair and return of function in the injured mammalian brain.

Description: Background: Ameloblast differentiation is the most critical stepwise process in amelogenesis, and it is controlled by precise molecular events. To better understand the mechanism controlling pre-ameloblasts (PABs) differentiation into secretory ameloblasts (SABs), a more precise identification of molecules and signaling networks will elucidate the mechanisms governing enamel formation and lay a foundation for enamel regeneration. Results: We analyzed transcriptional profiles of human PABs and SABs. From a total of 28,869 analyzed transcripts, we identified 923 differentially expressed genes (DEGs) with p   2. Among the DEGs, 647 genes showed elevated expression in PABs compared to SABs. Notably, 38 DEGs displayed greater than eight-fold changes. Comparative analysis revealed that highly expressed genes in PABs were involved in cell cycle control, DNA damage repair and apoptosis, while highly expressed genes in SABs were related to cell adhesion and extracellular matrix. Moreover, coexpression network analysis uncovered two highly conserved sub-networks contributing to differentiation, containing transcription regulators (RUNX2, ETV1 and ETV5), solute carrier family members (SLC15A1 and SLC7A11), enamel matrix protein (MMP20), and a polymodal excitatory ion channel (TRPA1). Conclusions: By combining comparative analysis and coexpression networks, this study provides novel biomarkers and research targets for ameloblast differentiation and the potential for their application in enamel regeneration.

Description: Background: Estimating the phylogenetic position of bacterial and archaeal organisms by genetic sequence comparisons is considered as the gold-standard in taxonomy. This is also a way to identify the species of origin of the sequence. The quality of the reference database used in such analyses is crucial: the database must reflect the up-to-date bacterial nomenclature and accurately indicate the species of origin of its sequences.DescriptionleBIBI QBPP is a web tool taking as input a series of nucleotide sequences belonging to one of a set of reference markers (e.g., SSU rRNA, rpoB, groEL2) and automatically retrieving closely related sequences, aligning them, and performing phylogenetic reconstruction using an approximate maximum likelihood approach. The system returns a set of quality parameters and, if possible, a suggested taxonomic assigment for the input sequences. The reference databases are extracted from GenBank and present four degrees of stringency, from the “superstringent” degree (one type strain per species) to the loosely parsed degree (“lax” database). A set of one hundred to more than a thousand sequences may be analyzed at a time. The speed of the process has been optimized through careful hardware selection and database design. Conclusion: leBIBI QBPP is a powerful tool helping biologists to position bacterial or archaeal sequence commonly used markers in a phylogeny. It is a diagnostic tool for clinical, industrial and environmental microbiology laboratory, as well as an exploratory tool for more specialized laboratories. Its main advantages, relatively to comparable systems are: i) the use of a broad set of databases covering diverse markers with various degrees of stringency; ii) the use of an approximate Maximum Likelihood approach for phylogenetic reconstruction; iii) a speed compatible with on-line usage; and iv) providing fully documented results to help the user in decision making.

Description: Background: Ossifying metaplasia is an unusual feature of urothelial carcinoma, with only a few cases reported. The largest series included 17 cases and was published in 1991. The mechanism of ossification is unknown and hypotheses of osteogenic precursor cells, inducing bone formation, are proposed.Case presentationA 75 year-old patient was treated for a high grade transitional cell carcinoma of the bladder by surgery, chemotherapy and radiotherapy. Histology showed foci of bone metaplasia, both at the periphery of the tumor, and in a lymph node metastasis. 1 year later, a heterotopic bone formation was discovered in the right retroperitoneal space, near the lumbar spine, increasing rapidly in size during follow-up. Several imaging exams were performed (2 CT, 1 MRI, 1 Pet-CT), but in the absence of typical features of sarcoma, diagnosis remained unclear.Histology of a CT-guided percutaneous biopsy showed urothelial carcinoma and mature lamellar bone. Integration of these findings with the radiological description of extraosseous localization was consistent with a diagnosis of osseous metaplasia of an urothelial carcinoma metastasis. The absence of bone atypia in both the primary and metastases argues against sarcomatoid urothelial carcinoma with osteosarcomatous differentiation. Conclusion: Osseous metaplasia of an urothelial carcinoma metastasis is unusual, and difficult to distinguish from radiotherapy induced sarcoma, or from sarcomatoid carcinoma. Rapid progression, sheathing of adjacent structures such as vessels (like inferior vena cava in our case) and nerves and bony feature of lymph node metastases necessitate histological confirmation and rapid treatment. Our case illustrates this disease and evaluates the imaging features. In addition we discuss the differential diagnosis of osseous retroperitoneal masses.

Description: Background: Herbicide resistance in agrestal weeds is a global problem threatening food security. Non-target-site resistance (NTSR) endowed by mechanisms neutralising the herbicide or compensating for its action is considered the most agronomically noxious type of resistance. Contrary to target-site resistance, NTSR mechanisms are far from being fully elucidated. A part of weed response to herbicide stress, NTSR is considered to be largely driven by gene regulation. Our purpose was to establish a transcriptome resource allowing investigation of the transcriptomic bases of NTSR in the major grass weed Alopecurus myosuroides L. (Poaceae) for which almost no genomic or transcriptomic data was available. Results: RNA-Seq was performed from plants in one F2 population that were sensitive or expressing NTSR to herbicides inhibiting acetolactate-synthase. Cloned plants were sampled over seven time-points ranging from before until 73 h after herbicide application. Assembly of over 159M high-quality Illumina reads generated a transcriptomic resource (ALOMYbase) containing 65,558 potentially active contigs (N50 = 1240 nucleotides) predicted to encode 32,138 peptides with 74 % GO annotation, of which 2017 were assigned to protein families presumably involved in NTSR. Comparison with the fully sequenced grass genomes indicated good coverage and correct representation of A. myosuroides transcriptome in ALOMYbase. The part of the herbicide transcriptomic response common to the resistant and the sensitive plants was consistent with the expected effects of acetolactate-synthase inhibition, with striking similarities observed with published Arabidopsis thaliana data. A. myosuroides plants with NTSR were first affected by herbicide action like sensitive plants, but ultimately overcame it. Analysis of differences in transcriptomic herbicide response between resistant and sensitive plants did not allow identification of processes directly explaining NTSR. Five contigs associated to NTSR in the F2 population studied were tentatively identified. They were predicted to encode three cytochromes P450 (CYP71A, CYP71B and CYP81D), one peroxidase and one disease resistance protein. Conclusions: Our data confirmed that gene regulation is at the root of herbicide response and of NTSR. ALOMYbase proved to be a relevant resource to support NTSR transcriptomic studies, and constitutes a valuable tool for future research aiming at elucidating gene regulations involved in NTSR in A. myosuroides.

Description: Background: Pineapple (Ananas comosus var. comosus) is known as the king of fruits for its crown and is the third most important tropical fruit after banana and citrus. The plant, which is indigenous to South America, is the most important species in the Bromeliaceae family and is largely traded for fresh fruit consumption. Here, we report the complete chloroplast sequence of the MD-2 pineapple that was sequenced using the PacBio sequencing technology. Results: In this study, the high error rate of PacBio long sequence reads of A. comosus’s total genomic DNA were improved by leveraging on the high accuracy but short Illumina reads for error-correction via the latest error correction module from Novocraft. Error corrected long PacBio reads were assembled by using a single tool to produce a contig representing the pineapple chloroplast genome. The genome of 159,636 bp in length is featured with the conserved quadripartite structure of chloroplast containing a large single copy region (LSC) with a size of 87,482 bp, a small single copy region (SSC) with a size of 18,622 bp and two inverted repeat regions (IRA and IRB) each with the size of 26,766 bp. Overall, the genome contained 117 unique coding regions and 30 were repeated in the IR region with its genes contents, structure and arrangement similar to its sister taxon, Typha latifolia. A total of 35 repeats structure were detected in both the coding and non-coding regions with a majority being tandem repeats. In addition, 205 SSRs were detected in the genome with six protein-coding genes contained more than two SSRs. Comparative chloroplast genomes from the subclass Commelinidae revealed a conservative protein coding gene albeit located in a highly divergence region. Analysis of selection pressure on protein-coding genes using Ka/Ks ratio showed significant positive selection exerted on the rps7 gene of the pineapple chloroplast with P less than 0.05. Phylogenetic analysis confirmed the recent taxonomical relation among the member of commelinids which support the monophyly relationship between Arecales and Dasypogonaceae and between Zingiberales to the Poales, which includes the A. comosus. Conclusions: The complete sequence of the chloroplast of pineapple provides insights to the divergence of genic chloroplast sequences from the members of the subclass Commelinidae. The complete pineapple chloroplast will serve as a reference for in-depth taxonomical studies in the Bromeliaceae family when more species under the family are sequenced in the future. The genetic sequence information will also make feasible other molecular applications of the pineapple chloroplast for plant genetic improvement.

Description: Background: Plant growth-promoting rhizobacteria are increasingly being seen as a way of complementing conventional inputs in agricultural systems. The effects on their host plants are diverse and include volatile-mediated growth enhancement. This study sought to assess the effects of bacterial volatiles on the biomass production and root system architecture of the model grass Brachypodium distachyon (L.) Beauv. Results: An in vitro experiment allowing plant-bacteria interaction throughout the gaseous phase without any physical contact was used to screen 19 bacterial strains for their growth-promotion ability over a 10-day co-cultivation period. Five groups of bacteria were defined and characterised based on their combined influence on biomass production and root system architecture. The observed effects ranged from unchanged to greatly increased biomass production coupled with increased root length and branching. Primary root length was increased only by the volatile compounds emitted by Enterobacter cloacae JM22 and Bacillus pumilus T4. Overall, the most significant results were obtained with Bacillus subtilis GB03, which induced an 81 % increase in total biomass, as well as enhancing total root length, total secondary root length and total adventitious root length by 88.5, 201.5 and 474.5 %, respectively. Conclusions: This study is the first report on bacterial volatile-mediated growth promotion of a grass plant. Contrasting modulations of biomass production coupled with changes in root system architecture were observed. Most of the strains that increased total plant biomass also modulated adventitious root growth. Under our screening conditions, total biomass production was strongly correlated with the length and branching of the root system components, except for primary root length. An analysis of the emission kinetics of the bacterial volatile compounds is being undertaken and should lead to the identification of the compounds responsible for the observed growth-promotion effects. Within the context of the inherent characteristics of our in vitro system, this paper identifies the next critical experimental steps and discusses them from both a fundamental and an applied perspective.

Description: Background: Species showing complex life cycles provide excellent opportunities to study the genetic associations between life cycle stages, as selective pressures may differ before and after metamorphosis. The European eel presents a complex life cycle with two metamorphoses, a first metamorphosis from larvae into glass eels (juvenile stage) and a second metamorphosis into silver eels (adult stage). We tested the hypothesis that different genes and gene pathways will be under selection at different life stages when comparing the genetic associations between glass eels and silver eels. Results: We used two sets of markers to test for selection: first, we genotyped individuals using a panel of 80 coding-gene single nucleotide polymorphisms (SNPs) developed in American eel; second, we investigated selection at the genome level using a total of 153,423 RAD-sequencing generated SNPs widely distributed across the genome. Using the RAD approach, outlier tests identified a total of 2413 (1.57 %) potentially selected SNPs. Functional annotation analysis identified signal transduction pathways as the most over-represented group of genes, including MAPK/Erk signalling, calcium signalling and GnRH (gonadotropin-releasing hormone) signalling. Many of the over-represented pathways were related to growth, while others could result from the different conditions that eels inhabit during their life cycle. Conclusions: The observation of different genes and gene pathways under selection when comparing glass eels vs. silver eels supports the adaptive decoupling hypothesis for the benefits of metamorphosis. Partitioning the life cycle into discrete morphological phases may be overall beneficial since it allows the different life stages to respond independently to their unique selection pressures. This might translate into a more effective use of food and niche resources and/or performance of phase-specific tasks (e.g. feeding in the case of glass eels, migrating and reproducing in the case of silver eels).

Description: Background: Avocado (Persea americana) is an economically important tropical fruit considered to be a good source of fatty acids. Despite its importance, the molecular and cellular characterization of biochemical and developmental processes in avocado is limited due to the lack of transcriptome and genomic information. Results: The transcriptomes of seeds, roots, stems, leaves, aerial buds and flowers were determined using different sequencing platforms. Additionally, the transcriptomes of three different stages of fruit ripening (pre-climacteric, climacteric and post-climacteric) were also analyzed. The analysis of the RNAseqatlas presented here reveals strong differences in gene expression patterns between different organs, especially between root and flower, but also reveals similarities among the gene expression patterns in other organs, such as stem, leaves and aerial buds (vegetative organs) or seed and fruit (storage organs). Important regulators, functional categories, and differentially expressed genes involved in avocado fruit ripening were identified. Additionally, to demonstrate the utility of the avocado gene expression atlas, we investigated the expression patterns of genes implicated in fatty acid metabolism and fruit ripening. Conclusions: A description of transcriptomic changes occurring during fruit ripening was obtained in Mexican avocado, contributing to a dynamic view of the expression patterns of genes involved in fatty acid biosynthesis and the fruit ripening process.

Description: Background: Deidentified newborn screening bloodspot samples (NBS) represent a valuable potential resource for genomic research if impediments to whole exome sequencing of NBS deoxyribonucleic acid (DNA), including the small amount of genomic DNA in NBS material, can be overcome. For instance, genomic analysis of NBS could be used to define allele frequencies of disease-associated variants in local populations, or to conduct prospective or retrospective studies relating genomic variation to disease emergence in pediatric populations over time. In this study, we compared the recovery of variant calls from exome sequences of amplified NBS genomic DNA to variant calls from exome sequencing of non-amplified NBS DNA from the same individuals. Results: Using a standard alignment-based Genome Analysis Toolkit (GATK), we find 62,000–76,000 additional variants in amplified samples. After application of a unique kmer enumeration and variant detection method (RUFUS), only 38,000–47,000 additional variants are observed in amplified gDNA. This result suggests that roughly half of the amplification-introduced variants identified using GATK may be the result of mapping errors and read misalignment. Conclusions: Our results show that it is possible to obtain informative, high-quality data from exome analysis of whole genome amplified NBS with the important caveat that different data generation and analysis methods can affect variant detection accuracy, and the concordance of variant calls in whole-genome amplified and non-amplified exomes.

Description: Background: Non-ribosomal peptide synthetase (NRPS) enzymes govern the assembly of amino acids and related monomers into peptide-like natural products. A key goal of the field is to develop methods to effective recombine NRPS domains or modules, and thereby generate modified or entirely novel products. We previously showed that substitution of the condensation (C) and adenylation (A) domains in module 2 of the pyoverdine synthetase PvdD from Pseudomonas aeruginosa led to synthesis of modified pyoverdines in a minority of cases, but that more often the recombinant enzymes were non-functional. One possible explanation was that the majority of introduced C domains were unable to effectively communicate with the thiolation (T) domain immediately upstream, in the first module of PvdD. Results: To test this we first compared the effectiveness of C-A domain substitution relative to T-C-A domain substitution using three different paired sets of domains. Having previously demonstrated that the PvdD A/T domain interfaces are tolerant of domain substitution, we hypothesised that T-C-A domain substitution would lead to more functional recombinant enzymes, by maintaining native T/C domain interactions. Although we successfully generated two recombinant pyoverdines, having a serine or a N5-formyl-N5-hydroxyornithine residue in place of the terminal threonine of wild type pyoverdine, in neither case did the T-C-A domain substitution strategy lead to substantially higher product yield. To more comprehensively examine the abilities of non-native T domains to communicate effectively with the C domain of PvdD module 2 we then substituted the module 1 T domain with 18 different T domains sourced from other pyoverdine NRPS enzymes. In 15/18 cases the recombinant NRPS was functional, including 6/6 cases where the introduced T domain was located upstream of a C domain in its native context. Conclusions: Our data indicate that T domains are generally able to interact effectively with non-native C domains, contrasting with previous findings that they are not generally portable upstream of epimerisation (E) or thioesterase (TE) domains. This offers promise for NRPS recombination efforts, but also raises the possibility that some C domains are unable to efficiently accept non-native peptides at their donor site due to steric constraints or other limitations.

Description: Background: Acanthus is a unique genus consisting of both true mangrove and terrestrial species; thus, it represents an ideal system for studying the origin and adaptive evolution of mangrove plants to intertidal environments. However, little is known regarding the two respects of mangrove species in Acanthus. In this study, we sequenced the transcriptomes of the pooled roots and leaves tissues for a mangrove species, Acanthus ilicifolius, and its terrestrial congener, A. leucostachyus, to illustrate the origin of the mangrove species in this genus and their adaptive evolution to harsh habitats. Results: We obtained 73,039 and 69,580 contigs with N50 values of 741 and 1557 bp for A. ilicifolius and A. leucostachyus, respectively. Phylogenetic analyses based on four nuclear segments and three chloroplast fragments revealed that mangroves and terrestrial species in Acanthus fell into different clades, indicating a single origin of the mangrove species in Acanthus. Based on 6634 orthologs, A. ilicifolius and A. leucostachyus were found to be highly divergent, with a peak of synonymous substitution rate (Ks) distribution of 0.145 and an estimated divergence time of approximately 16.8 million years ago (MYA). The transgression in the Early to Middle Miocene may be the major reason for the entry of the mangrove lineage of Acanthus into intertidal environments. Gene ontology (GO) classifications of the full transcriptomes did not show any apparent differences between A. ilicifolius and A. leucostachyus, suggesting the absence of gene components specific to the mangrove transcriptomes. A total of 99 genes in A. ilicifolius were identified with signals of positive selection. Twenty-three of the 99 positively selected genes (PSGs) were found to be involved in salt, heat and ultraviolet stress tolerance, seed germination and embryo development under periodic inundation. These stress-tolerance related PSGs may be crucial for the adaptation of the mangrove species in this genus to stressful marine environments and may contribute to speciation in Acanthus. Conclusions: We characterized the transcriptomes of one mangrove species of Acanthus, A. ilicifolius, and its terrestrial relative, A. leucostachyus, and provided insights into the origin of the mangrove Acanthus species and their adaptive evolution to abiotic stresses in intertidal environments.

Description: Background: Visualization and analysis of molecular profiling data together with biological networks are able to provide new mechanistic insights into biological functions. Currently, it is possible to visualize high-throughput data on top of pre-defined network layouts, but they are not always adapted to a given data analysis task. A network layout based simultaneously on the network structure and the associated multidimensional data might be advantageous for data visualization and analysis in some cases. Results: We developed a Cytoscape app, which allows constructing biological network layouts based on the data from molecular profiles imported as values of node attributes. DeDaL is a Cytoscape 3 app, which uses linear and non-linear algorithms of dimension reduction to produce data-driven network layouts based on multidimensional data (typically gene expression). DeDaL implements several data pre-processing and layout post-processing steps such as continuous morphing between two arbitrary network layouts and aligning one network layout with respect to another one by rotating and mirroring. The combination of all these functionalities facilitates the creation of insightful network layouts representing both structural network features and correlation patterns in multivariate data. We demonstrate the added value of applying DeDaL in several practical applications, including an example of a large protein-protein interaction network. Conclusions: DeDaL is a convenient tool for applying data dimensionality reduction methods and for designing insightful data displays based on data-driven layouts of biological networks, built within Cytoscape environment. DeDaL is freely available for downloading at http://bioinfo-out.curie.fr/projects/dedal/.

Description: Background: All organisms must synthesize the enzymatic cofactor coenzyme A (CoA) from the precursor pantothenate. Most bacteria can synthesize pantothenate de novo by the condensation of pantoate and β-alanine. The synthesis of β-alanine is catalyzed by l-aspartate-α-decarboxylase (PanD), a pyruvoyl enzyme that is initially synthesized as a zymogen (pro-PanD). Active PanD is generated by self-cleavage of pro-PanD at Gly24-Ser25 creating the active-site pyruvoyl moiety. In Salmonella enterica, this cleavage requires PanM, an acetyl-CoA sensor related to the Gcn5-like N-acetyltransferases. PanM does not acetylate pro-PanD, but the recent publication of the three-dimensional crystal structure of the PanM homologue PanZ in complex with the PanD zymogen of Escherichia coli provides validation to our predictions and provides a framework in which to further examine the cleavage mechanism. In contrast, PanD from bacteria lacking PanM efficiently cleaved in the absence of PanM in vivo. Results: Using phylogenetic analyses combined with in vivo phenotypic investigations, we showed that two classes of bacterial l-aspartate-α-decarboxylases exist. This classification is based on their posttranslational activation by self-cleavage of its zymogen. Class I l-aspartate-α-decarboxylase zymogens require the acetyl-CoA sensor PanM to be cleaved into active PanD. This class is found exclusively in the Gammaproteobacteria. Class II l-aspartate-α-decarboxylase zymogens self cleave efficiently in the absence of PanM, and are found in a wide number of bacterial phyla. Several members of the Euryarchaeota and Crenarchaeota also contain Class II l-aspartate-α-decarboxylases. Phylogenetic and amino acid conservation analyses of PanM revealed a conserved region of PanM distinct from conserved regions found in related Gcn5-related acetyltransferase enzymes (Pfam00583). This conserved region represents a putative domain for interactions with l-aspartate-α-decarboxylase zymogens. This work may inform future biochemical and structural studies of pro-PanD-PanM interactions. Conclusions: Experimental results indicate that S. enterica and C. glutamicum l-aspartate-α-decarboxylases represent two different classes of homologues of these enzymes. Class I homologues require PanM for activation, while Class II self cleave in the absence of PanM. Computer modeling of conserved amino acids using structure coordinates of PanM and l-aspartate-α-decarboxylase available in the protein data bank (RCSB PDB) revealed a putative site of interactions, which may help generate models to help understand the molecular details of the self-cleavage mechanism of l-aspartate-α-decarboxylases.

Description: Background: Heterochromatin is essential for chromosome segregation, gene silencing and genome integrity. The fission yeast Schizosaccharomyces pombe contains heterochromatin at centromeres, subtelomeres, and mating type genes, as well as at small islands of meiotic genes dispersed across the genome. This heterochromatin is generated by partially redundant mechanisms, including the production of small interfering RNAs (siRNAs) that are incorporated into the RITS protein complex (RNAi-Induced Transcriptional Silencing). The assembly of heterochromatin islands requires the function of the RNA-binding protein Mmi1, which recruits RITS to its mRNA targets and to heterochromatin islands. In addition, Mmi1 directs its targets to an exosome-dependent RNA elimination pathway. Results: Ccr4-Not is a conserved multiprotein complex that regulates gene expression at multiple levels, including RNA degradation and translation. We show here that Ccr4-Not is recruited by Mmi1 to its RNA targets. Surprisingly, Ccr4 and Caf1 (the mRNA deadenylase catalytic subunits of the Ccr4-Not complex) are not necessary for the degradation or translation of Mmi1 RNA targets, but are essential for heterochromatin integrity at Mmi1-dependent islands and, independently of Mmi1, at subtelomeric regions. Both roles require the deadenylase activity of Ccr4 and the Mot2/Not4 protein, a ubiquitin ligase that is also part of the complex. Genetic evidence shows that Ccr4-mediated silencing is essential for normal cell growth, indicating that this novel regulation is physiologically relevant. Moreover, Ccr4 interacts with components of the RITS complex in a Mmi1-independent manner. Conclusions: Taken together, our results demonstrate that the Ccr4-Not complex is required for heterochromatin integrity in both Mmi1-dependent and Mmi1-independent pathways.

Description: Background: In higher eukaryotes, the genome is partitioned into large "Topologically Associating Domains" (TADs) in which the chromatin displays favoured long-range contacts. While a crumpled/fractal globule organization has received experimental supports at higher-order levels, the organization principles that govern chromatin dynamics within these TADs remain unclear. Using simple polymer models, we previously showed that, in mouse liver cells, gene-rich domains tend to adopt a statistical helix shape when no significant locus-specific interaction takes place. Results: Here, we use data from diverse 3C-derived methods to explore chromatin dynamics within mouse and Drosophila TADs. In mouse Embryonic Stem Cells (mESC), that possess large TADs (median size of 840 kb), we show that the statistical helix model, but not globule models, is relevant not only in gene-rich TADs, but also in gene-poor and gene-desert TADs. Interestingly, this statistical helix organization is considerably relaxed in mESC compared to liver cells, indicating that the impact of the constraints responsible for this organization is weaker in pluripotent cells. Finally, depletion of histone H1 in mESC alters local chromatin flexibility but not the statistical helix organization. In Drosophila, which possesses TADs of smaller sizes (median size of 70 kb), we show that, while chromatin compaction and flexibility are finely tuned according to the epigenetic landscape, chromatin dynamics within TADs is generally compatible with an unconstrained polymer configuration. Conclusions: Models issued from polymer physics can accurately describe the organization principles governing chromatin dynamics in both mouse and Drosophila TADs. However, constraints applied on this dynamics within mammalian TADs have a peculiar impact resulting in a statistical helix organization.

Description: Background: Bone resorption takes place within the basic multicellular units (BMU), and the surface to be resorbed is isolated from adjacent bone surfaces by a sealing zone between osteoclast membrane and bone matrix, which defines the limits of the resorption lacuna. Considering that the extracellular fluid (ECF) in both BMU and the resorption lacuna can be isolated from its surroundings, I hypothesize that flow and ion composition of the bone ECF in these sites might contribute to the regulation of osteoclast H + secretion. To investigate this hypothesis, I evaluated the H + secretion properties of individual osteoclasts and osteoclast-like cells (OCL-cells) and investigated whether changes in flow or chloride content of the extracellular solution modify the H + secretion properties in vitro. Results: The results show that 1) osteoclasts are unable to secrete H + and regulate intracellular pH (pH i ) under continuous flow conditions and exhibit progressive intracellular acidification; 2) the cessation of flow coincides with the onset of H + secretion and subsequent progressive intracellular alkalinization of osteoclasts and OCL-cells; 3) osteoclasts exhibit spontaneous rhythmic oscillations of pH i in non-flowing ECF, 4) pH i oscillations are not abolished by concanamycin, NPPB, or removal of extracellular Na + or Cl − ; 5) extracellular Cl − removal modifies the pattern of oscillations, by diminishing H + secretion; 6) pH i oscillations are abolished by continuous flowing of ECF over osteoclasts and OCL-cells. Conclusions: The data suggest, for the first time, that ECF flow and Cl − content have direct effects on osteoclast H + secretion and could be part of a mechanism determining the onset of osteoclast H + secretion required for bone resorption.

Description: Background: Field pea (Pisum sativum L.) is a cool-season grain legume that is cultivated world-wide for both human consumption and stock-feed purposes. Enhancement of genetic and genomic resources for field pea will permit improved understanding of the control of traits relevant to crop productivity and quality. Advances in second-generation sequencing and associated bioinformatics analysis now provide unprecedented opportunities for the development of such resources. The objective of this study was to perform transcriptome sequencing and characterisation from two genotypes of field pea that differ in terms of seed and plant morphological characteristics. Results: Transcriptome sequencing was performed with RNA templates from multiple tissues of the field pea genotypes Kaspa and Parafield. Tissue samples were collected at various growth stages, and a total of 23 cDNA libraries were sequenced using Illumina high-throughput sequencing platforms. A total of 407 and 352 million paired-end reads from the Kaspa and Parafield transcriptomes, respectively were assembled into 129,282 and 149,272 contigs, which were filtered on the basis of known gene annotations, presence of open reading frames (ORFs), reciprocal matches and degree of coverage. Totals of 126,335 contigs from Kaspa and 145,730 from Parafield were subsequently selected as the reference set. Reciprocal sequence analysis revealed that c. 87 % of contigs were expressed in both cultivars, while a small proportion were unique to each genotype. Reads from different libraries were aligned to the genotype-specific assemblies in order to identify and characterise expression of contigs on a tissue-specific basis, of which 87 % were expressed in more than one tissue, while others showed distinct expression patterns in specific tissues, providing unique transcriptome signatures. Conclusion: This study provided a comprehensive assembled and annotated transcriptome set for field pea that can be used for development of genetic markers, in order to assess genetic diversity, construct linkage maps, perform trait-dissection and implement whole-genome selection strategies in varietal improvement programs, as well to identify target genes for genetic modification approaches on the basis of annotation and expression analysis. In addition, the reference field pea transcriptome will prove highly valuable for comparative genomics studies and construction of a finalised genome sequence.

Description: Background: Continental-scale aboveground biomass maps are increasingly available, but their estimates vary widely, particularly at high resolution. A comprehensive understanding of map discrepancies is required to improve their effectiveness in carbon accounting and local decision-making. To this end, we compare four continental-scale maps with a recent high-resolution lidar-derived biomass map over Maryland, USA. We conduct detailed comparisons at pixel-, county-, and state-level. Results: Spatial patterns of biomass are broadly consistent in all maps, but there are large differences at fine scales (RMSD 48.5–92.7 Mg ha −1 ). Discrepancies reduce with aggregation and the agreement among products improves at the county level. However, continental scale maps exhibit residual negative biases in mean (33.0–54.6 Mg ha −1 ) and total biomass (3.5–5.8 Tg) when compared to the high-resolution lidar biomass map. Three of the four continental scale maps reach near-perfect agreement at ~4 km and onward but do not converge with the high-resolution biomass map even at county scale. At the State level, these maps underestimate biomass by 30–80 Tg in forested and 40–50 Tg in non-forested areas. Conclusions: Local discrepancies in continental scale biomass maps are caused by factors including data inputs, modeling approaches, forest/non-forest definitions and time lags. There is a net underestimation over high biomass forests and non-forested areas that could impact carbon accounting at all levels. Local, high-resolution lidar-derived biomass maps provide a valuable bottom-up reference to improve the analysis and interpretation of large-scale maps produced in carbon monitoring systems.

Description: Background: Phalaenopsis is one of the important commercial orchids in the world. Members of the P. amabilis species complex represent invaluable germplasm for the breeding program. However, the phylogeny of the P. amabilis species complex is still uncertain. The Phalaenopsis amabilis species complex (Orchidaceae) consists of subspecies amabilis, moluccana, and rosenstromii of P. amabilis, as well as P. aphrodite ssp. aphrodite, P. ap. ssp. formosana, and P. sanderiana. The aims of this study were to reconstruct the phylogeny and biogeographcial patterns of the species complex using Neighbor Joining (NJ), Maxinum Parsimony (MP), Bayesian Evolutionary Analysis Sampling Trees (BEAST) and Reconstruct Ancestral State in Phylogenies (RASP) analyses based on sequences of internal transcribed spacers 1 and 2 from the nuclear ribosomal DNA and the trnH-psbA spacer from the plastid DNA. Results: A pattern of vicariance, dispersal, and vicariance + dispersal among disjunctly distributed taxa was uncovered based on RASP analysis. Although two subspecies of P. aphrodite could not be differentiated from each other in dispersal state, they were distinct from P. amabilis and P. sanderiana. Within P. amabilis, three subspecies were separated phylogenetically, in agreement with the vicariance or vicariance + dispersal scenario, with geographic subdivision along Huxley’s, Wallace’s and Lydekker’s Lines. Molecular dating revealed such subdivisions among taxa of P. amabilis complex dating back to the late Pleistocene. Population-dynamic analyses using a Bayesian skyline plot suggested that the species complex experienced an in situ range expansion and population concentration during the late Last Glacial Maximum (LGM). Conclusions: Taxa of the P. amabilis complex with disjunct distributions were differentiated due to vicariance or vicariance + dispersal, with events likely occurring in the late Pleistocene. Demographic growth associated with the climatic oscillations in the Würm glacial period followed the species splits. Nevertheless, a subsequent population slowdown occurred in the late LGM due to extinction of regional populations. The reduction of suitable habitats resulted in geographic fragmenttation of the remaining taxa.

Description: Background: Despite the tremendous drop in the cost of nucleotide sequencing in recent years, many research projects still utilize sequencing of pools containing multiple samples for the detection of sequence variants as a cost saving measure. Various software tools exist to analyze these pooled sequence data, yet little has been reported on the relative accuracy and ease of use of these different programs. Results: In this manuscript we evaluate five different variant detection programs—The Genome Analysis Toolkit (GATK), CRISP, LoFreq, VarScan, and SNVer—with regard to their ability to detect variants in synthetically pooled Illumina sequencing data, by creating simulated pooled binary alignment/map (BAM) files using single-sample sequencing data from varying numbers of previously characterized samples at varying depths of coverage per sample. We report the overall runtimes and memory usage of each program, as well as each program’s sensitivity and specificity to detect known true variants. Conclusions: GATK, CRISP, and LoFreq all gave balanced accuracy of 80 % or greater for datasets with varying per-sample depth of coverage and numbers of samples per pool. VarScan and SNVer generally had balanced accuracy lower than 80 %. CRISP and LoFreq required up to four times less computational time and up to ten times less physical memory than GATK did, and without filtering, gave results with the highest sensitivity. VarScan and SNVer had generally lower false positive rates, but also significantly lower sensitivity than the other three programs.

Description: For the third time, BMC Ecology is delighted to announce the winners of our Image Competition. This year featured entries from all over the world and showcased not only the creativity and talent of the participants, but also the exquisite beauty and diversity of our planet. We are pleased to present the winning selections of the editorial board of the journal and guest judge Dr. Ana Luz Porzecanski, as well as some highly commended images that are sure to impress.

Description: Background: Acid soil is a serious limitation to crop production all over the world. Toxic aluminium (Al) cations in acid soil inhibit root growth and reduce yield. Although a gene tolerant to acid soil has been identified, it has not been used in malting barley breeding, which is partly due to the acid soil tolerance gene being linked to unfavorable malting quality traits. Results: A Brazilian malting barley variety Br2 was identified as tolerant to acid soil. A doubled haploid (DH) population was developed from a cross between Br2 and the Australian acid-sensitive cultivar Hamelin. The DH population was tested for acid soil tolerance in native acid soil and a hydroponic system with pH 4.2, pH 4.2 + Al or pH 6.5, and genotyped using SSR, DArT and gene-specific markers. A single QTL was detected for all parameters related to acid soil tolerance. The QTL was mapped to the known HvMATE location on chromosome 4H. Sequence alignment of the HvMATE gene identified 13 INDELs and 87 SNPs, where one SNP coded for a single amino acid difference between the two varieties. A gene-specific marker was developed to detect the single nucleotide polymorphism between Hamelin and Br2. This marker co-segregated with aluminium tolerance and accounted for 79 % of phenotypic variation for acid soil tolerance. Conclusion: The present study identified a novel source of acid soil/Al tolerance from a Brazilian malting barley cultivar Br2. This variety tolerated Al toxicity but was sensitive to low pH which is similar to most other Al-tolerant varieties. A gene-specific marker Cit7 was developed based on the HvMATE gene sequence. Cit7 will improve the efficiency of molecular-assisted selection of new barley varieties with tolerance to acid soil. Multiple alleles exist for the acid soil tolerance gene on chromosome 4H, so a malting barley variety that tolerates acid soil could be developed by selecting suitable tolerant alleles. Tolerance to low pH may play an important role for barley to adapt to acid soils.

Description: Background: Formalin-fixed, paraffin-embedded (FFPE) tissues represent the most abundant resource of archived human specimens in pathology. Such tissue specimens are emerging as a highly valuable resource for translational proteomic studies. In quantitative proteomic analysis, reductive di-methylation of primary amines using stable isotopic formaldehyde variants is increasingly used due to its robustness and cost-effectiveness. Results: In the present study we show for the first time that isotopic amine dimethylation can be used in a straightforward manner for the quantitative proteomic analysis of FFPE specimens without interference from formalin employed in the FFPE process. Isotopic amine dimethylation of FFPE specimens showed equal labeling efficiency as for cryopreserved specimens. For both FFPE and cryopreserved specimens, differential labeling of identical samples yielded highly similar ratio distributions within the expected range for dimethyl labeling. In an initial application, we profiled proteome changes in clear cell renal cell carcinoma (ccRCC) FFPE tissue specimens compared to adjacent non–malignant renal tissue. Our findings highlight increased levels of glyocolytic enzymes, annexins as well as ribosomal and proteasomal proteins. Conclusion: Our study establishes isotopic amine dimethylation as a versatile tool for quantitative proteomic analysis of FFPE specimens and underlines proteome alterations in ccRCC.

Description: Background: Cultivated peanut (Arachis hypogaea. L) is one of the most important oil crops in the world. After flowering, the peanut plant forms aboveground pegs (gynophores) that penetrate the soil, giving rise to underground pods. This means of reproduction, referred to as geocarpy, distinguishes peanuts from most other plants. The molecular mechanism underlying geocarpic pod development in peanut is poorly understood. Results: To gain insight into the mechanism of geocarpy, we extracted proteins from aerial gynophores, subterranean unswollen gynophores, and gynophores that had just started to swell into pods. We analyzed the protein profiles in each of these samples by combining 1 DE with nanoLC-MS/MS approaches. In total, 2766, 2518, and 2280 proteins were identified from the three samples, respectively. An integrated analysis of proteome and transcriptome data revealed specifically or differentially expressed genes in the different developmental stages at both the mRNA and protein levels. A total of 69 proteins involved in the gravity response, light and mechanical stimulus, hormone biosynthesis, and transport were identified as being involved in geocarpy. Furthermore, we identified 91 genes that were specifically or abundantly expressed in aerial gynophores, including pectin methylesterase and expansin, which were presumed to promote the elongation of aerial gynophores. In addition, we identified 35 proteins involved in metabolism, defense, hormone biosynthesis and signal transduction, nitrogen fixation, and transport that accumulated in subterranean unswellen gynophores. Furthermore, 26 specific or highly abundant proteins related to fatty acid metabolism, starch synthesis, and lignin synthesis were identified in the early swelling pods. Conclusions: We identified thousands of proteins in the aerial gynophores, subterranean gynophores, and early swelling pods of peanut. This study provides the basis for examining the molecular mechanisms underlying peanut geocarpy pod development.

Description: Background: Mounting evidence shows that smartphone usage heavily disrupts our work life and social activities. Moreover, it is possible that overuse could resemble addictive tendencies. A key contributing factor to smartphone overuse seems to be usage of the messaging application WhatsApp. Although WhatsApp is one of the most commonly used communication applications on smartphones, research in this area is scarce. Given the huge societal debate on the impact of smartphone usage on our daily lives, the present study undertook a large-scale investigation in order to provide numbers on smartphone usage generally—and use of WhatsApp in particular, with the aim of providing a basis for a scientific debate. Methods: In a large sample of N = 2,418 users, we recorded WhatsApp behaviour over a 4 week period. Results: Our data show that use of WhatsApp accounted for 19.83% (= 32.11 min) of all smartphone behaviour (compare: Facebook only 9.38% = 15.19 min). The mean of general daily smartphone usage was 161.95 min. Females used WhatsApp for significantly longer periods of time than males and younger age was associated with longer duration of WhatsApp use. While the personality trait Extraversion was positively associated with daily WhatsApp use, Conscientiousness showed an inverse correlation with the length of daily WhatsApp use. Conclusions: The numbers on smartphone usage in the present study show that the smartphone dominates our daily life. In particular WhatsApp is a driving force, here. Given the length of daily smartphone and WhatsApp usage, more studies need to be conducted to better understand smartphone usage.

Description: Background: Interaction and genetic control for traits influencing the adaptation of the rice crop to varying environments was studied in a mapping population derived from parents (Moroberekan and Swarna) contrasting for drought tolerance, yield potential, lodging resistance, and adaptation to dry direct seeding. A BC 2 F 3 -derived mapping population for traits related to these four trait groups was phenotyped to understand the interactions among traits and to map and align QTLs using composite interval mapping (CIM). The study also aimed to identify QTLs for the four trait groups as composite traits using multivariate least square interval mapping (MLSIM) to further understand the genetic control of these traits. Results: Significant correlations between drought- and yield-related traits at seedling and reproductive stages respectively with traits for adaptation to dry direct-seeded conditions were observed. CIM and MLSIM methods were applied to identify QTLs for univariate and composite traits. QTL clusters showing alignment of QTLs for several traits within and across trait groups were detected at chromosomes 3, 4, and 7 through CIM. The largest number of QTLs related to traits belonging to all four trait groups were identified on chromosome 3 close to the qDTY 3.2 locus. These included QTLs for traits such as bleeding rate, shoot biomass, stem strength, and spikelet fertility. Multivariate QTLs were identified at loci supported by univariate QTLs such as on chromosomes 3 and 4 as well as at distinctly different loci on chromosome 8 which were undetected through CIM. Conclusion: Rice requires better adaptation across a wide range of environments and cultivation practices to adjust to climate change. Understanding the genetics and trade-offs related to each of these environments and cultivation practices thus becomes highly important to develop varieties with stability of yield across them. This study provides a wider picture of the genetics and physiology of adaptation of rice to wide range of environments. With a complete understanding of the processes and relationships between traits and trait groups, marker-assisted breeding can be used more efficiently to develop plant types that can combine all or most of the beneficial traits and show high stability across environments, ecosystems, and cultivation practices.

Description: Background: Lobosphaera incisa, formerly known as Myrmecia incisa and then Parietochloris incisa, is an oleaginous unicellular green alga belonging to the class Trebouxiophyceae (Chlorophyta). It is the richest known plant source of arachidonic acid, an ω-6 poly-unsaturated fatty acid valued by the pharmaceutical and baby-food industries. It is therefore an organism of high biotechnological interest, and we recently reported the sequence of its chloroplast genome. Results: We now report the complete sequence of the mitochondrial genome of L. incisa from high-throughput Illumina short-read sequencing. The circular chromosome of 69,997 bp is predicted to encode a total of 64 genes, some harboring specific self-splicing group I and group II introns. Overall, the gene content is highly similar to that of the mitochondrial genomes of other Trebouxiophyceae, with 34 protein-coding, 3 rRNA, and 27 tRNA genes. Genes are distributed in two clusters located on different DNA strands, a bipartite arrangement that suggests expression from two divergent promoters yielding polycistronic primary transcripts. The L. incisa mitochondrial genome contains families of intergenic dispersed DNA repeat sequences that are not shared with other known mitochondrial genomes of Trebouxiophyceae. The most peculiar feature of the genome is a repetitive palindromic repeat, the LIMP (L. Incisa Mitochondrial Palindrome), found 19 times in the genome. It is formed by repetitions of an AACCA pentanucleotide, followed by an invariant 7-nt loop and a complementary repeat of the TGGTT motif. Analysis of the genome sequencing reads indicates that the LIMP can be a substrate for large-scale genomic rearrangements. We speculate that LIMPs can act as origins of replication. Deep sequencing of the L. incisa transcriptome also suggests that the LIMPs with long stems are sites of transcript processing. The genome also contains five copies of a related palindromic repeat, the HyLIMP, with a 10-nt motif related to that of the LIMP. Conclusions: The mitochondrial genome of L. incisa encodes a unique type of repetitive palindromic repeat sequence, the LIMP, which can mediate genome rearrangements and play a role in mitochondrial gene expression. Experimental studies are needed to confirm and further characterize the functional role(s) of the LIMP.

Description: Background: The cereal rust fungi are destructive pathogens that affect grain production worldwide. Although the genomic and transcript sequences for three Puccinia species that attack wheat have been released, the functions of large repertories of genes from Puccinia still need to be addressed to understand the infection process of these obligate parasites. Host-induced gene silencing (HIGS) has emerged a useful tool to examine the importance of rust fungus genes while growing within host plants. In this study, HIGS was used to test genes from Puccinia with transcripts enriched in haustoria for their ability to interfere with full development of the rust fungi. Results: Approximately 1200 haustoria enriched genes from Puccinia graminis f. sp. tritici (Pgt) were identified by comparative RNA sequencing. Virus-induced gene silencing (VIGS) constructs with fragments of 86 Puccinia genes, were tested for their ability to interfere with full development of these rust fungi. Most of the genes tested had no noticeable effects, but 10 reduced Pgt development after co-inoculation with the gene VIGS constructs and Pgt. These included a predicted glycolytic enzyme, two other proteins that are probably secreted and involved in carbohydrate or sugar metabolism, a protein involved in thiazol biosynthesis, a protein involved in auxin biosynthesis, an amino acid permease, two hypothetical proteins with no conserved domains, a predicted small secreted protein and another protein predicted to be secreted with similarity to bacterial proteins involved in membrane transport. Transient silencing of four of these genes reduced development of P. striiformis (Pst), and three of also caused reduction of P. triticina (Pt) development. Conclusions: Partial suppression of transcripts involved in a large variety of biological processes in haustoria cells of Puccinia rusts can disrupt their development. Silencing of three genes resulted in suppression of all three rust diseases indicating that it may be possible to engineer durable resistance to multiple rust pathogens with a single gene in transgenic wheat plants for sustainable control of cereal rusts.

Description: Background: We present a genome-wide messenger RNA (mRNA) sequencing technique that converts small amounts of RNA from many samples into molecular phenotypes. It encompasses all steps from sample preparation to sequence analysis and is applicable to baseline profiling or perturbation measurements. Results: Multiplex sequencing of transcript 3′ ends identifies differential transcript abundance independent of gene annotation. We show that increasing biological replicate number while maintaining the total amount of sequencing identifies more differentially abundant transcripts. Conclusions: This method can be implemented on polyadenylated RNA from any organism with an annotated reference genome and in any laboratory with access to Illumina sequencing.

Description: Background: Verticillium wilt (VW) and Fusarium wilt (FW), caused by the soil-borne fungi Verticillium dahliae and Fusarium oxysporum f. sp. vasinfectum, respectively, are two most destructive diseases in cotton production worldwide. Root-knot nematodes (Meloidogyne incognita, RKN) and reniform nematodes (Rotylenchulus reniformis, RN) cause the highest yield loss in the U.S. Planting disease resistant cultivars is the most cost effective control method. Numerous studies have reported mapping of quantitative trait loci (QTLs) for disease resistance in cotton; however, very few reliable QTLs were identified for use in genomic research and breeding. Results: This study first performed a 4-year replicated test of a backcross inbred line (BIL) population for VW resistance, and 10 resistance QTLs were mapped based on a 2895 cM linkage map with 392 SSR markers. The 10 VW QTLs were then placed to a consensus linkage map with other 182 VW QTLs, 75 RKN QTLs, 27 FW QTLs, and 7 RN QTLs reported from 32 publications. A meta-analysis of QTLs identified 28 QTL clusters including 13, 8 and 3 QTL hotspots for resistance to VW, RKN and FW, respectively. The number of QTLs and QTL clusters on chromosomes especially in the A-subgenome was significantly correlated with the number of nucleotide-binding site (NBS) genes, and the distribution of QTLs between homeologous A- and D- subgenome chromosomes was also significantly correlated. Conclusions: Ten VW resistance QTL identified in a 4-year replicated study have added useful information to the understanding of the genetic basis of VW resistance in cotton. Twenty-eight disease resistance QTL clusters and 24 hotspots identified from a total of 306 QTLs and linked SSR markers provide important information for marker-assisted selection and high resolution mapping of resistance QTLs and genes. The non-overlapping of most resistance QTL hotspots for different diseases indicates that their resistances are controlled by different genes.

Description: Background: Genomic imprinting is an epigenetic mechanism that can lead to differential gene expression depending on the parent-of-origin of a received allele. While most studies on imprinting address its underlying molecular mechanisms or attempt at discovering genomic regions that might be subject to imprinting, few have focused on the amount of phenotypic variation contributed by such epigenetic process. In this report, we give a brief review of a one-locus imprinting model in a quantitative genetics framework, and provide a decomposition of the genetic variance according to this model. Analytical deductions from the proposed imprinting model indicated a non-negligible contribution of imprinting to genetic variation of complex traits. Also, we performed a whole-genome scan analysis on mouse body mass index (BMI) aiming at revealing potential consequences when existing imprinting effects are ignored in genetic analysis. Results: 10,021 SNP markers were used to perform a whole-genome single marker regression on mouse BMI using an additive and an imprinting model. Markers significant for imprinting indicated that BMI is subject to imprinting. Marked variance changed from 1.218 ×10 −4 to 1.842 ×10 −4 when imprinting was considered in the analysis, implying that one third of marked variance would be lost if existing imprinting effects were not accounted for. When both marker and pedigree information were used, estimated heritability increased from 0.176 to 0.195 when imprinting was considered. Conclusions: When a complex trait is subject to imprinting, using an additive model that ignores this phenomenon may result in an underestimate of additive variability, potentially leading to wrong inferences about the underlying genetic architecture of that trait. This could be a possible factor explaining part of the missing heritability commonly observed in genome-wide association studies (GWAS).

Description: Background: Enset (Ensete ventricosum (Welw.) Cheesman; Musaceae) is a multipurpose drought-tolerant food security crop with high conservation and improvement concern in Ethiopia, where it supplements the human calorie requirements of around 20 million people. The crop also has an enormous potential in other regions of Sub-Saharan Africa, where it is known only as a wild plant. Despite its potential, genetic and genomic studies supporting breeding programs and conservation efforts are very limited. Molecular methods would substantially improve current conventional approaches. Here we report the development of the first set of SSR markers from enset, their cross-transferability to Musa spp., and their application in genetic diversity, relationship and structure assessments in wild and cultivated enset germplasm. Results: SSR markers specific to E. ventricosum were developed through pyrosequencing of an enriched genomic library. Primer pairs were designed for 217 microsatellites with a repeat size 〉 20 bp from 900 candidates. Primers were validated in parallel by in silico and in vitro PCR approaches. A total of 67 primer pairs successfully amplified specific loci and 59 showed polymorphism. A subset of 34 polymorphic SSR markers were used to study 70 both wild and cultivated enset accessions. A large number of alleles were detected along with a moderate to high level of genetic diversity. AMOVA revealed that intra-population allelic variations contributed more to genetic diversity than inter-population variations. UPGMA based phylogenetic analysis and Discriminant Analysis of Principal Components show that wild enset is clearly separated from cultivated enset and is more closely related to the out-group Musa spp. No cluster pattern associated with the geographical regions, where this crop is grown, was observed for enset landraces. Our results reaffirm the long tradition of extensive seed-sucker exchange between enset cultivating communities in Southern Ethiopia. Conclusion: The first set of genomic SSR markers were developed in enset. A large proportion of these markers were polymorphic and some were also transferable to related species of the genus Musa. This study demonstrated the usefulness of the markers in assessing genetic diversity and structure in enset germplasm, and provides potentially useful information for developing conservation and breeding strategies in enset.

Description: Background: Mollusca is an extremely diverse animal phylum that includes the aculiferans (worm-like aplacophorans and eight-shelled polyplacophorans) and their sister group, the conchiferans, comprising monoplacophorans, bivalves (clams, mussels), gastropods (snails, slugs), scaphopods (tusk shells) and cephalopods (squids, octopuses). Studies on mollusks have revealed an overall number of 11 Hox genes in seven out of eight molluscan “class”-level taxa, but expression data of key developmental regulators such as homeotic genes are only available for three gastropod and two cephalopod species. These show that Hox genes are involved in the formation of specific features including shell, foot, funnel or tentacles and not in antero-posterior body plan patterning as in most other bilaterian animals. The role of Hox genes in non-conchiferan (i.e., aculiferan) mollusks remains entirely unknown. Results: Here we present the first data on the expression of seven Hox genes in apolyplacophoran mollusk, Acanthochitona crinita. In A. crinita the Hox genes Acr-Hox1-5, Hox7 and Post2 are expressed in a co-linear pattern along the antero-posterior axis, but not in molluscan-specific features such as the shell or the foot. The expression pattern is restricted to the post-trochal region and the transcripts are present in ecto-, endo- and mesodermal cell layers. Contrary to the situation in gastropods and cephalopods, we did neither find Hox gene expression in distinct neural subsets of A. crinita, nor in its developing shell plates. Conclusions: Our analysis and comparison with other lophotrochozoans indicate that the basal role of Hox genes is in antero-posterior axis patterning in mollusks, similar to the vast majority of bilaterian animals, and that this role has been conserved in polyplacophorans, while co-option into patterning of evolutionary novelties emerged either at the base of Conchifera or independently in gastropods and cephalopods. These morphological innovations most likely contributed to the evolutionary success of its representatives, as exemplified by, e.g., the wide ecological range and species richness of gastropods.

Description: Background: The use of transcriptomic and genomic datasets for phylogenetic reconstruction has become increasingly common as researchers attempt to resolve recalcitrant nodes with increasing amounts of data. The large size and complexity of these datasets introduce significant phylogenetic noise and conflict into subsequent analyses. The sources of conflict may include hybridization, incomplete lineage sorting, or horizontal gene transfer, and may vary across the phylogeny. For phylogenetic analysis, this noise and conflict has been accommodated in one of several ways: by binning gene regions into subsets to isolate consistent phylogenetic signal; by using gene-tree methods for reconstruction, where conflict is presumed to be explained by incomplete lineage sorting (ILS); or through concatenation, where noise is presumed to be the dominant source of conflict. The results provided herein emphasize that analysis of individual homologous gene regions can greatly improve our understanding of the underlying conflict within these datasets. Results: Here we examined two published transcriptomic datasets, the angiosperm group Caryophyllales and the aculeate Hymenoptera, for the presence of conflict, concordance, and gene duplications in individual homologs across the phylogeny. We found significant conflict throughout the phylogeny in both datasets and in particular along the backbone. While some nodes in each phylogeny showed patterns of conflict similar to what might be expected with ILS alone, the backbone nodes also exhibited low levels of phylogenetic signal. In addition, certain nodes, especially in the Caryophyllales, had highly elevated levels of strongly supported conflict that cannot be explained by ILS alone. Conclusion: This study demonstrates that phylogenetic signal is highly variable in phylogenomic data sampled across related species and poses challenges when conducting species tree analyses on large genomic and transcriptomic datasets. Further insight into the conflict and processes underlying these complex datasets is necessary to improve and develop adequate models for sequence analysis and downstream applications. To aid this effort, we developed the open source software phyparts (https://bitbucket.org/blackrim/phyparts), which calculates unique, conflicting, and concordant bipartitions, maps gene duplications, and outputs summary statistics such as internode certainy (ICA) scores and node-specific counts of gene duplications.

Description: Background: Bangladesh is facing serious shortage of trained health professionals. In the pluralistic healthcare system of Bangladesh, formal health care providers constitute only 5 % of the total workforce; the rest are informal health care providers. Information Communication Technologies (ICTs) are increasingly seen as a powerful tool for linking the community with formal healthcare providers. Our study assesses an intervention that linked village doctors (a cadre of informal health care providers practising modern medicine) to formal doctors through call centres from the perspective of the village doctors who participated in the intervention. Methods: The study was conducted in Chakaria, a remote rural area in south-eastern Bangladesh during April–May 2013. Twelve village doctors were selected purposively from a pool of 55 village doctors who participated in the mobile health (mHealth) intervention. In depth interviews were conducted to collect data. The data were manually analysed using themes that emerged.ResultThe village doctors talked about both business benefits (access to formal doctors, getting support for decision making, and being entitled to call trained doctors) and personal benefits (both financial and non-financial). Some of the major barriers mentioned were technical problems related to accessing the call centre, charging consultation fees, and unfamiliarity with the call centre physicians. Conclusion: Village doctors saw many benefits to having a business relationship with the trained doctors that the mHealth intervention provided. mHealth through call centres has the potential to ensure consultation services to populations through existing informal healthcare providers in settings with a shortage of qualified healthcare providers.

Description: Background: Prostate cancer is the most common form of cancer and the second leading cause of cancer death in North America. Auto-detection of prostate cancer can play a major role in early detection of prostate cancer, which has a significant impact on patient survival rates. While multi-parametric magnetic resonance imaging (MP-MRI) has shown promise in diagnosis of prostate cancer, the existing auto-detection algorithms do not take advantage of abundance of data available in MP-MRI to improve detection accuracy. The goal of this research was to design a radiomics-based auto-detection method for prostate cancer via utilizing MP-MRI data. Methods: In this work, we present new MP-MRI texture feature models for radiomics-driven detection of prostate cancer. In addition to commonly used non-invasive imaging sequences in conventional MP-MRI, namely T2-weighted MRI (T2w) and diffusion-weighted imaging (DWI), our proposed MP-MRI texture feature models incorporate computed high-b DWI (CHB-DWI) and a new diffusion imaging modality called correlated diffusion imaging (CDI). Moreover, the proposed texture feature models incorporate features from individual b-value images. A comprehensive set of texture features was calculated for both the conventional MP-MRI and new MP-MRI texture feature models. We performed feature selection analysis for each individual modality and then combined best features from each modality to construct the optimized texture feature models. Results: The performance of the proposed MP-MRI texture feature models was evaluated via leave-one-patient-out cross-validation using a support vector machine (SVM) classifier trained on 40,975 cancerous and healthy tissue samples obtained from real clinical MP-MRI datasets. The proposed MP-MRI texture feature models outperformed the conventional model (i.e., T2w+DWI) with regard to cancer detection accuracy. Conclusions: Comprehensive texture feature models were developed for improved radiomics-driven detection of prostate cancer using MP-MRI. Using a comprehensive set of texture features and a feature selection method, optimal texture feature models were constructed that improved the prostate cancer auto-detection significantly compared to conventional MP-MRI texture feature models.

Description: Background: Choosing the most appropriate family physician (FP) for the individual, plays a fundamental role in primary care. The aim of this study is to determine the selection criteria for the patients in choosing their family doctors and priority ranking of these criteria by using the multi-criteria decision-making method of the Analytic Hierarchy Process (AHP) model. Methods: The study was planned and conducted in two phases. In the first phase, factors affecting the patients’ decisions were revealed with a qualitative research. In the next phase, the priorities of FP selection criteria were determined by using AHP model. Criteria were compared in pairs. 96 patient were asked to fill the information forms which contains comparison scores in the Family Health Centres. Results: According to the analysis of focus group discussions FP selection criteria were congregated in to five groups: Individual Characteristics, Patient-Doctor relationship, Professional characteristics, the Setting, and Ethical Characteristics.For each of the 96 participants, comparison matrixes were formed based on the scores of their information forms. Of these, models of only 5 (5.2 %) of the participants were consistent, in other words, they have been able to score consistent ranking. The consistency ratios (CR) were found to be smaller than 0.10. Therefore the comparison matrix of this new model, which was formed based on the medians of scores only given by these 5 participants, was consistent (CR = 0.06 

Description: Background: The incidence of nonalcoholic fatty liver disease (NAFLD) has been increasing worldwide in parallel with the obesity epidemic. This study aims to investigate the effects of the total flavonoids in Stellera chamaejasme L. (TFSC) on the experimental NAFLD in high fat diet fed (HFD) rats. Methods: NAFLD model was induced in male Wistar rats by high-fat diet, and the rats in NAFLD group were randomized into NAFLD group (n = 20) and TFSC-treated group (n = 60). Both groups were given high-fat diet, and the normal group (n = 20) was given normal diet. In addition, the TFSC treated group was administered TFSC orally once a day at a low dose of 100 mg/kg (n = 20), medium dose of 200 mg/kg (n = 20), and high dose of 400 mg/kg (n = 20) for 6 weeks. Subsequently, the rats were sacrificed and body weight changes, lipid profiles in plasma and liver pathology were examined. The relative levels of fatty acid synthesis and β-oxidation gene expression in hepatic tissues were measured by quantitative real-time polymerase chain reaction (RT-PCR). Results: After the HFD administration for 4 weeks, the body weight,serum TC and TG levels in the rat of model group were significantly higher than in normal group (P 

Description: No description available

Description: Shear stress in arteries, which is a measure of the force exerted by blood flow on the arterial wall, is associated with the location of lipid plaques that cause heart disease. In this study, a mathematical model of shear stress was combined with cross-sectional x-ray images of an artery taken using Computed Tomography (CT) scanning, allowing the authors to explore patterns of shearing stress and shed light on the role of arterial architecture in heart disease.

Description: Background: Female fertility is an important trait in cattle breeding programs. In the Nordic countries selection is based on a fertility index (FTI). The fertility index is a weighted combination of four female fertility traits estimated breeding values for number of inseminations per conception (AIS), 56-day non-return rate (NRR), number of days from first to last insemination (IFL), and number of days between calving and first insemination (ICF). The objective of this study was to identify associations between sequence variants and fertility traits in Jersey cattle based on 1,225 Jersey sires from Denmark with official breeding values for female fertility traits. The association analyses were carried out in two steps: first the cattle genome was scanned for quantitative trait loci using a sire model for FTI using imputed whole genome sequence variants; second the significant quantitative trait locus regions were re-analyzed using a linear mixed model (animal model) for both FTI and its component traits AIS, NRR, IFL and ICF. The underlying traits were analyzed separately for heifers (first parity cows) and cows (later parity cows) for AIS, NRR, and IFL. Results: In the first step 6 QTL were detected for FTI: one QTL on each of BTA7, BTA20, BTA23, BTA25, and two QTL on BTA9 (QTL9–1 and QTL9–2). In the second step, ICF showed association with the QTL regions on BTA7, QTL9–2 QTL2 on BTA9, and BTA25, AIS for cows on BTA20 and BTA23, AIS for heifers on QTL9–2 on BTA9, IFL for cows on BTA20, BTA23 and BTA25, IFL for heifers on BTA7 and QTL9-2 on BTA9, NRR for heifers on BTA7 and BTA23, and NRR for cows on BTA23. Conclusion: The genome wide association study presented here revealed 6 genomic regions associated with FTI. Screening these 6 QTL regions for the underlying female fertility traits revealed that different female fertility traits showed associations with different subsets of the individual FTI QTL peaks. The result of this study contributed to a better insight into the genetic control of FTI in the Danish Jersey.

Description: Background: Major secondary metabolites, including flavonoids, caffeine, and theanine, are important components of tea products and are closely related to the taste, flavor, and health benefits of tea. Secondary metabolite biosynthesis in Camellia sinensis is differentially regulated in different tissues during growth and development. Until now, little was known about the expression patterns of genes involved in secondary metabolic pathways or their regulatory mechanisms. This study aimed to generate expression profiles for C. sinensis tissues and to build a gene regulation model of the secondary metabolic pathways. Results: RNA sequencing was performed on 13 different tissue samples from various organs and developmental stages of tea plants, including buds and leaves of different ages, stems, flowers, seeds, and roots. A total of 43.7 Gbp of raw sequencing data were generated, from which 347,827 unigenes were assembled and annotated. There were 46,693, 8446, 3814, 10,206, and 4948 unigenes specifically expressed in the buds and leaves, stems, flowers, seeds, and roots, respectively. In total, 1719 unigenes were identified as being involved in the secondary metabolic pathways in C. sinensis, and the expression patterns of the genes involved in flavonoid, caffeine, and theanine biosynthesis were characterized, revealing the dynamic nature of their regulation during plant growth and development. The possible transcription factor regulation network for the biosynthesis of flavonoid, caffeine, and theanine was built, encompassing 339 transcription factors from 35 families, namely bHLH, MYB, and NAC, among others. Remarkably, not only did the data reveal the possible critical check points in the flavonoid, caffeine, and theanine biosynthesis pathways, but also implicated the key transcription factors and related mechanisms in the regulation of secondary metabolite biosynthesis. Conclusions: Our study generated gene expression profiles for different tissues at different developmental stages in tea plants. The gene network responsible for the regulation of the secondary metabolic pathways was analyzed. Our work elucidated the possible cross talk in gene regulation between the secondary metabolite biosynthetic pathways in C. sinensis. The results increase our understanding of how secondary metabolic pathways are regulated during plant development and growth cycles, and help pave the way for genetic selection and engineering for germplasm improvement.

Description: Background: Evidence is accumulating that perturbation of early life microbial colonization of the gut induces long-lasting adverse health effects in individuals. Understanding the mechanisms behind these effects will facilitate modulation of intestinal health. The objective of this study was to identify biological processes involved in these long lasting effects and the (molecular) factors that regulate them. We used an antibiotic and the same antibiotic in combination with stress on piglets as an early life perturbation. Then we used host gene expression data from the gut (jejunum) tissue and community-scale analysis of gut microbiota from the same location of the gut, at three different time-points to gauge the reaction to the perturbation. We analysed the data by a new combination of existing tools. First, we analysed the data in two dimensions, treatment and time, with quadratic regression analysis. Then we applied network-based data integration approaches to find correlations between host gene expression and the resident microbial species. Results: The use of a new combination of data analysis tools allowed us to identify significant long-lasting differences in jejunal gene expression patterns resulting from the early life perturbations. In addition, we were able to identify potential key gene regulators (hubs) for these long-lasting effects. Furthermore, data integration also showed that there are a handful of bacterial groups that were associated with temporal changes in gene expression. Conclusion: The applied systems-biology approach allowed us to take the first steps in unravelling biological processes involved in long lasting effects in the gut due to early life perturbations. The observed data are consistent with the hypothesis that these long lasting effects are due to differences in the programming of the gut immune system as induced by the temporary early life changes in the composition and/or diversity of microbiota in the gut.

Description: Background: Trifoliate orange (Poncirus trifoliata (L.) Raf.) is extremely cold hardy after a full acclimation; however the underlying molecular mechanisms underlying this economically valuable trait remain poorly understood. In this study, global transcriptome profiles of trifoliate orange under cold conditions (4 °C) over a time course were generated by high-throughput sequencing. Results: More than 68 million high-quality reads were produced and assembled into a non-redundant data of 77,292 unigenes with an average length of 1112 bp (N50 = 1778 bp). Of these, 23,846 had significant sequence similarity to known genes and these were assigned to 61 gene ontology (GO) categories and 25 clusters of orthologous groups (COG) involved in 128 KEGG pathways. Sequences derived from cold-treated and control plants were mapped to the assembled transcriptome, resulting in the identification of 5549 differentially expressed genes (DEGs). These comprised 600 (462 up-regulated, 138 down-regulated), 2346 (1631 up-regulated, 715 down-regulated), and 5177 (2702 up-regulated, 2475 down-regulated) genes from the cold-treated samples at 6, 24 and 72 h, respectively. The accuracy of the RNA-seq derived transcript expression data was validated by analyzing the expression patterns of 17 DEGs by qPCR. Plant hormone signal transduction, plant-pathogen interaction, and secondary metabolism were the most significantly enriched GO categories amongst in the DEGs. A total of 60 transcription factors were shown to be cold responsive. In addition, a number of genes involved in the catabolism and signaling of hormones, such as abscisic acid, ethylene and gibberellin, were affected by the cold stress. Meanwhile, levels of putrescine progressively increased under cold, which was consistent with up-regulation of an arginine decarboxylase gene. Conclusions: This dataset provides valuable information regarding the trifoliate orange transcriptome changes in response to cold stress and may help guide future identification and functional analysis of genes that are importnatn for enhancing cold hardiness.

Description: Background: About 1,000 endangered Patagonian huemul deer (Hippocamelus bisulcus) remain in Chile and 350–500 in Argentina. Most groups (〉100) are not recovering, and prevalence of osteopathology in Argentina was at least 57%. Here I describe relevant cases of osteopathology from a Chilean population which, however, recently also provided data on trace mineral status, supporting the initial hypothesis that nutrition may be a primary etiologic factor. Additionally, recent data on bone chemical composition of Argentine cases and soil analyses are discussed. Results: Fluoride levels in Argentine cases with osteopathology were low and fluorosis was discarded as an etiological factor. Selenium deficiency occurred in 73% of huemul from the Chilean population which exhibited several cases with osteopathology. The pathophysiognomy included extensive erosion; tooth loss; porosification; perforations of palate, maxillar and mandibular bone with frequent exposure of tooth roots; and fractured mandibula. Areas currently used by remaining huemul have mainly acidic volcanic soils, which reduces selenium bioavailability: mean soil selenium levels from areas typically used by extant huemul were very deficient (0.19 ppm), corroborating documented overt selenium deficiency in local livestock and plants. The area of extant huemul is known to result in primary iodine deficiency in livestock which is aggravated by selenium deficiency. Conclusions: Currently the most parsimonious explanation for frequent osteopathology and lack of numerical recovery are the combined effects of selenium and iodine deficiencies based on: osteopathology in a population of selenium deficient huemul; selenium deficient livestock, plants and soils; acidic soils; and regional primary iodine deficiency. The nexus between mineral nutrition and population dynamics of huemul may be due to constraints on their movements to fertile lowlands, including the elimination of historic migratory traditions, and concomitant elimination of source populations.

Description: Adult body size is controlled by the mechanisms that stop growth when a species-characteristic size has been reached. The mechanisms by which size is sensed and by which this information is transduced to the growth regulating system are beginning to be understood in a few species of insects. Two rather different strategies for control have been discovered; one favors large body size and the other favors rapid development.

Description: Background: Multiple hypothesis testing is a pervasive problem in genomic data analysis. The conventional Bonferroni method which controls the family-wise error rate is conservative and with low power. The current paradigm is to control the false discovery rate. Results: We characterize the variability of the false discovery rate indices (local false discovery rates, q-value and false discovery proportion) using the bootstrapped method. A colon cancer gene-expression data and a visual refractive errors genome-wide association study data are analyzed as demonstration. We found a high variability in false discovery rate controls for typical genomic studies. Conclusions: We advise researchers to present the bootstrapped standard errors alongside with the false discovery rate indices.

Description: Background: Chickens are susceptible to infection with a limited number of Influenza A viruses and are a potential source of a human influenza pandemic. In particular, H5 and H7 haemagglutinin subtypes can evolve from low to highly pathogenic strains in gallinaceous poultry. Ducks on the other hand are a natural reservoir for these viruses and are able to withstand most avian influenza strains. Results: Transcriptomic sequencing of lung and ileum tissue samples from birds infected with high (H5N1) and low (H5N2) pathogenic influenza viruses has allowed us to compare the early host response to these infections in both these species. Chickens (but not ducks) lack the intracellular receptor for viral ssRNA, RIG-I and the gene for an important RIG-I binding protein, RNF135. These differences in gene content partly explain the differences in host responses to low pathogenic and highly pathogenic avian influenza virus in chicken and ducks. We reveal very different patterns of expression of members of the interferon-induced transmembrane protein (IFITM) gene family in ducks and chickens. In ducks, IFITM1, 2 and 3 are strongly up regulated in response to highly pathogenic avian influenza, where little response is seen in chickens. Clustering of gene expression profiles suggests IFITM1 and 2 have an anti-viral response and IFITM3 may restrict avian influenza virus through cell membrane fusion. We also show, through molecular phylogenetic analyses, that avian IFITM1 and IFITM3 genes have been subject to both episodic and pervasive positive selection at specific codons. In particular, avian IFITM1 showed evidence of positive selection in the duck lineage at sites known to restrict influenza virus infection. Conclusions: Taken together these results support a model where the IFITM123 protein family and RIG-I all play a crucial role in the tolerance of ducks to highly pathogenic and low pathogenic strains of avian influenza viruses when compared to the chicken.

Description: Background: The co-chaperone Hop [ h eat shock protein (HSP) o rganizing p rotein] has been shown to act as an adaptor for protein folding and maturation, in concert with Hsp70 and Hsp90. The hop gene is of eukaryotic origin. Likewise, the chloroplast elongation factor G (cEF-G) catalyzes the translocation step in chloroplast protein synthesis. The chl-fus gene, which encodes the cEF-G protein, is of plastid origin. Both proteins, Hop and cEF-G, derived from domain duplications. It was demonstrated that the nuclear chl-fus gene locates in opposite orientation to a hop gene in Glycine max. We explored 53 available plant genomes from Chlorophyta to higher plants, to determine whether the chl-fus gene was transferred directly downstream of the primordial hop in the proto-eukaryote host cell. Since both genes came from exon/module duplication events, we wanted to explore the involvement of introns in the early origin and the ensuing evolutionary changes in gene structure. Results: We reconstructed the evolutionary history of the two convergent plant genes, on the basis of their gene structure, microsynteny and microcolinearity, from 53 plant nuclear genomes. Despite a high degree (72 %) of microcolinearity among vascular plants, our results demonstrate that their adjacency was a product of chromosomal rearrangements. Based on predicted exon − intron structures, we inferred the molecular events giving rise to the current form of genes. Therefore, we propose a simple model of exon/module shuffling by intronic recombinations in which phase-0 introns were essential for domain duplication, and a phase-1 intron for transit peptide recruiting. Finally, we demonstrate a natural susceptibility of the intergenic region to recombine or delete, seriously threatening the integrity of the chl-fus gene for the future. Conclusions: Our results are consistent with the interpretation that the chl-fus gene was transferred from the chloroplast to a chromosome different from that of hop, in the primitive photosynthetic eukaryote, and much later before the appearance of angiosperms, it was recombined downstream of hop. Exon/module shuffling mediated by symmetric intron phases (i.e., phase-0 introns) was essential for gene evolution. The intergenic region is prone to recombine, risking the integrity of both genes.

Description: Background: Soil salinity is one of the most significant abiotic stresses affecting plant shoots and roots growth. The adjustment of root architecture to spatio-temporal heterogeneity in salinity is particularly critical for plant growth and survival. Bermudagrass (Cynodon dactylon) is a widely used turf and forage perennial grass with a high degree of salinity tolerance. Salinity appears to stimulate the growth of roots and decrease their mortality in tolerant bermudagrass. To estimate a broad spectrum of genes related to root elongation affected by salt stress and the molecular mechanisms that control the positive response of root architecture to salinity, we analyzed the transcriptome of bermudagrass root tips in response to salinity. Results: RNA-sequencing was performed in root tips of two bermudagrass genotypes contrasting in salt tolerance. A total of 237,850,130 high quality clean reads were generated and 250,359 transcripts were assembled with an average length of 1115 bp. Totally, 103,324 unigenes obtained with 53,765 unigenes (52 %) successfully annotated in databases. Bioinformatics analysis indicated that major transcription factor (TF) families linked to stress responses and growth regulation (MYB, bHLH, WRKY) were differentially expressed in root tips of bermudagrass under salinity. In addition, genes related to cell wall loosening and stiffening (xyloglucan endotransglucosylase/hydrolases, peroxidases) were identified. Conclusions: RNA-seq analysis identified candidate genes encoding TFs involved in the regulation of lignin synthesis, reactive oxygen species (ROS) homeostasis controlled by peroxidases, and the regulation of phytohormone signaling that promote cell wall loosening and therefore root growth under salinity.

Description: Background: In recent years the emergence of multidrug resistant Klebsiella pneumoniae strains has been an increasingly common event. This opportunistic species is one of the five main bacterial pathogens that cause hospital infections worldwide and multidrug resistance has been associated with the presence of high molecular weight plasmids. Plasmids are generally acquired through horizontal transfer and therefore is possible that systems that prevent the entry of foreign genetic material are inactive or absent. One of these systems is CRISPR/Cas. However, little is known regarding the clustered regularly interspaced short palindromic repeats and their associated Cas proteins (CRISPR/Cas) system in K. pneumoniae. The adaptive immune system CRISPR/Cas has been shown to limit the entry of foreign genetic elements into bacterial organisms and in some bacteria it has been shown to be involved in regulation of virulence genes. Thus in this work we used bioinformatics tools to determine the presence or absence of CRISPR/Cas systems in available K. pneumoniae genomes. Results: The complete CRISPR/Cas system was identified in two out of the eight complete K. pneumoniae genomes sequences and in four out of the 44 available draft genomes sequences. The cas genes in these strains comprises eight cas genes similar to those found in Escherichia coli, suggesting they belong to the type I-E group, although their arrangement is slightly different. As for the CRISPR sequences, the average lengths of the direct repeats and spacers were 29 and 33 bp, respectively. BLAST searches demonstrated that 38 of the 116 spacer sequences (33%) are significantly similar to either plasmid, phage or genome sequences, while the remaining 78 sequences (67%) showed no significant similarity to other sequences. The region where the CRISPR/Cas systems were located is the same in all the Klebsiella genomes containing it, it has a syntenic architecture, and is located among genes encoding for proteins likely involved in metabolism and resistance to antibiotics. Conclusions: The CRISPR/Cas system is not widely distributed in K. pneumoniae genomes, those present most likely belong to type I-E with few differences from the arrangement of the cse3 gene and most of the spacers have not been are not described yet. Given that the CRISPR/Cas system is scarcely distributed among K. pneumoniae genomes it is not clear whether it is involved in either immunity against foreign genetic material or virulence. We consider that this study represents a first step to understand the role of CRISPR/Cas in K. pneumoniae.

Description: Background: Gorham’s syndrome is a rare illness of unknown etiology. It is characterized by a local proliferation of blood or lymphatic vessels that in bones leads to progressive resorption and destruction. The cause of the disease is not elucidated, and therapeutic options remain limited.Case presentationWe report herein the case of a young female Caucasian patient aged 18 years with diffuse Gorham syndrome. In tissue specimens angiogenesis and massive lymphangiogenesis as well as the expression of vascular endothelial growth factor-A (VEGF-A) and neuropilins was observed. Lymphangiogenesis is a prominent feature of the disease and a number of lymphatic markers were found to be expressed, however only VEGF-A, but not vascular endothelial growth factor-C (VEGF-C) was found to be elevated in the circulation. Circulating levels of soluble VEGF receptor-1 were also not elevated. Furthermore, the patient responded favorably and the disease was stabilized following treatment with the beta-blocking agent Propranolol alone which acts on VEGF-A alone, but not on soluble VEGF receptor-1 levels. Conclusion: This suggests that the disease is dependent on VEGF-A, but on neither VEGF-C, the major driver of lymphangiogenesis, nor FLT1. Furthermore, Propranolol acts on VEGF-A but not FLT1 expression.

Description: Background: Spirulina (Arthrospira) platensis is the only cyanobacterium that in addition to being studied at the molecular level and subjected to gene manipulation, can also be mass cultivated in outdoor ponds for commercial use as a food supplement. Thus, encountering environmental changes, including temperature stresses, is common during the mass production of Spirulina. The use of cyanobacteria as an experimental platform, especially for photosynthetic gene manipulation in plants and bacteria, is becoming increasingly important. Understanding the mechanisms and protein-protein interaction networks that underlie low- and high-temperature responses is relevant to Spirulina mass production. To accomplish this goal, high-throughput techniques such as OMICs analyses are used. Thus, large datasets must be collected, managed and subjected to information extraction. Therefore, databases including (i) proteomic analysis and protein-protein interaction (PPI) data and (ii) domain/motif visualization tools are required for potential use in temperature response models for plant chloroplasts and photosynthetic bacteria.DescriptionsA web-based repository was developed including an embedded database, SpirPro, and tools for network visualization. Proteome data were analyzed integrated with protein-protein interactions and/or metabolic pathways from KEGG. The repository provides various information, ranging from raw data (2D-gel images) to associated results, such as data from interaction and/or pathway analyses. This integration allows in silico analyses of protein-protein interactions affected at the metabolic level and, particularly, analyses of interactions between and within the affected metabolic pathways under temperature stresses for comparative proteomic analysis. The developed tool, which is coded in HTML with CSS/JavaScript and depicted in Scalable Vector Graphics (SVG), is designed for interactive analysis and exploration of the constructed network. SpirPro is publicly available on the web at http://spirpro.sbi.kmutt.ac.th. Conclusions: SpirPro is an analysis platform containing an integrated proteome and PPI database that provides the most comprehensive data on this cyanobacterium at the systematic level. As an integrated database, SpirPro can be applied in various analyses, such as temperature stress response networking analysis in cyanobacterial models and interacting domain-domain analysis between proteins of interest.

Description: Background: The knowledge of the spatial organisation of the chromatin fibre in cell nuclei helps researchers to understand the nuclear machinery that regulates dna activity. Recent experimental techniques of the type Chromosome Conformation Capture (3c, or similar) provide high-resolution, high-throughput data consisting in the number of times any possible pair of dna fragments is found to be in contact, in a certain population of cells. As these data carry information on the structure of the chromatin fibre, several attempts have been made to use them to obtain high-resolution 3d reconstructions of entire chromosomes, or even an entire genome. The techniques proposed treat the data in different ways, possibly exploiting physical-geometric chromatin models. One popular strategy is to transform contact data into Euclidean distances between pairs of fragments, and then solve a classical distance-to-geometry problem. Results: We developed and tested a reconstruction technique that does not require translating contacts into distances, thus avoiding a number of related drawbacks. Also, we introduce a geometrical chromatin chain model that allows us to include sound biochemical and biological constraints in the problem. This model can be scaled at different genomic resolutions, where the structures of the coarser models are influenced by the reconstructions at finer resolutions. The search in the solution space is then performed by a classical simulated annealing, where the model is evolved efficiently through quaternion operators. The presence of appropriate constraints permits the less reliable data to be overlooked, so the result is a set of plausible chromatin configurations compatible with both the data and the prior knowledge. Conclusions: To test our method, we obtained a number of 3d chromatin configurations from hi-c data available in the literature for the long arm of human chromosome 1, and validated their features against known properties of gene density and transcriptional activity. Our results are compatible with biological features not introduced a priori in the problem: structurally different regions in our reconstructions highly correlate with functionally different regions as known from literature and genomic repositories.

Description: Background: The filamentous cyanobacterium Nostoc sp. strain PCC 7120 can fix N 2 when combined nitrogen is not available. Furthermore, it has to cope with reactive oxygen species generated as byproducts of photosynthesis and respiration. We have previously demonstrated the synthesis of Ser/Thr kinase Pkn22 as an important survival response of Nostoc to oxidative damage. In this study we wished to investigate the possible involvement of this kinase in signalling peroxide stress and nitrogen deprivation. Results: Quantitative RT-PCR experiments revealed that the pkn22 gene is induced in response to peroxide stress and to combined nitrogen starvation. Electrophoretic motility assays indicated that the pkn22 promoter is recognized by the global transcriptional regulators FurA and NtcA. Transcriptomic analysis comparing a pkn22-insertion mutant and the wild type strain indicated that this kinase regulates genes involved in important cellular functions such as photosynthesis, carbon metabolism and iron acquisition. Since metabolic changes may lead to oxidative stress, we investigated whether this is the case with nitrogen starvation. Our results rather invalidate this hypothesis thereby suggesting that the function of Pkn22 under nitrogen starvation is independent of its role in response to peroxide stress. Conclusions: Our analyses have permitted a more complete functional description of Ser/Thr kinase in Nostoc. We have decrypted the transcriptional regulation of the pkn22 gene, and analysed the whole set of genes under the control of this kinase in response to the two environmental changes often encountered by cyanobacteria in their natural habitat: oxidative stress and nitrogen deprivation.

Description: Background: Despite a number of recent reports of insect resistance to transgenic crops expressing insecticidal toxins from Bacillus thuringiensis (Bt), little is known about the mechanism of resistance to these toxins. The purpose of this study is to identify genes associated with the mechanism of Cry1F toxin resistance in European corn borer (Ostrinia nubilalis Hübner). For this, we compared the global transcriptomic response of laboratory selected resistant and susceptible O. nubilalis strain to Cry1F toxin. We further identified constitutive transcriptional differences between the two strains. Results: An O. nubilalis midgut transcriptome of 36,125 transcripts was assembled de novo from 106 million Illumina HiSeq and Roche 454 reads and used as a reference for estimation of differential gene expression analysis. Evaluation of gene expression profiles of midgut tissues from the Cry1F susceptible and resistant strains after toxin exposure identified a suite of genes that responded to the toxin in the susceptible strain (n = 1,654), but almost 20-fold fewer in the resistant strain (n = 84). A total of 5,455 midgut transcripts showed significant constitutive expression differences between Cry1F susceptible and resistant strains. Transcripts coding for previously identified Cry toxin receptors, cadherin and alkaline phosphatase and proteases were also differentially expressed in the midgut of the susceptible and resistant strains. Conclusions: Our current study provides a valuable resource for further molecular characterization of Bt resistance and insect response to Cry1F toxin in O. nubilalis and other pest species.

Description: Background: In gene expression analysis, overlapping genes, splice variants, and fusion transcripts are potential sources of data analysis artefacts, depending on how the observed intensity is assigned to one, or more genes. We here exemplify this by an in-depth analysis of the INS-IGF2 fusion transcript, which has recently been reported to be among the highest expressed transcripts in human pancreatic beta cells and its protein indicated as a novel autoantigen in Type 1 Diabetes. Results: Through RNA sequencing and variant specific qPCR analyses we demonstrate that the true abundance of INS-IGF2 is 〉20,000 fold lower than INS in human beta cells, and we suggest an explanation to the nature of the artefacts which have previously led to overestimation of the gene expression level in selected studies. We reinvestigated the previous reported findings of detection of INS-IGF2 using antibodies both in Western blotting and immunohistochemistry. We found that the one available commercial antibody (BO1P) raised against recombinant INS-IGF2 show strong cross-reaction to native proinsulin, and we did not detect INS-IGF2 protein in the human beta cell line EndoC-βH1. Furthermore, using highly sensitive proteomics analysis we could not demonstrate INS-IGF2 protein in samples of human islets nor in EndoC-βH1. Conclusions: Sequence features, such as fusion transcripts spanning multiple genes can lead to unexpected results in gene expression analysis, and care must be taken in generating and interpreting the results. For the specific case of INS-IGF2 we conclude that the abundance of the fusion transcript/protein is exceedingly lower than previously reported, and that current immuno-reagents available for detecting INS-IGF2 protein have a strong cross-reaction to native human proinsulin. Finally, we were unable to detect INS-IGF2 protein by proteomics analysis.

Description: Probiotic treatments are thought to increase the levels of commensal bacterial species that populate the human gut, causing no harm to their host and playing an important role in maintaining gut health. This study is an investigation of the effect of a probiotic treatment on the level of a known commensal bacterium in the guts of healthy human subjects, which was significantly increased with probiotic treatment compared with a control. The authors concluded that the probiotic may thus help to promote gut health.

Description: Background: Equations based on single-frequency bioelectrical impedance analysis at 50 kHz for determination of total body water content (TBW) have been previously validated in healthy non-sedated beagle dogs. We investigated whether these equations are predictive of TBW in various canine breeds by comparing the results of these equations with TBW values evaluated directly by deuterium oxide (D 2 O) dilution. Methods: Total body water content of 13 healthy adult pet dogs of various breeds was determined directly using D 2 O dilution and indirectly using previous equations based on values obtained with a portable bioelectric impedance device. Paired Student’s t-tests were used to compare TBW obtained by single-frequency bioelectrical impedance analysis and D 2 O dilution. A p-value of