Long non-coding RNA (lncRNA) is a class of non-coding RNA with important regulatory roles in biological process of organisms. The systematic comparison of lncRNAs with protein coding mRNAs in population expression and their response to environmental change are still poorly understood. Here we identified 17,610 lncRNAs and calculated their expression levels based on RNA-seq of 80 individuals of Miscanthus lutarioriparius from two environments, the nearly native habitats and transplanted field, respectively. LncRNAs had significantly higher expression diversity and lower expression frequency in population than protein coding mRNAs in both environments, which suggested that lncRNAs may experience more relaxed selection or divergent evolution in population compared with protein coding RNAs. In addition, the increase of expression diversity for lncRNAs was always significantly higher and the magnitude of fold change of expression in new stress environment was significantly larger than protein-coding mRNAs. These results suggested that lncRNAs may be more sensitive to environmental change than protein-coding mRNAs. Analysis of environment-robust and environment-specific lncRNA-mRNA co-expression network between two environments revealed the characterization of lncRNAs in response to environmental change. Furthermore, candidate lncRNAs contributing to water use efficiency (WUE) identified based on the WUE-lncRNA-mRNA co-expression network suggested the roles of lncRNAs in response to environmental change. Our study provided a comprehensive understanding of expression characterization of lncRNAs in population for M. lutarioriparius under field condition, which would be useful to explore the roles of lncRNAs and could accelerate the process of adaptation in new environment for many plants.

Understanding the genetic basis of water use efficiency (WUE) and its roles in plant adaptation to a drought environment is essential for the production of second-generation energy crops in water-deficit marginal land. In this study, RNA-Seq and WUE measurements were performed for 78 individuals of Miscanthus lutarioriparius grown in two common gardens, one located in warm and wet Central China near the native habitats of the species and the other located in the semiarid Loess Plateau, the domestication site of the energy crop. The field measurements showed that WUE of M. lutarioriparius in the semiarid location was significantly higher than that in the wet location. A matrix correlation analysis was conducted between gene expression levels and WUE to identify candidate genes involved in the improvement of WUE from the native to the domestication site. A total of 48 candidate genes were identified and assigned to functional categories, including photosynthesis, stomatal regulation, protein metabolism, and abiotic stress responses. Of these genes, nearly 73% were up-regulated in the semiarid site. It was also found that the relatively high expression variation of the WUE-related genes was affected to a larger extent by environment than by genetic variation. The study demonstrates that transcriptome-wide correlation between physiological phenotypes and expression levels offers an effective means for identifying candidate genes involved in the adaptation to environmental changes.

As a key strategy for mitigating global climate change, bioenergy production by reducing CO2 emissions plays an important role in ensuring sustainable development. However, land-use change by converting natural ecosystems into energy crop field could create a carbon debt at the beginning. Thus, the potential carbon debt calculation is necessary for determining a promising bioenergy crop production, especially in the region rich of marginal land. Here, we used high-resolution historical land-use data to identify the marginal land available and to evaluate the carbon debt of planting Miscanthus in the Loess Plateau, China. We found that there were 27.6 Mha for energy production and 9.7 Mha for ecological restoration, with total annual production of 0.41 billion tons of biomass. We also found that soil carbon sequestration and total CO2 mitigation were 9.3 Mt C year−1 and 542 Mt year−1, respectively. More importantly, the result showed that planting Miscanthus on marginal land in the Loess Plateau only took 0.97 years on average to repay the carbon debt. Our study demonstrated that Miscanthus production in suitable marginal land in the Loess Plateau can offer considerable renewable energy and mitigate climate change with little carbon debt. These results suggested that bioenergy production in the similar arid and semiarid region worldwide would contribute to carbon sequestration in the context of rapid climate change.

As a potential energy crop with high biomass yield, Miscanthus lutarioriparius (M. lutarioriparius), endemic to the Long River Range in central China, needs to be investigated for its acclimation to stressful climatic and soil conditions often found on the marginal land. In this study, traits related to acclimation and yield, including survival rates, plant height (PH), stem diameter (SD), tiller number (TN), water use efficiency (WUE), and photosynthetic rates (A), were examined for 41 M. lutarioriparius populations that transplanted to the arid and cold Loess Plateau of China. The results showed that the average survival rate of M. lutarioriparius populations was only 4.16% over the first winter but the overwinter rate increased to 35.03% after the second winter, suggesting that plants having survived the first winter could have acclaimed to the low temperature. The strikingly high survival rates over the second winter were found to be 95.83% and 80.85%, respectively, for HG18 and HG39 populations. These populations might be especially valuable for the selection of energy crops for such an area. Those individuals surviving for the two consecutive winters showed significantly higher WUE than those measured after the first winter. The high WUE and low stomatal conductance (g_s) observed in survived individuals could have been responsible for their acclimation to this new and harsh environment. A total of 61 individuals with productive growth traits and strong resistance to cold and drought were identified for further energy crop development. This study showed that the variation of M. lutarioriparius held great potential for developing energy crops following continuous field selection.

As a promising candidate for the second-generation C4 energy crop, Miscanthus lutarioriparius has well acclimated to the water-limited and high-light Loess Plateau in China by improving photosynthesis rate and water use efficiency (WUE) compared to its native habitat along Yangtze River. Photosynthetic genes were demonstrated as one major category of the candidate genes underlying the physiological superiority. To further study how photosynthetic genes interact to improve the acclimation potential of M. lutarioriparius, population expression patterns within photosynthesis pathway were explored between one mild environment and one harsh environment. We found that 108 transcripts in assembled transcriptome of M. lutarioriparius were highly similar to genes in three Kyoto Encyclopedia of Genes and Genomes (KEGG) photosynthesis pathways of sorghum and maize. Phylogenetic analyses using sorghum, maize, rice and Arabidopsis genes of dark reaction identified 23 orthologs and 30 paralogs of M. lutarioriparius photosynthetic genes. These genes were also clustered into two kinds of expression pattern. 87% of transcripts in dark reaction were up-regulated and all 14 chloroplast-encoded transcripts in light reaction increased degradation in the harsh environment compared to the mild environment. Moreover, 80.8 % of photosynthetic transcripts were coordinated at transcription level under the two environments. Interestingly, LHCI and PSI were significantly correlated with F-ATPase and C4 cycle. Overall, this study indicates the coordinated expression between cyclic electron transport (consisting of LHCI, PSI and ATPase) and CO2-concentrating mechanism (C4 cycle) could account for photosynthesis plasticity on M. lutarioriparius acclimation potential.

Adaptation is a characteristic that enhances the survival or reproduction of organisms; selection is the critical process leading to adaptive evolution. Therefore, detecting selection is important in studying evolutionary biology. Changes in allele frequency are fundamental to adaptive evolution. The allele frequency of entire genes at the genomic scale is more intensive and precise for analyzing selection effects, compared with simple sequence repeat and single nucleotide polymorphism (SNP) alleles from nuclear gene fragments. Here, we analyzed 29,094 SNPs derived from 80 individuals of 14 Miscanthus lutarioriparius L. Liou ex S.L. Chen & Renvoize populations planted near their native habitat (Jiangxia, Hubei Province, JH) and a stressful environment (Qingyang, Gansu Province, QG) to detect selection during initial adaptation. The nucleotide diversity of over 60% of genes was decreased in QG compared with JH, suggesting that most genes were undergoing selection in the stressful environment. We explored a new approach based on haplotype data inferred from RNA‐seq data to analyze the change in frequency between two sites and to detect selection signals. In total, 402 and 51 genes were found to be targets of positive and negative selection, respectively. Among these candidate genes, the enrichment of abiotic stress‐response genes and photosynthesis‐related genes might have been responsible for establishment in the stressful environment. This is the first study assessing the change in allele frequency at the genomic level during adaptation. The method in which allele frequency detects selection during initial adaptation using population RNA‐seq data would be useful for developing evolutionary biology.