Elevated levels of serum saturated fatty acid palmitate have been shown to promote insulin resistance, increase cellular ROS production, and trigger cell apoptosis in hepatocytes during the development of obesity. However, it remains unclear whether palmitate directly impacts the circadian clock in hepatocytes, which coordinates nutritional inputs and hormonal signaling with downstream metabolic outputs. Here we presented evidence that the molecular clock is a novel target of palmitate in hepatocytes. Palmitate exposure at low dose inhibits the molecular clock activity and suppresses the cyclic expression of circadian targets including Dbp, Nr1d1 and Per2 in hepatocytes. Palmitate treatment does not seem to alter localization or reduce protein expression of BMAL1 and CLOCK, the two core components of the molecular clock in hepatocytes. Instead, palmitate destabilizes the protein-protein interaction between BMAL1-CLOCK in a dose and time-dependent manner. Furthermore, we showed that SIRT1 activators could reverse the inhibitory action of palmitate on BMAL1-CLOCK interaction and the clock gene expression, whereas inhibitors of NAD synthesis mimic the palmitate effects on the clock function. In summary, our findings demonstrated that palmitate inhibits the clock function by suppressing SIRT1 function in hepatocytes.

Xyloglucan endo-transglycosylases (XETs) modify the xyloglucan-cellulose framework of plant cell walls and, thus, affect cell wall expansion and strength. Dissecting the mechanism by which natural variation in XETs affects wood properties can inform breeding efforts to improve wood quality and yield traits. To this end, we isolated a full-length PtoXET16A cDNA clone from Populus tomentosa. Real-time PCR analysis showed that PtoXET16A was maximally expressed in the root, followed by phloem, cambium, and developing xylem, suggesting that PtoXET16A plays important roles in the development of vascular tissues. Nucleotide diversity and linkage disequilibrium analysis revealed that PtoXET16A has high single nucleotide polymorphism (SNP) diversity (π = 0.01266 and θ(w) = 0.01392) and low linkage disequilibrium (r² ≥ 0.1, within 900 bp). SNP- and haplotype-based association analyses of 426 individuals from a natural population indicated that nine SNPs (including two non-synonymous markers and one splicing variant) (p ≤ 0.05, false discovery rate Q ≤ 0.01), and nine haplotypes (p ≤ 0.05) were significantly associated with growth and wood properties, each explaining from 3.40%-10.95% of phenotypic variance. This work shows that examination of allelic variation and linkage disequilibrium by a candidate-gene-based approach can help to decipher the genetic basis of wood formation. Moreover, the SNP markers identified in this study can potentially be applied for marker-assisted selection to improve growth and wood-property traits in Populus.

It is thought that methylcytosine can be inherited through meiosis and mitosis, and that epigenetic variation may be under genetic control or correlation may be caused by neutral drift. However, DNA methylation also varies with tissue, developmental stage, and environmental factors. Eliminating these factors, we analyzed the levels and patterns, diversity and structure of genomic methylcytosine in the xylem of nine natural populations of Chinese white poplar.

Although the molecular basis of poplar sex-specific flower development remains largely unknown, increasing evidence indicates an essential role for microRNAs (miRNAs). The specific miRNA types and precise miRNA expression patterns in dioecious plant flower development remain unclear. Here, we used andromonoecious poplar, an exceptional model system, to eliminate the confounding effects of genetic background of dioecious plants. This system, combined with high-throughput sequencing and computational analysis, allowed us to characterize sex-specific miRNAomes from female and male flowers. Comparative miRNAome analysis combined with quantitative real-time PCR revealed the expression patterns of 27 miRNAs in poplar flower and showed that the targets of these miRNAs are involved in flower organogenesis, Ca(2+) transport, phytohormone synthesis and metabolism, and DNA methylation. This paper describes a complex regulatory network consisting of these miRNAs expressed in sex-specific flower development in a dioecious plant. The conserved and novel miRNA locations were annotated in the Populus trichocarpa genome. Among these, miRNA Pto-F70 and 4 targets are located in the sex-determination regions of chromosome XIX. Furthermore, two novel miRNAs, Pto-F47 and Pto-F68, were shown for the first time to be regulatory factors in phytohormone interactions. To our knowledge, this report is the first systematic investigation of sex-specific flower-related miRNAs and their targets in poplar, and it deepens our understanding of the important regulatory functions of miRNAs in female and male flower development in this dioecious plant.

Transcription factors (TFs) regulate gene expression and can strongly affect phenotypes. However, few studies have examined TF variants and TF interactions with their targets in plants. Here, we used genetic association in 435 unrelated individuals of Populus tomentosa to explore the variants in Pto-Wuschela and its targets to decipher the genetic regulatory network of Pto-Wuschela. Our bioinformatics and co-expression analysis identified 53 genes with the motif TCACGTGA as putative targets of Pto-Wuschela. Single-marker association analysis showed that Pto-Wuschela was associated with wood properties, which is in agreement with the observation that it has higher expression in stem vascular tissues in Populus. Also, SNPs in the 53 targets were associated with growth or wood properties under additive or dominance effects, suggesting these genes and Pto-Wuschela may act in the same genetic pathways that affect variation in these quantitative traits. Epistasis analysis indicated that 75.5% of these genes directly or indirectly interacted Pto-Wuschela, revealing the coordinated genetic regulatory network formed by Pto-Wuschela and its targets. Thus, our study provides an alternative method for dissection of the interactions between a TF and its targets, which will strength our understanding of the regulatory roles of TFs in complex traits in plants.

In trees, xylem tissues play a key role in the formation of woody tissues, which have important uses for pulp and timber production; also DNA methylation plays an important part in gene regulation during xylogenesis in trees. In our study, methylation-sensitive amplified polymorphism (MSAP) analysis was used to analyze the role cytosine methylation plays in wood formation in the commercially important tree species Populus tomentosa. This analysis compared the methylation patterns between xylem tissues (developing xylem and mature xylem) and non-xylem tissues (cambium, shoot apex, young leaf, mature leaf, phloem, root, male catkin, and female catkin) and found 10,316 polymorphic methylation sites. MSAP identified 132 candidate genes with the same methylation patterns in xylem tissues, including seven wood-related genes. The expression of these genes differed significantly between xylem and non-xylem tissue types (P < 0.01). This indicated that the difference of expression of specific genes with unique methylation patterns, rather than relative methylation levels between the two tissue types plays a critical role in wood biosynthesis. However, 46.2% of candidate genes with the same methylation pattern in vascular tissues (cambium, phloem, and developing xylem) did not have distinct expression patterns in xylem and non-xylem tissue. Also, bisulfite sequencing and transcriptome sequencing of MYB, NAC and FASCICLIN-LIKE AGP 13 revealed that the location of cytosine methylation in the gene might affect the expression of different transcripts from the corresponding gene. The expression of different transcripts that produce distinct proteins from a single gene might play an important role in the regulation of xylogenesis.

Gene-derived simple sequence repeats (genic SSRs), also known as functional markers, are often preferred over random genomic markers because they represent variation in gene coding and/or regulatory regions. We characterized 544 genic SSR loci derived from 138 candidate genes involved in wood formation, distributed throughout the genome of Populus tomentosa, a key ecological and cultivated wood production species. Of these SSRs, three-quarters were located in the promoter or intron regions, and dinucleotide (59.7%) and trinucleotide repeat motifs (26.5%) predominated. By screening 15 wild P. tomentosa ecotypes, we identified 188 polymorphic genic SSRs with 861 alleles, 2-7 alleles for each marker. Transferability analysis of 30 random genic SSRs, testing whether these SSRs work in 26 genotypes of five genus Populus sections (outgroup, Salix matsudana), showed that 72% of the SSRs could be amplified in Turanga and 100% could be amplified in Leuce. Based on genotyping of these 26 genotypes, a neighbour-joining analysis showed the expected six phylogenetic groupings. In silico analysis of SSR variation in 220 sequences that are homologous between P. tomentosa and Populus trichocarpa suggested that genic SSR variations between relatives were predominantly affected by repeat motif variations or flanking sequence mutations. Inheritance tests and single-marker associations demonstrated the power of genic SSRs in family-based linkage mapping and candidate gene-based association studies, as well as marker-assisted selection and comparative genomic studies of P. tomentosa and related species.

The tetramer is the most important form for acetylcholinesterase in physiological conditions, i.e., in the neuromuscular junction and the nervous system. It is important to study the diffusion of acetylcholine to the active sites of the tetrameric enzyme to understand the overall signal transduction process in these cellular components. Crystallographic studies revealed two different forms of tetramers, suggesting a flexible tetramer model for acetylcholinesterase. Using a recently developed finite element solver for the steady-state Smoluchowski equation, we have calculated the reaction rate for three mouse acetylcholinesterase tetramers using these two crystal structures and an intermediate structure as templates. Our results show that the reaction rates differ for different individual active sites in the compact tetramer crystal structure, and the rates are similar for different individual active sites in the other crystal structure and the intermediate structure. In the limit of zero salt, the reaction rates per active site for the tetramers are the same as that for the monomer, whereas at higher ionic strength, the rates per active site for the tetramers are approximately 67%-75% of the rate for the monomer. By analyzing the effect of electrostatic forces on ACh diffusion, we find that electrostatic forces play an even more important role for the tetramers than for the monomer. This study also shows that the finite element solver is well suited for solving the diffusion problem within complicated geometries.

This article describes the numerical solution of the time-dependent Smoluchowski equation to study diffusion in biomolecular systems. Specifically, finite element methods have been developed to calculate ligand binding rate constants for large biomolecules. The resulting software has been validated and applied to the mouse acetylcholinesterase (mAChE) monomer and several tetramers. Rates for inhibitor binding to mAChE were calculated at various ionic strengths with several different time steps. Calculated rates show very good agreement with experimental and theoretical steady-state studies. Furthermore, these finite element methods require significantly fewer computational resources than existing particle-based Brownian dynamics methods and are robust for complicated geometries. The key finding of biological importance is that the rate accelerations of the monomeric and tetrameric mAChE that result from electrostatic steering are preserved under the non-steady-state conditions that are expected to occur in physiological circumstances.

Complex RNA transcription and processing produces a diverse range catalog of long noncoding RNAs (lncRNAs), important biological regulators that have been implicated in osmotic stress responses in plants. Promoter upstream transcript (PROMPT) lncRNAs share some regulatory elements with the promoters of their neighbouring protein-coding genes. However, their function remains unknown. Here, using strand-specific RNA sequencing, we identified 209 differentially regulated osmotic-responsive PROMPTs in poplar (Populus simonii). PROMPTs are transcribed bidirectionally and are more stable than other lncRNAs. Co-expression analysis of PROMPTs and protein-coding genes divided the regulatory network into five independent subnetworks including 27 network modules. Significantly enriched PROMPTs in the network were selected to validate their regulatory roles. We used delaminated layered double hydroxide lactate nanosheets (LDH-lactate-NS) to transport synthetic nucleic acids into live tissues to mimic overexpression and interference of a specific PROMPT. The altered expression of PROMPT_1281 induced the expression of its cis and trans targets, and this interaction was governed by its secondary structure rather than just its primary sequence. Based on this example, we proposed a model that a concentration gradient of PROMPT_1281 is established, which increases the probability of its interaction with targets near its transcription site that shares common motifs. Our results firstly demonstrated that PROMPT_1281 act as carriers of MYB transcription factors to induce the expression of target genes under osmotic stress. In sum, our study identified and validated a set of poplar PROMPTs that likely have regulatory functions in osmotic responses.

Long non-coding RNAs (lncRNAs) are important regulatory factors for plant growth and development, but little is known about the allelic interactions of lncRNAs with mRNA in perennial plants. Here, we analyzed the interaction of the NERD (Needed for RDR2-independent DNA methylation) Populus tomentosa gene PtoNERD with its putative regulator, the lncRNA NERDL (NERD-related lncRNA), which partially overlaps with the promoter region of this gene. Expression analysis in eight tissues showed a positive correlation between NERDL and PtoNERD (r = 0.62), suggesting that the interaction of NERDL with its putative target might be involved in wood formation. We conducted association mapping in a natural population of P. tomentosa (435 unrelated individuals) to evaluate genetic variation and the interaction of the lncRNA NERDL with PtoNERD. Using additive and dominant models, we identified 30 SNPs (P < 0.01) associated with five tree growth and wood property traits. Each SNP explained 3.90-8.57% of phenotypic variance, suggesting that NERDL and its putative target play a common role in wood formation. Epistasis analysis uncovered nine SNP-SNP association pairs between NERDL and PtoNERD, with an information gain of -7.55 to 2.16%, reflecting the strong interactions between NERDL and its putative target. This analysis provides a powerful method for deciphering the genetic interactions of lncRNAs with mRNA and dissecting the complex genetic network of quantitative traits in trees.

Impairment of insulin clearance is being increasingly recognized as a critical step in the development of insulin resistance and metabolic disease. The carcinoembryonic antigen-related cell adhesion molecule 1 (CEACAM1) promotes insulin clearance. Null deletion or liver-specific inactivation of Ceacam1 in mice causes a defect in insulin clearance, insulin resistance, steatohepatitis, and visceral obesity. Immunohistological analysis revealed reduction of hepatic CEACAM1 in obese subjects with fatty liver disease. Thus, we aimed to determine whether this occurs at the hepatocyte level in response to systemic extrahepatic factors and whether this holds across species. Northern and Western blot analyses demonstrate that CEACAM1 mRNA and protein levels are reduced in liver tissues of obese individuals compared to their lean age-matched counterparts. Furthermore, Western analysis reveals a comparable reduction of CEACAM1 protein in primary hepatocytes derived from the same obese subjects. Similar to humans, Ceacam1 mRNA level, assessed by quantitative RT-PCR analysis, is significantly reduced in the livers of obese Zucker (fa/fa, ZDF) and Koletsky (f/f) rats relative to their age-matched lean counterparts. These studies demonstrate that the reduction of hepatic CEACAM1 in obesity occurs at the level of hepatocytes and identify the reduction of hepatic CEACAM1 as a common denominator of obesity across multiple species.

As an essential component in reducing anthropogenic CO2 emissions to the atmosphere, tree planting is the key to keeping carbon dioxide emissions under control. In 1992, the United Nations agreed to take action at the Earth Summit to stabilize and reduce net zero global anthropogenic CO2 emissions. Tree planting was identified as an effective method to offset CO2 emissions. A high net photosynthetic rate (Pn) with fast-growing trees could efficiently fulfill the goal of CO2 emission reduction. Net photosynthetic rate model can provide refernece for plant's stability of photosynthesis productivity.

New kinase has emerged throughout evolution, but how new kinase evolve while maintaining their functions and acquiring new functions remains unclear. Fructokinase (FRK), the gateway kinase to fructose metabolism, plays essential roles in plant development, and stress tolerance. Here, we explored the evolution of FRK gene family in 20 plant species (from green algae to angiosperms) and their functional roles in Populus. We identified 125 putative FRK genes in the 20 plant species with an average of 6 members per species. Phylogenetic analysis separated these 125 genes into 8 clades including 3 conserved clades and 5 specific clades, the 5 of which only exist in green algae or angiosperms. Evolutionary analysis revealed that FRK genes in ancient land plants have the largest number of functional domains with the longest amino acid sequences, and the length of FRK genes became shorter during the transition to vascular plants. This was accompanied by loss, acquisition, and diversification of functional domains. In Populus, segmental duplication appears to be the main mechanism for the expansion of FRK genes. Specially, most FRK genes duplicated in salicoids are regulated by Populus-specific microRNAs. Furthermore, compared with common FRKs, Populus-specific FRKs have showed higher expression specificity and are associated with fewer growth and wood property traits, which suggests that these FRKs may have undergone functional divergence. Our study explores the specific roles of FRKs in the Populus genome and provides new insights for functional investigation of this gene family.

Idiopathic pulmonary fibrosis (IPF) is a fatal and age-related pulmonary disease. Nintedanib is a receptor tyrosine kinase inhibitor, and one of the only two listed drugs against IPF. Regorafenib is a novel, orally active, multi-kinase inhibitor that has similar targets to nintedanib and is applied to treat colorectal cancer and gastrointestinal stromal tumors in patients. In this study, we first identified that regorafenib could alleviate bleomycin-induced pulmonary fibrosis in mice. The in vivo experiments indicated that regorafenib suppresses collagen accumulation and myofibroblast activation. Further in vitro mechanism studies showed that regorafenib inhibits the activation and migration of myofibroblasts and extracellular matrix production, mainly through suppressing the transforming growth factor (TGF)-β1/Smad and non-Smad signaling pathways. In vitro studies have also indicated that regorafenib could augment autophagy in myofibroblasts by suppressing TGF-β1/mTOR (mechanistic target of rapamycin) signaling, and could promote apoptosis in myofibroblasts. In conclusion, regorafenib attenuates bleomycin-induced pulmonary fibrosis by suppressing the TGF-β1 signaling pathway.

This study investigated the different blood pressure patterns that were evaluated by ambulatory blood pressure monitoring (ABPM) among elderly patients and explored the effect of pressure patterns on cognitive impairment and mortality.

Complete and high-resolution temperature observation data are important input parameters for agrometeorological disaster monitoring and ecosystem modelling. Due to the limitation of field meteorological observation conditions, observation data are commonly missing, and an appropriate data imputation method is necessary in meteorological data applications. In this paper, we focus on filling long gaps in meteorological observation data at field sites. A deep learning-based model, BiLSTM-I, is proposed to impute missing half-hourly temperature observations with high accuracy by considering temperature observations obtained manually at a low frequency. An encoder-decoder structure is adopted by BiLSTM-I, which is conducive to fully learning the potential distribution pattern of data. In addition, the BiLSTM-I model error function incorporates the difference between the final estimates and true observations. Therefore, the error function evaluates the imputation results more directly, and the model convergence error and the imputation accuracy are directly related, thus ensuring that the imputation error can be minimized at the time the model converges. The experimental analysis results show that the BiLSTM-I model designed in this paper is superior to other methods. For a test set with a time interval gap of 30 days, or a time interval gap of 60 days, the root mean square errors (RMSEs) remain stable, indicating the model's excellent generalization ability for different missing value gaps. Although the model is only applied to temperature data imputation in this study, it also has the potential to be applied to other meteorological dataset-filling scenarios.

Compared with other forest systems, research interest in the potential for a stronger ecosystem carbon sequestration of evergreen forests throughout subtropical China has greatly increased. The eddy covariance technique is widely employed to determine accurate forest-atmosphere carbon dioxide (CO2) flux, which is subsequently used to determine forest ecosystem carbon exchange characteristics. The Dinghushan Biosphere Reserve, a subtropical monsoon evergreen broad-leaved forest, is a suitable study area due to its warm and humid climate (compared with other regions within the same latitude), consequently playing a role in the carbon cycle in the region. For this study, we hypothesized that the forest land in this region generally acts as a carbon sink, and that its carbon sequestration capacity increases over time despite the influence of climatic factors. Here, we compared net CO2 flux data derived from the eddy covariance technique over an 8-year study window. Additionally, we ascertained the effects of various environmental factors on net CO2 flux, while also using the Michaelis-Menten model and a physiologically based process model to track and report on ecosystem carbon exchange characteristics. We observed seasonal trends in daily ecosystem flux, indicative of sensitivity to climatic factors, such as air temperature, precipitation, and sunlight. The carbon sequestration capacity of the region exhibited seasonal variability, increasing from October to March (-264 g C m-2 year-1, i.e., 48.4%) while weakening from April to September (-150 g C m-2 year-1, i.e., 40.4%) on average. The net ecosystem exchange (NEE) rate varied from -518 to -211 g C m-2 year-1; ecosystem respiration (Re) varied from 1,142 to 899 g C m-2 year-1; and gross primary production (GPP) varied from 1,552 to 1,254 g C m-2 year-1. This study found that even though the Dinghushan Biosphere Reserve generally acts as a carbon sink, its carbon sequestration capacity did not increase significantly throughout the study period. The techniques (models) used in this study are suitable for application in other ecosystems globally, which can aid in their management and conservation. Finally, the Dinghushan Biosphere Reserve is both an exemplary and a model forest system useful in exploring CO2 absorption and sequestration from the atmosphere.

Economically important traits in many species generally show polygenic, quantitative inheritance. The components of genetic variation (additive, dominant and epistatic effects) of these traits conferred by multiple genes in shared biological pathways remain to be defined. Here, we investigated 11 full-length genes in cellulose biosynthesis, on 10 growth and wood-property traits, within a population of 460 unrelated Populus tomentosa individuals, via multi-gene association. To validate positive associations, we conducted single-marker analysis in a linkage population of 1,200 individuals. We identified 118, 121, and 43 associations (P< 0.01) corresponding to additive, dominant, and epistatic effects, respectively, with low to moderate proportions of phenotypic variance (R(2)). Epistatic interaction models uncovered a combination of three non-synonymous sites from three unique genes, representing a significant epistasis for diameter at breast height and stem volume. Single-marker analysis validated 61 associations (false discovery rate, Q ≤ 0.10), representing 38 SNPs from nine genes, and its average effect (R(2) = 3.8%) nearly 2-fold higher than that identified with multi-gene association, suggesting that multi-gene association can capture smaller individual variants. Moreover, a structural gene-gene network based on tissue-specific transcript abundances provides a better understanding of the multi-gene pathway affecting tree growth and lignocellulose biosynthesis. Our study highlights the importance of pathway-based multiple gene associations to uncover the nature of genetic variance for quantitative traits and may drive novel progress in molecular breeding.

Variation in regulatory factors, including microRNAs (miRNAs), contributes to variation in quantitative and complex traits. However, in plants, variants in miRNAs and their target genes that contribute to natural phenotypic variation, and the underlying regulatory networks, remain poorly characterized. We investigated the associations and interactions of single-nucleotide polymorphisms (SNPs) in miRNAs and their target genes with phenotypes in 435 individuals from a natural population of Populus. We used RNA-seq to identify 217 miRNAs differentially expressed in a tension wood system, and identified 1196 candidate target genes; degradome sequencing confirmed 60 of the target sites. In addition, 72 miRNA-target pairs showed significant co-expression. Gene ontology (GO) term analysis showed that most of the genes in the co-regulated pairs participate in biological regulation. Genome resequencing found 5383 common SNPs (frequency ≥ 0.05) in 139 miRNAs and 31 037 SNPs in 819 target genes. Single-SNP association analyses identified 232 significant associations between wood traits (P ≤ 0.05) and SNPs in 102 miRNAs and 1387 associations with 478 target genes. Among these, 102 miRNA-target pairs associated with the same traits. Multi-SNP associations found 102 epistatic pairs associated with traits. Furthermore, a reconstructed regulatory network contained 12 significantly co-expressed pairs, including eight miRNAs and nine targets associated with traits. Lastly, both expression and genetic association showed that miR156i, miR156j, miR396a and miR6445b were involved in the formation of tension wood. This study shows that variants in miRNAs and target genes contribute to natural phenotypic variation and annotated roles and interactions of miRNAs and their target genes by genetic association analysis.