Genome-wide identification of long non-coding RNAs and their potential functions in radish response to salt stress.

Long non-coding RNAs (lncRNAs) are increasingly recognized as cis- and trans-acting regulators of protein-coding genes in plants, particularly in response to abiotic stressors. Among these stressors, high soil salinity poses a significant challenge to crop productivity. Radish (Raphanus sativus L.) is a prominent root vegetable crop that exhibits moderate susceptibility to salt stress, particularly during the seedling stage. Nevertheless, the precise regulatory mechanisms through which lncRNAs contribute to salt response in radish remain largely unexplored. In this study, we performed genome-wide identification of lncRNAs using strand-specific RNA sequencing on radish fleshy root samples subjected to varying time points of salinity treatment. A total of 7,709 novel lncRNAs were identified, with 363 of them displaying significant differential expression in response to salt application. Furthermore, through target gene prediction, 5,006 cis- and 5,983 trans-target genes were obtained for the differentially expressed lncRNAs. The predicted target genes of these salt-responsive lncRNAs exhibited strong associations with various plant defense mechanisms, including signal perception and transduction, transcription regulation, ion homeostasis, osmoregulation, reactive oxygen species scavenging, photosynthesis, phytohormone regulation, and kinase activity. Notably, this study represents the first comprehensive genome-wide analysis of salt-responsive lncRNAs in radish, to the best of our knowledge. These findings provide a basis for future functional analysis of lncRNAs implicated in the defense response of radish against high salinity, which will aid in further understanding the regulatory mechanisms underlying radish response to salt stress.

Influence of the area of the aqueduct and region of interest on quantification of stroke volume in healthy volunteers using phase-contrast cine magnetic resonance imaging.

BACKGROUND: Phase-contrast cine magnetic resonance imaging (PC-MRI) has been used to measure cerebrospinal fluid (CSF) flow dynamics, but the influence of the area of the aqueduct and region of interest (ROI) on quantification of stroke volume (SV) has not been assessed. PURPOSE: To assess the influence of the area of the ROI in quantifying the aqueductal SV measured with PC-MRI within the cerebral aqueduct. MATERIAL AND METHODS: Nine healthy volunteers (mean age = 29.6 years) were enrolled in the study, and brain MRI examinations were performed on a 3.0-T system. Quantitative analysis of the aqueductal CSF flow was performed using manual ROI placement. ROIs were separately drawn for each of the 12 phases of the cardiac cycle, and changes in aqueduct size during the cardiac cycle were determined. The SV was calculated using 12 different aqueductal ROIs and compared with the SV calculated using a fixed ROI size. RESULTS: There was variation in the size of the aqueduct during the cardiac cycle. In addition, the measured SV increased with a greater area of the ROI. A significant difference in the calculated SVs with the 12 variable ROIs was observed compared with that using a fixed ROI throughout the cardiac cycle. CONCLUSION: To establish reliable reference values for the SV in future studies, a variable ROI should be considered.

Chromosome-level pepino genome provides insights into genome evolution and anthocyanin biosynthesis in Solanaceae.

Pepino (Solanum muricatum, 2n = 2x = 24), a member of the Solanaceae family, is an important globally grown fruit. Herein, we report high-quality, chromosome-level pepino genomes. The 91.67% genome sequence is anchored to 12 chromosomes, with a total length of 1.20 Gb and scaffold N50 of 87.03 Mb. More than half the genome comprises repetitive sequences. In addition to the shared ancient whole-genome triplication (WGT) event in eudicots, an additional new WGT event was present in the pepino. Our findings suggest that pepinos experienced chromosome rearrangements, fusions, and gene loss after a WGT event. The large number of gene removals indicated the instability of Solanaceae genomes, providing opportunities for species divergence and natural selection. The paucity of disease-resistance genes (NBS) in pepino and eggplant has been explained by extensive loss and limited generation of genes after WGT events in Solanaceae. The outbreak of NBS genes was not synchronized in Solanaceae species, which occurred before the Solanaceae WGT event in pepino, tomato, and tobacco, whereas it was almost synchronized with WGT events in the other four Solanaceae species. Transcriptome and comparative genomic analyses revealed several key genes involved in anthocyanin biosynthesis. Although an extra WGT event occurred in Solanaceae, CHS genes related to anthocyanin biosynthesis in grapes were still significantly expanded compared with those in Solanaceae species. Proximal and tandem duplications contributed to the expansion of CHS genes. In conclusion, the pepino genome and annotation facilitate further research into important gene functions and comparative genomic analysis in Solanaceae.

Genome-Wide Identification and Expression Profiling of CONSTANS-Like Genes in Pepper (Capsicum annuum): Gaining an Insight to Their Phylogenetic Evolution and Stress-Specific Roles.

CONSTANS-like (COL) genes play important regulatory roles in multiple growth and development processes of plants but have rarely been studied in Capsicum annuum. This study explored the evolutionary relationship and expression patterns of COL genes from C. annuum. A total of 10 COL genes were identified in the genome of the cultivated pepper Zunla-1 and were named CaCOL01-10. These genes were unequally distributed among five chromosomes and could be divided into three groups based on differences in gene structure characteristics. During evolutionary history, duplications and retentions were divergent among different groups of COL genes. Tandem duplication caused amplification of group I genes. Genetic distance among COL genes was the largest in group III, suggesting that group III genes undergo more relaxed selection pressure compared with the other groups. Expression patterns of CaCOLs in tissues were significantly different, with CaCOL08 exhibiting the highest expression in stem and leaf. Some COL orthologous genes showed markedly different expression patterns in pepper compared with tomato, such as COL_1 orthologs, which may be involved in fruit development in pepper. In addition, CaCOLs participated in the regulation of abiotic stresses to varying degrees. Five CaCOL genes were induced by cold, and CaCOL02 and CaCOL03 were specifically upregulated by cold and downregulated by heat. This study provides a theoretical basis for the in-depth understanding of the functions of COL genes in pepper and their molecular mechanisms involved in growth and development and responses to abiotic stresses.

Genome­wide characterization of the Gα subunit gene family in Rosaceae and expression analysis of PbrGPAs under heat stress.

The Gα subunit is an important component of the heterotrimeric G-protein complex and an integral component of several signal transduction pathways. It plays crucial roles in the diverse processes of plant growth and development, including the response to abiotic stress, regulation of root development, involvement in stomatal movement, and participation in hormone responses, which have been well investigated in many species. However, no comprehensive analysis has identified and explored the evolution, expression pattern characteristics and heat stress response of the Gα subunit genes in Rosaceae. In this study, 52 Gα subunit genes were identified in eight Rosaceae species; these genes were divided into three subfamilies (I, II, and III) based on their phylogenetic, conserved motif, and structural characteristics. Whole genome and dispersed duplication events were found to have contributed significantly to the expansion of the Gα subunit gene family, and purifying selection to have played a key role in the evolution of Gα subunit genes. An expression analysis identified some PbrGPA genes that were highly expressed in leaf, root, and fruit, and exhibited diverse spatiotemporal expression models in pear. Under abiotic stress conditions, the mRNA transcript levels of PbrGPA genes were up-regulated in response to high temperature treatment in leaves. Furthermore, three Gα subunit genes were shown to be located in the plasma membrane and nucleus in pear. In conclusion, the study of the Gα subunit gene family will help us to better understand its evolutionary history and expression patterns, while facilitating further investigations into the function of the Gα subunit gene in response to heat stress.

Genome-wide survey of Gγ subunit gene family in eight Rosaceae and expression analysis of PbrGGs in pear (Pyrus bretschneideri).

BACKGROUND: Heterotrimeric G-proteins, composed of Gα, Gß and Gγ subunits, are important signal transmitters, mediating the cellular response to multiple stimuli in animals and plants. The Gγ subunit is an essential component of the G-protein, providing appropriate functional specificity to the heterotrimer complex and has been well studied in many species. However, the evolutionary history, expression pattern and functional characteristics of Gγ subunits has not been explored in the Rosaceae, representing many important fruit crops. RESULTS: In this study, 35 Gγ subunit genes were identified from the eight species belonging to the Rosaceae family. Based on the structural gene characteristics, conserved protein motifs and phylogenetic analysis of the Gγ subunit genes, the genes were classified into three clades. Purifying selection was shown to play an important role in the evolution of Gγ subunit genes, while a recent whole-genome duplication event was the principal force determining the expansion of the Gγ subunit gene family in the subfamily Maloideae. Gγ subunit genes exhibited diverse spatiotemporal expression patterns in Chinese white pear, including fruit, root, ovary and bud, and under abiotic stress conditions, the relative expression of Gγ subunit genes were up-regulated or down-regulated. In addition, seven of the Gγ subunit proteins in pear were located on the plasma membrane, in the cytoplasm or nucleus. CONCLUSION: Overall, this study of the Gγ subunit gene family in eight Rosaceae species provided useful information to better understand the evolution and expression of these genes and facilitated further exploration of their functions in these important crop plants.

Ventriculoperitoneal shunt malfunction diagnosis based on substance dilution.

OBJECTIVE: Current methods for the diagnosis of ventriculoperitoneal (VP) shunt malfunction lack specific standards; therefore, it may be missed or misdiagnosed. Hence, providing a reliable diagnostic method will help improve the accuracy of preoperative decision-making. Therefore, the aim of the study was to provide a new method for the diagnosis of VP shunt malfunction. METHODS: After in vitro testing, we enrolled a total of 12 patients with VP shunt malfunction. Before revision surgery, 0.1 mL of a 5% sodium valproate (SV) solution was injected into the reservoir; 0.1 mL of the cerebrospinal fluid (CSF) was withdrawn 20 minutes later from the reservoir to measure the SV concentration. The process was repeated on the seventh day after surgery and compared with the preoperative results. RESULTS: The mean ±â€Šstandard deviation preoperative SV concentration in the cerebrospinal fluid was greater than the postoperative concentration (5967.8 ±â€Š1281.3 vs 391.1 ±â€Š184.6 µg/mL, P = .001). CONCLUSION: The proposed method is a reliable, safe, and relatively simple alternative for the diagnosis of VP shunt malfunction and further provides a reference for treatment.

Comprehensive analysis of SSRs and database construction using all complete gene-coding sequences in major horticultural and representative plants.

Simple sequence repeats (SSRs) are one of the most important genetic markers and widely exist in most species. Here, we identified 249,822 SSRs from 3,951,919 genes in 112 plants. Then, we conducted a comprehensive analysis of these SSRs and constructed a plant SSR database (PSSRD). Interestingly, more SSRs were found in lower plants than in higher plants, showing that lower plants needed to adapt to early extreme environments. Four specific enriched functional terms in the lower plant Chlamydomonas reinhardtii were detected when it was compared with seven other higher plants. In addition, Guanylate_cyc existed in more genes of lower plants than of higher plants. In our PSSRD, we constructed an interactive plotting function in the chart interface, and users can easily view the detailed information of SSRs. All SSR information, including sequences, primers, and annotations, can be downloaded from our database. Moreover, we developed Web SSR Finder and Batch SSR Finder tools, which can be easily used for identifying SSRs. Our database was developed using PHP, HTML, JavaScript, and MySQL, which are freely available at http://www.pssrd.info/ . We conducted an analysis of the Myb gene families and flowering genes as two applications of the PSSRD. Further analysis indicated that whole-genome duplication and whole-genome triplication played a major role in the expansion of the Myb gene families. These SSR markers in our database will greatly facilitate comparative genomics and functional genomics studies in the future.

Genome-wide survey of sucrose non-fermenting 1-related protein kinase 2 in Rosaceae and expression analysis of PbrSnRK2 in response to ABA stress.

BACKGROUND: The members of the sucrose non-fermenting 1-related protein kinase 2 (SnRK2) family are specific serine/threonine protein kinases in plants that play important roles in stress signal transduction and adaptation. Because of their positive regulatory roles in response to adverse conditions, the genes encoding thes proteins are considered potential candidates for breeding of plants for disease resistance and genetic improvement. However, there is far less information about this kinase family, and the function of these genes has not been explored in Rosaceae. RESULTS: A genome-wide survey and analysis of the genes encoding members of the SnRK2 family were performed in pear (Pyrus bretschneideri) and seven other Rosaceae species. A total of 71 SnRK2 genes were identified from the eight Rosaceae species and classified into three subgroups based on phylogenetic analysis and structural characteristics. Purifying selection played a crucial role in the evolution of SnRK2 genes, and whole-genome duplication and dispersed duplication were the primary forces underlying the characteristics of the SnRK2 gene family in Rosaceae. Transcriptome data and qRT-PCR assay results revealed that the distribution of PbrSnRK2s was very extensive, including across the roots, leaves, pollen, styles, and flowers, although most of them were mainly expressed in leaves. In addition, under stress conditions, the transcript levels of some of the genes were upregulated in leaves in response to ABA treatment. CONCLUSIONS: This study provides useful information and a theoretical introduction for the study of the evolution, expression, and functions of the SnRK2 gene family in plants.

Molecular Evolutionary and Expression Pattern Analysis of AKR Genes Shed New Light on GalUR Functional Characteristics in Brassica rapa.

The aldo-keto reductase (AKR) superfamily plays a major role in oxidation-reduction in plants. D-galacturonic acid reductase (GalUR), an ascorbic acid (AsA) biosynthetic enzyme, belongs to this superfamily. However, the phylogenetic relationship and evolutionary history of the AKR gene family in plants has not yet been clarified. In this study, a total of 1268 AKR genes identified in 36 plant species were used to determine this phylogenetic relationship. The retention, structural characteristics, and expression patterns of AKR homologous genes in Brassica rapa and Arabidopsis thaliana were analyzed to further explore their evolutionary history. We found that the AKRs originated in algae and could be divided into A and B groups according to the bootstrap value; GalURs belonged to group A. Group A AKR genes expanded significantly before the origin of angiosperms. Two groups of AKR genes demonstrated functional divergence due to environmental adaptability, while group A genes were more conservative than those in group B. All 12 candidate GalUR genes were cloned, and their expression patterns under stress were analyzed, in Pak-choi. These genes showed an obvious expression divergence under multiple stresses, and BrcAKR22 exhibited a positive correlation between its expression trend and AsA content. Our findings provide new insights into the evolution of the AKR superfamily and help build a foundation for further investigations of GalUR's functional characteristics.

Evaluation of an artificial intelligent hydrocephalus diagnosis model based on transfer learning.

To design and develop artificial intelligence (AI) hydrocephalus (HYC) imaging diagnostic model using a transfer learning algorithm and evaluate its application in the diagnosis of HYC by non-contrast material-enhanced head computed tomographic (CT) images.A training and validation dataset of non-contrast material-enhanced head CT examinations that comprised of 1000 patients with HYC and 1000 normal people with no HYC accumulating to 28,500 images. Images were pre-processed, and the feature variables were labeled. The feature variables were extracted by the neural network for transfer learning. AI algorithm performance was tested on a separate dataset containing 250 examinations of HYC and 250 of normal. Resident, attending and consultant in the department of radiology were also tested with the test sets, their results were compared with the AI model.Final model performance for HYC showed 93.6% sensitivity (95% confidence interval: 77%, 97%) and 94.4% specificity (95% confidence interval: 79%, 98%), with area under the characteristic curve of 0.93. Accuracy rate of model, resident, attending, and consultant were 94.0%, 93.4%, 95.6%, and 97.0%.AI can effectively identify the characteristics of HYC from CT images of the brain and automatically analyze the images. In the future, AI can provide auxiliary diagnosis of image results and reduce the burden on junior doctors.

Identification and Expression Analysis of Cold Shock Protein 3 (BcCSP3) in Non-Heading Chinese Cabbage (Brassica rapa ssp. chinensis).

A cold-related protein, cold shock protein 3 (BcCSP3), was isolated from non-heading Chinese cabbage in this study. BcCSP3 can encode 205 amino acids (aa) with an open reading frame (ORF) of 618 base pairs (bp). Multiple sequence alignment and phylogenetic tree analyses showed that BcCSP3 contains an N-terminal cold shock domain and is highly similar to AtCSP2, their kinship is recent. Real-time quantitative polymerase chain reaction (RT-qPCR) showed that the expression level of BcCSP3 in stems and leaves is higher than that in roots. Compared with other stress treatments, the change in BcCSP3 expression level was most pronounced under cold stress. In addition, a BcCSP3-GFP fusion protein was localized to the nucleus and cytoplasm. These results indicated that BcCSP3 may play an important role in response to cold stress in non-heading Chinese cabbage. This work may provide a reference for the identification and expression analysis of other CSP genes in non-heading Chinese cabbage.

Enhanced photosynthetic activity in pak choi hybrids is associated with increased grana thylakoids in chloroplasts.

Increased photosynthetic activity is closely linked to heterosis in plants, but the underlying molecular mechanisms remain elusive. Pak choi (Brassica rapa ssp. chinensis) is a widely grown vegetable in Asia, and the most commercial cultivars are F1 hybrids. Here, the inbred pak choi lines WTC and 2Q, and their reciprocal F1 hybrids WQ and QW, were used to characterize the increased photosynthetic activity in these hybrids at the physiological, cellular and molecular levels. We found that the hybrids had larger leaves, with more grana thylakoids. Additionally, these hybrids had significantly increased net photosynthetic rates (Pn ) under both saturating and low irradiance conditions. These data indicate that the increased photosynthetic activity in pak choi hybrids was associated with an improved photosynthetic mechanism and larger leaves. Next, we obtained genome-wide data using transcriptome and bisulfite sequencing. Gene ontology (GO) analysis showed that the differentially expressed genes among the parents and hybrids were mostly enriched in the 'photosynthesis', 'thylakoid', and 'chloroplast' categories, indicating that the increased number of grana thylakoids contributes to the enhanced photosynthetic capacity in hybrids. Furthermore, we found that the increased number of grana thylakoids was associated with the upregulation of light-harvesting complex of photosystem II 1 (BrLhcb1). Yeast one-hybrid and transient assay showed that the BrLhcb1 promoter was directly bound by CIRCADIAN CLOCK ASSOCIATED 1 (BrCCA1), resulting in increased BrLhcb1 expression and enhanced carbon fixation in hybrids. Finally, our findings provide new insight into molecular mechanisms underlying enhanced photosynthesis in pak choi hybrids.

BcMAF2 activates BcTEM1 and represses flowering in Pak-choi (Brassica rapa ssp. chinensis).

KEY MESSAGE: BcMAF2 plays a key role in flowering regulation by controlling BcTEM1, BcSOC1 and BCSPL15 in Pak-choi. Flowering is a key event in the life cycle of plants. Flowering time shows an extensive variation from different Pak-choi (Brassica rapa ssp. chinensis) cultivars. However, the regulation mechanism of flowering in Pak-choi remains rarely known. In this study, a systematic identification and functional analysis of a Pak-choi MADS Affecting Flowering (MAF) gene, BcMAF2, was carried out. BcMAF2 encoded a protein containing a conserved MADS-box domain, which was localized in the nucleus. QPCR analysis indicated that the expression of BcMAF2 was higher in the leaves and flowers. Overexpression of BcMAF2 in Arabidopsis showed that BcMAF2 repressed flowering, which was further confirmed by silencing endogenous BcMAF2 in Pak-choi. In addition, Tempranillo 1 (TEM1) expression was up-regulated and MAF2 expression was down-regulated in the BcMAF2-overexpressing Arabidopsis. The expression of BcMAF2 and BcTEM1 was down-regulated in BcMAF2-silencing Pak-choi plants. The yeast one-hybrid, dual luciferase and qPCR results revealed that BcMAF2 protein could directly bind to BcTEM1 promoter and activate its expression, which was not reported in Arabidopsis. Meanwhile, a self-inhibition was found in BcMAF2. Taken together, this work suggested that BcMAF2 could repress flowering by directly activating BcTEM1.

Genome-Wide Identification, Evolution, and Expression Analysis of the ATP-Binding Cassette Transporter Gene Family in Brassica rapa.

ATP-binding cassette (ABC) proteins can act as transporters of different substrates across biological membranes by hydrolyzing ATP. However, little information is available about ABC transporters in Brassica rapa, an important leafy vegetable. In the present study, we carried out genome-wide identification, characterization and molecular evolution analyses of ABC gene family in B. rapa and 9 other plant species. A total of 179 B. rapa ABC genes (BraABCs) were identified. Among them, 173 BraABCs were identified on 10 chromosomes. Based on phylogenetic analysis and domain organization, the BraABC family could be grouped into eight subfamilies. BraABCs in the same subfamily showed similar motif composition and exon-intron organization. Common and unique cis-elements involved in the transcriptional regulation were also identified in the promoter regions of BraABCs. Tissue-expression analysis of BraABCs demonstrated their diverse spatiotemporal expression profiles. Influences of the whole genome triplication (WGT) on the evolution of BraABCs were studied in detail. BraABCs were preferentially retained compared with their neighboring genes during diploidization after WGT. Synteny analysis identified 76 pairs of syntenic BraABC paralogs among the three subgenomes of B. rapa, and 10 paralog pairs underwent positive selection with ω (= Ka/Ks) ratios greater than 1. Analyses of the expression patterns of syntenic BraABC paralogs pairs across five tissues and under stress treatments revealed their functional conservation, sub-functionalization, neo-functionalization and pseudogenization during evolution. Our study presents a comprehensive overview of the ABC gene family in B. rapa and will be helpful for the further functional study of BraABCs in plant growth, development, and stress responses.

Comprehensive Analysis of the CDPK-SnRK Superfamily Genes in Chinese Cabbage and Its Evolutionary Implications in Plants.

The CDPK-SnRK (calcium-dependent protein kinase/Snf1-related protein kinase) gene superfamily plays important roles in signaling pathways for disease resistance and various stress responses, as indicated by emerging evidence. In this study, we constructed comparative analyses of gene structure, retention, expansion, whole-genome duplication (WGD) and expression patterns of CDPK-SnRK genes in Brassica rapa and their evolution in plants. A total of 49 BrCPKs, 14 BrCRKs, 3 BrPPCKs, 5 BrPEPRKs, and 56 BrSnRKs were identified in B. rapa. All BrCDPK-SnRK proteins had highly conserved kinase domains. By statistical analysis of the number of CDPK-SnRK genes in each species, we found that the expansion of the CDPK-SnRK gene family started from angiosperms. Segmental duplication played a predominant role in CDPK-SnRK gene expansion. The analysis showed that PEPRK was more preferentially retained than other subfamilies and that CPK was retained similarly to SnRK. Among the CPKs and SnRKs, CPKIII and SnRK1 genes were more preferentially retained than other groups. CRK was closest to CPK, which may share a common evolutionary origin. In addition, we identified 196 CPK genes and 252 SnRK genes in 6 species, and their different expansion and evolution types were discovered. Furthermore, the expression of BrCDPK-SnRK genes is dynamic in different tissues as well as in response to abiotic stresses, demonstrating their important roles in development in B. rapa. In summary, this study provides genome-wide insight into the evolutionary history and mechanisms of CDPK-SnRK genes following whole-genome triplication in B. rapa.

Cold acclimation alters DNA methylation patterns and confers tolerance to heat and increases growth rate in Brassica rapa.

Epigenetic modifications are implicated in plant adaptations to abiotic stresses. Exposure of plants to one stress can induce resistance to other stresses, a process termed cross-adaptation, which is not well understood. In this study, we aimed to unravel the epigenetic basis of elevated heat-tolerance in cold-acclimated Brassica rapa by conducting a genome-wide DNA methylation analysis of leaves from control (CK) and cold-acclimated (CA) plants. We found that both methylation and demethylation occurred during cold acclimation. Two significantly altered pathways, malate dehydrogenase activity and carbon fixation, and 1562 differentially methylated genes, including BramMDH1, BraKAT2, BraSHM4, and Bra4CL2, were identified in CA plants. Genetic validation and treatment of B. rapa with 5-aza-2-deoxycytidine (Aza) suggested that promoter demethylation of four candidate genes increased their transcriptional activities. Physiological analysis suggested that elevated heat-tolerance and high growth rate were closely related to increases in organic acids and photosynthesis, respectively. Functional analyses demonstrated that the candidate gene BramMDH1 (mMDH: mitochondrial malate dehydrogenase) directly enhances organic acids and photosynthesis to increase heat-tolerance and growth rate in Arabidopsis. However, Aza-treated B. rapa, which also has elevated BramMDH1 levels, did not exhibit enhanced heat-tolerance. We therefore suggest that DNA demethylation alone is not sufficient to increase heat-tolerance. This study demonstrates that altered DNA methylation contributes to cross-adaptation.

Role of vernalization-mediated demethylation in the floral transition of Brassica rapa.

MAIN CONCLUSION: Vernalization-mediated demethylation of BrCKA2 (casein kinase II α-subunit) and BrCKB4 (casein kinase II ß-subunit) shorten the period of the clock gene BrCCA1 (circadian clock associated 1) in Brassica rapa. Photoperiod and vernalization are two environmental cues involved in the regulation of floral transition, but the ways in which they interact remain unclear. DNA methylation is one of the main mechanisms involved in controlling the functional state of chromatin and gene expression in response to environmental signals. To study the interaction between photoperiod and vernalization in floral transition, we carried out a comparative genomic analysis of genome-wide DNA methylation profiles in normal (CK) and vernalized (CA) leaves from Brassica rapa using methylated-DNA immunoprecipitation sequencing (MeDIP-seq). Two subunits of casein kinase II (CK2), BrCKA2 (catalytic α-subunit of CK2) and BrCKB4 (regulatory ß-subunit of CK2), exhibited gradual DNA demethylation and increased expression in vernalized B. rapa. DNA methylation-defective plants demonstrated the causal link between DNA demethylation changes and changes in gene expression. Virus-induced gene silencing (VIGS) of BrCKA2 and BrCKB4 in B. rapa resulted in no change to the period of BrCCA1 (circadian clock associated 1) and a 1-week late flowering time. Finally, we demonstrated that increased levels of BrCKA2 and BrCKB4 in vernalized B. rapa confer elevated CK2 activity, resulting in a shortened period of the clock gene BrCCA1, which plays an important role in perceiving photoperiod in plants. Thus, our results suggest that there is a direct interaction between photoperiod and vernalization through DNA methylation mechanisms.

Conservation and Expression Patterns Divergence of Ascorbic Acid d-mannose/l-galactose Pathway Genes in Brassica rapa.

Ascorbic acid (AsA) participates in diverse biological processes, is regulated by multiple factors and is a potent antioxidant and cellular reductant. The D-Mannose/L-Galactose pathway is a major plant AsA biosynthetic pathway that is highly connected within biosynthetic networks, and generally conserved across plants. Previous work has shown that, although most genes of this pathway are expressed under standard growth conditions in Brassica rapa, some paralogs of these genes are not. We hypothesize that regulatory evolution in duplicate AsA pathway genes has occurred as an adaptation to environmental stressors, and that gene retention has been influenced by polyploidation events in Brassicas. To test these hypotheses, we explored the conservation of these genes in Brassicas and their expression patterns divergence in B. rapa. Similar retention and a high degree of gene sequence similarity were identified in B. rapa (A genome), B. oleracea (C genome) and B. napus (AC genome). However, the number of genes that encode the same type of enzymes varied among the three plant species. With the exception of GMP, which has nine genes, there were one to four genes that encoded the other enzymes. Moreover, we found that expression patterns divergence widely exists among these genes. (i) VTC2 and VTC5 are paralogous genes, but only VTC5 is influenced by FLC. (ii) Under light treatment, PMI1 co-regulates the AsA pool size with other D-Man/L-Gal pathway genes, whereas PMI2 is regulated only by darkness. (iii) Under NaCl, Cu(2+), MeJA and wounding stresses, most of the paralogs exhibit different expression patterns. Additionally, GME and GPP are the key regulatory enzymes that limit AsA biosynthesis in response to these treatments. In conclusion, our data support that the conservative and divergent expression patterns of D-Man/L-Gal pathway genes not only avoid AsA biosynthesis network instability but also allow B. rapa to better adapt to complex environments.