Genome-wide identification of cytosine-5 DNA methyltransferases and demethylases in Solanum lycopersicum

Recent studies have reported that decreased level of DNA cytosine methylation in the global genome was closely related to the initiation of tomato (Solanum lycopersicum) fruit ripening. However, genome-scale analysis of cytosine-5 DNA methyltransferases (C5-MTases) and demethylases in tomato has not been engaged. In this study, 7 C5-MTases and 3 demethylases were identified in tomato genome, which probably contributed to DNA cytosine methylation level in tomato. The 7 C5-MTases were categorized into 4 subgroups, and the 3 demethylases were classified into 2 subgroups based on phylogenetic analyses. Comprehensive analysis of their structure and genomic localization was also performed in this paper. According to online RNA-seq data, 4 S. lycopersicum C5-MTase (SlC5-MTase) genes (SlMET, SlDRM1L1, SlDRM5, SlMET3L) were expressed higher than others, and one DNA demethylase gene (SlDML) was significantly changed during tomato fruit development and ripening. Furthermore, all these five gene expressions at breaker (BK) stage changed with 1-methylcyclopropene (1-MCP) treatment, indicating that they were regulated by ethylene directly or indirectly in tomato fruit. In addition, subcellular localization analysis indicated that SlDRM1L1 and SlDRM5 located in the nucleus might have responsibility for RNA-directed DNA methylation (RdDM). Collectively, this paper provided a framework for gene discovery and functional characterization of C5-MTases and DNA demethylases in other Solanaceae species.     

Exact size and organization of DNA target-recognizing domains of multispecific DNA-(cytosine-C5)-methyltransferases.

A large portion of the sequences of type II DNA-(cytosine-C5)-methyltransferases (C5-MTases) represent highly conserved blocks of amino acids. General steps in the methylation reaction performed by C5-MTases have been found to be mediated by some of these domains. C5-MTases carry, in addition at the same relative location, a region variable in size and amino acid composition, part of which is associated with the capacity of each C5-MTase to recognize its characteristic target. Individual target-recognizing domains (TRDs) for the targets CCGG (M), CC(A/T)GG (E), GGCC (H), GCNGC (F) and G(G/A/T)GC(C/A/T)C (B) could be identified in the C-terminal part of the variable region of multispecific C5-MTases. With experiments reported here, we have established the organization of the variable regions of the multispecific MTases M.SPRI, M.phi3TI, M.H2I and M.rho 11SI at the resolution of individual amino acids. These regions comprise 204, 175, 268 and 268 amino acids, respectively. All variable regions are bipartite. They contain at their N-terminal side a very similar sequence of 71 amino acids. The integrity of this sequence must be assured to provide enzyme activity. Bracketed by 6-10 'linker' amino acids, they have, depending on the enzyme studied, towards their C-terminal end ensembles of individual TRDs of 38 (M), 39 (E), 40 (H), 44 (F) and 54 (B) amino acids. TRDs of different enzymes with equal specificity have the same size. TRDs do not overlap but are either separated by linker amino acids or abut each other.     

A directed evolution design of a GCG-specific DNA hemimethylase

DNA cytosine-5 methyltransferases (C5-MTases) are valuable models to study sequence-specific modification of DNA and are becoming increasingly important tools for biotechnology. Here we describe a structure-guided rational protein design combined with random mutagenesis and selection to change the specificity of the HhaI C5-MTase from GCGC to GCG. The specificity change was brought about by a five-residue deletion and introduction of two arginine residues within and nearby one of the target recognizing loops. DNA protection assays, bisulfite sequencing and enzyme kinetics showed that the best selected variant is comparable to wild-type M.HhaI in terms of sequence fidelity and methylation efficiency, and supersedes the parent enzyme in transalkylation of DNA using synthetic cofactor analogs. The designed C5-MTase can be used to produce hemimethylated CpG sites in DNA, which are valuable substrates for studies of mammalian maintenance MTases.     

A second type of N7-guanine RNA cap methyltransferase in an unusual locus of a large RNA virus genome

The order Nidovirales is a diverse group of (+)RNA viruses, with a common genome organization and conserved set of replicative and editing enzymes. In particular, RNA methyltransferases play a central role in mRNA stability and immune escape. However, their presence and distribution in different Nidovirales families is not homogeneous. In Coronaviridae, the best characterized family, two distinct methytransferases perform methylation of the N7-guanine and 2′-OH of the RNA-cap to generate a cap-1 structure (m7GpppNm). The genes of both of these enzymes are located in the ORF1b genomic region. While 2′-O-MTases can be identified for most other families based on conservation of both sequence motifs and genetic loci, identification of the N7-guanine methyltransferase has proved more challenging. Recently, we identified a putative N7-MTase domain in the ORF1a region (N7-MT-1a) of certain members of the large genome Tobaniviridae family. Here, we demonstrate that this domain indeed harbors N7-specific methyltransferase activity. We present its structure as the first N7-specific Rossmann-fold (RF) MTase identified for (+)RNA viruses, making it remarkably different from that of the known Coronaviridae ORF1b N7-MTase gene. We discuss the evolutionary implications of such an appearance in this unexpected location in the genome, which introduces a split-off in the classification of Tobaniviridae.     

Identification of a subdomain within DNA-(cytosine-C5)-methyltransferases responsible for the recognition of the 5' part of their DNA target.

In previous work on DNA-(cytosine-C5)-methyltransferases (C5-MTases), domains had been identified which are responsible for the sequence specificity of the different enzymes (target-recognizing domains, TRDs). Here we have analyzed the DNA methylation patterns of two C5-MTases containing reciprocal chimeric TRDs, consisting of the N- and C-terminal parts derived from two different parental TRDs specifying the recognition of 5'-CC(A/T)GG-3' and 5'-GCNGC-3'. Sequences recognized by these engineered MTases were non-symmetrical and degenerate, but contained at their 5' part a consensus sequence which was very similar to the 5' part of the target recognized by the parental TRD which contributed the N-terminal moiety of the chimeric TRD. The results are discussed in connection with the present understanding of the mechanism of DNA target recognition by C5-MTases. They demonstrate the possibility of designing C5-MTases with novel DNA methylation specificities.     

DNA methyltransferases: mechanistic models derived from kinetic analysis

The sequence-specific transfer of methyl groups from donor S-adenosyl-L-methionine (AdoMet) to certain positions of DNA-adenine or -cytosine residues by DNA methyltransferases (MTases) is a major form of epigenetic modification. It is virtually ubiquitous, except for some notable exceptions. Site-specific methylation can be regarded as a means to increase DNA information capacity and is involved in a large spectrum of biological processes. The importance of these functions necessitates a deeper understanding of the enzymatic mechanism(s) of DNA methylation. DNA MTases fall into one of two general classes; viz. amino-MTases and [C5-cytosine]-MTases. Amino-MTases, common in prokaryotes and lower eukaryotes, catalyze methylation of the exocyclic amino group of adenine ([N6-adenine]-MTase) or cytosine ([N4-cytosine]-MTase). In contrast, [C5-cytosine]-MTases methylate the cyclic carbon-5 atom of cytosine. Characteristics of DNA MTases are highly variable, differing in their affinity to their substrates or reaction products, their kinetic parameters, or other characteristics (order of substrate binding, rate limiting step in the overall reaction). It is not possible to present a unifying account of the published kinetic analyses of DNA methylation because different authors have used different substrate DNAs and/or reaction conditions. Nevertheless, it would be useful to describe those kinetic data and the mechanistic models that have been derived from them. Thus, this review considers in turn studies carried out with the most consistently and extensively investigated [N6-adenine]-, [N4-cytosine]- and [C5-cytosine]-DNA MTases.     

M.phi 3TII: a new monospecific DNA (cytosine-C5) methyltransferase with pronounced amino acid sequence similarity to a family of adenine-N6-DNA-methyltransferases.

The temperate B.subtilis phages phi 3T and rho 11s code, in addition to the multispecific DNA (cytosine-C5) methyltransferases (C5-MTases) M. phi 3TI and M. rho 11sI, which were previously characterized, for the identical monospecific C5-MTases M. phi 3TII and M. rho 11sII. These enzymes modify the C of TCGA sites, a novel target specificity among C5-MTases. The primary sequence of M. phi 3TII (326 amino acids) shows all conserved motifs typical of the building plan of C5-MTases. The degree of relatedness between M. phi 3TII and all other mono- or multispecific C5-MTases ranges from 30-40% amino acid identity. Particularly M. phi 3TII does not show pronounced similarity to M. phi 3TI indicating that both MTase genes were not generated from one another but were acquired independently by the phage. The amino terminal part of the M. phi 3TII (preceding the variable region 'V'), which predominantly constitutes the catalytic domain of the enzyme, exhibits pronounced sequence similarity to the amino termini of a family of A-N6-MTases, which--like M.TaqI--recognize the general sequence TNNA. This suggests that recently described similarities in the general three dimensional organization of C5- and A-N6-MTases imply divergent evolution of these enzymes originating from a common molecular ancestor.     

M.phi 3TII: a new monospecific DNA (cytosine-C5) methyltransferase with pronounced amino acid sequence similarity to a family of adenine-N6-DNA-methyltransferases.

The temperate B.subtilis phages phi 3T and rho 11s code, in addition to the multispecific DNA (cytosine-C5) methyltransferases (C5-MTases) M.phi 3TI and M.rho 11sI, which were previously characterized, for the identical monospecific C5-MTases M.phi 3TII and M.rho 11sII. These enzymes modify the C to TCGA sites, a novel target specificity among C5-MTases. The primary sequence of M.phi 3TII (326 amino acids) shows all conserved motifs typical of the building plan of C5-MTases. The degree of relatedness between M.phi 3TII and all other mono- or multispecific C5-MTases ranges from 30-40% amino acid identity. Particularly M.phi 3TII does not show pronounced similarity to M.phi 3TI indicating that both MTase genes were not generated from one another but were acquired independently by the phage. The amino terminal part of the M.phi 3TII (preceding the variable region 'V'), which predominantly constitutes the catalytic domain of the enzyme, exhibits pronounced sequence similarity to the amino termini of a family of A-N6-MTases, which--like M.Taql--recognize the general sequence TNNA. This suggests that recently described similarities in the general three dimensional organization of C5- and A-N6-MTases imply divergent evolution of these enzymes originating from a common molecular ancestor.     

M.BssHII,a Multispecific Cytosine-C5-DNA- methyltransferase with Unusual Target Recognizing Properties

A new multispecific cytosine-C5-DNA-methyltransferase (C5-MTase), M.BssHII, was identified inBacillus stearothermophilusH3. The M.BssHII gene was cloned and sequenced. The amino acid sequence deduced shows the characteristic building plan of a C5-MTase. By sequencing bisulfite- treated DNA methylated by M.BssHII and by restriction enzyme analysis, we defined the following methylation targets of M.BssHII: ACGCGT/ CCGCGG (MluI/SacII), PuGCGCPy (HaeII), PuCCGGPy (Cfr10I) and GCGCGC (BssHII). The relative location of the specificity determinants in the C5-MTase was derived from the analysis of M.BssHII derivatives carrying deletions within the variable region “V” and chimeric C5-Mtases constructed between M.BssHII and the related monospecific enzyme M.ψ3TII. Four of the M.BssHII specificities (MluI,SacII,Cfr10I andBssHII) could be associated with amino acid segments within the variable region “V”. The determinant forHaeII activity had to be assigned to sequences defining the enzyme core, the first example of a C5-MTase in which a sequence-specific methylation potential is mediated by structures outside of the variable region. Another intriguing result came from the analysis of one particular chimera made between M.BssHII and M.ψ3TII. This construct showed a relaxation of the methylation capacity, both with respect to the target recognized and the targeting of methylation within this sequence.     

Inhibition of HhaI DNA (Cytosine-C5) methyltransferase by oligodeoxyribonucleotides containing 5-aza-2'-deoxycytidine: examination of the intertwined roles of co-factor,target, transition state structure and enzyme conformation

The presence of 5-azacytosine (ZCyt) residues in DNA leads to potent inhibition of DNA (cytosine-C5) methyltranferases (C5-MTases) in vivo and in vitro. Enzymatic methylation of cytosine in mammalian DNA is an epigenetic modification that can alter gene activity and chromosomal stability, influencing both differentiation and tumorigenesis. Thus, it is important to understand the critical mechanistic determinants of ZCyt's inhibitory action. Although several DNA C5-MTases have been reported to undergo essentially irreversible binding to ZCyt in DNA, there is little agreement as to the role of AdoMet and/or methyl transfer in stabilizing enzyme interactions with ZCyt. Our results demonstrate that formation of stable complexes between HhaI methyltransferase (M.HhaI) and oligodeoxyribonucleotides containing ZCyt at the target position for methylation (ZCyt-ODNs) occurs in both the absence and presence of co-factors, AdoMet and AdoHcy. Both binary and ternary complexes survive SDS-PAGE under reducing conditions and take on a compact conformation that increases their electrophoretic mobility in comparison to free M.HhaI. Since methyl transfer can occur only in the presence of AdoMet, these results suggest (1) that the inhibitory capacity of ZCyt in DNA is based on its ability to induce a stable, tightly closed conformation of M.HhaI that prevents DNA and co-factor release and (2) that methylation of ZCyt in DNA is not required for inhibition of M.HhaI.     

植物miR159家族成员分子特性及其进化规律研究

为了解植物miR159家族成员的分子特性及其进化规律,该研究对miRBase数据库中登录的miR159家族成员进行分类统计、进化树构建、科间比较、二级结构预测及靶基因分析。结果表明:miR159家族在植物界分布非常广泛,蕨类植物可能是miR159家族的进化祖先;系统发育进化树分析显示,植物miR159家族成员间存在多个进化分支,且进化关系与植物属性有关,即植物亲缘关系越近的成员更易成枝,且具有相同进化方向的成员序列高度同源;Mfold预测显示,pre-miR159均会自发形成典型、稳定的茎环二级结构,并包含19~21个碱基为单位的功能片段,每个单位均有可能形成miR159成熟体;靶基因分析发现,miR159家族成员主要作用于MYB转录因子、转座因子和假定蛋白等,但在不同物种间及相同物种的不同成员间作用的靶标种类与靶基因ID数量均存在差异,尤其是miR159-3p与miR159-5p间的差异最为显著。     

Dicotyledons of cretaceous, palaeogene, and neogene. Adaptogenesis of the terminal phloem

A database containing information on the trait complex for terminal phloem of dicots leaf having been created during the last thirty years was analyzed. The typology of terminal phloem was verified. The monotypy of plant taxa up to the family level was confirmed. For orders and subclasses, the wide distribution of structural parallelisms was demonstrated. The lists of plant families with similar structures of the phloem terminals was elaborated. As a result of comparative analysis of these lists and paleobotanic data on the age of plant families, the specificity of the types of terminal phloem and the associated morphological and functional characteristics of the groups of families of the same evolutionary, age were revealed. Ancestral type prevails among the families belonging to Prodicots, which were emerged during Cretaceous. Symplastic type is characteristic to the families of Eudicots and biomes appeared in Paleogene, while the apoplastic type is characteristic to the families and biomes of Neogene. For each group to taxa, plant stem and root parasitism is characteristic for the crown taxa completing the structural sequences of trees and herbs. All these data allow discussion of the questions of climatic adaptogenesis of the character complex studied here, the impossibility of structural iterations in morphological evolution, and the reasons of similarities and differences in the evolutionary and contemporary zonal series of the terminal phloem in Dicots.     

Tree population studies in low-diversity forests,Guyana. I. Floristic composition and stand structure

Studies were undertaken on the floristic composition and stand structure of four 1 hectare plots in the lowland forests of Kurupukari, Guyana. A total of 3897 trees, covering 153 species and 31 plant families were recorded at greater than 5 cm diameter at breast height (dbh). The number of species per hectare ranged from 61 to 84 (>5.0 cm dbh) and 50–71 (>10.0 cm dbh). The total number of trees per hectare varied two-fold between study plots, with 45–50% of the trees within the 5–10 cm size-class. Mean total basal area varied from 32.39–34.63 m² per 100 m². The four most dominant plant families represented 43.8% of the total number of trees, while representing only 11.2% of the species. No one plant family dominated in more than one of the four study plots, and all four plots held at least one plant family with more than 20% of the total number of trees. Although 14 tree species were common to all four plots, only 26%–35% of the species were represented by a single tree. Between three and seven species represented 50% of the trees within all size-classes, with species dominance occurring within the highest density plot.These tropical forest types of central Guyana may represent some of the lowest diversity forests in the neotropics, whereby the total number of tree species is relatively limited, typically with six dominant canopy species, but the relative abundance of these species is highly variable between the forest types. Mechanisms influencing the competitive interactions associated with species dominance are discussed in relation to the importance of mycorrhizae and the persistence of species dominance.     

Sequence and Phylogenetic Analyses of 4 TMS Junctional Proteins of Animals: Connexins,Innexins, Claudins and Occludins

Connexins and probably innexins are the principal constituents of gap junctions, while claudins and occludins are principal tight junctional constituents. All have similar topologies with four alpha-helical transmembrane segments (TMSs), and all exhibit well-conserved extracytoplasmic cysteines that either are known to or potentially can form disulfide bridges. We have conducted sequence, topological and phylogenetic analyses of the proteins that comprise the connexin, innexin, claudin and occludin families. A multiple alignment of the sequences of each family was used to derive average hydropathy and similarity plots as well as phylogenetic trees. Analyses of the data generated led to the following evolutionary and functional suggestions: (1) In all four families, the most conserved regions of the proteins from each family are the four TMSs although the extracytoplasmic loops between TMSs 1 and 2, and TMSs 3 and 4 are usually well conserved. (2) The phylogenetic trees revealed sets of orthologues except for the innexins where phylogeny primarily reflects organismal source, probably due to a lack of relevant organismal sequence data. (3) The two halves of the connexins exhibit similarities suggesting that they were derived from a common origin by an internal gene duplication event. (4) Conserved cysteyl residues in the connexins and innexins may point to a similar extracellular structure involved in the docking of hemichannels to create intercellular communication channels. (5) We suggest a similar role in homomeric interactions for conserved extracellular residues in the claudins and occludins. The lack of sequence or motif similarity between the four different families indicates that, if they did evolve from a common ancestral gene, they have diverged considerably to fulfill separate, novel functions. We suggest that internal duplication was a general evolutionary strategy used to generate new families of channels and junctions with unique functions. These findings and suggestions should serve as guides for future studies concerning the structures, functions and evolutionary origins of junctional proteins.     

Morphological types of arbuscular mycorrhizas in pioneer woody plants growing in an oil palm farm in Sumatra, Indonesia

Morphological types of arbuscular mycorrhizas (AM) in pioneer woody plants were examined in an oil palm farm 1 year after reclamation. In total, 18 plant species belonging to 15 genera in 9 families were examined, and the Arum type was found in 12 species belonging to 9 genera in 5 families. In contrast, the Paris type was found in 5 species belonging to 5 genera in 4 families. The number of plant species with Arum-type AM was always higher than that with Paris-type AM in all four sampling plots. All the examined species in Euphorbiaceae, most of which are fast-growing pioneer trees and were dominant in forming the canopy in secondary forests, were found to have Arum-type AM. These results suggest that Arum-type colonization is beneficial for fast-growing woody plant species.     

蔷薇科果树类受体激酶的研究进展

类受体激酶(receptor like kinase,RLK)参与调控植物几乎所有的生命活动,是植物生长发育和环境适应的“中央处理器”。该文对近年来国内外有关蔷薇科果树RLK基因鉴定、进化特征及其在各器官生长发育、非生物和生物逆境中的作用及调控机制等方面的研究进展进行了综述。蔷薇科果树基因组中存在数目庞大的RLKs,不同树种间的RLK数目和各亚家族成员数目都存在较大差异,而且蔷薇科果树RLK存在极为普遍的部分重复和串联重复现象,是导致家族成员迅速变化的重要原因。有研究发现,一些RLKs调控蔷薇科果树器官发育和对环境的适应性。在器官发育方面,LRR RLK亚家族成员调控根系发育,CrRLK1L、LysM RLK和LRR RLK亚家族部分成员参与调控果实发育,CrRLK1L亚家族成员参与调控花粉管发育,LRR RLK、LysM RLK、L LEC RLK和B Lectin RLK亚家族部分成员调控蔷薇科果树对生物逆境的适应。今后RLK功能研究可侧重于蔷薇科果树特色性状,通过提高目标基因的筛选和验证的效率,加速主效RLKs的筛选进程,并通过筛选主效RLKs诱导方式和加速分子育种进程等途径,将研究成果应用于实际生产。     

Global angiosperm family richness revisited: linking ecology and evolution to climate

Aim The global richness gradient of angiosperm families is correlated with current climate, and it has been claimed that historical processes are not necessary to understand patterns of plant family richness. This claim has drawn criticism, and there have been doubts about the quality of the data used to quantify the pattern. We revisit this issue using the Angiosperm Phylogeny Group (APG) III classification and revised range maps, and we incorporate an evolutionary variable, family age, to explore covariation between evolution and ecology and their links to climate via the tropical conservatism hypothesis (TCH). Location Global. Methods The richness pattern for 408 families was derived from range maps, and family ages were derived from a dated angiosperm phylogeny. Patterns were generated for all families, 143 families composed of trees, and 149 families composed of herbs. We also examined family range size patterns to test the extent to which extratropical floras are nested subsets of tropical floras. Ordinary least squares (OLS) multiple and partial regressions were used to generate climate models for richness, mean range size and mean age for each plant dataset and to evaluate the covariation between contemporary climate and clade age as correlates of family richness. Results We confirmed the strong association between contemporary climate and family richness. Age patterns predicted by TCH were also found for families comprising trees. The richness of herbaceous families, in contrast, was correlated with climate but the age pattern was not as predicted by TCH. Floras in cold and dry areas are strongly nested within richer tropical floras. Main conclusions Phylogenetic niche conservatism at the family level offers a likely explanation for the global diversity gradient of trees, but not for non‐desert herbs, probably because of the faster evolutionary rates for herbs and less constrained evolutionary responses to climate change. Thus, it appears that multiple processes account for the overall angiosperm family gradient. Our analysis also demonstrates that even very strong associations of taxon richness and climate do not preclude evolutionary processes, as has been widely argued, and that climatic and evolutionary hypotheses for richness gradients are not mutually exclusive.     

Genotype and environment interact to influence acceptability and suitability of white spruce for a specialist herbivore, Zeiraphera canadensis

Abstract.

• 1 The independent and interacting effects of plant genotype and site (i.e. environment) on the acceptability of white spruce, Picea gluaca (Moench) Voss, to the spruce bud moth, Zeiraphera canadensis Mut. & Free. and on plant suitability for egg development, were studied at four sites in New Brunswick, Canada.
• 2 A greater proportion of shoots on trees in two half-sib families, previously designated as highly susceptible, were partially eaten by spruce bud moths than shoots on trees in two half-sib families with low susceptibility.
• 3 At the site with the highest bud moth population, oviposition was highest on trees in susceptible families and on branches damaged by bud moth larvae. Oviposition was not higher on trees in susceptible families at the other three sites, resulting in a strong tree genotype × site interaction for oviposition.
• 4 Although there was a significant tree genotype × site interaction for egg predation, egg survival was higher on trees in susceptible families at all sites, due to lower levels of egg parasitism and predation.
• 5 Egg densities were positively but weakly correlated to shoot length and diameter. There were no consistent relationships between shoot length, shoot diameter, needle length or needle density and per cent egg survival, parasitism or predation.
• 6 Because egg survival was higher on trees in susceptible families at all sites but egg densities were only higher on trees in susceptible families at one site, host plant acceptability and suitability were positively related at only one site. We speculate that Z.canadensis does not distinguish between hosts of different suitability until trees are heavily damaged.