Search

Sorting-related receptor with A-type repeats (SORLA) is an important receptor for regulating normal cellular functions via protein sorting. Here, we determined the structures of the full-length SORLA and identified two distinct conformations of apo-SORLA using single-particle cryogenic electron microscopy. In contrast to homologous proteins, both monomer and dimer forms of SORLA existed in a neutral solution. Only three hydrogen bonds in the vicinity of the dimer interface implied the involvement in dimerization. The orientation of residue R490 was a key point for ligand binding. These results suggest a unique mechanism of SORLA dimerization for protein trafficking.

Bearing is an important part of rotary machinery, which usually operates under a variety of complicated and severe conditions, and is prone to break down. To reduce and prevent the loss caused by the fault of bearings, a method based on multi‐band filtering (MBF) is proposed and applied to bearing fault diagnosis in this paper. Bearing signal is decomposed into multiple sub‐band signals by an MBF constructed from a specific prototype filter (finite impulse response filter) through cosine modulation. Then, the required sub‐band can be selected adaptively according to bearing fault frequency. To make the fault frequency more prominent, an adaptive filtering method is exploited to reduce the noise contained in the selected sub‐band. Finally, bearing fault diagnosis is realized by envelope spectrum analysis. Experimental results show that better performance in both simulated and real data are achieved by the proposed algorithm, which indicates that the proposed method can realize bearing fault diagnosis efficiently when compared to other state of the art methods.

Zika virus (ZIKV) infection during pregnancy can cause microcephaly in newborns, yet the underlying mechanisms remain largely unexplored. Here, we reveal extensive and large-scale metabolic reprogramming events in ZIKV-infected mouse brains by performing a multi-omics study comprising transcriptomics, proteomics, phosphoproteomics and metabolomics approaches. Our proteomics and metabolomics analyses uncover dramatic alteration of nicotinamide adenine dinucleotide (NAD+)-related metabolic pathways, including oxidative phosphorylation, TCA cycle and tryptophan metabolism. Phosphoproteomics analysis indicates that MAPK and cyclic GMP-protein kinase G signaling may be associated with ZIKV-induced microcephaly. Notably, we demonstrate the utility of our rich multi-omics datasets with follow-up in vivo experiments, which confirm that boosting NAD+ by NAD+ or nicotinamide riboside supplementation alleviates cell death and increases cortex thickness in ZIKV-infected mouse brains. Nicotinamide riboside supplementation increases the brain and body weight as well as improves the survival in ZIKV-infected mice. Our study provides a comprehensive resource of biological data to support future investigations of ZIKV-induced microcephaly and demonstrates that metabolic alterations can be potentially exploited for developing therapeutic strategies.

Although institutional pressures have huge strategic implications for organizational activities, this certainly does not mean that organizations under institutional pressures can improve environmental performance automatically. Institutional pressures are critical but not sufficient to affect environmental performance directly. Therefore, additional research is needed to explore the missing link between institutional pressures and environmental performance. Based on the “pressure-response-performance” framework, this study integrates perspectives of institutional theory and organizational learning to argue the mediating role of organizational learning in the relationship between institutional pressures and environmental performance. Data were collected via 268 valid questionnaires from construction firms located in Shanxi Province in central China. Hypotheses in the conceptual model were tested with structural equation modeling. Empirical results reveal that both coercive and mimetic pressures have significantly positive effects on organizational learning, whereas normative pressures have a non-significant effect on organizational learning. Besides that, organizational learning has a significantly positive effect on environmental performance. In addition, organizational learning partially mediates the relationship between coercive pressures and environmental performance and completely mediates the relationship between mimetic pressures and environmental performance. By exploring the mediating role of organizational learning, the article uncovers the missing link in the relationship between institutional pressures and environmental performance.

Sonic hedgehog (Shh) signaling plays a critical role in regulating cerebellum development by maintaining the physiological proliferation of granule neuron precursors (GNPs), and its dysregulation leads to the oncogenesis of medulloblastoma. O-GlcNAcylation (O-GlcNAc) of proteins is an emerging regulator of brain function that maintains normal development and neuronal circuitry. Here, we demonstrate that O-GlcNAc transferase (OGT) in GNPs mediate the cerebellum development, and the progression of the Shh subgroup of medulloblastoma. Specifically, OGT regulates the neurogenesis of GNPs by activating the Shh signaling pathway via O-GlcNAcylation at S355 of GLI family zinc finger 2 (Gli2), which in turn promotes its deacetylation and transcriptional activity via dissociation from p300, a histone acetyltransferases. Inhibition of OGT via genetic ablation or chemical inhibition improves survival in a medulloblastoma mouse model. These data uncover a critical role for O-GlcNAc signaling in cerebellar development, and pinpoint a potential therapeutic target for Shh-associated medulloblastoma.

There have been large foodborne outbreaks related to Enterohemorrhagic

Summary Cellular responses to transforming growth factor β (TGF-β) depend on cell context. Here, we explored how TGF-β/nodal signaling crosstalks with the epigenome to promote mesendodermal differentiation. We find that expression of a group of mesendodermal genes depends on both TRIM33 and nodal signaling in embryoid bodies (EBs) but not in embryonic stem cells (ESCs). Only in EBs, TRIM33 binds these genes in the presence of expanded H3K18ac marks. Furthermore, the H3K18ac landscape at mesendodermal genes promotes TRIM33 recruitment. We reveal that HDAC1 binds to active gene promoters and interferes with TRIM33 recruitment to mesendodermal gene promoters. However, the TRIM33-interacting protein p300 deposits H3K18ac and further enhances TRIM33 recruitment. ATAC-seq data demonstrate that TRIM33 primes mesendodermal genes for activation by maintaining chromatin accessibility at their regulatory regions. Altogether, our study suggests that HDAC1 and p300 are key factors linking the epigenome through TRIM33 to the cell context-dependent nodal response during mesendodermal differentiation., Highlights • The H3K18ac landscape changes at gene promoters during mesendodermal differentiation • Histone acetylation facilitates TRIM33 recruitment to the chromatin • HDAC1 and p300 are the key factors for nodal signaling crosstalk with epigenome • Chromatin accessibility at mesendodermal genes depends on TRIM33, Xi and colleagues explored how nodal signaling crosstalks with epigenome to specifically promote mesendodermal differentiation. They showed that TRIM33 specifically regulates the expression of a group of mesendodermal genes upon nodal signaling by maintaining the H3K18ac mark, and that HDAC1 and p300 are the key factors linking epigenome to the cell context-dependent nature of nodal responses during mesendodermal differentiation.

The extensive collection of glossy (gl) and eceriferum (cer) mutants of maize and Arabidopsis have proven invaluable in dissecting the branched metabolic pathways that support cuticular lipid deposition. This bifurcated pathway integrates a fatty acid elongation‐decarbonylative branch and a fatty acid elongation‐reductive branch, which collectively has the capacity to generate hundreds of cuticular lipid metabolites. In this study, a combined transgenic and biochemical strategy was implemented to explore and compare the physiological function of three homologous genes, Gl2, Gl2‐like, and CER2, in the context of this branched pathway. These biochemical characterizations integrated new extraction chromatographic procedures with high spatial resolution mass spectrometric imaging methods to profile the cuticular lipids on developing floral tissues transgenically expressing these transgenes in wild‐type or cer2 mutant lines of Arabidopsis. Collectively, these datasets establish that both the maize Gl2 and Gl2‐like genes are functional homologs of the Arabidopsis CER2 gene. In addition, the dynamic distribution of cuticular lipid deposition follows distinct floral organ localization patterns indicating that the fatty acid elongation‐decarbonylative branch of the pathway is differentially localized from the fatty acid elongation‐reductive branch of the pathway.

EcoHIV/NDK is a recombinant virus modified from HIV-1 that uses mice as a host. NanoLuc (Nluc) luciferase is a genetically engineered small molecule enzyme (19.1kDa). It is currently the best bioluminescence reporter gene which uses furimazine as a luminescent substrate to produce detectable fluorescence. In order to quantify and locate the levels of EcoHIV/NDK infection and expression, the EcoHIV/NDK plasmid and the Nluc luciferase gene were integrated together to construct EcoHIV/NDK-NanoLuc virus. To verify whether EcoHIV/NDK-NanoLuc is effective in in vivo and in vitro, we conducted the following experiments. First, the existence of the nanoluc fragment in the recombinant plasmid was verified by PCR amplification, and then the correctness of the plasmid construction was verified by restriction enzyme digestion. As a result, the correct recombinant and plasmid were constructed. The Nluc luciferase reporter gene was used to indicate the infection of the successfully packaged virus on Hela-mCAT cells and the three strains of Kunming mice, BALB/c nude mice, and 129S2/SvPasCrl mice, respectively. Our results show that EcoHIV/NDK-NanoLuc can effectively infect Hela-mCAT cells and three strains of mice, in which viral replication site and distribution was detected through the luciferase reporter gene.

Currently, the laboratory diagnostic tests available for HIV-1 viral infection are mainly based on serological testing which relies on enzyme-linked immunosorbent assay (ELISA) for blood HIV antigen detection and reverse transcription polymerase chain reaction (RT-PCR) for HIV specific RNA sequence identification. However, these methods are expensive and time-consuming, and suffer from false positive and/or false negative results. Thus, there is an urgent need for developing a cost effective, rapid and accurate diagnostic method for HIV-1 infection. In order to reduce the barriers for effective diagnosis, a near-infrared spectroscopy (NIR) method was used to detect the HIV-1 virus in human serum, specifically, three absorption peaks with dose-dependent at 1582nm, 1810nm and 2363nm were found by multiple FBiPLSR test analysis for HIV-nano and HIV-EGFP, but not for MLV. Therefore, we recommend the use of 1582nm, 1810nm and 2363nm as the characteristic spectrum peak, for early screening and rapid diagnosis of serum HIV.

The advanced aero-engine rotor is composed of different materials and connected by the joints such as bolted structure and spline, which shows discontinuity. The rotor always operates at high speed, and its bending deformation changes the contact status of joint interface, which leads to the stiffness loss and the change of the rotor dynamics. For the high-speed rotor system with interface, taking the spline structure as an example, the bend stiffness loss model is established, and the cause of bend stiffness loss at the joint is explained. On the basis of theory, the stiffness loss correction coefficient is introduced, and the influence of stiffness loss on the critical speed and bearing load are analyzed. In order to optimize the rotor dynamic performance, the robust design methods for modal properties and dynamic response behavior of rotor systems are developed with strain energy distribution on the rotor, aiming to reduce the stiffness loss and the sensitivity of the rotor dynamics. A numerical example is presented, the results show that optimizing the bearing support position and stiffness contributes to the reduction of the strain energy ratio on the joint and the sensitivity of the rotor dynamics, thus optimizing the rotor dynamic performance and its robustness.

The extensive collection of glossy (gl) and eceriferum (cer) mutants of maize and Arabidopsis have proven invaluable in dissecting the branched metabolic pathways that support cuticular lipid deposition. This branched pathway integrates the fatty acid elongation-decarbonylative branch and the fatty acid elongation-reductive branch that has the capacity to generate hundreds of cuticular lipid metabolites. In this study a combined transgenic and biochemical strategy was implemented to explore and compare the physiological function of three homologous genes, Gl2, Gl2-like and CER2 in the context of this branched pathway. These biochemical characterizations integrated new extraction-chromatographic procedures with high-spatial resolution mass spectrometric imaging methods to profile the cuticular lipids on developing floral tissues transgenically expressing these transgenes in wild-type or cer2 mutant lines of Arabidopsis. Collectively, these datasets establish that both the maize Gl2 and Gl2-like genes are functional homologs of the Arabidopsis CER2 gene. In addition, the dynamic distribution of cuticular lipid deposition follows distinct floral organ localization patterns indicating that the fatty acid elongation-decarbonylative branch of the pathway is differentially localized from the fatty acid elongation-reductive branch of the pathway.

Zika virus (ZIKV) infection during pregnancy can cause microcephaly in newborns, yet the underlying mechanisms remain largely unexplored. Here, we reveal extensive and large-scale metabolic reprogramming events in ZIKV-infected mouse brains by performing a multi-omics study comprising transcriptomics, proteomics, phosphoproteomics and metabolomics approaches. Our proteomics and metabolomics analyses uncover dramatic alteration of nicotinamide adenine dinucleotide (NAD

Background Poly C Binding Protein 1 (PCBP1) is an RNA-binding protein that binds and regulates translational activity of subsets of cellular mRNAs. Depletion of PCBP1 is implicated in various carcinomas, but the underlying mechanism in tumorigenesis remains elusive. Methods We performed a transcriptome-wide screen to identify novel bounding mRNA of PCBP1. The bind regions between PCBP1 with target mRNA were investigated by using point mutation and luciferase assay. Cell proliferation, cell cycle, tumorigenesis and cell apoptosis were also evaluated in ovary and colon cancer cell lines. The mechanism that PCBP1 affects p27 was analyzed by mRNA stability and ribosome profiling assays. We analyzed PCBP1 and p27 expression in ovary, colon and renal tumor samples and adjacent non-tumor tissues using RT-PCR, Western Blotting and immunohistochemistry. The prognostic significance of PCBP1 and p27 also analyzed using online databases. Results We identified cell cycle inhibitor p27Kip1 (p27) as a novel PCBP1-bound transcript. We then demonstrated that binding of PCBP1 to p27 3’UTR via its KH1 domain mainly stabilizes p27 mRNA, while enhances its translation to fuel p27 expression, prior to p27 protein degradation. The upregulated p27 consequently inhibits cell proliferation, cell cycle progression and tumorigenesis, whereas promotes cell apoptosis under paclitaxel treatment. Conversely, knockdown of PCBP1 in turn compromises p27 mRNA stability, leading to lower p27 level and tumorigenesis in vivo. Moreover, forced depletion of p27 counteracts the tumor suppressive ability of PCBP1 in the same PCBP1 over-expressing cells. Physiologically, we showed that decreases of both p27 mRNA and its protein expressions are well correlated to PCBP1 depletion in ovary, colon and renal tumor samples, independent of the p27 ubiquitin ligase Skp2 level. Correlation of PCBP1 with p27 is also found in the tamoxifen, doxorubincin and lapatinib resistant breast cancer cells of GEO database. Conclusion Our results thereby indicate that loss of PCBP1 expression firstly attenuates p27 expression at post-transcriptional level, and subsequently promotes carcinogenesis. PCBP1 could be used as a diagnostic marker to cancer patients. Electronic supplementary material The online version of this article (10.1186/s13046-018-0840-1) contains supplementary material, which is available to authorized users.

Rapid urbanization has great potential to adversely impact the local hydrological cycle, the environment, and ecosystems, as well as trigger severe problems, such as urban flooding, waterlogging, and water contamination. Understanding runoff generation mechanisms in urban areas helps identify the impacts of urbanization on runoff response and better simulate urban floods. However, in most urban hydrological models, urban surfaces are coarsely classified into impervious and pervious ones, which likely ignores differential rainfall-runoff responses of different urban surfaces and leads to large discrepancies between simulated and observed runoff. Here we developed the TVGM_Urban model, a new urban hydrological model based on the time variant gain model (TVGM), that represents nonlinear rainfall-runoff relationships for different urban surfaces. We applied the model to the Fenghuangcheng region of Shenzhen City, China and compared the results to those of two commonly used urban hydrological models: the Horton infiltration and Soil Conservation Service Curve Number (SCSCN) model. For each of these models, we conducted model uncertainty analysis using the GLUE method. The results showed that the TVGM_Urban model outperformed the other two models in terms of total runoff and peak flow. This could be due to the consideration of different land covers and both saturation excess and infiltration excess runoff generation in the TVGM_Urban model. In addition, the uncertainty analysis indicated better performance of the TVGM_Urban model in reducing structural uncertainty and prediction uncertainty. This study highlights the need to account for detailed land covers and different runoff generation mechanisms in urban hydrological modeling.

Understanding drought dynamics spatially and temporally is essential for environmental and socio-economic systems, especially under climate change with expected increasing drought. This study aims to update the knowledge about drought spatiotemporal changes over mainland China in the 21st century by utilizing the latest Coupled Model Intercomparison Project Phase 6 (CMIP6) datasets. Precipitation from the nine CMIP6 models is first bias-corrected using a quantile mapping approach and then used to calculate the 3-month Standardized Precipitation Index (SPI) for the historical period and future three scenarios (SSP1-2.6, SSP2-4.5 and SSP5-8.5). We conduct a clustering algorithm to identify three-dimensional (i.e., latitude, longitude, and time) drought events to access the drought dynamics in both space and time. Results show that the drought events occur more frequently during June to September than the other months for both historical and future periods. The drought frequency is significantly increasing during Sep–Dec in the 21st century (increase by about 8.5%, p-value

Germination is a highly complex process by which seeds begin to develop and establish themselves as viable organisms. In this study, we utilize a combination of gas chromatography-mass spectrometry, liquid chromatography-fluorescence, and mass spectrometry imaging approaches to profile and visualize the metabolic distributions of germinating seeds from two different inbreds of maize (Zea mays) seeds, B73 and Mo17. Gas chromatography and liquid chromatography analyses demonstrate that the two inbreds are highly differentiated in their metabolite profiles throughout the course of germination, especially with regard to amino acids, sugar alcohols, and small organic acids. Crude dissection of the seed followed by gas chromatography-mass spectrometry analysis of polar metabolites also revealed that many compounds were highly sequestered among the various seed tissue types. To further localize compounds, matrix-assisted laser desorption/ionization mass spectrometry imaging was utilized to visualize compounds in fine detail in their native environments over the course of germination. Most notably, the fatty acyl chain-dependent differential localization of phospholipids and triacylglycerols was observed within the embryo and radicle, showing correlation with the heterogeneous distribution of fatty acids. Other interesting observations include unusual localization of ceramides on the endosperm/scutellum boundary and subcellular localization of ferulate in the aleurone.

The authors describe a quartz crystal microbalance based aptasensor for the determination of Pb2+. In order to enhance its response, oligonucleotide-embellished gold nanoparticles (AuNPs) were used to amplify the frequency changes. The method is based on the use of specific aptamers immobilized on the surface of the quartz crystal microbalance (QCM) and the binding of Pb2+, which prevents the self-assembly of the AuNPs on the QCM. Trace concentrations of Pb2+ can be determined by monitoring the change in the response frequency of the QCM. The method has a 4 nmol∙L−1 detection limit and works in the 5 to 200 nmol∙L−1 Pb2+ concentration range. This aptasensor also shows adequate selectivity for Pb2+ over potentially interfering metal ions.

Given the growing need for biofuel production but the lack of suitable land for producing biomass feedstock, development of stress-tolerant energy crops will be increasingly important. We used comparative transcriptomics to reveal differential responses to long-term salt stress among five populations of Miscanthus lutarioriparius grown in the natural habitats and salinity experimental site. A total of 59 genes were found to be potentially responsive to the high salinity conditions shared by the five populations, including those involved in detoxification, plant defense, photosynthesis, and signal transduction. Of these genes, about 70% were related to abiotic stress response. Among five populations, the most contrasting performance between relatively high survival rates and the relatively weak growing traits were in accordance with the down regulation of genes involved in growth and upregulation of genes related to plant stress tolerance in one of the populations. These results might reveal a potential tolerance-productivity tradeoff, where resources were allocated from growth to stress resistance. The comparative transcriptomics of different populations among different environments will provide a basis for breeding and domestication of energy crops. This article is protected by copyright. All rights reserved.

The regional frequency analysis (RFA) is a widely used method in analyzing the changes of extreme precipitation (EP). The uncertainties in the identification of homogeneous subregions and the selection of optimal regional frequency distributions can largely influence the results of the RFA. In this study, the fuzzy c-means method combined with the extended Xie-Benn index (FCXB) is used to help determine the optimal division of subregions. In addition, we introduce a new comprehensive index (CI) to overcome the shortcomings of present measures and reduce the uncertainty in regional frequency distribution selection. The changes of EP at 93 meteorological stations in the Songliao River Basin (SRB) during 1960 to 2016 is analyzed. The main results show that: (1) FCXB can effectively identify the optimum number of homogeneous subregions automatically, and the corresponding subregion division is proven to be reasonable and reliable; (2) compared with the single goodness-of-fit measure, the developed CI can reduce the uncertainties in distribution selection and determine the optimal regional distribution in a reliable way; (3) the estimated EP under different return periods both decrease from the south to the north of the SRB, which indicates the risk of high-intensity EP events in the southern SRB is relatively higher. These findings can provide technical support for local policymakers to formulate effective measures to lessen the damages of the EP on ecosystem and society.

Highlight We uncovered significant structural and compositional differences between leaf cuticles of drought-tolerant and drought-sensitive wheat cultivars. Specific MYB factors drive cuticular responses to drought and may underlie genotypes with contrasting drought tolerance., A plant cuticle forms a hydrophobic layer covering plant organs, and plays an important role in plant development and protection from environmental stresses. We examined epicuticular structure, composition, and a MYB-based regulatory network in two Australian wheat cultivars, RAC875 and Kukri, with contrasting cuticle appearance (glaucousness) and drought tolerance. Metabolomics and microscopic analyses of epicuticular waxes revealed that the content of β-diketones was the major compositional and structural difference between RAC875 and Kukri. The content of β-diketones remained the same while those of alkanes and primary alcohols were increased by drought in both cultivars, suggesting that the interplay of all components rather than a single one defines the difference in drought tolerance between cultivars. Six wheat genes encoding MYB transcription factors (TFs) were cloned; four of them were regulated in flag leaves of both cultivars by rapid dehydration and/or slowly developing cyclic drought. The involvement of selected MYB TFs in the regulation of cuticle biosynthesis was confirmed by a transient expression assay in wheat cell culture, using the promoters of wheat genes encoding cuticle biosynthesis-related enzymes and the SHINE1 (SHN1) TF. Two functional MYB-responsive elements, specifically recognized by TaMYB74 but not by other MYB TFs, were localized in the TdSHN1 promoter. Protein structural determinants underlying the binding specificity of TaMYB74 for functional DNA cis-elements were defined, using 3D protein molecular modelling. A scheme, linking drought-induced expression of the investigated TFs with downstream genes that participate in the synthesis of cuticle components, is proposed.

The widely used probability based drought indices for drought characterization usually underlie the assumption of the considered hydrometeorological variables to be stationary, which is largely challenged under the current global change. Incorporating the parameters of probability distribution with covariates, such as time or climate indices, has become a common way to develop drought index considering nonstationarity. However most previous drought studies emphasized more on the linear relationship between distribution parameters and related covariates, which may neglect the inherent nonlinear relationship and lost useful information for drought assessment. Our study aims to develop a Nonstationary Standardized Precipitation Index (NSPI) with distribution parameters nonlinearly varying with potential influencing climate indices. The middle and lower reaches of the Yangtze River Basin (MLRYRB) of China are selected as the case study area to examine the performance of the NSPI with the performance of the traditional standardized precipitation index (SPI). Results indicate that the NSPI (with climate indices as covariates) is more robust than both the NSPI (with time as the covariate) and the traditional SPI. Furthermore, the nonlinear dependence of distribution parameters on climate indices considered in the nonstationary model can be more suitable for drought assessment than the linear dependence of distribution parameters. Besides, there exists a higher frequency of extreme drought by the NSPI than the SPI in recent decades in most regions of the MLRYRB. Although the NSPI is a little computationally expensive with more required inputs of climate variables, it can depict the influence of changing environment on the drought occurrence, thus can be a feasible alternative for drought assessment considering nonstationarity in view of future change and provide valuable support for further studies.

The authors describe a piezoelectric aptasensor for arsenite. A self assembeled monolayer (SAM) of mercaptoethylamine was prepared to immobilize arsenite on the surface of a quartz crystal microbalance. Gold nanoparticles were modified with arsenite aptamer to amplify the response frequency of the biosensor. Arsenite first binds to the SAM on the gold surface of the QCM. On addition of gold nanoparticles with aptamer (DNA-AuNp), the SAM-As(III)-aptamer sandwich is formed. This increases the resonance frequency of the sensor and allows trace concentration of arsenite to be determined. The aptasensor can detect arsenite in the 8 to 1000 nmol·L−1 concentration range with a 4.4 nmol·L−1 lower detection limit (at S/N = 3). The sandwich structure improves the specificity of the aptasensor without considering the conformational transition of the aptamer. The strategy described here conceivably has a large potential as it shows that small molecules can be sensed by using aptamers with unknown working mechanism.

Coexpression network revealing genes with Co-variation Expression pattern (CE) and those with Top rank of Expression fold change (TE) played different roles in responding to new environment of Miscanthus lutarioriparius. Variation in gene expression level, the product of genetic and/or environmental perturbation, determines the robustness-to-plasticity spectrum of a phenotype in plants. Understanding how expression variation of plant population response to a new field is crucial to domesticate energy crops. Weighted Gene Coexpression Network Analysis (WGCNA) was used to explore the patterns of expression variation based on 72 Miscanthus lutarioriparius transcriptomes from two contrasting environments, one near the native habitat and the other in one harsh domesticating region. The 932 genes with Co-variation Expression pattern (CE) and other 932 genes with Top rank of Expression fold change (TE) were identified and the former were strongly associated with the water use efficiency (r ≥ 0.55, P ≤ 10−7). Functional enrichment of CE genes were related to three organelles, which well matched the annotation of twelve motifs identified from their conserved noncoding sequence; while TE genes were mostly related to biotic and/or abiotic stress. The expression robustness of CE genes with high genetic diversity kept relatively stable between environments while the harsh environment reduced the expression robustness of TE genes with low genetic diversity. The expression plasticity of CE genes was increased less than that of TE genes. These results suggested that expression variation of CE genes and TE genes could account for the robustness and plasticity of acclimation ability of Miscanthus, respectively. The patterns of expression variation revealed by transcriptomic network would shed new light on breeding and domestication of energy crops.