Transcriptome analysis reveals antioxidant defense mechanisms in the red swamp crayfish Procambarus clarkia after exposure to chromium.

Chromium (Cr) as a chromate anion has a strong redox capacity that seriously threatens the ecological environment and human health. Cr can contaminate water and impart toxicity to aquatic species. Procambarus clarkii is an important food source that once represented a large proportion of the aquaculture industry due to its rapid reproduction and high economic value. However, there have been reports on the death of P. clarkii due to heavy metal pollution. The underlying mechanism regarding heavy metal toxicity was studied in this paper. The transcriptome data of hemocytes extracted from P. clarkii injected with Cr were analyzed by high-throughput sequencing and compared to the control group. In total, 48,128,748 clean reads were obtained in the treatment group and 56,480,556 clean reads were obtained in the control group. The reads were assembled using Trinity and the identified unigenes were then annotated. Then, 421 differentially-expressed genes (DEGs) were found, 170 of which were upregulated and 251 downregulated. Many of these genes were found to be related to glutathione metabolism and transportation. The glutathione metabolic pathway of P. clarkii was thus activated by Cr exposure to detoxify and maintain body function. Validation of DEGs with quantitative real-time PCR confirms the changes in gene expression. Thus, this study provides data supporting a glutathione-focused response of P. clarkii to exposure to heavy metals.

Ancient Gene Duplications, Rather Than Polyploidization, Facilitate Diversification of Petal Pigmentation Patterns in Clarkia gracilis (Onagraceae).

It has been suggested that gene duplication and polyploidization create opportunities for the evolution of novel characters. However, the connections between the effects of polyploidization and morphological novelties have rarely been examined. In this study, we investigated whether petal pigmentation patterning in an allotetraploid Clarkia gracilis has evolved as a result of polyploidization. Clarkia gracilis is thought to be derived through a recent polyploidization event with two diploid species, C. amoena huntiana and an extinct species that is closely related to C. lassenensis. We reconstructed phylogenetic relationships of the R2R3-MYBs (the regulators of petal pigmentation) from two subspecies of C. gracilis and the two purported progenitors, C. a. huntiana and C. lassenensis. The gene tree reveals that these R2R3-MYB genes have arisen through duplications that occurred before the divergence of the two progenitor species, that is, before polyploidization. After polyploidization and subsequent gene loss, only one of the two orthologous copies inherited from the progenitors was retained in the polyploid, turning it to diploid inheritance. We examined evolutionary changes in these R2R3-MYBs and in their expression, which reveals that the changes affecting patterning (including expression domain contraction, loss-of-function mutation, cis-regulatory mutation) occurred after polyploidization within the C. gracilis lineages. Our results thus suggest that polyploidization itself is not necessary in producing novel petal color patterns. By contrast, duplications of R2R3-MYB genes in the common ancestor of the two progenitors have apparently facilitated diversification of petal pigmentation patterns.

R2R3-MYB genes control petal pigmentation patterning in Clarkia gracilis ssp. sonomensis (Onagraceae).

Petal pigmentation patterning is widespread in flowering plants. The genetics of these pattern elements has been of great interest for understanding the evolution of phenotypic diversification. Here, we investigate the genetic changes responsible for the evolution of an unpigmented petal element on a colored background. We used transcriptome analysis, gene expression assays, cosegregation in F2 plants and functional tests to identify the gene(s) involved in petal coloration in Clarkia gracilis ssp. sonomensis. We identified an R2R3-MYB transcription factor (CgsMYB12) responsible for anthocyanin pigmentation of the basal region ('cup') in the petal of C. gracilis ssp. sonomensis. A functional mutation in CgsMYB12 creates a white cup on a pink petal background. Additionally, we found that two R2R3-MYB genes (CgsMYB6 and CgsMYB11) are also involved in petal background pigmentation. Each of these three R2R3-MYB genes exhibits a different spatiotemporal expression pattern. The functionality of these R2R3-MYB genes was confirmed through stable transformation of Arabidopsis. Distinct spatial patterns of R2R3-MYB expression have created the possibility that pigmentation in different sections of the petal can evolve independently. This finding suggests that recent gene duplication has been central to the evolution of petal pigmentation patterning in C. gracilis ssp. sonomensis.

QM/MM free energy simulations of salicylic acid methyltransferase: effects of stabilization of TS-like structures on substrate specificity.

Salicylic acid methyltransferases (SAMTs) synthesize methyl salicylate (MeSA) using salicylate as the substrate. MeSA synthesized in plants may function as an airborne signal to activate the expression of defense-related genes and could also be a critical mobile signaling molecule that travels from the site of plant infection to establish systemic immunity in the induction of disease resistance. Here the results of QM/MM free energy simulations for the methyl transfer process in Clarkia breweri SAMT (CbSAMT) are reported to determine the origin of the substrate specificity of SAMTs. The free energy barrier for the methyl transfer from S-adenosyl-L-methionine (AdoMet) to 4-hydroxybenzoate in CbSAMT is found to be about 5 kcal/mol higher than that from AdoMet to salicylate, consistent with the experimental observations. It is suggested that the relatively high efficiency for the methylation of salicylate compared to 4-hydroxybenzoate is due, at least in part, to the reason that a part of the stabilization of the transition state (TS) configuration is already reflected in the reactant complex, presumably, through the binding. The results seem to indicate that the creation of the substrate complex (e.g., through mutagenesis and substrate modifications) with its structure closely resembling TS might be fruitful for improving the catalytic efficiency for some enzymes. The results show that the computer simulations may provide important insights into the origin of the substrate specificity for the SABATH family and could be used to help experimental efforts in generating engineered enzymes with altered substrate specificity.

Mating system and ploidy influence levels of inbreeding depression in Clarkia (Onagraceae).

Inbreeding depression is the reduction in offspring fitness associated with inbreeding and is thought to be one of the primary forces selecting against the evolution of self-fertilization. Studies suggest that most inbreeding depression is caused by the expression of recessive deleterious alleles in homozygotes whose frequency increases as a result of self-fertilization or mating among relatives. This process leads to the selective elimination of deleterious alleles such that highly selfing species may show remarkably little inbreeding depression. Genome duplication (polyploidy) has also been hypothesized to influence levels of inbreeding depression, with polyploids expected to exhibit less inbreeding depression than diploids. We studied levels of inbreeding depression in allotetraploid and diploid species of Clarkia (Onagraceae) that vary in mating system (each cytotype was represented by an outcrossing and a selfing species). The outcrossing species exhibited more inbreeding depression than the selfing species for most fitness components and for two different measures of cumulative fitness. In contrast, though inbreeding depression was generally lower for the polyploid species than for the diploid species, the difference was statistically significant only for flower number and one of the two measures of cumulative fitness. Further, we detected no significant interaction between mating system and ploidy in determining inbreeding depression. In sum, our results suggest that a taxon's current mating system is more important than ploidy in influencing levels of inbreeding depression in natural populations of these annual plants.

Structural basis for substrate recognition in the salicylic acid carboxyl methyltransferase family.

Recently, a novel family of methyltransferases was identified in plants. Some members of this newly discovered and recently characterized methyltransferase family catalyze the formation of small-molecule methyl esters using S-adenosyl-L-Met (SAM) as a methyl donor and carboxylic acid-bearing substrates as methyl acceptors. These enzymes include SAMT (SAM:salicylic acid carboxyl methyltransferase), BAMT (SAM:benzoic acid carboxyl methyltransferase), and JMT (SAM:jasmonic acid carboxyl methyltransferase). Moreover, other members of this family of plant methyltransferases have been found to catalyze the N-methylation of caffeine precursors. The 3.0-A crystal structure of Clarkia breweri SAMT in complex with the substrate salicylic acid and the demethylated product S-adenosyl-L-homocysteine reveals a protein structure that possesses a helical active site capping domain and a unique dimerization interface. In addition, the chemical determinants responsible for the selection of salicylic acid demonstrate the structural basis for facile variations of substrate selectivity among functionally characterized plant carboxyl-directed and nitrogen-directed methyltransferases and a growing set of related proteins that have yet to be examined biochemically. Using the three-dimensional structure of SAMT as a guide, we examined the substrate specificity of SAMT by site-directed mutagenesis and activity assays against 12 carboxyl-containing small molecules. Moreover, the utility of structural information for the functional characterization of this large family of plant methyltransferases was demonstrated by the discovery of an Arabidopsis methyltransferase that is specific for the carboxyl-bearing phytohormone indole-3-acetic acid.

The 5' leader of plant PgiC has an intron: the leader shows both the loss and maintenance of constraints compared with introns and exons in the coding region.

PgiC, a complex gene with 23 coding exons and 22 intervening introns, encodes the cytosolic isozyme of phosphoglucose isomerase (EC 5.3.1.9) in higher plants. Here, we report RNA ligase-mediated rapid amplification of cDNA ends experiments that showed that PgiC in Clarkia (Onagraceae) and Arabidopsis thaliana has an intron in the 5' leader. Comparison of the EMBL accessions of the cDNA and genomic sequences showed that this is also the case in rice (Oryza sativa), suggesting that a leader intron is generally present in higher plant PgiC. The intron is bounded by consensus 5'-GT and AG-3' splice sites but showed alternative splicing in Clarkia, resulting in mature transcripts that differ by 8-19 nt in length. The intron is located 18 or 10 nt upstream of the start codon in Clarkia, 2 nt upstream in Arabidopsis, and 9 nt in rice. PgiC in Clarkia was duplicated before the divergence of the extant species, many of which have two expressed genes PgiC1 and PgiC2. Full-length transcripts of both genes identified the transcription start and made it possible to identify the leader intron and leader exon (between the transcription start and leader intron) from previously obtained genomic sequences of both genes in other Clarkia species. These data permit the comparison of evolution in the leader exon and intron with the exons and introns of the coding region, a topic that has not been studied previously. Both the leader exon and the leader intron resemble introns of the coding region in base substitution rate and accumulation of gaps. But the leader intron splice junctions are not strictly conserved in position as are those of the coding region introns. Also, in base composition, the leader intron resembles the other introns, whereas the leader exon more nearly resembles the coding exons. A difference in base composition between coding exons and flanking introns is known to be important for the recognition of splice sites. Thus, the marked difference in base composition between the leader exon and leader intron is probably maintained by selection despite a high rate of sequence divergence.

Single mutations silence PGiC2 genes in two very recent allotetraploid species of Clarkia.

Our understanding of how polyploidy influences gene evolution is limited by the fact there have been few molecular descriptions of particular genes and their expression in polyploid plants and their diploid progenitors. Here we use evidence from sequencing of genomic DNA and cDNA obtained by reverse transcriptase-polymerase chain reaction and 3' rapid amplification of cDNA ends to describe PgiC genes and their expression in two allotetraploid species of the wildflower genus Clarkia, C. delicata and C. similis. PgiC encodes the cytosolic isozyme of phosphoglucose isomerase (EC 5.3.1.9) and was duplicated in the ancestral stock of Clarkia, giving rise to paralogous genes PgiC1 and PgiC2. The active form of the PGIC enzyme is a dimer of like subunits. The electrophoretic patterns in the parent species show three bands of activity, representing two homodimers and a heterodimer of intermediate mobility, and are encoded by two genes. The electrophoretic patterns in the tetraploids also show three bands, but the tetraploids were expected to have multiple PGIC isozymes encoded by four genes. Our molecular studies demonstrated that each tetraploid has two PgiC1 and two PgiC2 genes, as predicted. One gene in each of them has been silenced by a single mutation, and a functional protein is no longer produced. In C. similis, PgiC2(mod) was silenced by a mutation of a single nucleotide in exon 5 that created a stop codon. In C. delicata, a polymorphism exists between a normal allele and a defective allele of PgiC2(epi) that has a deletion of a splice junction in intron 19 that results in the synthesis of a transcript lacking an entire exon, an example of exon skipping. The three-banded PGIC electrophoretic pattern of both tetraploid species arises because isozymes encoded by two or three of the genes comigrate. A very recent origin for both tetraploids is suggested by the near identity of several of their PgiC genes to their corresponding diploid orthologues and the absence of any acceleration in mutation rates. The problem of assessing genetic redundancy in tetraploids is discussed.