Genome-Wide Analysis Indicates Lineage-Specific Gene Loss during Papilionoideae Evolution

Gene loss is the driving force for changes in genome and morphology; however, this particular evolutionary event has been poorly investigated in leguminous plants. Legumes (Fabaceae) have some lineage-specific and diagnostic characteristics that are distinct from other angiosperms. To understand the potential role of gene loss in the evolution of legumes, we compared six genome-sequenced legume species of Papilionoideae, the largest representative clade of Fabaceae, such as Glycine max, with 34 nonlegume plant species, such as Arabidopsis thaliana. The results showed that the putative orthologs of the 34 Arabidopsis genes belonging to 29 gene families were absent in these legume species but these were conserved in the sequenced nonlegume angiosperm lineages. Further evolutionary analyses indicated that the orthologs of these genes were almost completely lost in the Papillionoideae ancestors, thus designated as the legume lost genes (LLGs), and these underwent purifying selection in nonlegume plants. Most LLGs were functionally unknown. In Arabidopsis, two LLGs were well-known genes that played a role in plant immunity such as HARMLESS TO OZONE LAYER 1 and HOPZ-ACTIVATED RESISTANCE 1, and 16 additional LLGs were predicted to participate in plant-pathogen interactions in in silico expression and protein-protein interaction network analyses. Most of these LLGs' orthologs in various plants were also found to be associated with biotic stress response, indicating the conserved role of these genes in plant defense. The evolutionary implication of LLGs during the development of the ability of symbiotic nitrogen fixation involving plant and bacterial interactions, which is a well-known characteristic of most legumes, is also discussed. Our work sheds light on the evolutionary implication of gene loss events in Papilionoideae evolution, as well as provides new insights into crop design to improve nitrogen fixation capacity

Genome-Wide Analysis Indicates Lineage-Specific Gene Loss during Papilionoideae Evolution

Gene loss is the driving force for changes in genome and morphology; however, this particular evolutionary event has been poorly investigated in leguminous plants. Legumes (Fabaceae) have some lineage-specific and diagnostic characteristics that are distinct from other angiosperms. To understand the potential role of gene loss in the evolution of legumes, we compared six genome-sequenced legume species of Papilionoideae, the largest representative clade of Fabaceae, such as Glycine max, with 34 nonlegume plant species, such as Arabidopsis thaliana. The results showed that the putative orthologs of the 34 Arabidopsis genes belonging to 29 gene families were absent in these legume species but these were conserved in the sequenced nonlegume angiosperm lineages. Further evolutionary analyses indicated that the orthologs of these genes were almost completely lost in the Papillionoideae ancestors, thus designated as the legume lost genes (LLGs), and these underwent purifying selection in nonlegume plants. Most LLGs were functionally unknown. In Arabidopsis, two LLGs were well-known genes that played a role in plant immunity such as HARMLESS TO OZONE LAYER 1 and HOPZ-ACTIVATED RESISTANCE 1, and 16 additional LLGs were predicted to participate in plant–pathogen interactions in in silico expression and protein–protein interaction network analyses. Most of these LLGs’ orthologs in various plants were also found to be associated with biotic stress response, indicating the conserved role of these genes in plant defense. The evolutionary implication of LLGs during the development of the ability of symbiotic nitrogen fixation involving plant and bacterial interactions, which is a well-known characteristic of most legumes, is also discussed. Our work sheds light on the evolutionary implication of gene loss events in Papilionoideae evolution, as well as provides new insights into crop design to improve nitrogen fixation capacity

Population transcriptomics uncovers the regulation of gene expression variation in adaptation to changing environment

Expression variation plays an important role in plant adaptation, but little is known about the factors impacting the expression variation when population adapts to changing environment. We used RNA-seq data from 80 individuals in 14 Miscanthus lutarioriparius populations, which were transplanted into a harsh environment from native habitat, to investigate the expression level, expression diversity and genetic diversity for genes expressed in both environments. The expression level of genes with lower expression level or without SNP tended to be more changeable in new environment, which suggested highly expressed genes experienced stronger purifying selection than those at lower level. Low proportion of genes with population effect confirmed the weak population structure and frequent gene flow in these populations. Meanwhile, the number of genes with environment effect was the most frequent compared with that with population effect. Our results showed that environment and genetic diversity were the main factors determining gene expression variation in population. This study could facilitate understanding the mechanisms of global gene expression variation when plant population adapts to changing environment

Population transcriptomics uncovers the regulation of gene expression variation in adaptation to changing environment

Expression variation plays an important role in plant adaptation, but little is known about the factors impacting the expression variation when population adapts to changing environment. We used RNA-seq data from 80 individuals in 14 Miscanthus lutarioriparius populations, which were transplanted into a harsh environment from native habitat, to investigate the expression level, expression diversity and genetic diversity for genes expressed in both environments. The expression level of genes with lower expression level or without SNP tended to be more changeable in new environment, which suggested highly expressed genes experienced stronger purifying selection than those at lower level. Low proportion of genes with population effect confirmed the weak population structure and frequent gene flow in these populations. Meanwhile, the number of genes with environment effect was the most frequent compared with that with population effect. Our results showed that environment and genetic diversity were the main factors determining gene expression variation in population. This study could facilitate understanding the mechanisms of global gene expression variation when plant population adapts to changing environment

Coexpression network revealing the plasticity and robustness of population transcriptome during the initial stage of domesticating energy crop Miscanthus lutarioriparius

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 ae 0.55, P ae 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

Expression of miR-652-3p and Effect on Apoptosis and Drug Sensitivity in Pediatric Acute Lymphoblastic Leukemia

MicroRNAs (miRNAs) expression profiles were screened in plasma samples from pediatric patients with acute lymphoblastic leukemia (ALL) and healthy controls, using qRT-PCR-based TaqMan low-density miRNA arrays. MiR-652-3p (a circulating miRNA) was downregulated in new diagnosis (ND) patients compared with healthy controls. The levels of miR652-3p were restored in complete remission (CR) but were downregulated again in disease relapse (RE). The expression pattern of miR-652-3p was validated in bone marrow (BM) samples from other pediatric ALL patients. MiR-652-3p was significantly upregulated in BM when the patients (n=86) achieved CR, as compared with the matched ND samples (p<0.001). Moreover, the miR-652-3p levels in BM decreased again in two patients at RE. In addition, the lymphoblastic leukemia cell lines Reh and RS4:11 were found to have lower levels of miR-625-3p than the normal B-cell line. Overexpression of miR-652-3p using agomir increased the sensitivity to vincristine and cytarabine (all p<0.05) and promoted apoptosis (both p<0.05) in Reh and RS4:11 cells. In conclusion, the results suggested that a low level of miR-652-3p might be involved in the pathogenesis of pediatric ALL. Overexpression of miR-652-3p might suppress lymphoblastic leukemia cells, promoting apoptosis and increasing sensitivity to chemotherapeutic drugs

Population transcriptomic characterization of the genetic and expression variation of a candidate progenitor of Miscanthus energy crops

The use of transcriptome data in the study of the population genetics of a species can capture faint signals of both genetic variation and expression variation and can provide a broad picture of a species' genomic response to environmental conditions. In this study, we characterized the genetic and expression diversity of Miscanthus lutarioriparius by comparing more than 16,225 transcripts obtained from 78 individuals, belonging to 10 populations distributed across the species' entire geographic range. We only observed a low level of nucleotide diversity (pi = 0.000434) among the transcriptome data of these populations, which is consistent with highly conserved sequences of functional elements and protein-coding genes captured with this method. Tests of population divergence using the transcriptome data were consistent with previous microsatellite data but proved to be more sensitive, particularly if gene expression variation was considered as well. For example, the analysis of expression data showed that genes involved in photosynthetic processes and responses to temperature or reactive oxygen species stimuli were significantly enriched in certain populations. This differential gene expression was primarily observed among populations and not within populations. Interestingly, nucleotide diversity was significantly negatively correlated with expression diversity within populations, while this correlation was positive among populations. This suggests that genetic and expression variation play separate roles in adaptation and population persistence. Combining analyses of genetic and gene expression variation represents a promising approach for studying the population genetics of wild species and may uncover both adaptive and nonadaptive processes

The coordination of gene expression within photosynthesis pathway for acclimation of C4 energy crop Miscanthus lutarioriparius

As a promising candidate for the second-generation C(4) 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 C(4) cycle. Overall, this study indicates the coordinated expression between cyclic electron transport (consisting of LHCI, PSI, and ATPase) and CO(2)-concentrating mechanism (C(4) cycle) could account for photosynthesis plasticity on M. lutarioriparius acclimation potential

Population transcriptomic characterization of the genetic and expression variation of a candidate progenitor of Miscanthus energy crops

The use of transcriptome data in the study of the population genetics of a species can capture faint signals of both genetic variation and expression variation, and can provide a broad picture of a species’ genomic response to environmental conditions. In this study, we characterized the genetic and expression diversity of Miscanthus lutarioriparius by comparing more than 16,225 transcripts obtained from 78 individuals, belonging to 10 populations distributed across the species’ entire geographic range. We only observed a low level of nucleotide diversity (π = 0.000434) among the transcriptome data of these populations, which is consistent with highly conserved sequences of functional elements and protein-coding genes captured with this method. Tests of population divergence using the transcriptome data were consistent with previous microsatellite data but proved to be more sensitive, particularly if gene expression variation was considered as well. For example, the analysis of expression data showed that genes involved in photosynthetic processes and responses to temperature or reactive oxygen species stimuli were significantly enriched in certain populations. This differential gene expression was primarily observed among populations and not within populations. Interestingly, nucleotide diversity was significantly negatively correlated with expression diversity within populations, while this correlation was positive among populations. This suggests that genetic and expression variation play separate roles in adaptation and population persistence. Combining analyses of genetic and gene expression variation represents a promising approach for studying the population genetics of wild species, and may uncover both adaptive and non-adaptive processes