Genome-Wide Analysis and Expression Profiling of the Glutathione Peroxidase-like Enzyme Gene Family in Solanum tuberosum.

Glutathione peroxidase-like enzyme is an important enzymatic antioxidant in plants. It is involved in scavenging reactive oxygen species, which can effectively prevent oxidative damage and improve resistance. GPXL has been studied in many plants but has not been reported in potatoes, the world's fourth-largest food crop. This study identified eight StGPXL genes in potatoes for the first time through genome-wide bioinformatics analysis and further studied the expression patterns of these genes using qRT-PCR. The results showed that the expression of StGPXL1 was significantly upregulated under high-temperature stress, indicating its involvement in potato defense against high-temperature stress, while the expression levels of StGPXL4 and StGPXL5 were significantly downregulated. The expression of StGPXL1, StGPXL2, StGPXL3, and StGPXL6 was significantly upregulated under drought stress, indicating their involvement in potato defense against drought stress. After MeJA hormone treatment, the expression level of StGPXL6 was significantly upregulated, indicating its involvement in the chemical defense mechanism of potatoes. The expression of all StGPXL genes is inhibited under biotic stress, which indicates that GPXL is a multifunctional gene family, which may endow plants with resistance to various stresses. This study will help deepen the understanding of the function of the potato GPXL gene family, provide comprehensive information for the further analysis of the molecular function of the potato GPXL gene family as well as a theoretical basis for potato molecular breeding.

In-silico genome wide analysis of Mitogen activated protein kinase kinase kinase gene family in C. sinensis.

Mitogen activated protein kinase kinase kinase (MAPKKK) form the upstream component of MAPK cascade. It is well characterized in several plants such as Arabidopsis and rice however the knowledge about MAPKKKs in tea plant is largely unknown. In the present study, MAPKKK genes of tea were obtained through a genome wide search using Arabidopsis thaliana as the reference genome. Among 59 candidate MAPKKK genes in tea, 17 genes were MEKK-like, 31 genes were Raf-like and 11 genes were ZIK- like. Additionally, phylogenetic relationships were established along with structural analysis, which includes gene structure, its location as well as conserved motifs, cis-acting regulatory elements and functional domain signatures that were systematically examined. Also, on the basis of one orthologous gene found between tea and Arabidopsis, functional interaction was carried out in C. sinensis based on an Arabidopsis association model. The expressional profiles indicated major involvement of MAPKKK genes from tea in response to various abiotic stress factors. Taken together, this study provides the targets for additional inclusive identification, functional study, and provides comprehensive knowledge for a better understanding of the MAPKKK cascade regulatory network in C. sinensis.

OPT gene family analysis of potato (Solanum tuberosum) responding to heavy metal stress: Comparative omics and co-expression networks revealed the underlying core templates and specific response patterns.

Oligopeptides transporter (OPT) can maintain intracellular metal homeostat, however, their evolutionary characteristics, as well as their expression patterns in heavy metal exposure, remain unclear. Compared with previous OPT family identification, we identified 94 OPT genes (including 21 in potato) in potato and 4 other plants by HMMER program based on OPT domain (PF03169) for the first time. Secondly, conserved and special OPTs were found through comprehensive analysis. Thirdly, spatio-temporal tissue specific expression patterns and co-expression frameworks of potato OPT genes under different heavy metal stress were constructed. These data can provide excellent gene resources for food security and soil remediation.

Customized optical mapping by CRISPR-Cas9 mediated DNA labeling with multiple sgRNAs.

Whole-genome mapping technologies have been developed as a complementary tool to provide scaffolds for genome assembly and structural variation analysis (1,2). We recently introduced a novel DNA labeling strategy based on a CRISPR-Cas9 genome editing system, which can target any 20bp sequences. The labeling strategy is specifically useful in targeting repetitive sequences, and sequences not accessible to other labeling methods. In this report, we present customized mapping strategies that extend the applications of CRISPR-Cas9 DNA labeling. We first design a CRISPR-Cas9 labeling strategy to interrogate and differentiate the single allele differences in NGG protospacer adjacent motifs (PAM sequence). Combined with sequence motif labeling, we can pinpoint the single-base differences in highly conserved sequences. In the second strategy, we design mapping patterns across a genome by selecting sets of specific single-guide RNAs (sgRNAs) for labeling multiple loci of a genomic region or a whole genome. By developing and optimizing a single tube synthesis of multiple sgRNAs, we demonstrate the utility of CRISPR-Cas9 mapping with 162 sgRNAs targeting the 2Mb Haemophilus influenzae chromosome. These CRISPR-Cas9 mapping approaches could be particularly useful for applications in defining long-distance haplotypes and pinpointing the breakpoints in large structural variants in complex genomes and microbial mixtures.

Evolution of a plant gene cluster in Solanaceae and emergence of metabolic diversity.

Plants produce phylogenetically and spatially restricted, as well as structurally diverse specialized metabolites via multistep metabolic pathways. Hallmarks of specialized metabolic evolution include enzymatic promiscuity and recruitment of primary metabolic enzymes and examples of genomic clustering of pathway genes. Solanaceae glandular trichomes produce defensive acylsugars, with sidechains that vary in length across the family. We describe a tomato gene cluster on chromosome 7 involved in medium chain acylsugar accumulation due to trichome specific acyl-CoA synthetase and enoyl-CoA hydratase genes. This cluster co-localizes with a tomato steroidal alkaloid gene cluster and is syntenic to a chromosome 12 region containing another acylsugar pathway gene. We reconstructed the evolutionary events leading to this gene cluster and found that its phylogenetic distribution correlates with medium chain acylsugar accumulation across the Solanaceae. This work reveals insights into the dynamics behind gene cluster evolution and cell-type specific metabolite diversity.

Plants produce a vast variety of different molecules known as secondary or specialized metabolites to attract pollinating insects, such as bees, or protect themselves against herbivores and pests. The secondary metabolites are made from simple building blocks that are readily available in plants, including amino acids, fatty acids and sugars. Different species of plant, and even different parts of the same plant, produce their own sets of secondary metabolites. For example, the hairs on the surface of tomatoes and other members of the nightshade family of plants make metabolites known as acylsugars. These chemicals deter herbivores and pests from damaging the plants. To make acylsugars, the plants attach long chains known as fatty acyl groups to molecules of sugar, such as sucrose. Some members of the nightshade family produce acylsugars with longer chains than others. In particular, acylsugars with long chains are only found in tomatoes and other closely-related species. It remained unclear how the nightshade family evolved to produce acylsugars with chains of different lengths. To address this question, Fan et al. used genetic and biochemical approaches to study tomato plants and other members of the nightshade family. The experiments identified two genes known as AACS and AECH in tomatoes that produce acylsugars with long chains. These two genes originated from the genes of older enzymes that metabolize fatty acids ­ the building blocks of fats ­ in plant cells. Unlike the older genes, AACS and AECH were only active at the tips of the hairs on the plant's surface. Fan et al. then investigated the evolutionary relationship between 11 members of the nightshade family and two other plant species. This revealed that AACS and AECH emerged in the nightshade family around the same time that longer chains of acylsugars started appearing. These findings provide insights into how plants evolved to be able to produce a variety of secondary metabolites that may protect them from a broader range of pests. The gene cluster identified in this work could be used to engineer other species of crop plants to start producing acylsugars as natural pesticides.

A Methionine Sulfoxide Reductase B Is Required for the Establishment of Astragalus sinicus-Mesorhizobium Symbiosis.

Methionine sulfoxide reductase B (MsrB) is involved in oxidative stress or defense responses in plants. However, little is known about its role in legume-rhizobium symbiosis. In this study, an MsrB gene was identified from Astragalus sinicus and its function in symbiosis was characterized. AsMsrB was induced under phosphorus starvation and displayed different expression patterns under symbiotic and nonsymbiotic conditions. Hydrogen peroxide or methyl viologen treatment enhanced the transcript level of AsMsrB in roots and nodules. Subcellular localization showed that AsMsrB was localized in the cytoplasm of onion epidermal cells and co-localized with rhizobia in nodules. Plants with AsMsrB-RNAi hairy roots exhibited significant decreases in nodule number, nodule nitrogenase activity and fresh weight of the aerial part, as well as an abnormal nodule and symbiosome development. Statistical analysis of infection events showed that plants with AsMsrB-RNAi hairy roots had significant decreases in the number of root hair curling events, infection threads and nodule primordia compared with the control. The content of hydrogen peroxide increased in AsMsrB-RNAi roots but decreased in AsMsrB overexpression roots at the early stage of infection. The transcriptome analysis showed synergistic modulations of the expression of genes involved in reactive oxygen species generation and scavenging, defense and pathogenesis and early nodulation. In addition, a candidate protein interacting with AsMsrB was identified and confirmed by bimolecular fluorescence complementation. Taken together, our results indicate that AsMsrB plays an essential role in nodule development and symbiotic nitrogen fixation by affecting the redox homeostasis in roots and nodules.

Genome-wide characterization of tea plant (Camellia sinensis) Hsf transcription factor family and role of CsHsfA2 in heat tolerance.

BACKGROUND: Heat stress factors (Hsfs) play vital roles in signal transduction pathways operating in responses to environmental stresses. However, Hsf gene family has not been thoroughly explored in tea plant (Camellia sinensis L.). RESULTS: In this study, we identified 25 CsHsf genes in C. sinensis that were separated by phylogenetic analysis into three sub-families (i.e., A, B, and C). Gene structures, conserved domains and motifs analyses indicated that the CsHsf members in each class were relatively conserved. Various cis-acting elements involved in plant growth regulation, hormone responses, stress responses, and light responses were located in the promoter regions of CsHsfs. Furthermore, degradome sequencing analysis revealed that 7 CsHsfs could be targeted by 9 miRNAs. The expression pattern of each CsHsf gene was significantly different in eight tissues. Many CsHsfs were differentially regulated by drought, salt, and heat stresses, as well as exogenous abscisic acid (ABA) and Ca2+. In addition, CsHsfA2 was located in the nucleus. Heterologous expression of CsHsfA2 improved thermotolerance in transgenic yeast, suggesting its potential role in the regulation of heat stress response. CONCLUSIONS: A comprehensive genome-wide analysis of Hsf in C. sinensis present the global identification and functional prediction of CsHsfs. Most of them were implicated in a complex gene regulatory network controlling various abiotic stress responses and signal transduction pathways in tea plants. Additionally, heterologous expression of CsHsfA2 increased thermotolerance of transgenic yeast. These findings provide new insights into the functional divergence of CsHsfs and a basis for further research on CsHsfs functions.

A GDSL-type esterase/lipase gene, GELP77, is necessary for pollen dissociation and fertility in Arabidopsis.

Pollen wall characteristics are dramatically changed during pollen maturation. Many genes have been identified as regulators of such changes in pollen wall characteristics, but mechanisms of such changes have not been completely understood. Here, a GDSL-type esterase/lipase gene, GELP77, is shown to regulate such changes in Arabidopsis thaliana. GELP77-deficient (gelp77) plants exhibited male sterility, and this phenotype was suppressed by introduction of a GELP77 genomic fragment. Mature pollen grains of wild-type Arabidopsis plants have an organized reticulate surface structure and are dissociated from each other. In contrast, pollen grains of gelp77 lacked such a structure and were shrunken and stuck to each other. Nuclei were not detectable in gelp77 microspores at a putative uninucleate stage, suggesting that GELP77 is required as early as this stage. In plants that have the GELP77 promoter-GELP77-GFP transgene, the GELP77-GFP fusion protein was detected in microspores, tapetal cells and middle layer cells in anthers at post-meiotic stages, whereas not anthers at pre-meiotic stages. Analysis of amino acid sequences suggests that GELP77 is phylogenetically distant from the other 104 GDSL-type esterase/lipase genes in Arabidopsis and that GELP77 orthologs are present in various plant species. Together, these results indicate that GELP77 regulates pollen wall characteristics in Arabidopsis.

Genome-wide investigation of superoxide dismutase (SOD) gene family and their regulatory miRNAs reveal the involvement in abiotic stress and hormone response in tea plant (Camellia sinensis).

Superoxide dismutases (SODs), as a family of metalloenzymes related to the removal of reactive oxygen species (ROS), have not previously been investigated at genome-wide level in tea plant. In this study, 10 CsSOD genes were identified in tea plant genome, including 7 Cu/Zn-SODs (CSDs), 2 Fe-SODs (FSDs) and one Mn-SOD (MSD), and phylogenetically classified in three subgroups, respectively. Physico-chemical characteristic, conserved motifs and potential protein interaction analyses about CsSOD proteins were carried out. Exon-intron structures and codon usage bias about CsSOD genes were also examined. Exon-intron structures analysis revealed that different CsSOD genes contained various number of introns. On the basis of the prediction of regulatory miRNAs of CsSODs, a modification 5' RNA ligase-mediated (RLM)-RACE was performed and validated that csn-miR398a-3p-1 directly cleaves CsCSD4. By prediction of cis-acting elements, the expression patterns of 10 CsSOD genes and their regulatory miRNAs were detected under cold, drought, exogenous methyl jasmonate (MeJA) and gibberellin (GA3) treatments. The results showed that most of CsSODs except for CsFSD2 were induced under cold stress and CsCSDs may play primary roles under drought stress; exogenous GA3 and MeJA could also stimulated/inhibited distinct CsSODs at different stages. In addition, we found that csn-miR398a-3p-1 negatively regulated the expression of CsCSD4 may be a crucial regulatory mechanism under cold stress. This study provides a certain basis for the studies about stress resistance in tea plants, even provide insight into comprehending the classification, evolution, diverse functions and influencing factors of expression patterns for CsSOD genes.

Genome-wide identification, characterization, and expression patterns of the BZR transcription factor family in sugar beet (Beta vulgaris L.).

BACKGROUND: BRASSINAZOLE-RESISTANT (BZR) family genes encode plant-specific transcription factors (TFs) that participate in brassinosteroid signal transduction. BZR TFs have vital roles in plant growth, including cell elongation. However, little is known about BZR genes in sugar beet (Beta vulgaris L.). RESULTS: Therefore, we performed a genome-wide investigation of BvBZR genes in sugar beet. Through an analysis of the BES1_N conserved domain, six BvBZR gene family members were identified in the sugar beet genome, which clustered into three subgroups according to a phylogenetic analysis. Each clade was well defined by the conserved motifs, implying that close genetic relationships could be identified among the members of each subfamily. According to chromosomal distribution mapping, 2, 1, 1, 1, and 1 genes were located on chromosomes 1, 4, 5, 6, and 8, respectively. The cis-acting elements related to taproot growth were randomly distributed in the promoter sequences of the BvBZR genes. Tissue-specific expression analyses indicated that all BvBZR genes were expressed in all three major tissue types (roots, stems, and leaves), with significantly higher expression in leaves. Subcellular localization analysis revealed that Bv1_fxre and Bv6_nyuw are localized in the nuclei, consistent with the prediction of Wolf PSORT. CONCLUSION: These findings offer a basis to predict the functions of BZR genes in sugar beet, and lay a foundation for further research of the biological functions of BZR genes in sugar beet.

Nondisjunction and unequal spindle organization accompany the drive of Aegilops speltoides B chromosomes.

B chromosomes (Bs) are supernumerary chromosomes, which are often preferentially inherited. When transmission rates of chromosomes are higher than 0.5, not obeying the Mendelian law of equal segregation, the resulting transmission advantage is collectively referred to as 'chromosome drive'. Here we analysed the drive mechanism of Aegilops speltoides Bs. The repeat AesTR-183 of A. speltoides Bs, which also can be detected on the Bs of Aegilops mutica and rye, was used to track Bs during pollen development. Nondisjunction of CENH3-positive, tubulin interacting B sister chromatids and an asymmetric spindle during first pollen grain mitosis are key for the accumulation process. A quantitative flow cytometric approach revealed that, independent of the number of Bs present in the mother plant, Bs accumulate in the generative nuclei to > 93%. Nine out of 11 tested (peri)centromeric repeats were shared by A and B chromosomes. Our findings provide new insights into the process of chromosome drive. Quantitative flow cytometry is a useful and reliable method to study the drive frequency of Bs. Nondisjunction and unequal spindle organization accompany during first pollen mitosis the drive of A. speltoides Bs. The prerequisites for the drive process seems to be common in Poaceae.

Genome-wide analysis of the MATE gene family in potato.

The multidrug and toxic compound extrusion (MATE) protein family is a newly discovered family of secondary transporters that extrude metabolic waste and a variety of antibiotics out of the cell using an electrochemical gradient of H+ or Na+ across the membrane. The main function of MATE gene family is to participate in the process of plant detoxification and morphogenesis. The genome-wide analysis of the MATE genes in potato genome was conducted. At least 48 genes were initially identified and classified into six subfamilies. The chromosomal localization of MATE gene family showed that they could be distributed on 11 chromosomes except chromosome 9. The number of amino acids is 145-616, the molecular weight of proteins is 15.96-66.13 KD, the isoelectric point is 4.97-9.17, and they were located on the endoplasmic reticulum with having 4-13 transmembrane segments. They contain only two parts of the exons and UTR without introns. Some members of the first subfamily of potato MATE gene family are clustered with At2g04070 and they may be related to the transport of toxic compounds such as alkaloids and heavy metal. The function of the members of the second subfamily may be similar to that of At3g23560, which is related to tetramethylammonium transport. Some members of the third subfamily are clustered with At3g59030 and they may be involved in the transport of flavonoids. The fifth subfamily may be related to the transport of iron ions. The function of the sixth subfamily may be similar to that of At4g39030, which is related to salicylic acid transport. There are three kinds of conserved motifs in potato MATE genes, including the motif 1, motif 2, and motif 3. Each motif has 50 amino acids. The number of each motif is different in the gene sequence, of which 45 MATE genes contain at least a motif, but there is no motif in ST0015301, ST0045283, and ST0082336. These results provide a reference for further research on the function of potato MATE genes.

Genome-wide identification and characterization of the CsSnRK2 family in Camellia sinensis.

The sucrose nonfermenting 1 (SNF1)-related protein kinase 2 (SnRK2) genes play central roles in plant stress signal transduction. In this study, 8 SnRK2 genes were identified from the tea plant genome database and named CsSnRK2.1-8. Phylogenetic analysis showed that the CsSnRK2 genes were classifiable into three groups, similar to those of Arabidopsis thaliana, Oryza sativa and maize. The coding sequences (CDSs) of all CsSnRK2s were separated by eight introns, and their exon-intron organizations exhibited high similarity to those of other plants. The fluorescence of GFP fused with CsSnRK2.3 was detected in only the cytoplasm, while the rest of the proteins showed GFP signal in both the nucleus and the cytoplasm. The results of the expression patterns of the CsSnRK2 genes showed that CsSnRK2s were differentially induced by salt, polyethylene glycol (PEG) and abscisic acid (ABA) stress. Interestingly, The expression of CsSnRK2.3 was inhibited by ABA, suggesting the complicated roles of CsSnRK2s in the ABA signal transduction pathway. Some CsSnRK2 gene pairs showed significant expression change correlations under stresses, indicating that CsSnRK2s might exhibit synergistic effects of signal regulation in response to various stresses. In summary, this comprehensive analysis will facilitate further studies of the SnRK2 family of Camellia sinensis and provide useful information for the functional validation of CsSnRK2s.

Genome-wide identification and characterization of ALTERNATIVE OXIDASE genes and their response under abiotic stresses in Camellia sinensis (L.) O. Kuntze.

MAIN CONCLUSION: Four typical ALTERNATIVE OXIDASE genes have been identified in tea plants, and their sequence features and gene expression profiles have provided useful information for further studies on function and regulation. Alternative oxidase (AOX) is a terminal oxidase located in the respiratory electron transport chain. AOX catalyzes the oxidation of quinol and the reduction of oxygen into water. In this study, a genome-wide search and subsequent DNA cloning were performed to identify and characterize AOX genes in tea plant (Camellia sinensis (L.) O. Kuntze cv. Longjing43). Our results showed that tea plant possesses four AOX genes, i.e., CsAOX1a, CsAOX1d, CsAOX2a and CsAOX2b. Gene structure and protein sequence analyses revealed that all CsAOXs share a four-exon/three-intron structure with highly conserved regions and amino acid residues, which are necessary for AOX secondary structures, catalytic activities and post-translational regulations. All CsAOX were shown to localize in mitochondria using the green fluorescent protein (GFP)-targeting assay. Both CsAOX1a and CsAOX1d were induced by cold, salt and drought stresses, and with different expression patterns in young and mature leaves. Reactive oxygen species (ROS) accumulated strongly after 72 and 96 h cold treatments in both young and mature leaves, while the polyphenol and total catechin decreased significantly only in mature leaves. In comparison to AtAOX1a in Arabidopsis thaliana, CsAOX1a lost almost all of the stress-responsive cis-acting regulatory elements in its promoter region (1500 bp upstream), but possesses a flavonoid biosynthesis-related MBSII cis-acting regulatory element. These results suggest a link between CsAOX1a function and the metabolism of some secondary metabolites in tea plant. Our studies provide a basis for the further elucidation of the biological function and regulation of the AOX pathway in tea plants.

Systematic analysis of DEMETER-like DNA glycosylase genes shows lineage-specific Smi-miR7972 involved in SmDML1 regulation in Salvia miltiorrhiza.

DEMETER-like DNA glycosylases (DMLs) initiate the base excision repair-dependent DNA demethylation to regulate a wide range of biological processes in plants. Six putative SmDML genes, termed SmDML1-SmDML6, were identified from the genome of S. miltiorrhiza, an emerging model plant for Traditional Chinese Medicine (TCM) studies. Integrated analysis of gene structures, sequence features, conserved domains and motifs, phylogenetic analysis and differential expression showed the conservation and divergence of SmDMLs. SmDML1, SmDML2 and SmDML4 were significantly down-regulated by the treatment of 5Aza-dC, a general DNA methylation inhibitor, suggesting involvement of SmDMLs in genome DNA methylation change. SmDML1 was predicted and experimentally validated to be target of Smi-miR7972. Computational analysis of forty whole genome sequences and almost all of RNA-seq data from Lamiids revealed that MIR7972s were only distributed in some plants of the three orders, including Lamiales, Solanales and Boraginales, and the number of MIR7972 genes varied among species. It suggests that MIR7972 genes underwent expansion and loss during the evolution of some Lamiids species. Phylogenetic analysis of MIR7972s showed closer evolutionary relationships between MIR7972s in Boraginales and Solanales in comparison with Lamiales. These results provide a valuable resource for elucidating DNA demethylation mechanism in S. miltiorrhiza.

Genome-wide identification of conserved and novel microRNAs in one bud and two tender leaves of tea plant (Camellia sinensis) by small RNA sequencing, microarray-based hybridization and genome survey scaffold sequences.

BACKGROUND: MicroRNAs (miRNAs) are important for plant growth and responses to environmental stresses via post-transcriptional regulation of gene expression. Tea, which is primarily produced from one bud and two tender leaves of the tea plant (Camellia sinensis), is one of the most popular non-alcoholic beverages worldwide owing to its abundance of secondary metabolites. A large number of miRNAs have been identified in various plants, including non-model species. However, due to the lack of reference genome sequences and/or information of tea plant genome survey scaffold sequences, discovery of miRNAs has been limited in C. sinensis. RESULTS: Using small RNA sequencing, combined with our recently obtained genome survey data, we have identified and analyzed 175 conserved and 83 novel miRNAs mainly in one bud and two tender leaves of the tea plant. Among these, 93 conserved and 18 novel miRNAs were validated using miRNA microarray hybridization. In addition, the expression pattern of 11 conserved and 8 novel miRNAs were validated by stem-loop-qRT-PCR. A total of 716 potential target genes of identified miRNAs were predicted. Further, Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that most of the target genes were primarily involved in stress response and enzymes related to phenylpropanoid biosynthesis. The predicted targets of 4 conserved miRNAs were further validated by 5'RLM-RACE. A negative correlation between expression profiles of 3 out of 4 conserved miRNAs (csn-miR160a-5p, csn-miR164a, csn-miR828 and csn-miR858a) and their targets (ARF17, NAC100, WER and MYB12 transcription factor) were observed. CONCLUSION: In summary, the present study is one of few such studies on miRNA detection and identification in the tea plant. The predicted target genes of majority of miRNAs encoded enzymes, transcription factors, and functional proteins. The miRNA-target transcription factor gene interactions may provide important clues about the regulatory mechanism of these miRNAs in the tea plant. The data reported in this study will make a huge contribution to knowledge on the potential miRNA regulators of the secondary metabolism pathway and other important biological processes in C. sinensis.

Characterization of the polyphenol oxidase gene family reveals a novel microRNA involved in posttranscriptional regulation of PPOs in Salvia miltiorrhiza.

Salvia miltiorrhiza is a well-known material of traditional Chinese medicine. Understanding the regulatory mechanisms of phenolic acid biosynthesis and metabolism are important for S. miltiorrhiza quality improvement. We report here that S. miltiorrhiza contains 19 polyphenol oxidases (PPOs), forming the largest PPO gene family in plant species to our knowledge. Analysis of gene structures and sequence features revealed the conservation and divergence of SmPPOs. SmPPOs were differentially expressed in plant tissues and eight of them were predominantly expressed in phloem and xylem, indicating that some SmPPOs are functionally redundant, whereas the others are associated with different physiological processes. Expression patterns of eighteen SmPPOs were significantly altered under MeJA treatment, and twelve were yeast extract and Ag+-responsive, suggesting the majority of SmPPOs are stress-responsive. Analysis of high-throughput small RNA sequences and degradome data showed that miR1444-mediated regulation of PPOs existing in P. trichocarpa is absent from S. miltiorrhiza. Instead, a subset of SmPPOs was posttranscriptionally regulated by a novel miRNA, termed Smi-miR12112. It indicates the specificity and significance of miRNA-mediated regulation of PPOs. The results shed light on the regulation of SmPPO expression and suggest the complexity of SmPPO-associated phenolic acid biosynthesis and metabolism.

Identification of Symmetrical RNA Editing Events in the Mitochondria of Salvia miltiorrhiza by Strand-specific RNA Sequencing.

Salvia miltiorrhiza is one of the most widely-used medicinal plants. Here, we systematically analyzed the RNA editing events in its mitochondria. We developed a pipeline using REDItools to predict RNA editing events from stand-specific RNA-Seq data. The predictions were validated using reverse transcription, RT-PCR amplification and Sanger sequencing experiments. Putative sequences motifs were characterized. Comparative analyses were carried out between S. miltiorrhiza, Arabidopsis thaliana and Oryza sativa. We discovered 1123 editing sites, including 225 "C to U" sites in the protein-coding regions. Fourteen of sixteen (87.5%) sites were validated. Three putative DNA motifs were identified around the predicted sites. The nucleotides on both strands at 115 of the 225 sites had undergone RNA editing, which we called symmetrical RNA editing (SRE). Four of six these SRE sites (66.7%) were experimentally confirmed. Re-examination of strand-specific RNA-Seq data from A. thaliana and O. sativa identified 327 and 369 SRE sites respectively. 78, 20 and 13 SRE sites were found to be conserved among A. thaliana, O. sativa and S. miltiorrhiza respectively. This study provides a comprehensive picture of RNA editing events in the mitochondrial genome of S. miltiorrhiza. We identified SREs for the first time, which may represent a universal phenomenon.

Transcriptome-Wide Identification and Prediction of miRNAs and Their Targets in Paris polyphylla var. yunnanensis by High-Throughput Sequencing Analysis.

Long dormancy period of seeds limits the large-scale artificial cultivation of the scarce Paris polyphylla var. yunnanensis, an important traditional Chinese medicine. Characterizing miRNAs and their targets is crucial to understanding the role of miRNAs during seed dormancy in this species. Considering the limited genome information of this species, we first sequenced and assembled the transcriptome data of dormant seeds and their seed coats as the reference genome. A total of 146,671 unigenes with an average length of 923 bp were identified and showed functional diversity based on different annotation methods. Two small RNA libraries from respective seeds and seed coats were sequenced and the combining data indicates that 263 conserved miRNAs belonging to at least 83 families and 768 novel miRNAs in 1174 transcripts were found. The annotations of the predicted putative targets of miRNAs suggest that these miRNAs were mainly involved in the cell, metabolism and genetic information processing by direct and indirect regulation patterns in dormant seeds of P. polyphylla var. yunnanensis. Therefore, we provide the first known miRNA profiles and their targets, which will assist with further study of the molecular mechanism of seed dormancy in P. polyphylla var. yunnanensis.

IRREGULAR POLLEN EXINE1 Is a Novel Factor in Anther Cuticle and Pollen Exine Formation.

Anther cuticle and pollen exine are protective barriers for pollen development and fertilization. Despite that several regulators have been identified for anther cuticle and pollen exine development in rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana), few genes have been characterized in maize (Zea mays) and the underlying regulatory mechanism remains elusive. Here, we report a novel male-sterile mutant in maize, irregular pollen exine1 (ipe1), which exhibited a glossy outer anther surface, abnormal Ubisch bodies, and defective pollen exine. Using map-based cloning, the IPE1 gene was isolated as a putative glucose-methanol-choline oxidoreductase targeted to the endoplasmic reticulum. Transcripts of IPE1 were preferentially accumulated in the tapetum during the tetrad and early uninucleate microspore stage. A biochemical assay indicated that ipe1 anthers had altered constituents of wax and a significant reduction of cutin monomers and fatty acids. RNA sequencing data revealed that genes implicated in wax and flavonoid metabolism, fatty acid synthesis, and elongation were differentially expressed in ipe1 mutant anthers. In addition, the analysis of transfer DNA insertional lines of the orthologous gene in Arabidopsis suggested that IPE1 and their orthologs have a partially conserved function in male organ development. Our results showed that IPE1 participates in the putative oxidative pathway of C16/C18 ω-hydroxy fatty acids and controls anther cuticle and pollen exine development together with MALE STERILITY26 and MALE STERILITY45 in maize.