Transcript Abundance Levels Research Articles

Identifying Relevant Covariates in RNA-seq Analysis by Pseudo-Variable Augmentation

AbstractRNA-sequencing (RNA-seq) technology allows for the identification of differentially expressed genes, which are genes whose mean transcript abundance levels vary across conditions. In practice, RNA-seq datasets often include covariates that are of primary interest in addition to a set of covariates that are subject to selection. Some of these covariates may be relevant to gene expression levels, while others may be irrelevant. Ignoring relevant covariates or attempting to adjust for the effect of irrelevant covariates can compromise the identification of differentially expressed genes. To address this issue, we propose a variable selection method that uses pseudo-variables to control the expected proportion of selected covariates that are irrelevant. Our method accurately selects relevant covariates while keeping the false selection rate below a specified level. We demonstrate that our method outperforms existing methods for detecting differentially expressed genes when working with available covariates. Our method is implemented in function of the R package , which is available at www.github.com/ntyet/csrnaseq. The analysis and simulation are available at www.github.com/ntyet/csrnaseq/tree/main/analysis.

Identification of a specific APOE transcript and functional elements associated with Alzheimer’s disease

BackgroundThe APOE gene is the strongest genetic risk factor for late-onset Alzheimer’s Disease (LOAD). However, the gene regulatory mechanisms at this locus remain incompletely characterized.MethodsTo identify novel AD-linked functional elements within the APOE locus, we integrated SNP variants with multi-omics data from human postmortem brains including 2,179 RNA-seq samples from 3 brain regions and two ancestries (European and African), 667 DNA methylation samples, and ChIP-seq samples. Additionally, we plotted the expression trajectory of APOE transcripts in human brains during development.ResultsWe identified an AD-linked APOE transcript (jxn1.2.2) particularly observed in the dorsolateral prefrontal cortex (DLPFC). The APOE jxn1.2.2 transcript is associated with brain neuropathological features, cognitive impairment, and the presence of the APOE4 allele in DLPFC. We prioritized two independent functional SNPs (rs157580 and rs439401) significantly associated with jxn1.2.2 transcript abundance and DNA methylation levels. These SNPs are located within active chromatin regions and affect brain-related transcription factor-binding affinities. The two SNPs shared effects on the jxn1.2.2 transcript between European and African ethnic groups.ConclusionThe novel APOE functional elements provide potential therapeutic targets with mechanistic insight into the disease etiology.

Tyrosine transfer RNA levels and modifications during blood-feeding and vitellogenesis in the mosquito, Aedes aegypti.

Mosquitoes such as Aedes aegypti must consume a blood meal for the nutrients necessary for egg production. Several transcriptome and proteome changes occur post-blood meal that likely corresponds with codon usage alterations. Transfer RNA (tRNA) is the adapter molecule that reads messenger RNA codons to add the appropriate amino acid during protein synthesis. Chemical modifications to tRNA enhance codon decoding, improving the accuracy and efficiency of protein synthesis. Here, we examined tRNA modifications and transcripts associated with the blood meal and subsequent periods of vitellogenesis in A. aegypti. More specifically, we assessed tRNA transcript abundance and modification levels in the fat body at critical times post blood-feeding. Based on a combination of alternative codon usage and identification of particular modifications, we discovered that increased transcription of tyrosine tRNAs is likely critical during the synthesis of egg yolk proteins in the fat body following a blood meal. Altogether, changes in both the abundance and modification of tRNA are essential factors in the process of vitellogenin production after blood-feeding in mosquitoes.

Neuronal expression of human amyloid-β and Tau drives global phenotypic and multi-omic changes in C. elegans.

Alzheimer's disease (AD) and Alzheimer's related diseases (ADRD) are prevalent age-related neurodegenerative disorders characterized by the accumulation of amyloid-β (Aβ) plaques and Tau neurofibrillary tangles. The nematode Caenorhabditis elegan s ( C. elegans ) serves as an invaluable model organism in diseases of old age-due to its rapid aging. Here we performed an unbiased systems analysis of a C. elegans strain expressing both Aβ and Tau proteins within neurons. We set out to determine if there was a phenotypic interaction between Aβ and Tau. In addition, we were interested in determining the temporal order of the phenotypic and multi-omic (geromic) outcomes. At an early stage of adulthood, we observed reproductive impairments and mitochondrial dysfunction consistent with disruptions in mRNA transcript abundance, protein solubility, and metabolite levels. Notably, the expression of these neurotoxic proteins exhibited a synergistic effect, leading to accelerated aging. Our findings shed light on the close relationship between normal aging and ADRD. Specifically, we demonstrate alterations to metabolic functions preceding age-related neurotoxicity, offering a resource for the development of new therapeutic strategies.

The integral role of de novo lipogenesis in the preparation for seasonal dormancy

Animals can alter their body compositions in anticipation of dormancy to endure seasons with limited food availability. Accumulation of lipid reserves, mostly in the form of triglycerides (TAGs), is observed during the preparation for dormancy in diverse animals, including insects (diapause) and mammals (hibernation). However, the mechanisms involved in the regulation of lipid accumulation and the ecological consequences of failure to accumulate adequate lipid stores in preparation for animal dormancy remain understudied. In the broadest sense, lipid reserves can be accumulated in two ways: the animal either receives lipids directly from the environment or converts the sugars and amino acids present in food to fatty acids through de novo lipogenesis and then to TAGs. Here, we show that preparation for diapause in the Colorado potato beetle (Leptinotarsa decemlineata) involves orchestrated upregulation of genes involved in lipid metabolism with a transcript peak in 8- and 10-d-old diapause-destined insects. Regulation at the transcript abundance level was associated with the accumulation of substantial fat stores. Furthermore, the knockdown of de novo lipogenesis enzymes (ACCase and FAS-1) prolonged the preparatory phase, while the knockdown of fatty acid transportation genes shortened the preparatory phase. Our findings suggest a model in which the insects dynamically decide when to transition from the preparation phase into diapause, depending on the progress in lipid accumulation through de novo lipogenesis.

Assessing the potential relevance of CEACAM6 as a blood transcriptional biomarker.

Changes in blood transcript abundance levels have been associated with pathogenesis in a wide range of diseases. While next generation sequencing technology can measure transcript abundance on a genome-wide scale, downstream clinical applications often require small sets of genes to be selected for inclusion in targeted panels. Here we set out to gather information from the literature and transcriptome datasets that would help researchers determine whether to include the gene CEACAM6 in such panels. We employed a workflow to systematically retrieve, structure, and aggregate information derived from both the literature and public transcriptome datasets. It consisted of profiling the CEACAM6 literature to identify major diseases associated with this candidate gene and establish its relevance as a biomarker. Accessing blood transcriptome datasets identified additional instances where CEACAM6 transcript levels differ in cases vs controls. Finally, the information retrieved throughout this process was captured in a structured format and aggregated in interactive circle packing plots. Although it is not routinely used clinically, the relevance of CEACAM6 as a biomarker has already been well established in the cancer field, where it has invariably been found to be associated with poor prognosis. Focusing on the blood transcriptome literature, we found studies reporting elevated levels of CEACAM6 abundance across a wide range of pathologies, especially diseases where inflammation plays a dominant role, such as asthma, psoriasis, or Parkinson's disease. The screening of public blood transcriptome datasets completed this picture, showing higher abundance levels in patients with infectious diseases caused by viral and bacterial pathogens. Targeted assays measuring CEACAM6 transcript abundance in blood may be of potential utility for the management of patients with diseases presenting with systemic inflammation and for the management of patients with cancer, where the assay could potentially be run both on blood and tumor tissues.

Elevated temperatures dampen the innate immune capacity of developing lake sturgeon (Acipenser fulvescens).

Chronic exposure to high temperatures may leave freshwater fishes vulnerable to opportunistic pathogens, particularly during early life stages. Lake sturgeon, Acipenser fulvescens, populations within the northern expanse of their range in Manitoba, Canada, may be susceptible to high temperature stress and pathogenic infection. We acclimated developing lake sturgeon for 21 days to two ecologically-relevant, summer temperatures (16 and 20°C). Individuals from both acclimation treatments were then exposed to 0, 30, and 60 µg.ml-1 bacterial lipopolysaccharides (endotoxins), as an immune stimulus, for 48 h and sampled 4 and 48 h during trial exposures and following a 7-day recovery period. We then measured whole body transcriptional (mRNA) responses involved in the innate immune, stress, and fatty acid responses following acute exposure to the bacterial endotoxins. Data revealed that overall levels of mRNA transcript abundance were higher in 20°C reared sturgeon under control conditions. However, following exposure to a bacterial stimulus, lake sturgeon acclimated to 16°C produced a more robust and persistent transcriptional response with higher mRNA transcript abundance across innate immune, stress, and fatty acid responses than their 20°C acclimated counterparts. Additional whole-animal performance metrics (critical thermal maximum, metabolic rate, cortisol concentration, and whole-body and mucosal lysozyme activity) demonstrated acclimation-specific responses, indicating compromised metabolic, stress, and enzymatic capacity following the initiation of immune related responses. Our study showed that acclimation to 20°C during early development impaired the immune capacity of developing lake sturgeon as well as the activation of molecular pathways involved in the immune, stress, and fatty acid responses. The present study highlights the effects of ecologically-relevant, chronic thermal stress on seasonal pathogen susceptibility in this endangered species.

Expression of Microcystis Biosynthetic Gene Clusters in Natural Populations Suggests Temporally Dynamic Synthesis of Novel and Known Secondary Metabolites in Western Lake Erie.

Microcystis spp. produce diverse secondary metabolites within freshwater cyanobacterial harmful algal blooms (cyanoHABs) around the world. In addition to the biosynthetic gene clusters (BGCs) encoding known compounds, Microcystis genomes harbor numerous BGCs of unknown function, indicating a poorly understood chemical repertoire. While recent studies show that Microcystis produces several metabolites in the lab and field, little work has focused on analyzing the abundance and expression of its broader suite of BGCs during cyanoHAB events. Here, we use metagenomic and metatranscriptomic approaches to track the relative abundance of Microcystis BGCs and their transcripts throughout the 2014 western Lake Erie cyanoHAB. The results indicate the presence of several transcriptionally active BGCs that are predicted to synthesize both known and novel secondary metabolites. The abundance and expression of these BGCs shifted throughout the bloom, with transcript abundance levels correlating with temperature, nitrate, and phosphorus concentrations and the abundance of co-occurring predatory and competitive eukaryotic microorganisms, suggesting the importance of both abiotic and biotic controls in regulating expression. This work highlights the need for understanding the chemical ecology and potential risks to human and environmental health posed by secondary metabolites that are produced but often unmonitored. It also indicates the prospects for identifying pharmaceutical-like molecules from cyanoHAB-derived BGCs. IMPORTANCE Microcystis spp. dominate cyanobacterial harmful algal blooms (cyanoHABs) worldwide and pose significant threats to water quality through the production of secondary metabolites, many of which are toxic. While the toxicity and biochemistry of microcystins and several other compounds have been studied, the broader suite of secondary metabolites produced by Microcystis remains poorly understood, leaving gaps in our understanding of their impacts on human and ecosystem health. We used community DNA and RNA sequences to track the diversity of genes encoding synthesis of secondary metabolites in natural Microcystis populations and assess patterns of transcription in western Lake Erie cyanoHABs. Our results reveal the presence of both known gene clusters that encode toxic secondary metabolites as well as novel ones that may encode cryptic compounds. This research highlights the need for targeted studies of the secondary metabolite diversity in western Lake Erie, a vital freshwater source to the United States and Canada.

A prognostic signature consisting of N6-methyladenosine modified mRNAs demonstrates clinical potential in prediction of biochemical recurrence and guidance on precision therapy in prostate cancer

Novel biomarkers are urgently needed to improve the prediction of clinical outcomes and guide personalized treatment for prostate cancer (PCa) patients. However, the role of N6-methyladenosine (m6A) modifications in PCa initiation and progression remains largely elusive. In our study, we collected benign Prostate Hyperplasia (BPH), localized PCa, and metastatic PCa samples from patients and performed methylated RNA immunoprecipitation sequencing (MeRIP-Seq) to map m6A-methylated mRNAs. Furthermore, we developed a prognostic signature based on 239 differentially methylated RNAs and the TCGA-PRAD dataset, which can be used to calculate an m6A-modified mRNA (MMM) score for a PCa patient, validated by independent multi-center cohorts. Our findings revealed that differential m6A modifications were positively correlated with altered expressions of mapped m6A-modified mRNAs. Higher MMM scores were associated with shorter times to biochemical recurrence (BCR) in PCa patients, and the MMM scoring system outperformed three well-established signatures in nine independent validation cohorts, as demonstrated by Kaplan-Meier survival analysis, C-index and ROC. Patients who did not respond to androgen receptor signaling inhibitor (ARSI) therapy and immunotherapy were found to have high MMM scores. Two hub genes, TLE1 and PFKL, were confirmed to have m6A sites through MeRIP-qPCR, and their knockdown promoted PCa cell invasion. Bioinformatics analysis of single-cell databases identified cell types with high transcript abundance levels of these two genes. In summary, our study is the first to perform transcriptome-wide m6A mapping in prostate tissues. The translational potential of a prognostic signature, comprising m6A-methylated mRNAs, in predicting clinical outcomes and therapy responses for PCa patients, is demonstrated.

Characterization and genetic diversity of MHC class II DRB genes in the Arabian camel (Camelus dromedarius).

This study investigated the MHC DRB genes in the Arabian camel (Camelus dromedarius). The results revealed the presence of - at least - two transcribed DRB-like genes in chromosome 20, designated MhcCadr-DRB1 and MhcCadr-DRB2. These genes are 155Kb apart, have similar gene structure, and are transcribed in opposite directions. Compared to DRB1, the DRB2 locus contains a deletion of 12 nucleotides in the second exon (270bp), exhibits lower transcript abundance, and is expressed as two splice variants differing by exon 2 skipping. This gene seems to be of minor functional relevance in the dromedary camel. Conversely, the DRB1 is thought to be the main gene in this species showing higher transcript abundance and polymorphism levels. A total of seven DRB1 exon 2 alleles were identified in the Tunisian dromedary camel resulting from 18 amino acid substitutions. Six full length alleles were characterized at the mRNA level. Although there is no clear evidence for balancing selection (i.e., heterozygote advantage), signals of weak historical positive selection acting on the DRB1 gene were detected, as indicated by the limited number of the sites being positively selected. This trend might be related to the low exposure to pathogens and to the demographic history of the species. Comparative analysis with Bactrian and wild camel genomes suggested occurrence of trans species polymorphism (TSP) in the Camelus genus. The results lay the foundation for the MHC DRB1 genetic diversity analysis in this genus since the developed genotyping protocols are fully applicable in the three Camelus species.

PpERF1b-like enhances lignin synthesis in pear (Pyrus pyrifolia) 'hard-end' fruit.

The hard-end is a disorder of pear fruit, however, the mechanisms underlying its development remain unknown. In this study, we found that the hard-end fruit contained a higher transcript abundance level of ethylene-response factor 1b-like (PpERF1b-like) and released more ethylene compared to normal pear. In the ethephon treated normal fruit, flesh tissues accumulated more lignin together with elevated expression of PpERF1b-like. Overexpressing PpERF1b-like transiently in fruit and stably in callus increased lignin accumulation and the expression of lignin biosynthesis genes; the opposite results were observed in fruit showing repressed expression of PpERF1b-like. These results confirmed the role of PpERF1b-like in promoting hard-end formation through promoting lignin synthesis. This study provided valuable information for further clarifying the regulation of hard-end formation in pear.

Assessing the potential relevance of CEACAM6 as a blood transcriptional biomarker

Background Changes in blood transcript abundance levels have been associated with pathogenesis in a wide range of diseases. While next generation sequencing technology can measure transcript abundance on a genome-wide scale, downstream clinical applications often require small sets of genes to be selected for inclusion in targeted panels. Here we set out to gather information from the literature and transcriptome datasets that would help researchers determine whether to include the gene CEACAM6 in such panels. Methods We employed a workflow to systematically retrieve, structure, and aggregate information derived from both the literature and public transcriptome datasets. It consisted of profiling the CEACAM6 literature to identify major diseases associated with this candidate gene and establish its relevance as a biomarker. Accessing blood transcriptome datasets identified additional instances where CEACAM6 transcript levels differ in cases vs controls. Finally, the information retrieved throughout this process was captured in a structured format and aggregated in interactive circle packing plots. Results Although it is not routinely used clinically, the relevance of CEACAM6 as a biomarker has already been well-established in the cancer field, where it has invariably been found to be associated with poor prognosis. Focusing on the blood transcriptome literature, we found studies reporting elevated levels of CEACAM6 abundance across a wide range of pathologies, especially diseases where inflammation plays a dominant role, such as asthma, psoriasis, or Parkinson’s disease. The screening of public blood transcriptome datasets completed this picture, showing higher abundance levels in patients with infectious diseases caused by viral and bacterial pathogens. Conclusions Targeted assays measuring CEACAM6 transcript abundance in blood may be of potential utility for the management of patients with diseases presenting with systemic inflammation and for the management of patients with cancer, where the assay could potentially be run both on blood and tumor tissues.

Investigation of the optimal culture time for warmed bovine ovarian tissues before transplantation†.

Ovarian tissue cryopreservation by vitrification is an effective technique, but there are still many unresolved issues related to the procedure. The aim of this study was to investigate the optimal culture time of postwarmed ovarian tissues and their viability before ovarian tissue transplantation. The bovine ovarian tissues were used to evaluate the effect of postwarming culture periods (0, 0.25, 0.5, 1, 2, 5, and 24h) in the levels of residual cryoprotectant, LDH release, ROS generation, gene and protein abundance, and follicle viability and its mitochondrial membrane potential. Residual cryoprotectant concentration decreased significantly after 1h of culture. The warmed ovarian tissues that underwent between 0 and 2h of culture time showed similar LDH and ROS levels compared with fresh nonfrozen tissues. The anti-Mullerian hormone transcript abundance did not differ in any of the groups. No increase in the relative transcript abundance and protein level of Caspase 3 and Cleaved-Caspase 3, respectively, in the first 2h of culture after warming. On the other hand, an increased protein level of double stranded DNA breaks (gamma-H2AX) was observed in postwarmed tissues disregarding the length of culture time, and a temporary reduction in pan-AKT was detected in postwarming tissues between 0 and 0.25h of culture time. Prolonged culture time lowered the percentage of viable follicles in warmed tissues, but it did not seem to affect the follicular mitochondrial membrane potential. In conclusion, 1-2h of culture time would be optimal for vitrified-warmed tissues before transplantation.

Identifying potential flavonoid biosynthesis regulator in Zanthoxylum bungeanum Maxim. by genome-wide characterization of the MYB transcription factor gene family

Plant MYB transcription factors (TFs) play crucial roles in regulating the biosynthesis of flavonoids but current analysis on their role in Zanthoxylum bungeanum Maxim. (ZBM) is far from comprehensive. In this study, we identified 270 MYB genes in ZBM and divided them into four subfamilies. The R2R3-MYB (ZbMYB) category contained 251 genes and was classified into 33 subfamilies according to their phylogenetic results and sequence similarity. These subfamilies included 24 subgroups containing both MYBs of ZBM plants and AtMYBs, and nine subgroups containing only ZBM MYBs or AtMYBs. ZbMYBs with similar functions clustered into the same subgroup, indicating functional conservation. The subcellular localization analysis predicted that most ZbMYB genes were found in the nucleus. The transposed duplications appeared to play a major role in the expansion of the MYB gene family in ZBM. Through phylogenetic analysis and transcriptome profiling, it was found that 28 ZbMYB genes may regulate the biosynthesis of flavonoids in ZBM, and these genes expression presented distinct temporal and spatial expression patterns. In different fruit development stages of ZBM, the expression patterns of EVM0042160 and EVM0033809 genes obtained by qRT-PCR analysis are very similar to the flavonoid and anthocyanin content curves in ZBM. Further correlation analysis showed that the content of flavonoids in different fruit development stages and the transcript abundance levels of 28 ZbMYB genes have different degrees of correlation relationship. These results indicated that the ZbMYB genes might be involved in the flavonoid metabolic pathway. This comprehensive and systematic analysis of MYB family genes provided a solid foundation for further functional analysis of MYB TFs in ZBM.

A 4 bp InDel in the Promoter of Wheat Gene TaAFP-B Affecting Seed Dormancy Confirmed in Transgenic Rice.

BackgroundWheat (Triticum aestivum L.) ABA insensitive five (ABI5) binding protein gene (TaAFP) is a homologue of the ABI5 binding protein (AFP) gene in Arabidopsis thaliana. It is well documented that AtAFP is a negative regulator of ABA signaling that regulates embryo germination and seed dormancy. TaABI5 was earlier shown to be expressed specifically in seed and its transcript accumulated during wheat grain maturation and acquisition of dormancy. It plays an important role in seed dormancy. In a previous study, we identified two allelic variants TaAFP-B1a and TaAFP-B1b of TaAFP on chromosome arm 2BS in common wheat, designated as, respectively. Sequence analysis revealed a 4 bp insertion in the promoter of TaAFP-B1a compared with TaAFP-B1b that affected mRNA transcription level, mRNA stability, GUS and tdTomatoER translation level, and GUS activity determining seed dormancy.ResultsThe transcription and translation levels of TaAFP-B were significantly reduced in TaAFP-Ba and TaAFP-Ba-GFP transgenic plants compared with TaAFP-Bb and TaAFP-Bb-GFP. The average GI (germination index) values of TaAFP-Ba and TaAFP-Ba-GFP were significantly lower than those of TaAFP-Bb and TaAFP-Bb-GFP in T1 and T2 transgenic rice seeds, whereas mature TaAFP-Ba and TaAFP-Ba-GFP transgenic seeds exhibited increased ABA sensitivity and content of endogenous ABA compared with TaAFP-Bb and TaAFP-Bb-GFP.ConclusionThe 4 bp insertion in the promoter of TaAFP-Ba decreased transcript abundance and translation level in transgenic rice. This insertion increased sensitivity to ABA and content of endogenous ABA in mature seeds, leading to a higher seed dormancy and pre-harvest sprouting tolerance in transgenic rice.

Engineered Fungus Thermothelomyces thermophilus Producing Plant Storage Proteins.

An efficient Agrobacterium-mediated genetic transformation based on the plant binary vector pPZP-RCS2 was carried out for the multiple heterologous protein production in filamentous fungus Thermothelomyces thermophilus F-859 (formerly Myceliophthora thermophila F-859). The engineered fungus Th. thermophilus was able to produce plant storage proteins of Zea mays (α-zein Z19) and Amaranthus hypochondriacus (albumin A1) to enrich fungal biomass by valuable nutritional proteins and improved amino acid content. The mRNA levels of z19 and a1 genes were significantly dependent on their driving promoters: the promoter of tryptophan synthase (PtrpC) was more efficient to express a1, while the promoter of translation elongation factor (Ptef) provided much higher levels of z19 transcript abundance. In general, the total recombinant proteins and amino acid contents were higher in the Ptef-containing clones. This work describes a new strategy to improve mycoprotein nutritive value by overexpression of plant storage proteins.

Genome-Wide Analysis of the Thioredoxin Gene Family in Gossypium hirsutum L. and the Role of the Atypical Thioredoxin Gene GhTRXL3-2 in Flowering

Thioredoxin (TRX) is a highly conserved low-molecular-weight protein and a ubiquitous antioxidant enzyme that plays key role in the regulation of plant growth and development. Here, using the whole-genome sequence, we performed a systematic analysis for the TRX gene family in upland cotton (Gossypium hirsutum L.) and analyzed their structural characteristics, evolution, and expression profiles during growth and development. At least 86 GhTRX members, 40 typical and 46 atypical, were identified in the cotton genome, and they were unevenly distributed on the 26 chromosomes. Conserved domains and phylogenic tree construction classified the typical TRX gene family into seven subfamilies and the atypical TRX into nine subfamilies. An evolutionary analysis revealed that the TRX gene family underwent purification selection during evolution. In addition, an RNA-Seq analysis showed that, during vegetative and reproductive development, the differences in transcript abundance levels and organ-specific expression patterns suggest functional diversity. Biochemical assays demonstrated that the atypical TRX protein GhTRXL3-2 interacted with the cotton FLOWERING LOCUS T protein GhFT. The overexpression of GhTRXL3-2 in Arabidopsis thaliana resulted in early flowering compared with control plants. Additionally, the silencing of GhTRXL3-2 in cotton delayed maturation, suggesting that it has important roles in cotton’s flowering regulation. These results help clarify the evolution of the TRX genes and elucidate their biological functions in cotton flowering regulation.

Sex-stratified genome-wide association study of multisite chronic pain in UK Biobank.

Chronic pain is highly prevalent worldwide and imparts a significant socioeconomic and public health burden. Factors influencing susceptibility to, and mechanisms of, chronic pain development, are not fully understood, but sex is thought to play a significant role, and chronic pain is more prevalent in women than in men. To investigate sex differences in chronic pain, we carried out a sex-stratified genome-wide association study of Multisite Chronic Pain (MCP), a derived chronic pain phenotype, in UK Biobank on 178,556 men and 209,093 women, as well as investigating sex-specific genetic correlations with a range of psychiatric, autoimmune and anthropometric phenotypes and the relationship between sex-specific polygenic risk scores for MCP and chronic widespread pain. We also assessed whether MCP-associated genes showed expression pattern enrichment across tissues. A total of 123 SNPs at five independent loci were significantly associated with MCP in men. In women, a total of 286 genome-wide significant SNPs at ten independent loci were discovered. Meta-analysis of sex-stratified GWAS outputs revealed a further 87 independent associated SNPs. Gene-level analyses revealed sex-specific MCP associations, with 31 genes significantly associated in females, 37 genes associated in males, and a single gene, DCC, associated in both sexes. We found evidence for sex-specific pleiotropy and risk for MCP was found to be associated with chronic widespread pain in a sex-differential manner. Male and female MCP were highly genetically correlated, but at an rg of significantly less than 1 (0.92). All 37 male MCP-associated genes and all but one of 31 female MCP-associated genes were found to be expressed in the dorsal root ganglion, and there was a degree of enrichment for expression in sex-specific tissues. Overall, the findings indicate that sex differences in chronic pain exist at the SNP, gene and transcript abundance level, and highlight possible sex-specific pleiotropy for MCP. Results support the proposition of a strong central nervous-system component to chronic pain in both sexes, additionally highlighting a potential role for the DRG and nociception.

Aedes aegypti miRNA-33 modulates permethrin induced toxicity by regulating VGSC transcripts

Aedes aegypti is a major vector of Zika, dengue, and other arboviruses. Permethrin adulticidal spraying, which targets the voltage-gated sodium channel (VGSC), is commonly done to reduce local mosquito populations and protect humans from exposure to arbovirus pathogens transmitted by this dangerous pest. Permethrin resistance, however, is a growing problem and understanding its underlying molecular basis may identify avenues to combat it. We identified a single G:C polymorphism in pre-miR-33 that was genetically associated with permethrin resistance; resulting isoforms had structural differences that may affect DICER-1/pre-miRNA processing rates. We then assessed the effects of overexpression of pre-miR-33 isoforms on permethrin toxicological phenotypes, VGSC transcript abundance and protein levels for two genetically related mosquito strains. One strain had its naturally high permethrin resistance levels maintained by periodic treatment, and the other was released from selection. VGSC protein levels were lower in the permethrin resistant strain than in the related permethrin-susceptible strain. Overexpression of the G-pre-miR-33 isoform reduced VGSC expression levels in both strains. To further elucidate changes in gene expression associated with permethrin resistance, exome-capture gDNA deep sequencing, genetic association mapping and subsequent gene set enrichment analysis revealed that transport genes, in particular, were selected in resistant versus susceptible mosquitoes. Collectively, these data indicate that miR-33 regulates VGSC expression as part of a nuanced system of neuronal regulation that contributes to a network of heritable features determining permethrin resistance.

Exogenous melatonin maintains leaf quality of postharvest Chinese flowering cabbage by modulating respiratory metabolism and energy status

The efficacy of melatonin on postharvest quality of Chinese flowering cabbage was investigated in this study. Treatment of Chinese cabbage with melatonin significantly delayed leaf yellowing and maintained higher contents of nutrients such as soluble sugars, protein, vitamin C, total flavonoids and phenols during a 7-day storage. Also, melatonin enhanced the activities and expression level of glucose-6-phosphate dehydrogenase (G6PDH) and 6-phosphogluconate dehydrogenase (6-PGDH). Melatonin treatment decreased tissue weight loss, respiration rate, activities of phosphohexose isomerase (PHI), succinate dehydrogenase (SDH), cytochrome C oxidase (CCO), and ascorbic acid oxidase (AAO) and the expression of their corresponding genes. Additionally, melatonin-treated leaves maintained high energy status, as evidenced by higher total ATPase and nicotinamide adenine dinucleotide kinase (NADK) activities, higher energy charge level, adenosine triphosphate (ATP) and adenosine diphosphate (ADP) contents, lower level of adenosine monophosphate (AMP) and transcript abundance of respiratory enzymes. These findings collectively suggest that melatonin maintains leaf quality of postharvest Chinese flowering cabbage might be, at least in part, due to the reduced ratio of Embden-Meyerhof-Parnas pathway (EMP), tricarboxylic acid (TCA) cycle, and cytochrome pathway (CCP), a relatively higher level of pentose phosphate pathway (PPP) and energy status.