Evolution of cyclohexadienyl dehydratase from an ancestral solute-binding protein.

The emergence of enzymes through the neofunctionalization of noncatalytic proteins is ultimately responsible for the extraordinary range of biological catalysts observed in nature. Although the evolution of some enzymes from binding proteins can be inferred by homology, we have a limited understanding of the nature of the biochemical and biophysical adaptations along these evolutionary trajectories and the sequence in which they occurred. Here we reconstructed and characterized evolutionary intermediate states linking an ancestral solute-binding protein to the extant enzyme cyclohexadienyl dehydratase. We show how the intrinsic reactivity of a desolvated general acid was harnessed by a series of mutations radiating from the active site, which optimized enzyme-substrate complementarity and transition-state stabilization and minimized sampling of noncatalytic conformations. Our work reveals the molecular evolutionary processes that underlie the emergence of enzymes de novo, which are notably mirrored by recent examples of computational enzyme design and directed evolution.

Identification of a small protein domain present in all plant lineages that confers high prephenate dehydratase activity.

l-Phenylalanine serves as a building block for the biosynthesis of proteins, but also as a precursor for a wide range of plant-derived compounds essential for plants and animals. Plants can synthesize Phe within the plastids using arogenate as a precursor; however, an alternative pathway using phenylpyruvate as an intermediate, described for most microorganisms, has recently been proposed. The functionality of this pathway requires the existence of enzymes with prephenate dehydratase (PDT) activity (EC 4.2.1.51) in plants. Using phylogenetic studies, functional complementation assays in yeast and biochemical analysis, we have identified the enzymes displaying PDT activity in Pinus pinaster. Through sequence alignment comparisons and site-directed mutagenesis we have identified a 22-amino acid region conferring PDT activity (PAC domain) and a single Ala314 residue critical to trigger this activity. Our results demonstrate that all plant clades include PAC domain-containing ADTs, suggesting that the PDT activity, and thus the ability to synthesize Phe using phenylpyruvate as an intermediate, has been preserved throughout the evolution of plants. Moreover, this pathway together with the arogenate pathway gives plants a broad and versatile capacity to synthesize Phe and its derived compounds. PAC domain-containing enzymes are also present in green and red algae, and glaucophytes, the three emerging clades following the primary endosymbiont event resulting in the acquisition of plastids in eukaryotes. The evolutionary prokaryotic origin of this domain is discussed.

Conserved Arginines at the P-Protein Stalk Binding Site and the Active Site Are Critical for Ribosome Interactions of Shiga Toxins but Do Not Contribute to Differences in the Affinity of the A1 Subunits for the Ribosome.

The A1 subunits of Shiga toxin 1 (Stx1A1) and Shiga toxin 2 (Stx2A1) interact with the conserved C termini of ribosomal-stalk P-proteins to remove a specific adenine from the sarcin/ricin loop. We previously showed that Stx2A1 has higher affinity for the ribosome and higher catalytic activity than Stx1A1. To determine if conserved arginines at the distal face of the active site contribute to the higher affinity of Stx2A1 for the ribosome, we mutated Arg172, Arg176, and Arg179 in both toxins. We show that Arg172 and Arg176 are more important than Arg179 for the depurination activity and toxicity of Stx1A1 and Stx2A1. Mutation of a single arginine reduced the depurination activity of Stx1A1 more than that of Stx2A1. In contrast, mutation of at least two arginines was necessary to reduce depurination by Stx2A1 to a level similar to that of Stx1A1. R176A and R172A/R176A mutations eliminated interaction of Stx1A1 and Stx2A1 with ribosomes and with the stalk, while mutation of Arg170 at the active site reduced the binding affinity of Stx1A1 and Stx2A1 for the ribosome, but not for the stalk. These results demonstrate that conserved arginines at the distal face of the active site are critical for interactions of Stx1A1 and Stx2A1 with the stalk, while a conserved arginine at the active site is critical for non-stalk-specific interactions with the ribosome. Arginine mutations at either site reduced ribosome interactions of Stx1A1 and Stx2A1 similarly, indicating that conserved arginines are critical for ribosome interactions but do not contribute to the higher affinity of Stx2A1 for the ribosome.

Characterization of two key enzymes for aromatic amino acid biosynthesis in symbiotic archaea.

Biosynthesis of L-tyrosine (L-Tyr) and L-phenylalanine (L-Phe) is directed by the interplay of three enzymes. Chorismate mutase (CM) catalyzes the rearrangement of chorismate to prephenate, which can be either converted to hydroxyphenylpyruvate by prephenate dehydrogenase (PD) or to phenylpyruvate by prephenate dehydratase (PDT). This work reports the first characterization of a trifunctional PD-CM-PDT from the smallest hyperthermophilic archaeon Nanoarchaeum equitans and a bifunctional CM-PD from its host, the crenarchaeon Ignicoccus hospitalis. Hexa-histidine tagged proteins were expressed in Escherichia coli and purified by affinity chromatography. Specific activities determined for the trifunctional enzyme were 21, 80, and 30 U/mg for CM, PD, and PDT, respectively, and 47 and 21 U/mg for bifunctional CM and PD, respectively. Unlike most PDs, these two archaeal enzymes were insensitive to regulation by L-Tyr and preferred NADP(+) to NAD(+) as a cofactor. Both the enzymes were highly thermally stable and exhibited maximal activity at 90 °C. N. equitans PDT was feedback inhibited by L-Phe (Ki = 0.8 µM) in a non-competitive fashion consistent with L-Phe's combination at a site separate from that of prephenate. Our results suggest that PD from the unique symbiotic archaeal pair encompass a distinct subfamily of prephenate dehydrogenases with regard to their regulation and co-substrate specificity.

Enhanced production of L-phenylalanine in Corynebacterium glutamicum due to the introduction of Escherichia coli wild-type gene aroH.

Metabolic engineering is a powerful tool which has been widely used for producing valuable products. For improving L-phenylalanine (L-Phe) accumulation in Corynebacterium glutamicum, we have investigated the target genes involved in the biosynthetic pathways. The genes involved in the biosynthesis of L-Phe were found to be strictly regulated genes by feedback inhibition. As a result, overexpression of the native wild-type genes aroF, aroG or pheA resulted in a slight increase of L-Phe. In contrast, overexpression of aroF (wt) or pheA (fbr) from E. coli significantly increased L-Phe production. Co-overexpression of aroF (wt) and pheA (fbr) improved the titer of L-Phe to 4.46 ± 0.06 g l⁻¹. To further analyze the target enzymes in the aromatic amino acid synthesis pathway between C. glutamicum and E. coli, the wild-type gene aroH from E. coli was overexpressed and evaluated in C. glutamicum. As predicted, upregulation of the wild-type gene aroH resulted in a remarkable increase of L-Phe production. Co-overexpression of the mutated pheA (fbr) and the wild-type gene aroH resulted in the production of L-Phe up to 4.64 ± 0.09 g l⁻¹. Based on these results we conclude that the wild-type gene aroH from E. coli is an appropriate target gene for pathway engineering in C. glutamicum for the production of aromatic amino acids.

Identification of a pheA gene associated with Streptococcus mitis by using suppression subtractive hybridization.

We performed suppression subtractive hybridization to identify genomic differences between Streptococcus mitis and Streptococcus pneumoniae. Based on the pheA gene, a primer set specific to S. mitis detection was found in 18 out of 103 S. mitis-specific clones. Our findings would be useful for discrimination of S. mitis from other closely related cocci in the oral environment.

The role of oligomerization in the optimization of cyclohexadienyl dehydratase conformational dynamics and catalytic activity.

The emergence of oligomers is common during the evolution and diversification of protein families, yet the selective advantage of oligomerization is often cryptic or unclear. Oligomerization can involve the formation of isologous head-to-head interfaces (e.g., in symmetrical dimers) or heterologous head-to-tail interfaces (e.g., in cyclic complexes), the latter of which is less well studied and understood. In this work, we retrace the emergence of the trimeric form of cyclohexadienyl dehydratase from Pseudomonas aeruginosa (PaCDT) by introducing residues that form the PaCDT trimer-interfaces into AncCDT-5 (a monomeric reconstructed ancestor of PaCDT). We find that single interface mutations can switch the oligomeric state of the variants and that trimerization corresponds with a reduction in the KM value of the enzyme from a promiscuous level to the physiologically relevant range. In addition, we find that removal of a C-terminal extension present in PaCDT leads to a variant with reduced catalytic activity, indicating that the C-terminal region has a role in tuning enzymatic activity. We show that these observations can be rationalized at the structural and dynamic levels, with trimerization and C-terminal extension leading to reduced sampling of non-catalytic conformational substates in molecular dynamics simulations. Overall, this work provides insight into how neutral sampling of distinct oligomeric states along an evolutionary trajectory can facilitate the evolution and optimization of enzyme function.

13C isotope effect on the reaction catalyzed by prephenate dehydratase.

The (13)C isotope effect for the conversion of prephenate to phenylpyruvate by the enzyme prephenate dehydratase from Methanocaldococcus jannaschii is 1.0334+/-0.0006. The size of this isotope effect suggests that the reaction is concerted. From the X-ray structure of a related enzyme, it appears that the only residue capable of acting as the general acid needed for removal of the hydroxyl group is threonine-172, which is contained in a conserved TRF motif. The more favorable entropy of activation for the enzyme-catalyzed process (25 eu larger than for the acid-catalyzed reaction) has been explained by a preorganized microenvironment that obviates the need for extensive solvent reorganization. This is consistent with forced planarity of the ring and side chain, which would place the leaving carboxyl and hydroxyl out of plane. Such distortion of the substrate may be a major contributor to catalysis.

Development of L-tryptophan production strains by defined genetic modification in Escherichia coli.

Construction and improvement of industrial strains play a central role in the commercial development of microbial fermentation processes. L-tryptophan producers have usually been developed by classical random mutagenesis due to its complicated metabolic network and regulatory mechanism. However, in the present study, an L-tryptophan overproducing Escherichia coli strain was developed by defined genetic modification methodology. Feedback inhibitions of 3-deoxy-D-arabinoheptulosonate 7-phosphate synthase (AroF) and anthranilate synthase (TrpED) were eliminated by site-directed mutagenesis. Expression of deregulated AroF and TrpED was achieved by using a temperature-inducible expression plasmid pSV. Transcriptional regulation of trp repressor was removed by deleting trpR. Pathway for L-Trp degradation was removed by deleting tnaA. L-phenylalanine and L-tyrosine biosynthesis pathways that compete with L-tryptophan biosynthesis were blocked by deleting their critical genes (pheA and tyrA). The final engineered E. coli can produce 13.3 g/l of L-tryptophan. Fermentation characteristics of the engineered strains were also analyzed.

Metabolic engineering of Escherichia coli for improving L-3,4-dihydroxyphenylalanine (L-DOPA) synthesis from glucose.

L-3,4-dihydroxyphenylalanine (L-DOPA) is an aromatic compound employed for the treatment of Parkinson's disease. Metabolic engineering was applied to generate Escherichia coli strains for the production of L-DOPA from glucose by modifying the phosphoenolpyruvate:sugar phosphotransferase system (PTS) and aromatic biosynthetic pathways. Carbon flow was directed to the biosynthesis of L-tyrosine (L-Tyr), an L-DOPA precursor, by transforming strains with compatible plasmids carrying genes encoding a feedback-inhibition resistant version of 3-deoxy-D-arabino-heptulosonate-7-phosphate synthase, transketolase, the chorismate mutase domain from chorismate mutase-prephenate dehydratase from E. coli and cyclohexadienyl dehydrogenase from Zymomonas mobilis. The effects on L-Tyr production of PTS inactivation (PTS(-) gluc(+) phenotype), as well as inactivation of the regulatory protein TyrR, were evaluated. PTS inactivation caused a threefold increase in the specific rate of L-Tyr production (q( L-Tyr)), whereas inactivation of TyrR caused 1.7- and 1.9-fold increases in q( L-Tyr) in the PTS(+) and the PTS(-) gluc(+) strains, respectively. An 8.6-fold increase in L-Tyr yield from glucose was observed in the PTS(-) gluc(+) tyrR (-) strain. Expression of hpaBC genes encoding the enzyme 4-hydroxyphenylacetate 3-hydroxylase from E. coli W in the strains modified for L-Tyr production caused the synthesis of L-DOPA. One of such strains, having the PTS(-) gluc(+) tyrR (-) phenotype, displayed the best production parameters in minimal medium, with a specific rate of L-DOPA production of 13.6 mg/g/h, L-DOPA yield from glucose of 51.7 mg/g and a final L-DOPA titer of 320 mg/l. In a batch fermentor culture in rich medium this strain produced 1.51 g/l of L-DOPA in 50 h.

The Komagataeibacter europaeus GqqA is the prototype of a novel bifunctional N-Acyl-homoserine lactone acylase with prephenate dehydratase activity.

Previously, we reported the isolation of a quorum quenching protein (QQ), designated GqqA, from Komagataeibacter europaeus CECT 8546 that is highly homologous to prephenate dehydratases (PDT) (Valera et al. in Microb Cell Fact 15, 88. https://doi.org/10.1186/s12934-016-0482-y , 2016). GqqA strongly interfered with N-acyl-homoserine lactone (AHL) quorum sensing signals from Gram-negative bacteria and affected biofilm formation in its native host strain Komagataeibacter europaeus. Here we present and discuss data identifying GqqA as a novel acylase. ESI-MS-MS data showed unambiguously that GqqA hydrolyzes the amide bond of the acyl side-chain of AHL molecules, but not the lactone ring. Consistent with this observation the protein sequence does not carry a conserved Zn2+ binding motif, known to be essential for metal-dependent lactonases, but in fact harboring the typical periplasmatic binding protein domain (PBP domain), acting as catalytic domain. We report structural details for the native structure at 2.5 Å resolution and for a truncated GqqA structure at 1.7 Å. The structures obtained highlight that GqqA acts as a dimer and complementary docking studies indicate that the lactone ring of the substrate binds within a cleft of the PBP domain and interacts with polar residues Y16, S17 and T174. The biochemical and phylogenetic analyses imply that GqqA represents the first member of a novel type of QQ family enzymes.

Expression of a bacterial bi-functional chorismate mutase/prephenate dehydratase modulates primary and secondary metabolism associated with aromatic amino acids in Arabidopsis.

Plants can synthesize the aromatic amino acid Phe via arogenate, but it is still not known whether they also use an alternative route for Phe biosynthesis via phenylpyruvate, like many micro-organisms. To examine this possibility, we expressed a bacterial bi-functional PheA (chorismate mutase/prephenate dehydratase) gene in Arabidopsis thaliana that converts chorismate via prephenate into phenylpyruvate. The PheA-expressing plants showed a large increase in the level of Phe, implying that they can convert phenylpyruvate into Phe. In addition, PheA expression rendered the plants more sensitive than wild-type plants to the Trp biosynthesis inhibitor 5-methyl-Trp, implying that Phe biosynthesis competes with Trp biosynthesis from their common precursor chorismate. Surprisingly, GC-MS, LC-MS and microarray analyses showed that this increase in Phe accumulation only had a very minor effect on the levels of other primary metabolites as well as on the transcriptome profile, implying little regulatory cross-interaction between the aromatic amino acid biosynthesis network and the bulk of the Arabidopsis transcriptome and primary metabolism. However, the levels of a number of secondary metabolites derived from all three aromatic amino acids (Phe, Trp and Tyr) were altered in the PheA plants, implying regulatory cross-interactions between the flux of aromatic amino acid biosynthesis from chorismate and their further metabolism into various secondary metabolites. Taken together, our results provide insights into the regulatory mechanisms of aromatic amino acid biosynthesis and their interaction with central primary metabolism, as well as the regulatory interface between primary and secondary metabolism.

Pathway Engineering for Phenethylamine Production in Escherichia coli.

In this study, the metabolic pathway of phenethylamine synthesis was reconstructed by chromosomal integration and overexpression of the Enterococcus faecium pdc gene encoding phenylalanine decarboxylase in Escherichia coli. The genes encoding 3-deoxy-d-arabinoheptulosonate-7-phosphate synthase (aroG), shikimate kinase II (aroL), chorismate mutase/prephenate dehydratase (pheA), and tyrosine aminotransferase (tyrB) in the phenethylamine synthetic pathway were sequentially chromosomally overexpressed. The phosphotransferase system was replaced by deleting the ptsH-ptsI-crr genes and chromosomally overexpressing the genes encoding galactose permease (galP) and glucokinase (glk). In addition, the zwf gene encoding glucose-6-phosphate dehydrogenase in the pentose phosphate pathway was chromosomally overexpressed, generating the final engineered E. coli strain AUD9. The AUD9 strain produced 2.65 g L-1 phenethylamine with a yield of 0.27 g of phenethylamine g-1 glucose in batch fermentation; fed-batch fermentation of AUD9 produced 38.82 g L-1 phenethylamine with a productivity of 1.08 g L-1 h-1 phenethylamine, demonstrating its potential for industrial fermentative production of phenethylamine.

Altered conformational sampling along an evolutionary trajectory changes the catalytic activity of an enzyme.

Several enzymes are known to have evolved from non-catalytic proteins such as solute-binding proteins (SBPs). Although attention has been focused on how a binding site can evolve to become catalytic, an equally important question is: how do the structural dynamics of a binding protein change as it becomes an efficient enzyme? Here we performed a variety of experiments, including propargyl-DO3A-Gd(III) tagging and double electron-electron resonance (DEER) to study the rigid body protein dynamics of reconstructed evolutionary intermediates to determine how the conformational sampling of a protein changes along an evolutionary trajectory linking an arginine SBP to a cyclohexadienyl dehydratase (CDT). We observed that primitive dehydratases predominantly populate catalytically unproductive conformations that are vestiges of their ancestral SBP function. Non-productive conformational states, including a wide-open state, are frozen out of the conformational landscape via remote mutations, eventually leading to extant CDT that exclusively samples catalytically relevant compact states. These results show that remote mutations can reshape the global conformational landscape of an enzyme as a mechanism for increasing catalytic activity.

Completion of the cytosolic post-chorismate phenylalanine biosynthetic pathway in plants.

In addition to being a vital component of proteins, phenylalanine is also a precursor of numerous aromatic primary and secondary metabolites with broad physiological functions. In plants phenylalanine is synthesized predominantly via the arogenate pathway in plastids. Here, we describe the structure, molecular players and subcellular localization of a microbial-like phenylpyruvate pathway for phenylalanine biosynthesis in plants. Using a reverse genetic approach and metabolic flux analysis, we provide evidence that the cytosolic chorismate mutase is responsible for directing carbon flux towards cytosolic phenylalanine production via the phenylpyruvate pathway. We also show that an alternative transcription start site of a known plastidial enzyme produces a functional cytosolic prephenate dehydratase that catalyzes the conversion of prephenate to phenylpyruvate, the intermediate step between chorismate mutase and phenylpyruvate aminotransferase. Thus, our results complete elucidation of phenylalanine biosynthesis via phenylpyruvate in plants, showing that this pathway splits from the known plastidial arogenate pathway at chorismate, instead of prephenate as previously thought, and the complete pathway is localized in the cytosol.

Metabolic engineering of Escherichia coli for L-tyrosine production by expression of genes coding for the chorismate mutase domain of the native chorismate mutase-prephenate dehydratase and a cyclohexadienyl dehydrogenase from Zymomonas mobilis.

The expression of the feedback inhibition-insensitive enzyme cyclohexadienyl dehydrogenase (TyrC) from Zymomonas mobilis and the chorismate mutase domain from native chorismate mutase-prephenate dehydratase (PheA(CM)) from Escherichia coli was compared to the expression of native feedback inhibition-sensitive chorismate mutase-prephenate dehydrogenase (CM-TyrA(p)) with regard to the capacity to produce l-tyrosine in E. coli strains modified to increase the carbon flow to chorismate. Shake flask experiments showed that TyrC increased the yield of l-tyrosine from glucose (Y(l-Tyr/Glc)) by 6.8-fold compared to the yield obtained with CM-TyrA(p). In bioreactor experiments, a strain expressing both TyrC and PheA(CM) produced 3 g/liter of l-tyrosine with a Y(l-Tyr/Glc) of 66 mg/g. These values are 46 and 48% higher than the values for a strain expressing only TyrC. The results show that the feedback inhibition-insensitive enzymes can be employed for strain development as part of a metabolic engineering strategy for l-tyrosine production.

Redesigning the PheA domain of gramicidin synthetase leads to a new understanding of the enzyme's mechanism and selectivity.

The PheA domain of gramicidin synthetase A, a non-ribosomal peptide synthetase, selectively binds phenylalanine along with ATP and Mg2+ and catalyzes the formation of an aminoacyl adenylate. In this study, we have used a novel protein redesign algorithm, K*, to predict mutations in PheA that should exhibit improved binding for tyrosine. Interestingly, the introduction of two predicted mutations to PheA did not significantly improve KD, as measured by equilibrium fluorescence quenching. However, the mutations improved the specificity of the enzyme for tyrosine (as measured by kcat/KM), primarily driven by a 56-fold improvement in KM, although the improvement did not make tyrosine the preferred substrate over phenylalanine. Using stopped-flow fluorometry, we examined binding of different amino acid substrates to the wild-type and mutant enzymes in the pre-steady state in order to understand the improvement in KM. Through these investigations, it became evident that substrate binding to the wild-type enzyme is more complex than previously described. These experiments show that the wild-type enzyme binds phenylalanine in a kinetically selective manner; no other amino acids tested appeared to bind the enzyme in the early time frame examined (500 ms). Furthermore, experiments with PheA, phenylalanine, and ATP reveal a two-step binding process, suggesting that the PheA-ATP-phenylalanine complex may undergo a conformational change toward a catalytically relevant intermediate on the pathway to adenylation; experiments with PheA, phenylalanine, and other nucleotides exhibit only a one-step binding process. The improvement in KM for the mutant enzyme toward tyrosine, as predicted by K*, may indicate that redesigning the side-chain binding pocket allows the substrate backbone to adopt productive conformations for catalysis but that further improvements may be afforded by modeling an enzyme:ATP:substrate complex, which is capable of undergoing conformational change.

Crystallization and preliminary X-ray diffraction analysis of prephenate dehydratase from Mycobacterium tuberculosis H37Rv.

Tuberculosis remains the leading cause of mortality arising from a bacterial pathogen (Mycobacterium tuberculosis). There is an urgent need for the development of new antimycobacterial agents. The aromatic amino-acid pathway is essential for the survival of this pathogen and represents a target for structure-based drug design. Accordingly, the M. tuberculosis prephenate dehydratase has been cloned, expressed, purified and crystallized by the hanging-drop vapour-diffusion method using PEG 400 as a precipitant. The crystal belongs to the orthorhombic space group I222 or I2(1)2(1)2(1), with unit-cell parameters a = 98.26, b = 133.22, c = 225.01 angstroms, and contains four molecules in the asymmetric unit. A complete data set was collected to 3.2 angstroms resolution using a synchrotron-radiation source.

Genome engineering Escherichia coli for L-DOPA overproduction from glucose.

Genome engineering has become a powerful tool for creating useful strains in research and industry. In this study, we applied singleplex and multiplex genome engineering approaches to construct an E. coli strain for the production of L-DOPA from glucose. We first used the singleplex genome engineering approach to create an L-DOPA-producing strain, E. coli DOPA-1, by deleting transcriptional regulators (tyrosine repressor tyrR and carbon storage regulator A csrA), altering glucose transport from the phosphotransferase system (PTS) to ATP-dependent uptake and the phosphorylation system overexpressing galactose permease gene (galP) and glucokinase gene (glk), knocking out glucose-6-phosphate dehydrogenase gene (zwf) and prephenate dehydratase and its leader peptide genes (pheLA) and integrating the fusion protein chimera of the downstream pathway of chorismate. Then, multiplex automated genome engineering (MAGE) based on 23 targets was used to further improve L-DOPA production. The resulting strain, E. coli DOPA-30N, produced 8.67 g/L of L-DOPA in 60 h in a 5 L fed-batch fermentation. This titer is the highest achieved in metabolically engineered E. coli having PHAH activity from glucose.

Genomic analysis of sequence variation in tandemly repeated DNA. Evidence for localized homogeneous sequence domains within arrays of alpha-satellite DNA.

As a model to examine the local distribution of sequence variation within large arrays of tandemly repeated DNA in complex genomes, the long-range organization of alpha-satellite DNA from human chromosome 17 was investigated. Three individual chromosomes, representing different alpha-satellite haplotypes, were segregated into mouse and human somatic cell hybrids and their arrays sized by pulse-field gel electrophoresis. An inventory of the higher-order repeat units found in multiple separate regions of these megabase arrays was obtained using cosmid mapping and two-dimensional gel electrophoresis, a technique that combines the large-scale resolution of pulsed-field gel electrophoresis with the small-scale resolution of conventional gel electrophoresis. These analyses show that alpha-satellite arrays are characterized by the presence of localized homogeneous domains containing only one distinct type of repeat unit. These domains, which consist of sequence variants and/or higher-order repeat length variants, can be up to at least several hundred thousands of bases in length. Both abundant and rare variant repeat units can be localized in these distinct domains, which may correspond to transition states in the evolution of tandem multicopy DNA families. This description of the organization of large arrays of tandem repeats provides insight into mechanisms involved in their homogenization.