Tandem gene arrays: a challenge for functional genomics

Genomic view of heavy-ion-induced deletions associated with distribution of essential genes in Arabidopsis thaliana

Heavy-ion beam, a type of ionizing radiation, has been applied to plant breeding as a powerful mutagen and is a promising tool to induce large deletions and chromosomal rearrangements. The effectiveness of heavy-ion irradiation can be explained by linear energy transfer (LET; keV µm-1). Heavy-ion beams with different LET values induce different types and sizes of mutations. It has been suggested that deletion size increases with increasing LET value, and complex chromosomal rearrangements are induced in higher LET radiations. In this study, we mapped heavy-ion beam-induced deletions detected in Arabidopsis mutants to its genome. We revealed that deletion sizes were similar between different LETs (100 to 290 keV μm-1), that their upper limit was affected by the distribution of essential genes, and that the detected chromosomal rearrangements avoid disrupting the essential genes. We also focused on tandemly arrayed genes (TAGs), where two or more homologous genes are adjacent to one another in the genome. Our results suggested that 100 keV µm-1 of LET is enough to disrupt TAGs and that the distribution of essential genes strongly affects the heritability of mutations overlapping them. Our results provide a genomic view of large deletion inductions in the Arabidopsis genome.

Enzymes Involved in Antioxidant and Detoxification Processes Present Changes in the Expression Levels of Their Coding Genes under the Stress Caused by the Presence of Antimony in Tomato

Currently, there is an increasing presence of heavy metals and metalloids in soils and water due to anthropogenic activities. However, the biggest problem caused by this increase is the difficulty in recycling these elements and their high permanence in soils. There are plants with great capacity to assimilate these elements or make them less accessible to other organisms. We analyzed the behavior of Solanum lycopersicum L., a crop with great agronomic interest, under the stress caused by antimony (Sb). We evaluated the antioxidant response throughout different exposure times to the metalloid. Our results showed that the enzymes involved in the AsA-GSH cycle show changes in their expression level under the stress caused by Sb but could not find a relationship between the NITROSOGLUTATHIONE REDUCTASE (GSNOR) expression data and nitric oxide (NO) content in tomato roots exposed to Sb. We hypothesize that a better understanding of how these enzymes work could be key to develop more tolerant varieties to this kind of abiotic stress and could explain a greater or lesser phytoremediation capacity. Moreover, we deepened our knowledge about Glutathione S-transferase (GST) and Glutathione Reductase (GR) due to their involvement in the elimination of the xenobiotic component.

Evolution and diversification of the O-methyltransferase (OMT) gene family in Solanaceae

O-methyltransferases (OMTs) are a group of enzymes involved in several fundamental biological processes in plants, including lignin biosynthesis, pigmentation, and aroma production. Despite the intensive investigation of the role of OMTs in plant secondary metabolism, the evolution and diversification of this gene family in Solanaceae remain poorly understood. Here, we conducted a genome-wide survey of OMT genes in six Solanaceae species, reconstructing gene phylogenetic trees, predicting the potential involvement in biological processes, and investigating the exon/intron structure and chromosomal location. We identified 57 caffeoyl-CoA OMTs (CCoAOMTs) and 196 caffeic acid OMTs (COMTs) in the studied species. We observed a conserved gene block on chromosome 2 that consisted of tandem duplicated copies of OMT genes. Our results suggest that the expansion of the OMT gene family in Solanaceae was driven by whole genome duplication, segmental duplication, and tandem duplication, with multiple genes being retained by neofunctionalization and subfunctionalization. This study represents an essential first step in unraveling the evolutionary history of OMTs in Solanaceae. Our findings deepen our understanding of the crucial role of OMTs in several biological processes and highlight their significance as potential biotechnological targets.

Hidden prevalence of deletion-inversion bi-alleles in CRISPR-mediated deletions of tandemly arrayed genes in plants

Tandemly arrayed genes (TAGs) with functional redundancy and chromosomal linkage constitute 14 ~ 35% in sequenced plant genomes. The multiplex CRISPR system is the tool of choice for creating targeted TAG deletions. Here, we show that up to ~80% of CRISPR-mediated TAG knockout alleles in Arabidopsis and rice are deletion-inversion (delinver) bi-alleles, which are easily misidentified as homozygous deletion alleles by routine PCR-based genotyping. This can lead to misinterpretation of experimental data and production of progenies with genetic heterogeneity in an unnoticed manner. In ~2,650 transgenic events, delinver mutation frequencies are predominantly correlated with deletion frequencies but unrelated to chromosomal locations or deletion sizes. Delinver mutations also occur frequently at genomic non-TAG loci during multiplexed CRISPR editing. Our work raises the alarm about delinver mutations as common unwanted products of targeted TAG deletions in plants and helps prevent false interpretation of plant TAG functions due to this hidden genotype issue.

Multiplex Editing of the Nucleoredoxin1 Tandem Array in Poplar: From Small Indels to Translocations and Complex Inversions

The CRISPR-Cas9 system has been deployed for precision mutagenesis in an ever-growing number of species, including agricultural crops and forest trees. Its application to closely linked genes with extremely high sequence similarities has been less explored. In this study, we used CRISPR-Cas9 to mutagenize a tandem array of seven Nucleoredoxin1 (NRX1) genes spanning ∼100 kb in Populus tremula × Populus alba. We demonstrated efficient multiplex editing with one single guide RNA in 42 transgenic lines. The mutation profiles ranged from small insertions and deletions and local deletions in individual genes to large genomic dropouts and rearrangements spanning tandem genes. We also detected complex rearrangements including translocations and inversions resulting from multiple cleavage and repair events. Target capture sequencing was instrumental for unbiased assessments of repair outcomes to reconstruct unusual mutant alleles. The work highlights the power of CRISPR-Cas9 for multiplex editing of tandemly duplicated genes to generate diverse mutants with structural and copy number variations to aid future functional characterization.

Genome-Wide Identification and Expression Profiles of 13 Key Structural Gene Families Involved in the Biosynthesis of Rice Flavonoid Scaffolds

Flavonoids are a class of key polyphenolic secondary metabolites with broad functions in plants, including stress defense, growth, development and reproduction. Oryza sativa L. (rice) is a well-known model plant for monocots, with a wide range of flavonoids, but the key flavonoid biosynthesis-related genes and their molecular features in rice have not been comprehensively and systematically characterized. Here, we identified 85 key structural gene candidates associated with flavonoid biosynthesis in the rice genome. They belong to 13 families potentially encoding chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonol synthase (FLS), leucoanthocyanidin dioxygenase (LDOX), anthocyanidin synthase (ANS), flavone synthase II (FNSII), flavanone 2-hydroxylase (F2H), flavonoid 3′-hydroxylase (F3′H), flavonoid 3′,5′-hydroxylase (F3′5′H), dihydroflavonol 4-reductase (DFR), anthocyanidin reductase (ANR) and leucoanthocyanidin reductase (LAR). Through structural features, motif analyses and phylogenetic relationships, these gene families were further grouped into five distinct lineages and were examined for conservation and divergence. Subsequently, 22 duplication events were identified out of a total of 85 genes, among which seven pairs were derived from segmental duplication events and 15 pairs were from tandem duplications, demonstrating that segmental and tandem duplication events play important roles in the expansion of key flavonoid biosynthesis-related genes in rice. Furthermore, these 85 genes showed spatial and temporal regulation in a tissue-specific manner and differentially responded to abiotic stress (including six hormones and cold and salt treatments). RNA-Seq, microarray analysis and qRT-PCR indicated that these genes might be involved in abiotic stress response, plant growth and development. Our results provide a valuable basis for further functional analysis of the genes involved in the flavonoid biosynthesis pathway in rice.

Segmental and tandem chromosome duplications led to divergent evolution of the chalcone synthase gene family in Phalaenopsis orchids

BACKGROUND AND AIMS

Orchidaceae is a large plant family, and its extraordinary adaptations may have guaranteed its evolutionary success. Flavonoids are a group of secondary metabolites that mediate plant acclimation to challenge environments. Chalcone synthase (CHS) catalyses the initial step in the flavonoid biosynthetic pathway. This is the first chromosome-level investigation of the CHS gene family in Phalaenopsis aphrodite and was conducted to elucidate if divergence of this gene family is associated with chromosome evolution.

METHODS

Complete CHS genes were identified from our whole-genome sequencing data sets and their gene expression profiles were obtained from our transcriptomic data sets. Fluorescence in situ hybridization (FISH) was conducted to position five CHS genes to high-resolution pachytene chromosomes.

KEY RESULTS

The five Phalaenopsis CHS genes can be classified into three groups, PaCHS1, PaCHS2 and the tandemly arrayed three-gene cluster, which diverged earlier than those of the orchid genera and species. Additionally, pachytene chromosome-based FISH mapping showed that the three groups of CHS genes are localized on three distinct chromosomes. Moreover, an expression analysis of RNA sequencing revealed that the five CHS genes had highly differentiated expression patterns and its expression pattern-based clustering showed high correlations between sequence divergences and chromosomal localizations of the CHS gene family in P. aphrodite.

CONCLUSIONS

Based on their phylogenetic relationships, expression clustering analysis and chromosomal distributions of the five paralogous PaCHS genes, we proposed that expansion of this gene family in P. aphrodite occurred through segmental duplications, followed by tandem duplications. These findings provide information for further studies of CHS functions and regulations, and shed light on the divergence of an important gene family in orchids.

A differentially regulated AP2/ERF transcription factor gene cluster acts downstream of a MAP kinase cascade to modulate terpenoid indole alkaloid biosynthesis in Catharanthus roseus

Catharanthus roseus produces bioactive terpenoid indole alkaloids (TIAs), including the chemotherapeutics, vincristine and vinblastine. Transcriptional regulation of TIA biosynthesis is not fully understood. The jasmonic acid (JA)-responsive AP2/ERF transcription factor (TF), ORCA3, and its regulator, CrMYC2, play key roles in TIA biosynthesis. ORCA3 forms a physical cluster with two uncharacterized AP2/ERFs, ORCA4 and 5. Here, we report that (1) the ORCA gene cluster is differentially regulated; (2) ORCA4, while overlapping functionally with ORCA3, modulates an additional set of TIA genes. Unlike ORCA3, ORCA4 overexpression resulted in dramatic increase of TIA accumulation in C. roseus hairy roots. In addition, CrMYC2 is capable of activating ORCA3 and co-regulating TIA pathway genes concomitantly with ORCA3. The ORCA gene cluster and CrMYC2 act downstream of a MAP kinase cascade that includes a previously uncharacterized MAP kinase kinase, CrMAPKK1. Overexpression of CrMAPKK1 in C. roseus hairy roots upregulated TIA pathways genes and increased TIA accumulation. This work provides detailed characterization of a TF gene cluster and advances our understanding of the transcriptional and post-translational regulatory mechanisms that govern TIA biosynthesis in C. roseus.

Targeted deletion and inversion of tandemly arrayed genes in Arabidopsis thaliana using zinc finger nucleases

Tandemly arrayed genes (TAGs) or gene clusters are prevalent in higher eukaryotic genomes. For example, approximately 17% of genes are organized in tandem in the model plant Arabidopsis thaliana. The genetic redundancy created by TAGs presents a challenge for reverse genetics. As molecular scissors, engineered zinc finger nucleases (ZFNs) make DNA double-strand breaks in a sequence-specific manner. ZFNs thus provide a means to delete TAGs by creating two double-strand breaks in the gene cluster. Using engineered ZFNs, we successfully targeted seven genes from three TAGs on two Arabidopsis chromosomes, including the well-known RPP4 gene cluster, which contains eight resistance (R) genes. The resulting gene cluster deletions ranged from a few kb to 55 kb with frequencies approximating 1% in somatic cells. We also obtained large chromosomal deletions of ~9 Mb at approximately one tenth the frequency, and gene cluster inversions and duplications also were achieved. This study demonstrates the ability to use sequence-specific nucleases in plants to make targeted chromosome rearrangements and create novel chimeric genes for reverse genetics and biotechnology.

Plant genome engineering with sequence-specific nucleases

Recent advances in genome engineering provide newfound control over a plant's genetic material. It is now possible for most bench scientists to alter DNA in living plant cells in a variety of ways, including introducing specific nucleotide substitutions in a gene that change a protein's amino acid sequence, deleting genes or chromosomal segments, and inserting foreign DNA at precise genomic locations. Such targeted DNA sequence modifications are enabled by sequence-specific nucleases that create double-strand breaks in the genomic loci to be altered. The repair of the breaks, through either homologous recombination or nonhomologous end joining, can be controlled to achieve the desired sequence modification. Genome engineering promises to advance basic plant research by linking DNA sequences to biological function. Further, genome engineering will enable plants' biosynthetic capacity to be harnessed to produce the many agricultural products required by an expanding world population.

Co-expression of soybean Dicer-like genes in response to stress and development

Regulation of gene transcription and post-transcriptional processes is critical for proper development, genome integrity, and stress responses in plants. Many genes involved in the key processes of transcriptional and post-transcriptional regulation have been well studied in model diploid organisms. However, gene and genome duplication may alter the function of the genes involved in these processes. To address this question, we assayed the stress-induced transcription patterns of duplicated gene pairs involved in RNAi and DNA methylation processes in the paleopolyploid soybean. Real-time quantitative PCR and Sequenom MassARRAY expression assays were used to profile the relative expression ratios of eight gene pairs across eight different biotic and abiotic stress conditions. The transcriptional responses to stress for genes involved in DNA methylation, RNAi processing, and miRNA processing were compared. The strongest evidence for pairwise co-expression in response to stresses was exhibited by non-paralogous Dicer-like (DCL) genes GmDCL2a-GmDCL3a and GmDCL1b-GmDCL2b, most profoundly in root tissues. Among homoeologous or paralogous DCL genes, the Dicer-like 2 (DCL2) gene pair exhibited the strongest response to stress and most conserved co-expression pattern. This was surprising because the DCL2 duplication event is more ancient than the other DCL duplications. Possible mechanisms that may be driving the DCL2 co-expression are discussed.

Using transgenic modulation of protein synthesis and accumulation to probe protein signaling networks in Arabidopsis thaliana

Deployment of new model species in the plant biology community requires the development and/or improvement of numerous genetic tools. Sequencing of the Arabidopsis thaliana genome opened up a new challenge of assigning biological function to each gene. As many genes exhibit spatiotemporal or other conditional regulation of biological processes, probing for gene function necessitates applications that can be geared toward temporal, spatial and quantitative functional analysis in vivo. The continuing quest to establish new platforms to examine plant gene function has resulted in the availability of numerous genomic and proteomic tools. Classical and more recent genome-wide experimental approaches include conventional mutagenesis, tagged DNA insertional mutagenesis, ectopic expression of transgenes, activation tagging, RNA interference and two-component transactivation systems. The utilization of these molecular tools has resulted in conclusive evidence for the existence of many genes, and expanded knowledge on gene structure and function. This review covers several molecular tools that have become increasingly useful in basic plant research. We discuss their advantages and limitations for probing cellular protein function while emphasizing the contributions made to lay the fundamental groundwork for genetic manipulation of crops using plant biotechnology.

Origin and evolutionary analysis of the plant-specific TIFY transcription factor family

A substantial number of transcription factor families have been investigated from all kingdoms of life, but a particular class of plant-specific TIFY transcription factors, characterized by a highly conserved TIFY domain, lacks a systemic analysis of its origin and evolutionary relationships among different plant species. After exhaustive genome-wide searches against 14 genomes, TIFY transcription factors were identified and classified into four subfamilies TIFY, PPD, JAZ and ZML according to their different domain architectures. Results show that the TIFY domain of the ZML subfamily possesses a core "TLS[F/Y]XG" motif rather than the "TIFYXG" motif that is dominant in the other three subfamilies. A comprehensive survey of the TIFY family allowed us to discover a new group within the JAZ subfamily and to identify several novel conserved motifs via phylogenetic analysis. Evolutional analysis indicates that whole genome duplication and tandem duplication contributed to the expansion of the TIFY family in plants.

Targeted mutagenesis in Arabidopsis using zinc-finger nucleases

We report here an efficient method for making targeted mutations in Arabidopsis thaliana genes. The approach uses zinc-finger nucleases (ZFNs) - enzymes engineered to create DNA double-strand breaks at specific target loci. Imprecise repair of double-strand breaks by nonhomologous end-joining generates small insertions or deletions at the cleavage site. In this protocol, constructs encoding ZFNs for specific loci are transformed into Arabidopsis by Agrobacterium-mediated transformation. ZFN expression is induced during germination to initiate mutagenesis of the target locus. Typically, more than 20% of the primary transgenics segregate loss-of-function mutations in the next generation. ZFNs make it possible to expand the range of Arabidopsis mutants available for study and to create mutations in genes missed by random mutagenesis approaches, such as those using T-DNA, transposons, or chemical mutagens.

Clustered transcription factor genes regulate nicotine biosynthesis in tobacco

Tobacco (Nicotiana tabacum) synthesizes nicotine and related pyridine alkaloids in the root, and their synthesis increases upon herbivory on the leaf via a jasmonate-mediated signaling cascade. Regulatory NIC loci that positively regulate nicotine biosynthesis have been genetically identified, and their mutant alleles have been used to breed low-nicotine tobacco varieties. Here, we report that the NIC2 locus, originally called locus B, comprises clustered transcription factor genes of an ethylene response factor (ERF) subfamily; in the nic2 mutant, at least seven ERF genes are deleted altogether. Overexpression, suppression, and dominant repression experiments using transgenic tobacco roots showed both functional redundancy and divergence among the NIC2-locus ERF genes. These transcription factors recognized a GCC-box element in the promoter of a nicotine pathway gene and specifically activated all known structural genes in the pathway. The NIC2-locus ERF genes are expressed in the root and upregulated by jasmonate with kinetics that are distinct among the members. Thus, gene duplication events generated a cluster of highly homologous transcription factor genes with transcriptional and functional diversity. The NIC2-locus ERFs are close homologs of ORCA3, a jasmonate-responsive transcriptional activator of indole alkaloid biosynthesis in Catharanthus roseus, indicating that the NIC2/ORCA3 ERF subfamily was recruited independently to regulate jasmonate-inducible secondary metabolism in distinct plant lineages.

High frequency targeted mutagenesis in Arabidopsis thaliana using zinc finger nucleases

We report here an efficient method for targeted mutagenesis of Arabidopsis genes through regulated expression of zinc finger nucleases (ZFNs)-enzymes engineered to create DNA double-strand breaks at specific target loci. ZFNs recognizing the Arabidopsis ADH1 and TT4 genes were made by Oligomerized Pool ENgineering (OPEN)-a publicly available, selection-based platform that yields high quality zinc finger arrays. The ADH1 and TT4 ZFNs were placed under control of an estrogen-inducible promoter and introduced into Arabidopsis plants by floral-dip transformation. Primary transgenic Arabidopsis seedlings induced to express the ADH1 or TT4 ZFNs exhibited somatic mutation frequencies of 7% or 16%, respectively. The induced mutations were typically insertions or deletions (1-142 bp) that were localized at the ZFN cleavage site and likely derived from imprecise repair of chromosome breaks by nonhomologous end-joining. Mutations were transmitted to the next generation for 69% of primary transgenics expressing the ADH1 ZFNs and 33% of transgenics expressing the TT4 ZFNs. Furthermore, approximately 20% of the mutant-producing plants were homozygous for mutations at ADH1 or TT4, indicating that both alleles were disrupted. ADH1 and TT4 were chosen as targets for this study because of their selectable or screenable phenotypes (adh1, allyl alcohol resistance; tt4, lack of anthocyanins in the seed coat). However, the high frequency of observed ZFN-induced mutagenesis suggests that targeted mutations can readily be recovered by simply screening progeny of primary transgenic plants by PCR and DNA sequencing. Taken together, our results suggest that it should now be possible to obtain mutations in any Arabidopsis target gene regardless of its mutant phenotype.

Deletion-based reverse genetics in Medicago truncatula

The primary goal of reverse genetics, the identification of null mutations in targeted genes, is achieved through screening large populations of randomly mutagenized plants. T-DNA and transposon-based mutagenesis has been widely employed but is limited to species in which transformation and tissue culture are efficient. In other species, TILLING (for Targeting Induced Local Lesions IN Genomes), based on chemical mutagenesis, has provided an efficient method for the identification of single base pair mutations, only 5% of which will be null mutations. Furthermore, the efficiency of inducing point mutations, like insertion-based mutations, is dependent on target size. Here, we describe an alternative reverse genetic strategy based on physically induced genomic deletions that, independent of target size, exclusively recovers knockout mutants. Deletion TILLING (De-TILLING) employs fast neutron mutagenesis and a sensitive polymerase chain reaction-based detection. A population of 156,000 Medicago truncatula plants has been structured as 13 towers each representing 12,000 M2 plants. The De-TILLING strategy allows a single tower to be screened using just four polymerase chain reaction reactions. Dual screening and three-dimensional pooling allows efficient location of mutants from within the towers. With this method, we have demonstrated the detection of mutants from this population at a rate of 29% using five targets per gene. This De-TILLING reverse genetic strategy is independent of tissue culture and efficient plant transformation and therefore applicable to any plant species. De-TILLING mutants offer advantages for crop improvement as they possess relatively few background mutations and no exogenous DNA.

An Ac/Ds-mediated gene trap system for functional genomics in barley

Background

Gene trapping is a powerful tool for gene discovery and functional genomics in both animals and plants. Upon insertion of the gene trap construct into an expressed gene, splice donor and acceptor sites facilitate the generation of transcriptional fusions between the flanking sequence and the reporter. Consequently, detection of reporter gene expression allows the identification of genes based on their expression pattern. Up to now rice is the only cereal crop for which gene trap approaches exist. In this study we describe a gene trap system in barley (Hordeum vulgare L.) based on the maize transposable elements Ac/Ds.

Results

We generated gene trap barley lines by crossing Ac transposase expressing plants with multiple independent transformants carrying the Ds based gene trap construct GTDsB. Upstream of the β-Glucuronidase start codon GTDsB carries splice donor and acceptor sites optimized for monocotyledonous plants. DNA blot analysis revealed GTDsB transposition frequencies of 11% and 26% in the F₁ and F₂ generation of gene trap lines and perpetuation of transposition activity in later generations. Furthermore, analysis of sequences flanking transposed GTDsB elements evidenced preferential insertion into expressed regions of the barley genome. We screened leaves, nodes, immature florets, pollinated florets, immature grains and seedlings of F₂ plants and detected GUS expression in 51% (72/141) of the plants. Thus, reporter gene expression was found in 24 of the 28 F₁ lines tested and in progeny of all GTDsB parental lines.

Conclusion

Due to the frequent transposition of GTDsB and the efficient expression of the GUS reporter gene, we conclude that this Ac/Ds-based gene trap system is an applicable approach for gene discovery in barley. The successful introduction of a gene trap construct optimized for monocots in barley contributes a novel functional genomics tool for this cereal crop.

Dormancy-associated MADS genes from the EVG locus of peach [Prunus persica (L.) Batsch] have distinct seasonal and photoperiodic expression patterns

Mapping and sequencing of the non-dormant evg mutant in peach [Prunus persica (L.) Batsch] identified six tandem-arrayed DAM (dormancy-associated MADS-box) genes as candidates for regulating growth cessation and terminal bud formation. To narrow the list of candidate genes, an attempt was made to associate bud phenology with the seasonal and environmental patterns of expression of the candidates in wild-type trees. The expression of the six peach DAM genes at the EVG locus of peach was characterized throughout an annual growing cycle in the field, and under controlled conditions in response to a long day-short day photoperiod transition. DAM1, 2, 4, 5, and 6 were responsive to a reduction in photoperiod in controlled conditions and the direction of response correlated with the seasonal timing of expression in field-grown trees. DAM3 did not respond to photoperiod and may be regulated by chilling temperatures. The DAM genes in peach appear to have at least four distinct patterns of expression. DAM1, 2, and 4 are temporally associated with seasonal elongation cessation and bud formation and are the most likely candidates for control of the evg phenotype.

Functional analysis of a predicted flavonol synthase gene family in Arabidopsis

The genome of Arabidopsis (Arabidopsis thaliana) contains five sequences with high similarity to FLAVONOL SYNTHASE1 (AtFLS1), a previously characterized flavonol synthase gene that plays a central role in flavonoid metabolism. This apparent redundancy suggests the possibility that Arabidopsis uses multiple isoforms of FLS with different substrate specificities to mediate the production of the flavonols, quercetin and kaempferol, in a tissue-specific and inducible manner. However, biochemical and genetic analysis of the six AtFLS sequences indicates that, although several of the members are expressed, only AtFLS1 encodes a catalytically competent protein. AtFLS1 also appears to be the only member of this group that influences flavonoid levels and the root gravitropic response in seedlings under nonstressed conditions. This study showed that the other expressed AtFLS sequences have tissue- and cell type-specific promoter activities that overlap with those of AtFLS1 and encode proteins that interact with other flavonoid enzymes in yeast two-hybrid assays. Thus, it is possible that these "pseudogenes" have alternative, noncatalytic functions that have not yet been uncovered.