From protein microarray to biology: The discovery of epigenetic regulatory molecules

Chromatin, composed of repeating nucleosome units, is the genetic polymer of life. To aid in DNA compaction and organized storage, the double helix wraps around a core complex of histone proteins to form the nucleosome, and is therefore no longer freely accessible to cellular proteins for the processes of transcription, replication and DNA repair. Over the course of evolution, DNA-based applications have developed routes to access DNA bound up in chromatin, and further, have actually utilized the chromatin structure to create another level of complexity and information storage. The histone molecules that DNA surrounds have free-floating tails that extend out of the nucleosome. These tails are post-translationally modified to create docking sites for the proteins involved in transcription, replication and repair, thus providing one prominent way that specific genomic sequences are accessed and manipulated. Adding another degree of information storage, histone tail-modifications paint the genome in precise manners to influence a state of transcriptional activity or repression, to generate euchromatin, containing gene-dense regions, or heterochromatin, containing repeat sequences and low-density gene regions. The work presented here is the study of histone tail modifications, how they are written and how they are read, divided into two projects. Both begin with protein microarray experiments where we discover the protein domains that can bind modified histone tails, and how multiple tail modifications can influence this binding. Project one then looks deeper into the enzymes that lay down the tail modifications. Specifically, we studied histone-tail arginine methylation by PRMT6. We found that methylation of a specific histone residue by PRMT6, arginine 2 of H3, can antagonize the binding of protein domains to the H3 tail and therefore affect transcription of genes regulated by the H3-tail binding proteins. Project two focuses on a protein we identified to bind modified histone tails, PHF20, and was an endeavor to discover the biological role of this protein. Thus, in total, we are looking at a complete process: (1) histone tail modification by an enzyme (here, PRMT6), (2) how this and other modifications are bound by conserved protein domains, and (3) by using PHF20 as an example, the functional outcome of binding through investigating the biological role of a chromatin reader. ^

Identification of a conserved cluster in the RH domain of GRK critical for activation by GPCRs

One of the most critical aspects of G Protein Coupled Receptors (GPCRs) regulation is their rapid and acute desensitization following agonist stimulation. Phosphorylation of these receptors by GPCR kinases (GRK) is a major mechanism of desensitization. Considerable evidence from studies of rhodopsin kinase and GRK2 suggests there is an allosteric docking site for the receptor distinct from the GRK catalytic site. While the agonist-activated GPCR appears crucial for GRK activation, the molecular details of this interaction remain unclear. Recent studies suggested an important role for the N- and C-termini and domains in the small lobe of the kinase domain in allosteric activation; however, neither the mechanism of action of that site nor the RH domain contributions have been elucidated. To search for the allosteric site, we first indentified evolutionarily conserved sites within the RH and kinase domains presumably deterministic of protein function employing evolutionary trace (ET) methodology and crystal structures of GRK6. Focusing on a conserved cluster centered on helices 3, 9, and 10 in the RH domain, key residues of GRK5 and 6 were targeted for mutagenesis and functional assays. We found that a number of double mutations within helices 3, 9, and 10 and the N-terminus markedly reduced (50–90%) the constitutive phosphorylation of the β-2 Adrenergic Receptor (β2AR) in intact cells and phosphorylation of light-activated rhodopsin (Rho*) in vitro as compared to wild type (WT) GRK5 or 6. Based on these results, we designed peptide mimetics of GRK5 helix 9 both computationally and through chemical modifications with the goal of both confirming the importance of helix 9 and developing a useful inhibitor to disrupt the GPCR-GRK interaction. Several peptides were found to block Rho* phosphorylation by GRK5 including the native helix 9 sequence, Peptide Builder designed-peptide preserving only the key ET residues, and chemically locked helices. Most peptidomimetics showed inhibition of GRK5 activity greater than 80 % with an IC50 of ∼ 30 µM. Alanine scanning of helix 9 has further revealed both essential and non-essential residues for inhibition. Importantly, substitution of Arg 169 by an alanine in the native helix 9-based peptide gave an almost complete inhibition at 30 µM with an IC50 of ∼ 10 µM. In summary we report a previously unrecognized crucial role for the RH domain of GRK5 and 6, and the subsequent identification of a lead peptide inhibitor of protein-protein interaction with potential for specific blockade of GPCR desensitization. ^

Functional Characterization of AtxA, the Bacillus anthracis Virulence Regulator

Coordinated expression of virulence genes in Bacillus anthracis occurs via a multi-faceted signal transduction pathway that is dependent upon the AtxA protein. Intricate control of atxA gene transcription and AtxA protein function have become apparent from studies of AtxA-induced synthesis of the anthrax toxin proteins and the poly-D-glutamic acid capsule, two factors with important roles in B. anthracis pathogenesis. The amino-terminal region of the AtxA protein contains winged-helix (WH) and helix-turn-helix (HTH) motifs, structural features associated with DNA-binding. Using filter binding assays, I determined that AtxA interacted non-specifically at a low nanomolar affinity with a target promoter (Plef) and AtxA-independent promoters. AtxA also contains motifs associated with phosphoenolpyruvate: sugar phosphotransferase system (PTS) regulation. These PTS-regulated domains, PRD1 and PRD2, are within the central amino acid sequence. Specific histidines in the PRDs serve as sites of phosphorylation (H199 and H379). Phosphorylation of H199 increases AtxA activity; whereas, H379 phosphorylation decreases AtxA function. For my dissertation, I hypothesized that AtxA binds target promoters to activate transcription and that DNA-binding activity is regulated via structural changes within the PRDs and a carboxy-terminal EIIB-like motif that are induced by phosphorylation and ligand binding. I determined that AtxA has one large protease-inaccessible domain containing the PRDs and the carboxy-terminal end of the protein. These results suggest that AtxA has a domain that is distinct from the putative DNA-binding region of the protein. My data indicate that AtxA activity is associated with AtxA multimerization. Oligomeric AtxA was detected when co-affinity purification, non-denaturing gel electrophoresis, and bis(maleimido)hexane (BMH) cross-linking techniques were employed. I exploited the specificity of BMH for cysteine residues to show that AtxA was cross-linked at C402, implicating the carboxy-terminal EIIB-like region in protein-protein interactions. In addition, higher amounts of the cross-linked dimeric form of AtxA were observed when cells were cultured in conditions that promote toxin gene expression. Based on the results, I propose that AtxA multimerization requires the EIIB-like motif and multimerization of AtxA positively impacts function. I investigated the role of the PTS in the function of AtxA and the impact of phosphomimetic residues on AtxA multimerization. B. anthracis Enzyme I (EI) and HPr did not facilitate phosphorylation of AtxA in vitro. Moreover, markerless deletion of ptsHI in B. anthracis did not perturb AtxA function. Taken together, these results suggest that proteins other than the PTS phosphorylate AtxA. Point mutations mimicking phosphohistidine (H to D) and non-phosphorylated histidine (H to A) were tested for an impact on AtxA activity and multimerization. AtxA H199D, AtxA H199A, and AtxA H379A displayed multimerization phenotypes similar to that of the native protein, whereas AtxA H379D was not susceptible to BMH cross-linking or co-affinity purification with AtxA-His. These data suggest that phosphorylation of H379 may decrease AtxA activity by preventing AtxA multimerization. Overall, my data support the following model of AtxA function. AtxA binds to target gene promoters in an oligomeric state. AtxA activity is increased in response to the host-related signal bicarbonate/CO2 because this signal enhances AtxA multimerization. In contrast, AtxA activity is decreased by phosphorylation at H379 because multimerization is inhibited. Future studies will address the interplay between bicarbonate/CO2 signaling and phosphorylation on AtxA function.

Characterization of the MLC box: A conserved {\it cis\/}-acting DNA element present in the chicken MLC1s/1c gene

The slow/cardiac alkali myosin light chain (MLC1s/1c) is a member of a multigene family whose protein products are essential for activation of the myosin ATPase. In the adult, the MLC1s/1c isoform is expressed in both cardiac and slow-twitch skeletal muscles, while it is expressed by all skeletal muscles during development.^ To elucidate the molecular mechanisms that underlie the transcriptional regulation of MLC1s/1c gene expression, the immediate 5$\sp\prime$ flanking region of the gene was isolated and shown to be capable of directing reporter gene expression. Analysis of this region revealed a 110 bp muscle-specific enhancer that includes a myocyte-specific enhancer-binding factor 2 (MEF-2) site, E-boxes, which are potential binding sites for the basic-helix-loop-helix proteins such as MyoD, and a MLC box. The focus of the thesis was to identify the role of the MLC box in expression of the MLC1s/1c gene.^ The MLC box is a member of the family of CArG box containing cis-acting DNA elements. Mutagenesis showed that the MLC box is necessary, but not sufficient, for the expression of a reporter gene linked to the 5$\sp\prime$ flanking region of the MLC1s/1c gene. Linker scanner and site-directed mutagenesis identified a number of potential sites within the 110 bp muscle-specific enhancer that may cooperate with the MLC box. These are the MEF-2 site, the E-box site, and a 10 bp element located upstream of the MEF-2 site that does not have sequence similarity with any known cis-acting element. The MLC box is capable of binding to factors present in muscle nuclear extracts, as well as to human recombinant serum response factor (SRF). Binding of SRF to the MLC box was correlated with the ability of the 5$\sp\prime$ flanking region of the MLC1s/1c gene to drive reporter gene expression. Results suggest a model in which binding of SRF to the MLC box activates expression of the MLC1s/1c gene while binding of the factors present in the nuclear extracts suppresses the expression of the gene. (Abstract shortened with permission of author.) ^

Functional studies of {\it twist\/} during mouse embryogenesis

A fundamental task in developmental biology is to understand the molecular mechanisms governing early embryogenesis. The aim of this study was to understand the developmental role of a putative basic helix-loop-helix (b-HLH) transcription factor, twist, during mouse embryogenesis.^ twist was originally identified in Drosophila as one of the zygotic genes, including snail, that were required for dorsal-ventral patterning. In Drosophila embryogenesis, twist is expressed in the cells of the ventral midline destined to form mesoderm. In embryos lacking twist expression, their ventral cells fail to form a ventral furrow and subsequently no mesoderm is formed.^ During mouse embryogenesis, twist is expressed after initial mesoderm formation in both mesoderm and cranial neural crest cell derivatives. To study the role of twist in vivo, twist-null embryos were generated by gene targeting. Embryos homozygous for the twist mutation die at midgestation. The most prominent phenotype in the present study was a failure of the cranial neural tube to close (exencephaly). twist-null embryos also showed defects in head mesenchyme, branchial arches, somites, and limb buds.^ To understand whether twist functions cell-autonomously and to investigate how twist-null cells interact with wild-type cells in vivo, twist chimeras composed of both twist-null and wild-type cells marked by the expression of the lacZgene were generated. Chimeric analysis revealed a correlation between the incidence of exencephaly and the contribution of the underlying twist-null head mesenchyme, thus strongly suggesting that twist-expressing head mesenchyme is required for the closure of the cranial neural tube. These studies have identified twist as a critical regulator for the mesenchymal fate determination within the cranial neural crest lineage. Most strikingly, twist-null head mesenchyme cells were always segregated from wild-type cells, indicating that the twist mutation altered the adhesive specificity of these cells. Furthermore, these results also indicated that twist functions cell-autonomously in the head, arch, and limb mesenchyme but non-cell-autonomously in the somites. Taken together, these studies have established the essential role of twist during mouse embryogenesis. ^

Cloning and molecular analysis of paraxis: A {\it trans\/}-acting factor required for somite maturation

During vertebrate embryogenesis, cells from the paraxial mesoderm coalesce in a rostral-to-caudal progression to form the somites. Subsequent compartmentalization of the somites yields the sclerotome, myotome and dermatome, which give rise to the axial skeleton, axial musculature, and dermis, respectively. Recently, we cloned a novel basic-Helix-Loop-Helix (bHLH) protein, called scleraxis, which is expressed in the sclerotome, in mesenchymal precursors of bone and cartilage, and in connective tissues. This dissertation focuses on the cloning, expression and functional analysis of a bHLH protein termed paraxis, which is nearly identical to scleraxis within the bHLH region but diverges in both its amino and carboxyl termini. During the process of mouse embryogenesis, paraxis transcripts are first detected at about day 7.5 post coitum within the primitive mesoderm lying posterior to the head and heart primordia. Subsequently, paraxis expression progresses caudally through the paraxial mesoderm, immediately preceding somite formation. Paraxis is expressed at high levels in newly formed somites before the first detectable expression of the myogenic bHLH genes, and as the somite becomes compartmentalized, paraxis becomes downregulated within the myotome.^ To determine the function of paraxis during mammalian embryogenesis, mice were generated with a null mutation in the paraxis locus. Paraxis null mice survived until birth, but exhibited severe foreshortening along the anteroposterior axis due to the absence of vertebrae caudal to the midthoracic region. The phenotype also included axial skeletal defects, particularly shortened bifurcated ribs which were detached from the vertebral column, fused vertebrae and extensive truncation and disorganization caudal to the hindlimbs. Mutant neonates also lacked normal levels of trunk muscle and exhibited defects in the dermis as well as the stratification of the epidermis. Analysis of paraxis -/- mutant embryos has revealed a failure of the somites to both properly epithelialize and compartmentalize, resulting in defects in somite-derived cell lineages. These results suggest that paraxis is an essential component of the genetic pathway regulating somitogenesis. ^

G -rich oligodeoxyribonucleotides: Triplexes and transcription

Formation of a triple helix resulting from oligonucleotide binding to the DNA double helix offers new possibilities to control gene expression at the transcriptional level. Purine-motif triplexes can be formed under physiological pH. Nevertheless, this formation was inhibited by certain monovalent cations during the association but not during dissociation. Since triplexes are very stable, it was possible to assemble them in the absence of KCl and have them survive throughout the course of an in vitro transcription reaction. As for the design of a better triplex-forming oligonucleotide, 12 nucleotides in length afforded the highest binding affinity. G/T-rich oligonucleotides can be very polymorphic in solution. The conditions for forming purine-motif triplexes, duplexes or G-quartets were determined. Understanding these parameters will be important for the practical use of G-rich oligonucleotides in the development of DNA aptamers where the structure of the oligonucleotide is paramount in dictating its function. Finally, purine-motif triplexes were demonstrated to significantly inhibit gene transcription in vitro. The optimal effect on this process was dependent on the location of triplexes within the promoter, i.e., whether upstream or proximally downstream of the transcription start site. The mechanism for the inhibition of transcription appeared to be interference with initiation through preventing engagement by RNA polymerase. This finding is revolutionary when compared to the conventional model where triplexes inhibit transcription only by occluding binding by trans-acting proteins. Our findings broaden the utility of triplexes and support a strategy for antigene therapy by triplexes. ^

Functional analysis of GNA -1, a Galphai superfamily member found in the filamentous fungus Neurospora crassa

Heterotrimeric GTP-binding proteins, G proteins, are integral components of eukaryotic signaling systems linking extracellular signals to intracellular responses. Through coupling to seven-transmembrane helix receptors, G proteins convey primary signaling events into multi-leveled cascades of intracellular activity by regulating downstream enzymes, collectively called effectors. The effector enzymes regulated by G proteins include adenylyl cyclase, cAMP phosphodiesterase, phospolipase C-β, mitogen-activated protein kinases, and ion channels. ^ Neurospora crassa is a multicellular, filamentous fungus that is capable of both asexual and sexual reproduction by elaboration of specialized, developmentally controlled structures that give rise to either asexual or sexual spores, respectively. N. crassa possesses at least three heterotrimeric Gα proteins (GNA-1–3) and one Gβ subunit (GNB-1). GNA-1 was the first microbial protein that could be classified in the Gαi superfamily based on its amino acid identity and demonstration that it is a substrate for ADP-ribosylation by pertussis toxin. ^ Experiments were designed to identify the signal transduction pathways and the effector enzymes regulated by GNA-1. Targeted gene-replacement of gna-1 revealed that GNA-1 controls multiple developmental pathways including both asexual and sexual reproduction, maintenance of growth, and resistance to osmotic stress. The Gαi and Gαz members of the Gαi superfamily negatively regulate adenylyl cyclase activity in mammalian cells; therefore, adenylyl cyclase and cAMP levels were measured in Δgna-1 strains and also in strains that were deleted for both gna-1 and gna-2, a second Gα in N. crassa shown to have overlapping functions with GNA-1. Direct measurements of adenylyl cyclase activity revealed that GNA-1, but not GNA-2, was responsible for GTP-stimulated adenylyl cyclase activity in N. crassa. Furthermore, anti-GNA-1 IgG could specifically inhibit GTP-stimulated adenylyl cyclase activity in wild-type strain extracts. These studies also provided evidence that N. crassa possesses feedback mechanisms that control steady-state cAMP levels through indirect regulation of cAMP-phosphodiesterase activity; mutations in gna-1 and gna-2 were additive in their effect on lowering cAMP-phosphodiesterase activity under growth conditions where steady-state cAMP levels were normal but GTP-stimulated adenylyl cyclase activity was reduced 90% in comparison to control strains. ^ Genetic and biochemical epistasis experiments utilizing a Δ gna-1 cr-1 mutant suggest that GNA-1 is essential for female fertility in a cAMP-independent pathway. Furthermore, deletion of gna-1 in a cr-1 background exacerbated many of the defects already observed in the cr-1 strain including more severe growth restriction and developmental defects. However, deletion of gna-1 had no effect on the increased thermotolerance of cr-1, which has been attributed to loss of cAMP. cr-1 possesses GNA-1 protein, and crude membrane fractions from this strain reconstituted GTP-stimulated adenylyl cyclase activity in Δgna-1 membrane fractions. These studies provide direct evidence for the involvement of Gα proteins in the regulation of adenylyl cyclase activity in eukaryotic microbes. ^

Heterotrimeric G protein -mediated signal transduction in Neurospora crassa

Heterotrimeric G protein-mediated signal transduction is one of numerous means that cells utilize to respond to external stimuli. G proteins consist of α, β andγ subunits. Extracellular ligands bind to seven-transmembrane helix receptors, triggering conformational changes. This is followed by activation of coupled G proteins through the exchange of GDP for GTP on the Gα subunit. Once activated, Gα-GTP dissociates from the βγ dimer. Both of these two moieties can interact with downstream effectors, such as adenylyl cyclase, phospholipase C, phosphodiesterases, or ion channels, leading to a series of changes in cellular metabolism and physiology. ^ Neurospora crassa is a eukaryotic multicellular filamentous fungus, with asexual/vegetative and sexual phases to its life cycle. Three Gα (GNA-1, GNA-2, GNA-3) and one Gβ (GNB-1) proteins have been identified in this organism. This dissertation investigates GNA-1 and GNB-1 mediated signaling pathways in N. crassa. ^ GNA-1 was the first identified microbial Gα that belongs to a mammalian superfamily (Gαi). Deletion of GNA-1 leads to multiple defects in N. crassa. During the asexual cycle, Δgna-1 strains display a slower growth rate and delayed conidiation on solid medium. In the sexual cycle, the Δgna-1 mutant is male-fertile but female-sterile. Biochemical studies have shown that Δ gna-1 strains have lower adenosine 3′–5 ′ cyclic monophosphate (cAMP) levels than wild type under conditions where phenotypic defects are observed. In this thesis work, strains containing one of two GTPase-deficient gna-1 alleles (gna-1 R178C, gna-1Q204L) leading to constitutive activation of GNA-1 have been constructed and characterized. Activation of GNA-1 causes uncontrolled aerial hyphae proliferation, elevated sensitivity to heat and oxidative stresses, and lower carotenoid synthesis. To further study the function of GNA-1, constructs to enable expression of mammalian Gαi superfamily members were transformed into a Δ gna-1 strain, and complementation of Δgna-1 defects investigated. Gαs, which is not a member of Gα i superfamily was used as a control. These mammalian Gα genes were able to rescue the vegetative growth rate defect of the Δ gna-1 strain in the following order: Gαz > Gα o > Gαs > Gαt > Gαi. In contrast, only Gαo was able to complement the sexual defect of a Δgna-1 strain. With regard to the thermotolerance phenotype, none of the mammalian Gα genes restored the sensitivity to a wild type level. These results suggest that GNA-1 regulates two independent pathways during the vegetative and sexual cycles in N. crassa. ^ GNB-1, a G protein β subunit from N. crassa, was identified and its functions investigated in this thesis work. The sequence of the gnb-1 gene predicts a polypeptide of 358 residues with a molecular mass of 39.7 kDa. GNB-1 exhibits 91% identity to Cryphonectria parasitica CPGB-1, and also displays significant homology with human and Dictyostelium Gβ genes (∼66%). A Δ gnb-1 strain was constructed and shown to exhibit defects in asexual spore germination, vacuole number and size, mass accumulation and female fertility. A novel role for GNB-1 in regulation of GNA-1 and GNA-2 protein levels was also demonstrated. ^

Constitutively active cAMP receptor mutants dominantly inhibit receptor-mediated signaling pathways in Dictyostelium discoideum

Dictyostelium, a soil amoeba, is able to develop from free-living cells to multicellular fruiting bodies upon starvation using extracellular cAMP to mediate cell-cell communication, chemotaxis and developmental gene expression. The seven transmembrane G protein-coupled cAMP receptor-1 (cAR1) mediated responses, such as the activation of adenylyl cyclase and guanylyl cyclase, are transient, due to the existence of poorly understood adaptation mechanisms. For this dissertation, the powerful genetics of the Dictyostelium system was employed to study the adaptation mechanism of cAR1-mediated cAMP signaling as well as mechanisms intrinsic to cAR1 that regulate its activation. ^ We proposed that constitutively active cAR1 would cause constant adaptation, thus inhibiting downstream pathways that are essential for aggregation and development. Therefore, a screen for dominant negative cAR1 mutants was undertaken to identify constitutively active receptor mutants. Three dominant negative cAR1 mutants were identified. All appear to be constitutively active receptor mutants because they are constitutively phosphorylated and possess high affinity for cAMP. Biochemical studies showed that these mutant receptors prevented the activation of downstream effectors, including adenylyl and guanylyl cyclases. In addition, these cells also were defective in cAMP chemotaxis and cAR1-mediated gene expression. These findings suggest that the mutant receptors block development by constantly activating multiple adaptation pathways. ^ Sequence analysis revealed that these mutations (I104N, L100H) are clustered in a conserved region of the third transmembrane helix (TM3) of cAR1. To investigate the role of this region in receptor activation, one of these residues, I104, was mutated to all the other 19 possible amino acids. We found that all but the most conservative substitutions increase the receptor's affinity about 20- to 70-fold. However, only highly polar substitutions of I104, particularly basic residues, resulted in receptors that are constitutively phosphorylated and dominantly inhibit development, suggesting that highly polar substitutions not only disrupt an interaction constraining the receptor in its low-affinity, inactive state but also promote an additional conformational change that resembles the ligand-bound conformation. Our findings suggest that I104 plays a specific role in constraining the receptor in its inactive state and that substituting it with highly polar residues results in constitutive activation. ^

Characterization of the NOP -1 opsin photoreceptor in the filamentous fungus Neurospora crassa

Light absorption is an important process for energy production and sensory perception in many organisms. In the filamentous fungus, Neurospora crassa, blue-light is an important regulator of both asexual and sexual development, but the identity of the blue-light receptor is unknown. The work presented in this dissertation initiated the characterization of the putative N. crassa opsin photoreceptor, NOP-1. Opsins were thought to exist only in the archaea and mammals until the discovery of nop-1. All opsins have the same conserved structure of seven transmembrane helical domains with a lysine residue in the seventh helix specific for forming a Schiff-base linkage with retinal. The predicted NOP-1 protein sequence is equally similar to archaeal rhodopsins and a newly identified fungal opsin-related protein group (ORPs). ORPs maintain the seven transmembrane helical structure of opsins, but lack the conserved lysine residue for binding retinal. An ORP gene, orp-1 was identified in N. crassa and this work includes the cloning and sequence analysis of this gene. Characterization of NOP-1 function in N. crassa development began with the construction of a Δnop-1 deletion mutant. Extensive phenotypic analysis of Δnop-1 mutants revealed only subtle defects during development primarily under environmental conditions that induce a stress response. NOP-1 was overexpressed in the heterologous system Pichia pastoris, and it was demonstrated that NOP-1 protein bound all-trans retinal to form a green-light absorbing pigment (λmax = 534 nm) with a photochemical reaction cycle similar to archaeal sensory rhodopsins. nop-1 gene expression was monitored during N. crassa development. nop-1 transcript is highly expressed during asexual sporulation (conidiation) and transcript levels are abundant in the later stages of conidial development. nop-1 expression is not regulated by blue-light or elevated temperatures. Potential functions for NOP-1 were discovered through the transcriptional analysis of conidiation-associated genes in Δnop-1 mutants. NOP-1 exhibits antagonistic transcriptional regulation of conidiation-associated genes late in conidial development, by enhancing the carotenogenic gene, al-2 and repressing the conidiation-specific genes, con-10 and con-13. ^