MeV+R: using MeV as a graphical user interface for Bioconductor applications in microarray analysis

The Dsup coordinates grain development and abiotic stress in rice

DNA damage is a serious threat to all living organisms and may be induced by environmental stressors. Previous studies have revealed that the tardigrade (Ramazzotius varieornatus) DNA damage suppressor protein Dsup has protective effects in human cells and tobacco. However, whether Dsup provides radiation damage protection more widely in crops is unclear. To explore the effects of Dsup in other crops, stable Dsup overexpression lines through Agrobacterium-mediated transformation were generated and their agronomic traits were deeply investigated. In this study, the overexpression of Dsup not only enhanced the DNA damage resistance at the seeds and seedlings stages, they also exhibited grain size enlargement and starch granule structure and cell size alteration by the scanning electron microscopy observation. Notably, the RNA-seq revealed that the Dsup plants increased radiation-related and abiotic stress-related gene expression in comparison to wild types, suggesting that Dsup is capable to coordinate normal growth and abiotic stress resistance in rice. Immunoprecipitation enrichment with liquid chromatography-tandem mass spectrometry (IP-LC-MS) assays uncovered 21 proteins preferably interacting with Dsup in plants, suggesting that Dsup binds to transcription and translation related proteins to regulate the homeostasis between DNA protection and plant development. In conclusion, our data provide a detailed agronomic analysis of Dsup plants and potential mechanisms of Dsup function in crops. Our findings provide novel insights for the breeding of crop radiation resistance.

A proteomic meta-analysis refinement of plasma extracellular vesicles

Extracellular vesicles play major roles in cell-to-cell communication and are excellent biomarker candidates. However, studying plasma extracellular vesicles is challenging due to contaminants. Here, we performed a proteomics meta-analysis of public data to refine the plasma EV composition by separating EV proteins and contaminants into different clusters. We obtained two clusters with a total of 1717 proteins that were depleted of known contaminants and enriched in EV markers with independently validated 71% true-positive. These clusters had 133 clusters of differentiation (CD) antigens and were enriched with proteins from cell-to-cell communication and signaling. We compared our data with the proteins deposited in PeptideAtlas, making our refined EV protein list a resource for mechanistic and biomarker studies. As a use case example for this resource, we validated the type 1 diabetes biomarker proplatelet basic protein in EVs and showed that it regulates apoptosis of β cells and macrophages, two key players in the disease development. Our approach provides a refinement of the EV composition and a resource for the scientific community.

Proteomic association with age-dependent sex differences in Wisconsin Card Sorting Test performance in healthy Thai subjects

Sex differences in cognitive function exist, but they are not stable and undergo dynamic change during the lifespan. However, our understanding of how sex-related neural information transmission evolves with age is still in its infancy. This study utilized the Wisconsin Card Sorting Test (WCST) and the label-free proteomics method with bioinformatic analysis to investigate the molecular mechanisms underlying age-related sex differences in cognitive performance in 199 healthy Thai subjects (aged 20-70 years), as well as explore the sex-dependent protein complexes for predicting cognitive aging. The results showed that males outperformed females in two of the five WCST sub-scores: %Corrects and %Errors. Sex differences in these scores were related to aging, becoming noticeable in those over 60. At the molecular level, differently expressed individual proteins and protein complexes between both sexes are associated with the potential N-methyl-D-aspartate type glutamate receptor (NMDAR)-mediated excitotoxicity, with the NMDAR complex being enriched exclusively in elderly female samples. These findings provided a preliminary indication that healthy Thai females might be more susceptible to such neurotoxicity, as evidenced by their cognitive performance. NMDAR protein complex enrichment in serum could be proposed as a potential indication for predicting cognitive aging in healthy Thai females.

Pyramiding Breeding of Low-Glutelin-Content Indica Rice with Good Quality and Resistance

Low-glutelin-content rice, a type of functional rice with glutelin levels below 4%, is an essential dietary supplement for chronic kidney disease (CKD) patients. Developing low-glutelin-content rice varieties is crucial to catering to the growing CKD population. In this study, we aimed to create a new low-glutelin indica rice variety with excellent agronomic traits. To achieve this, we employed a combination of molecular-marker-assisted selection and traditional breeding techniques. The cultivars W3660, Wushansimiao (WSSM), and Nantaixiangzhan (NTXZ) were crossbred, incorporating the Lgc-1, Pi-2, Xa23, and fgr alleles into a single line. The result of this breeding effort was "Yishenxiangsimiao", a new indica rice variety that inherits the desirable characteristics of its parent lines. Yishenxiangsimiao (YSXSM) possesses not only a low glutelin content but also dual resistance to blast and bacterial blight (BB). It exhibits high-quality grains with a fragrant aroma. This new low-glutelin indica cultivar not only ensures a stable food supply for CKD patients but also serves as a healthy dietary option for the general public. We also performed RNA-seq of these rice varieties to investigate their internal gene expression differences. The YSXSM exhibited a higher biotic-resistance gene expression in comparison to NTXZ. In summary, we successfully developed a novel low-glutelin indica rice variety, "Yishenxiangsimiao", with superior agronomic traits. This rice variety addresses the dietary needs of CKD patients and offers a nutritious choice for all consumers.

Cholinergic-estrogen interaction is associated with the effect of education on attenuating cognitive sex differences in a Thai healthy population

The development of human brain is shaped by both genetic and environmental factors. Sex differences in cognitive function have been found in humans as a result of sexual dimorphism in neural information transmission. Numerous studies have reported the positive effects of education on cognitive functions. However, little work has investigated the effect of education on attenuating cognitive sex differences and the neural mechanisms behind it based on healthy population. In this study, the Wisconsin Card Sorting Test (WCST) was employed to examine sex differences in cognitive function in 135 Thai healthy subjects, and label-free quantitative proteomic method and bioinformatic analysis were used to study sex-specific neurotransmission-related protein expression profiles. The results showed sex differences in two WCST sub-scores: percentage of Total corrects and Total errors in the primary education group (Bayes factor>100) with males performed better, while such differences eliminated in secondary and tertiary education levels. Moreover, 11 differentially expressed proteins (DEPs) between men and women (FDR<0.1) were presented in both education groups, with majority of them upregulated in females. Half of those DEPs interacted directly with nAChR3, whereas the other DEPs were indirectly connected to the cholinergic pathways through interaction with estrogen. These findings provided a preliminary indication that a cholinergic-estrogen interaction relates to, and might underpin, the effect of education on attenuating cognitive sex differences in a Thai healthy population.

Deficiency of the circadian clock gene Rev-erbα induces mood disorder-like behaviours and dysregulation of the serotonergic system in mice

Mood disorders such as depression, anxiety, and bipolar disorder are highly associated with disrupted daily rhythms of activity, which are often observed in shift work and sleep disturbance in humans. Recent studies have proposed the REV-ERBα protein as a key circadian nuclear receptor that links behavioural rhythms to mood regulation. However, how the Rev-erbα gene participates in the regulation of mood remains poorly understood. Here, we show that the regulation of the serotonergic (5-HTergic) system, which plays a central role in stress-induced mood behaviours, is markedly disrupted in Rev-erbα^-/- mice. Rev-erbα^-/- mice exhibit both negative and positive behavioural phenotypes, including anxiety-like and mania-like behaviours, when subjected to a stressful environment. Importantly, Rev-erbα^-/- mice show a significant decrease in the expression of a gene that encodes the rate-limiting enzyme of serotonin (5-HT) synthesis in the raphe nuclei (RN). In addition, 5-HT levels in Rev-erbα^-/- mice are significantly reduced in the prefrontal cortex, which receives strong inputs from the RN and controls stress-related behaviours. Our findings indicate that Rev-erbα plays an important role in controlling the 5-HTergic system and thus regulates mood and behaviour.

A review of bioinformatics tools and web servers in different microarray platforms used in cancer research

Over the past decade, conventional lab work strategies have gradually shifted from being limited to a laboratory setting towards a bioinformatics era to help manage and process the vast amounts of data generated by omics technologies. The present work outlines the latest contributions of bioinformatics in analyzing microarray data and their application to cancer. We dissect different microarray platforms and their use in gene expression in cancer models. We highlight how computational advances empowered the microarray technology in gene expression analysis. The study on protein-protein interaction databases classified into primary, derived, meta-database, and prediction databases describes the strategies to curate and predict novel interaction networks in silico. In addition, we summarize the areas of bioinformatics where neural graph networks are currently being used, such as protein functions, protein interaction prediction, and in silico drug discovery and development. We also discuss the role of deep learning as a potential tool in the prognosis, diagnosis, and treatment of cancer. Integrating these resources efficiently, practically, and ethically is likely to be the most challenging task for the healthcare industry over the next decade; however, we believe that it is achievable in the long term.

Proteomic landscape subtype and clinical prognosis of patients with the cognitive impairment by Japanese encephalitis infection

BACKGROUND

Cognitive impairment is one of the primary sequelae affecting the quality of life of patients with Japanese encephalitis (JE). The clinical treatment is mainly focused on life support, lacking of targeted treatment strategy.

METHODS

A cerebrospinal fluid (CSF) proteomic profiling study was performed including 26 patients with JE in Gansu province of China from June 2017 to October 2018 and 33 other concurrent hospitalized patients who were excluded central nervous system (CNS) organic or CNS infection diseases. The clinical and proteomics data of patients with JE were undergoing combined analysis for the first time.

RESULTS

Two subtypes of JE associated with significantly different prognoses were identified. Compared to JE1, the JE2 subtype is associated with lower overall survival rate and a higher risk of cognitive impairment. The percentages of neutrophils (N%), lymphocyte (L%), and monocytes (M%) decreased in JE2 significantly.

CONCLUSIONS

The differences in proteomic landscape between JE subgroups have specificity for the prognosis of cognitive impairment. The data also provided some potential target proteins for treatment of cognitive impairments caused by JE. Trial registration ChiCTR, ChiCTR2000030499. Registered 1st June 2017, http://www.medresman.org.cn/pub/cn/proj/projectshow.aspx?proj=6333.

Microscopic and Transcriptomic Analyses of Dalbergoid Legume Peanut Reveal a Divergent Evolution Leading to Nod-Factor-Dependent Epidermal Crack-Entry and Terminal Bacteroid Differentiation

Root nodule symbiosis (RNS) is the pillar behind sustainable agriculture and plays a pivotal role in the environmental nitrogen cycle. Most of the genetic, molecular, and cell-biological knowledge on RNS comes from model legumes that exhibit a root-hair mode of bacterial infection, in contrast to the Dalbergoid legumes exhibiting crack-entry of rhizobia. As a step toward understanding this important group of legumes, we have combined microscopic analysis and temporal transcriptome to obtain a dynamic view of plant gene expression during Arachis hypogaea (peanut) nodule development. We generated comprehensive transcriptome data by mapping the reads to A. hypogaea, and two diploid progenitor genomes. Additionally, we performed BLAST searches to identify nodule-induced yet-to-be annotated peanut genes. Comparison between peanut, Medicago truncatula, Lotus japonicus, and Glycine max showed upregulation of 61 peanut orthologs among 111 tested known RNS-related genes, indicating conservation in mechanisms of nodule development among members of the Papilionoid family. Unlike model legumes, recruitment of class 1 phytoglobin-derived symbiotic hemoglobin (SymH) in peanut indicates diversification of oxygen-scavenging mechanisms in the Papilionoid family. Finally, the absence of cysteine-rich motif-1-containing nodule-specific cysteine-rich peptide (NCR) genes but the recruitment of defensin-like NCRs suggest a diverse molecular mechanism of terminal bacteroid differentiation. In summary, our work describes genetic conservation and diversification in legume-rhizobia symbiosis in the Papilionoid family, as well as among members of the Dalbergoid legumes.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Genome-wide identification of ZmSnRK2 genes and functional analysis of ZmSnRK2.10 in ABA signaling pathway in maize (Zea mays L)

BACKGROUND

Phytohormone abscisic acid (ABA) is involved in the regulation of a wide range of biological processes. In Arabidopsis, it has been well-known that SnRK2s are the central components of the ABA signaling pathway that control the balance between plant growth and stress response, but the functions of ZmSnRK2 in maize are rarely reported. Therefore, the study of ZmSnRK2 is of great importance to understand the ABA signaling pathways in maize.

RESULTS

In this study, 14 ZmSnRK2 genes were identified in the latest version of maize genome database. Phylogenetic analysis revealed that ZmSnRK2s are divided into three subclasses based on their diversity of C-terminal domains. The exon-intron structures, phylogenetic, synteny and collinearity analysis indicated that SnRK2s, especially the subclass III of SnRK2, are evolutionally conserved in maize, rice and Arabidopsis. Subcellular localization showed that ZmSnRK2 proteins are localized in the nucleus and cytoplasm. The RNA-Seq datasets and qRT-PCR analysis showed that ZmSnRK2 genes exhibit spatial and temporal expression patterns during the growth and development of different maize tissues, and the transcript levels of some ZmSnRK2 genes in kernel are significantly induced by ABA and sucrose treatment. In addition, we found that ZmSnRK2.10, which belongs to subclass III, is highly expressed in kernel and activated by ABA. Overexpression of ZmSnRK2.10 partially rescued the ABA-insensitive phenotype of snrk2.2/2.3 double and snrk2.2/2.3/2.6 triple mutants and led to delaying plant flowering in Arabidopsis.

CONCLUSION

The SnRK2 gene family exhibits a high evolutionary conservation and has expanded with whole-genome duplication events in plants. The ZmSnRK2s expanded in maize with whole-genome and segmental duplication, not tandem duplication. The expression pattern analysis of ZmSnRK2s in maize offers important information to study their functions. Study of the functions of ZmSnRK.10 in Arabidopsis suggests that the ABA-dependent members of SnRK2s are evolutionarily conserved in plants. Our study elucidated the structure and evolution of SnRK2 genes in plants and provided a basis for the functional study of ZmSnRK2s protein in maize.

Oocyte-specific linker histone H1foo interacts with Esrrb to induce chromatin decondensation at specific gene loci

Linker histone H1 is mainly localized in the linker DNA region, between two nucleosome cores, and regulates chromatin structures linking gene expression. Mammalian oocytes contain the histone H1foo, a distinct member with low sequence similarity to other members in the H1 histone family. Although, from various previous studies, evidence related to H1foo function in chromatin structures is being accumulated, the distribution of H1foo at the target gene loci in a genome-wide manner and the molecular mechanism of H1foo-dependent chromatin architecture remain unclear. In this study, we aimed to identify the target loci and the physiological factor bound to H1foo at the loci. Chromatin immunoprecipitation sequencing analysis of H1foo-overexpressing mouse embryonic stem cells showed that H1foo is enriched around the transcriptional start sites of genes such as oocyte-specific genes and that the chromatin structures at these regions were relaxed. We demonstrated that H1foo was physiologically bound to the nuclear receptor estrogen-related receptor beta (Esrrb), and Esrrb was necessary for H1foo activity of chromatin decondensation at the target loci. The specific localization and interaction with Esrrb were validated in endogenous H1foo of oocytes. Overall, H1foo induces chromatin decondensation in a locus-specific manner and this function is achieved by interacting with Esrrb.

Site-specific glycoproteomic analysis revealing increased core-fucosylation on FOLR1 enhances folate uptake capacity of HCC cells to promote EMT

Rationale: Epithelial-mesenchymal transition (EMT) has been recognized as an important step toward high invasion and metastasis of many cancers including hepatocellular carcinoma (HCC), while the mechanism for EMT promotion is still ambiguous.

Methods: The dynamic alterations of site-specific glycosylation during HGF/TGF-β1-induced EMT process of three HCC cell lines were systematically investigated using precision glycoproteomic methods. The possible roles of EMT-related glycoproteins and site-specific glycans were further confirmed by various molecular biological approaches.

Results: Using mass spectrometry-based glycoproteomic methods, we totally identified 2306 unique intact glycopeptides from SMMC-7721 and HepG2 cell lines, and found that core-fucosylated glycans were accounted for the largest proportion of complex N-glycans. Through quantification analysis of intact glycopeptides, we found that the majority of core-fucosylated intact glycopeptides from folate receptor α (FOLR1) were up-regulated in the three HGF-treated cell lines. Similarly, core-fucosylation of FOLR1 were up-regulated in SMMC-7721 and Hep3B cells with TGF-β1 treatment. Using molecular approaches, we further demonstrated that FUT8 was a driver for HGF/TGF-β1-induced EMT. The silencing of FUT8 reduced core-fucosylation and partially blocked the progress of HGF-induced EMT. Finally, we confirmed that the level of core-fucosylation on FOLR1 especially at the glycosite Asn-201 positively regulated the cellular uptake capacity of folates, and enhanced uptake of folates could promote the EMT of HCC cells.

Conclusions: Based on the results, we proposed a potential pathway for HGF or TGF-β1-induced EMT of HCC cells: HGF or TGF-β1 treatment of HCC cells can increase the expression of glycosyltransferase FUT8 to up-regulate the core-fucosylation of N-glycans on glycoproteins including the FOLR1; core-fucosylation on FOLR1 can then enhance the folate uptake capacity to finally promote the EMT progress of HCC cells.

SARS-CoV-2 Proteome Microarray for Mapping COVID-19 Antibody Interactions at Amino Acid Resolution

Comprehensive profiling of humoral antibody response to severe acute respiratory syndrome (SARS) coronavirus-2 (CoV-2) proteins is essential in understanding the host immunity and in developing diagnostic tests and vaccines. To address this concern, we developed a SARS-CoV-2 proteome peptide microarray to analyze antibody interactions at the amino acid resolution. With the array, we demonstrate the feasibility of employing SARS-CoV-1 antibodies to detect the SARS-CoV-2 nucleocapsid phosphoprotein. The first landscape of B-cell epitopes for SARS-CoV-2 IgM and IgG antibodies in the serum of 10 coronavirus disease of 2019 (COVID-19) patients with early infection is also constructed. With array data and structural analysis, a peptide epitope for neutralizing antibodies within the SARS-CoV-2 spike receptor-binding domain's interaction interface with the angiotensin-converting enzyme 2 receptor was predicted. All the results demonstrate the utility of our microarray as a platform to determine the changes of antibody responses in COVID-19 patients and animal models as well as to identify potential targets for diagnosis and treatment.

Phosphate-Starvation-Inducible S-Like RNase Genes in Rice Are Involved in Phosphate Source Recycling by RNA Decay

The fine-tuning of inorganic phosphate (Pi) for enhanced use efficiency has long been a challenging subject in agriculture, particularly in regard to rice as a major crop plant. Among ribonucleases (RNases), the RNase T2 family is broadly distributed across kingdoms, but little has been known on its substrate specificity compared to RNase A and RNase T1 families. Class I and class II of the RNase T2 family are defined as the S-like RNase (RNS) family and have showed the connection to Pi recycling in Arabidopsis. In this study, we first carried out a phylogenetic analysis of eight rice and five Arabidopsis RNS genes and identified mono-specific class I and dicot-specific class I RNS genes, suggesting the possibility of functional diversity between class I RNS family members in monocot and dicot species through evolution. We then compared the in silico expression patterns of all RNS genes in rice and Arabidopsis under normal and Pi-deficient conditions and further confirmed the expression patterns of rice RNS genes via qRT-PCR analysis. Subsequently, we found that most of the OsRNS genes were differentially regulated under Pi-deficient treatment. Association of Pi recycling by RNase activity in rice was confirmed by measuring total RNA concentration and ribonuclease activity of shoot and root samples under Pi-sufficient or Pi-deficient treatment during 21 days. The total RNA concentrations were decreased by < 60% in shoots and < 80% in roots under Pi starvation, respectively, while ribonuclease activity increased correspondingly. We further elucidate the signaling pathway of Pi starvation through upregulation of the OsRNS genes. The 2-kb promoter region of all OsRNS genes with inducible expression patterns under Pi deficiency contains a high frequency of P1BS cis-acting regulatory element (CRE) known as the OsPHR2 binding site, suggesting that the OsRNS family is likely to be controlled by OsPHR2. Finally, the dynamic transcriptional regulation of OsRNS genes by overexpression of OsPHR2, ospho2 mutant, and overexpression of OsPT1 lines involved in Pi signaling pathway suggests the molecular basis of OsRNS family in Pi recycling via RNA decay under Pi starvation.

Serum Protein Profiling Reveals a Landscape of Inflammation and Immune Signaling in Early-stage COVID-19 Infection

Coronavirus disease 2019 (COVID-19) is a highly contagious infection and threating the human lives in the world. The elevation of cytokines in blood is crucial to induce cytokine storm and immunosuppression in the transition of severity in COVID-19 patients. However, the comprehensive changes of serum proteins in COVID-19 patients throughout the SARS-CoV-2 infection is unknown. In this work, we developed a high-density antibody microarray and performed an in-depth proteomics analysis of serum samples collected from early COVID-19 (n = 15) and influenza (n = 13) patients. We identified a large set of differentially expressed proteins (n = 132) that participate in a landscape of inflammation and immune signaling related to the SARS-CoV-2 infection. Furthermore, the significant correlations of neutrophil and lymphocyte with the CCL2 and CXCL10 mediated cytokine signaling pathways was identified. These information are valuable for the understanding of COVID-19 pathogenesis, identification of biomarkers and development of the optimal anti-inflammation therapy.

Identification of TYROBP and C1QB as Two Novel Key Genes With Prognostic Value in Gastric Cancer by Network Analysis

Background: Gastric cancer (GC) is the fifth most frequently diagnosed malignancy, and the third leading cause of tumor-related mortalities worldwide. Due to a high heterogeneity in GC, its treatment and prognosis are challenging, necessitating urgent identification of novel prognostic predictors for GC patients. Methods: We downloaded RNA sequence data, from the Cancer Genome Atlas and microarray data from Gene Expression Omnibus database, then identified common differentially-expressed genes (DEGs) between GC and normal gastric tissues across four datasets. We then used a combination of protein-protein interaction (PPI) network and weighted gene co-expression network analysis (WGCNA) to identify key genes with prognostic value in GC. Thereafter, we used quantitative real time polymerase chain reaction (qRT-PCR) to validate expression of the identified key genes in the Zhejiang University (ZJU) cohort. Finally, we evaluated the relationships between gene expression and immune factors, including immune cells and biomarkers of immunotherapy. Results: Among 426 common DEGs screened, 333 and 93 were upregulated and downregulated, respectively. PPI network and WGCNA successfully identified the top 30 hub genes, among which PTPRC, TYROBP, CCR1, CYBB, LCP2, and C1QB were common. Furthermore, TYROBP and C1QB were negatively associated with prognosis of GC patients, implying that they were key GC predictors. Interestingly, TYROBP and C1QB were positively correlated with predictive biomarkers for GC immunotherapy, including PD-L1 expression, CD8⁺ T cells infiltration, and EBV status. Conclusions: TYROBP and C1QB were identified as two novel key genes with prognostic value in GC by network analysis.

ZEITLUPE facilitates the rhythmic movements of Nicotiana attenuata flowers

Circadian organ movements are ubiquitous in plants. These rhythmic outputs are thought to be regulated by the circadian clock and auxin signalling, but the underlying mechanisms have not been clarified. Flowers of Nicotiana attenuata change their orientation during the daytime through a 140° arc to balance the need for pollinators and the protection of their reproductive organs. This rhythmic trait is under the control of the circadian clock and results from bending and re-straightening movements of the pedicel, stems that connect flowers to the inflorescence. Using an explant system that allowed pedicel growth and curvature responses to be characterized with high spatial and temporal resolution, we demonstrated that this movement is organ autonomous and mediated by auxin. Changes in the growth curvature of the pedicel are accompanied by an auxin gradient and dorsiventral asymmetry in auxin-dependent transcriptional responses; application of auxin transport inhibitors influenced the normal movements of this organ. Silencing the expression of the circadian clock component ZEITLUPE (ZTL) arrested changes in the growth curvature of the pedicel and altered auxin signalling and responses. IAA19-like, an Aux/IAA transcriptional repressor that is circadian regulated and differentially expressed between opposite tissues of the pedicel, and therefore possibly involved in the regulation of changes in organ curvature, physically interacted with ZTL. Together, these results are consistent with a direct link between the circadian clock and the auxin signalling pathway in the regulation of this rhythmic floral movement.

Systematic re-evaluation of the long-used standard protocol of urease-dependent metabolome sample preparation

In the urinary metabolomics for finding biomarkers in urine, owing to high concentrations of urea, for chromatography-based metabolomic analysis, urea needed to be degraded by urease. This urease pretreatment has been the key step of sample preparation for standard urinary metabolomics until today even for mass spectrometry-based analysis. The urease pretreatment involving incubation of urine with urease contradicts the concept of metabolome sampling, which should immediately arrest metabolic reactions to prevent alterations of a metabolite profile. Nonetheless, the impact of urease pretreatment has not been clearly elucidated yet. We found that activities of urease and endogenous urinary enzymes and metabolite contaminants from the urease preparations introduce artefacts into metabolite profiles, thus leading to misinterpretation.

Temporal Profile of Brain Gene Expression After Prey Catching Conditioning in an Anuran Amphibian

A key goal in modern neurobiology is to understand the mechanisms underlying learning and memory. To that end, it is essential to identify the patterns of gene expression and the temporal sequence of molecular events associated with learning and memory processes. It is also important to ascertain if and how these molecular events vary between organisms. In vertebrates, learning and memory processes are characterized by distinct phases of molecular activity involving gene transcription, structural change, and long-term maintenance of such structural change in the nervous system. Utilizing next generation sequencing techniques, we profiled the temporal expression patterns of genes in the brain of the fire-bellied toad Bombina orientalis after prey catching conditioning. The fire-bellied toad is a basal tetrapod whose neural architecture and molecular pathways may help us understand the ancestral state of learning and memory mechanisms in tetrapods. Differential gene expression following conditioning revealed activity in molecular pathways related to immediate early genes (IEG), cytoskeletal modification, axon guidance activity, and apoptotic processes. Conditioning induced early IEG activity coinciding with transcriptional activity and neuron structural modification, followed by axon guidance and cell adhesion activity, and late neuronal pruning. While some of these gene expression patterns are similar to those found in mammals submitted to conditioning, some interesting divergent expression profiles were seen, and differential expression of some well-known learning-related mammalian genes is missing altogether. These results highlight the importance of using a comparative approach in the study of the mechanisms of leaning and memory and provide molecular resources for a novel vertebrate model in the relatively poorly studied Amphibia.

Kynurenine, 3-OH-kynurenine, and anthranilate are nutrient metabolites that alter H3K4 trimethylation and H2AS40 O-GlcNAcylation at hypothalamus-related loci

Epigenetic mechanisms can establish and maintain mitotically stable patterns of gene expression while retaining the DNA sequence. These mechanisms can be affected by environmental factors such as nutrients. The importance of intracellular dosages of nutrient metabolites such as acetyl coenzyme A and S-adenosylmethionine, which are utilized as donors for post-translational modifications, is well-known in epigenetic regulation; however, the significance of indirect metabolites in epigenetic regulation is not clear. In this study, we screened for metabolites that function as epigenetic modulators. Because the expression of genes related to hypothalamic function is reportedly affected by nutritional conditions, we used a neural cell culture system and evaluated hypothalamic-linked loci. We supplemented the culture medium with 129 metabolites separately during induction of human-iPS-derived neural cells and used high-throughput ChIP-qPCR to determine the epigenetic status at 37 hypothalamus-linked loci. We found three metabolites (kynurenine, 3-OH-kynurenine, and anthranilate) from tryptophan pathways that increased H3K4 trimethylation and H2AS40 O-GlcNAcylation, resulting in upregulated gene expression at most loci, except those encoding pan-neural markers. Dietary supplementation of these three metabolites and the resulting epigenetic modification were important for stability in gene expression. In conclusion, our findings provide a better understanding of how nutrients play a role in epigenetic mechanisms.

An Integrative Approach to Analyze Seed Germination in Brassica napus

Seed germination is a complex trait determined by the interaction of hormonal, metabolic, genetic, and environmental components. Variability of this trait in crops has a big impact on seedling establishment and yield in the field. Classical studies of this trait in crops have focused mainly on the analyses of one level of regulation in the cascade of events leading to seed germination. We have carried out an integrative and extensive approach to deepen our understanding of seed germination in Brassica napus by generating transcriptomic, metabolic, and hormonal data at different stages upon seed imbibition. Deep phenotyping of different seed germination-associated traits in six winter-type B. napus accessions has revealed that seed germination kinetics, in particular seed germination speed, are major contributors to the variability of this trait. Metabolic profiling of these accessions has allowed us to describe a common pattern of metabolic change and to identify the levels of malate and aspartate metabolites as putative metabolic markers to estimate germination performance. Additionally, analysis of seed content of different hormones suggests that hormonal balance between ABA, GA, and IAA at crucial time points during this process might underlie seed germination differences in these accessions. In this study, we have also defined the major transcriptome changes accompanying the germination process in B. napus. Furthermore, we have observed that earlier activation of key germination regulatory genes seems to generate the differences in germination speed observed between accessions in B. napus. Finally, we have found that protein-protein interactions between some of these key regulator are conserved in B. napus, suggesting a shared regulatory network with other plant species. Altogether, our results provide a comprehensive and detailed picture of seed germination dynamics in oilseed rape. This new framework will be extremely valuable not only to evaluate germination performance of B. napus accessions but also to identify key targets for crop improvement in this important process.

KIAA1199, a Target of MicoRNA-486-5p, Promotes Papillary Thyroid Cancer Invasion by Influencing Epithelial-Mesenchymal Transition (EMT)

Background

KIAA1199 has been reported to be associated with malignant progression and poor clinical outcomes in various human malignancies. However, its clinical role and molecular function remain unknown in papillary thyroid cancer (PTC).

Material/Methods

The Cancer Genome Atlas (TCGA) was used to investigate the expression profiles of KIAA1199 and miR-486-5p in PTC. Immunohistochemistry was used to validate the protein expression of KIAA1199 in PTC. The Weighted Gene Co-expression Network Analysis (WGCNA) and Gene Set Enrichment Analysis (GSEA) were used to explore the potential pathway underling KIAA1199 in PTC. In vitro and in vivo experiments were performed to investigate the biological role of KIAA1199 in PTC progression. Luciferase reporter assays and Western blot analysis were performed to determine whether KIAA1199 is a downstream target of miR-486-5p.

Results

We found that KIAA1199 was aberrantly elevated in PTC tissues compared with normal tissues, and upregulation of KIAA1199 was positively correlated with more advanced clinical variables. Additionally, bioinformatic analysis indicated that KIAA1199 was involved in cell migration and invasion. KIAA1199 silencing inhibited the invasive ability of PTC cells by affecting epithelial-mesenchymal transition (EMT) in vitro and in vivo. Furthermore, miR-486-5p was identified as an upstream microRNA that directly targets the 3′-UTR region of KIAA1199.

Conclusions

The miR-486-5p/KIAA1199/EMT axis might play a critical role in PTC invasion and metastasis and offers a potential therapeutic strategy for PTC.

miR-126-3p down-regulation contributes to dabrafenib acquired resistance in melanoma by up-regulating ADAM9 and VEGF-A

Background

Development of resistance to inhibitors of BRAF (BRAFi) and MEK (MEKi) remains a great challenge for targeted therapy in patients with BRAF-mutant melanoma. Here, we explored the role of miRNAs in melanoma acquired resistance to BRAFi.

Methods

miRNA expression in two BRAF-mutant melanoma cell lines and their dabrafenib-resistant sublines was determined using Affymetrix GeneChip® miRNA 3.1 microarrays and/or qRT-PCR. The effects of miR-126-3p re-expression on proliferation, apoptosis, cell cycle, ERK1/2 and AKT phosphorylation, dabrafenib sensitivity, invasiveness and VEGF-A secretion were evaluated in the dabrafenib-resistant sublines using MTT assays, flow cytometry, immunoblotting, invasion assays in Boyden chambers and ELISA. ADAM9, PIK3R2, MMP7 and CXCR4 expression in the sensitive and dabrafenib-resistant cells was determined by immunoblotting. Small RNA interference was performed to investigate the consequence of VEGFA or ADAM9 silencing on proliferation, invasiveness or dabrafenib sensitivity of the resistant sublines. Long-term proliferation assays were carried out in dabrafenib-sensitive cells to assess the effects of enforced miR-126-3p expression or ADAM9 silencing on resistance development. VEGF-A serum levels in melanoma patients treated with BRAFi or BRAFi+MEKi were evaluated at baseline (T0), after two months of treatment (T2) and at progression (TP) by ELISA.

Results

miR-126-3p was significantly down-regulated in the dabrafenib-resistant sublines as compared with their parental counterparts. miR-126-3p replacement in the drug-resistant cells inhibited proliferation, cell cycle progression, phosphorylation of ERK1/2 and/or AKT, invasiveness, VEGF-A and ADAM9 expression, and increased dabrafenib sensitivity. VEGFA or ADAM9 silencing impaired proliferation and invasiveness of the drug-resistant sublines. ADAM9 knock-down in the resistant cells increased dabrafenib sensitivity, whereas miR-126-3p enforced expression or ADAM9 silencing in the drug-sensitive cells delayed the development of resistance. At T0 and T2, statistically significant differences were observed in VEGF-A serum levels between patients who responded to therapy and patients who did not. In responder patients, a significant increase of VEGF-A levels was observed at TP versus T2.

Conclusions

Strategies restoring miR-126-3p expression or targeting VEGF-A or ADAM9 could restrain growth and metastasis of dabrafenib-resistant melanomas and increase their drug sensitivity. Circulating VEGF-A is a promising biomarker for predicting patients’ response to BRAFi or BRAFi+MEKi and for monitoring the onset of resistance.

Electronic supplementary material

The online version of this article (10.1186/s13046-019-1238-4) contains supplementary material, which is available to authorized users.

Metabolite profile changes and increased antioxidative and antiinflammatory activities of mixed vegetables after fermentation by Lactobacillus plantarum

Fermented vegetables have emerged as prebiotics with various health benefits. However, the possible mechanisms behind their health benefits are unclear. To relate the metabolite profile changes in fermented mixed vegetables with associated health benefits of fermented vegetables, we analyzed the metabolite profiles of mixed vegetables, before and after fermentation by Lactobacillus plantarum, using gas chromatography/time-of-flight-mass spectrometry (GC/TOF-MS). To analyze health benefits of fermented vegetables, antioxidative and antiinflammatory activities were measured using RAW 264.7 cells. Among 78 metabolites identified by GC/TOF-MS in this study, those significantly increased after fermentation include antioxidative and/or antiinflammatory agents such as lactate, 3-phennyllactate, indole-3-lactate, β-hydroxybutyrate, γ-aminobutyrate, and glycerol. These metabolites may have been either newly synthesized or depolymerized from high molecular weight polymers from vegetables during fermentation. This is the first metabolomics study to relate metabolite profile changes with increased health benefits of fermented vegetables.

Aspergillus fumigatus Challenged by Human Dendritic Cells: Metabolic and Regulatory Pathway Responses Testify a Tight Battle

Dendritic cells (DCs) are antigen presenting cells which serve as a passage between the innate and the acquired immunity. Aspergillosis is a major lethal condition in immunocompromised patients caused by the adaptable saprophytic fungus Aspergillus fumigatus. The healthy human immune system is capable to ward off A. fumigatus infections however immune-deficient patients are highly vulnerable to invasive aspergillosis. A. fumigatus can persist during infection due to its ability to survive the immune response of human DCs. Therefore, the study of the metabolism specific to the context of infection may allow us to gain insight into the adaptation strategies of both the pathogen and the immune cells. We established a metabolic model of A. fumigatus central metabolism during infection of DCs and calculated the metabolic pathway (elementary modes; EMs). Transcriptome data were used to identify pathways activated when A. fumigatus is challenged with DCs. In particular, amino acid metabolic pathways, alternative carbon metabolic pathways and stress regulating enzymes were found to be active. Metabolic flux modeling identified further active enzymes such as alcohol dehydrogenase, inositol oxygenase and GTP cyclohydrolase participating in different stress responses in A. fumigatus. These were further validated by qRT-PCR from RNA extracted under these different conditions. For DCs, we outlined the activation of metabolic pathways in response to the confrontation with A. fumigatus. We found the fatty acid metabolism plays a crucial role, along with other metabolic changes. The gene expression data and their analysis illuminate additional regulatory pathways activated in the DCs apart from interleukin regulation. In particular, Toll-like receptor signaling, NOD-like receptor signaling and RIG-I-like receptor signaling were active pathways. Moreover, we identified subnetworks and several novel key regulators such as UBC, EGFR, and CUL3 of DCs to be activated in response to A. fumigatus. In conclusion, we analyze the metabolic and regulatory responses of A. fumigatus and DCs when confronted with each other.

miRNA profiling of urinary exosomes to assess the progression of acute kidney injury

Because exosomes have gained attention as a source of biomarkers, we investigated if miRNAs in exosomes (exo-miRs) can report the disease progression of organ injury. Using rat renal ischemia-reperfusion injury (IRI) as a model of acute kidney injury (AKI), we determined temporally-released exo-miRs in urine during IRI and found that these exo-miRs could reliably mirror the progression of AKI. From the longitudinal measurements of miRNA expression in kidney and urine, we found that release of exo- miRs was a regulated sorting process. In the injury state, miR-16, miR-24, and miR-200c were increased in the urine. Interestingly, expression of target mRNAs of these exo-miRs was significantly altered in renal medulla. Next, in the early recovery state, exo-miRs (miR-9a, miR-141, miR-200a, miR-200c, miR-429), which share Zeb1/2 as a common target mRNA, were upregulated together, indicating that they reflect TGF-β-associated renal fibrosis. Finally, release of exo-miRs (miR-125a, miR-351) was regulated by TGF-β1 and was able to differentiate the sham and IRI even after the injured kidneys were recovered. Altogether, these data indicate that exo-miRs released in renal IRI are associated with TGF-β signaling. Temporal release of exo-miRs which share targets might be a regulatory mechanism to control the progression of AKI.

Toward Minimal Residual Disease-Directed Therapy in Melanoma

Many patients with advanced cancers achieve dramatic responses to a panoply of therapeutics yet retain minimal residual disease (MRD), which ultimately results in relapse. To gain insights into the biology of MRD, we applied single-cell RNA sequencing to malignant cells isolated from BRAF mutant patient-derived xenograft melanoma cohorts exposed to concurrent RAF/MEK-inhibition. We identified distinct drug-tolerant transcriptional states, varying combinations of which co-occurred within MRDs from PDXs and biopsies of patients on treatment. One of these exhibited a neural crest stem cell (NCSC) transcriptional program largely driven by the nuclear receptor RXRG. An RXR antagonist mitigated accumulation of NCSCs in MRD and delayed the development of resistance. These data identify NCSCs as key drivers of resistance and illustrate the therapeutic potential of MRD-directed therapy. They also highlight how gene regulatory network architecture reprogramming may be therapeutically exploited to limit cellular heterogeneity, a key driver of disease progression and therapy resistance.

A Large Open Pangenome and a Small Core Genome for Giant Pandoraviruses

Giant viruses of amoebae are distinct from classical viruses by the giant size of their virions and genomes. Pandoraviruses are the record holders in size of genomes and number of predicted genes. Three strains, P. salinus, P. dulcis, and P. inopinatum, have been described to date. We isolated three new ones, namely P. massiliensis, P. braziliensis, and P. pampulha, from environmental samples collected in Brazil. We describe here their genomes, the transcriptome and proteome of P. massiliensis, and the pangenome of the group encompassing the six pandoravirus isolates. Genome sequencing was performed with an Illumina MiSeq instrument. Genome annotation was performed using GeneMarkS and Prodigal softwares and comparative genomic analyses. The core genome and pangenome were determined using notably ProteinOrtho and CD-HIT programs. Transcriptomics was performed for P. massiliensis with the Illumina MiSeq instrument; proteomics was also performed for this virus using 1D/2D gel electrophoresis and mass spectrometry on a Synapt G2Si Q-TOF traveling wave mobility spectrometer. The genomes of the three new pandoraviruses are comprised between 1.6 and 1.8 Mbp. The genomes of P. massiliensis, P. pampulha, and P. braziliensis were predicted to harbor 1,414, 2,368, and 2,696 genes, respectively. These genes comprise up to 67% of ORFans. Phylogenomic analyses showed that P. massiliensis and P. braziliensis were more closely related to each other than to the other pandoraviruses. The core genome of pandoraviruses comprises 352 clusters of genes, and the ratio core genome/pangenome is less than 0.05. The extinction curve shows clearly that the pangenome is still open. A quarter of the gene content of P. massiliensis was detected by transcriptomics. In addition, a product for a total of 162 open reading frames were found by proteomic analysis of P. massiliensis virions, including notably the products of 28 ORFans, 99 hypothetical proteins, and 90 core genes. Further analyses should allow to gain a better knowledge and understanding of the evolution and origin of these giant pandoraviruses, and of their relationships with viruses and cellular microorganisms.

Concentration-dependent protein loading of extracellular vesicles released by Histoplasma capsulatum after antibody treatment and its modulatory action upon macrophages

Diverse pathogenic fungi secrete extracellular vesicles (EV) that contain macromolecules, including virulence factors that can modulate the host immune response. We recently demonstrated that the binding of monoclonal antibodies (mAb) modulates how Histoplasma capsulatum load and releases its extracellular vesicles (EV). In the present paper, we addressed a concentration-dependent impact on the fungus' EV loading and release with different mAb, as well as the pathophysiological role of these EV during the host-pathogen interaction. We found that the mAbs differentially regulate EV content in concentration-dependent and independent manners. Enzymatic assays demonstrated that laccase activity in EV from H. capsulatum opsonized with 6B7 was reduced, but urease activity was not altered. The uptake of H. capsulatum by macrophages pre-treated with EV, presented an antibody concentration-dependent phenotype. The intracellular killing of yeast cells was potently inhibited in macrophages pre-treated with EV from 7B6 (non-protective) mAb-opsonized H. capsulatum and this inhibition was associated with a decrease in the reactive-oxygen species generated by these macrophages. In summary, our findings show that opsonization quantitatively and qualitatively modifies H. capsulatum EV load and secretion leading to distinct effects on the host's immune effector mechanisms, supporting the hypothesis that EV sorting and secretion are dynamic mechanisms for a fine-tuned response by fungal cells.

Streamlined Tandem Mass Tag (SL-TMT) Protocol: An Efficient Strategy for Quantitative (Phospho)proteome Profiling Using Tandem Mass Tag-Synchronous Precursor Selection-MS3

Mass spectrometry (MS) coupled toisobaric labeling has developed rapidly into a powerful strategy for high-throughput protein quantification. Sample multiplexing and exceptional sensitivity allow for the quantification of tens of thousands of peptides and, by inference, thousands of proteins from multiple samples in a single MS experiment. Accurate quantification demands a consistent and robust sample-preparation strategy. Here, we present a detailed workflow for SPS-MS3-based quantitative abundance profiling of tandem mass tag (TMT)-labeled proteins and phosphopeptides that we have named the streamlined (SL)-TMT protocol. We describe a universally applicable strategy that requires minimal individual sample processing and permits the seamless addition of a phosphopeptide enrichment step ("mini-phos") with little deviation from the deep proteome analysis. To showcase our workflow, we profile the proteome of wild-type Saccharomyces cerevisiae yeast grown with either glucose or pyruvate as the carbon source. Here, we have established a streamlined TMT protocol that enables deep proteome and medium-scale phosphoproteome analysis.

Nucleosomes of polyploid trophoblast giant cells mostly consist of histone variants and form a loose chromatin structure

Trophoblast giant cells (TGCs) are one of the cell types that form the placenta and play multiple essential roles in maintaining pregnancy in rodents. TGCs have large, polyploid nuclei resulting from endoreduplication. While previous studies have shown distinct gene expression profiles of TGCs, their chromatin structure remains largely unknown. An appropriate combination of canonical and non-canonical histones, also known as histone variants, allows each cell to exert its cell type-specific functions. Here, we aimed to reveal the dynamics of histone usage and chromatin structure during the differentiation of trophoblast stem cells (TSCs) into TGCs. Although the expression of most genes encoding canonical histones was downregulated, the expression of a few genes encoding histone variants such as H2AX, H2AZ, and H3.3 was maintained at a relatively high level in TGCs. Both the micrococcal nuclease digestion assay and nucleosome stability assay using a microfluidic device indicated that chromatin became increasingly loose as TSCs differentiated. Combinatorial experiments involving H3.3-knockdown and -overexpression demonstrated that variant H3.3 resulted in the formation of loose nucleosomes in TGCs. In conclusion, our study revealed that TGCs possessed loose nucleosomes owing to alterations in their histone composition during differentiation.

6.0 K microarray reveals differential transcriptomic responses in the dinoflagellate Prorocentrum minimum exposed to polychlorinated biphenyl (PCB)

Endocrine disrupting chemicals (EDCs) have toxic effects on algae; however, their molecular genomic responses have not been sufficiently elucidated. Here, we evaluated genome-scaled responses of the dinoflagellate alga Prorocentrum minimum exposed to an EDC, polychlorinated biphenyl (PCB), using a 6.0 K microarray. Based on two-fold change cut-off, we identified that 609 genes (∼10.2%) responded to the PCB treatment. KEGG pathway analysis showed that differentially expressed genes (DEGs) were related to ribosomes, biosynthesis of amino acids, spliceosomes, and cellular processes. Many DEGs were involved in cell cycle progression, apoptosis, signal transduction, ion binding, and cellular transportation. In contrast, only a few genes related to photosynthesis and oxidative stress were expressed in response to PCB exposure. This was supported by that fact that there were no obvious changes in the photosynthetic efficiency and reactive oxygen species (ROS) production. These results suggest that PCB might not cause chloroplast and oxidative damage, but could lead to cell cycle arrest and apoptosis. In addition, various signal transduction and transport pathways might be disrupted in the cells, which could further contribute to cell death. These results expand the genomic understanding of the effects of EDCs on this dinoflagellate protist.

Meta-expression analysis of unannotated genes in rice and approaches for network construction to suggest the probable roles

This work suggests 2020 potential candidates in rice for the functional annotation of unannotated genes using meta-analysis of anatomical samples derived from microarray and RNA-seq technologies and this information will be useful to identify novel morphological agronomic traits. Although the genome of rice (Oryza sativa) has been sequenced, 14,365 genes are considered unannotated because they lack putative annotation information. According to the Rice Genome Annotation Project Database ( http://rice.plantbiology.msu.edu/ ), the proportion of functionally characterized unannotated genes (0.35%) is quite limited when compared with the approximately 3.9% of annotated genes with assigned putative functions. Researchers require additional information to help them investigate the molecular mechanisms associated with those unannotated genes. To determine which of them might regulate morphological or physiological traits in the rice genome, we conducted a meta-analysis of expression data that covered a wide range of tissue/organ samples. Overall, 2020 genes showed cultivar-, tissue-, or organ-preferential patterns of expression. Representative candidates from featured groups were validated by RT-PCR, and the GUS reporter system was used to validate the expression of genes that were clustered according to their leaf or root preference. Taking a molecular and genetics approach, we examined meta-expression data and found that 127 genes were differentially expressed between japonica and indica rice cultivars. This is potentially significant for future agronomic applications. We also used a T-DNA insertional mutant and performed a co-expression network analysis of Sword shape dwarf1 (SSD1), a gene that regulates cell division. This network was refined via RT-PCR analysis. Our results suggested that SSD1 represses the expression of four genes related to the processes of DNA replication or cell division and provides insight into possible molecular mechanisms. Together, these strategies present a valuable tool for in-depth characterization of currently unannotated genes.

Gain-of-function miRNA signature by mutant p53 associates with poor cancer outcome

Missense mutation of p53 not only impairs its tumor suppression function, but also causes oncogenic gain of function (GOF). The molecular underpinning of mutant p53 (mutp53) GOF is not fully understood, especially for the potential roles of non-coding genes. Here we identify the microRNA expression profile (microRNAome) of mutp53 on Arg282 by controlled microarray experiments, and clarify the prognostic significance of mutp53-regulated miRNAs in cancers. A predominant repression effect on miRNA expression was found for mutant p53, with 183 significantly downregulated and only 12 upregulated miRNAs. Mutp53 and wild-type (wtp53) commonly upregulate let-7i, and other two miRNAs were upregulated by wtp53 but repressed by mutp53 (miR-610 and miR-3065–3p). Based the mutp53-regulated miRNA signature, a non-negative matrix factorization (NMF) model classified gastric cancer (GC) cases into subgroups with significantly different Disease-free survival (Kaplan-Meier test, P = 0.013). In contrast, the NMF model based on all miRNAs did not associate with cancer outcome. The mutp53 miRNA signature associated with the outcomes of breast cancer (P = 0.024) and hepatocellular cancer (P = 0.012). The miRPath analysis revealed that mutp53-suppressed miRNAs associate with Hippo, TGF-β and stem cell signaling pathways. Taken together, our results highlight a miRNA-mediated GOF mechanism of mutant p53 on Arg282, and suggest the prognostic potential of mutp53-associated miRNA signature.

Quantitative trait locus mapping of Populus bark features and stem diameter

BACKGROUND

Bark plays important roles in photosynthate transport and storage, along with physical and chemical protection. Bark texture varies extensively among species, from smooth to fissured to deeply furrowed, but its genetic control is unknown. This study sought to determine the main genomic regions associated with natural variation in bark features and stem diameter. Quantitative trait loci (QTL) were mapped using an interspecific pseudo-backcross pedigree (Populus trichocarpa x P. deltoides and P. deltoides) for bark texture, bark thickness and diameter collected across three years, two sites and three biological replicates per site.

RESULTS

QTL specific to bark texture were highly reproducible in shared intervals across sites, years and replicates. Significant positive correlations and co-localization between trait QTL suggest pleiotropic regulators or closely linked genes. A list of candidate genes with related putative function, location close to QTL maxima and with the highest expression level in the phloem, xylem and cambium was identified.

CONCLUSION

Candidate genes for bark texture included an ortholog of Arabidopsis ANAC104 (PopNAC128), which plays a role in lignified fiber cell and ray development, as well as Pinin and Fasciclin (PopFLA) genes with a role in cell adhesion, cell shape and migration. The results presented in this study provide a basis for future genomic characterization of genes found within the QTL for bark texture, bark thickness and diameter in order to better understand stem and bark development in Populus and other woody perennial plants. The QTL mapping approach identified a list of prime candidate genes for further validation using functional genomics or forward genetics approaches.

Mouse Cutaneous Melanoma Induced by Mutant BRaf Arises from Expansion and Dedifferentiation of Mature Pigmented Melanocytes

To identify the cells at the origin of melanoma, we combined single-cell lineage-tracing and transcriptomics approaches with time-lapse imaging. A mouse model that recapitulates key histopathological features of human melanomagenesis was created by inducing a BRafV600E-driven melanomagenic program in tail interfollicular melanocytes. Most targeted mature, melanin-producing melanocytes expanded clonally within the epidermis before losing their differentiated features through transcriptional reprogramming and eventually invading the dermis. Tumors did not form within interscales, which contain both mature and dormant amelanotic melanocytes. The hair follicle bulge, which contains melanocyte stem cells, was also refractory to melanomagenesis. These studies identify varying tumor susceptibilities within the melanocytic lineage, highlighting pigment-producing cells as the melanoma cell of origin, and indicate that regional variation in tumor predisposition is dictated by microenvironmental cues rather than intrinsic differences in cellular origin. Critically, this work provides in vivo evidence that differentiated somatic cells can be reprogrammed into cancer initiating cells.

H2A O-GlcNAcylation at serine 40 functions genomic protection in association with acetylated H2AZ or γH2AX

Background

We have previously reported a novel O-GlcNAc modification at serine 40 (S40) of H2A (H2AS40Gc). S40-type H2A isoforms susceptible to O-GlcNAcylation are evolutionarily new and restricted to the viviparous animals; however, the biological function of H2AS40Gc is largely unknown. H2A isoforms are consisted of S40 and alanine 40 (A40) type and this residue on H2A is located in the L1 of the globular domain, which is also known as a variable portion that distinguishes between the canonical and non-canonical H2A variants. In this study, by considering the similarity between the S40-type H2A and histone H2A variants, we explored the function of H2AS40Gc in mouse embryonic stem cells (mESCs).

Results

We found several similarities between the S40-type H2A isoforms and histone H2A variants such H2AZ and H2AX. mRNA of S40-type H2A isoforms (H2A1 N and H2A3) had a poly(A) tail and was produced throughout the cell cycle in contrast to that of A40-type. Importantly, H2AS40Gc level increased owing to chemical-induced DNA damage, similar to phosphorylated H2AX (γH2AX) and acetylated H2AZ (AcH2AZ). H2AS40Gc was accumulated at the restricted area (± 1.5 kb) of DNA damage sites induced by CRISPR/CAS9 system in contrast to accumulation of γH2AX, which was widely scattered. Overexpression of the wild-type (WT) H2A3, but not the S40 to A40 mutation (S40A-mutant), protected the mESC genome against chemical-induced DNA damage. Furthermore, 3 h after the DNA damage treatment, the genome was almost recovered in WT mESCs, whereas the damage advanced further in the S40A-mutant mESCs, suggesting functions of H2AS40Gc in the DNA repair mechanism. Furthermore, the S40A mutant prevented the accumulation of the DNA repair apparatus such as DNA-PKcs and Rad51 at the damage site. Co-immunoprecipitation experiment in WT and S40A-mutant mESCs revealed that H2AS40Gc physiologically bound to AcH2AZ at the initial phase upon DNA damage, followed by binding with γH2AX during the DNA damage repair process.

Conclusions

These data suggest that H2AS40Gc functions to maintain genome integrity through the DNA repair mechanism in association with AcH2AZ and γH2AX.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0157-x) contains supplementary material, which is available to authorized users.

Genome-wide analysis of basic helix-loop-helix (bHLH) transcription factors in Brachypodium distachyon

BACKGROUND

As a superfamily of transcription factors (TFs), the basic helix-loop-helix (bHLH) proteins have been characterized functionally in many plants with a vital role in the regulation of diverse biological processes including growth, development, response to various stresses, and so on. However, no systemic analysis of the bHLH TFs has been reported in Brachypodium distachyon, an emerging model plant in Poaceae.

RESULTS

A total of 146 bHLH TFs were identified in the Brachypodium distachyon genome and classified into 24 subfamilies. BdbHLHs in the same subfamily share similar protein motifs and gene structures. Gene duplication events showed a close relationship to rice, maize and sorghum, and segment duplications might play a key role in the expansion of this gene family. The amino acid sequence of the bHLH domains were quite conservative, especially Leu-27 and Leu-54. Based on the predicted binding activities, the BdbHLHs were divided into DNA binding and non-DNA binding types. According to the gene ontology (GO) analysis, BdbHLHs were speculated to function in homodimer or heterodimer manner. By integrating the available high throughput data in public database and results of quantitative RT-PCR, we found the expression profiles of BdbHLHs were different, implying their differentiated functions.

CONCLUSION

One hundred fourty-six BdbHLHs were identified and their conserved domains, sequence features, phylogenetic relationship, chromosomal distribution, GO annotations, gene structures, gene duplication and expression profiles were investigated. Our findings lay a foundation for further evolutionary and functional elucidation of BdbHLH genes.

Preeclampsia: novel insights from global RNA profiling of trophoblast subpopulations

BACKGROUND

The maternal signs of preeclampsia, which include the new onset of high blood pressure, can occur because of faulty placentation. We theorized that transcriptomic analyses of trophoblast subpopulations in situ would lend new insights into the role of these cells in preeclampsia pathogenesis.

OBJECTIVE

Our goal was to enrich syncytiotrophoblasts, invasive cytotrophoblasts, or endovascular cytotrophoblasts from the placentas of severe preeclampsia cases. Total RNA was subjected to global transcriptional profiling to identify RNAs that were misexpressed compared with controls.

STUDY DESIGN

This was a cross-sectional analysis of placentas from women who had been diagnosed with severe preeclampsia. Gestational age-matched controls were placentas from women who had a preterm birth with no signs of infection. Laser microdissection enabled enrichment of syncytiotrophoblasts, invasive cytotrophoblasts, or endovascular cytotrophoblasts. After RNA isolation, a microarray approach was used for global transcriptional profiling. Immunolocalization identified changes in messenger RNA expression that carried over to the protein level. Differential expression of non-protein-coding RNAs was confirmed by in situ hybridization. A 2-way analysis of variance of non-coding RNA expression identified particular classes that distinguished trophoblasts in cases vs controls. Cajal body foci were visualized by coilin immunolocalization.

RESULTS

Comparison of the trophoblast subtype data within each group (severe preeclampsia or noninfected preterm birth) identified many highly differentially expressed genes. They included molecules that are known to be expressed by each subpopulation, which is evidence that the method worked. Genes that were expressed differentially between the 2 groups, in a cell-type-specific manner, encoded a combination of molecules that previous studies associated with severe preeclampsia and those that were not known to be dysregulated in this pregnancy complication. Gene ontology analysis of the syncytiotrophoblast data highlighted the dysregulation of immune functions, morphogenesis, transport, and responses to vascular endothelial growth factor and progesterone. The invasive cytotrophoblast data provided evidence of alterations in cellular movement, which is consistent with the shallow invasion often associated with severe preeclampsia. Other dysregulated pathways included immune, lipid, oxygen, and transforming growth factor-beta responses. The data for endovascular cytotrophoblasts showed disordered metabolism, signaling, and vascular development. Additionally, the transcriptional data revealed the differential expression in severe preeclampsia of 2 classes of non-coding RNAs: long non-coding RNAs and small nucleolar RNAs. The long non-coding RNA, urothelial cancer associated 1, was the most highly up-regulated in this class. In situ hybridization confirmed severe preeclampsia-associated expression in syncytiotrophoblasts. The small nucleolar RNAs, which chemically modify RNA structure, also correlated with severe preeclampsia. Thus, we enumerated Cajal body foci, sites of small nucleolar RNA activity, in primary cytotrophoblasts that were isolated from control and severe preeclampsia placentas. In severe preeclampsia, cytotrophoblasts had approximately double the number of these foci as the control samples.

CONCLUSION

A laser microdissection approach enabled the identification of novel messenger RNAs and non-coding RNAs that were misexpressed by various trophoblast subpopulations in severe preeclampsia. The results suggested new avenues of investigation, in particular, the roles of PRG2, Kell blood group determinants, and urothelial cancer associated 1 in syncytiotrophoblast diseases. Additionally, many of the newly identified dysregulated molecules might have clinical utility as biomarkers of severe preeclampsia.

BIG LEAF is a regulator of organ size and adventitious root formation in poplar

Here we report the discovery through activation tagging and subsequent characterization of the BIG LEAF (BL) gene from poplar. In poplar, BL regulates leaf size via positively affecting cell proliferation. Up and downregulation of the gene led to increased and decreased leaf size, respectively, and these phenotypes corresponded to increased and decreased cell numbers. BL function encompasses the early stages of leaf development as native BL expression was specific to the shoot apical meristem and leaf primordia and was absent from the later stages of leaf development and other organs. Consistently, BL downregulation reduced leaf size at the earliest stages of leaf development. Ectopic expression in mature leaves resulted in continued growth most probably via sustained cell proliferation and thus the increased leaf size. In contrast to the positive effect on leaf growth, ectopic BL expression in stems interfered with and significantly reduced stem thickening, suggesting that BL is a highly specific activator of growth. In addition, stem cuttings from BL overexpressing plants developed roots, whereas the wild type was difficult to root, demonstrating that BL is a positive regulator of adventitious rooting. Large transcriptomic changes in plants that overexpressed BL indicated that BL may have a broad integrative role, encompassing many genes linked to organ growth. We conclude that BL plays a fundamental role in control of leaf size and thus may be a useful tool for modifying plant biomass productivity and adventitious rooting.

shinyheatmap: Ultra fast low memory heatmap web interface for big data genomics

BACKGROUND

Transcriptomics, metabolomics, metagenomics, and other various next-generation sequencing (-omics) fields are known for their production of large datasets, especially across single-cell sequencing studies. Visualizing such big data has posed technical challenges in biology, both in terms of available computational resources as well as programming acumen. Since heatmaps are used to depict high-dimensional numerical data as a colored grid of cells, efficiency and speed have often proven to be critical considerations in the process of successfully converting data into graphics. For example, rendering interactive heatmaps from large input datasets (e.g., 100k+ rows) has been computationally infeasible on both desktop computers and web browsers. In addition to memory requirements, programming skills and knowledge have frequently been barriers-to-entry for creating highly customizable heatmaps.

RESULTS

We propose shinyheatmap: an advanced user-friendly heatmap software suite capable of efficiently creating highly customizable static and interactive biological heatmaps in a web browser. shinyheatmap is a low memory footprint program, making it particularly well-suited for the interactive visualization of extremely large datasets that cannot typically be computed in-memory due to size restrictions. Also, shinyheatmap features a built-in high performance web plug-in, fastheatmap, for rapidly plotting interactive heatmaps of datasets as large as 105-107 rows within seconds, effectively shattering previous performance benchmarks of heatmap rendering speed.

CONCLUSIONS

shinyheatmap is hosted online as a freely available web server with an intuitive graphical user interface: http://shinyheatmap.com. The methods are implemented in R, and are available as part of the shinyheatmap project at: https://github.com/Bohdan-Khomtchouk/shinyheatmap. Users can access fastheatmap directly from within the shinyheatmap web interface, and all source code has been made publicly available on Github: https://github.com/Bohdan-Khomtchouk/fastheatmap.

The protein phosphatase 2C clade A protein OsPP2C51 positively regulates seed germination by directly inactivating OsbZIP10

Protein phosphatase 2C clade A members are major signaling components in the ABA-dependent signaling cascade that regulates seed germination. To elucidate the role of PP2CA genes in germination of rice seed, we selected OsPP2C51, which shows highly specific expression in the embryo compared with other protein phosphatases based on microarray data. GUS histochemical assay confirmed that OsPP2C51 is expressed in the seed embryo and that this expression pattern is unique compared with those of other OsPP2CA genes. Data obtained from germination assays and alpha-amylase assays of OsPP2C51 knockout and overexpression lines suggest that OsPP2C51 positively regulates seed germination in rice. The expression of alpha-amylase synthesizing genes was high in OsPP2C51 overexpressing plants, suggesting that elevated levels of OsPP2C51 might enhance gene expression related to higher rates of seed germination. Analysis of protein interactions between ABA signaling components showed that OsPP2C51 interacts with OsPYL/RCAR5 in an ABA-dependent manner. Furthermore, interactions were observed between OsPP2C51 and SAPK2, and between OsPP2C51 and OsbZIP10 and we found out that OsPP2C51 can dephosphorylates OsbZIP10. These findings suggest the existence of a new branch in ABA signaling pathway consisting of OsPYL/RCAR-OsPP2C-bZIP apart from the previously reported OsPYL/RCAR-OsPP2C-SAPK-bZIP. Overall, our result suggests that OsPP2C51 is a positive regulator of seed germination by directly suppressing active phosphorylated OsbZIP10.

Alteration of N-glycan expression profile and glycan pattern of glycoproteins in human hepatoma cells after HCV infection

BACKGROUND

Hepatitis C virus (HCV) infection causes chronic liver diseases, liver fibrosis and even hepatocellular carcinoma (HCC). However little is known about any information of N-glycan pattern in human liver cell after HCV infection.

METHODS

The altered profiles of N-glycans in HCV-infected Huh7.5.1 cell were analyzed by using mass spectrometry. Then, lectin microarray, lectin pull-down assay, reverse transcription-quantitative real time PCR (RT-qPCR) and western-blotting were used to identify the altered N-glycosylated proteins and glycosyltransferases.

RESULTS

Compared to uninfected cells, significantly elevated levels of fucosylated, sialylated and complex N-glycans were found in HCV infected cells. Furthermore, Lens culinaris agglutinin (LCA)-binding glycoconjugates were increased most. Then, the LCA-agarose was used to precipitate the specific glycosylated proteins and identify that fucosylated modified annexin A2 (ANXA2) and heat shock protein 90 beta family member 1 (HSP90B1) was greatly increased in HCV-infected cells. However, the total ANXA2 and HSP90B1 protein levels remained unchanged. Additionally, we screened the mRNA expressions of 47 types of different glycosyltransferases and found that α1,6-fucosyltransferase 8 (FUT8) was the most up-regulated and contributed to strengthen the LCA binding capability to fucosylated modified ANXA2 and HSP90B1 after HCV infection.

CONCLUSIONS

HCV infection caused the altered N-glycans profiles, increased expressions of FUT8, fucosylated ANXA2 and HSP90B1 as well as enhanced LCA binding to Huh7.5.1.

GENERAL SIGNIFICANCE

Our results may lay the foundation for clarifying the role of N-glycans and facilitate the development of novel diagnostic biomarkers and therapeutic targets based on the increased FUT8, fucosylated ANXA2 and HSP90B1 after HCV infection.

Genome-Wide Identification and Analysis of Genes, Conserved between japonica and indica Rice Cultivars, that Respond to Low-Temperature Stress at the Vegetative Growth Stage

Cold stress is very detrimental to crop production. However, only a few genes in rice have been identified with known functions related to cold tolerance. To meet this agronomic challenge more effectively, researchers must take global approaches to select useful candidate genes and find the major regulatory factors. We used five Gene expression omnibus series data series of Affymetrix array data, produced with cold stress-treated samples from the NCBI Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/), and identified 502 cold-inducible genes common to both japonica and indica rice cultivars. From them, we confirmed that the expression of two randomly chosen genes was increased by cold stress in planta. In addition, overexpression of OsWRKY71 enhanced cold tolerance in 'Dongjin,' the tested japonica cultivar. Comparisons between japonica and indica rice, based on calculations of plant survival rates and chlorophyll fluorescence, confirmed that the japonica rice was more cold-tolerant. Gene Ontology enrichment analysis indicate that the 'L-phenylalanine catabolic process,' within the Biological Process category, was the most highly overrepresented under cold-stress conditions, implying its significance in that response in rice. MapMan analysis classified 'Major Metabolic' processes and 'Regulatory Gene Modules' as two other major determinants of the cold-stress response and suggested several key cis-regulatory elements. Based on these results, we proposed a model that includes a pathway for cold stress-responsive signaling. Results from our functional analysis of the main signal transduction and transcription regulation factors identified in that pathway will provide insight into novel regulatory metabolism(s), as well as a foundation by which we can develop crop plants with enhanced cold tolerance.

Integrated and Quantitative Proteomics of Human Tumors

Quantitative proteomics represents a powerful approach for the comprehensive analysis of proteins expressed under defined conditions. These properties have been used to investigate the proteome of disease states, including cancer. It has become a major subject of studies to apply proteomics for biomarker and therapeutic target identification. In the last decades, technical advances in mass spectrometry have increased the capacity of protein identification and quantification. Moreover, the analysis of posttranslational modification (PTM), especially phosphorylation, has allowed large-scale identification of biological mechanisms. Even so, increasing evidence indicates that global protein quantification is often insufficient for the explanation of biology and has shown to pose challenges in identifying new and robust biomarkers. As a consequence, to improve the accuracy of the discoveries made using proteomics in human tumors, it is necessary to combine (i) robust and reproducible methods for sample preparation allowing statistical comparison, (ii) PTM analyses in addition to global proteomics for additional levels of knowledge, and (iii) use of bioinformatics for decrypting protein list. Herein, we present technical specificities for samples preparation involving isobaric tag labeling, TiO₂-based phosphopeptides enrichment and hydrazyde-based glycopeptides purification as well as the key points for the quantitative analysis and interpretation of the protein lists. The method is based on our experience with tumors analysis derived from hepatocellular carcinoma, chondrosarcoma, human embryonic intervertebral disk, and chordoma experiments.

MicroScope: ChIP-seq and RNA-seq software analysis suite for gene expression heatmaps

Background

Heatmaps are an indispensible visualization tool for examining large-scale snapshots of genomic activity across various types of next-generation sequencing datasets. However, traditional heatmap software do not typically offer multi-scale insight across multiple layers of genomic analysis (e.g., differential expression analysis, principal component analysis, gene ontology analysis, and network analysis) or multiple types of next-generation sequencing datasets (e.g., ChIP-seq and RNA-seq). As such, it is natural to want to interact with a heatmap’s contents using an extensive set of integrated analysis tools applicable to a broad array of genomic data types.

Results

We propose a user-friendly ChIP-seq and RNA-seq software suite for the interactive visualization and analysis of genomic data, including integrated features to support differential expression analysis, interactive heatmap production, principal component analysis, gene ontology analysis, and dynamic network analysis.

Conclusions

MicroScope is hosted online as an R Shiny web application based on the D3 JavaScript library: http://microscopebioinformatics.org/. The methods are implemented in R, and are available as part of the MicroScope project at: https://github.com/Bohdan-Khomtchouk/Microscope.