Deletion of Tmem268 in mice suppresses anti-infectious immune responses by downregulating CD11b signaling.

Transmembrane protein 268 (TMEM268) is a novel, tumor growth-related protein first reported by our laboratory. It interacts with the integrin subunit ß4 (ITGB4) and plays a positive role in the regulation of the ITGB4/PLEC signaling pathway. Here, we investigated the effects and mechanism of TMEM268 in anti-infectious immune response in mice. Tmem268 knockout in mice aggravated cecal ligation and puncture-induced sepsis, as evidenced by higher bacterial burden in various tissues and organs, congestion, and apoptosis. Moreover, Tmem268 deficiency in mice inhibited phagocyte adhesion and migration, thus decreasing phagocyte infiltration at the site of infection and complement-dependent phagocytosis. Further findings indicated that TMEM268 interacts with CD11b and inhibits its degradation via the endosome-lysosome pathway. Our results reveal a positive regulatory role of TMEM268 in ß2 integrin-associated anti-infectious immune responses and signify the potential value of targeting the TMEM268-CD11b signaling axis for the maintenance of immune homeostasis and immunotherapy for sepsis and related immune disorders.

Verticillium dahliae Vta3 promotes ELV1 virulence factor gene expression in xylem sap, but tames Mtf1-mediated late stages of fungus-plant interactions and microsclerotia formation.

Verticillium transcription activator of adhesion 3 (Vta3) is required for plant root colonization and pathogenicity of the soil-borne vascular fungus Verticillium dahliae. RNA sequencing identified Vta3-dependent genetic networks required for growth in tomato xylem sap. Vta3 affects the expression of more than 1,000 transcripts, including candidates with predicted functions in virulence and morphogenesis such as Egh16-like virulence factor 1 (Elv1) and Master transcription factor 1 (Mtf1). The genes encoding Elv1 and Mtf1 were deleted and their functions in V. dahliae growth and virulence on tomato (Solanum lycopersicum) plants were investigated using genetics, plant infection experiments, gene expression studies and phytohormone analyses. Vta3 contributes to virulence by promoting ELV1 expression, which is dispensable for vegetative growth and conidiation. Vta3 decreases disease symptoms mediated by Mtf1 in advanced stages of tomato plant colonization, while Mtf1 induces the expression of fungal effector genes and tomato pathogenesis-related protein genes. The levels of pipecolic and salicylic acids functioning in tomato defense signaling against (hemi-) biotrophic pathogens depend on the presence of MTF1, which promotes the formation of resting structures at the end of the infection cycle. In summary, the presence of VTA3 alters gene expression of virulence factors and tames the Mtf1 genetic subnetwork for late stages of plant disease progression and subsequent survival of the fungus in the soil.

Solution Combustion Synthesis of Submicron-Sized Titanium Niobium Oxide Anodes for High-Rate and Ultrastable Lithium-Ion Batteries.

Recently, the development of high-rate performance lithium-ion batteries is crucial for the development of next-generation energy storage systems. Nanoarchitecturing of the electrode material is a common strategy to improve the effective Li+ diffusion transport rate. However, this method often results in a reduction of volumetric energy density and battery stability. In this work, we propose a different strategy by synthesizing submicron-sized Ti2Nb10O29 (s-TNO) as a durable high-rate anode material using a facile and scalable solution combustion method, eliminating the dependence nanoarchitectures. The s-TNO electrode material exhibits a large tunnel structure and an excellent pseudocapacitive performance. The results show that this electrode material delivers a commendable reversible capacity of 238.7 mAh g-1 at 0.5 C and retains 78.2% of its capacity after 10,000 cycles at 10 C. This work provides a valuable guide for the synthesis of submicron-structured electrode materials using the solution combustion method, particularly for high-capacity, high-rate, and high-stability electrode materials.

Endosome associated trafficking regulator 1 promotes tumor growth and invasion of glioblastoma multiforme via inhibiting TNF signaling pathway.

PURPOSE: Endosome associated trafficking regulator 1 (ENTR1) is a novel endosomal protein, which can affect multiple cellular biological behavior by remodeling plasma membrane structures. However, little is known regarding its function and underlying mechanisms in glioblastoma multiforme. METHODS: Expression profile and clinical signature were obtained from The Public Database of human tumor. Immunohistochemical staining and western blotting assays were used to measure ENTR1 expression level. Human primary GBM tumor cells and human GBM cell lines A172, U87 and U251 were used to clarify the precise role of ENTR1. CCK-8 assays, wound healing and transwell invasion assays were designed to investigate cell viability, invasion and migration of GBM cells, respectively. Underlying molecular mechanisms of ENTR1 were determined via RNA-seq analysis. Tumor formation assay was used to validate the influence of ENTR1 in vivo. RESULTS: Compared with normal brain tissues, ENTR1 was highly expressed in gliomas and correlated with malignant grades of gliomas and poor overall survival time. The proliferation and invasion of GBM cells could be weaken and the sensitivity to temozolomide (TMZ) chemotherapy increased after knocking down ENTR1. Overexpression of ENTR1 could reverse this effect. RNA-seq analysis showed that tumor necrosis factor (TNF) signaling pathway might be a putative regulatory target of ENTR1. Tumor formation assay validated that ENTR1 was a significant factor in tumor growth. CONCLUSION: Our results indicated that ENTR1 played an important role in cell proliferation, invasion and chemotherapeutic sensitivity of GBM, suggesting that ENTR1 might be a novel prognostic marker and significant therapeutic target for GBM.

Emodin derivative E35 and its combination with autophagy inhibitor against acute myeloid leukemia cells in vitro and in vivo.

Acute myeloid leukemia (AML) is an aggressive hematopoietic malignancy with poor prognosis and high recurrence rate. The discovery of more effective therapeutic strategies for AML plays a crucial role. The present work showed that E35, a novel derivative of emodin, significantly inhibited cell proliferation and induced autophagy and apoptosis in AML cells. Treatment with E35 markedly induced Beclin-1, LC3-II, cleaved Caspase-9 and PARP, and suppressed mitogen-activated protein kinase (MAPK) pathway. E35 exposure evoked autophagic activity prior to apoptosis induction, and autophagy inhibition by 3-methyladenine (3-MA) dramatically increased E35-induced apoptosis in both AML cell lines and patient-derived AML cells. Nevertheless, study on AML xenograft model showed that the combination E35 with 3-MA exhibited much more inhibitory effects on leukemia cell growth in vivo. No obvious adverse reactions occurred in the xenograft animals administered E35 alone or its cotreatment with 3-MA. These findings suggest that E35 could exert anti-leukemia effects, and that the combination of E35 and autophagy inhibitor might prove a more highly efficient strategy for AML treatment.

The Vici Syndrome Protein EPG5 Is a Rab7 Effector that Determines the Fusion Specificity of Autophagosomes with Late Endosomes/Lysosomes.

Mutations in the human autophagy gene EPG5 cause the multisystem disorder Vici syndrome. Here we demonstrated that EPG5 is a Rab7 effector that determines the fusion specificity of autophagosomes with late endosomes/lysosomes. EPG5 is recruited to late endosomes/lysosomes by direct interaction with Rab7 and the late endosomal/lysosomal R-SNARE VAMP7/8. EPG5 also binds to LC3/LGG-1 (mammalian and C. elegans Atg8 homolog, respectively) and to assembled STX17-SNAP29 Qabc SNARE complexes on autophagosomes. EPG5 stabilizes and facilitates the assembly of STX17-SNAP29-VAMP7/8 trans-SNARE complexes, and promotes STX17-SNAP29-VAMP7-mediated fusion of reconstituted proteoliposomes. Loss of EPG5 activity causes abnormal fusion of autophagosomes with various endocytic vesicles, in part due to elevated assembly of STX17-SNAP25-VAMP8 complexes. SNAP25 knockdown partially suppresses the autophagy defect caused by EPG5 depletion. Our study reveals that EPG5 is a Rab7 effector involved in autophagosome maturation, providing insight into the molecular mechanism underlying Vici syndrome.

Omic horizon expression: a database of gene expression based on RNA sequencing data.

BACKGROUND: Gene expression profiles have important significance for gene expression characteristics and further functional studies. More attention has been given to the expression databases in humans and mice, but less attention has been given to rats, while rat models also play an irreplaceable role in biomedical experiments. RESULTS: To depict the rat gene expression profiles in mRNA expression levels, we analyzed over 2,700 RNA sequencing (RNA-Seq) samples from 48 tissues, 40 primary cell types and 25 cell lines; and then mapped them to the latest version of the rat genome reference, mRatBN7.2. Based on these datasets and reanalysis, we constructed a new database, the Omic Horizon Expression Database ( http://immudb.bjmu.edu.cn/expression.html ), which allows expressional profile query of over 25,000 rat genes based on non-redundant gene symbols. The database supports requests using gene symbols (or alias), Ensemble and Entrez gene IDs. Gene expression profiles can be queried in three categories: tissues, primary cells and cell lines. Application examples including expression profiling and comparison, as well as identification of novel rat genes, were illustrated to show the utility of the database. CONCLUSIONS: As an omic resource, the Omic Horizon Expression Database provides horizons of gene expression profiles across various tissues and cells, which greatly facilitates the identification of rat genes as well as functional clues.

Gut microbiome variations in Rhinopithecus roxellanae caused by changes in the environment.

BACKGROUND: The snub-nosed monkey (Rhinopithecus roxellanae) is an endangered animal species mainly distributed in China and needs to be protected. Gut microbiome is an important determinant of animal health and population survival as it affects the adaptation of the animals to different foods and environments under kinetic changes of intrinsic and extrinsic factors. Therefore, this study aimed to elucidate gut fecal microbiome profiles of snub-nosed monkeys affected by several extrinsic and intrinsic factors, including raising patterns (captive vs. wild), age, sex, and diarrheal status to provide a reference for making protection strategies. RESULTS: The 16S rRNA gene sequencing was firstly used to pre-check clustering of 38 fecal samples from the monkeys including 30 wild and 8 captive (5 healthy and 3 diarrheal) from three Regions of Shennongjia Nature Reserve, Hubei Province, China. Then the 24 samples with high-quality DNA from 18 wild and 6 captive (4 healthy and 2 diarrheal) monkeys were subjected to shotgun metagenomic sequencing to characterize bacterial gut microbial communities. We discovered that the raising pattern (captive and wild) rather than age and sex was the predominant factor attributed to gut microbiome structure and proportionality. Wild monkeys had significantly higher bacterial diversity and lower Bacteroidetes/Firmicutes ratios than captive animals. Moreover, the gut microbiomes in wild healthy monkeys were enriched for the genes involved in fatty acid production, while in captive animals, genes were enriched for vitamin biosynthesis and metabolism and amino acid biosynthesis from carbohydrate intermediates. Additionally, a total of 37 antibiotic resistant genes (ARG) types were detected. Unlike the microbiome diversity, the captive monkeys have a higher diversity of ARG than the wild animals. CONCLUSION: Taken together, we highlight the importance of self-reprogramed metabolism in the snub-nosed monkey gut microbiome to help captive and wild monkeys adapt to different intrinsic and extrinsic environmental change.

Environmental quality mediates the ecological dominance of cooperatively breeding birds.

Although social species as diverse as humans and ants are among the most abundant organisms on Earth, animals cooperate and form groups for many reasons. How these different reasons for grouping affect a species' ecological dominance remains unknown. Here we use a theoretical model to demonstrate that the different fitness benefits that animals receive by forming groups depend on the quality of their environment, which in turn impacts their ecological dominance and resilience to global change. We then test the model's key predictions using phylogenetic comparative analysis of >6500 bird species. As predicted, we find that cooperative breeders occurring in harsh and fluctuating environments have larger ranges and greater abundances than non-cooperative breeders, but cooperative breeders occurring in benign and stable environments do not. Using our model, we further show that social species living in harsh and fluctuating environments will be less vulnerable to climate change than non-social species.

Distribution, causal agents, and infection dynamic of emerging ink disease of sweet chestnut in Southern Switzerland.

Emerging diseases caused by both native and exotic pathogens represent a main threat to forest ecosystems worldwide. The two invasive soilborne pathogens Phytophthora cinnamomi and Phytophthora × cambivora are the causal agents of ink disease, which has been threatening Castanea sativa in Europe for several centuries and seems to be re-emerging in recent years. Here, we investigated the distribution, causal agents, and infection dynamics of ink disease in southern Switzerland. A total of 25 outbreaks were identified, 19 with only P. cinnamomi, 5 with only P. × cambivora, and 1 with both species. Dendrochronological analyses showed that the disease emerged in the last 20-30 years. Infected trees either died rapidly within 5-15 years post-infection or showed a prolonged state of general decline until death. Based on a generalized linear model, the local risk of occurrence of ink disease was increased by an S-SE aspect of the chestnut stand, the presence of a pure chestnut stand, management activities, the proximity of roads and buildings, and increasing annual mean temperature and precipitation. The genetic structure of the local P. cinnamomi population suggests independent introductions and local spread of the pathogen.

Unique Fluorescent Protein-Encoded Microbeads with Supramaximal Encoding Capacity and High Sensitivity for Multiplex Bioassays.

Fluorescence-encoded microbeads (FEBs) have been widely used as a critical component in multiplexed biomolecular assays. Here, we propose a simple, sustainable, low-cost, and safe strategy for preparing FEBs by assembling fluorescent proteins (FPs) onto magnetic microbeads (MBs) via chemical coupling. Combining the type of FP, the concentration of FP, and the size of the magnetic microbeads as encoding elements, an ultralarge encoding capacity with 506 barcodes was obtained. We demonstrate that the FP-based FEBs have good stability during long-term storage and tolerate the use of an organic solution. Multiplex detection of femtomolar ssDNA molecules was achieved via flow cytometry, and the detection procedure is simple and fast because it does not require amplification and washing strategies. The advantages of this advanced method for multiplex detections including high sensitivity, specificity, accuracy, repeatability, rapidity, and cost-effectiveness show a broad application prospect in basic and applied research fields such as disease diagnosis, food safety, environmental protection, proteomics, genomics, and drug screening.

Multi-parameter estimation resistant to the random polarization effect for coherent optical receivers.

Optical parameter estimation based on the data obtained by coherent optical receivers is critical for optical performance monitoring (OPM) and the stable operation of the receiver digital signal processing (DSP). A robust multi-parameter estimation is intricate due to the interference of various system effects. By resorting to the cyclostationary theory, we are able to formulate a chromatic dispersion (CD), frequency offset (FO), and optical signal-to-noise ratio (OSNR) joint estimation strategy that is resistant to the random polarization effect, including polarization mode dispersion (PMD) and polarization rotation. The method uses data directly after the DSP resampling and matched filtering. Both numerical simulation and field optical cable experiment validate our method.

Soluble expression and purification of recombinant bovine ferritin H-chain.

Ferritin is a potential medicine delivery vehicle and vaccine platform, and its efficient expression is a prerequisite for widespread application. This study introduces a soluble expression strategy for recombinant bovine ferritin heavy chain (rFTH) in a prokaryotic system and an improved protein purification method. The amplified rFTH gene was ligated into the prokaryotic expression vector pET30a. The recombinant vectors with the N-terminal His-tag(N-His) or C-terminal His-tag(C-His) were translated and expressed separately. The results showed that the solubility of rFTH with C-His was significantly higher than that with N-His. The expression of rFTH with C-His was attempted at 37 °C and 16 °C, respectively. The results showed that the proportion of soluble protein expressed at 37 °C was more than 90%, higher than that expressed at 16 °C. Then rFTH with C-His was purified successfully using anion exchange chromatography, modified PEG precipitation, and dialysis. The rFTH protein was characterized using SDS-PAGE, Native-PAGE, Western blot, transmission electron microscopy, and dynamic light scattering. The results demonstrated that the purified rFTH protein self-assembled into ferritin nanoparticles with a regular shape and uniform size. This study sheds new light on the soluble expression of ferritin and provides a foundation for the construction of bovine ferritin nanoparticle production platforms.

The complete mitochondrial genome of Ogmocotyle ailuri: gene content, composition and rearrangement and phylogenetic implications.

Trematodes of the genus Ogmocotyle are intestinal flukes that can infect a variety of definitive hosts, resulting in significant economic losses worldwide. However, there are few studies on molecular data of these trematodes. In this study, the mitochondrial (mt) genome of Ogmocotyle ailuri isolated from red panda (Ailurus fulgens) was determined and compared with those from Pronocephalata to investigate the mt genome content, genetic distance, gene rearrangements and phylogeny. The complete mt genome of O. ailuri is a typical closed circular molecule of 14 642 base pairs, comprising 12 protein-coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes and 2 non-coding regions. All genes are transcribed in the same direction. In addition, 23 intergenic spacers and 2 locations with gene overlaps were determined. Sequence identities and sliding window analysis indicated that cox1 is the most conserved gene among 12 PCGs in O. ailuri mt genome. The sequenced mt genomes of the 48 Plagiorchiida trematodes showed 5 types of gene arrangement based on all mt genome genes, with the gene arrangement of O. ailuri being type I. Phylogenetic analysis using concatenated amino acid sequences of 12 PCGs revealed that O. ailuri was closer to Ogmocotyle sikae than to Notocotylus intestinalis. These data enhance the Ogmocotyle mt genome database and provide molecular resources for further studies of Pronocephalata taxonomy, population genetics and systematics.

Identification of Differential Circular RNA Expression Profiles and Functional Networks in Human Macrophages Induced by Virulent and Avirulent Mycobacterium tuberculosis Strains.

Circular RNAs (circRNAs) are noncoding RNAs with diverse functions. However, most Mycobacterium tuberculosis (M.tb)-related circRNAs remain undiscovered. In this study, we infected THP-1 cells with virulent and avirulent M.tb strains and then sequenced the cellular circRNAs. Bioinformatic analysis predicted 58,009 circRNAs in all the cells. In total, 2035 differentially expressed circRNAs were identified between the M.tb-infected and uninfected THP-1 cells and 1258 circRNAs were identified in the virulent and avirulent M.tb strains. Further, the top 10 circRNAs were confirmed by Sanger sequencing, among which four circRNAs, namely circSOD2, circCHSY1, circTNFRSF21, and circDHTKD1, which were highly differentially expressed in infected cells compared with those in uninfected cells, were further confirmed by ring formation, specific primers, and RNase R digestion. Next, circRNA-miRNA-mRNA subnetworks were constructed, such as circDHTKD1/miR-660-3p/IL-12B axis. Some of the individual downstream genes, such as miR-660-3p and IL-12B, were previously reported to be associated with cellular defense against pathological processes induced by M.tb infection. Because macrophages are important immune cells and the major host cells of M.tb, these findings provide novel ideas for exploring the M.tb pathogenesis and host defense by focusing on the regulation of circRNAs during M.tb infection.

The LncRNA AK018453 regulates TRAP1/Smad signaling in IL-17-activated astrocytes: A potential role in EAE pathogenesis.

The pro-inflammatory cytokine interleukin 17 (IL-17), that is mainly produced by Th17 cells, has been recognized as a key regulator in multiple sclerosis (MS) and experimental autoimmune encephalomyelitis (EAE). Reactive astrocytes stimulated by proinflammatory cytokines including IL-17 are involved in blood brain barrier destruction, inflammatory cells infiltration and spinal cord injury. However, the role of long non-coding RNAs (lncRNAs) induced by IL-17 in the pathogenesis of MS and EAE remains unknown. Herein, we found that an IL-17-induced lncRNA AK018453 promoted TGF-ß receptor-associated protein 1 (TRAP1) expression and Smad-dependent signaling in mouse primary astrocytes. Knockdown of AK018453 significantly suppressed astrocytosis, attenuated the phosphorylation of Smad2/3, reduced NF-κB p65 and CBP/P300 binding to the TRAP1 promoter, and diminished pro-inflammatory cytokine production in the IL-17-treated astrocytes. AK018453 knockdown in astrocytes by a lentiviral vector in vivo dramatically inhibited inflammation and prevented the mice from demyelination in the spinal cord during the progression of EAE. Together, these results suggest that AK018453 regulates IL-17-dependent inflammatory response in reactive astrocytes and potentially promotes the pathogenesis of EAE via the TRAP1/Smad pathway. Targeting this pathway may have a therapeutic potential for intervening inflammatory demyelinating diseases.

Robust immune response stimulated by in situ injection of CpG/αOX40/cGAMP in αPD-1-resistant malignancy.

Recently, the emergence of immunotherapy has revolutionized traditional tumour treatment. However, effective treatments for patients exhibiting αPD-1 resistance are still lacking. In our study, a combination of cytosine-phosphate-guanine oligodeoxynucleotides (CpG-ODNs), anti-OX40 and cyclic guanosine monophosphate-adenosine monophosphate (cGAMP) injection in situ systematically generated a robust antitumour immune response in TC1 and B16 cells, which are αPD-1-resistant malignancies. More precisely, this method activates both adaptive and innate immunity. Additionally, in situ vaccination with CpG/αOX40/cGAMP fully activates the production of cytokines. However, the combination of αPD-1 does not improve the efficacy of triple therapy, prompting further questions. Collectively, the combination of CpG/αOX40/cGAMP causes the regression of various αPD-1-resistant tumours through the full mobilization of innate and adaptive immunity. In addition, we explored the therapeutic effect of triple therapy on the αPD-1-sensitive cell line CT26. The results showed that triple therapy could significantly enhance the therapeutic effect of αPD-1, and some mice even achieved complete tumour regression after the combined application of αPD-1 and triple treatment.

An emerging role for cyclic dinucleotide phosphodiesterase and nanoRNase activities in Mycoplasma bovis: Securing survival in cell culture.

Mycoplasmas are host-restricted prokaryotes with a nearly minimal genome. To overcome their metabolic limitations, these wall-less bacteria establish intimate interactions with epithelial cells at mucosal surfaces. The alarming rate of antimicrobial resistance among pathogenic species is of particular concern in the medical and veterinary fields. Taking advantage of the reduced mycoplasma genome, random transposon mutagenesis was combined with high-throughput screening in order to identify key determinants of mycoplasma survival in the host-cell environment and potential targets for drug development. With the use of the ruminant pathogen Mycoplasma bovis as a model, three phosphodiesterases of the DHH superfamily were identified as essential for the proliferation of this species under cell culture conditions, while dispensable for axenic growth. Despite a similar domain architecture, recombinant Mbov_0327 and Mbov_0328 products displayed different substrate specificities. While rMbovP328 protein exhibited activity towards cyclic dinucleotides and nanoRNAs, rMbovP327 protein was only able to degrade nanoRNAs. The Mbov_0276 product was identified as a member of the membrane-associated GdpP family of phosphodiesterases that was found to participate in cyclic dinucleotide and nanoRNA degradation, an activity which might therefore be redundant in the genome-reduced M. bovis. Remarkably, all these enzymes were able to convert their substrates into mononucleotides, and medium supplementation with nucleoside monophosphates or nucleosides fully restored the capacity of a Mbov_0328/0327 knock-out mutant to grow under cell culture conditions. Since mycoplasmas are unable to synthesize DNA/RNA precursors de novo, cyclic dinucleotide and nanoRNA degradation are likely contributing to the survival of M. bovis by securing the recycling of purines and pyrimidines. These results point toward proteins of the DHH superfamily as promising targets for the development of new antimicrobials against multidrug-resistant pathogenic mycoplasma species.

Prevalence of Theileria in cattle in China: A systematic review and meta-analysis.

Theileria, one of the causative agents of blood protozoan, has brought a huge economic loss to the cattle industry worldwide. However, the epidemiology of Theileria in Chinese cattle has not been systematically investigated. This comprehensive review aimed at investigating the prevalence of Theileria infection in cattle in China. A total of 48 published papers on Theileria infection in cattle in China (including data from 21,366 animals) from inception to October 8, 2021 met the inclusion standard after searching in five databases (Technology Periodical Database, Wan Fang Database, China National Knowledge Infrastructure, PubMed, and ScienceDirect). The pooled prevalence of Theileria in cattle in China was 32.4% identified by using a random effects model. The prevalence in Northeastern China (45.3%) was higher than that in other regions. In the sex subgroup, the prevalence of Theileria was higher in females (48.9%) than that in males (45.8%). The prevalence of Theileria was higher in cattle of free range (34.4%) compared with that of intensive farming (22.3%). The prevalence prior to 2013 (36.1%) was higher than that after 2013 (33.6%). Among three cattle species, dairy cows had the lowest prevalence (21.5%). The prevalence of Theileria (T.) annulata (22.2%) and T. sergenti (26.2%) was higher than other species of Theileria (T. buffeli: 17.5%, T. luwenshuni: 0.9%, T. orientalis: 15.5%, T. ovis: 0.21%, T. sinensis: 20.2%, T. uilenbergi: 6.2%, Others: 0.9%). We also analyzed the impact of different geographic factor subgroups (longitude, latitude, precipitation, temperature, humidity, and altitude) on the prevalence of Theileria in cattle. Among them, climatic factors of longitude, latitude, precipitation, humidity, temperature were associated with the prevalence of Theileria. These analyses suggested that Theileria was common in cattle in China. Targeted prevention programs based on geographic and climatic conditions in different areas may play an important role in reducing Theileria infection among cattle.

Detect differentially methylated regions using non-homogeneous hidden Markov model for bisulfite sequencing data.

DNA methylation plays an important role in many biological processes and diseases. With the rise of the whole genome bisulfite sequencing technique, aberrant methylation patterns can now be detected by comparing paired normal and disease samples at the single nucleotide level. We develop a novel Bayesian method for detecting differentially methylated regions from paired bisulfite sequencing data, and implement it as a R package called BSDMR. Based on a non-homogeneous hidden Markov model, BSDMR provides a better modeling strategy for the spatial correlation between CpG sites and takes into consideration the relationship between methylation signals from normal and disease samples. Simulations show that BSDMR performs well even under low read depth and has a smaller false discovery rates than existing methods. We also apply BSDMR to the colon cancer data from Gene Expression Omnibus. The detected DMRs are well supported by existing biomedical literatures.