蔷薇类的分子系统学研究—兼论系统发育基因组学的原理和方法(以SET基因超家族为例)

分子系统发育分析的主要任务包括：（1）帮助建立生命之树（tree of life）;（2）追踪基因和基因家族（gene family）的起源和进化, 以获知基因在进化过程中的功能分化和伴随发生的重要分子事件（key molecular events）和形态性状的关键创新（key innovation）。这两个方面在本研究中都有所涉及。对于前者，选用植物线粒体matR基因重建被子植物蔷薇类群的系统发育关系;对于后者，则以SET基因超家族为例，探讨其在真核生物中的进化分类以及与功能多样性的关系。 I 蔷薇类的分子系统学 蔷薇类（rosids）是基于分子数据建立的被子植物的主要分支之一，包含13个目，大约三分之一的被子植物物种。两个主要蔷薇类内部分支是豆类fabids（包含7个目）和锦葵类malvids（包含3个目）。然而，这两个分支内部，以及这两个分支与蔷薇类基部类群，包括牻牛儿苗目（Geraniales）、桃金娘目（Myrtales）和流苏子目（Crossosomatales）之间的关系大多是不清楚的。本研究中，我们选取174个物种来代表72个蔷薇类（rosids）的科，利用两个数据集，即线粒体matR单基因数据集和包括线粒体matR基因、两个质体基因(rbcL、 atpB)和一个核基因(18S rDNA) 的4基因数据集，重建蔷薇类在科以上分类阶元水平的系统发育关系。同时，还对线粒体matR基因的进化特征和用于大尺度系统发育分析的适合度和潜力进行了评价。 线粒体matR单基因数据支持malvids和大多数蔷薇类目的单系性质，然而，豆类（fabids）成员没有形成一个分支，其COM亚支，包括卫矛目（Celastrales）、酢浆草目（Oxalidales）、金虎尾目（Malpighiales）和蒜树科（Huaceae），分辨为锦葵类（malvids）的姐妹群。这个关系在最近根据花结构特征曾被提出过，但从未在之前的分子系统发育分析中得到分辨。4基因数据集支持首先是牻牛儿苗目（Geraniales），接着是桃金娘目（Myrtales）作为蔷薇类（rosids）的最基部的分支;流苏子目（Crossosomatales）是锦葵类（malvids）姐妹群，以及蔷薇类（rosids）的核心部分包括豆类（fabids），锦葵类（malvids）和流苏子目（Crossosomatales）。线粒体matR基因的进化特征分析显示，与两个叶绿体基因(rbcL 和atpB)比较，同义替代速率约是它们的1/4，而非同义替代速率接近于自身的同义替代速率，表明matR 基因具有松弛的选择压力。线粒体matR基因相对慢速的进化使非同源相似（homoplasious）突变减少，提高了系统发育信息的质量，同时，松弛的选择压力使非同义替代数量增加，弥补了慢速进化导致的系统发育信息数量不足的缺陷，这两个方面的结合使线粒体matR基因非常适用于被子植物在科以上水平的系统发育研究。 II SET基因超家族的系统发育基因组学分析 SET基因超家族基因编码含有SET结构域的蛋白，在真核生物中，SET-domain蛋白一般是多结构域（multi-domain）的。SET-domain蛋白具有对组蛋白H3和H4的N末端尾部进行赖氨酸残基甲基化修饰的酶活性;从异染色质形成到基因转录，甲基化的组蛋白广泛影响染色质水平的基因调控。依据SET结构域一级序列的相似性和结构域组织（domain architecture）特征，目前，SET-domain基因超家族被划分为4-7个家族。由于这些划分或者使用动物或者使用植物SET基因，只有少数其它类群的物种加入分析，因此这样的划分可能是不完整的。本研究采用系统发育基 因组学方法（phylogenomic approach），在真核生物范围内广泛取样，期望获得相对完整的SET-domain基因家族的 进化分类方案，在此基础上加深理解SET-domain基因的进化机制和功能多样性。 在提取了17个物种，代表5个真核超群的SET蛋白序列基础上，系统发育分析结合“结构域组织特征”鉴别了9个SET基因家族，其中一个是新的SET基因家族。以前的SET8和Class VI家族,及SMYD和SUV4-20家族分别合并为一个家族。大部分家族在进化过程中发生了2次以上的基因重复事件，通过获得不同的结构域产生具有不同功能的新基因。一个SET基因家族在进化过程中推测发生了从脊椎动物祖先向盘基网柄菌（Dictyostelium discoideum）的水平基因转移。

Molecular studies in gorgonian alloimmunity : search for gene homologs of the immunoglobulin gene superfamily in Swiftia exserta

Humoral and cells surface molecules of the mammalian immune system, grouped into the Immunoglobulin Gene Superfamily, share protein structure and gene sequence homologies with molecules found among diverse phylogenetic groups. In histocompatibility studies, the gorgonian coral Swiftia exserta has recently demonstrated specific alloimmunity with memory (Salter-Cid and Bigger, 1991. Biological Bulletin Vol 181). In an attempt to shed light on the origins of this gene family and the evolution of the vertebrate immune response, genomic DNA from Swiftia exserta was isolated, purified, and analyzed by Southern blot hybridization with mouse gene probes corresponding to two molecules of the Immunoglobulin Gene Superfamily, the Thy-1 antigen, and the alpha-3 domain of the MHC Class I histocompatibility marker. Hybridizations were conducted under low to non-stringent conditions to allow binding of mismatched homologs that may exist between the mouse gene probes and the Swiftia DNA. Removal of non-specific binding (sequences less than 70% homologous) occurred in washing steps. Results show that with the probes selected, the method chosen, and the conditions applied, no evidence of sequences of 70% or greater homology to the mouse Thy-1 or MHC Class I alpha-3 genes exist in Swiftia exserta genome.

Ensuring the statistical soundness of competitive gene set approaches: gene filtering and genome-scale coverage are essential

In this article, we focus on the analysis of competitive gene set methods for detecting the statistical significance of pathways from gene expression data. Our main result is to demonstrate that some of the most frequently used gene set methods, GSEA, GSEArot and GAGE, are severely influenced by the filtering of the data in a way that such an analysis is no longer reconcilable with the principles of statistical inference, rendering the obtained results in the worst case inexpressive. A possible consequence of this is that these methods can increase their power by the addition of unrelated data and noise. Our results are obtained within a bootstrapping framework that allows a rigorous assessment of the robustness of results and enables power estimates. Our results indicate that when using competitive gene set methods, it is imperative to apply a stringent gene filtering criterion. However, even when genes are filtered appropriately, for gene expression data from chips that do not provide a genome-scale coverage of the expression values of all mRNAs, this is not enough for GSEA, GSEArot and GAGE to ensure the statistical soundness of the applied procedure. For this reason, for biomedical and clinical studies, we strongly advice not to use GSEA, GSEArot and GAGE for such data sets.

Genome dynamics explain the evolution of flowering time CCT domain gene families in the Poaceae

Numerous CCT domain genes are known to control flowering in plants. They belong to the CONSTANS-like (COL) and PREUDORESPONSE REGULATOR (PRR) gene families, which in addition to a CCT domain possess B-box or response-regulator domains, respectively. Ghd7 is the most recently identified COL gene to have a proven role in the control of flowering time in the Poaceae. However, as it lacks B-box domains, its inclusion within the COL gene family, technically, is incorrect. Here, we show Ghd7 belongs to a larger family of previously uncharacterized Poaceae genes which possess just a single CCT domain, termed here CCT MOTIF FAMILY (CMF) genes. We molecularly describe the CMF (and related COL and PRR) gene families in four sequenced Poaceae species, as well as in the draft genome assembly of barley (Hordeum vulgare). Genetic mapping of the ten barley CMF genes identified, as well as twelve previously unmapped HvCOL and HvPRR genes, finds the majority map to colinear positions relative to their Poaceae orthologues. Combined inter-/intra-species comparative and phylogenetic analysis of CMF, COL and PRR gene families indicates they evolved prior to the monocot/dicot divergence ~200 mya, with Poaceae CMF evolution described as the interplay between whole genome duplication in the ancestral cereal, and subsequent clade-specific mutation, deletion and duplication events. Given the proven role of CMF genes in the modulation of cereals flowering, the molecular, phylogenetic and comparative analysis of the Poaceae CMF, COL and PRR gene families presented here provides the foundation from which functional investigation can be undertaken.

Systematic Reverse Genetics of Transfer-DNA-Tagged Lines of Arabidopsis1 : Isolation of Mutations in the Cytochrome P450 Gene Superfamily

We have developed an efficient reverse-genetics protocol that uses expedient pooling and hybridization strategies to identify individual transfer-DNA insertion lines from a collection of 6000 independently transformed lines in as few as 36 polymerase chain reactions. We have used this protocol to systematically isolate Arabidopsis lines containing insertional mutations in individual cytochrome P450 genes. In higher plants P450 genes encode enzymes that perform an exceptionally wide range of functions, including the biosynthesis of primary metabolites necessary for normal growth and development, the biosynthesis of secondary products, and the catabolism of xenobiotics. Despite their importance, progress in assigning enzymatic function to individual P450 gene products has been slow. Here we report the isolation of the first 12 such lines, including one (CYP83B1-1) that displays a runt phenotype (small plants with hooked leaves), and three insertions in abundantly expressed genes. The DNAs used in this study are publicly available and can be used to systematically isolate mutants in Arabidopsis.

Targeted deletion of the mouse POU domain gene Brn-3a causes selective loss of neurons in the brainstem and trigeminal ganglion, uncoordinated limb movement, and impaired suckling.

The Brn-3 subfamily of POU domain genes are expressed in sensory neurons and in select brainstem nuclei. Earlier work has shown that targeted deletion of the Brn-3b and Brn-3c genes produce, respectively, defects in the retina and in the inner ear. We show herein that targeted deletion of the Brn-3a gene results in defective suckling and in uncoordinated limb and trunk movements, leading to early postnatal death. Brn-3a (-/-) mice show a loss of neurons in the trigeminal ganglia, the medial habenula, the red nucleus, and the caudal region of the inferior olivary nucleus but not in the retina and dorsal root ganglia. In the trigeminal and dorsal root ganglia, but not in the retina, there is a marked decrease in the frequency of neurons expressing Brn-3b and Brn-3c, suggesting that Brn-3a positively regulates Brn-3b and Brn-3c expression in somatosensory neurons. Thus, Brn-3a exerts its major developmental effects in somatosensory neurons and in brainstem nuclei involved in motor control. The pheno-types of Brn-3a, Brn-3b, and Brn-3c mutant mice indicate that individual Brn-3 genes have evolved to control development in the auditory, visual, or somatosensory systems and that despite differences between these systems in transduction mechanisms, sensory organ structures, and central information processing, there may be fundamental homologies in the genetic regulatory events that control their development.

Analysis of MLL2 gene in the first Brazilian family with Kabuki syndrome

Most patients with Kabuki syndrome (KS) are the only person in their family with the condition. However, familial cases of KS have been described showing evidence that this syndrome can be inherited as a dominant trait with variable expressivity. We report on two related individuals with facial findings characteristic of KS. The proposita had arched eyebrows, long and upward slanting palpebral fissures, cleft lip and palate, retromicrognathia, brachydactyly of hands and feet, stubby fingers, nail hypoplasia, and prominent finger pads. Her mother had eyebrows with dispersed lateral half, long and upward slanting palpebral fissures, retrognathia, abnormal and posteriorly rotated ears, prominent finger pads, brachydactyly of feet, learning difficulties, and psychomotor development delay. DNA sequencing revealed a novel missense mutation in the MLL2 gene in both the proposita and her mother. The mutation (p.R5432Q) was found in the exon 51, within the SET domain of the gene, which confers methyltransferase activity on the protein. Therefore, the epigenetic and transcriptional regulatory properties of this protein may be altered and this suggests that the mutation is the cause of phenotype observed in both the patient and her mother. The clinical signs and the molecular evidence in this family further support the notion that KS is an autosomal dominant condition with variable expressivity. To our knowledge this is the first report of a Brazilian family with recurrence of this syndrome. (C) 2012 Wiley Periodicals, Inc.

Control of fertilization-independent endosperm development by the MEDEA polycomb gene in Arabidopsis

Higher plant reproduction is unique because two cells are fertilized in the haploid female gametophyte. Egg and sperm nuclei fuse to form the embryo. A second sperm nucleus fuses with the central cell nucleus that replicates to generate the endosperm, a tissue that supports embryo development. To understand mechanisms that initiate reproduction, we isolated a mutation in Arabidopsis, f644, that allows for replication of the central cell and subsequent endosperm development without fertilization. When mutant f644 egg and central cells are fertilized by wild-type sperm, embryo development is inhibited, and endosperm is overproduced. By using a map-based strategy, we cloned and sequenced the F644 gene and showed that it encodes a SET-domain polycomb protein. Subsequently, we found that F644 is identical to MEDEA (MEA), a gene whose maternal-derived allele is required for embryogenesis [Grossniklaus, U., Vielle-Calzada, J.-P., Hoeppner, M. A. & Gagliano, W. B. (1998) Science 280, 446–450]. Together, these results reveal functions for plant polycomb proteins in the suppression of central cell proliferation and endosperm development. We discuss models to explain how polycomb proteins function to suppress endosperm and promote embryo development.

The C-terminal SET domains of ALL-1 and TRITHORAX interact with the INI1 and SNR1 proteins, components of the SWI/SNF complex

The ALL-1 gene was discovered by virtue of its involvement in human acute leukemia. Its Drosophila homolog trithorax (trx) is a member of the trx-Polycomb gene family, which maintains correct spatial expression of the Antennapedia and bithorax complexes during embryogenesis. The C-terminal SET domain of ALL-1 and TRITHORAX (TRX) is a 150-aa motif, highly conserved during evolution. We performed yeast two hybrid screening of Drosophila cDNA library and detected interaction between a TRX polypeptide spanning SET and the SNR1 protein. SNR1 is a product of snr1, which is classified as a trx group gene. We found parallel interaction in yeast between the SET domain of ALL-1 and the human homolog of SNR1, INI1 (hSNF5). These results were confirmed by in vitro binding studies and by demonstrating coimmunoprecipitation of the proteins from cultured cells and/or transgenic flies. Epitope-tagged SNR1 was detected at discrete sites on larval salivary gland polytene chromosomes, and these sites colocalized with around one-half of TRX binding sites. Because SNR1 and INI1 are constituents of the SWI/SNF complex, which acts to remodel chromatin and consequently to activate transcription, the interactions we observed suggest a mechanism by which the SWI/SNF complex is recruited to ALL-1/trx targets through physical interactions between the C-terminal domains of ALL-1 and TRX and INI1/SNR1.

Plastid-localized acetyl-CoA carboxylase of bread wheat is encoded by a single gene on each of the three ancestral chromosome sets

5′-End fragments of two genes encoding plastid-localized acetyl-CoA carboxylase (ACCase; EC 6.4.1.2) of wheat (Triticum aestivum) were cloned and sequenced. The sequences of the two genes, Acc-1,1 and Acc-1,2, are 89% identical. Their exon sequences are 98% identical. The amino acid sequence of the biotin carboxylase domain encoded by Acc-1,1 and Acc-1,2 is 93% identical with the maize plastid ACCase but only 80–84% identical with the cytosolic ACCases from other plants and from wheat. Four overlapping fragments of cDNA covering the entire coding region were cloned by PCR and sequenced. The wheat plastid ACCase ORF contains 2,311 amino acids with a predicted molecular mass of 255 kDa. A putative transit peptide is present at the N terminus. Comparison of the genomic and cDNA sequences revealed introns at conserved sites found in the genes of other plant multifunctional ACCases, including two introns absent from the wheat cytosolic ACCase genes. Transcription start sites of the plastid ACCase genes were estimated from the longest cDNA clones obtained by 5′-RACE (rapid amplification of cDNA ends). The untranslated leader sequence encoded by the Acc-1 genes is at least 130–170 nucleotides long and is interrupted by an intron. Southern analysis indicates the presence of only one copy of the gene in each ancestral chromosome set. The gene maps near the telomere on the short arm of chromosomes 2A, 2B, and 2D. Identification of three different cDNAs, two corresponding to genes Acc-1,1 and Acc-1,2, indicates that all three genes are transcriptionally active.

Probing the transmembrane topology of cyclic nucleotide-gated ion channels with a gene fusion approach.

Cyclic nucleotide-gated (CNG) cation channels contain two short sequence motifs--a residual voltage-sensor (S4) and a pore-forming (P) segment--that are reminiscent of similar segments in voltage-activated Shaker-type K+ channels. It has been tacitly assumed that CNG channels and this K+ channel subfamily share a common overall topology, characterized by a hydrophobic domain comprising six membrane-spanning segments. We have systematically investigated the topology of CNG channels from bovine rod photoreceptor and Drosophila melanogaster by a gene fusion approach using the bacterial reporter enzymes alkaline phosphatase and beta-galactosidase, which are active only in the periplasm and only in the cytoplasm, respectively. Enzymatic activity was determined after expression of fusion constructs in Escherichia coli. CNG channels were found to have six membrane-spanning segments, suggesting that CNG and Shaker-type K+ channels, albeit distant relatives within a gene superfamily of ion channels, share a common topology.

G9a histone methyltransferase contributes to imprinting in the mouse placenta

Whereas DNA methylation is essential for genomic imprinting, the importance of histone methylation in the allelic expression of imprinted genes is unclear. Imprinting control regions (ICRs), however, are marked by histone H3-K9 methylation on their DNA-methylated allele. In the placenta, the paternal silencing along the Kcnq1 domain on distal chromosome 7 also correlates with the presence of H3-K9 methylation, but imprinted repression at these genes is maintained independently of DNA methylation. To explore which histone methyltransferase (HMT) could mediate the allelic H3-K9 methylation on distal chromosome 7, and at ICRs, we generated mouse conceptuses deficient for the SET domain protein G9a. We found that in the embryo and placenta, the differential DNA methylation at ICRs and imprinted genes is maintained in the absence of G9a. Accordingly, in embryos, imprinted gene expression was unchanged at the domains analyzed, in spite of a global loss of H3-K9 dimethylation (H3K9me2). In contrast, the placenta-specific imprinting of genes on distal chromosome 7 is impaired in the absence of G9a, and this correlates with reduced levels of H3K9me2 and H3K9me3. These findings provide the first evidence for the involvement of an HMT and suggest that histone methylation contributes to imprinted gene repression in the trophoblast.

A biochemical and structural characterization of Drosophila neuroglian

A number of cell-cell interactions in the nervous system are mediated by immunoglobulin gene superfamily members. For example, neuroglian, a homophilic neural cell adhesion molecule in Drosophila, has an extracellular portion comprising six C- 2 type immunoglobulin-like domains followed by five fibronectin type III (FnIII) repeats. Neuroglian shares this domain organization and significant sequence identity with Ll, a murine neural adhesion molecule that could be a functional homologue. Here I report the crystal structure of a proteolytic fragment containing the first two FnIII repeats of neuroglian (NgFn 1,2) at 2.0Å. The interpretation of photomicrographs of rotary shadowed Ng, the entire extracellular portion of neuroglian, and NgFnl-5, the five neuroglian Fn III domains, is also discussed.

The structure of NgFn 1,2 consists of two roughly cylindrical β-barrel structural motifs arranged in a head-to-tail fashion with the domains meeting at an angle of ~120, as defined by the cylinder axes. The folding topology of each domain is identical to that previously observed for single FnIII domains from tenascin and fibronectin. The domains of NgFn1,2 are related by an approximate two fold screw axis that is nearly parallel to the longest dimension of the fragment. Assuming this relative orientation is a general property of tandem FnIII repeats, the multiple tandem FnIII domains in neuroglian and other proteins are modeled as thin straight rods with two domain zig-zag repeats. When combined with the dimensions of pairs of tandem immunoglobulin-like domains from CD4 and CD2, this model suggests that neuroglian is a long narrow molecule (20 - 30 Å in diameter) that extends up to 370Å from the cell surface.

In photomicrographs, rotary shadowed Ng and NgFn1-5 appear to be highly flexible rod-like molecules. NgFn 1-5 is observed to bend in at least two positions and has a mean total length consistent with models generated from the NgFn1,2 structure. Ng molecules have up to four bends and a mean total length of 392 Å, consistent with a head-to-tail packing of neuroglian's C2-type domains.

Rubisco in complex with Rubisco large subunit methyltransferase.

SET domain protein lysine methyltransferases (PKMT) are a structurally unique class of enzymes that catalyze the specific methylation of lysine residues in a number of different substrates. Especially histone-specific SET domain PKMTs have received widespread attention because of their roles in the regulation of epigenetic gene expression and the development of some cancers. Rubisco large subunit methyltransferase (RLSMT) is a chloroplast-localized SET domain PKMT responsible for the formation of trimethyl-lysine-14 in the large subunit of Rubisco, an essential photosynthetic enzyme. Here, we have used cryoelectron microscopy to produce an 11-A density map of the Rubisco-RLSMT complex. The atomic model of the complex, obtained by fitting crystal structures of Rubisco and RLSMT into the density map, shows that the extensive contact regions between the 2 proteins are mainly mediated by hydrophobic residues and leucine-rich repeats. It further provides insights into potential conformational changes that may occur during substrate binding and catalysis. This study presents the first structural analysis of a SET domain PKMT in complex with its intact polypeptide substrate.

Involvement of histone H3 methylation in Tup1-Ssn6 repression

The Tup1-Ssn6 complex regulates the expression of diverse classes of genes in Saccharomyces cerevisiae including those regulated by mating type, DNA damage, glucose, and anaerobic stress. The complex is recruited to target genes by sequence-specific repressor proteins. Once recruited to particular promoters, it is not completely clear how it functions to block transcription. Repression probably occurs through interactions with both the basal transcriptional machinery and components of chromatin. Tup1 interactions with chromatin are strongly influenced by acetylation of histories H3 and H4. Tup1 binds to underacetylated histone tails and requires multiple histone deacetylases (HDACs) for its repressive functions. Like acetylation, histone methylation is involved in regulation of gene expression. The possible role of histone methylation in Tup1 repression is not known. Here we examined possible roles of histone methyltransferases in Tup1-Ssn6 functions. We found that like other genes, Tup1-Ssn6 target genes exhibit increases in the levels of histone H3 lysine 4 methylation upon activation. However, deletion of individual or multiple histone methyltransferases (HMTs) and other SET-domain containing genes has no apparent effect on Tup1-Ssn6 mediated repression of a number of well-defined targets. Interestingly, we discovered that Ssn6 interacts with Set2. Since deletion of SET2 does not affect Tup1-Ssn6 repression, Ssn6 may utilize Set2 in other contexts to regulate gene repression. In order examine if the two components of the Tup1-Ssn6 complex have independent functions in the cell, we identified genes differentially expressed in tup1Δ and ssn6Δ mutants using DNA microarrays. Our data indicate that ∼4% of genes in the cell are regulated by Ssn6 independently of Tup1. In addition, expression of genes regulated by Tup1-Ssn6 seems to be differently affected by deletion of Ssn6 and deletion of Tup1, which indicates that these components might have separate functions. Our data shed new light on the classical view of Tup1-Ssn6 functions, and indicate that Ssn6 might have repressive functions as well. ^