Water-Accelerated Photooxidation of CH3NH3PbI3 Perovskite.

Halide perovskites have the potential to disrupt the photovoltaics market based on their high performance and low cost. However, the decomposition of perovskites under moisture, oxygen, and light raises concerns about service lifetime, especially because degradation mechanisms and the corresponding rate laws that fit the observed data have thus far eluded researchers. Here, we report a water-accelerated photooxidation mechanism dominating the degradation kinetics of archetypal perovskite CH3NH3PbI3 in air under >1% relative humidity at 25 °C. From this mechanism, we develop a kinetic model that quantitatively predicts the degradation rate as a function of temperature, ambient O2 and H2O levels, and illumination. Because water is a possible product of dry photooxidation, these results highlight the need for encapsulation schemes that rigorously block oxygen ingress, as product water may accumulate beneath the encapsulant and initiate the more rapid water-accelerated photooxidative decomposition.

Mutation landscape of TSC1/TSC2 in Chinese patients with tuberous sclerosis complex.

Genetic testing of TSC1 and TSC2 is important for the diagnosis of tuberous sclerosis complex (TSC), an autosomal dominant neurocutaneous disease. This study retrospectively reviewed 347 samples from patients with clinically suspected TSC being tested for mutations in TSC1 and TSC2 genes using next-generation sequencing and multiplex ligation-dependent probe amplification. Two hundred eighty-one patients (80.98%) were classified as definite/possible/uncertain diagnosis of TSC and the mutational spectrum of TSC1/TSC2 was described. Two hundred eighteen unique nonsynonymous SNVs/Indels (64 in TSC1, 154 in TSC2) and 13 copy number variants (CNVs) were identified in 241 samples (85.77%), including 82 novel variants. CNVs involving 12 large deletions and one duplication were detected exclusively in TSC2. Both TSC1 and TSC2 mutations were nearly uniformly distributed in their protein-coding regions. Furthermore, a string of non-TSC1/TSC2 deleterious variants in 12 genes was identified in the patients, especially overwhelmingly present in the patients with no mutation identified (NMI) in TSC1/TSC2. Our study provides a comprehensive TSC1/TSC2 mutation landscape and reveal some potential risk non-TSCs variants present in patients with NMI.

A prognostic 11 long noncoding RNA expression signature for breast invasive carcinoma.

Breast cancer, the most common cancer in women worldwide, is associated with high mortality. The long non-coding RNAs (lncRNAs) with a little capacity of coding proteins is playing an increasingly important role in the cancer paradigm. Accumulating evidences demonstrate that lncRNAs have crucial connections with breast cancer prognosis while the studies of lncRNAs in breast cancer are still in its primary stage. In this study, we collected 1052 clinical patient samples, a comparatively large sample size, including 13 159 lncRNA expression profiles of breast invasive carcinoma (BRCA) from The Cancer Genome Atlas database to identify prognosis-related lncRNAs. We randomly separated all of these clinical patient samples into training and testing sets. In the training set, we performed univariable Cox regression analysis for primary screening and played the model for Robust likelihood-based survival for 1000 times. Then 11 lncRNAs with a frequency more than 600 were selected for prediction of the prognosis of BRCA. Using the analysis of multivariate Cox regression, we established a signature risk-score formula for 11 lncRNA to identify the relationship between lncRNA signatures and overall survival. The 11 lncRNA signature was validated both in the testing and the complete set and could effectively classify the high-/low-risk group with different OS. We also verified our results in different stages. Moreover, we analyzed the connection between the 11 lncRNAs and the genes of ESR1, PGR, and Her2, of which protein products (ESR, PGR, and HER2) were used to classify the breast cancer subtypes widely. The results indicated correlations between 11 lncRNAs and the gene of PGR and ESR1. Thus, a prognostic model for 11 lncRNA expression was developed to classify the BRAC clinical patient samples, providing new avenues in understanding the potential therapeutic methods of breast cancer.

Integrated analysis of transcription factors and targets co-expression profiles reveals reduced correlation between transcription factors and target genes in cancer.

Transcription factors are recognized as the key regulators of gene expression. However, the changes in the correlation of transcription factors and their target genes between normal and tumor tissues are usually ignored. In this research, we used mRNA expression profile data from The Cancer Genome Atlas which included 5726 samples across 11 major human cancers to perform co-expression analysis by the Pearson correlation coefficients. Then, integrating 81,357 pairs of transcription factors and target genes from transcription factors databases to find out the changes in the co-expression correlation of these gene pairs from normal to tumor tissues. Based on the changes in the number of co-expressed TF-TG pairs and changes in the level of co-expression, we found the generally reduced correlation between transcription factors and their target genes in cancer. Additionally, we screened out universal and specific transcription factors-target genes pairs which may significant influence particular cancer. Then, we obtained 423 cancer cell line expression profiles from Broad Institute Cancer Cell Line Encyclopedia to verify our results. Some of these pairs like XRCC5-XRCC6 have been reported to involve in multiple cancers, while pairs like IRF1-PSMB9 without any previous articles related to tumor but involve in the biological processes of cancer, which are of great potential to be therapeutic targets. Our research may provide insights to better understand the tumor development mechanisms and find potential therapeutic targets.

Controlled Synthesis and Photocatalytic Performance of Au@ZnO Nanospheres with Core-Shell and Yolk-Shell Structures Assisted by Carbonaceous Layers as Intermediate.

Au@ZnO nanospheres with controlled core-shell and yolk-shell features were successfully fabricated through a facile template-based solvothermal method using carbonaceous layers as intermediate. The active carbonaceous layers containing numerous surface functional groups such as OH and C═O play key roles in fabrication of the Au@ZnO nanospheres. The structures of Au@ZnO can easily controlled by adjusting the concentration of Zn ions aqueous solution and/or the size of Au@C microspheres. The synthesis strategy proposed in this paper is also verified to be suitable for fabricating some other Au@metallic oxides nanospheres such as Au@SnO2 and Au@Fe2O3 with controlled structures. Furthermore, due to the synergic effect of the surface plasma resonance effect of Au cores and the large surface area of 122.67 m2 g-1, the Au@ZnO nanospheres with yolk-shell feature exhibite good photocatalytic activity and stability under visible-light irradiation.

Ancestry, Plasmodium cynomolgi prevalence and rhesus macaque admixture in cynomolgus macaques (Macaca fascicularis) bred for export in Chinese breeding farms.

BACKGROUND: Most cynomolgus macaques (Macaca fascicularis) used in the United States as animal models are imported from Chinese breeding farms without documented ancestry. Cynomolgus macaques with varying rhesus macaque ancestry proportions may exhibit differences, such as susceptibility to malaria, that affect their suitability as a research model. METHODS: DNA of 400 cynomolgus macaques from 10 Chinese breeding farms was genotyped to characterize their regional origin and rhesus ancestry proportion. A nested PCR assay was used to detect Plasmodium cynomolgi infection in sampled individuals. RESULTS: All populations exhibited high levels of genetic heterogeneity and low levels of inbreeding and genetic subdivision. Almost all individuals exhibited an Indochinese origin and a rhesus ancestry proportion of 5%-48%. The incidence of P. cynomolgi infection in cynomolgus macaques is strongly associated with proportion of rhesus ancestry. CONCLUSIONS: The varying amount of rhesus ancestry in cynomolgus macaques underscores the importance of monitoring their genetic similarity in malaria research.

Multi-omics data integration reveals the molecular network of dysregulation IQGAP2-mTOR promotes cell proliferation.

IQGAP2 as a tumor suppressor gene can influence cell proliferation in multiple tumor cell lines. However, the regulation network of cell proliferation resulting solely from the deficiency of IQGAP2 in cells was still unclear. Here, we integrated transcriptome, proteome, and phosphoproteome analyses to investigate the regulatory network of cell proliferation in IQGAP2 knockdown HaCaT and HEK293 cells. Our findings revealed that the dysregulation of the IQGAP2-mTOR molecular network led to increased cell proliferation. We demonstrated that IQGAP2 knockdown enhanced the phosphorylation levels of AKT and S6K, leading to increased cell proliferation. Additionally, we found that AKT and mTOR inhibitors partially rescued abnormal cell proliferation by reducing hyperphosphorylation. Our data suggest a potential connection between the mTOR signaling pathway and aberrant cell proliferation in IQGAP2 knockdown cells. These findings offer a new therapeutic strategy for patients with IQGAP2 deficiency.

RNA atlas and competing endogenous RNA regulation in tissue-derived exosomes from luminal B and triple-negative breast cancer patients.

Background: Luminal B and triple-negative breast cancer (TNBC) are malignant subtypes of breast cancer (BC), which can be attributed to the multifaceted roles of tissue-derived exosomes (T-exos). Competing endogenous RNA (ceRNA) networks can regulate gene expression post-transcriptionally. Methods: RNAs in T-exos from luminal B BC (n=8) and TNBC (n=8) patients were compared with those from persons with benign breast disease (n=8). The differentially expressed (DE) mRNA, microRNA (miRNA), and long noncoding RNA (lncRNA) target genes were annotated using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) to reveal the relevant biological processes.The ceRNA networks were constructed to show distinct regulation, and the mRNAs involved were annotated. The miRNAs involved in the ceRNA networks were screened with the Kaplan-Meier Plotter database to identify dysregulated ceRNAs with prognostic power. Results: In total, 802 DE mRNAs, 441 DE lncRNAs, and 104 DE miRNAs were identified in luminal B BC T-exos, while 1699 DE mRNAs, 590 DE lncRNAs, and 277 DE miRNAs were identified in TNBC T-exos. Gene annotation revealed that the RAS-MAPK pathway was the primary biological process in luminal B BC T-exos, while endocrine system development and growth were the main processes in TNBC T-exos. Survival analysis established seven survival-related lncRNA/miRNA/mRNA regulations in luminal B BC T-exos, and nineteen survival-related lncRNA/miRNA/mRNA regulations in TNBC T-exos. Conclusion: In addition to survival-related ceRNA regulations, ceRNA regulation of RAS-MAPK in luminal B and endocrine system development and growth regulation in TNBC might contribute to the tumorigenesis of BC.

MYD88-Mutated Chronic Lymphocytic Leukaemia/Small Lymphocytic Lymphoma as a Distinctive Molecular Subgroup Is Associated with Atypical Immunophenotypes in Chinese Patients.

Chronic lymphocytic leukaemia/small lymphocytic lymphoma (CLL/SLL) is a heterogeneous disease in Western and Chinese populations, and it is still not well characterized in Chinese patients. Based on a large cohort of newly diagnosed CLL/SLL patients from China, we investigated immunophenotypes, genetic abnormalities, and their correlations. Eighty-four percent of the CLL/SLL patients showed typical immunophenotypes with scores of 4 or 5 points in the Royal Marsden Hospital (RMH) scoring system (classic group), and the remaining 16% of patients were atypical with scores lower than 4 points (atypical group). Trisomy 12 and variants of TP53, NOTCH1, SF3B1, ATM, and MYD88 were the most recurrent genetic aberrations. Additionally, unsupervised genomic analysis based on molecular genetics revealed distinctive characteristics of MYD88 variants in CLL/SLL. By overlapping different correlation grouping analysis from genetics to immunophenotypes, the results showed MYD88 variants to be highly related to atypical CLL/SLL immunophenotypes. Furthermore, compared with mantle cell lymphoma (MCL), the genetic landscape showed potential value in clinical differential diagnosis of atypical CLL/SLL and MCL patients. These results reveal immunophenotypic and genetic features, and may provide insights into the tumorigenesis and clinical management of Chinese CLL/SLL patients.

Machine learning models-based on integration of next-generation sequencing testing and tumor cell sizes improve subtype classification of mature B-cell neoplasms.

Background: Next-generation sequencing (NGS) panels for mature B-cell neoplasms (MBNs) are widely applied clinically but have yet to be routinely used in a manner that is suitable for subtype differential diagnosis. This study retrospectively investigated newly diagnosed cases of MBNs from our laboratory to investigate mutation landscapes in Chinese patients with MBNs and to combine mutational information and machine learning (ML) into clinical applications for MBNs, especially for subtype classification. Methods: Samples from the Catalogue Of Somatic Mutations In Cancer (COSMIC) database were collected for ML model construction and cases from our laboratory were used for ML model validation. Five repeats of 10-fold cross-validation Random Forest algorithm was used for ML model construction. Mutation detection was performed by NGS and tumor cell size was confirmed by cell morphology and/or flow cytometry in our laboratory. Results: Totally 849 newly diagnosed MBN cases from our laboratory were retrospectively identified and included in mutational landscape analyses. Patterns of gene mutations in a variety of MBN subtypes were found, important to investigate tumorigenesis in MBNs. A long list of novel mutations was revealed, valuable to both functional studies and clinical applications. By combining gene mutation information revealed by NGS and ML, we established ML models that provide valuable information for MBN subtype classification. In total, 8895 cases of 8 subtypes of MBNs in the COSMIC database were collected and utilized for ML model construction, and the models were validated on the 849 MBN cases from our laboratory. A series of ML models was constructed in this study, and the most efficient model, with an accuracy of 0.87, was based on integration of NGS testing and tumor cell sizes. Conclusions: The ML models were of great significance in the differential diagnosis of all cases and different MBN subtypes. Additionally, using NGS results to assist in subtype classification of MBNs by method of ML has positive clinical potential.

Iodine Promoted Ultralow Zn Nucleation Overpotential and Zn-Rich Cathode for Low-Cost, Fast-Production and High-Energy Density Anode-Free Zn-Iodine Batteries.

The anode-free design is a promising strategy to increase the energy density of aqueous Zn metal batteries (AZMBs). However, the scarcity of Zn-rich cathodes and the rapid loss of limited Zn greatly hinder their commercial applications. To address these issues, a novel anode-free Zn-iodine battery (AFZIB) was designed via a simple, low-cost and scalable approach. Iodine plays bifunctional roles in improving the AFZIB overall performance: enabling high-performance Zn-rich cathode and modulating Zn deposition behavior. On the cathode side, the ZnI2 serves as Zn-rich cathode material. The graphene/polyvinyl pyrrolidone heterostructure was employed as an efficient host for ZnI2 to enhance electron conductivity and suppress the shuttle effect of iodine species. On the anode side, trace I3- additive in the electrolyte creates surface reconstruction on the commercial Cu foil. The in situ formed zincophilic Cu nanocluster allows ultralow-overpotential and uniform Zn deposition and superior reversibility (average coulombic efficiency > 99.91% over 7,000 cycles). Based on such a configuration, AFZIB exhibits significantly increased energy density (162 Wh kg-1) and durable cycle stability (63.8% capacity retention after 200 cycles) under practical application conditions. Considering the low cost and simple preparation methods of the electrode materials, this work paves the way for the practical application of AZMBs.

Analyses of Long-Term Epidemic Trends and Evolution Characteristics of Haplotype Subtypes Reveal the Dynamic Selection on SARS-CoV-2.

The scale of SARS-CoV-2 infection and death is so enormous that further study of the molecular and evolutionary characteristics of SARS-CoV-2 will help us better understand and respond to SARS-CoV-2 outbreaks. The present study analyzed the epidemic and evolutionary characteristics of haplotype subtypes or regions based on 1.8 million high-quality SARS-CoV-2 genomic data. The estimated ratio of the rates of non-synonymous to synonymous changes (Ka/Ks) in North America and the United States were always more than 1.0, while the Ka/Ks in other continents and countries showed a sharp decline, then a slow increase to 1.0, and a dramatic increase over time. H1 (B.1) with the highest substitution rate has become the most dominant haplotype subtype since March 2020 and has evolved into multiple haplotype subtypes with smaller substitution rates. Many evolutionary characteristics of early SARS-CoV-2, such as H3 being the only early haplotype subtype that existed for the shortest time, the global prevalence of H1 and H1-5 (B.1.1) within a month after being detected, and many high divergent genome sequences early in February 2020, indicate the missing of early SARS-CoV-2 genomic data. SARS-CoV-2 experienced dynamic selection from December 2019 to August 2021 and has been under strong positive selection since May 2021. Its transmissibility and the ability of immune escape may be greatly enhanced over time. This will bring greater challenges to the control of the pandemic.

Dilution effect for highly efficient multiple-component organic solar cells.

Although the multiple-component (MC) blend strategy has been frequently used as a very effective way to improve the performance of organic solar cells (OSCs), there is a strong need to understand the fundamental working mechanism and material selection rule for achieving optimal MC-OSCs. Here we present the 'dilution effect' as the mechanism for MC-OSCs, where two highly miscible components are molecularly intermixed. Contrary to the aggregation-induced non-radiative decay, the dilution effect enables higher luminescence quantum efficiencies and open-circuit voltages (VOC) in MC-OSCs via suppressed electron-vibration coupling. The continuously broadened bandgap together with reduced electron-vibration coupling also explains the composition-dependent VOC in ternary blends well. Moreover, we show that electrons can transfer between different acceptors, depending on the energy offset between them, which contributes to the largely unperturbed charge transport and high fill factors in MC-OSCs. The discovery of the dilution effect enables the demonstration of a high power conversion efficiency of 18.31% in an MC-OSC.

AutoVEM: An automated tool to real-time monitor epidemic trends and key mutations in SARS-CoV-2 evolution.

With the global epidemic of SARS-CoV-2, it is important to effectively monitor the variation, haplotype subgroup epidemic trends and key mutations of SARS-CoV-2 over time. This is of great significance to the development of new vaccines, the update of therapeutic drugs, and the improvement of detection methods. The AutoVEM tool developed in the present study could complete all mutations detections, haplotypes classification, haplotype subgroup epidemic trends and candidate key mutations analysis for 131,576 SARS-CoV-2 genome sequences in 18 h on a 1 core CPU and 2 GB RAM computer. Through haplotype subgroup epidemic trends analysis of 131,576 genome sequences, the great significance of the previous 4 specific sites (C241T, C3037T, C14408T and A23403G) was further revealed, and 6 new mutation sites of highly linked (T445C, C6286T, C22227T, G25563T, C26801G and G29645T) were discovered for the first time that might be related to the infectivity, pathogenicity or host adaptability of SARS-CoV-2. In brief, we proposed an integrative method and developed an efficient automated tool to monitor haplotype subgroup epidemic trends and screen for the candidate key mutations in the evolution of SARS-CoV-2 over time for the first time, and all data could be updated quickly to track the prevalence of previous key mutations and new candidate key mutations because of high efficiency of the tool. In addition, the idea of combinatorial analysis in the present study can also provide a reference for the mutation monitoring of other viruses.

Comprehensive evolution and molecular characteristics of a large number of SARS-CoV-2 genomes reveal its epidemic trends.

OBJECTIVES: To further reveal the phylogenetic evolution and molecular characteristics of the whole genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) based on a large number of genomes and provide a basis for the prevention and treatment of SARS-CoV-2. METHODS: Various evolution analysis methods were employed. RESULTS: The estimated ratio of the rates of non-synonymous to synonymous changes (Ka/Ks) of SARS-CoV-2 was 1.008 or 1.094 based on 622 or 3624 SARS-CoV-2 genomes and nine key specific sites of high linkage, and four major haplotypes were found: H1, H2, H3 and H4. The results of Ka/Ks, detected population size and development trends of each major haplotype showed that H3 and H4 subgroups were going through a purify evolution and almost disappeared after detection, indicating that they might have existed for a long time. The H1 and H2 subgroups were going through a near neutral or neutral evolution and globally increased with time, and the frequency of H1 was generally high in Europe and correlated with the death rate (r >0.37), suggesting that these two haplotypes might relate to the infectivity or pathogenicity of SARS-CoV-2. CONCLUSIONS: Several key specific sites and haplotypes related to the infectivity or pathogenicity of SARS-CoV-2, and the possible earlier origin time and place of SARS-CoV-2 were indicated based on the evolution and epidemiology of 16,373 SARS-CoV-2 genomes.

Interface Modified Flexible Printed Conductive Films via Ag2O Nanoparticle Decorated Ag Flake Inks.

A new approach to stable, low resistance inexpensive printed flexible conductive inks is proposed. Silver inks have been extensively studied and commercialized for applications in printed electronics due to the inherent high conductivity and stability of silver, even in particulate-based percolation networks processed at temperatures compatible with low-cost polymer films such as poly(ethylene terephthalate) (PET). Recent interest in flexible and even stretchable circuits, however, has presented new challenges for particle-based inks as mechanical strains can result in the opening of critical particle-to-particle contacts. Here we report a facile, low-cost method for the single-step synthesis of stable, printable nanoscale Ag2O decorated Ag flake inks which can be converted to highly conductive Ag films at 150 °C curing temperature without the use of limited shelf life organometallics or low metal loading nanoparticles to modify the interface between silver flakes. Analysis indicates that decoration of Ag flakes with Ag2O nanoparticles (NPs) during ink synthesis improves the conductivity and flexibility of printed silver films by forming bridging interconnections between Ag flakes after low temperature reduction of the Ag2O NPs. In this work, printed nanodecorated silver conductors with starting oxide to metal weight ratios of 5:95 exhibited lateral resistivities lower than 1.5 × 10-5 Ω·cm, which was 35% less than films derived from undecorated Ag flake inks of the same total Ag loading and binder system. This resistivity difference increased to 45% after cyclic bend testing showing increased resilience to repeated flexing for the nanodecorated inks. Through detailed compositional and morphological characterizations, we demonstrate that such improved conductivity and flexibility are due to a more effective bridging afforded by the in situ synthesized Ag NPs on the surface of Ag flakes. These properties, combined with the simplified syntheses method of the nanoink, make the material a viable, advantageous alternative to the limited number of stretchable conductors currently available.

Systematic comparison of germline variant calling pipelines cross multiple next-generation sequencers.

The development and innovation of next generation sequencing (NGS) and the subsequent analysis tools have gain popularity in scientific researches and clinical diagnostic applications. Hence, a systematic comparison of the sequencing platforms and variant calling pipelines could provide significant guidance to NGS-based scientific and clinical genomics. In this study, we compared the performance, concordance and operating efficiency of 27 combinations of sequencing platforms and variant calling pipelines, testing three variant calling pipelines-Genome Analysis Tool Kit HaplotypeCaller, Strelka2 and Samtools-Varscan2 for nine data sets for the NA12878 genome sequenced by different platforms including BGISEQ500, MGISEQ2000, HiSeq4000, NovaSeq and HiSeq Xten. For the variants calling performance of 12 combinations in WES datasets, all combinations displayed good performance in calling SNPs, with their F-scores entirely higher than 0.96, and their performance in calling INDELs varies from 0.75 to 0.91. And all 15 combinations in WGS datasets also manifested good performance, with F-scores in calling SNPs were entirely higher than 0.975 and their performance in calling INDELs varies from 0.71 to 0.93. All of these combinations manifested high concordance in variant identification, while the divergence of variants identification in WGS datasets were larger than that in WES datasets. We also down-sampled the original WES and WGS datasets at a series of gradient coverage across multiple platforms, then the variants calling period consumed by the three pipelines at each coverage were counted, respectively. For the GIAB datasets on both BGI and Illumina platforms, Strelka2 manifested its ultra-performance in detecting accuracy and processing efficiency compared with other two pipelines on each sequencing platform, which was recommended in the further promotion and application of next generation sequencing technology. The results of our researches will provide useful and comprehensive guidelines for personal or organizational researchers in reliable and consistent variants identification.

Transcriptomic responses of hypothalamus to acute exercise in type 2 diabetic Goto-Kakizaki rats.

The hypothalamus has an integral role in energy homeostasis regulation, and its dysfunctions lead to the development of type 2 diabetes (T2D). Physical activity positively affects the prevention and treatment of T2D. However, there is not much information on the adaptive mechanisms of the hypothalamus. In this study, RNA sequencing was used to determine how acute exercise affects hypothalamic transcriptome from both type 2 diabetic Goto-Kakizaki (GK) and control Wistar rats with or without a single session of running (15 m/min for 60 min). Through pairwise comparisons, we identified 957 differentially expressed genes (DEGs), of which 726, 197, and 98 genes were found between GK and Wistar, exercised GK and GK, and exercised Wistar and Wistar, respectively. The results of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment revealed that lipid metabolism-related terms and pathways were enriched in GK and exercised GK rats, and nervous system related terms and pathways were enriched in exercised GK and Wistar rats. Furthermore, 45 DEGs were associated with T2D and related phenotypes according to the annotations in the Rat Genome Database. Among these 45 DEGs, several genes (Plin2, Cd36, Lpl, Wfs1, Cck) related to lipid metabolism or the nervous system are associated with the exercise-induced benefits in the hypothalamus of GK rats. Our findings might assist in identifying potential therapeutic targets for T2D prevention and treatment.

Transcriptome Changes of Skeletal Muscle RNA-Seq Speculates the Mechanism of Postprandial Hyperglycemia in Diabetic Goto-Kakizaki Rats During the Early Stage of T2D.

To address how skeletal muscle contributes to postprandial hyperglycemia, we performed skeletal muscle transcriptome analysis of diabetic Goto-Kakizaki (GK) and control Wistar rats by RNA sequencing (RNA-Seq). We obtained 600 and 1785 differentially expressed genes in GK rats compared to those Wistar rats at three and four weeks of age, respectively. Specifically, Tbc1d4, involved in glucose uptake, was significantly downregulated in the skeletal muscle of GK aged both three and four weeks compared to those of age-matched Wistar rats. Pdk4, related to glucose uptake and oxidation, was significantly upregulated in the skeletal muscle of GK aged both three and four weeks compared to that of age-matched Wistar rats. Genes (Acadl, Acsl1 and Fabp4) implicated in fatty acid oxidation were significantly upregulated in the skeletal muscle of GK aged four weeks compared to those of age-matched Wistar rats. The overexpression or knockout of Tbc1d4, Pdk4, Acadl, Acsl1 and Fabp4 has been reported to change glucose uptake and fatty acid oxidation directly in rodents. By taking the results of previous studies into consideration, we speculated that dysregulation of key dysregulated genes (Tbc1d4, Pdk4, Acadl, Acsl1 and Fabp4) may lead to a decrease in glucose uptake and oxidation, and an increase in fatty acid oxidation in GK skeletal muscle at three and four weeks, which may, in turn, contribute to postprandial hyperglycemia. Our research revealed transcriptome changes in GK skeletal muscle at three and four weeks. Tbc1d4, Acadl, Acsl1 and Fabp4 were found to be associated with early diabetes in GK rats for the first time, which may provide a new scope for pathogenesis of postprandial hyperglycemia.