For most sequenced flowering plants, multiple whole-genome duplications (WGDs) are found. Duplicated genes following WGD often have different fates that can quickly disappear again, be retained for long(er) periods, or subsequently undergo small-scale duplications. However, how different expression, epigenetic regulation, and functional constraints are associated with these different gene fates following a WGD still requires further investigation due to successive WGDs in angiosperms complicating the gene trajectories. In this study, we investigate lotus (Nelumbo nucifera), an angiosperm with a single WGD during the K–pg boundary. Based on improved intraspecific-synteny identification by a chromosome-level assembly, transcriptome, and bisulfite sequencing, we explore not only the fundamental distinctions in genomic features, expression, and methylation patterns of genes with different fates after a WGD but also the factors that shape post-WGD expression divergence and expression bias between duplicates. We found that after a WGD genes that returned to single copies show the highest levels and breadth of expression, gene body methylation, and intron numbers, whereas the long-retained duplicates exhibit the highest degrees of protein–protein interactions and protein lengths and the lowest methylation in gene flanking regions. For those long-retained duplicate pairs, the degree of expression divergence correlates with their sequence divergence, degree in protein–protein interactions, and expression level, whereas their biases in expression level reflecting subgenome dominance are associated with the bias of subgenome fractionation. Overall, our study on the paleopolyploid nature of lotus highlights the impact of different functional constraints on gene fate and duplicate divergence following a single WGD in plant.

Background: Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNAmediated gene inhibition and alternative-splicing of pre-mRNAs. However, the coordination of interaction between miRNAs with different splicing isoforms, and the change of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus (Nelumbo nucifera), an economically important flowering plant. Results: Through RNA-seq analyses on miRNAs and their target genes (isoforms) among six lotus tissues, expression of most miRNAs seem to be negatively correlated with their targets and tend to be tissue-specific. Further, our results showed that preferential interactions between miRNAs and hub gene isoforms in one coexpression module which is highly correlated with leaf. Intriguingly, for many genes, their corresponding isoforms were assigned to different coexpressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. Further detailed functional enrichment analysis of miRNA targets revealed that miRNAs are involved in the regulation of lotus growth and development by regulating plant hormone-related pathway genes. Conclusions: Taken together, our comprehensive analyses of miRNA and mRNA transcriptome elucidate the coordination of interaction between miRNAs and different splicing isoforms, and highlight the functional divergence of many transcript isoforms from the same locus in lotus.

No abstract

Background: Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNA-mediated gene inhibition and alternative-splicing of pre-mRNAs. However, the coordination of interaction between miRNAs with different splicing isoforms, and the change of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus (Nelumbo nucifera), an economically important flowering plant.Results: Through RNA-seq analyses on miRNAs and their target genes (isoforms) among six lotus tissues, expression of most miRNAs seem to be negatively correlated with their targets and tend to be tissue-specific. Further, our results showed that preferential interactions between miRNAs and hub gene isoforms in one coexpression module which is highly correlated with leaf. Intriguingly, for many genes, their corresponding isoforms were assigned to different co-expressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. Further detailed functional enrichment analysis of miRNA targets revealed that miRNAs are involved in the regulation of lotus growth and development by regulating plant hormone-related pathway genes.Conclusions: Taken together, our comprehensive analyses of miRNA and mRNA transcriptome elucidate the coordination of interaction between miRNAs and different splicing isoforms, and highlight the functional divergence of many transcript isoforms from the same locus in lotus.

Background Brassinosteroid (BR) signaling regulates plant growth and development in concert with other signaling pathways. Although many genes have been identified that play a role in BR signaling, the biological and functional consequences of disrupting those key BR genes still require detailed investigation. Results Here we performed phenotypic and transcriptomic comparisons of A. thaliana lines carrying a loss-of-function mutation in BRI1 gene, bri1–5, that exhibits a dwarf phenotype and its three activation-tag suppressor lines that were able to partially revert the bri1–5 mutant phenotype to a WS2 phenotype, namely bri1–5/bri1–1D, bri1–5/brs1–1D, and bri1–5/bak1–1D. From the three investigated bri1–5 suppressors, bri1–5/bak1–1D was the most effective suppressor at the transcriptional level. All three bri1–5 suppressors showed altered expression of the genes in the abscisic acid (ABA signaling) pathway, indicating that ABA likely contributes to the partial recovery of the wild-type phenotype in these bri1–5 suppressors. Network analysis revealed crosstalk between BR and other phytohormone signaling pathways, suggesting that interference with one hormone signaling pathway affects other hormone signaling pathways. In addition, differential expression analysis suggested the existence of a strong negative feedback from BR signaling on BR biosynthesis and also predicted that BRS1, rather than being directly involved in signaling, might be responsible for providing an optimal environment for the interaction between BRI1 and its ligand. Conclusions Our study provides insights into the molecular mechanisms and functions of key brassinosteroid (BR) signaling genes, especially BRS1.

BackgroundBrassinosteroid (BR) signaling regulates plant growth and development in concert with other signaling pathways. Although many genes have been identified that play a role in Brassinosteroid (BR) signaling, the biological and functional consequences of disrupting those key BR genes still requires detailed investigation.ResultsHere we performed phenotypic and transcriptomic comparisons of A. thaliana lines carrying a loss-of-function mutation in BRI1 gene, bri1-5, that exhibits a dwarf phenotype along with its three activation-tag suppressor lines that were able to partially revert the bri1-5 mutant phenotype to a WT phenotype, namely bri1-5/bri1-1D, bri1-5/brs1-1D, bri1-5/bak1-1D. From the three investigated bri1-5 suppressors, bri1-5/bak1-1D was the most effective suppressor at the transcriptional level. All three bri1-5 suppressors showed altered expression of the genes in the abscisic acid (ABA signaling) pathway, indicating that ABA likely contributes to the partial recovery of the wild type phenotype in these bri1-5 suppressors. Network analysis revealed crosstalk between BR and other phytohormone signaling pathways, suggesting that interference with one hormone signaling pathway affects other hormone signaling pathways. In addition, differential expression analysis suggested the existence of a strong negative feedback from BR signaling on BR biosynthesis and also predicted that BRS1, rather than being directly involved in signaling, is likely responsible for providing an optimal environment for the interaction between BRI1 and its ligand. ConclusionsOur study provides insights into the molecular mechanisms and functions of key brassinosteroid (BR) signaling genes, especially BRS1.

With the decreasing cost of sequencing and availability of larger numbers of sequenced genomes, comparative genomics is becoming increasingly attractive to complement experimental techniques for the task of transcription factor binding site identification. In this study, we redesigned BLSSpeller, a motif discovery algorithm, to cope with larger sequence datasets. BLSSpeller was used to identify novel motifs in Zea mays in a comparative genomics setting with 16 monocot lineages. We discovered 61 motifs of which 20 matched previously described motif models in Arabidopsis. In addition, novel, yet uncharacterized motifs were detected, several of which are supported by available sequence-based and/or functional data. Instances of the predicted motifs were enriched around transcription start sites and contained signatures of selection. Moreover, the enrichment of the predicted motif instances in open chromatin and transcription factor binding sites indicates their functionality, supported by the fact that genes carrying instances of these motifs were often found to be coexpressed and/or enriched in similar GO functions. Overall, our study unveiled several novel candidate motifs that might help our understanding of the genotype to phenotype association in crops.

BackgroundA better understanding of the molecular effects of salinity stress is key to improving salt tolerance in Zea mays. In this study, we combined phenotyping with transcript profiling and network analysis to study genotype-specific differences in salt tolerance in Zea mays. ResultAn extensive phenotypic screening identified two genotypes with an extreme phenotypic difference in tolerance towards salt stress. RNA-seq analysis of the selected salt-tolerant (R9) and salt-sensitive (S46) genotype was performed to unveil the molecular mechanism underlying the difference in salt tolerance. GO enrichment and network analysis on the results of the expression analysis identified phosphorylation-dependent signaling processes, ion transportation, oxidation-reduction, glutathione and tryptophan metabolism as the main processes different between the selected tolerant and sensitive genotypes. Genes belonging to the subnetwork enriched for phosphorylation and kinase activity shared a common regulatory element in their promoter region, which matched the binding site of an Arabidopsis TF with known role in salt-stress response.ConclusionNetwork-based transcriptome analysis of two maize genotypes identified pathways associated with differences in genotype-specific salt tolerance and identified a link between transcriptional and posttranslational regulation of salt tolerance.