euAP2a, a key gene that regulates flowering time in peach (Prunus persica) by modulating thermo-responsive transcription programming.

Frequent spring frost damage threatens temperate fruit production, and breeding of late-flowering cultivars is an effective strategy for preventing such damage. However, this effort is often hampered by the lack of specific genes and markers and a lack of understanding of the mechanisms. We examined a Late-Flowering Peach (LFP) germplasm and found that its floral buds require a longer chilling period to release from their dormancy and a longer warming period to bloom than the control cultivar, two key characteristics associated with flowering time. We discovered that a 983-bp deletion in euAP2a, an APETALA2 (AP2)-related gene with known roles in regulating floral organ identity and flowering time, was primarily responsible for late flowering in LFP. This deletion disrupts an miR172 binding site, resulting in a gain-of-function mutation in euAP2a. Transcriptomic analyses revealed that at different stages of floral development, two chilling-responsive modules and four warm-responsive modules, comprising approximately 600 genes, were sequentially activated, forming a unique transcription programming. Furthermore, we found that euAP2a was transiently downregulated during the activation of these thermal-responsive modules at various stages. However, the loss of such transient, stage-specific downregulation of euAP2a caused by the deletion of miR172 binding sites resulted in the deactivation or delay of these modules in the LFP flower buds, suggesting that euAP2a acts as a transcription repressor to control floral developmental pace in peaches by modulating the thermo-responsive transcription programming. The findings shed light on the mechanisms behind late flowering in deciduous fruit trees, which is instrumental for breeding frost-tolerant cultivars.

Regulatory frameworks involved in the floral induction, formation and developmental programming of woody horticultural plants: a case study on blueberries.

Flowering represents a crucial stage in the life cycles of plants. Ensuring strong and consistent flowering is vital for maintaining crop production amidst the challenges presented by climate change. In this review, we summarized key recent efforts aimed at unraveling the complexities of plant flowering through genetic, genomic, physiological, and biochemical studies in woody species, with a special focus on the genetic control of floral initiation and activation in woody horticultural species. Key topics covered in the review include major flowering pathway genes in deciduous woody plants, regulation of the phase transition from juvenile to adult stage, the roles of CONSTANS (CO) and CO-like gene and FLOWERING LOCUS T genes in flower induction, the floral regulatory role of GA-DELLA pathway, and the multifunctional roles of MADS-box genes in flowering and dormancy release triggered by chilling. Based on our own research work in blueberries, we highlighted the central roles played by two key flowering pathway genes, FLOWERING LOCUS T and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1, which regulate floral initiation and activation (dormancy release), respectively. Collectively, our survey shows both the conserved and diverse aspects of the flowering pathway in annual and woody plants, providing insights into the potential molecular mechanisms governing woody plants. This paves the way for enhancing the resilience and productivity of fruit-bearing crops in the face of changing climatic conditions, all through the perspective of genetic interventions.

CW198 acts as a genetic insulator to block enhancer-promoter interaction in plants.

Insulators in vertebrates play a role in genome architecture and orchestrate temporo-spatial enhancer-promoter interactions. In plants, insulators and their associated binding factors have not been documented as of yet, largely as a result of a lack of characterized insulators. In this study, we took a comprehensive strategy to identify and validate the enhancer-blocking insulator CW198. We show that a 1.08-kb CW198 fragment from Arabidopsis can, when interposed between an enhancer and a promoter, efficiently abrogate the activation function of both constitutive and floral organ-specific enhancers in transgenic Arabidopsis and tobacco plants. In plants, both transcriptional crosstalk and spreading of histone modifications were rarely detectable across CW198, which resembles the insulation property observed across the CTCF insulator in the mammalian genome. Taken together, our findings support that CW198 acts as an enhancer-blocking insulator in both Arabidopsis and tobacco. The significance of the present findings and their relevance to the mitigation of mutual interference between enhancers and promoters, as well as multiple promoters in transgenes, is discussed.

An Improved Procedure for Agrobacterium-Mediated Transformation of 'Carrizo' Citrange.

Although several protocols for genetic transformation of citrus have been published, it is highly desirable to further improve its efficiency. Here we report treatments of Agrobacterium cells and citrus explants prior to and during co-cultivation process to enhance transformation efficiency using a commercially used rootstock 'Carrizo' citrange [Citrus sinensis (L.) Osb. × Poncirius trifoliata (L.) Raf.] as a model plant. We found explants from light-grown seedlings exhibited higher transformation efficiency than those from etiolated seedlings. We pre-cultured Agrobacterium cells in a 1/10 MS, 0.5 g/L 2-(N-morpholino) ethanesulfonic acid (MES) and 100 µM acetosyringone liquid medium for 6 h at 25 °C before used to infect citrus explants. We incubated epicotyl segments in an MS liquid medium containing 13.2 µM 6-BA, 4.5 µM 2,4-D, 0.5 µM NAA for 3 h at 25 °C prior to Agrobacterium infection. In the co-cultivation medium, we added 30 µM paclobutrazol and 10 µM lipoic acid. Each of these treatments significantly increased the efficiencies of transformation up to 30.4% (treating Agrobacterium with acetosyringone), 31.8% (treating explants with cytokinin and auxin), 34.9% (paclobutrazol) and 38.6% (lipoic acid), respectively. When the three treatments were combined, we observed that the transformation efficiency was enhanced from 11.5% to 52.3%. The improvement of genetic transformation efficiency mediated by these three simple treatments may facilitate more efficient applications of transgenic and gene editing technologies for functional characterization of citrus genes and for genetic improvement of citrus cultivars.

Rice RS2-9, which is bound by transcription factor OSH1, blocks enhancer-promoter interactions in plants.

Insulators characterized in Drosophila and mammals have been shown to play a key role in the restriction of promiscuous enhancer-promoter interactions, as well as reshaping the topological landscape of chromosomes. Yet the role of insulators in plants remains poorly understood, in large part because of a lack of well-characterized insulators and binding factor(s). In this study, we isolated a 1.2-kb RS2-9 insulator from the Oryza sativa (rice) genome that can, when interposed between an enhancer and promoter, efficiently block the activation function of both constitutive and floral organ-specific enhancers in transgenic Arabidopsis and Nicotiana tabacum (tobacco). In the rice genome, the genes flanking RS2-9 exhibit an absence of mutual transcriptional interactions, as well as a lack of histone modification spread. We further determined that O. sativa Homeobox 1 (OSH1) bound two regions of RS2-9, as well as over 50 000 additional sites in the rice genome, the majority of which resided in intergenic regions. Mutation of one of the two OSH1-binding sites in RS2-9 impaired insulation activity by up to 60%, whereas the mutation of both binding sites virtually abolished insulator function. We also demonstrated that OSH1 binding sites were associated with 72% of the boundaries of topologically associated domains (TADs) identified in the rice genome, which is comparable to the 77% of TAD boundaries bound by the insulator CCCTC-binding factor (CTCF) in mammals. Taken together, our findings indicate that OSH1-RS2-9 acts as a true insulator in plants, and highlight a potential role for OSH1 in gene insulation and topological organization in plant genomes.

Genome-Wide Changes of Regulatory Non-Coding RNAs Reveal Pollen Development Initiated at Ecodormancy in Peach.

Bud dormancy is under the regulation of complex mechanisms including genetic and epigenetic factors. To study the function of regulatory non-coding RNAs in winter dormancy release, we analyzed the small RNA and long non-coding RNA (lncRNA) expression from peach (Prunus persica) floral buds in endodormancy, ecodormancy and bud break stages. Small RNAs underwent a major shift in expression primarily between dormancy and flowering with specific pairs of microRNAs and their mRNA target genes undergoing coordinated differential expression. From endodormancy to ecodormancy, ppe-miR6285 was significantly upregulated while its target gene, an ASPARAGINE-RICH PROTEIN involved in the regulation of abscisic acid signaling, was downregulated. At ecodormancy, ppe-miR2275, a homolog of meiosis-specific miR2275 across angiosperms, was significantly upregulated, supporting microsporogenesis in anthers at a late stage of dormancy. The expression of 785 lncRNAs, unlike the overall expression pattern in the small RNAs, demonstrated distinctive expression signatures across all dormancy and flowering stages. We predicted that a subset of lncRNAs were targets of microRNAs and found 18 lncRNA/microRNA target pairs with both differentially expressed across time points. The genome-wide differential expression and network analysis of non-coding RNAs and mRNAs from the same tissues provide new candidate loci for dormancy regulation and suggest complex noncoding RNA interactions control transcriptional regulation across these key developmental time points.

Thermal-responsive genetic and epigenetic regulation of DAM cluster controlling dormancy and chilling requirement in peach floral buds.

The Dormancy-associated MADS-box (DAM) gene cluster in peach serves as a key regulatory hub on which the seasonal temperatures act and orchestrate dormancy onset and exit, chilling response and floral bud developmental pace. Yet, how different temperature regimes interact with and regulate the six linked DAM genes remains unclear. Here, we demonstrate that chilling downregulates DAM1 and DAM3-6 in dormant floral buds with distinct patterns and identify DAM4 as the most abundantly expressed one. We reveal multiple epigenetic events, with tri-methyl histone H3 lysine 27 (H3K27me3) induced by chilling specifically in DAM1 and DAM5, a 21-nt sRNA in DAM3 and a ncRNA induced in DAM4. Such induction is inversely correlated with downregulation of their cognate DAMs. We also show that the six DAMs were hypermethylated, associating with the production of 24-nt sRNAs. Hence, the chilling-responsive dynamic of the different epigenetic elements and their interactions likely define distinct expression abundance and downregulation pattern of each DAM. We further show that the expression of the five DAMs remains steadily unchanged or continuously downregulated at the ensuing warm temperature after chilling, and this state of regulation correlates with robust increase of sRNA expression, H3K27me3 and CHH methylation, which is particularly pronounced in DAM4. Such robust increase of repressive epigenetic marks may irreversibly reinforce the chilling-imposed repression of DAMs to ensure flower-developmental programming free from any residual DAM inhibition. Taken together, we reveal novel information about genetic and epigenetic regulation of the DAM cluster in peach, which will be of fundamental significance in understanding of the regulatory mechanisms underlying chilling requirement and dormancy release, and of practical application for improvement of plasticity of flower time and bud break in fruit trees to adapt changing climates.

Distinctive Gene Expression Patterns Define Endodormancy to Ecodormancy Transition in Apricot and Peach.

Dormancy is a physiological state that plants enter for winter hardiness. Environmental-induced dormancy onset and release in temperate perennials coordinate growth cessation and resumption, but how the entire process, especially chilling-dependent dormancy release and flowering, is regulated remains largely unclear. We utilized the transcriptome profiles of floral buds from fall to spring in apricot (Prunus armeniaca) genotypes with contrasting bloom dates and peach (Prunus persica) genotypes with contrasting chilling requirements (CR) to explore the genetic regulation of bud dormancy. We identified distinct gene expression programming patterns in endodormancy and ecodormancy that reproducibly occur between different genotypes and species. During the transition from endo- to eco-dormancy, 1,367 and 2,102 genes changed in expression in apricot and peach, respectively. Over 600 differentially expressed genes were shared in peach and apricot, including three DORMANCY ASSOCIATED MADS-box (DAM) genes (DAM4, DAM5, and DAM6). Of the shared genes, 99 are located within peach CR quantitative trait loci, suggesting these genes as candidates for dormancy regulation. Co-expression and functional analyses revealed that distinctive metabolic processes distinguish dormancy stages, with genes expressed during endodormancy involved in chromatin remodeling and reproduction, while the genes induced at ecodormancy were mainly related to pollen development and cell wall biosynthesis. Gene expression analyses between two Prunus species highlighted the conserved transcriptional control of physiological activities in endodormancy and ecodormancy and revealed genes that may be involved in the transition between the two stages.

Association of the phenylpropanoid pathway with dormancy and adaptive trait variation in apricot (Prunus armeniaca).

Trees use many mechanisms to adapt and respond to stressful conditions. The phenylpropanoid pathway in particular is known to be associated with a diverse suite of plant stress responses. In this study, we explored the relationship between the phenylpropanoid pathway metabolite production, gene expression and adaptive trait variation associated with floral bud reactivation during and following dormancy in Prunus armeniaca L. (apricot). Concentrations of eight phenylpropanoid metabolites were measured during chill accumulation and at developmental stages corresponding to the emergence of sepals and petals in floral buds of varieties that differ phenotypically in bloom date (BD). A significant interaction effect of chill hours and BD phenotype on the concentration of each of the compounds was observed (mixed analysis of variance, P < 0.05), with the concentration of most phenylpropanoid metabolites dropping precipitously when sepals and petals emerged. While phenylpropanoid biosynthetic gene expression patterns were more variable in general, expression changed over time and was impacted, although to a lesser degree, by BD phenotype. Furthermore, separation of BD phenotypic groups was most pronounced when early and late BD varieties were at different developmental stages, i.e., 800 chill hours. Taken together, these results suggest that the phenylpropanoid pathway is associated with floral bud reactivation in apricot. Furthermore, we show that the phenylpropanoid pathway is also impacted by phenological trait variation associated with dormancy. A better understanding of how apricot and other perennial tree species respond and adapt to environmental perturbations will be critical for improvement programs aimed at identifying and breeding trees more suitable for rapidly changing environments.

Small RNAs, emerging regulators critical for the development of horticultural traits.

Small RNAs (sRNAs) have been recently recognized as key genetic and epigenetic regulators in various organisms, ranging from the modification of DNA and histone methylations to the modulation of the abundance of coding or non-coding RNAs. In plants, major regulatory sRNAs are classified as respective microRNA (miRNA) and small interfering RNA (siRNA) species, with the former primarily engaging in posttranscriptional regulation while the latter in transcriptional one. Many of these characterized sRNAs are involved in regulation of diverse biological programs, processes, and pathways in response to developmental cues, environmental signals/stresses, pathogen infection, and pest attacks. Recently, sRNAs-mediated regulations have also been extensively investigated in horticultural plants, with many novel mechanisms unveiled, which display far more mechanistic complexity and unique regulatory features compared to those studied in model species. Here, we review the recent progress of sRNA research in horticultural plants, with emphasis on mechanistic aspects as well as their relevance to trait regulation. Given that major and pioneered sRNA research has been carried out in the model and other plants, we also discuss ongoing sRNA research on these plants. Because miRNAs and phased siRNAs (phasiRNAs) are the most studied sRNA regulators, this review focuses on their biogenesis, conservation, function, and targeted genes and traits as well as the mechanistic relation between them, aiming at providing readers comprehensive information instrumental for future sRNA research in horticulture crops.

Author Correction: Dynamic Changes of Genome-Wide DNA Methylation during Soybean Seed Development.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Characterization of a new apple luteovirus identified by high-throughput sequencing.

BACKGROUND: 'Rapid Apple Decline' (RAD) is a newly emerging problem of young, dwarf apple trees in the Northeastern USA. The affected trees show trunk necrosis, cracking and canker before collapse in summer. In this study, we discovered and characterized a new luteovirus from apple trees in RAD-affected orchards using high-throughput sequencing (HTS) technology and subsequent Sanger sequencing. METHODS: Illumina NextSeq sequencing was applied to total RNAs prepared from three diseased apple trees. Sequence reads were de novo assembled, and contigs were annotated by BLASTx. RT-PCR and 5'/3' RACE sequencing were used to obtain the complete genome of a new virus. RT-PCR was used to detect the virus. RESULTS: Three common apple viruses and a new luteovirus were identified from the diseased trees by HTS and RT-PCR. Sequence analyses of the complete genome of the new virus show that it is a new species of the genus Luteovirus in the family Luteoviridae. The virus is graft transmissible and detected by RT-PCR in apple trees in a couple of orchards. CONCLUSIONS: A new luteovirus and/or three known viruses were found to be associated with RAD. Molecular characterization of the new luteovirus provides important information for further investigation of its distribution and etiological role.

Exposure in vitro to an Environmentally Isolated Strain TC09 of Cladosporium sphaerospermum Triggers Plant Growth Promotion, Early Flowering, and Fruit Yield Increase.

A growing number of bacteria and fungi have been found to promote plant growth through mutualistic interactions involving elements such as volatile organic compounds (VOCs). Here, we report the identification of an environmentally isolated strain of Cladosporium sphaerospermum (herein named TC09), that substantially enhances plant growth after exposure in vitro beyond what has previously been reported. When cultured on Murashige and Skoog (MS) medium under in vitro conditions, tobacco seedlings (Nicotiana tabacum) exposed to TC09 cultures for 20 days increased stem height and whole plant biomass up to 25- and 15-fold, respectively, over controls without exposure. TC09-mediated growth promotion required >5 g/L sucrose in the plant culture medium and was influenced by the duration of exposure ranging from one to 10 days, beyond which no differences were detected. When transplanted to soil under greenhouse conditions, TC09-exposed tobacco plants retained higher rates of growth. Comparative transcriptome analyses using tobacco seedlings exposed to TC09 for 10 days uncovered differentially expressed genes (DEGs) associated with diverse biological processes including cell expansion and cell cycle, photosynthesis, phytohormone homeostasis and defense responses. To test the potential efficacy of TC09-mediated growth promotion on agricultural productivity, pepper plants (Capsicum annuum L.) of two different varieties, Cayenne and Minisweet, were pre-exposed to TC09 and planted in the greenhouse to monitor growth, flowering, and fruit production. Results showed that treated pepper plants flowered 20 days earlier and yielded up to 213% more fruit than untreated controls. Altogether the data suggest that exposure of young plants to C. sphaerospermum produced VOCs may provide a useful tool to improve crop productivity.

Dynamic Changes of Genome-Wide DNA Methylation during Soybean Seed Development.

Seed development is programmed by expression of many genes in plants. Seed maturation is an important developmental process to soybean seed quality and yield. DNA methylation is a major epigenetic modification regulating gene expression. However, little is known about the dynamic nature of DNA methylation and its effects on gene expression during plant development. Through whole-genome bisulfite sequencing, we showed that DNA methylation went through dynamic changes during seed maturation. An average of 66% CG, 45% CHG and 9% CHH contexts was methylated in cotyledons. CHH methylation levels in cotyledons changed greatly from 6% at the early stage to 11% at the late stage. Transcribed genes were approximately two-fold more likely to be differentially methylated than non-transcribed genes. We identified 40, 66 and 2136 genes containing differentially methylated regions (DMRs) with negative correlation between their expression and methylation in the CG, CHG and CHH contexts, respectively. The majority of the DMR genes in the CHH context were transcriptionally down-regulated as seeds mature: 99% of them during early maturation were down-regulated, and preferentially associated with DNA replication and cell division. The results provide novel insights into the dynamic nature of DNA methylation and its relationship with gene regulation in seed development.

SacB-SacR Gene Cassette As the Negative Selection Marker to Suppress Agrobacterium Overgrowth in Agrobacterium-Mediated Plant Transformation.

Agrobacterium overgrowth is a common problem in Agrobacterium-mediated plant transformation. To suppress the Agrobacterium overgrowth, various antibiotics have been used during plant tissue culture steps. The antibiotics are expensive and may adversely affect plant cell differentiation and reduce plant transformation efficiency. The SacB-SacR proteins are toxic to most Agrobacterium tumefaciens strains when they are grown on culture medium supplemented with sucrose. Therefore, SacB-SacR genes can be used as negative selection markers to suppress the overgrowth of A. tumefaciens in the plant tissue culture process. We generated a mutant A. tumefaciens strain GV2260 (recA-SacB/R) that has the SacB-SacR cassette inserted into the bacterial genome at the recA gene locus. The mutant Agrobacterium strain is sensitive to sucrose but maintains its ability to transform plant cells in both transient and stable transformation assays. We demonstrated that the mutant strain GV2260 (recA-SacB/R) can be inhibited by sucrose that reduces the overgrowth of Agrobacterium and therefore improves the plant transformation efficiency. We employed GV2260 (recA-SacB/R) to generate stable transgenic N. benthamiana plants expressing a CRISPR-Cas9 for knocking out a WRKY transcription factor.

A Novel Two-Step Method for Screening Shade Tolerant Mutant Plants via Dwarfism.

When subjected to shade, plants undergo rapid shoot elongation, which often makes them more prone to disease and mechanical damage. Shade-tolerant plants can be difficult to breed; however, they offer a substantial benefit over other varieties in low-light areas. Although perennial ryegrass (Lolium perenne L.) is a popular species of turf grasses because of their good appearance and fast establishment, the plant normally does not perform well under shade conditions. It has been reported that, in turfgrass, induced dwarfism can enhance shade tolerance. Here we describe a two-step procedure for isolating shade tolerant mutants of perennial ryegrass by first screening for dominant dwarf mutants, and then screening dwarf plants for shade tolerance. The two-step screening process to isolate shade tolerant mutants can be done efficiently with limited space at early seedling stages, which enables quick and efficient isolation of shade tolerant mutants, and thus facilitates development of shade tolerant new cultivars of turfgrasses. Using the method, we isolated 136 dwarf mutants from 300,000 mutagenized seeds, with 65 being shade tolerant (0.022%). When screened directly for shade tolerance, we recovered only four mutants from a population of 150,000 (0.003%) mutagenized seeds. One shade tolerant mutant, shadow-1, was characterized in detail. In addition to dwarfism, shadow-1 and its sexual progeny displayed high degrees of tolerance to both natural and artificial shade. We showed that endogenous gibberellin (GA) content in shadow-1 was higher than wild-type controls, and shadow-1 was also partially GA insensitive. Our novel, simple and effective two-step screening method should be applicable to breeding shade tolerant cultivars of turfgrasses, ground covers, and other economically important crop plants that can be used under canopies of existing vegetation to increase productivity per unit area of land.

Novel and Recently Evolved MicroRNA Clusters Regulate Expansive F-BOX Gene Networks through Phased Small Interfering RNAs in Wild Diploid Strawberry.

The wild strawberry (Fragaria vesca) has recently emerged as an excellent model for cultivated strawberry (Fragaria × ananassa) as well as other Rosaceae fruit crops due to its short seed-to-fruit cycle, diploidy, and sequenced genome. Deep sequencing and parallel analysis of RNA ends were used to identify F. vesca microRNAs (miRNAs) and their target genes, respectively. Thirty-eight novel and 31 known miRNAs were identified. Many known miRNAs targeted not only conserved mRNA targets but also developed new target genes in F. vesca. Significantly, two new clusters of miRNAs were found to collectively target 94 F-BOX (FBX) genes. One of the miRNAs in the new cluster is 22 nucleotides and triggers phased small interfering RNA production from six FBX genes, which amplifies the silencing to additional FBX genes. Comparative genomics revealed that the main novel miRNA cluster evolved from duplications of FBX genes. Finally, conserved trans-acting siRNA pathways were characterized and confirmed with distinct features. Our work identified novel miRNA-FBX networks in F. vesca and shed light on the evolution of miRNAs/phased small interfering RNA networks that regulate large gene families in higher plants.

Highly interactive nature of flower-specific enhancers and promoters, and its potential impact on tissue-specific expression and engineering of multiple genes or agronomic traits.

Molecular stacking enables multiple traits to be effectively engineered in crops using a single vector. However, the co-existence of distinct plant promoters in the same transgenic unit might, like their mammalian counterparts, interfere with one another. In this study, we devised a novel approach to investigate enhancer-promoter and promoter-promoter interactions in transgenic plants and demonstrated that three of four flower-specific enhancer/promoters were capable of distantly activating a pollen- and stigma-specific Pps promoter (fused to the cytotoxic DT-A gene) in other tissues, as revealed by novel tissue ablation phenotypes in transgenic plants. The NtAGI1 enhancer exclusively activated stamen- and carpel-specific DT-A expression, thus resulting in tissue ablation in an orientation-independent manner; this activation was completely abolished by the insertion of an enhancer-blocking insulator (EXOB) between the NtAGI1 enhancer and Pps promoter. Similarly, AGL8 and AP1Lb1, but not AP1La, promoters also activated distinct tissue-specific DT-A expression and ablation, with the former causing global growth retardation and the latter ablating apical inflorescences. While the tissue specificity of the enhancer/promoters generally defined their activation specificities, the strength of their activity in particular tissues or developmental stages appeared to determine whether activation actually occurred. Our findings provide the first evidence that plant-derived enhancer/promoters can distantly interact/interfere with one another, which could pose potential problems for the tissue-specific engineering of multiple traits using a single-vector stacking approach. Therefore, our work highlights the importance of adopting enhancer-blocking insulators in transformation vectors to minimize promoter-promoter interactions. The practical and fundamental significance of these findings will be discussed.

Conserved transcriptional regulatory programs underlying rice and barley germination.

Germination is a biological process important to plant development and agricultural production. Barley and rice diverged 50 million years ago, but share a similar germination process. To gain insight into the conservation of their underlying gene regulatory programs, we compared transcriptomes of barley and rice at start, middle and end points of germination, and revealed that germination regulated barley and rice genes (BRs) diverged significantly in expression patterns and/or protein sequences. However, BRs with higher protein sequence similarity tended to have more conserved expression patterns. We identified and characterized 316 sets of conserved barley and rice genes (cBRs) with high similarity in both protein sequences and expression patterns, and provided a comprehensive depiction of the transcriptional regulatory program conserved in barley and rice germination at gene, pathway and systems levels. The cBRs encoded proteins involved in a variety of biological pathways and had a wide range of expression patterns. The cBRs encoding key regulatory components in signaling pathways often had diverse expression patterns. Early germination up-regulation of cell wall metabolic pathway and peroxidases, and late germination up-regulation of chromatin structure and remodeling pathways were conserved in both barley and rice. Protein sequence and expression pattern of a gene change quickly if it is not subjected to a functional constraint. Preserving germination-regulated expression patterns and protein sequences of those cBRs for 50 million years strongly suggests that the cBRs are functionally significant and equivalent in germination, and contribute to the ancient characteristics of germination preserved in barley and rice. The functional significance and equivalence of the cBR genes predicted here can serve as a foundation to further characterize their biological functions and facilitate bridging rice and barley germination research with greater confidence.

Systems and evolutionary characterization of microRNAs and their underlying regulatory networks in soybean cotyledons.

MicroRNAs (miRNAs) are an emerging class of small RNAs regulating a wide range of biological processes. Soybean cotyledons evolved as sink tissues to synthesize and store seed reserves which directly affect soybean seed yield and quality. However, little is known about miRNAs and their regulatory networks in soybean cotyledons. We sequenced 292 million small RNA reads expressed in soybean cotyledons, and discovered 130 novel miRNA genes and 72 novel miRNA families. The cotyledon miRNAs arose at various stages of land plant evolution. Evolutionary analysis of the miRNA genes in duplicated genome segments from the recent Glycine whole genome duplication revealed that the majority of novel soybean cotyledon miRNAs were young, and likely arose after the duplication event 13 million years ago. We revealed the evolutionary pathway of a soybean cotyledon miRNA family (soy-miR15/49) that evolved from a neutral invertase gene through an inverted duplication and a series of DNA amplification and deletion events. A total of 304 miRNA genes were expressed in soybean cotyledons. The miRNAs were predicted to target 1910 genes, and form complex miRNA networks regulating a wide range of biological pathways in cotyledons. The comprehensive characterization of the miRNAs and their underlying regulatory networks at gene, pathway and system levels provides a foundation for further studies of miRNAs in cotyledons.