Reverse transcriptase (RT) mutations contribute to hepatitis B virus resistance during antiviral therapy with nucleos(t)ide analogs. However, the composition of the RT quasispecies and their interactions during antiviral treatment have not yet been thoroughly defined. In this report, 10 patients from each of 3 different virological response groups, i.e., complete virological response, partial virological response and virological breakthrough, were selected from a multicenter trial of Telbivudine treatment. Variations in the drug resistance-related critical RT regions in 107 serial serum samples from the 30 patients were examined by ultra-deep sequencing. A total of 496,577 sequence reads were obtained, with an average sequencing coverage of 4,641X per sample. The phylogenies of the quasispecies revealed the independent origins of two critical quasispecies, i.e., the rtA181T and rtM204I mutants. Data analyses and theoretical modeling showed a cooperative-competitive interplay among the quasispecies. In particular, rtM204I mutants compete against other quasispecies, which eventually leads to virological breakthrough. However, in the absence of rtM204I mutants, synergistic growth of the drug-resistant rtA181T mutants with the wild-type quasispecies could drive the composition of the viral population into a state of partial virological response. Furthermore, we demonstrated that the frequency of drug-resistant mutations in the early phase of treatment is important for predicting the virological response to antiviral therapy.

In this study we examined the relationship between autophagy and apoptosis in diabetic rats after spinal cord injury (SCI), also we determined the role of autophagy in diabetes-aggravated neurological injury in vivo and in vitro. Our results showed that diabetes decreased the survival of neurons, promoted astrocytes proliferation, increased inflammatory cells infiltration and inhibited functional recovery after SCI. Diabetes was shown to confer increased activation of apoptotic pathways, along with an increase in autophagy; similar effects were also observed in vitro in neuronal PC12 cells. Treatment with rapamycin, an autophagy activator, partially abolished the adverse effect of diabetes, suggesting that diabetes may enhance neurological damage and suppress locomotor recovery after SCI, in addition to its effects on apoptosis and autophagy. In contrast, further stimulation of autophagy improved neurological function via inhibition of apoptosis. These results explained how diabetes exacerbates SCI in cellular level and suggested autophagy stimulation to be a new therapeutic strategy for diabetic SCI.

Leptospirosis, caused by pathogenic Leptospira spp., has recently been recognized as an emerging infectious disease worldwide. Despite its severity and global importance, knowledge about the molecular pathogenesis and virulence evolution of Leptospira spp. remains limited. Here we sequenced and analyzed 102 isolates representing global sources. A high genomic variability were observed among different Leptospira species, which was attributed to massive gene gain and loss events allowing for adaptation to specific niche conditions and changing host environments. Horizontal gene transfer and gene duplication allowed the stepwise acquisition of virulence factors in pathogenic Leptospira evolved from a recent common ancestor. More importantly, the abundant expansion of specific virulence-related protein families, such as metalloproteases-associated paralogs, were exclusively identified in pathogenic species, reflecting the importance of these protein families in the pathogenesis of leptospirosis. Our observations also indicated that positive selection played a crucial role on this bacteria adaptation to hosts. These novel findings may lead to greater understanding of the global diversity and virulence evolution of Leptospira spp.

Coxsackievirus A16 (CA16) is one of the major pathogens associated with human hand, foot, and mouth disease (HFMD) in the Asia-pacific region. Although CA16 infections are generally mild, severe neurological manifestations or even death has been reported. Studies on CA16 pathogenesis and vaccine development are severely hampered because the small animal models that are currently available show major limitations. In this study, gerbils (Meriones unguiculatus) were investigated for their suitability as an animal model to study CA16 pathogenesis and vaccine development. Our results showed that gerbils up to the age of 21 days were fully susceptible to CA16 and all died within five days post-infection. CA16 showed a tropism towards the skeletal muscle, spinal cord and brainstem of gerbils, and severe lesions, including necrosis, were observed. In addition, an inactivated CA16 whole-virus vaccine administrated to gerbils was able to provide full protection to the gerbils against lethal doses of CA16 strains. These results demonstrate that gerbils are a suitable animal model to study CA16 infection and vaccine development.

The fungal kingdom potentially has the most complex chitin synthase (CHS) gene family, but evolution of the fungal CHS gene family and its diversification to fulfill multiple functions remain to be elucidated. Here, we identified the full complement of CHSs from 231 fungal species. Using the largest dataset to date, we characterized the evolution of the fungal CHS gene family using phylogenetic and domain structure analysis. Gene duplication, domain recombination and accretion are major mechanisms underlying the diversification of the fungal CHS gene family, producing at least 7 CHS classes. Contraction of the CHS gene family is morphology-specific, with significant loss in unicellular fungi, whereas family expansion is lineage-specific with obvious expansion in early-diverging fungi. ClassV and ClassVII CHSs with the same domain structure were produced by the recruitment of domains PF00063 and PF08766 and subsequent duplications. Comparative analysis of their functions in multiple fungal species shows that the emergence of ClassV and ClassVII CHSs is important for the morphogenesis of filamentous fungi, development of pathogenicity in pathogenic fungi, and heat stress tolerance in Pezizomycotina fungi. This work reveals the evolution of the fungal CHS gene family, and its correlation with fungal morphogenesis and adaptation to ecological niches.

Cassava (Manihot esculenta Crantz) is an important tropical starchy root crop that is adapted to drought but extremely cold sensitive. A cold-tolerant, high-quality, and robust supply of cassava is urgently needed. Here, we clarify genome-wide distribution and classification of CCGG hemi-methylation and full-methylation, and detected 77 much candidate QTLsepi for cold stress and 103 much candidate QTLsepi for storage root quality and yield in 186 cassava population, generated by crossing two non-inbred lines with female parent KU50 and male parent SC124 (KS population). We developed amplified-fragment single nucleotide polymorphism and methylation (AFSM) genetic map in this population. We also constructed the AFSM QTL map, identified 260 much candidate QTL genes for cold stress and 301 much candidate QTL genes for storage root quality and yield, based on the years greenhouse and field trials. This may accounted for a significant amount of the variation in the key traits controlling cold tolerance and the high quality and yield of cassava.

More than 150 ginsenosides have been isolated and identified from Panax plants. Ginsenosides with different glycosylation degrees have demonstrated different chemical properties and bioactivity. In this study, we systematically cloned and characterized 46 UGT94 family UDP-glycosyltransferases (UGT94s) from a mixed Panax ginseng/callus cDNA sample with high amino acid identity. These UGT94s were found to catalyze sugar chain elongation at C3-O-Glc and/or C20-O-Glc of protopanaxadiol (PPD)-type, C20-O-Glc or C6-O-Glc of protopanaxatriol (PPT)-type or both C3-O-Glc of PPD-type and C6-O-Glc of PPT-type or C20-O-Glc of PPD-type and PPT-type ginsenosides with different efficiencies. We also cloned 26 and 51 UGT94s from individual P. ginseng and P. notoginseng plants, respectively; our characterization results suggest that there is a group of UGT94s with high amino acid identity but diverse functions or catalyzing activities even within individual plants. These UGT94s were classified into three clades of the phylogenetic tree and consistent with their catalytic function. Based on these UGT94s, we elucidated the biosynthetic pathway of a group of ginsenosides. Our present results reveal a series of UGTs involved in second sugar chain elongation of saponins in Panax plants, and provide a scientific basis for understanding the diverse evolution mechanisms of UGT94s among plants.

The association and potential role of the protein hormone adiponectin in autoimmune diseases causing musculoskeletal disorders, including rheumatoid arthritis (RA), are controversial. Conflicting results may arise from the influences of confounding factors linked to genetic backgrounds, disease stage, disease-modifying anti-rheumatic drugs and patients' metabolic characteristics. Here, we examined serum level of adiponectin and its relationship with disease activity score 28 with erythrocytes sedimentation rate (DAS28[ESR]) and Sharp score in a treatment-naïve Han Chinese RA population. This cross-sectional study enrolled 125 RA patients. Serum level of total adiponectin was assessed by enzyme-linked immunosorbent assay (ELISA). Other important clinical and laboratory parameters were collected from the hospital database. DAS28(ESR) was calculated according to the equation previously published. Sharp score was evaluated based on hands radiographs by an independent radiologist. The correlation between serum adiponectin level and DAS28(ESR) or the Sharp score was investigated by univariate and multivariable linear regression analyses, respectively. Multiple imputation by chained equations was used to account for missing data. Univariate analyses showed a significant positive correlation between DAS28(ESR) and age or C-reactive protein (CRP) (both p = 0.003), while serum adiponectin level was negatively correlated with DAS28(ESR) (p = 0.015). The negative correlation between adiponectin level and DAS28(ESR) remained true in multivariable analyses adjusted for confounders. In addition, the univariate analyses revealed positive correlations of Sharp score to disease duration (p < 0.001), CRP (p = 0.023) and ESR (p < 0.001). In the multivariable model adjusted for confounders, adiponectin was negatively correlated with Sharp score (p = 0.013). In this single-institution cross-sectional study, serum adiponectin level in treatment-naive RA patients is negatively correlated with DAS28(ESR) and the Sharp score after adjustment for prominent identified confounders. Serum adiponectin may be potentially useful for assessing disease activity and radiographic progression of RA.

Mycoplasma mycoides subsp. mycoides is the causative agent of contagious bovine pleuropneumonia. A pathogenic strain BEN-1 was isolated from bovine lung and underwent continuous passages in rabbits for 468 generations. During this process, the strain's strong virulence became weak and, gradually, it lost the ability to confer protective immunity in cattle but developed virulence in rabbits. In order to gain insight into the mechanisms behind the reduction in virulence and the loss of immunogenicity, we sequenced five representative strains of the BEN series, including the original strain (BEN-1), the strain generation that first acquired virulence in rabbits (BEN-50), the two vaccine strain generations (BEN-181 and BEN-326), and the strain generation showing the greatest loss of immunogenicity (BEN-468). The gene mutation rate in the four different propagation stages varied greatly, and over half of variations observed in each generation were removed during the propagation process. However, the variation maintained in the BEN-468 generation might contribute to its changes in virulence and immunogenicity. We thus identified 18 genes associated with host adaptation, six genes contributing to virulence in cattle, and 35 genes participating in conferring immunity in cattle. These findings might help us optimize the vaccine to obtain more effective immunization results.

The dysregulation of transcription factors has an important impact on the oncogenesis and tumor progression. Nonetheless, its functions in colorectal cancer metastasis are still unclear. In this study, four transcription factors (HNF4A, HSF1, MECP2 and RAD21) were demonstrated to be associated with the metastasis of colorectal cancer in both RNA and protein levels. To comprehensively explore the intrinsic mechanisms, we profiled the molecular landscape of these metastasis-related transcription factors from multiple perspectives. In particular, as the crucial factors affecting genome stability, both copy number variation and DNA methylation exerted their strengths on the expression of these transcription factors (except MECP2). Additionally, based on a series of bioinformatics analyses, putative long non-coding RNAs were identified as functional regulators. Besides that, rely on the ATAC-Seq and ChIP-Seq profiles, we detected the target genes regulated by each transcription factor in the active chromatin zones. Finally, we inferred the associations between the target genes by Bayesian networks and identified LMO7 and ARL8A as potential clinical biomarkers. Taken together, our research systematically characterized the regulatory cascades of HNF4A, HSF1, MECP2 and RAD21 in colorectal cancer metastasis.

Milk microbiota has a great influence on the safety and quality of dairy products. However, few studies have investigated the variations of bacterial composition in raw milk. In this study, raw milk samples were collected in 12 successive months, and their bacterial compositions were determined by 16 S rRNA gene sequencing. The highest diversity of bacterial composition was detected in June, while the lowest was in December. Firmicutes, Proteobacteria and Actinobacteria were the most abundant phyla and exhibited a counter-balanced relationship. Pseudomonas, Lactococcus and Acinetobacter were the most prevalent genera (>1%), and a tiny core microbiota (Acinetobacter and Pseudomonas) was observed. Temperature and humidity were the determining factors for most variation in bacterial compositions at both the phylum and genus levels. Higher abundances of Pseudomonas, Propionibacterium and Flavobacterium were correlated with low temperature. Furthermore, Pseudomonas/Propionibacterium and Lactobacillus/Bifidobacterium were two pairs of genera that had synergistic effects. Associations between the microbiota and milk quality parameters were analyzed. The abundances of Propionibacterium and Pseudoalteromonas were negatively correlated to total bacterial count, which meant that they helped to maintain milk quality, while a series of environmental microorganisms contributed to the spoilage of raw milk.

While chromatin remodeling mediated by post-translational modification of histone is extensively studied in carcinogenesis and cancer cell's response to chemotherapy and radiotherapy, little is known about the role of histone expression in chemoresistance. Here we report a novel chemoresistance mechanism involving histone H4 expression. Extended from our previous studies showing that concurrent blockage of the NF-κB and Akt signaling pathways sensitizes lung cancer cells to cisplatin-induced apoptosis, we for the first time found that knockdown of Akt1 and the NF-κB-activating kinase IKKβ cooperatively downregulated histone H4 expression, which increased cisplatin-induced apoptosis in lung cancer cells. The enhanced cisplatin cytotoxicity in histone H4 knockdown cells was associated with proteasomal degradation of RIP1, accumulation of cellular ROS and degradation of IAPs (cIAP1 and XIAP). The cisplatin-induced DNA-PK activation was suppressed in histone H4 knockdown cells, and inhibiting DNA-PK reduced expression of RIP1 and IAPs in cisplatin-treated cells. These results establish a novel mechanism by which NF-κB and Akt contribute to chemoresistance involving a signaling pathway consisting of histone H4, DNA-PK, RIP1 and IAPs that attenuates ROS-mediated apoptosis, and targeting this pathway may improve the anticancer efficacy of platinum-based chemotherapy.

Beauveria bassiana is an environmentally friendly alternative to chemical insecticides against various agricultural insect pests and vectors of human diseases. However, its application has been limited due to slow kill and sensitivity to abiotic stresses. Understanding of the molecular pathogenesis and physiological characteristics would facilitate improvement of the fungal performance. Loss-of-function mutagenesis is the most powerful tool to characterize gene functions, but it is hampered by the low rate of homologous recombination and the limited availability of selectable markers. Here, by combining the use of uridine auxotrophy as recipient and donor DNAs harboring auxotrophic complementation gene ura5 as a selectable marker with the blastospore-based transformation system, we established a highly efficient, low false-positive background and cost-effective CRISPR/Cas9-mediated gene editing system in B. bassiana. This system has been demonstrated as a simple and powerful tool for targeted gene knock-out and/or knock-in in B. bassiana in a single gene disruption. We further demonstrated that our system allows simultaneous disruption of multiple genes via homology-directed repair in a single transformation. This technology will allow us to study functionally redundant genes and holds significant potential to greatly accelerate functional genomics studies of B. bassiana.

Cutaneous melanoma (CM, hereafter referred to as melanoma) is a highly malignant tumor that typically undergoes early metastasis. Pyroptosis, as a special programmed cell death process that releases inflammatory factors and has been widely studied in tumors, but its role in melanoma has not been fully elucidated. In this study, we examined the relationship between pyroptosis and the prognosis of melanoma through bioinformatic analysis of RNA-sequencing data. Our results demonstrated that pyroptosis is a protective factor associated with melanoma prognosis. A higher pyroptosis score was associated with a more favorable overall survival. We used weighted gene co-expression networks analysis (WGCNA) to establish an effective prognosis model based on 12 pyroptosis-related genes. We then validated it in two independent cohorts. Furthermore, a nomogram combining clinicopathological characteristics and a pyroptosis-related gene signature (PGS) score was designed to effectively evaluate the prognosis of melanoma. Additionally, we analyzed the potential roles of pyroptosis in the tumor immune microenvironment and drug response. Interestingly, we found that the elevated infiltration of multiple immune cells, such as CD4+ T cells, CD8+ T cells, dendritic cells, and M1 macrophages, may be associated with the occurrence of pyroptosis. Pyroptosis was also related to a better response of melanoma to interferon-α, paclitaxel, cisplatin and imatinib. Through Spearman correlation analysis of the 12 pyroptosis-related genes and 135 chemotherapeutic agents in the Genomics of Drug Sensitivity in Cancer database, we identified solute carrier family 31 member 2 (SLC31A2) and collagen type 4 alpha 5 chain (COL4A5) as being associated with resistance to most of these drugs. In conclusion, this PGS is an effective and novelty prognostic indicator in melanoma, and also has an association with the melanoma immune microenvironment and melanoma treatment decision-making.

The Asian honeybee Apis cerana is one of two bee species that have been commercially kept with immense economic value. Here we present the analysis of genomic sequence and transcriptomic exploration for A. cerana as well as the comparative genomic analysis of the Asian honeybee and the European honeybee A. mellifera. The genome and RNA-seq data yield new insights into the behavioral and physiological resistance to the parasitic mite Varroa the evolution of antimicrobial peptides, and the genetic basis for labor division in A. cerana. Comparison of genes between the two sister species revealed genes specific to A. cerana, 54.5% of which have no homology to any known proteins. The observation that A. cerana displayed significantly more vigilant grooming behaviors to the presence of Varroa than A. mellifera in conjunction with gene expression analysis suggests that parasite-defensive grooming in A. cerana is likely triggered not only by exogenous stimuli through visual and olfactory detection of the parasite, but also by genetically endogenous processes that periodically activates a bout of grooming to remove the ectoparasite. This information provides a valuable platform to facilitate the traits unique to A. cerana as well as those shared with other social bees for health improvement.

The peptide hormone acyl-ghrelin and its receptor, GHSR1a, represent intriguing therapeutic targets due to their actions in metabolic homeostasis and reward activity. However, this pleotropic activity makes it difficult to intervene in this system without inducing unwanted effects. Thus, it is desirable to identify passive and active regulatory mechanisms that allow differentiation between functional domains. Anatomical restriction by the blood brain barrier represents one major passive regulatory mechanism. However, it is likely that the ghrelin system is subject to additional passive mechanisms that promote independent regulation of orexigenic behavior and reward processing. By applying acyl-ghrelin sequestering antibodies, it was determined that peripheral sequestration of acyl-ghrelin is sufficient to blunt weight gain, but not cocaine rewarding effects. However, both weight gain and reward-associated behaviors were shown to be blocked by direct antagonism of GHSR1a. Overall, these data indicate that GHSR1a effects on reward are independent from peripheral acyl-ghrelin binding, whereas centrally-mediated alteration of energy storage requires peripheral acyl-ghrelin binding. This demonstration of variable ligand-dependence amongst functionally-distinct GHSR1a populations is used to generate a regulatory model for functional manipulation of specific effects when attempting to therapeutically target the ghrelin system.

MicroRNAs (miRNAs) play important roles in the fibrosis of systemic sclerosis (SSc). However, the underlying miRNA-mRNA regulatory network is not fully understood. A systemic investigation of the role of miRNAs would be very valuable for increasing our knowledge of the pathogenesis of SSc. Here, we combined miRNA and mRNA expression profiles and bioinformatics analyses and then performed validation experiments. we identified 21 miRNAs and 2698 mRNAs that were differentially expressed in SSc. Among these, 17 miRNAs and their 33 target mRNAs (55 miRNA-mRNA pairs) were involved in Toll-like receptor, transforming growth factor β and Wnt signalling pathways. Validation experiments revealed that miR-146b, miR-130b, miR-21, miR-31 and miR-34a levels were higher whereas miR-145 levels were lower in SSc skin tissues and fibroblasts, normal fibroblasts and endothelial cells that were stimulated with SSc serum. ACVR2B, FZD2, FZD5 and SOX2 levels were increased in SSc skin fibroblasts, normal fibroblasts and endothelial cells that were stimulated with SSc serum. We did not identify any negative correlations among these miRNA-mRNA pairs. miR-21 was specifically expressed at higher levels in SSc serum. Six miRNAs and 4 mRNAs appear to play important roles in the pathogenesis of SSc are worth investigating in future functional studies.

Many infertile men are the victims of spermatogenesis disorder. However, conventional clinical test could not provide efficient information on the causes of spermatogenesis disorder and guide the doctor how to treat it. More effective diagnosis and treating methods could be developed if the key genes that regulate spermatogenesis were determined. Many works have been done on animal models, while there are few works on human beings due to the limited sample resources. In current work, testis tissues were obtained from 27 patients with obstructive azoospermia via surgery. The combination of Fluorescence Activated Cell Sorting and Magnetic Activated Cell Sorting was chosen as the efficient method to sort typical germ cells during spermatogenesis. RNA Sequencing was carried out to screen the change of transcriptomic profile of the germ cells during spermatogenesis. Differential expressed genes were clustered according to their expression patterns. Gene Ontology annotation, pathway analysis, and Gene Set Enrichment Analysis were carried out on genes with specific expression patterns and the potential key genes such as HOXs, JUN, SP1, and TCF3 which were involved in the regulation of spermatogenesis, with the potential value serve as molecular tools for clinical purpose, were predicted.

The maintenance of the heart rhythm and the conduction of excitatory signals require changing excitatory signals via electrical activity and coordination by communication between working and conductive cardiomyocytes. Understanding how the ventricular conduction system is established provides novel insights into the pathophysiological progress of cardiac arrhythmias. However, the major hurdle in this field is the lack of a specific genetic tool that targets the Purkinje fibres of the ventricular conduction system and no other types of cardiomyocytes or coronary vessels. Here, we generated a Sema3a-CreERT2 knock-in mouse line to test its specificity for genetically labelled Purkinje fibres. We found that Sema3a was expressed in the subendocardial layer of the trabecular myocardium in the embryonic heart and was restricted to the Purkinje fibres in the adult heart. A fate mapping study based on the Sema3a-CreERT2 line revealed that the Sema3a+ cardiomyocytes were restricted to the fate of Purkinje fibres in the perinatal but not the embryonic stage. Collectively, our study provides a new genetic tool, i.e., Sema3a-CreERT2, for studying the molecular mechanisms that regulate the function of Purkinje fibres.

The need for high-affinity, SARS-CoV-2-specific monoclonal antibodies (mAbs) is critical in the face of the global COVID-19 pandemic, as such reagents can have important diagnostic, research, and therapeutic applications. Of greatest interest is the ~ 300 amino acid receptor binding domain (RBD) within the S1 subunit of the spike protein because of its key interaction with the human angiotensin converting enzyme 2 (hACE2) receptor present on many cell types, especially lung epithelial cells. We report here the development and functional characterization of 29 nM-affinity mouse SARS-CoV-2 mAbs created by an accelerated immunization and hybridoma screening process. Differing functions, including binding of diverse protein epitopes, viral neutralization, impact on RBD-hACE2 binding, and immunohistochemical staining of infected lung tissue, were correlated with variable gene usage and sequence.