Surgical resection remains a critical treatment option for many patients with primary and secondary hepatic neoplasms. Extended hepatectomy (eHx) may be required for some patients with large tumors, which may cause liver failure and death. Partial hepatectomy (pHx) and eHx mouse models were constructed, liver tissues were sampled at 18, 36, and 72 h posthepatectomy. Transcriptome and metabolome analyses were employed to explore the different potential mechanisms in regeneration and injury between pHx and eHx. The results showed that eHx was associated with more severe liver injury and lower survival rates than pHx. Transcriptomics data showed there were 1842, 2129, and 1277 differentially expressed genes (DEGs) in eHx and 962, 1305, and 732 DEGs in pHx at 18, 36, and 72 h posthepatectomy, respectively, compared with the those in the sham groups. Compared with pHx, the number of DEGs in the eHx group reached a maximum of 230 at 18 h after surgery and decreased sequentially to 87 and 43 at 36 and 72 h. Metabolomics analysis identified a total of 1399 metabolites, and 48 significant differentially produced metabolites (DPMs) were screened between eHx and pHx. Combined analysis of DEGs and DPMs indicated that cholesterol metabolism and insulin resistance may be two important pathways for liver regeneration and mouse survival postextended hepatectomy. Our results showed the global influence of pHx and eHx on the transcriptome and metabolome in mouse liver, and revealed cholesterol metabolism and insulin resistance pathways might be involved in regeneration post-pHx and -eHx.

BACKGROUND: Ligilactobacillus salivarius has been frequently isolated from the gut microbiota of humans and domesticated animals and has been studied as a candidate probiotic. Badger (Meles meles) is known as a "generalist" species that consumes complex foods and exhibits tolerance and resistance to certain pathogens, which can be partly attributed to the beneficial microbes such as L. salivarius in the gut microbiota. However, our understanding of the beneficial traits and genomic features of badger-originated L. salivarius remains elusive. RESULTS: In this study, nine L. salivarius strains were isolated from wild badgers' feces, one of which exhibited good probiotic properties. Complete genomes of the nine L. salivarius strains were generated, and comparative genomic analysis was performed with the publicly available complete genomes of L. salivarius obtained from humans and domesticated animals. The strains originating from badgers harbored a larger genome, a higher number of protein-coding sequences, and functionally annotated genes than those originating from humans and chickens. The pan-genome phylogenetic tree demonstrated that the strains originating from badgers formed a separate clade, and totally 412 gene families (12.6% of the total gene families in the pan-genome) were identified as genes gained by the last common ancestor of the badger group. The badger group harbored significantly more gene families responsible for the degradation of complex carbohydrate substrates and production of polysaccharides than strains from other hosts; many of these were acquired by gene gain events. CONCLUSIONS: A candidate probiotic and nine L. salivarius complete genomes were obtained from the badgers' gut microbiome, and several beneficial genes were identified to be specifically present in the badger-originated strains that were gained in the evolution. Our study provides novel insights into the adaptation of L. salivarius to the intestinal habitat of wild badgers and provides valuable strain and genome resources for the development of L. salivarius as a probiotic.

BACKGROUND: Due to the influence of extreme weather, the environment in China's main cotton-producing areas is prone to drought stress conditions, which affect the growth and development of cotton and lead to a decrease in cotton yield. RESULTS: In this study, 188 upland cotton germplasm resources were phenotyped for data of 8 traits (including 3 major yield traits) under drought conditions in three environments for two consecutive years. Correlation analysis revealed significant positive correlations between the three yield traits. Genetic analysis showed that the estimated heritability of the seed cotton index (SC) under drought conditions was the highest (80.81%), followed by that of boll weight (BW) (80.64%) and the lint cotton index (LC) (70.49%) With genome-wide association study (GWAS) analysis, a total of 75 quantitative trait loci (QTLs) were identified, including two highly credible new QTL hotspots. Three candidate genes (Gh_D09G064400, Gh_D10G261000 and Gh_D10G254000) located in the two new QTL hotspots, QTL51 and QTL55, were highly expressed in the early stage of fiber development and showed significant correlations with SC, LC and BW. The expression of three candidate genes in two extreme materials after drought stress was analyzed by qRT-PCR, and the expression of these two materials in fibers at 15, 20 and 25 DPA. The expression of these three candidate genes was significantly upregulated after drought stress and was significantly higher in drought-tolerant materials than in drought-sensitive materials. In addition, the expression levels of the three candidate genes were higher in the early stage of fiber development (15 DPA), and the expression levels in drought-tolerant germplasm were higher than those in drought-sensitive germplasm. These three candidate genes may play an important role in determining cotton yield under drought conditions. CONCLUSIONS: This study is helpful for understanding the regulatory genes affecting cotton yield under drought conditions and provides germplasm and candidate gene resources for breeding high-yield cotton varieties under these conditions.

BACKGROUND: Colletotrichum camelliae, one of the most important phytopathogenic fungi infecting tea plants (Camellia sinensis), causes brown blight disease resulting in significant economic losses in yield of some sensitive cultivated tea varieties. To better understand its phytopathogenic mechanism, the genetic information is worth being resolved. RESULTS: Here, a high-quality genomic sequence of C. camelliae (strain LT-3-1) was sequenced using PacBio RSII sequencing platform, one of the most advanced Three-generation sequencing platforms and assembled. The result showed that the fungal genomic sequence is 67.74 Mb in size (with the N50 contig 5.6 Mb in size) containing 14,849 putative genes, of which about 95.27% were annotated. The data revealed a large class of genomic clusters potentially related to fungal pathogenicity. Based on the Pathogen Host Interactions database, a total of 1698 genes (11.44% of the total ones) were annotated, containing 541 genes related to plant cell wall hydrolases which is remarkably higher than those of most species of Colletotrichum and others considered to be hemibiotrophic and necrotrophic fungi. It's likely that the increase in cell wall-degrading enzymes reflects a crucial adaptive characteristic for infecting tea plants. CONCLUSION: Considering that C. camelliae has a specific host range and unique morphological and biological traits that distinguish it from other species of the genus Colletotrichum, characterization of the fungal genome will improve our understanding of the fungus and its phytopathogenic mechanism as well.

To combat drought stress in rice, a major threat to global food security, three major quantitative trait loci for 'yield under drought stress' (qDTYs) were successfully exploited in the last decade. However, their molecular basis still remains unknown. To understand the role of secondary regulation by miRNA in drought stress response and their relation, if any, with the three qDTYs, the miRNA dynamics under drought stress was studied at booting stage in two drought tolerant (Sahbaghi Dhan and Vandana) and one drought sensitive (IR 20) cultivars. In total, 53 known and 40 novel differentially expressed (DE) miRNAs were identified. The primary drought responsive miRNAs were Osa-MIR2919, Osa-MIR3979, Osa-MIR159f, Osa-MIR156k, Osa-MIR528, Osa-MIR530, Osa-MIR2091, Osa-MIR531a, Osa-MIR531b as well as three novel ones. Sixty-one target genes that corresponded to 11 known and 4 novel DE miRNAs were found to be co-localized with the three qDTYs, out of the 1746 target genes identified. We could validate miRNA-mRNA expression under drought for nine known and three novel miRNAs in eight different rice genotypes showing varying degree of tolerance. From our study, Osa-MIR2919, Osa-MIR3979, Osa-MIR528, Osa-MIR2091-5p and Chr01_11911S14Astr and their target genes LOC_Os01g72000, LOC_Os01g66890, LOC_Os01g57990, LOC_Os01g56780, LOC_Os01g72834, LOC_Os01g61880 and LOC_Os01g72780 were identified as the most promising candidates for drought tolerance at booting stage. Of these, Osa-MIR2919 with 19 target genes in the qDTYs is being reported for the first time. It acts as a negative regulator of drought stress tolerance by modulating the cytokinin and brassinosteroid signalling pathway.

BACKGROUND: Environmental changes are expected to intensify in the future. The invasion of toxic plants under environmental changes may change herbivore feeding environments. Herbivores living long-term in toxic plant-feeding environments will inevitably ingest plant secondary metabolites (PSMs), and under different feeding environments are likely to have unique protection mechanisms that support improved adaptation to PSMs in their habitat. We aimed to compare different subterranean herbivore population responses and adaptations to toxic plants to unveil their feeding challenges. RESULTS: Here, we investigated the adaptive capacity of the liver in two geographically separated populations of plateau zokors (Eospalax baileyi) before and after exposure to the toxic plant Stellera chamaejasme (SC), at the organ, biochemical, and transcriptomic levels. The results showed no significant liver granules or inflammatory reactions in the Tianzhu (TZ) population after the SC treatment. The transaminase level in the TZ population was significantly lower than that in the Luqu population. Transcriptome analysis revealed that the TZ population exhibited interactions with other detoxification metabolic pathways by oxytocin pathway-associated genes, including diacylglycerol lipase alpha (Dagla), calcium/calmodulin dependent protein kinase II Alpha (Camk2a), and CD38 molecule (Cd38). The phase II process of liver drug metabolism increased to promote the rate of metabolism. We found that alternative splicing (AS) and the expression of the cyclin D (Ccnd1) gene interact-a TZ population hallmark-reduced liver inflammatory responses. CONCLUSION: Our study supports the detoxification limitation hypothesis that differences in liver detoxification metabolism gene expression and AS are potential factors in herbivore adaptation to PSMs and may be a strategy of different herbivore populations to improve toxic plant adaptability.

BACKGROUND: The harsh conditions of high-altitude environments are known to drive the evolution of physiological and morphological traits in endothermic animals. These conditions are expected to result in the adaptive evolution of protein coding genes encoded in mitochondrial genomes that are vital for the oxidative phosphorylation pathway. In this study, we formally tested for signatures of adaptive evolution on mitochondrial protein coding genes in Tapirus pinchaque and other odd-toed ungulates inhabiting high-elevation environments. RESULTS: The AT-rich mitochondrial genome of T. pinchaque is 16,750 bp long. A phylomitogenomic analysis supports the monophyly of the genus Tapirus and families in the Perissodactyla. The ratio of non-synonymous to synonymous substitutions demonstrated that all mitochondrial genes undergo purifying selection in T. pinchaque and other odd ungulates living at high elevations. Over this negative background selection, Branch Models suggested that cox3 and nad6 might be undergoing stronger purifying selection than other mitochondrial protein coding genes. Furthermore, Site Models suggested that one and four sites in nad2 and nad5, respectively, could be experiencing positive selection. However, these results were supported by Likelihood Ratio Tests but not Bayesian Empirical Bayes posterior probabilities. Additional analyses (in DataMonkey) indicated a relaxation of selection strength in nad6, evidence of episodic diversifying selection in cob, and revealed episodic positive/diversifying selection signatures for two sites in nad1, and one site each in nad2 and nad4. CONCLUSION: The mitochondrial genome of T. pinchaque is an important genomic resource for conservation of this species and this study contributes to the understanding of adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates inhabiting high-altitude environments.

BACKGROUND: The incidence of kidney disease caused by thyroid cancer is rising worldwide. Observational studies cannot recognize whether thyroid cancer is independently associated with kidney disease. We performed the Mendelian randomization (MR) approach to genetically investigate the causality of thyroid cancer on immunoglobulin A nephropathy (IgAN). METHODS AND RESULTS: We explored the causal effect of thyroid cancer on IgAN by MR analysis. Fifty-two genetic loci and single nucleotide polymorphisms were related to thyroid cancer. The primary approach in this MR analysis was the inverse variance weighted (IVW) method, and MR‒Egger was the secondary method. Weighted mode and penalized weighted median were used to analyze the sensitivity. In this study, the random-effect IVW models showed the causal impact of genetically predicted thyroid cancer across the IgAN risk (OR, 1.191; 95% CI, 1.131-1.253, P < 0.001). Similar results were also obtained in the weighted mode method (OR, 1.048; 95% CI, 0.980-1.120, P = 0.179) and penalized weighted median (OR, 1.185; 95% CI, 1.110-1.264, P < 0.001). However, the MR‒Egger method revealed that thyroid cancer decreased the risk of IgAN, but this difference was not significant (OR, 0.948; 95% CI, 0.855-1.051, P = 0.316). The leave-one-out sensitivity analysis did not reveal the driving influence of any individual SNP on the association between thyroid cancer and IgAN. CONCLUSION: The IVW model indicated a significant causality of thyroid cancer with IgAN. However, MR‒Egger had a point estimation in the opposite direction. According to the MR principle, the evidence of this study did not support a stable significant causal association between thyroid cancer and IgAN. The results still need to be confirmed by future studies.

BACKGROUND: Gut microbiota plays a significant role in host survival, health, and diseases; however, compared to other livestock, research on the gut microbiome of donkeys is limited. RESULTS: In this study, a total of 30 donkey samples of rectal contents from six regions, including Shigatse, Changdu, Yunnan, Xinjiang, Qinghai, and Dezhou, were collected for metagenomic sequencing. The results of the species annotation revealed that the dominant phyla were Firmicutes and Bacteroidetes, and the dominant genera were Bacteroides, unclassified_o_Clostridiales (short for Clostridiales) and unclassified_f_Lachnospiraceae (short for Lachnospiraceae). The dominant phyla, genera and key discriminators were Bacteroidetes, Clostridiales and Bacteroidetes in Tibet donkeys (Shigatse); Firmicutes, Clostridiales and Clostridiales in Tibet donkeys (Changdu); Firmicutes, Fibrobacter and Tenericutes in Qinghai donkeys; Firmicutes, Clostridiales and Negativicutes in Yunnan donkeys; Firmicutes, Fibrobacter and Fibrobacteres in Xinjiang donkeys; Firmicutes, Clostridiales and Firmicutes in Dezhou donkeys. In the functional annotation, it was mainly enriched in the glycolysis and gluconeogenesis of carbohydrate metabolism, and the abundance was the highest in Dezhou donkeys. These results combined with altitude correlation analysis demonstrated that donkeys in the Dezhou region exhibited strong glucose-conversion ability, those in the Shigatse region exhibited strong glucose metabolism and utilization ability, those in the Changdu region exhibited a strong microbial metabolic function, and those in the Xinjiang region exhibited the strongest ability to decompose cellulose and hemicellulose. CONCLUSION: According to published literature, this is the first study to construct a dataset with multi-regional donkey breeds. Our study revealed the differences in the composition and function of gut microbes in donkeys from different geographic regions and environmental settings and is valuable for donkey gut microbiome research.

BACKGROUND: Substantial variation in male fertility is regularly observed in farmed Arctic charr. However, detailed investigations of its fluctuation during a reproductive season and across years are lacking. Furthermore, information about the effect of underlying genetic factors influencing sperm quality is scarce. The current study focused on seasonal and age-related factors that may affect sperm quality characteristics in males reared in natural and delayed photoperiods. Animals were sampled three times for two consecutive years, and sperm quality parameters were recorded using a computer-assisted sperm analysis (CASA) system. Thereafter, high-throughput sequencing technologies were applied, aiming to identify genomic regions related to the variation of sperm quality throughout the reproductive season. RESULTS: An across-season variation in the recorded sperm quality parameters was evident. Overall, 29% and 42% of males from the natural and delayed spawning groups had a highly variable total progressive motility. Males at four years of age showed significantly higher sperm motility and velocities during the early October and November recordings compared to the following year when the same animals were five years of age. On the other hand, the opposite was observed regarding sperm concentration during the last sampling. A genome-wide FST scan detected SNP differentiation among males with high and low variability in total progressive motility (PM) on eight chromosomes (FST > 0.17), Genome wide windows with the highest FST contained SNPs in proximity (within 250 kb up- and downstream distance) to 16 genes with sperm quality biological functions in mammalian species. CONCLUSION: Our findings provide a detailed view of seasonal, age-related, and genetic effects on sperm quality and can be used to guide decisions on broodstock selection and hatchery management.

BACKGROUND: Ubiquitination controls almost all cellular processes. The dysregulation of ubiquitination signals is closely associated with the initiation and progression of multiple diseases. However, there is little comprehensive research on the interaction and potential function of ubiquitination regulators (UBRs) in spermatogenesis and cancer. METHODS: We systematically characterized the mRNA and protein expression of UBRs across tissues and further evaluated their roles in testicular development and spermatogenesis. Subsequently, we explored the genetic alterations, expression perturbations, cancer hallmark-related pathways, and clinical relevance of UBRs in pan-cancer. RESULTS: This work reveals heterogeneity in the expression patterns of UBRs across tissues, and the expression pattern in testis is the most distinct. UBRs are dynamically expressed during testis development, which are critical for normal spermatogenesis. Furthermore, UBRs have widespread genetic alterations and expression perturbations in pan-cancer. The expression of 79 UBRs was identified to be closely correlated with the activity of 32 cancer hallmark-related pathways, and ten hub genes were screened for further clinical relevance analysis by a network-based method. More than 90% of UBRs can affect the survival of cancer patients, and hub genes have an excellent prognostic classification for specific cancer types. CONCLUSIONS: Our study provides a comprehensive analysis of UBRs in spermatogenesis and pan-cancer, which can build a foundation for understanding male infertility and developing cancer drugs in the aspect of ubiquitination.

BACKGROUND: Evident adolescent idiopathic scoliosis (AIS) incurs high treatment costs, low quality of life, and many complications. Early screening of AIS is essential to avoid progressing to an evident stage. However, there is no valid serum biomarker for AIS for early screening. METHODS: Antibody-based array is a large-scale study of proteins, which is expected to reveal a serum protein signature as biomarker for AIS. There are two segments of the research, including biomarkers screening and validation. In the biomarkers screening group, a total of 16 volunteers participated in this study, and we carried out differentially expressed proteins screening via protein array assay between No-AIS group and the AIS group, through which GeneSet enrichment analysis was performed. In the validation group with a total of 62 volunteers, the differentially expressed proteins from screening group were verified by Enzyme-Linked immunosorbent assay (ELISA), and then multiple regression analysis. RESULTS: In our study, there were twenty-nine differentially expressed proteins in AIS, through Protein array assay and GeneSet enrichment analysis in the biomarkers screening group. Then the expression of FAP, CD23 and B2M decreased as the degree of AIS increased via ELISA in validation group (FAP, p < 0.0001; CD23, p = 0.0002; B2M, p < 0.0001). Further, the results of multiple regression analysis showed that FAP, CD23 are linked to Cobb angle, whereas B2M were excluded because of multicollinearity. CONCLUSIONS: Altogether, we found that serum protein FAP and CD23 are intimately related to AIS, suggesting FAP and CD23 are expected to serve as the serum biomarkers, which significantly facilitate frequent longitudinal monitoring as to keep track of disease progression and tailor treatment accordingly.

BACKGROUND: Safety of probiotic strains that are used in human and animal trials is a prerequisite. Genome based safety assessment of probiotics has gained popularity due its cost efficiency and speed, and even became a part of national regulation on foods containing probiotics in Indonesia. However, reliability of the safety assessment based only on a full genome sequence is not clear. Here, for the first time, we sequenced, assembled, and analysed the genome of the probiotic strain Lactiplantibacillus plantarum IS-10506, that was isolated from dadih, a traditional fermented buffalo milk. The strain has already been used as a probiotic for more than a decade, and in several clinical trials proven to be completely safe. METHODS: The genome of the probiotic strain L. plantarum IS-10506 was sequenced using Nanopore sequencing technology, assembled, annotated and screened for potential harmful (PH) and beneficial genomic features. The presence of the PH features was assessed from general annotation, as well as with the use of specialised tools. In addition, PH regions in the genome were compared to all other probiotic and non-probiotic L. plantarum strains available in the NCBI RefSeq database. RESULTS: For the first time, a high-quality complete genome of L. plantarum IS-10506 was obtained, and an extensive search for PH and a beneficial signature was performed. We discovered a number of PH features within the genome of L. plantarum IS-10506 based on the general annotation, including various antibiotic resistant genes (AMR); however, with a few exceptions, bioinformatics tools specifically developed for AMR detection did not confirm their presence. We further demonstrated the presence of the detected PH genes across multiple L. plantarum strains, including probiotics, and overall high genetic similarities between strains. CONCLUSION: The genome of L. plantarum IS-10506 is predicted to have several PH features. However, the strain has been utilized as a probiotic for over a decade in several clinical trials without any adverse effects, even in immunocompromised children with HIV infection and undernourished children. This implies the presence of PH feature signatures within the probiotic genome does not necessarily indicate their manifestation during administration. Importantly, specialized tools for the search of PH features were found more robust and should be preferred over manual searches in a general annotation.

The autosomal dominant form of polycystic kidney disease (ADPKD) is the most common hereditary disease that causes late-onset renal cyst development and end-stage renal disease. Preimplantation genetic testing for monogenic disease (PGT-M) has emerged as an effective strategy to prevent pathogenic mutation transmission rely on SNP linkage analysis between pedigree members. Yet, it remains challenging to establish reliable PGT-M methods for ADPKD cases or other monogenic diseases with de novo mutations or without a family history. Here we reported the application of long-read sequencing for direct haplotyping in a female patient with de novo PKD1 c.11,526 G > C mutation and successfully established the high-risk haplotype. Together with targeted short-read sequencing of SNPs for the couple and embryos, the carrier status for embryos was identified. A healthy baby was born without the PKD1 pathogenic mutation. Our PGT-M strategy based on long-read sequencing for direct haplotyping combined with targeted SNP haplotype can be widely applied to other monogenic disease carriers with de novo mutation.

BACKGROUND: Previous preliminary work found that Latroeggtoxin-VI (LETX-VI), a proteinaceous neurotoxin from the eggs of spider Latrodectus tredecimguttatus, could promote the synthesis and release of dopamine in PC12 cells. However, the underlying mechanisms have not been fully clear. Here, the effects of LETX-VI on the gene expression profile and dopamine in PC12 cells were analyzed with the differential transcriptome-based strategies. RESULTS: After treatment of PC12 cells with LETX-VI for 24 h, a total of 356 differentially expressed transcripts were identified. Of them 165 were up-regulated and 191 down-regulated. Relevant GO analysis indicated that LETX-VI modulated the expression of certain genes and thereby affected multiple biological processes in PC12 cells, including protein metabolism, nucleic acid metabolism, substance transport, signaling, neurotransmitter metabolism and release. When western blot analysis was employed to confirm the abundance levels of synaptojanin 1 and synuclein alpha interacting protein, the representatives of highly up- and down-regulated transcript-encoded proteins that are closely related with dopamine respectively, it was found that the level of synaptojanin 1 in the PC12 cells treated with LETX-VI was increased, whereas that of synuclein alpha interacting protein was not obviously altered, suggesting that synaptojanin 1 may be much more involved in the effects of LETX-VI on dopamine. After synaptojanin 1 level was knocked down using siRNA, the levels of both total and released dopamine were significantly decreased, indicating that synaptojanin 1 is a protein positively modulating the synthesis and secretion of dopamine. When the PC12 cells with knocked down synaptojanin 1 were treated by LETX-VI, the adverse effects of synaptojanin 1 knockdown on dopamine were attenuated, confirming that LETX-VI promotes the synthesis and secretion of dopamine at least partially by enhancing the expression of the gene SYNJ1 encoding synaptojanin 1. CONCLUSIONS: This work demonstrates that LETX-VI exerts multiple regulatory effects on the cellular processes in PC12 cells by altering the gene expression profile. LETX-VI modulates the expression of the genes closely related to the synthesis, transport and release of neurotransmitters especially dopamine in PC12 cells, with the gene SYNJ1 encoding synaptojanin 1 as a main target.

BACKGROUND: Symbiotic nitrogen fixation differs among Bradyrhizobium japonicum strains. Soybean inoculated with USDA123 has a lower yield than strains known to have high nitrogen fixation efficiency, such as USDA110. In the main soybean-producing area in the Midwest of the United States, USDA123 has a high nodule incidence in field-grown soybean and is competitive but inefficient in nitrogen fixation. In this study, a high-throughput system was developed to characterize nodule number among 1,321 Glycine max and 69 Glycine soja accessions single inoculated with USDA110 and USDA123. RESULTS: Seventy-three G. max accessions with significantly different nodule number of USDA110 and USDA123 were identified. After double inoculating 35 of the 73 accessions, it was observed that PI189939, PI317335, PI324187B, PI548461, PI562373, and PI628961 were occupied by USDA110 and double-strain nodules but not by USDA123 nodules alone. PI567624 was only occupied by USDA110 nodules, and PI507429 restricted all strains. Analysis showed that 35 loci were associated with nodule number in G. max when inoculated with strain USDA110 and 35 loci with USDA123. Twenty-three loci were identified in G. soja when inoculated with strain USDA110 and 34 with USDA123. Only four loci were common across two treatments, and each locus could only explain 0.8 to 1.5% of phenotypic variation. CONCLUSIONS: High-throughput phenotyping systems to characterize nodule number and occupancy were developed, and soybean germplasm restricting rhizobium strain USDA123 but preferring USDA110 was identified. The larger number of minor effects and a small few common loci controlling the nodule number indicated trait genetic complexity and strain-dependent nodulation restriction. The information from the present study will add to the development of cultivars that limit USDA123, thereby increasing nitrogen fixation efficiency and productivity.

BACKGROUND: Anatabine, although being one of four major tobacco alkaloids, is never accumulated in high quantity in any of the naturally occurring species from the Nicotiana genus. Previous studies therefore focused on transgenic approaches to synthetize anatabine, most notably by generating transgenic lines with suppressed putrescine methyltransferase (PMT) activity. This led to promising results, but the global gene expression of plants with such distinct metabolism has not been analyzed. In the current study, we describe how these plants respond to topping and the downstream effects on alkaloid biosynthesis. RESULTS: The surge in anatabine accumulation in PMT transgenic lines after topping treatment and its effects on gene expression changes were analyzed. The results revealed increases in expression of isoflavone reductase-like (A622) and berberine bridge-like enzymes (BBLs) oxidoreductase genes, previously shown to be crucial for the final steps of nicotine biosynthesis. We also observed significantly higher methylputrescine oxidase (MPO) expression in all plants subjected to topping treatment. In order to investigate if MPO suppression would have the same effects as that of PMT, we generated transgenic plants. These plants with suppressed MPO expression showed an almost complete drop in leaf nicotine content, whereas leaf anatabine was observed to increase by a factor of ~ 1.6X. CONCLUSION: Our results are the first concrete evidence that suppression of MPO leads to decreased nicotine in favor of anatabine in tobacco roots and that this anatabine is successfully transported to tobacco leaves. Alkaloid transport in plants remains to be investigated to higher detail due to high variation of its efficiency among Nicotiana species and varieties of tobacco. Our research adds important step to better understand pyrrolidine ring biosynthesis and its effects on gene expression and subsequent accumulation of anatabine.

BACKGROUND: Increasing grain nutritional value in sorghum (Sorghum bicolor) is a paramount breeding objective, as is increasing drought resistance (DR), because sorghum is grown mainly in drought-prone areas. The genetic basis of grain nutritional traits remains largely unknown. Marker-assisted selection using significant loci identified through genome-wide association study (GWAS) shows potential for selecting desirable traits in crops. This study assessed natural variation available in sorghum accessions from around the globe to identify novel genes or genomic regions with potential for improving grain nutritional value, and to study associations between DR traits and grain weight and nutritional composition. RESULTS: We dissected the genetic architecture of grain nutritional composition, protein content, thousand-kernel weight (TKW), and plant height (PH) in sorghum through GWAS of 163 unique African and Asian accessions under irrigated and post-flowering drought conditions. Several QTLs were detected. Some were significantly associated with DR, TKW, PH, protein, and Zn, Mn, and Ca contents. Genomic regions on chromosomes 1, 2, 4, 8, 9, and 10 were associated with TKW, nutritional, and DR traits; colocalization patterns of these markers indicate potential for simultaneous improvement of these traits. In African accessions, markers associated with TKW were mapped to six regions also associated with protein, Zn, Ca, Mn, Na, and DR, suggesting the potential for simultaneous selection for higher grain nutrition and TKW. Our results indicate that it may be possible to select for increased DR on the basis of grain nutrition and weight potential. CONCLUSIONS: This study provides a valuable resource for selecting landraces for use in plant breeding programs and for identifying loci that may contribute to grain nutrition and weight with the hope of producing cultivars that combine improved yield traits, nutrition, and DR.

BACKGROUND: 5-methylcytosine (m5C) modification is widely associated with many biological and pathological processes. However, knowledge of m5C modification in osteoarthritis (OA) remains lacking. Thus, our study aimed to identify common m5C features in OA. RESULTS: In the present study, we identified 1395 differentially methylated genes (DMGs) and 1673 differentially expressed genes (DEGs) using methylated RNA immunoprecipitation next-generation sequencing (MeRIP-seq) and RNA-sequencing. A co-expression analysis of DMGs and DEGs showed that the expression of 133 genes was significantly affected by m5C methylation. A protein-protein interaction network of the 133 genes was constructed using the STRING database, and the cytoHubba plug-in of Cytoscape was used to hub genes were screen out 11 hub genes, including MMP14, VTN, COL15A1, COL6A2, SPARC, COL5A1, COL6A3, COL6A1, COL8A2, ADAMTS2 and COL7A1. The Pathway enrichment analysis by the ClueGO and CluePedia plugins in Cytoscape showed that the hub genes were significantly enriched in collagen degradation and extracellular matrix degradation. CONCLUSIONS: Our study indicated that m5C modification might play an important role in OA pathogenesis, and the present study provides worthwhile insight into identifying m5C-related therapeutic targets in OA.

BACKGROUND: Intramuscular fat (IMF) content is the major indicator for evaluating chicken meat quality due to its positive correlation with tenderness, juiciness, and flavor. An increasing number of studies are focusing on the functions of microRNAs (miRNAs) in intramuscular adipocyte differentiation. However, little is known about the association of miR-128-3p with intramuscular adipocyte differentiation. Our previous RNA-seq results indicated that miR-128-3p was differentially expressed at different periods in chicken intramuscular adipocytes, revealing a possible association with intramuscular adipogenesis. The purpose of this research was to investigate the biological functions and regulatory mechanism of miR-128-3p in chicken intramuscular adipogenesis. RESULTS: The results of a series of assays confirmed that miR-128-3p could promote the proliferation and inhibit the differentiation of intramuscular adipocytes. A total of 223 and 1,050 differentially expressed genes (DEGs) were identified in the mimic treatment group and inhibitor treatment group, respectively, compared with the control group. Functional enrichment analysis revealed that the DEGs were involved in lipid metabolism-related pathways, such as the MAPK and TGF-β signaling pathways. Furthermore, target gene prediction analysis showed that miR-128-3p can target many of the DEGs, such as FDPS, GGT5, TMEM37, and ASL2. The luciferase assay results showed that miR-128-3p targeted the 3' UTR of FDPS. The results of subsequent functional assays demonstrated that miR-128-3p acted as an inhibitor of intramuscular adipocyte differentiation by targeting FDPS. CONCLUSION: miR-128-3p inhibits chicken intramuscular adipocyte differentiation by downregulating FDPS. Our findings provide a theoretical basis for the study of lipid metabolism and reveal a potential target for molecular breeding to improve meat quality.

Rhodophyta are among the closest known relatives of green plants. Studying the codons of their genomes can help us understand the codon usage pattern and characteristics of the ancestor of green plants. By studying the codon usage pattern of all available red algae, it was found that although there are some differences among species, high-bias genes in most red algae prefer codons ending with GC. Correlation analysis, Nc-GC3s plots, parity rule 2 plots, neutrality plot analysis, differential protein region analysis and comparison of the nucleotide content of introns and flanking sequences showed that the bias phenomenon is likely to be influenced by local mutation pressure and natural selection, the latter of which is the dominant factor in terms of translation accuracy and efficiency. It is worth noting that selection on translation accuracy could even be detected in the low-bias genes of individual species. In addition, we identified 15 common optimal codons in seven red algae except for G. sulphuraria for the first time, most of which were found to be complementary and bound to the tRNA genes with the highest copy number. Interestingly, tRNA modification was found for the highly degenerate amino acids of all multicellular red algae and individual unicellular red algae, which indicates that highly biased genes tend to use modified tRNA in translation. Our research not only lays a foundation for exploring the characteristics of codon usage of the red algae as green plant ancestors, but will also facilitate the design and performance of transgenic work in some economic red algae in the future.

BACKGROUND: The mechanism of grain development in elite maize breeding lines has not been fully elucidated. Grain length, grain width and grain weight are key components of maize grain yield. Previously, using the Chinese elite maize breeding line Chang7-2 and its large grain mutant tc19, we characterized the grain size developmental difference between Chang7-2 and tc19 and performed transcriptomic analysis. RESULTS: In this paper, using Chang7-2 and tc19, we performed comparative transcriptomic, proteomic and metabolomic analyses at different grain development stages. Through proteomics analyses, we found 2884, 505 and 126 differentially expressed proteins (DEPs) at 14, 21 and 28 days after pollination, respectively. Through metabolomics analysis, we identified 51, 32 and 36 differentially accumulated metabolites (DAMs) at 14, 21 and 28 days after pollination, respectively. Through multiomics comparative analysis, we showed that the phenylpropanoid pathways are influenced at transcriptomic, proteomic and metabolomic levels in all the three grain developmental stages. CONCLUSION: We identified several genes in phenylpropanoid biosynthesis, which may be related to the large grain phenotype of tc19. In summary, our results provided new insights into maize grain development.

BACKGROUND: Biological denitrification has been commonly adopted for the removal of nitrogen from sewage effluents. However, due to the low temperature during winter, microorganisms in the wastewater biological treatment unit usually encounter problems such as slow cell growth and low enzymatic efficiency. Hence, the isolation and screening of cold-tolerant aerobic denitrifying bacteria (ADB) have recently drawn attention. In our previous study, two Pseudomonas strains PMCC200344 and PMCC200367 isolated from Arctic soil demonstrated strong denitrification ability at low temperatures. The two Arctic strains show potential for biological nitrogen removal from sewage in cold environments. However, the genome sequences of these two organisms have not been reported thus far. RESULTS: Here, the basic characteristics and genetic diversity of strains PMCC200344 and PMCC200367 were described, together with the complete genomes and comparative genomic results. The genome of Pseudomonas sp. PMCC200344 was composed of a circular chromosome of 6,478,166 bp with a G + C content of 58.60% and contained a total of 5,853 genes. The genome of Pseudomonas sp. PMCC200367 was composed of a circular chromosome of 6,360,061 bp with a G + C content of 58.68% and contained 5,801 genes. Not only prophages but also genomic islands were identified in the two Pseudomonas strains. No plasmids were observed. All genes of a complete set of denitrification pathways as well as various putative cold adaptation and heavy metal resistance genes in the genomes were identified and analyzed. These genes were usually detected on genomic islands in bacterial genomes. CONCLUSIONS: These analytical results provide insights into the genomic basis of microbial denitrification in cold environments, indicating the potential of Arctic Pseudomonas strains in nitrogen removal from sewage effluents at low temperatures.

BACKGROUND: Cytochrome P450 (Cytochrome P450s) genes are involved in the catalysis of various reactions, including growth, development, and secondary metabolite biosynthetic pathways. However, little is known about the characteristics and functions of the P450 gene family in Camellia sinensis (C. sinensis). RESULTS: To reveal the mechanisms of tea plant P450s coping with abiotic stresses, analyses of the tea plant P450 gene family were conducted using bioinformatics-based methods. In total, 273 putative P450 genes were identified from the genome database of C. sinensis. The results showed that P450s were well-balanced across the chromosomes I to XV of entire genome, with amino acid lengths of 268-612 aa, molecular weights of 30.95-68.5 kDa, and isoelectric points of 4.93-10.17. Phylogenetic analysis divided CsP450s into 34 subfamilies, of which CYP71 was the most abundant. The predicted subcellular localization results showed that P450 was distributed in a variety of organelles, with chloroplasts, plasma membrane,,and cytoplasm localized more frequently. The promoter region of CsP450s contained various cis-acting elements related to phytohormones and stress responses. In addition, ten conserved motifs (Motif1-Motif10) were identified in the CsP450 family proteins, with 27 genes lacking introns and only one exon. The results of genome large segment duplication showed that there were 37 pairs of genes with tandem duplication. Interaction network analysis showed that CsP450 could interact with multiple types of target genes, and there are protein interactions within the family. Tissue expression analysis showed that P450 was highly expressed in roots and stems. Moreover, qPCR analysis of the relative expression level of the gene under drought and cold stress correlated with the sequencing results. CONCLUSIONS: This study lays the foundation for resolving the classification and functional study of P450 family genes and provides a reference for the molecular breeding of C. sinensis.

BACKGROUND: Paenarthrobacter nicotinovorans ATCC 49919 uses the pyridine-pathway to degrade nicotine and could provide a renewable source of precursors from nicotine-containing waste as well as a model for studying the molecular evolution of catabolic pathways and their spread by horizontal gene transfer via soil bacterial plasmids. RESULTS: In the present study, the strain was sequenced using the Illumina NovaSeq 6000 and Oxford Nanopore Technology (ONT) MinION platforms. Following hybrid assembly with Unicycler, the complete genome sequence of the strain was obtained and used as reference for whole-genome-based phylogeny analyses. A total of 64 related genomes were analysed; five Arthrobacter strains showed both digital DNA-DNA hybridization and average nucleotide identity values over the species threshold when compared to P. nicotinovorans ATCC 49919. Five plasmids and two contigs belonging to Arthrobacter and Paenarthrobacter strains were shown to be virtually identical with the pAO1 plasmid of Paenarthrobacter nicotinovorans ATCC 49919. Moreover, a highly syntenic nic-genes cluster was identified on five plasmids, one contig and three chromosomes. The nic-genes cluster contains two major locally collinear blocks that appear to form a putative catabolic transposon. Although the origins of the nic-genes cluster and the putative transposon still elude us, we hypothesise here that the ATCC 49919 strain most probably evolved from Paenarthrobacter sp. YJN-D or a very closely related strain by acquiring the pAO1 megaplasmid and the nicotine degradation pathway. CONCLUSIONS: The data presented here offers another snapshot into the evolution of plasmids harboured by Arthrobacter and Paenarthrobacter species and their role in the spread of metabolic traits by horizontal gene transfer among related soil bacteria.

BACKGROUND: Atabecestat, a potent brain penetrable BACE1 inhibitor that reduces CSF amyloid beta (Aβ), was developed as an oral treatment for Alzheimer's disease (AD). Elevated liver enzyme adverse events were reported in three studies although only one case met Hy's law criteria to predict serious hepatotoxicity. METHOD: We performed a case-control genome-wide association study (GWAS) to identify genetic risk variants associated with liver enzyme elevation using 42 cases with alanine transaminase (ALT) above three times the upper limit of normal (ULN) and 141 controls below ULN. Additionally, we performed a GWAS using continuous maximal ALT/ULN (expressed as times the ULN) upon exposure to atabecestat as the outcome measure (n = 285). RESULTS: No variant passed the genome-wide significance threshold (p = 5 × 10- 8) in the case-control GWAS. We identified suggestive association signals in genes (NLRP1, SCIMP, and C1QBP) implicated in the inflammatory processes. Among the genes implicated by position mapping using variants suggestively associated (p < 1 × 10- 5) with ALT elevation case-control status, gene sets involved in innate immune response (adjusted p-value = 0.05) and regulation of cytokine production (adjusted p-value = 0.04) were enriched. One genomic region in the intronic region of GABRG3 passed the genome-wide significance threshold in the continuous max(ALT/ULN) GWAS, and this variant was nominally associated with ALT elevation case status (p = 0.009). CONCLUSION: The suggestive GWAS signals in the case-control GWAS analysis suggest the potential role of inflammation in atabecestat-induced liver enzyme elevation.

BACKGROUND: The cellular and molecular dynamics of human prepuce are crucial for understanding its biological and physiological functions, as well as the prevention of related genital diseases. However, the cellular compositions and heterogeneity of human prepuce at single-cell resolution are still largely unknown. Here we systematically dissected the prepuce of children and adults based on the single-cell RNA-seq data of 90,770 qualified cells. RESULTS: We identified 15 prepuce cell subtypes, including fibroblast, smooth muscle cells, T/natural killer cells, macrophages, vascular endothelial cells, and dendritic cells. The proportions of these cell types varied among different individuals as well as between children and adults. Moreover, we detected cell-type-specific gene regulatory networks (GRNs), which could contribute to the unique functions of related cell types. The GRNs were also highly dynamic between the prepuce cells of children and adults. Our cell-cell communication network analysis among different cell types revealed a set of child-specific (e.g., CD96, EPO, IFN-1, and WNT signaling pathways) and adult-specific (e.g., BMP10, NEGR, ncWNT, and NPR1 signaling pathways) signaling pathways. The variations of GRNs and cellular communications could be closely associated with prepuce development in children and prepuce maintenance in adults. CONCLUSIONS: Collectively, we systematically analyzed the cellular variations and molecular changes of the human prepuce at single-cell resolution. Our results gained insights into the heterogeneity of prepuce cells and shed light on the underlying molecular mechanisms of prepuce development and maintenance.

The nucleolus is a large nuclear body that serves as the primary site for ribosome biogenesis. Recent studies have suggested that it also plays an important role in organizing chromatin architecture. However, to establish a causal relationship between nucleolar ribosome assembly and chromatin architecture, genetic tools are required to disrupt nucleolar ribosome biogenesis. In this study, we used ATAC-seq to investigate changes in chromatin accessibility upon specific depletion of two ribosome biogenesis components, RPOA-2 and GRWD-1, in the model organism Caenorhabditis elegans. To facilitate the analysis of ATAC-seq data, we introduced two tools: SRAlign, an extensible NGS data processing workflow, and SRAtac, a customizable end-to-end ATAC-seq analysis pipeline. Our results revealed highly comparable changes in chromatin accessibility following both RPOA-2 and GRWD-1 perturbations. However, we observed a weak correlation between changes in chromatin accessibility and gene expression. While our findings corroborate the idea of a feedback mechanism between ribosomal RNA synthesis, nucleolar ribosome large subunit biogenesis, and chromatin structure during the L1 stage of C. elegans development, they also prompt questions regarding the functional impact of these alterations on gene expression.

BACKGROUND: As the prepubertal stage is a crucial point for the proper development of the mammary gland and milk production, this study aims to evaluate how protein restriction at this stage can affect methylation marks in milk somatic cells. Here, 28 Assaf ewes were subjected to 42.3% nutritional protein restriction (14 animals, NPR) or fed standard diets (14 animals, C) during the prepubertal stage. During the second lactation, the milk somatic cells of these ewes were sampled, and the extracted DNA was subjected to whole-genome bisulfite sequencing. RESULTS: A total of 1154 differentially methylated regions (DMRs) were identified between the NPR and C groups. Indeed, the results of functional enrichment analyses of the genes harboring these DMRs suggested their relevant effects on the development of the mammary gland and lipid metabolism in sheep. The additional analysis of the correlations of the mean methylation levels within these DMRs with fat, protein, and dry extract percentages in the milk and milk somatic cell counts suggested associations between several DMRs and milk production traits. However, there were no phenotypic differences in these traits between the NPR and C groups. CONCLUSION: In light of the above, the results obtained in the current study might suggest potential candidate genes for the regulation of milk production traits in the sheep mammary gland. Further studies focusing on elucidating the genetic mechanisms affected by the identified DMRs may help to better understand the biological mechanisms modified in the mammary gland of dairy sheep as a response to nutritional challenges and their potential effects on milk production.