Extensive Phenotypic Variation among Allelic T-DNA Inserts in Arabidopsis thaliana

T-DNA insertion mutants are a tool used widely in Arabidopsis thaliana to disrupt gene function. We phenotyped multiple homozygous T-DNA A. thaliana mutants at each of two loci (AT1G11060 and AT4G00210). We measured life history traits, including germination, size at reproduction and fruit production. Allelic T-DNA lines differed for most traits at AT1G11060 but not at AT4G00210. However, insertions in exons differed from other insertion positions in AT4G00210 but not in AT1G11060. We found evidence for additional insertions in approximately half of the lines, but found few phenotypic consequences. In general, our results suggest that a cautious interpretation of T-DNA phenotypes is warranted.     

Computational Prediction of Candidate Proteins for S-Nitrosylation in Arabidopsis thaliana

Nitric oxide (NO) is an important signaling molecule that regulates many physiological processes in plants. One of the most important regulatory mechanisms of NO is S-nitrosylation—the covalent attachment of NO to cysteine residues. Although the involvement of cysteine S-nitrosylation in the regulation of protein functions is well established, its substrate specificity remains unknown. Identification of candidates for S-nitrosylation and their target cysteine residues is fundamental for studying the molecular mechanisms and regulatory roles of S-nitrosylation in plants. Several experimental methods that are based on the biotin switch have been developed to identify target proteins for S-nitrosylation. However, these methods have their limits. Thus, computational methods are attracting considerable attention for the identification of modification sites in proteins. Using GPS-SNO version 1.0, a recently developed S-nitrosylation site-prediction program, a set of 16,610 candidate proteins for S-nitrosylation containing 31,900 S-nitrosylation sites was isolated from the entire Arabidopsis proteome using the medium threshold. In the compartments “chloroplast,” “CUL4-RING ubiquitin ligase complex,” and “membrane” more than 70% of the proteins were identified as candidates for S-nitrosylation. The high number of identified candidates in the proteome reflects the importance of redox signaling in these compartments. An analysis of the functional distribution of the predicted candidates showed that proteins involved in signaling processes exhibited the highest prediction rate. In a set of 46 proteins, where 53 putative S-nitrosylation sites were already experimentally determined, the GPS-SNO program predicted 60 S-nitrosylation sites, but only 11 overlap with the results of the experimental approach. In general, a computer-assisted method for the prediction of targets for S-nitrosylation is a very good tool; however, further development, such as including the three dimensional structure of proteins in such analyses, would improve the identification of S-nitrosylation sites.     

Herkogamy and Its Effects on Mating Patterns in Arabidopsis thaliana

The evolution of mating systems, which exhibit an extraordinary diversity in flowering plants, is of central interest in plant biology. Herkogamy, the spatial separation of sexual organs within flowers, is a widespread floral mechanism that is thought to be an adaptive trait reducing self-pollination in hermaphroditic plants. In contrast with previous studies of herkogamy that focused on plants with relatively large floral displays, we here characterized herkogamy in Arabidopsis thaliana, a model plant with a strong selfing syndrome. Developmental features, reproductive consequences, and genetic architecture of herkogamy were exploited using naturally variable A. thaliana accessions, under both greenhouse and natural conditions. Our results demonstrate that the degree of herkogamy can strongly influence the mating patterns of A. thaliana: approach herkogamy can effectively promote outcrossing, no herkogamy is also capable of enhancing the opportunity for outcrossing, and reverse herkogamy facilitates efficient self-pollination. In addition, we found that the expression of herkogamy in A. thaliana was environment-dependent and regulated by multiple quantitative trait loci. This study reveals how minor modifications in floral morphology may cause dramatic changes in plant mating patterns, provides new insights into the function of herkogamy, and suggests the way for dissecting the genetic basis of this important character in a model plant.     

Natural Variation Identifies Multiple Loci Controlling Petal Shape and Size in Arabidopsis thaliana

Natural variation in organ morphologies can have adaptive significance and contribute to speciation. However, the underlying allelic differences responsible for variation in organ size and shape remain poorly understood. We have utilized natural phenotypic variation in three Arabidopsis thaliana ecotypes to examine the genetic basis for quantitative variation in petal length, width, area, and shape. We identified 23 loci responsible for such variation, many of which appear to correspond to genes not previously implicated in controlling organ morphology. These analyses also demonstrated that allelic differences at distinct loci can independently affect petal length, width, area or shape, suggesting that these traits behave as independent modules. We also showed that ERECTA (ER), encoding a leucine-rich repeat (LRR) receptor-like serine-threonine kinase, is a major effect locus determining petal shape. Allelic variation at the ER locus was associated with differences in petal cell proliferation and concomitant effects on petal shape. ER has been previously shown to be required for regulating cell division and expansion in other contexts; the ER receptor-like kinase functioning to also control organ-specific proliferation patterns suggests that allelic variation in common signaling components may nonetheless have been a key factor in morphological diversification.     

Structural and Functional Characterization of the Protein Kinase Mps1 in Arabidopsis thaliana

In eukaryotes, protein kinases catalyze the transfer of a gamma-phosphate from ATP (or GTP) to specific amino acids in protein targets. In plants, protein kinases have been shown to participate in signaling cascades driving responses to environmental stimuli and developmental processes. Plant meristems are undifferentiated tissues that provide the major source of cells that will form organs throughout development. However, non-dividing specialized cells can also dedifferentiate and re-initiate cell division if exposed to appropriate conditions. Mps1 (Monopolar spindle) is a dual-specificity protein kinase that plays a critical role in monitoring the accuracy of chromosome segregation in the mitotic checkpoint mechanism. Although Mps1 functions have been clearly demonstrated in animals and fungi, its role in plants is so far unclear. Here, using structural and biochemical analyses here we show that Mps1 has highly similar homologs in many plant genomes across distinct lineages (e.g. AtMps1 in Arabidopsis thaliana). Several structural features (i.e. catalytic site, DFG motif and threonine triad) are clearly conserved in plant Mps1 kinases. Structural and sequence analysis also suggest that AtMps1 interact with other cell cycle proteins, such as Mad2 and MAPK1. By using a very specific Mps1 inhibitor (SP600125) we show that compromised AtMps1 activity hampers the development of A. thaliana seedlings in a dose-dependent manner, especially in secondary roots. Moreover, concomitant administration of the auxin IAA neutralizes the AtMps1 inhibition phenotype, allowing secondary root development. These observations let us to hypothesize that AtMps1 might be a downstream regulator of IAA signaling in the formation of secondary roots. Our results indicate that Mps1 might be a universal component of the Spindle Assembly Checkpoint machinery across very distant lineages of eukaryotes.     

"Missing" G x E Variation Controls Flowering Time in Arabidopsis thaliana

Understanding how genetic variation interacts with the environment is essential for understanding adaptation. In particular, the life cycle of plants is tightly coordinated with local environmental signals through complex interactions with the genetic variation (G x E). The mechanistic basis for G x E is almost completely unknown. We collected flowering time data for 173 natural inbred lines of Arabidopsis thaliana from Sweden under two growth temperatures (10°C and 16°C), and observed massive G x E variation. To identify the genetic polymorphisms underlying this variation, we conducted genome-wide scans using both SNPs and local variance components. The SNP-based scan identified several variants that had common effects in both environments, but found no trace of G x E effects, whereas the scan using local variance components found both. Furthermore, the G x E effects appears to be concentrated in a small fraction of the genome (0.5%). Our conclusion is that G x E effects in this study are mostly due to large numbers of allele or haplotypes at a small number of loci, many of which correspond to previously identified flowering time genes.     

Diuretics Prime Plant Immunity in Arabidopsis thaliana

Plant activators are agrochemicals that activate the plant immune system, thereby enhancing disease resistance. Due to their prophylactic and durable effects on a wide spectrum of diseases, plant activators can provide synergistic crop protection when used in combination with traditional pest controls. Although plant activators have achieved great success in wet-rice farming practices in Asia, their use is still limited. To isolate novel plant activators applicable to other crops, we screened a chemical library using a method that can selectively identify immune-priming compounds. Here, we report the isolation and characterization of three diuretics, bumetanide, bendroflumethiazide and clopamide, as immune-priming compounds. These drugs upregulate the immunity-related cell death of Arabidopsis suspension-cultured cells induced with an avirulent strain of Pseudomonas syringae pv. tomato in a concentration-dependent manner. The application of these compounds to Arabidopsis plants confers disease resistance to not only the avirulent but also a virulent strain of the pathogen. Unlike salicylic acid, an endogenous phytohormone that governs disease resistance in response to biotrophic pathogens, the three diuretic compounds analyzed here do not induce PR1 or inhibit plant growth, showing potential as lead compounds in a practical application.     

Regulation of Active DNA Demethylation by a Methyl-CpG-Binding Domain Protein in Arabidopsis thaliana

Active DNA demethylation plays crucial roles in the regulation of gene expression in both plants and animals. In Arabidopsis thaliana, active DNA demethylation is initiated by the ROS1 subfamily of 5-methylcytosine-specific DNA glycosylases via a base excision repair mechanism. Recently, IDM1 and IDM2 were shown to be required for the recruitment of ROS1 to some of its target loci. However, the mechanism(s) by which IDM1 is targeted to specific genomic loci remains to be determined. Affinity purification of IDM1- and IDM2- associating proteins demonstrated that IDM1 and IDM2 copurify together with two novel components, methyl-CpG-binding domain protein 7 (MBD7) and IDM2-like protein 1 (IDL1). IDL1 encodes an α-crystallin domain protein that shows high sequence similarity with IDM2. MBD7 interacts with IDM2 and IDL1 in vitro and in vivo and they form a protein complex associating with IDM1 in vivo. MBD7 directly binds to the target loci and is required for the H3K18 and H3K23 acetylation in planta. MBD7 dysfunction causes DNA hypermethylation and silencing of reporter genes and a subset of endogenous genes. Our results suggest that a histone acetyltransferase complex functions in active DNA demethylation and in suppression of gene silencing at some loci in Arabidopsis.     

The Structure of the Karrikin-Insensitive Protein (KAI2) in Arabidopsis thaliana

KARRIKIN INSENSITIVE 2 (KAI2) is an α/β hydrolase involved in seed germination and seedling development. It is essential for plant responses to karrikins, a class of butenolide compounds derived from burnt plant material that are structurally similar to strigolactone plant hormones. The mechanistic basis for the function of KAI2 in plant development remains unclear. We have determined the crystal structure of Arabidopsis thaliana KAI2 in space groups P2₁ 2₁ 2₁ (a  = 63.57 Å, b  = 66.26 Å, c  = 78.25 Å) and P2₁ (a  = 50.20 Å, b  = 56.04 Å, c  = 52.43 Å, β  = 116.12°) to 1.55 and 2.11 Å respectively. The catalytic residues are positioned within a large hydrophobic pocket similar to that of DAD2, a protein required for strigolactone response in Petunia hybrida. KAI2 possesses a second solvent-accessible pocket, adjacent to the active site cavity, which offers the possibility of allosteric regulation. The structure of KAI2 is consistent with its designation as a serine hydrolase, as well as previous data implicating the protein in karrikin and strigolactone signalling.     

Identification of Nitrogen Starvation-Responsive MicroRNAs in Arabidopsis thaliana

microRNAs (miRNAs) are a class of negative regulators that take part in many processes such as growth and development, stress responses, and metabolism in plants. Recently, miRNAs were shown to function in plant nutrient metabolism. Moreover, several miRNAs were identified in the response to nitrogen (N) deficiency. To investigate the functions of other miRNAs in N deficiency, deep sequencing technology was used to detect the expression of small RNAs under N-sufficient and -deficient conditions. The results showed that members from the same miRNA families displayed differential expression in response to N deficiency. Upon N starvation, the expression of miR169, miR171, miR395, miR397, miR398, miR399, miR408, miR827, and miR857 was repressed, whereas those of miR160, miR780, miR826, miR842, and miR846 were induced. miR826, a newly identified N-starvation-induced miRNA, was found to target the AOP2 gene. Among these N-starvation-responsive miRNAs, several were involved in cross-talk among responses to different nutrient (N, P, S, Cu) deficiencies. miR160, miR167, and miR171 could be responsible for the development of Arabidopsis root systems under N-starvation conditions. In addition, twenty novel miRNAs were identified and nine of them were significantly responsive to N-starvation. This study represents comprehensive expression profiling of N-starvation-responsive miRNAs and advances our understanding of the regulation of N homeostasis mediated by miRNAs.     

CpLEPA Is Critical for Chloroplast Protein Synthesis Under Suboptimal Conditions in Arabidopsis thaliana

LEPA is one of the most conserved translation factors and is found from bacteria to higher plants. However, the physiological function of the chloroplast LEPA homolog in higher plants remains unknown. Herein, we demonstrate the physiological role of cpLEPA in enabling efficient photosynthesis in higher plants. The cplepa-1 mutant displays slightly high chlorophyll fluorescence and pale green phenotypes under normal growth conditions. The growth of the cplepa-1 mutant is reduced when grown on soil, and greater reduction is observed under intense light illumination. Photosynthetic activity is impaired in the cplepa-1 mutants, which is reflected in the decreased steady-state levels of chloroplast proteins. In vivo protein labeling experiments explained the decrease in the steady-state levels of chloroplast proteins. An abnormal association of the chloroplast-encoded mRNAs with ribosomes suggests that the protein synthesis deficiencies in cplepa-1 are due to defects in translation initiation in the chloroplasts. The cpLEPA protein appears to be an essential translation factor that promotes the efficiency of chloroplast protein synthesis.     

A Genome-Wide Association Analysis Reveals Epistatic Cancellation of Additive Genetic Variance for Root Length in Arabidopsis thaliana

Efforts to identify loci underlying complex traits generally assume that most genetic variance is additive. Here, we examined the genetics of Arabidopsis thaliana root length and found that the genomic narrow-sense heritability for this trait in the examined population was statistically zero. The low amount of additive genetic variance that could be captured by the genome-wide genotypes likely explains why no associations to root length could be found using standard additive-model-based genome-wide association (GWA) approaches. However, as the broad-sense heritability for root length was significantly larger, and primarily due to epistasis, we also performed an epistatic GWA analysis to map loci contributing to the epistatic genetic variance. Four interacting pairs of loci were revealed, involving seven chromosomal loci that passed a standard multiple-testing corrected significance threshold. The genotype-phenotype maps for these pairs revealed epistasis that cancelled out the additive genetic variance, explaining why these loci were not detected in the additive GWA analysis. Small population sizes, such as in our experiment, increase the risk of identifying false epistatic interactions due to testing for associations with very large numbers of multi-marker genotypes in few phenotyped individuals. Therefore, we estimated the false-positive risk using a new statistical approach that suggested half of the associated pairs to be true positive associations. Our experimental evaluation of candidate genes within the seven associated loci suggests that this estimate is conservative; we identified functional candidate genes that affected root development in four loci that were part of three of the pairs. The statistical epistatic analyses were thus indispensable for confirming known, and identifying new, candidate genes for root length in this population of wild-collected A. thaliana accessions. We also illustrate how epistatic cancellation of the additive genetic variance explains the insignificant narrow-sense and significant broad-sense heritability by using a combination of careful statistical epistatic analyses and functional genetic experiments.     

Structural Insights of the ssDNA Binding Site in the Multifunctional Endonuclease AtBFN2 from Arabidopsis thaliana

The multi S1/P1 nuclease AtBFN2 (EC 3.1.30.1) encoded by the Arabidopsis thaliana At1g68290 gene is a glycoprotein that digests RNA, ssDNA, and dsDNA. AtBFN2 depends on three zinc ions for cleaving DNA and RNA at 3′-OH to yield 5′-nucleotides. In addition, AtBFN2′s enzymatic activity is strongly glycan dependent. Plant Zn²⁺-dependent endonucleases present a unique fold, and belong to the Phospholipase C (PLC)/P1 nuclease superfamily. In this work, we present the first complete, ligand-free, AtBFN2 crystal structure, along with sulfate, phosphate and ssDNA co-crystal structures. With these, we were able to provide better insight into the glycan structure and possible enzymatic mechanism. In comparison with other nucleases, the AtBFN2/ligand-free and AtBFN2/PO₄ models suggest a similar, previously proposed, catalytic mechanism. Our data also confirm that the phosphate and vanadate can inhibit the enzyme activity by occupying the active site. More importantly, the AtBFN2/A₅T structure reveals a novel and conserved secondary binding site, which seems to be important for plant Zn²⁺-dependent endonucleases. Based on these findings, we propose a rational ssDNA binding model, in which the ssDNA wraps itself around the protein and the attached surface glycan, in turn, reinforces the binding complex.     

Accuracy of Protein-Protein Binding Sites in High-Throughput Template-Based Modeling

The accuracy of protein structures, particularly their binding sites, is essential for the success of modeling protein complexes. Computationally inexpensive methodology is required for genome-wide modeling of such structures. For systematic evaluation of potential accuracy in high-throughput modeling of binding sites, a statistical analysis of target-template sequence alignments was performed for a representative set of protein complexes. For most of the complexes, alignments containing all residues of the interface were found. The full interface alignments were obtained even in the case of poor alignments where a relatively small part of the target sequence (as low as 40%) aligned to the template sequence, with a low overall alignment identity (<30%). Although such poor overall alignments might be considered inadequate for modeling of whole proteins, the alignment of the interfaces was strong enough for docking. In the set of homology models built on these alignments, one third of those ranked 1 by a simple sequence identity criteria had RMSD<5 Å, the accuracy suitable for low-resolution template free docking. Such models corresponded to multi-domain target proteins, whereas for single-domain proteins the best models had 5 Å<RMSD<10 Å, the accuracy suitable for less sensitive structure-alignment methods. Overall, ∼50% of complexes with the interfaces modeled by high-throughput techniques had accuracy suitable for meaningful docking experiments. This percentage will grow with the increasing availability of co-crystallized protein-protein complexes.     

Accurate Prediction of Peptide Binding Sites on Protein Surfaces

Many important protein–protein interactions are mediated by the binding of a short peptide stretch in one protein to a large globular segment in another. Recent efforts have provided hundreds of examples of new peptides binding to proteins for which a three-dimensional structure is available (either known experimentally or readily modeled) but where no structure of the protein–peptide complex is known. To address this gap, we present an approach that can accurately predict peptide binding sites on protein surfaces. For peptides known to bind a particular protein, the method predicts binding sites with great accuracy, and the specificity of the approach means that it can also be used to predict whether or not a putative or predicted peptide partner will bind. We used known protein–peptide complexes to derive preferences, in the form of spatial position specific scoring matrices, which describe the binding-site environment in globular proteins for each type of amino acid in bound peptides. We then scan the surface of a putative binding protein for sites for each of the amino acids present in a peptide partner and search for combinations of high-scoring amino acid sites that satisfy constraints deduced from the peptide sequence. The method performed well in a benchmark and largely agreed with experimental data mapping binding sites for several recently discovered interactions mediated by peptides, including RG-rich proteins with SMN domains, Epstein-Barr virus LMP1 with TRADD domains, DBC1 with Sir2, and the Ago hook with Argonaute PIWI domain. The method, and associated statistics, is an excellent tool for predicting and studying binding sites for newly discovered peptides mediating critical events in biology.     

基于蛋白质序列的泛素化位点预测研究进展

泛素化是目前广受关注的一种翻译后修饰过程,对蛋白质降解、DNA修复等多种细胞过程都具有重要的调控作用。本文根据国内外蛋白质泛素化位点预测的研究,分析了预测泛素化位点的特征属性,总结了对这些特征进行优化的特征选择方法,并对预测过程中所使用的各种机器学习分类器进行了概述。