Stochastic simulation is a key tool in population genetics, since the models involved are often analytically intractable and simulation is usually the only way of obtaining ground-truth data to evaluate inferences. Because of this, a large number of specialized simulation programs have been developed, each filling a particular niche, but with largely overlapping functionality and a substantial duplication of effort. Here, we introduce msprime version 1.0, which efficiently implements ancestry and mutation simulations based on the succinct tree sequence data structure and the tskit library. We summarize msprime’s many features, and show that its performance is excellent, often many times faster and more memory efficient than specialized alternatives. These high-performance features have been thoroughly tested and validated, and built using a collaborative, open source development model, which reduces duplication of effort and promotes software quality via community engagement.

We reported HIVID (high-throughput Viral Integration Detection), a novel experimental and computational method to detect the location of Hepatitis B Virus (HBV) integration breakpoints in Hepatocellular Carcinoma (HCC) genome. In this method, the fragments with HBV sequence were enriched by a set of HBV probes and then processed to high-throughput sequencing. In order to evaluate the performance of HIVID, we compared the results of HIVID with that of whole genome sequencing method (WGS) in 28 HCC tumors. We detected a total of 246 HBV integration breakpoints in HCC genome, 113 out of which were within 400bp upstream or downstream of 125 breakpoints identified by WGS method, covering 89.3% (125/140) of total breakpoints. The integration was located in the gene TERT, MLL4, and CCNE1. In addition, we discovered 133 novel breakpoints missed by WGS method, with 66.7% (10/15) of validation rate. Our study shows HIVID is a cost-effective methodology with high specificity and sensitivity to identify viral integration in human genome.

No abstract

Stochastic simulation is a key tool in population genetics, since the models involved are often analytically intractable and simulation is usually the only way of obtaining ground-truth data to evaluate inferences. Because of this necessity, a large number of specialised simulation programs have been developed, each filling a particular niche, but with largely overlapping functionality and a substantial duplication of effort. Here, we introduce msprime version 1.0, which efficiently implements ancestry and mutation simulations based on the succinct tree sequence data structure and tskit library. We summarise msprime's many features, and show that its performance is excellent, often many times faster and more memory efficient than specialised alternatives. These high-performance features have been thoroughly tested and validated, and built using a collaborative, open source development model, which reduces duplication of effort and promotes software quality via community engagement.

Melanisation has been considered to be an important virulence factor of Fonsecaea monophora. However, the biosynthetic mechanisms of melanisation remain unknown. We therefore used next generation sequencing technology to investigate the transcriptome and digital gene expression data, which are valuable resources to better understand the molecular and biological mechanisms regulating melanisation in F. monophora. We performed de novo transcriptome assembly and digital gene expression (DGE) profiling analyses of parent (CBS 122845) and albino (CBS 125194) strains using the Illumina RNA-seq system. A total of 17 352 annotated unigenes were found by BLAST search of NR, Swiss-Prot, Gene Ontology, Clusters of Orthologous Groups and Kyoto Encyclopedia of Genes and Genomes (KEGG) (E-value <1e‒5). A total of 2 283 unigenes were judged to be the differentially expressed between the two genotypes. We identified most of the genes coding for key enzymes involved in melanin biosynthesis pathways, including polyketide synthase (pks), multicopper oxidase (mco), laccase, tyrosinase and homogentisate 1,2-dioxygenase (hmgA). DEG analysis showed extensive down-regulation of key genes in the DHN pathway, while up-regulation was noted in the DOPA pathway of the albino mutant. The transcript levels of partial genes were confirmed by real time RT-PCR, while the crucial role of key enzymes was confirmed by either inhibitor or substrate tests in vitro. Meanwhile, numbers of genes involved in light sensing, cell wall synthesis, morphology and environmental stress were identified in the transcriptome of F. monophora. In addition, 3 353 SSRs (Simple Sequence Repeats) markers were identified from 21 600 consensus sequences. Blocking of the DNH pathway is the most likely reason of melanin deficiency in the albino strain, while the production of pheomelanin and pyomelanin were probably regulated by unknown transcription factors on upstream of both pathways. Most of genes involved in environmental tolerance to oxidants, irradiation and extreme temperatures were also assembled and annotated in transcriptomes of F. monophora. In addition, thousands of identified cSSR (combined SSR) markers will favour further genetic linkage studies. In conclusion, these data will contribute to understanding the regulation of melanin biosynthesis and help to improve the studies of pathogenicity of F. monophora.

No abstract

Background The woodwasp Sirex noctilio Fabricius is a major quarantine pest worldwide that was first discovered in China in 2013 and mainly harms Pinus sylvestris var. mongolica Litv.. S. nitobei Matsumura is a native species in China and is closely related to S. noctilio. Recently, the two woodwasps species were found attacking the P. sylvestris var. mongolica Litv in succession. The olfactory system is the foundation of insect behavior. Olfactory genes were identified through antennal transcriptome analysis. The expression profiles odorant binding proteins (OBPs) were analyzed with RT-qPCR. Results From our transcriptome analysis, 16 OBPs, 7 chemosensory proteins (CSPs), 41 odorant receptors (ORs), 8 gustatory receptors (GRs), 13 ionotropic receptors (IRs), and one sensory neuron membrane protein (SNMP) were identified in S. noctilio, while 15 OBPs, 6 CSPs, 43 ORs, 10 GRs, 16 IRs, and 1 SNMP were identified in S. nitobei. Most of the olfactory genes identified in two species were homologous. However, some species-specific olfactory genes were identified from the antennal transcriptomes, including SnocOBP13, SnocCSP6, SnocOR26, SnocGR2, SnocIR7 in S. noctilio and SnitGR9, SnitGR11, SnitIR17 in S. nitobei. In total, 14 OBPs were expressed primarily in the antennae. SnocOBP9 and SnitOBP9, identified as PBP homologues, were sex-biased expression in two siricid, but with different pattern. SnocOBP11 and SnitOBP11 were highly expressed in antennae and clearly expressed in external genitalia. SnocOBP7 and SnitOBP7 were highly expressed in male genitalia. SnocOBP3 and SnocOBP10 were highly expressed in female genitalia and male heads, while SnitOBP3 and SnitOBP10 did not show obvious tissue bias. Conclusion We analyzed 86 and 91 olfactory genes from S. noctilio and S. nitobei, respectively. Most of the olfactory genes identified were homologous, but also some species-specific olfactory genes were identified, which indicated the similarities and differences of the molecular mechanisms between the two closely-related species. Different expression in the antennae, external genitals or heads, exhibiting an obvious sex bias, suggested their different role in recognizing sex pheromones or plant volatiles. Species-specific expression for several OBPs genes may suggest that they strengthened or lost their original function during species differentiation, resulting in olfactory differences between the two species.

Bone morphogenetic protein-2 (BMP-2) is a multi-functional growth factor belonging to the transforming growth factor β superfamily that has a broad range of activities that affect many different cell types. BMP-2 induces odontoblastic differentiation of human dental pulp cells (DPCs), but the underlying mechanism remains unclear. In this study, we investigated the potential role of the JNK mitogen-activated protein kinases (MAPK) pathway in BMP-2-induced odontoblastic differentiation of DPCs. The levels of phosphorylated and unphosphorylated JNK MAPK were quantified by Western blot analysis following treatment with BMP-2 and the JNK inhibitor SP600125. The role of JNK MAPK in the BMP-2-induced odontoblastic differentiation of DPCs was determined by measuring alkaline phosphatase (ALP) activity and by examining the expression of odontoblastic markers using quantitative real-time polymerase chain reaction analysis. The effect of JNK MAPK silencing on odontoblastic differentiation was also investigated. BMP-2 upregulated the phosphorylation of JNK in DPCs in a dose- and time-dependent manner. Early markers of odontoblastic differentiation, including ALP activity, osteopontin and dentin matrix protein-1, were not inhibited by the JNK inhibitor. However, the JNK inhibitor, SP600125, significantly inhibited late-stage differentiation of odontoblasts, including the gene expression of osteocalcin, dentin sialophosphoprotein and bone sialoprotein, and also reduced the formation of mineralized nodules in BMP-2-treated DPCs. Consistent with this observation, silencing of JNK MAPK also decreased late-stage odontoblastic differentiation. Taken together, these findings suggest that JNK activity is required for late-stage odontoblastic differentiation induced by BMP-2.