UDP-glucuronosyltransferases (UGTs), being multifunctional detoxification enzymes, play a major role in the process of resistance to various pesticides in insects. However, the mechanism underlying the molecular regulation of pesticide resistance remains unclear, especially in Apis cerana cerana. In this study, all of the UGTs in Apis cerana cerana (AccUGT) have been identified through the multiple alignment and phylogenetic analysis. Expression of AccUGT genes under different pesticides, and antioxidant genes after silencing of AccUGT2B20-like, were detected by qRT-PCR. The resistance of overexpressed AccUGT2B20-like to oxidative stress was investigated by an Escherichia coli overexpression system. Also, antioxidant-related enzyme activity was detected after silencing of the AccUGT2B20-like gene. Expression pattern analysis showed that almost all UGT genes were upregulated under different pesticide treatments. This result indicated that AccUGTs participate in the detoxification process of pesticides. AccUGT2B20-like was the major gene because it was more highly induced than the others. Overexpression of AccUGT2B20-like in E. coli could effectively improve oxidative stress resistance. Specifically, silencing the AccUGT2B20-like gene increased oxidative stress by repressing the expression of oxidation-related genes, decreasing antioxidant-related enzyme activity, and increasing malondialdehyde concentration. Taken together, our results indicate that AccUGTs are involved in pesticide resistance, among which, AccUGT2B20-like contributes to the detoxification of pesticides by eliminating oxidative stress in Apis cerana cerana. This study explains the molecular basis for the resistance of bees to pesticides and provides an important safeguard for maintaining ecological balance.

There are dozens of recognized indigenous dog breeds in China. However, these breeds have not had extensive studies to describe their population structure, genomic linkage disequilibrium (LD) patterns, and selection signatures. Here, we systematically surveyed the genomes of 157 unrelated dogs that were from 15 diverse Chinese dog breeds. Canine 170K SNP chips were used to compare the genomic structures of Chinese and Western dogs. The genotyping data of 170K SNP chips in Western dogs were downloaded from the LUPA (a European initiative of canine genome project) database. Chinese indigenous dogs had lower LD and shorter accumulative runs of homozygosity (ROH) in the genome. The genetic distances between individuals within each Chinese breed were larger than those within Western breeds. Chinese indigenous and Western dog breeds were clearly differentiated into two separate clades revealed by the PCA and NJ-tree. We found evidence for historical introgression of Western dogs into Chinese Kazakhstan shepherd and Mongolia Xi dogs. We suggested that Greenland sledge dog, Papillon, and European Eurasier have Chinese dog lineages. Selection sweep analysis identified genome-wide selection signatures of each Chinese breed and three breed groups. We highlighted several genes including EPAS1 and DNAH9 that show signatures of natural selection in Qinghai-Tibetan plateau dogs and are likely important for genetic adaptation to high altitude. Comparison of our findings with previous reports suggested RBP7, NMNAT1, SLC2A5, and H6PD that exhibit signatures of natural selection in Chinese mountain hounds as promising candidate genes for the traits of endurance and night vision, and NOL8, KRT9, RORB, and CAMTA1 that show signals of selection in Xi dogs might be candidate genes influencing dog running speed. The results about genomic and population structures, and selection signatures of Chinese dog breeds reinforce the conclusion that Chinese indigenous dogs with great variations of phenotypes are important resources for identifying genes responsible for complex traits.

The characteristics of head and neck squamous cell carcinoma (HNSCC) across different anatomic sites in the Chinese population have not been studied. To determine the genomic abnormalities underlying HNSCC across different anatomic sites, the alterations of selected cancer-related genes were evaluated.

Hepatocellular carcinoma (HCC) is ranked fifth among the most common cancer worldwide. Hypoxia can induce tumor growth, but the relationship with HCC prognosis remains unclear. Our study aims to construct a hypoxia-related multigene model to predict the prognosis of HCC.

Forsythiae Fructus (Lianqiao in Chinese) is widely used in traditional Chinese medicine. The lipid components in Forsythiae Fructus are the basis of plant growth and active metabolism. Samples were collected at two growth stages for a comprehensive study. Transcriptome and lipidomics were performed by using the RNA-seq and UPLC-Q-TOF-MS techniques separately. For the first time, it was reported that there were 5802 lipid components in Lianqiao comprised of 31.7% glycerolipids, 16.57% phospholipids, 13.18% sphingolipids, and 10.54% fatty acids. Lipid components such as terpenes and flavonoids have pharmacological activity, but their content was low. Among these lipids which were isolated from Forsythiae Fructus, 139 showed significant differences from the May and July harvest periods. The lipids of natural products are mainly concentrated in pregnenolones and polyvinyl lipids. RNA-Seq analysis revealed 92,294 unigenes, and 1533 of these were differentially expressed. There were 551 differential genes enriched in 119 KEGG pathways. The de novo synthesis pathways of terpenoids and flavonoids were explored. Combined with the results of lipidomics and transcriptomics, it is hypothesized that in the synthesis of abscisic acid, a terpenoid, may be under the dynamic regulation of genes EC: 1.1.1.288, EC: 1.14.14.137 and EC: 1.13.11.51 in balanced state. In the synthesis of gibberellin, GA20-oxidase (GA20ox, EC: 1.14.11.12), and GA3-oxidase (GA3ox, EC: 1.14.11.15) catalyze the production of active GAs, and EC: 1.14.11.13 is the metabolic enzymes of active GAs. In the synthesis of flavonoids, MF (multifunctional), PAL (phenylalanine ammonia-lyase), CHS (chalcone synthase), ANS (anthocyanidin synthase), FLS (flavonol synthase) are all key enzymes. The results of the present study provide valuable reference information for further research on the metabolic pathways of the secondary metabolites of Forsythia suspensa.

Bacteroides ovatus ELH-B2 is considered as a potential next-generation probiotic due to its preventive effects on lipopolysaccharides-associated inflammation and intestinal microbiota disorders in mice. To study safety issues associated with B. ovatus ELH-B2, we conducted comprehensive and systematic experiments, including in vitro genetic assessments of potential virulence and antimicrobial resistance genes, and an in vivo acute toxicity study of both immunocompetent and immunosuppressed mice via cyclophosphamide treatment. The results indicated that this novel strain is non-toxigenic, fragilysin is not expressed, and most of potential virulence genes are correlated with cellular structures such as capsular polysaccharide and polysaccharide utilizations. The antibiotic resistance features are unlikely be transferred to other intestinal microorganisms as no plasmids nor related genomic islands were identified. Side effects were not observed in mice. B. ovatus ELH-B2 also alleviated the damages caused by cyclophosphamide injection.

Background: Ephrin A3 (EFNA3), like most genes in the ephrin family, plays a central role in embryonic development and can be dysregulated in a variety of tumors. However, the relationship between EFNA3 and gastric cancer (GC) prognosis and tumor-infiltrating lymphocytes remains unclear. Methods: Tumor Immune Estimation Resource (TIMER) and Gene Expression Profiling Interactive Analysis 2 (GEPIA2) were used to analyze the expression of EFNA3. Kaplan-Meier plots and GEPIA2 were used to evaluate the relationship between EFNA3 expression and GC prognosis. Univariable survival and multivariate Cox analyses were used to compare various clinical characteristics with survival. LinkedOmics database was used for gene set enrichment analysis (GSEA). TIMER database and CIBERSORT algorithm were used to examine the relationship between EFNA3 expression and immune infiltration in GC and to explore cumulative survival in GC. The relationship between EFNA3 and immune checkpoints was examined using cBioPortal genomics analysis. Finally, EFNA3 expression in GC cells and tissues was assayed using quantitative real-time polymerase chain reaction. Results: EFNA3 expression differs in a variety of cancers, and EFNA3 expression was higher in GC tissue than normal gastric tissue. GC patients with high expression of EFNA3 had worse overall survival, disease-free survival, and first progression. Multivariate analysis identified EFNA3 as an independent prognostic factor for GC. GSEA identified ribosome, cell cycle, ribosome biogenesis in eukaryotes, and aminoacyl-tRNA biosynthesis pathways as differentially enriched in patients with high EFNA3 expression. B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells were significantly negatively correlated with a variety of immune markers. EFNA3 participates in changes in GC immune checkpoint markers in a collinear manner. EFNA3 expression in HGC-27, AGS, MKN45, and NCI-N87 was cell lines higher than that in GES-1, and patients with high expression of EFNA3 had a worse prognosis. Conclusion: EFNA3 can be used as a prognostic and immune infiltration and checkpoint marker in GC patients.

Autosomal dominant mental retardation-7 (MRD7) is a rare anomaly, characterized by severe intellectual disability, feeding difficulties, behavior abnormalities, and distinctive facial features, including microcephaly, deep-set eyes, large simple ears, and a pointed or bulbous nasal tip. Some studies show that the disorder has a close correlation with variants in DYRK1A. Herein we described a Chinese girl presenting typical clinical features diagnosed at 4 years old. Whole-exome sequencing of the familial genomic DNA identified a novel mutation c.930C > A (p.Tyr310*) in exon 7 of DYRK1A in the proband. The nonsense mutation was predicted to render the truncation of the protein. Our results suggested that the de novo heterozygous mutation in DYRK1A was responsible for the MRD7 in this Chinese family, which both extended the knowledge of mutation spectrum in MRD7 patients and highlighted the clinical application of exome sequencing.

Background: Chemotherapy resistance based on fluorouracil and cisplatin is one of the most encountered postoperative clinical problems in patients diagnosed with gastric cancer (GC), resulting in poor prognosis. Aim of the Study: This study aimed to combine autophagy-related genes (ARGs) to investigate the susceptibility patients with GC to postoperative chemotherapy. Methods: Based on The Cancer Genome Atlas (TCGA) database, gene expression data for GC patients undergoing chemotherapy were integrated and analyzed. Prognostic genes were screened based on univariate and multivariate analysis regression analysis. Subjects were divided into high-risk and low-risk groups according to the median risk score. Kaplan-Meier method was used to evaluate OS and DFS. The accuracy of the prediction was determined by the subject operating characteristic curve analysis. In addition, stratified analyses based on different clinical variables was performed to assess the correlation between risk scores and clinical variables. Quantitative real-time (qRT) PCR was used to verify the expression of CXCR4 in GC tissues and cell lines. Results: A total of nine ARGs related to the prognosis of chemotherapy patients were screened out. Compared with normal gastric mucosa cell, CXCR4 showed elevated expression in GC and was significantly associated with survival. Based on GEO and TCGA databases, the model accurately predicted DFS and OS after chemotherapy. Conclusion: This study established prognostic markers based on nine genes, predicting that ARGs are related to chemotherapy susceptibility of GC patients, which can provide better individualized treatment regimens for clinical practice.

Clinical data mining and bioinformatics analysis can be employed effectively to elucidate the function and underlying mechanisms of the gene of interest. Here, we have proposed a framework for the identification and validation of independent biomarkers in human cancer and for mechanistic profiling using gene sets enrichment analysis and pathway analysis. This is followed by validation with in vitro experiments. Using this framework to analyze the clinical relevance of SEC23A, we have discovered the prognostic potential of SEC23A in different cancers and identified SEC23A as an independent prognostic factor for poor prognosis in bladder cancer, which implicates SEC23A, for the first time, as an oncogene. Bioinformatic analyses have elucidated an association between SEC23A expression and the upregulation of the MAPK signaling pathway. Using the T24 human bladder cell line, we confirmed that knockdown of SEC23A expression could effectively impact the MAPK signaling pathway. Further, through PCR verification, we showed that MEF2A, one of the key genes of the MAPK signaling pathway, might be a downstream factor of the SEC23A gene.

Aegilops tauschii is the diploid progenitor of the D subgenome of hexaploid wheat (Triticum aestivum L.). Here, the phenotypic data of kernel length (KL), kernel width (KW), kernel volume (KV), kernel surface area (KSA), kernel width to length ratio (KWL), and hundred-kernel weight (HKW) for 223 A. tauschii accessions were gathered across three continuous years. Based on population structure analysis, 223 A. tauschii were divided into two subpopulations, namely T-group (mainly included A. tauschii ssp. tauschii accessions) and S-group (mainly included A. tauschii ssp. strangulata). Classifications based on cluster analysis were highly consistent with the population structure results. Meanwhile, the extent of linkage disequilibrium decay distance (r 2 = 0.5) was about 110 kb and 290 kb for T-group and S-group, respectively. Furthermore, a genome-wide association analysis was performed on these kernel traits using 6,723 single nucleotide polymorphism (SNP) markers. Sixty-six significant markers, distributed on all seven chromosomes, were identified using a mixed linear model explaining 4.82-13.36% of the phenotypic variations. Among them, 15, 28, 22, 14, 21, and 13 SNPs were identified for KL, KW, KV, KSA, KWL, and HKW, respectively. Moreover, six candidate genes that may control kernel traits were identified (AET2Gv20774800, AET4Gv20799000, AET5Gv20005900, AET5Gv20084100, AET7Gv20644900, and AET5Gv21111700). The transfer of beneficial genes from A. tauschii to wheat using marker-assisted selection will broaden the wheat D subgenome improve the efficiency of breeding.

Besides being one of the most prevalent cancers among women, incidence and mortality rates of endometrial cancer (EC) are still increasing. The E2F family of transcriptional factors is involved in cell differentiation, apoptosis, and inhibition of DNA damage response, thus affecting growth and invasion of tumor cells.

Spinal muscular atrophy (SMA) is a severe motor neuron degenerative disease caused by loss-of-function mutations in the survival motor neuron gene SMN1. It is widely posited that defective gene expression underlies SMA. However, the identities of these affected genes remain to be elucidated. By analyzing the transcriptome of a Caenorhabditis elegans SMA model at the pre-symptomatic stage, we found that the expression of numerous nuclear encoded mitochondrial genes and vacuolar H+-ATPase genes was significantly down-regulated, while that of histone genes was significantly up-regulated. We previously showed that the uaf-1 gene, encoding key splicing factor U2AF large subunit, could affect the behavior and lifespan of smn-1 mutants. Here, we found that smn-1 and uaf-1 interact to affect the recognition of 3' and 5' splice sites in a gene-specific manner. Altogether, our results suggest a functional interaction between smn-1 and uaf-1 in affecting RNA splicing and a potential effect of smn-1 on the expression of mitochondrial and histone genes.

The Duroc × (Landrace × Large White) hybrid pig (DLY) is the most popular commercial pig used in the Chinese pig industry. DLY pigs are usually white but sometimes show colored phenotypes. Colored DLY pigs are not favored by slaughterhouses and retailers, thus causing certain economic losses to farmers in China. In this study, we first conducted a genome-wide association study and RNA sequencing to demonstrate that KIT variants are responsible for diversifying coat color phenotypes segregating in a DLY population. We then defined the precise sizes and locations of four duplications (DUP1-4), four candidate causative mutations at the KIT locus, in the pig reference genome using the whole-genome sequence data of representative colored individuals. The sequence data also enabled us to identify a list of new KIT alleles. By investigating the association between these new alleles and coat color phenotypes, we provide further evidence that DUP2 is another causative mutation for the solid white coat color in pigs. DUP1 (the KIT gene duplication), DUP2 and the splice mutation are all required for the manifestation of a solid white coat color. DUP4 had a more significant effect on the formation of the belt phenotype compared with DUP3. Given the necessity of DUP2 for the solid white coat color, we detected IN /IN homozygotes lacking DUP2 in Large White and Landrace pigs and found that French Landrace pigs had the highest frequency (8.98%) of IN /IN individuals. This study not only advances our understanding of the molecular mechanism of the color phenotype in pigs, but also establishes a simple and accurate method for the screening of KIT IN /IN homozygotes in Large White and Landrace that would cause colored DLY pigs.

DNA methylation plays an essential role in the pathogenesis of coronary artery disease (CAD) through regulating mRNA expressions. This study aimed to identify hub genes regulated by DNA methylation as biomarkers of CAD. Gene expression and methylation datasets of peripheral blood leukocytes (PBLs) of CAD were downloaded from the Gene Expression Omnibus (GEO) database. Subsequently, multiple computational approaches were performed to analyze the regulatory networks and to recognize hub genes. Finally, top hub genes were verified in a case-control study, based on their differential expressions and methylation levels between CAD cases and controls. In total, 535 differentially expressed-methylated genes (DEMGs) were identified and partitioned into 4 subgroups. TSS200 and 5'UTR were confirmed as high enrichment areas of differentially methylated CpGs sites (DMCs). The function of DEMGs is enriched in processes of histone H3-K27 methylation, regulation of post-transcription and DNA-directed RNA polymerase activity. Pathway enrichment showed DEMGs participated in the VEGF signaling pathway, adipocytokine signaling pathway, and PI3K-Akt signaling pathway. Besides, expressions of hub genes fibronectin 1 (FN1), phosphatase (PTEN), and tensin homolog and RNA polymerase III subunit A (POLR3A) were discordantly expressed between CAD patients and controls and related with DNA methylation levels. In conclusion, our study identified the potential biomarkers of PBLs for CAD, in which FN1, PTEN, and POLR3A were confirmed.

Gastric cancer (GC) is a product of multiple genetic abnormalities, including genetic and epigenetic modifications. This study aimed to integrate various biomolecules, such as miRNAs, mRNA, and DNA methylation, into a genome-wide network and develop a nomogram for predicting the overall survival (OS) of GC.

Homologous recombination (HR), the most significant event in meiosis, has important implications for genetic diversity and evolution in organisms. Heteroduplex DNA (hDNA), the product of HR, can be captured by artificially induced chromosome doubling during the development of the embryo sac to inhibit postmeiotic segregation, subsequently, and hDNAs are directly detected using codominant simple sequence repeat (SSR) markers. In the present study, two hybrid triploid populations derived from doubling the chromosomes of the embryo sac induced by high temperature in Populus tomentosa served as starting materials. Eighty-seven, 62, and 79 SSR markers on chromosomes 01, 04, and 19, respectively, that were heterozygous in the maternal parent and different from the paternal parent were screened to detect and characterize the hDNA in P. tomentosa. The results showed that the hDNA frequency patterns on chromosomes changed slightly when the number of SSR primers increased. The highest hDNA frequency occurred at the adjacent terminal on chromosomes, which was slightly higher than those at the terminals in the two genotypic individuals, and the hDNA frequency gradually decreased as the locus-centromere distance decreased. With the increase in the number of SSR markers employed for detection, the number of recombination events (REs) detected significantly increased. In regions with high methylation or long terminal repeat (LTR) retrotransposon enrichment, the frequency of hDNA was low, and high frequencies were observed in regions with low sequence complexity and high gene density. High-frequency recombination occurring at high gene density regions strongly affected the association between molecular markers and quantitative trait loci (QTLs), which was an important factor contributing to the difficulty encountered by MAS in achieving the expected breeding results.

Cannabinoid receptor 1 activation by the major psychoactive component in cannabis, Δ9-tetrahydrocannabinol (THC), produces motor impairments, hypothermia, and analgesia upon acute exposure. In previous work, we demonstrated significant sex and strain differences in acute responses to THC following administration of a single dose (10 mg/kg, i.p.) in C57BL/6J (B6) and DBA/2J (D2) inbred mice. To determine the extent to which these differences are heritable, we quantified acute responses to a single dose of THC (10 mg/kg, i.p.) in males and females from 20 members of the BXD family of inbred strains derived by crossing and inbreeding B6 and D2 mice. Acute THC responses (initial sensitivity) were quantified as changes from baseline for: 1. spontaneous activity in the open field (mobility), 2. body temperature (hypothermia), and 3. tail withdrawal latency to a thermal stimulus (antinociception). Initial sensitivity to the immobilizing, hypothermic, and antinociceptive effects of THC varied substantially across the BXD family. Heritability was highest for mobility and hypothermia traits, indicating that segregating genetic variants modulate initial sensitivity to THC. We identified genomic loci and candidate genes, including Ndufs2, Scp2, Rps6kb1 or P70S6K, Pde4d, and Pten, that may control variation in THC initial sensitivity. We also detected strong correlations between initial responses to THC and legacy phenotypes related to intake or response to other drugs of abuse (cocaine, ethanol, and morphine). Our study demonstrates the feasibility of mapping genes and variants modulating THC responses in the BXDs to systematically define biological processes and liabilities associated with drug use and abuse.

CircRNAs have been reported to play essential roles in regulating immunity and inflammation, which may be an important regulatory factor in the development of vitiligo. However, the expression profile of circRNAs and their potential biological functions in vitiligo have not been reported so far. In our study we found there are 64 dysregulated circRNAs and 14 dysregulated miRNAs in the patients with vitiligo. Through the correlation analysis, we obtained 12 dysregulated circRNAs and 5 dysregulated miRNAs, forming 48 relationships in the circRNA-miRNA-mRNA regulatory network. Gene Ontology analysis indicated dysregulated circRNAs in vitiligo is closely related to the disorder of the metabolic pathway. The KEGG pathway of dysregulation of circRNAs mainly enriched in the biological processes such as ubiquitin mediated proteolysis, endocytosis and RNA degradation, and in Jak-STAT signaling pathway. Therefore, we found the circRNA-miRNA-mRNA regulatory network are involved in the regulation of numerous melanocyte functions, and these dysregulated circRNAs may closely related to the melanocyte metabolism. Our study provides a theoretical basis for studying the vitiligo pathogenesis from the perspective of circRNA-miRNA-mRNA network.

Silicosis is a fatal occupational lung disease which currently has no effective clinical cure. Recent studies examining the underlying mechanism of silicosis have primarily examined experimental models, which may not perfectly reflect the nature of human silicosis progression. A comprehensive profiling of the molecular changes in human silicosis lungs is urgently needed. Here, we conducted RNA sequencing (RNA-seq) on the lung tissues of 10 silicosis patients and 7 non-diseased donors. A total of 2,605 differentially expressed genes (DEGs) and critical pathway changes were identified in human silicosis lungs. Further, the DEGs in silicosis were compared with those in idiopathic pulmonary fibrosis (IPF) and chronic obstructive pulmonary diseases (COPD), to extend current knowledge about the disease mechanisms and develop potential drugs. This analysis revealed both common and specific regulations in silicosis, along with several critical genes (e.g., MUC5AC and FGF10), which are potential drug targets for silicosis treatment. Drugs including Plerixafor and Retinoic acid were predicted as potential candidates in treating silicosis. Overall, this study provides the first transcriptomic fingerprint of human silicosis lungs. The comparative transcriptome analyses comprehensively characterize pathological regulations resulting from silicosis, and provide valuable cues for silicosis treatment.