Centromeric DNA characterization in the model grass Brachypodium distachyon provides insights on the evolution of the genus.

Brachypodium distachyon is a well-established model monocot plant, and its small and compact genome has been used as an accurate reference for the much larger and often polyploid genomes of cereals such as Avena sativa (oats), Hordeum vulgare (barley) and Triticum aestivum (wheat). Centromeres are indispensable functional units of chromosomes and they play a core role in genome polyploidization events during evolution. As the Brachypodium genus contains about 20 species that differ significantly in terms of their basic chromosome numbers, genome size, ploidy levels and life strategies, studying their centromeres may provide important insight into the structure and evolution of the genome in this interesting and important genus. In this study, we isolated the centromeric DNA of the B. distachyon reference line Bd21 and characterized its composition via the chromatin immunoprecipitation of the nucleosomes that contain the centromere-specific histone CENH3. We revealed that the centromeres of Bd21 have the features of typical multicellular eukaryotic centromeres. Strikingly, these centromeres contain relatively few centromeric satellite DNAs; in particular, the centromere of chromosome 5 (Bd5) consists of only ~40 kb. Moreover, the centromeric retrotransposons in B. distachyon (CRBds) are evolutionarily young. These transposable elements are located both within and adjacent to the CENH3 binding domains, and have similar compositions. Moreover, based on the presence of CRBds in the centromeres, the species in this study can be grouped into two distinct lineages. This may provide new evidence regarding the phylogenetic relationships within the Brachypodium genus.

Critical function of DNA methyltransferase 1 in tomato development and regulation of the DNA methylome and transcriptome.

DNA methylation confers epigenetic regulation on gene expression and thereby on various biological processes. Tomato has emerged as an excellent system to study the function of DNA methylation in plant development. To date, regulation and function of DNA methylation maintenance remains unclear in tomato plants. Here, we report the critical function of tomato (Solanum lycopersicum) Methyltransferase 1 (SlMET1) in plant development and DNA methylome and transcriptome regulation. Using CRISPR-Cas9 gene editing, we generated slmet1 mutants and observed severe developmental defects with a frame-shift mutation, including small and curly leaves, defective inflorescence, and parthenocarpy. In leaf tissues, mutations in SlMET1 caused CG hypomethylation and CHH hypermethylation on a whole-genome scale, leading to a disturbed transcriptome including ectopic expression of many RIN target genes such as ACC2 in leaf tissues, which are normally expressed in fruits. Neither the CG hypomethylation nor CHH hypermethylation in the slmet1 mutants is related to tissue culture. Meanwhile, tissue culture induces non-CG hypomethylation, which occurs more frequently at gene regions than at TE regions. Our results depict SlMET1- and tissue culture-dependent tomato DNA methylomes, and that SlMET1 is required for maintaining a normal transcriptome and normal development of tomato.

Comparative analysis of plant MKK gene family reveals novel expansion mechanism of the members and sheds new light on functional conservation.

BACKGROUND: Mitogen-activated protein kinase (MAPK) cascades play critical functions in almost every aspect of plant growth and development, which regulates many physiological and biochemical processes. As a middle nodal point of the MAPK cascades, although evolutionary analysis of MKK from individual plant families had some reports, their evolutionary history in entire plants is still not clear. RESULTS: To better understand the evolution and function of plant MKKs, we performed systematical molecular evolutionary analysis of the MAPKK gene family and also surveyed their gene organizations, sequence features and expression patterns in different subfamilies. Phylogenetic analysis showed that plant MAPKK fall into five different groups (Group A-E). Majority orthology groups seemed to be a single or low-copy genes in all plant species analyzed in Group B, C and D, whereas group A MKKs undergo several duplication events, generating multiple gene copies. Further analysis showed that these duplication events were on account of whole genome duplications (WGDs) in plants and the duplicate genes maybe have undergone functional divergence. We also found that group E MKKs had mutation with one change of serine or theronine might lead to inactivity originated through the ancient tandem duplicates in monocots. Moreover, we also identified MKK3 integrated NTF2 domain that might have gradually lost the cytoplasmic-nuclear trafficking activity, which suggests that they may involve with the gene function more and more sophistication in the evolutionary process. Moreover, expression analyses indicated that plant MKK genes play probable roles in UV-B signaling. CONCLUSION: In general, ancient gene and genome duplications are significantly conducive to the expansion of the plant MKK gene family. Our study reveals two distinct evolutionary patterns for plant MKK proteins and sheds new light on the functional evolution of this gene family.

Genome-wide identification and evolutionary analyses of the PP2C gene family with their expression profiling in response to multiple stresses in Brachypodium distachyon.

BACKGROUND: The type-2C protein phosphatases (PP2Cs), negatively regulating ABA responses and MAPK cascade pathways, play important roles in stress signal transduction in plants. Brachypodium distachyon is a new model plant for exploring the functional genomics of temperate grasses, cereals and biofuel crops. To date, genome-wide identification and analysis of the PP2C gene family in B. distachyon have not been investigated. RESULTS: In this study, 86 PP2C genes in B. distachyon were identified. Domain-based analyses of PP2C proteins showed that they all contained the phosphatase domains featured as 11 conserved signature motifs. Although not all phosphatase domains of BdPP2C members included all 11 motifs, tertiary structure analysis showed that four residues contributing to magnesium/manganese ions (Mg(2+)/Mn(2+)) coordination were conserved, except for two noncanonical members. The analysis of their chromosomal localizations showed that most of the BdPP2C genes were located within the low CpG density region. Phylogenetic tree and synteny blocks analyses among B. distachyon, Arabidopsis thaliana and Oryza sativa revealed that all PP2C members from the three species can be phylogenetically categorized into 13 subgroups (A-M) and BdPP2Cs were evolutionarily more closely related to OsPP2Cs than to AtPP2Cs. Segmental duplications contributed particularly to the expansion of the BdPP2C gene family and all duplicated BdPP2Cs evolved mainly from purifying selection. Real-time quantitative reverse transcription PCR (qRT-PCR) analysis showed that BdPP2Cs were broadly expressed in disparate tissues. We also found that almost all members displayed up-regulation in response to abiotic stresses such as cold, heat, PEG and NaCl treatments, but down-regulation to biotic stresses such as Ph14, Guy11 and F0968 infection. CONCLUSIONS: In the present study, a comprehensive analysis of genome-wide identification and characterization of protein domains, phylogenetic relationship, gene and protein structure, chromosome location and expression pattern of the PP2C gene family was carried out for the first time in a new model monocot, i.e., B. distachyon. Our results provide a reference for genome-wide identification of the PP2C gene family of other species and also provide a foundation for future functional research on PP2C genes in B. distachyon.

Genome-wide evolutionary characterization and analysis of bZIP transcription factors and their expression profiles in response to multiple abiotic stresses in Brachypodium distachyon.

BACKGROUND: Plant basic leucine zipper (bZIP) transcription factors are one of the largest and most diverse gene families and play key roles in regulating diverse stress processes. Brachypodium distachyon is emerging as a widely recognized model plant for the temperate grass family and the herbaceous energy crops, however there is no comprehensive analysis of bZIPs in B. distachyon, especially those involved in stress tolerances. RESULTS: In this study, 96 bZIP genes (BdbZIPs) were identified distributing unevenly on each chromosome of B. distachyon, and most of them were scattered in the low CpG content regions. Gene duplications were widespread throughout B. distachyon genome. Evolutionary comparisons suggested B. distachyon and rice's bZIPs had the similar evolutionary patterns. The exon splicing in BdbZIP motifs were more complex and diverse than those in other plant species. We further revealed the potential close relationships between BdbZIP gene expressions and items including gene structure, exon splicing pattern and dimerization features. In addition, multiple stresses expression profile demonstrated that BdbZIPs exhibited significant expression patterns responding to 14 stresses, and those responding to heavy metal treatments showed opposite expression pattern comparing to the treatments of environmental factors and phytohormones. We also screened certain up- and down-regulated BdbZIP genes with fold changes ≥2, which were more sensitive to abiotic stress conditions. CONCLUSIONS: BdbZIP genes behaved diverse functional characters and showed discrepant and some regular expression patterns in response to abiotic stresses. Comprehensive analysis indicated these BdbZIPs' expressions were associated not only with gene structure, exon splicing pattern and dimerization feature, but also with abiotic stress treatments. It is possible that our findings are crucial for revealing the potentialities of utilizing these candidate BdbZIPs to improve productivity of grass plants and cereal crops.

Genome-wide exploration of the molecular evolution and regulatory network of mitogen-activated protein kinase cascades upon multiple stresses in Brachypodium distachyon.

BACKGROUND: Brachypodium distachyon is emerging as a widely recognized model plant that has very close relations with several economically important Poaceae species. MAPK cascade is known to be an evolutionarily conserved signaling module involved in multiple stresses. Although the gene sequences of MAPK and MAPKK family have been fully identified in B. distachyon, the information related to the upstream MAPKKK gene family especially the regulatory network among MAPKs, MAPKKs and MAPKKKs upon multiple stresses remains to be understood. RESULTS: In this study, we have identified MAPKKKs which belong to the biggest gene family of MAPK cascade kinases. We have systematically investigated the evolution of whole MAPK cascade kinase gene family in terms of gene structures, protein structural organization, chromosomal localization, orthologs construction and gene duplication analysis. Our results showed that most BdMAPK cascade kinases were located at the low-CpG-density region, and the clustered members in each group shared similar structures of the genes and proteins. Synteny analysis showed that 62 or 21 pairs of duplicated orthologs were present between B. distachyon and Oryza sativa, or between B. distachyon and Arabidopsis thaliana respectively. Gene expression data revealed that BdMAPK cascade kinases were rapidly regulated by stresses and phytohormones. Importantly, we have constructed a regulation network based on co-expression patterns of the expression profiles upon multiple stresses performed in this study. CONCLUSIONS: BdMAPK cascade kinases were involved in the signaling pathways of multiple stresses in B. distachyon. The network of co-expression regulation showed the most of duplicated BdMAPK cascade kinase gene orthologs demonstrated their convergent function, whereas few of them developed divergent function in the evolutionary process. The molecular evolution analysis of identified MAPK family genes and the constructed MAPK cascade regulation network under multiple stresses provide valuable information for further investigation of the functions of BdMAPK cascade kinase genes.

Transcriptional profiling analysis of individual kinase-deletion strains of fission yeast in response to nitrogen starvation.

Nitrogen starvation (NS) induces sexual development when mating partners are available or enter into quiescent state (G0) in heterothallic background in fission yeast. However, little is known whether the two processes share common signaling molecules or cells defective in the two processes share common transcriptional signatures. To address these questions, we first assessed 77 kinase-deletion strains for NS-induced G0-arrest phenotypes. Our result indicated that 10 out of 77 kinase-deletion strains exhibited defect in G0-arrest, only 3 of which were defective in sexual development based on a previous study, suggesting that the two processes hardly share common signaling components. We subsequently performed transcriptional profiling analysis. Our result indicated that NS-induced transcriptional change was so robust that it prevailed the alteration by individual kinase-deletion alleles. Based on comparison between kinase-deletion strains proficient and deficient in sexual development or G0-arrest, we identified subsets of genes that were associated with sexual development-deficient or G0-arrest-deficient kinase-deletion strains. Multiple pairing analyses allowed grouping of functional related kinases. Furthermore, we showed that Pka1-mediated pathways were required for upregulation of NS-induced genes upon NS and downregulation of the same set of genes under the N-replete conditions. Taken together, our analyses indicate that sexual development and NS-induced G0-arrest are unrelated; and sexual development-deficient and G0-arrest-deficient kinase-deletion strains possess distinct transcriptional signatures. We propose that Pka1 is a key regulator of nitrogen metabolic pathways and Pka1-mediated signaling pathways play roles in regulation of NS-induced genes under both N-depleted and N-replete conditions.

Deconvolution of chromatin immunoprecipitation-microarray (ChIP-chip) analysis of MBF occupancies reveals the temporal recruitment of Rep2 at the MBF target genes.

MBF (or DSC1) is known to regulate transcription of a set of G(1)/S-phase genes encoding proteins involved in regulation of DNA replication. Previous studies have shown that MBF binds not only the promoter of G(1)/S-phase genes, but also the constitutive genes; however, it was unclear if the MBF bindings at the G(1)/S-phase and constitutive genes were mechanistically distinguishable. Here, we report a chromatin immunoprecipitation-microarray (ChIP-chip) analysis of MBF binding in the Schizosaccharomyces pombe genome using high-resolution genome tiling microarrays. ChIP-chip analysis indicates that the majority of the MBF occupancies are located at the intragenic regions. Deconvolution analysis using Rpb1 ChIP-chip results distinguishes the Cdc10 bindings at the Rpb1-poor loci (promoters) from those at the Rpb1-rich loci (intragenic sequences). Importantly, Res1 binding at the Rpb1-poor loci, but not at the Rpb1-rich loci, is dependent on the Cdc10 function, suggesting a distinct binding mechanism. Most Cdc10 promoter bindings at the Rpb1-poor loci are associated with the G(1)/S-phase genes. While Res1 or Res2 is found at both the Cdc10 promoter and intragenic binding sites, Rep2 appears to be absent at the Cdc10 promoter binding sites but present at the intragenic sites. Time course ChIP-chip analysis demonstrates that Rep2 is temporally accumulated at the coding region of the MBF target genes, resembling the RNAP-II occupancies. Taken together, our results show that deconvolution analysis of Cdc10 occupancies refines the functional subset of genomic binding sites. We propose that the MBF activator Rep2 plays a role in mediating the cell cycle-specific transcription through the recruitment of RNAP-II to the MBF-bound G(1)/S-phase genes.

R gene expression induced by a type-III effector triggers disease resistance in rice.

Disease resistance (R) genes in plants encode products that specifically recognise incompatible pathogens and trigger a cascade of events leading to disease resistance in the host plant. R-gene specificity is dictated by both host R genes and cognate avirulence (avr) genes in pathogens. However, the basis of gene-for-gene specificity is not well understood. Here, we report the cloning of the R gene Xa27 from rice and the cognate avr gene avrXa27 from Xanthomonas oryzae pv. oryzae. Resistant and susceptible alleles of Xa27 encode identical proteins. However, expression of only the resistant allele occurs when a rice plant is challenged by bacteria harbouring avrXa27, whose product is a nuclear localized type-III effector. Induction of Xa27 occurs only in the immediate vicinity of infected tissue, whereas ectopic expression of Xa27 resulted in resistance to otherwise compatible strains of the pathogen. Thus Xa27 specificity towards incompatible pathogens involves the differential expression of the R gene in the presence of the AvrXa27 effector.

Modulation of cell cycle-specific gene expressions at the onset of S phase arrest contributes to the robust DNA replication checkpoint response in fission yeast.

Fission yeast replication checkpoint kinases Rad3p and Cds1p are essential for maintaining cell viability after transient treatment with hydroxyurea (HU), an agent that blocks DNA replication. Although current studies have focused on the cyclin-dependent protein kinase Cdc2p that is regulated by these checkpoint kinases, other aspects of their functions at the onset of S phase arrest have not been fully understood. In this study, we use genome-wide DNA microarray analyses to show that HU-induced change of expression profiles in synchronized G(2) cells occurs specifically at the onset of S phase arrest. Induction of many core environmental stress response genes and repression of ribosomal genes happen during S phase arrest. Significantly, peak expression level of the MluI-like cell cycle box (MCB)-cluster (G(1)) genes is maintained at the onset of S phase arrest in a Rad3p- and Cds1p-dependent manner. Expression level maintenance of the MCB-cluster is mediated through the accumulation of Rep2p, a putative transcriptional activator of the MBF complex. Conversely, the FKH-cluster (M) genes are repressed during the onset of S phase arrest in a Rad3p-dependent manner. Repression of the FKH-cluster genes is mediated through the decreased levels of one of the putative forkhead transcription factors, Sep1p, but not Fkh2p. Together, our results demonstrate that Rad3p and Cds1p modulate transcriptional response during the onset of S phase arrest.

Integrated bioinformatics analysis of As, Au, Cd, Pb and Cu heavy metal responsive marker genes through Arabidopsis thaliana GEO datasets.

BACKGROUND: Current environmental pollution factors, particularly the distribution and diffusion of heavy metals in soil and water, are a high risk to local environments and humans. Despite striking advances in methods to detect contaminants by a variety of chemical and physical solutions, these methods have inherent limitations such as small dimensions and very low coverage. Therefore, identifying novel contaminant biomarkers are urgently needed. METHODS: To better track heavy metal contaminations in soil and water, integrated bioinformatics analysis to identify biomarkers of relevant heavy metal, such as As, Cd, Pb and Cu, is a suitable method for long-term and large-scale surveys of such heavy metal pollutants. Subsequently, the accuracy and stability of the results screened were experimentally validated by quantitative PCR experiment. RESULTS: We obtained 168 differentially expressed genes (DEGs) which contained 59 up-regulated genes and 109 down-regulated genes through comparative bioinformatics analyses. Subsequently, the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichments of these DEGs were performed, respectively. GO analyses found that these DEGs were mainly related to responses to chemicals, responses to stimulus, responses to stress, responses to abiotic stimulus, and so on. KEGG pathway analyses of DEGs were mainly involved in the protein degradation process and other biologic process, such as the phenylpropanoid biosynthesis pathways and nitrogen metabolism. Moreover, we also speculated that nine candidate core biomarker genes (namely, NILR1, PGPS1, WRKY33, BCS1, AR781, CYP81D8, NR1, EAP1 and MYB15) might be tightly correlated with the response or transport of heavy metals. Finally, experimental results displayed that these genes had the same expression trend response to different stresses as mentioned above (Cd, Pb and Cu) and no mentioned above (Zn and Cr). CONCLUSION: In general, the identified biomarker genes could help us understand the potential molecular mechanisms or signaling pathways responsive to heavy metal stress in plants, and could be applied as marker genes to track heavy metal pollution in soil and water through detecting their expression in plants growing in those environments.

Partial defoliation of Brachypodium distachyon plants grown in petri dishes under low light increases P and other nutrient levels concomitantly with transcriptional changes in the roots.

BACKGROUND: There have been few studies on the partial defoliation response of grass. It has been unclear how partial defoliation may affect roots at the levels of nutrient accumulation and transcriptional regulation. Hereby we report a comprehensive investigation on molecular impacts of partial defoliation by using a model grass species, Brachypodium distachyon. RESULTS: Our Inductively Coupled Plasma Mass Spectrometry analyses of B. distachyon revealed shoot- and root-specific accumulation patterns of a group of macronutrients including potassium (K), Phosphorus (P), Calcium (Ca), Magnesium (Mg), and micronutrients including Sodium (Na), iron (Fe), and Manganese (Mn). Meanwhile, our genome-wide profiling of gene expression patterns depicts transcriptional impacts on B. distachyon roots by cutting the aerial portion. The RNAseq analyses identified a total of 1,268 differentially expressed genes in B. distachyon with partial defoliation treatment. Our comprehensive analyses by means of multiple approaches, including Gene Ontology, InterPro and Pfam protein classification, KEGG pathways, and Plant TFDB, jointly highlight the involvement of hormone-mediated wounding response, primary and secondary metabolites, and ion homeostasis, in B. distachyon after the partial defoliation treatment. In addition, evidence is provided that roots respond to partial defoliation by modifying nutrient uptake and rhizosphere acidification rate, indicating that an alteration of the root/soil interaction occurs in response to this practice. CONCLUSIONS: This study reveals how partial defoliation alters ion accumulation levels in shoots and roots, as well as partial defoliation-induced transcriptional reprogramming on a whole-genome scale, thereby providing insight into the molecular mechanisms underlying the recovery process of grass after partial defoliation.

In silico-prediction of protein-protein interactions network about MAPKs and PP2Cs reveals a novel docking site variants in Brachypodium distachyon.

Protein-protein interactions (PPIs) underlie the molecular mechanisms of most biological processes. Mitogen-activated protein kinases (MAPKs) can be dephosphorylated by MAPK-specific phosphatases such as PP2C, which are critical to transduce extracellular signals into adaptive and programmed responses. However, the experimental approaches for identifying PPIs are expensive, time-consuming, laborious and challenging. In response, many computational methods have been developed to predict PPIs. Yet, these methods have inherent disadvantages such as high false positive and negative results. Thus, it is crucial to develop in silico approaches for predicting PPIs efficiently and accurately. In this study, we identified PPIs among 16 BdMAPKs and 86 BdPP2Cs in B. distachyon using a novel docking approach. Further, we systematically investigated the docking site (D-site) of BdPP2C which plays a vital role for recognition and docking of BdMAPKs. D-site analysis revealed that there were 96 pairs of PPIs including all BdMAPKs and most BdPP2Cs, which indicated that BdPP2C may play roles in other signaling networks. Moreover, most BdPP2Cs have a D-site for BdMAPKs in our prediction results, which suggested that our method can effectively predict PPIs, as confirmed by their 3D structure. In addition, we validated this methodology with known Arabidopsis and yeast phosphatase-MAPK interactions from the STRING database. The results obtained provide a vital research resource for exploring an accurate network of PPIs between BdMAPKs and BdPP2Cs.

Genome Sequence of Bacillus megaterium Strain YC4-R4, a Plant Growth-Promoting Rhizobacterium Isolated from a High-Salinity Environment.

Here, we report the complete genome sequence for Bacillus megaterium strain YC4-R4, a highly salt-tolerant rhizobacterium that promotes growth in plants. The sequencing process was performed by combining pyrosequencing and single-molecule sequencing techniques. The complete genome is estimated to be approximately 5.44 Mb, containing a total of 5,673 predicted protein-coding DNA sequences (CDSs).

Complete Genome Sequence of Bacillus megaterium Strain TG1-E1, a Plant Drought Tolerance-Enhancing Bacterium.

Based on a combination of next-generation sequencing and single-molecule sequencing, we obtained the whole-genome sequence of Bacillus megaterium strain TG1-E1, which is a highly salt-tolerant rhizobacterium that enhances plant tolerance to drought stress. The complete genome is estimated to be approximately 5.48 Mb containing a total of 5,858 predicted protein-coding DNA sequences.

Adaptive expression responses in the Pol-gamma null strain of S. pombe depleted of mitochondrial genome.

BACKGROUND: DNA polymerase gamma(Pol-gamma) has been shown to be essential for maintenance of the mitochondrial genome (mtDNA) in the petite-positive budding yeast Saccharomyces cerevisiae. Budding yeast cells lacking mitochondria exhibit a slow-growing or petite-colony phenotype. Petite strains fail to grow on non-fermentable carbon sources. However, it is not clear whether the Pol-gamma is required for mtDNA maintenance in the petite-negative fission yeast Schizosaccharomyces pombe. RESULTS: We show that disruption of the nuclear gene pog1+ that encodes Pol-gamma is sufficient to deplete mtDNA in S. pombe. Cells bearing pog1Delta allele require substantial growth periods to form petite colonies. Mitotracker assays indicate that pog1Delta cells are defective in mitochondrial function and EM analyses suggest that pog1Delta cells lack normal mitochondrial structures. Depletion of mtDNA in pog1Delta cells is evident from quantitative real-time PCR assays. Genome-wide expression profiles of pog1Delta and other mtDNA-less cells reveal that many genes involved in response to stimulus, energy derivation by oxidation of organic compounds, cellular carbohydrate metabolism, and energy reserve metabolism are induced. Conversely, many genes encoding proteins involved in amino acid metabolism and oxidative phosphorylation are repressed. CONCLUSION: By showing that Pol-gamma is essential for mtDNA maintenance and disruption of pog1+ alters the genome-wide expression profiles, we demonstrated that cells lacking mtDNA exhibit adaptive nuclear gene expression responses in the petite-negative S. pombe.

Roles of Nuclear Pores and Nucleo-cytoplasmic Trafficking in Plant Stress Responses.

The nuclear pore complex (NPC) is a large protein complex that controls the exchange of components between the nucleus and the cytoplasm. In plants, the NPC family components play critical roles not only in essential growth and developmental processes, but also in plant responses to various environmental stress conditions. The involvement of NPC components in plant stress responses is mainly attributed to different mechanisms including control of mRNA/protein nucleo-cytoplasmic trafficking and transcriptional gene regulation. This mini review summarizes current knowledge of the NPC-mediated plant stress responses and provides an overview of the underlying molecular mechanisms.

Methylated-antibody affinity purification to improve proteomic identification of plant RNA polymerase Pol V complex and the interacting proteins.

Affinity purification followed by enzymatic digestion and mass spectrometry has been widely utilized for the sensitive detection of interacting proteins and protein complexes in various organisms. In plants, the method is technically challenging due to the low abundance proteins, non-specific binding and difficulties of eluting interacting proteins from antibody beads. In this report, we describe a strategy to modify antibodies by reductive methylation of lysines without affecting their binding properties, followed by on-bead digestion of bound proteins with endoproteinase Lys-C. By this method, the antibody remains intact and does not interfere with the downstream identification of interacting proteins. Non-specific binding proteins were excluded using 14N/15N-metabolic labeling of wild-type and the transgenic plant counterparts. The method was employed to identify 12 co-immunoprecipitated protein subunits in Pol V complex and to discover 17 potential interacting protein targets in Arabidopsis. Our results demonstrated that the modification of antibodies by reductive dimethylation can improve the reliability and sensitivity of identifying low-abundance proteins through on-bead digestion and mass spectrometry. We also show that coupling this technique with chemical crosslinking enables in-depth characterization of endogenous protein complexes and the protein-protein interaction networks including mapping the surface topology and post-translational modifications of interacting proteins.

Genome-wide evolutionary characterization and expression analyses of WRKY family genes in Brachypodium distachyon.

Members of plant WRKY gene family are ancient transcription factors that function in plant growth and development and respond to biotic and abiotic stresses. In our present study, we have investigated WRKY family genes in Brachypodium distachyon, a new model plant of family Poaceae. We identified a total of 86 WRKY genes from B. distachyon and explored their chromosomal distribution and evolution, domain alignment, promoter cis-elements, and expression profiles. Combining the analysis of phylogenetic tree of BdWRKY genes and the result of expression profiling, results showed that most of clustered gene pairs had higher similarities in the WRKY domain, suggesting that they might be functionally redundant. Neighbour-joining analysis of 301 WRKY domains from Oryza sativa, Arabidopsis thaliana, and B. distachyon suggested that BdWRKY domains are evolutionarily more closely related to O. sativa WRKY domains than those of A. thaliana. Moreover, tissue-specific expression profile of BdWRKY genes and their responses to phytohormones and several biotic or abiotic stresses were analysed by quantitative real-time PCR. The results showed that the expression of BdWRKY genes was rapidly regulated by stresses and phytohormones, and there was a strong correlation between promoter cis-elements and the phytohormones-induced BdWRKY gene expression.

Validation of a reference gene (BdFIM) for quantifying transgene copy numbers in Brachypodium distachyon by real-time PCR.

Brachypodium distachyon has been proposed as a new model system for gramineous plants with a sequenced genome and an efficient transformation system. Many transgenic B. distachyon plants have been generated in recent years. To develop a reliable fast method for detecting transgenic B. distachyon and quantifying its transgene copy numbers, a species-specific reference gene is of great priority to be validated both in qualitative PCR and quantitative real-time PCR detection. In this study, we first proved that the BdFIM (B. distachyon fimbrin-like protein) gene is a suitable reference gene in qualitative PCR and quantitative real-time PCR for B. distachyon. Fourteen different B. distachyon varieties were tested by both qualitative and quantitative PCRs, and identical amplification products of BdFIM were obtained with all of them, while no amplification products were observed with samples from 14 other plant species, suggesting that BdFIM gene was specific to B. distachyon. The results of Southern blot analysis revealed that the BdFIM gene was low copy number in seven tested B. distachyon varieties. In conclusion, the BdFIM gene can be used as a reference gene, since it had species specificity, low heterogeneity, and low copy number among the tested B. distachyon varieties. Furthermore, the copy number of inserted sequences from transgenic B. distachyon obtained by real-time PCR methods and Southern blot confirmed that the BdFIM gene was an applicable reference gene in B. distachyon.