Search

Abstract The nucleosome is one of the hallmarks of eukaryotes, a dynamic platform that supports many critical functions in eukaryotic cells. Here, we engineer the in vivo assembly of the nucleosome core in the model bacterium Escherichia coli. We show that bacterial chromosome DNA and eukaryotic histones can assemble in vivo to form nucleosome complexes with many features resembling those found in eukaryotes. The formation of nucleosomes in E. coli was visualized with atomic force microscopy and using tripartite split green fluorescent protein. Under a condition that moderate histones expression was induced at 1 µM IPTG, the nucleosome-forming bacterium is viable and has sustained growth for at least 110 divisions in longer-term growth experiments. It exhibits stable nucleosome formation, a consistent transcriptome across passages, and reduced growth fitness under stress conditions. In particular, the nucleosome arrays in E. coli genic regions have profiles resembling those in eukaryotic cells. The observed compatibility between the eukaryotic nucleosome and the bacterial chromosome machinery may reflect a prerequisite for bacteria-archaea union, providing insight into eukaryogenesis and the origin of the nucleosome.

Abstract Background Crohn's disease (CD) is a disease characterized by intestinal immune dysfunction, often accompanied by metabolic abnormalities. Disturbances in lactate metabolism have been found in the intestine of patients with CD, but studies on the role of lactate and related Lactylation in the pathogenesis of CD are still unknown. Methods We identified the core genes associated with Lactylation by downloading and merging three CD-related datasets (GSE16879, GSE75214, and GSE112366) from the GEO database, and analyzed the functions associated with the hub genes and the correlation between their expression levels and immune infiltration through comprehensive analysis. We explored the Lactylation levels of different immune cells using single-cell data and further analyzed the differences in Lactylation levels between inflammatory and non-inflammatory sites. Results We identified six Lactylation-related hub genes that are highly associated with CD. Further analysis revealed that these six hub genes were highly correlated with the level of immune cell infiltration. To further clarify the effect of Lactylation on immune cells, we analyzed single-cell sequencing data of immune cells from inflammatory and non-inflammatory sites in CD patients and found that there were significant differences in the levels of Lactylation between different types of immune cells, and that the levels of Lactylation were significantly higher in immune cells from inflammatory sites. Conclusions These results suggest that Lactylation-related genes and their functions are closely associated with changes in inflammatory cells in CD patients.

Abstract Cancer immunotherapies have achieved unprecedented success in clinic, but they remain largely ineffective in some major types of cancer, such as colorectal cancer with microsatellite stability (MSS CRC). It is therefore important to study tumor microenvironment of resistant cancers for developing new intervention strategies. In this study, we identify a metabolic cue that determines the unique immune landscape of MSS CRC. Through secretion of distal cholesterol precursors, which directly activate RORγt, MSS CRC cells can polarize T cells toward Th17 cells that have well-characterized pro-tumor functions in colorectal cancer. Analysis of large human cancer cohorts revealed an asynchronous pattern of the cholesterol biosynthesis in MSS CRC, which is responsible for the abnormal accumulation of distal cholesterol precursors. Inhibiting the cholesterol biosynthesis enzyme Cyp51, by pharmacological or genetic interventions, reduced the levels of intratumoral distal cholesterol precursors and suppressed tumor progression through a Th17-modulation mechanism in preclinical MSS CRC models. Our study therefore reveals a novel mechanism of cancer–immune interaction and an intervention strategy for the difficult-to-treat MSS CRC.

Abstract Single-molecule detection and phasing of A-to-I RNA editing events remain an unresolved problem. Long-read and PCR-free nanopore native RNA sequencing offers a great opportunity for direct RNA editing detection. Here, we develop a neural network model, DeepEdit, that not only recognizes A-to-I editing events in single reads of Oxford Nanopore direct RNA sequencing, but also resolves the phasing of RNA editing events on transcripts. We illustrate the robustness of DeepEdit by applying it to Schizosaccharomyces pombe and Homo sapiens transcriptome data. We anticipate DeepEdit to be a powerful tool for the study of RNA editing from a new perspective.

Coxsackievirus A10 (CVA10) has recently emerged as one of the major causative agents of hand, foot, and mouth disease. CVA10 may also cause a variety of complications. No approved vaccine or drug is currently available for CVA10. The residues of CVA10 critical for viral attachment, infectivity and in vivo pathogenicity have not been identified by experiment. Here, we report the identification of CVA10 residues important for binding to cellular receptor KREMEN1. We identified VP2 N142 as a key receptor-binding residue by screening of CVA10 mutants resistant to neutralization by soluble KREMEN1 protein. The receptor-binding residue N142 is exposed on the canyon rim but highly conserved in all naturally occurring CVA10 strains, which provides a counterexample to the canyon hypothesis. Residue N142 when mutated drastically reduced receptor-binding activity, resulting in decreased viral attachment and infection in cell culture. More importantly, residue N142 when mutated reduced viral replication in limb muscle and spinal cord of infected mice, leading to lower mortality and less severe clinical symptoms. Additionally, residue N142 when mutated could decrease viral binding affinity to anti-CVA10 polyclonal antibodies and a neutralizing monoclonal antibody and render CVA10 resistant to neutralization by the anti-CVA10 antibodies. Overall, our study highlights the essential role of VP2 residue N142 of CVA10 in the interactions with KREMEN1 receptor and neutralizing antibodies and viral virulence in mice, facilitating the understanding of the molecular mechanisms of CVA10 infection and immunity. Our study also provides important information for rational development of antibody-based treatment and vaccines against CVA10 infection.

Anti-melanoma differentiation-associated gene 5-positive dermatomyositis is a severe autoimmune disease with largely unknown aetiology. Here authors identify the key immune cell types that contribute to the disease phenotype and implicate type I interferons signalling in the patho-mechanism.

Abstract Models developed using Nanopore direct RNA sequencing data from in vitro synthetic RNA with all adenosine replaced by N 6-methyladenosine (m6A) are likely distorted due to superimposed signals from saturated m6A residues. Here, we develop a neural network, DENA, for m6A quantification using the sequencing data of in vivo transcripts from Arabidopsis. DENA identifies 90% of miCLIP-detected m6A sites in Arabidopsis and obtains modification rates in human consistent to those found by SCARLET, demonstrating its robustness across species. We sequence the transcriptome of two additional m6A-deficient Arabidopsis, mtb and fip37-4, using Nanopore and evaluate their single-nucleotide m6A profiles using DENA.

Abstract Background Virus screening and viral genome reconstruction are urgent and crucial for the rapid identification of viral pathogens, i.e., tracing the source and understanding the pathogenesis when a viral outbreak occurs. Next-generation sequencing (NGS) provides an efficient and unbiased way to identify viral pathogens in host-associated and environmental samples without prior knowledge. Despite the availability of software, data analysis still requires human operations. A mature pipeline is urgently needed when thousands of viral pathogen and viral genome reconstruction samples need to be rapidly identified. Results In this paper, we present a rapid and accurate workflow to screen metagenomics sequencing data for viral pathogens and other compositions, as well as enable a reference-based assembler to reconstruct viral genomes. Moreover, we tested our workflow on several metagenomics datasets, including a SARS-CoV-2 patient sample with NGS data, pangolins tissues with NGS data, Middle East Respiratory Syndrome (MERS)-infected cells with NGS data, etc. Our workflow demonstrated high accuracy and efficiency when identifying target viruses from large scale NGS metagenomics data. Our workflow was flexible when working with a broad range of NGS datasets from small (kb) to large (100 Gb). This took from a few minutes to a few hours to complete each task. At the same time, our workflow automatically generates reports that incorporate visualized feedback (e.g., metagenomics data quality statistics, host and viral sequence compositions, details about each of the identified viral pathogens and their coverages, and reassembled viral pathogen sequences based on their closest references). Conclusions Overall, our system enabled the rapid screening and identification of viral pathogens from metagenomics data, providing an important piece to support viral pathogen research during a pandemic. The visualized report contains information from raw sequence quality to a reconstructed viral sequence, which allows non-professional people to screen their samples for viruses by themselves (Additional file 1).

Abstract Background Single-cell RNA sequencing (scRNA-seq) is the most widely used technique to obtain gene expression profiles from complex tissues. Cell subsets and developmental states are often identified via differential gene expression patterns. Most of the single-cell tools utilized highly variable genes to annotate cell subsets and states. However, we have discovered that a group of genes, which sensitively respond to environmental stimuli with high coefficients of variation (CV), might impose overwhelming influences on the cell type annotation. Result In this research, we developed a method, based on the CV-rank and Shannon entropy, to identify these noise genes, and termed them as “sensitive genes”. To validate the reliability of our methods, we applied our tools in 11 single-cell data sets from different human tissues. The results showed that most of the sensitive genes were enriched pathways related to cellular stress response. Furthermore, we noticed that the unsupervised result was closer to the ground-truth cell labels, after removing the sensitive genes detected by our tools. Conclusion Our study revealed the prevalence of stochastic gene expression patterns in most types of cells, compared the differences among cell marker genes, housekeeping genes (HK genes), and sensitive genes, demonstrated the similarities of functions of sensitive genes in various scRNA-seq data sets, and improved the results of unsupervised clustering towards the ground-truth labels. We hope our method would provide new insights into the reduction of data noise in scRNA-seq data analysis and contribute to the development of better scRNA-seq unsupervised clustering algorithms in the future.

This work presents the application of both uniaxial and bi-axial in-plane loads to unidirectional carbon fiber reinforced polymer laminates using a newly developed Arcan fixture, which is a reliable experimental set-up to obtain a uniform shear stress field in a butterfly specimen. The set-up can be used for both damage model validation and parameters identification at various fiber orientations while using the same specimens. A sequential damage study was completed to highlight the influence of diffused damage induced in pure shear on the fiber direction tensile behavior of the laminate. This was accomplished by applying load on the specimens in two steps: (i) the pure shear step and then unloading at approximately 70% of the shear failure strength, (ii) in the tensile step until final failure. A clear drop in the tensile behavior of the laminate was observed by the diffused damage induced in the first loading step of pure shear. The experimental study is also supplemented with numerical simulations using a nonlinear elasto-plastic coupled damage constitutive law by employing Puck’s failure theory for mesodamage activation. In addition to the damage pattern, the non-linear mechanical behavior in shear is predicted and found in good correspondence with the experimental results.

In view of the problems that existing underground positioning algorithm based on UWB is complex and solution value is not global optimal, an improved trilateral positioning algorithm based on UWB was proposed. The algorithm is based on DW1000 chip and trilateral positioning algorithm, adopts bilateral two-way ranging method. Optimization objective function is set by solving the optimal solution of binary quadratic equation, and coordinate is two-dimensional, so as to reduce complexity of the algorithm. The invalid data generated in measurement process were screened out, effective search area is determined, and global traversal is conducted in the effective area to find global optimal positioning coordinates, which solves the problem that solution value is not global optimal. At the same time, in order to further improve positioning accuracy, appropriate measures are taken to increase the number of reference measurement nodes. The experimental results show that the average error of the proposed positioning algorithm is lower than that of the weighted trilateral positioning algorithm and the Taylor algorithm, and the positioning accuracy is improved obviously. Meanwhile, the algorithm has strong practicability, and measurement accuracy can be effectively improved by increasing the number of reference nodes.

Here, the authors show that Zika virus (ZIKV) infection induces Th1-like Tfh cells that depend on T-bet for their development and are essential for class switching of ZIKV-specific IgG2c antibodies and maintenance of long-term neutralizing antibody responses.

Abstract Background Apical periodontitis (AP) is essentially an inflammatory disease of microbial etiology primarily caused by infection of the pulp and root canal system. Variation of the bacterial communities caused by AP, as well as their changes responding to dental therapy, are of utmost importance to understand the pathogensis of the apical periodontitis and establishing effective antimicrobial therapeutic strategies. This study aims to uncover the composition and diversity of microbiota associated to the root apex to identify the relevant bacteria highly involved in AP, with the consideration of root apex samples from the infected teeth (with/without root canal treatment), healthy teeth as well as the healthy oral. Methods Four groups of specimens are considered, the apical part of root from diseased teeth with and without root canal treatment, and wisdom teeth extracted to avoid being impacted (tooth healthy control), as well as an additional healthy oral control from biofilm of the buccal mucosa. DNA was extracted from these specimens and the microbiome was examined through focusing on the V3-V4 hypervariable region of the 16S rRNA gene using sequencing on Illumina MiSeq platform. Composition and diversity of the bacterial community were tested for individual samples, and between-group comparisons were done through differential analysis to identify the significant changes. Results We observed reduced community richness and diversity in microbiota samples from diseased teeth compared to healthy controls. Through differential analysis between AP teeth and healthy teeth, we identified 49 OTUs significantly down-regulated as well as 40 up-regulated OTUs for AP. Conclusion This study provides a global view of the microbial community of the AP associated cohorts, and revealed that AP involved not only bacteria accumulated with a high abundance, but also those significantly reduced ones due to microbial infection.

Abstract Background This study aimed to explore the molecular mechanism and clinical relevance of iguratimod in the regulation of human B cell terminal differentiation. Methods An in vitro human antibody-secreting cell (ASC) differentiation system was established to test the effect of iguratimod. B cell phenotype and key transcription factors (TFs) relevant to ASC differentiation were analyzed through flow cytometry and qPCR. The COX-2 activity was measured by enzyme immunoassay (EIA). RNA sequencing was used to identify potential targets of iguratimod. We enrolled six treatment-naive rheumatoid arthritis (RA) patients whose blood samples were collected for phenotypic and molecular studies along with 12-week iguratimod monotherapy. Results Iguratimod inhibited human ASC generation without affecting B cell activation and proliferation. Iguratimod showed only weak COX-2 activity. Gene set enrichment analysis (GSEA) identified that protein kinase C (PKC) pathway was targeted by iguratimod which was confirmed by PKC activity detection. Furthermore, early growth response 1 (EGR1), a target of PKC and a non-redundant TF for ASC differentiation, was found to be the most downregulated gene in iguratimod-treated B cells. Lastly, iguratimod monotherapy decreased peripheral ASCs and was associated with improved disease activity. The expression of major ASC-related TFs, including EGR1, was similarly downregulated in patient blood samples. Conclusions Iguratimod inhibits ASC differentiation both in vitro and in RA patients. Our study suggests that PKC/EGR1 axis, rather than COX-2, is critically involved in the inhibitory effect by iguratimod on human ASC differentiation. Iguratimod could have a broader application to treat B cell-related autoimmune diseases in clinics.

Though an increasing number of studies have indicated the existence of an interaction between the virome and human health or disease, the specific role of these plasma viral components remains largely unsolved. We provide evidence here that an altered plasma virome profile is associated with different immune statuses of HIV-1 infection.

BackgroundThe pathogenesis of COVID-19 emerges as complex, with multiple factors leading to injury of different organs. Some of the studies on aspects of SARS-CoV-2 cell entry and innate immunity have produced seemingly contradictory claims. In this situation, a comprehensive comparative analysis of a large number of related datasets from several studies could bring more clarity, which is imperative for therapy development.MethodsWe therefore performed a comprehensive comparative study, analyzing RNA-Seq data of infections with SARS-CoV-2, SARS-CoV and MERS-CoV, including data from different types of cells as well as COVID-19 patients. Using these data, we investigated viral entry routes and innate immune responses.Results and ConclusionFirst, our analyses support the existence of cell entry mechanisms for SARS and SARS-CoV-2 other than the ACE2 route with evidence of inefficient infection of cells without expression of ACE2; expression of TMPRSS2/TPMRSS4 is unnecessary for efficient SARS-CoV-2 infection with evidence of efficient infection of A549 cells transduced with a vector expressing human ACE2. Second, we find that innate immune responses in terms of interferons and interferon simulated genes are strong in relevant cells, for example Calu3 cells, but vary markedly with cell type, virus dose, and virus type.

Servo test system was used to obtain stress-strain curves of Beishan granite under different confining pressure, and the loading process was jointly controlled by means of axial pressure and circumferential deformation.Variation characteristics of energy and characteristic stress of Beishan granite were researched.Characteristic stress are linear increase with confining pressure increasing, crack initiation stress, damage stress and peak strength increased 3.30 MPa, 6.74 MPa, 7.68 MPa respectively along with confining pressure increased 1 MPa.Total energy and elastic energy corresponding to characteristic stress are greatly increase with confining pressure increasing.When confining pressure increased form 2 MPa to 30 MPa, elastic energy corresponding to initiation stress, damage stress and peak stress increased 5.2 times, 5.8 times and 5.1 times respectively, corresponding total energy increased 4.9 times, 5.7 times and 4.7 times respectively.The increment of dissipation energy increase with increasing of confining pressure, and energy at the crack unstable propagation stage was significant affected by confining pressure.Energy storage coefficient was defined to estimate rock stability around the peak strength, effect of stress state coefficient on energy storage was further analyzed.

Hematopoietic stem cells (HSCs) in the bone marrow are quiescent, but are activated in response to stress. Here, the authors show that loss of Med23 leads to greater activation and enhanced myeloid potential of HSCs in response to stress, also Med23 maintains stemness gene expression and suppresses myeloid genes.

Abstract Background RNA editing is an important mechanism that expands the diversity and complexity of genetic codes. The conversions of adenosine (A) to inosine (I) and cytosine (C) to uridine (U) are two prominent types of RNA editing in animals. The roles of RNA editing events have been implicated in important biological pathways. Cellular RNA editing activity in response to influenza A virus infection has not been fully characterized in human and avian hosts. This study was designed as a big data analysis to investigate the role and response of RNA editing in epithelial cells during the course of infection with various subtypes of influenza A viruses. Results Using a bioinformatics pipeline modified from our previous study, we characterized the profiles of A-to-I and C-to-U RNA editing events in human epithelial cells during the course of influenza A virus infection. Our results revealed a striking diversity of A-to-I RNA editing activities in human epithelial cells in responses to different subtypes of influenza A viruses. The infection of H1N1 and H3N2 significantly up-regulated normalized A-to-I RNA editing levels in human epithelial cells, whereas that of H5N1 did not change it and H7N9 infection significantly down-regulated normalized A-to-I editing level in A549 cells. Next, the expression levels of ADAR and APOBEC enzymes responsible for A-to-I and C-to-U RNA editing during the course of virus infection were examined. The increase of A-to-I RNA editing activities in infections with some influenza A viruses (H1N1 and H3N2) is linked to the up-regulation of ADAR1 but not ADAR2. Further, the pattern recognition receptors of human epithelial cells infected with H1N1, H3N2, H5N1 and H7N9 were examined. Variable responsive changes in gene expression were observed with RIG-I like receptors and Toll like receptors. Finally, the effect of influenza A virus infection on cellular RNA editing activity was also analyzed in avian hosts. Conclusion This work represents the first comprehensive study of cellular RNA editing activity in response to different influenza A virus infections in human and avian hosts, highlighting the critical role of RNA editing in innate immune response and the pathogenicity of different subtypes of influenza A viruses.

Abstract Background Adenosine-to-Inosine (A-to-I) RNA editing is catalyzed by the adenosine deaminase acting on RNA (ADAR) family of enzymes, which induces alterations in mRNA sequence. It has been shown that A-to-I RNA editing events are of significance in the cell’s innate immunity and cellular response to viral infections. However, whether RNA editing plays a role in cellular response to microorganism/fungi infection has not been determined. Candida albicans, one of the most prevalent human pathogenic fungi, usually act as a commensal on skin and superficial mucosal, but has been found to cause candidiasis in immunosuppression patients. Previously, we have revealed the up-regulation of A-to-I RNA editing activity in response to different types of influenza virus infections. The current work is designed to study the effect of microorganism/fungi infection on the activity of A-to-I RNA editing in infected hosts. Results We first detected and characterized the A-to-I RNA editing events in oral epithelial cells (OKF6) and primary human umbilical vein endothelial cells (HUVEC), under normal growth condition or with C. albicans infection. Eighty nine thousand six hundred forty eight and 60,872 A-to-I editing sites were detected in normal OKF6 and HUVEC cells, respectively. They were validated against the RNA editing databases, DARNED, RADAR, and REDIportal with 50, 80, and 80% success rates, respectively. While over 95% editing sites were detected in Alu regions, among the rest of the editing sites in non repetitive regions, the majority was located in introns and UTRs. The distributions of A-to-I editing activity and editing depth were analyzed during the course of C. albicans infection. While the normalized editing levels of common editing sites exhibited a significant increase, especially in Alu regions, no significant change in the expression of ADAR1 or ADAR2 was observed. Second, we performed further analysis on data from in vivo mouse study with C. albicans infection. One thousand one hundred thirty three and 955 A-to-I editing sites were identified in mouse tongue and kidney tissues, respectively. The number of A-to-I editing events was much smaller than in human epithelial or endothelial cells, due to the lack of Alu elements in mouse genome. Furthermore, during the course of C. albicans infection we observed stable level of A-to-I editing activity in 131 and 190 common editing sites in the mouse tongue and kidney tissues, and found no significant change in ADAR1 or ADAR2 expression (with the exception of ADAR2 displaying a significant increase at 12 h after infection in mouse kidney tissue before returning to normal). Conclusions This work represents the first comprehensive analysis of A-to-I RNA editome in human epithelial and endothelial cells. C. albicans infection of human epithelial and endothelial cells led to the up-regulation of A-to-I editing activities, through a mechanism different from that of viral infections in human hosts. However, the in vivo mouse model with C. albicans infection did not show significant changes in A-to-I editing activities in tongue and kidney tissues. The different results in the mouse model were likely due to the presence of more complex in vivo environments, e.g. circulation and mixed cell types.

Abstract Background Innate immunity provides first line of defense against viral infections. The interactions between hosts and influenza A virus and the response of host innate immunity to viral infection are critical determinants for the pathogenicity or virulence of influenza A viruses. This study was designed to investigate global changes of gene expression and detailed responses of innate immune systems in human and avian hosts during the course of infection with various subtypes of influenza A viruses, using collected and self-generated transcriptome sequencing data from human bronchial epithelial (HBE), human tracheobronchial epithelial (HTBE), and A549 cells infected with influenza A virus subtypes, namely H1N1, H3N2, H5N1 HALo mutant, and H7N9, and from ileum and lung of chicken and quail infected with H5N1, or H5N2. Results We examined the induction of various cytokines and chemokines in human hosts infected with different subtypes of influenza A viruses. Type I and III interferons were found to be differentially induced with each subtype. H3N2 caused abrupt and the strongest response of IFN-β and IFN-λ, followed by H1N1 (though much weaker), whereas H5N1 HALo mutant and H7N9 induced very minor change in expression of type I and III interferons. Similarly, differential responses of other innate immunity-related genes were observed, including TMEM173, MX1, OASL, IFI6, IFITs, IFITMs, and various chemokine genes like CCL5, CX3CL1, and chemokine (C-X-C motif) ligands, SOCS (suppressors of cytokine signaling) genes. Third, the replication kinetics of H1N1, H3N2, H5N1 HALo mutant and H7N9 subtypes were analyzed, H5N1 HALo mutant was found to have the highest viral replication rate, followed by H3N2, and H1N1, while H7N9 had a rate similar to that of H1N1 or H3N2 though in different host cell type. Conclusion Our study illustrated the differential responses of innate immunity to infections of different subtypes of influenza A viruses. We found the influenza viruses which induced stronger innate immune responses replicate slower than those induces weaker innate immune responses. Our study provides important insight into links between the differential innate immune responses from hosts and the pathogenicity/ virulence of different subtypes of influenza A viruses.

The properties and data structure of distance histogram algorithm of three-dimensional space data were analyzed, and a general computing method based on graphics processor and a high performance computing method based on FPGA were proposed. The calculation method based on graphics processor is used to implement single instruction multiple data parallel optimization of distance histogram algorithm of three-dimensional space data; the calculation method based on FPGA can realize hardware block optimization, and achieve the optimal matching of the hardware structure of the algorithm. The experimental results show that using the calculation method based on graphics processor can make the algorithm reach performance acceleration with an average of 18 times, and the computing method based on FPGA can make the algorithm reach performance acceleration with an average of 30 times, which greatly improve data processing ability of the algorithm.

Orifice throttle is a common way of throttling in the air static pressure guide rail,In order to study the gas film pressure distribution of the air static pressure guide rail,a gas thrust bearing gas film pressure distribution test bench was developed. In this paper,the orifice throttling is an example,based on the law of the gas film pressure decline,the gas film is divided into the pressure driven area and the Newton friction area. The gas flow pattern in the gas film support area is simulated by Fluent and the flow rate and pressure are calculated. The experimental results are in good agreement with the numerical analysis of the gas pressure distribution. It is concluded that the gas film is divided into the pressure driven area and the Newton friction zone along the length of the gas film,which verifies the validity of the pressure partition theory.

Regulatory T (Treg) cells play central roles in maintaining immune homeostasis and self-tolerance. However, the molecular mechanisms underlying Treg cell homeostasis and suppressive function are still not fully understood. Here, we report that the deletion of another P subfamily members of the forkhead box (Foxp) subfamily member Foxp1 in Treg cells led to increased numbers of activated Treg (aTreg) cells at the expense of quiescent Treg cells, and also resulted in impaired Treg suppressive function. Mice with Foxp1-deficient Treg cells developed spontaneous inflammatory disease with age; they also had more severe inflammatory disease in colitis and experimental autoimmune encephalomyelitis (EAE) models. Mechanistically, we found that Foxp1 bound to the conserved noncoding sequence 2 (CNS2) element of the Foxp3 locus and helped maintain Treg suppressive function by stabilizing the Foxp3 expression. Furthermore, we found that Foxp1 and Foxp3 coordinated the regulation of cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) expression levels. Taken together, our study demonstrates that Foxp1 plays critical roles in both maintaining Treg cell quiescence during homeostasis and regulating Treg suppressive function.

Bainitangceng in the Yunnan and Guizhou area has a close relationship with manganese deposits. Previous geochemical studies focused on the whole rock of siliceous limestones rather than siliceous components. The question of whether or not whole rocks can really reflect the source of siliceous components needs further study. In this study, Bainitangceng' siliceous limestone in Yunnan and Guizhou area was selected and HCl was used to remove the carbonate minerals in siliceous limestone to acquire simple composition siliceous components. Trace elements in siliceous components were determined by ICP-MS/OES. The origin of siliceous components is discussed by comparing the trace elemental composition of the whole rock of siliceous limestone and siliceous components after removing carbonate. Results show that Th/U ratios (0.02-0.05) in siliceous components are lower than those (0.08-0.2) of the whole rock of siliceous limestone. Siliceous components after removing the carbonate minerals have Eu/Eu* values of 0.86 to 1.54 with positive Eu anomaly, different from Eu/Eu* values (0.87-0.93) of the whole rock of siliceous limestone with negative Eu anomaly. Both whole rock of siliceous limestone and siliceous components are characterized by enriching the light rare earth elements, but the proportion of light rare earth elements in siliceous components is larger than that in whole rock. The above characteristics indicate that the siliceous component of ‘Bainitangceng’ may have been derived from hydrothermal fluid, which is important for explaining the formation of Permian manganese deposits.

When lubricating gas film of aerostatic guide way worked in micron level,internal airflow was affected by rarefied gas,the pressure and flow rate along flow direction was different to the result of calculating. In the gas film height direction,due to the characteristic scale decreased,the effect of gas molecule was intensified,traditional assumption of ignored vertical pressure gradient no longer founded,not only interior air film present obvious stratified trend,but velocity slip phenomenon was more significant. Simulating flow patterns and calculating its velocity and pressure in gas film by LAMMPS（ Large-scale Molecular Massively Parallel Simulator）,LBM（ Lattice Boltzmann Method） and Fluent in this paper,came to the relevant conclusion: the gas film in the longitudinal direction was divided into pressure-driven area and Newton friction area,and it also was divided into near wall layer,thin layer,continuous flow layer along its height direction,velocity slip mainly occurred in the near wall layer and thin layer of pressure-driven area. Gas stratification and velocity slip present in pressure-driven and Newton friction areas;Gas Stratification and velocity slip was enhanced with increasing gas flow rate,gas film stratification in Newton friction zone was gradually weakening,slippage speed gradually degenerated into near-wall boundary layer.

Elucidating protective immunity against HCV is important for the development of a preventative vaccine. We hypothesize that spontaneous resolution of acute HCV infection offers clue to protective immune responses, and that DAA therapy affects the quality and quantity of HCV-specific T cell responses. To test these hypotheses, we performed T cell epitope mapping in 111 HCV-infected individuals including 61 chronically HCV-1b (CHC-1b) infected, 24 chronically HCV-2a (CHC-2a) infected and 26 spontaneously recovered (SPR) patients with 376 overlapping peptides covering the entire HCV polyprotein. Selected T cell epitopes were then used to evaluate T cell responses in another 22 chronically HCV-1b infected patients on DAA therapy. Results showed that SPR had better HCV-specific T cell responses than CHC, as manifested by higher response rate, greater magnitude and broader epitope coverage. In addition, SPR recognized novel epitopes in Core, E1, E2, NS4B, NS5A regions that were not present in the CHC. Furthermore, during the first 24 weeks of DAA therapy, there was no functional immune reconstitution of HCV-specific T cells. These results indicate that T cell responses may be a correlate of protection. Therefore, effective preventative vaccines should elicit a robust T cell response. Although various DAA regimens efficiently cleared viruses from the blood of HCV-infected patients, there was no contemporaneous early T cell immune reconstitution, suggesting that early treatment is needed for preserving the functions of HCV-specific T cells.

The hydrolytic deamination of adenosine to inosine (A-to-I editing) in precursor mRNA induces variable gene products at the post-transcription level. How and to what extent A-to-I RNA editing diversifies transcriptome is not fully characterized in the evolution, and very little is known about the selective constraints that drive the evolution of RNA editing events. Here we present a study on A-to-I RNA editing, by generating a global profile of A-to-I editing for a phylogeny of seven Drosophila species, a model system spanning an evolutionary timeframe of approximately 45 million years. Of totally 9281 editing events identified, 5150 (55.5%) are located in the coding sequences (CDS) of 2734 genes. Phylogenetic analysis places these genes into 1,526 homologous families, about 5% of total gene families in the fly lineages. Based on conservation of the editing sites, the editing events in CDS are categorized into three distinct types, representing events on singleton genes (type I), and events not conserved (type II) or conserved (type III) within multi-gene families. While both type I and II events are subject to purifying selection, notably type III events are positively selected, and highly enriched in the components and functions of the nervous system. The tissue profiles are documented for three editing types, and their critical roles are further implicated by their shifting patterns during holometabolous development and in post-mating response. In conclusion, three A-to-I RNA editing types are found to have distinct evolutionary dynamics. It appears that nervous system functions are mainly tested to determine if an A-to-I editing is beneficial for an organism. The coding plasticity enabled by A-to-I editing creates a new class of binary variations, which is a superior alternative to maintain heterozygosity of expressed genes in a diploid mating system.

Bacterial pathogens continue to threaten public health worldwide today. Identification of bacterial virulence factors can help to find novel drug/vaccine targets against pathogenicity. It can also help to reveal the mechanisms of the related diseases at the molecular level. With the explosive growth in protein sequences generated in the postgenomic age, it is highly desired to develop computational methods for rapidly and effectively identifying virulence factors according to their sequence information alone. In this study, based on the protein-protein interaction networks from the STRING database, a novel network-based method was proposed for identifying the virulence factors in the proteomes of UPEC 536, UPEC CFT073, P. aeruginosa PAO1, L. pneumophila Philadelphia 1, C. jejuni NCTC 11168 and M. tuberculosis H37Rv. Evaluated on the same benchmark datasets derived from the aforementioned species, the identification accuracies achieved by the network-based method were around 0.9, significantly higher than those by the sequence-based methods such as BLAST, feature selection and VirulentPred. Further analysis showed that the functional associations such as the gene neighborhood and co-occurrence were the primary associations between these virulence factors in the STRING database. The high success rates indicate that the network-based method is quite promising. The novel approach holds high potential for identifying virulence factors in many other various organisms as well because it can be easily extended to identify the virulence factors in many other bacterial species, as long as the relevant significant statistical data are available for them.

Influenza virus, which spreads around the world in seasonal epidemics and leads to large numbers of deaths every year, has several ribonucleoproteins in the central core of the viral particle. These viral ribonucleoproteins can specifically bind the conserved 3' and 5' caps of the viral RNAs with responsibility for replication and transcription of the viral RNA in the nucleus of infected cells. A fundamental question of most importance is that how the cap-binding proteins in the influenza virus discriminates between capped RNAs and non-capped ones. To get an answer, we performed molecular dynamics simulations and free energy calculations on the influenza A virus PB2 subunit, an important component of the RNP complexes, with a cap analog m7GTP. Our calculations showed that some key residues in the active site, such as Arg355, His357, Glu361 as well as Gln406, could offer significant hydrogen bonding and hydrophobic interactions with the guanine ring of the cap analog m7GTP to form an aromatic sandwich mechanism for the cap recognition and positioning in the active site. Subsequently, we applied this idea to a virtual screening procedure and identified 5 potential candidates that might be inhibitors against the PB2 subunit. Interestingly, 2 candidates Cpd1 and Cpd2 have been already reported to have inhibitory activities to the influenza virus cap-binding proteins. Further calculation also showed that they had comparatively higher binding affinities to the PB2 subunit than that of m7GTP. We believed that our findings could give an atomic insight into the deeper understanding of the cap recognition and binding mechanism, providing useful information for searching or designing novel drugs against influenza viruses.

Lactobacillus delbrueckii subsp. bulgaricus (Lb. bulgaricus) is an important species of Lactic Acid Bacteria (LAB) used for cheese and yogurt fermentation. The genome of Lb. bulgaricus 2038, an industrial strain mainly used for yogurt production, was completely sequenced and compared against the other two ATCC collection strains of the same subspecies. Specific physiological properties of strain 2038, such as lysine biosynthesis, formate production, aspartate-related carbon-skeleton intermediate metabolism, unique EPS synthesis and efficient DNA restriction/modification systems, are all different from those of the collection strains that might benefit the industrial production of yogurt. Other common features shared by Lb. bulgaricus strains, such as efficient protocooperation with Streptococcus thermophilus and lactate production as well as well-equipped stress tolerance mechanisms may account for it being selected originally for yogurt fermentation industry. Multiple lines of evidence suggested that Lb. bulgaricus 2038 was genetically closer to the common ancestor of the subspecies than the other two sequenced collection strains, probably due to a strict industrial maintenance process for strain 2038 that might have halted its genome decay and sustained a gene network suitable for large scale yogurt production.

Pyrrolidone carboxylic acid (PCA) is formed during a common post-translational modification (PTM) of extracellular and multi-pass membrane proteins. In this study, we developed a new predictor to predict the modification sites of PCA based on maximum relevance minimum redundancy (mRMR) and incremental feature selection (IFS). We incorporated 727 features that belonged to 7 kinds of protein properties to predict the modification sites, including sequence conservation, residual disorder, amino acid factor, secondary structure and solvent accessibility, gain/loss of amino acid during evolution, propensity of amino acid to be conserved at protein-protein interface and protein surface, and deviation of side chain carbon atom number. Among these 727 features, 244 features were selected by mRMR and IFS as the optimized features for the prediction, with which the prediction model achieved a maximum of MCC of 0.7812. Feature analysis showed that all feature types contributed to the modification process. Further site-specific feature analysis showed that the features derived from PCA's surrounding sites contributed more to the determination of PCA sites than other sites. The detailed feature analysis in this paper might provide important clues for understanding the mechanism of the PCA formation and guide relevant experimental validations.