Duplication and amplification of toxin genes in Vibrio cholerae

Vibrio cholerae strains of the classical biotype all contain two widely separated copies of the cholera toxin operon ctxAB. In contrast, EI Tor strains containing multiple copies of ctx have their copies arranged on large tandem repeats which are either 7 or 9.7 kb in length. The variation in size among these large tandem duplications was due to a difference in the copy number of a smaller, 2.7 kb, tandemly repeated sequence (RS1) that is located at the novel joint of these duplications, as well as upstream and downstream of ctx. Southern blot hybridization analysis indicated that amplification of a DNA region carrying ctx and flanked by direct repeats of RS1 may be responsible for the hypertoxinogenic phenotype of EI Tor variants selected by intraintestinal growth in rabbits.     

Molecular characterization of a repeat element causing large-scale size variation in the mitochondrial DNA of the sea scallop Placopecten magellanicus

The scallop Placopecten magellanicus has the largest reported animal mitochondrial DNA (average 35 kb) and exhibits large inter- and intraindividual length variation owing to the varying copy number of a repeated element. We have characterized the repeat array by using restriction mapping and sequence analysis. The repeated element consists of 1,442 bp flanked on either side by the sequence ACTTTCC in a direct orientation. The array contains two to eight copies of the repeated element arranged in a direct orientation and in tandem. Only complete copies of the element are present in the array. The repeat element contains three regions with characteristic nucleotide sequences: a 10-bp inverted repeat shown to extrude into a cruciform in a supercoiled DNA plasmid, a 120-bp tract rich in G/C (70%) and adjacent to the inverted repeat, and periodically interspersed homopolymer runs of A and T occurring near the middle of the element which induce DNA curvature in dimeric constructs of the element. The element appears to be unique to P. magellanicus. The structural properties of the repeat element and its organization in an array of repeats may be important in explaining the generation and maintenance of large-scale mitochondrial DNA size variation observed in many animal species.     

Copy Number Variation at the APOL1 Locus

Two coding variants in the APOL1 gene (G1 and G2) explain most of the high rate of kidney disease in African Americans. APOL1-associated kidney disease risk inheritance follows an autosomal recessive pattern: The relative risk of kidney disease associated with inheritance of two high-risk variants is 7–30 fold, depending on the specific kidney phenotype. We wished to determine if the variability in phenotype might in part reflect structural differences in APOL1 gene. We analyzed sequence coverage from 1000 Genomes Project Phase 3 samples as well as exome sequencing data from African American kidney disease cases for copy number variation. 8 samples sequenced in the 1000 Genomes Project showed increased coverage over a ~100kb region that includes APOL2, APOL1 and part of MYH9, suggesting the presence of APOL1 copy number greater than 2. We reasoned that such duplications should be enriched in apparent G1 heterozygotes with kidney disease. Using a PCR-based assay, we observed the presence of this duplication in additional samples from apparent G0G1 or G0G2 individuals. The frequency of this APOL1 duplication was compared among cases (n = 123) and controls (n = 255) with apparent G0G1 heterozygosity. The presence of APOL1 duplication was observed in 4.06% of cases and 0.78% controls, preliminary evidence that this APOL1 duplication may alter susceptibility to kidney disease (p = 0.03). Taqman-based copy number assays confirmed the presence of 3 APOL1 copies in individuals positive for this specific duplication by PCR assay, but also identified a small number of individuals with additional APOL1 copies of presumably different structure. These observations motivate further studies to better assess the contribution of APOL1 copy number on kidney disease risk and on APOL1 function. Investigators and clinicians genotyping APOL1 should also consider whether the particular genotyping platform used is subject to technical errors when more than two copies of APOL1 are present.     

Generation of Tandem Direct Duplications by Reversed-Ends Transposition of Maize Ac Elements

Tandem direct duplications are a common feature of the genomes of eukaryotes ranging from yeast to human, where they comprise a significant fraction of copy number variations. The prevailing model for the formation of tandem direct duplications is non-allelic homologous recombination (NAHR). Here we report the isolation of a series of duplications and reciprocal deletions isolated de novo from a maize allele containing two Class II Ac/Ds transposons. The duplication/deletion structures suggest that they were generated by alternative transposition reactions involving the termini of two nearby transposable elements. The deletion/duplication breakpoint junctions contain 8 bp target site duplications characteristic of Ac/Ds transposition events, confirming their formation directly by an alternative transposition mechanism. Tandem direct duplications and reciprocal deletions were generated at a relatively high frequency (∼0.5 to 1%) in the materials examined here in which transposons are positioned nearby each other in appropriate orientation; frequencies would likely be much lower in other genotypes. To test whether this mechanism may have contributed to maize genome evolution, we analyzed sequences flanking Ac/Ds and other hAT family transposons and identified three small tandem direct duplications with the structural features predicted by the alternative transposition mechanism. Together these results show that some class II transposons are capable of directly inducing tandem sequence duplications, and that this activity has contributed to the evolution of the maize genome.     

Large Inverted Duplications in the Human Genome Form via a Fold-Back Mechanism

Inverted duplications are a common type of copy number variation (CNV) in germline and somatic genomes. Large duplications that include many genes can lead to both neurodevelopmental phenotypes in children and gene amplifications in tumors. There are several models for inverted duplication formation, most of which include a dicentric chromosome intermediate followed by breakage-fusion-bridge (BFB) cycles, but the mechanisms that give rise to the inverted dicentric chromosome in most inverted duplications remain unknown. Here we have combined high-resolution array CGH, custom sequence capture, next-generation sequencing, and long-range PCR to analyze the breakpoints of 50 nonrecurrent inverted duplications in patients with intellectual disability, autism, and congenital anomalies. For half of the rearrangements in our study, we sequenced at least one breakpoint junction. Sequence analysis of breakpoint junctions reveals a normal-copy disomic spacer between inverted and non-inverted copies of the duplication. Further, short inverted sequences are present at the boundary of the disomic spacer and the inverted duplication. These data support a mechanism of inverted duplication formation whereby a chromosome with a double-strand break intrastrand pairs with itself to form a “fold-back” intermediate that, after DNA replication, produces a dicentric inverted chromosome with a disomic spacer corresponding to the site of the fold-back loop. This process can lead to inverted duplications adjacent to terminal deletions, inverted duplications juxtaposed to translocations, and inverted duplication ring chromosomes.     

Unequal Crossing Over Is the Principal Pathway of Homologous Recombination in Tandem Duplications of <Emphasis Type="Italic">Escherichia coli</Emphasis>

Homologous recombination between direct DNA repeats in tandem duplications usually leads to their dissociation. An even number of crossovers between two copies of a duplication should lead to the formation of diploid segregants, i.e., to the preservation of the duplication. However, in studies of the genotype of diploid segregants in heterozygous tandem duplications of Escherichia coli, it was shown that they arise by unequal exchanges between sister chromosomes rather than by intrachromosomal exchanges. Generally, these exchanges lead to the establishment of the homozygous state of (heterozygous) duplications. Since the available data suggest that the exchange between sister chromosomes may be coupled with DNA replication, it is supposed that unequal exchanges between direct DNA repeats occur in the process of DNA replication.__________Translated from Genetika, Vol. 41, No. 8, 2005, pp. 1038–1044.Original Russian Text Copyright © 2005 by Prokop’ev, Sukhodolets.     

Plastidic trnFUUC pseudogenes in North American genus Boechera (Brassicaceae): mechanistic aspects of evolution

The origin and maintenance of a plastidic tandem repeat next to the TRNF (UUC) gene were analyzed in the genus BOECHERA in a phylogenetic context and were compared to published analogous examples that emerged in parallel in the Asteraceae and Juncaceae, respectively. Although we identified some features common to these taxonomic groups with respect to structure and origin of the region, obvious differences were encountered, which argue against a specific mechanism or evolutionary principle underlying the parallel origin and maintenance of the TRNF-tandem repeats in those families. In contrast to the situation in the Asteraceae, no reciprocal recombinant repeat types have been observed in the Brassicaceae. Forty copy types, classified into three groups, were isolated from 103 chloroplast haplotypes of BOECHERA and it was demonstrated that they are composed of four subregions of various origins. We discuss various mutation mechanisms such as DNA replication slippage, and inter- and intrachromosomal recombination which were reported to mediate variation in copy numbers and other types of observed sequence length polymorphism. It is shown that the observed molecular structure of the tandem repeat region did not fully fit the particular patterns expected under a scenario of evolution including any of the known mechanisms. Nevertheless, it appeared that intermolecular unequal crossing-over is most likely the driving force in the evolution of this tandem repeat. However, it remains to be explained, why no reciprocal recombinant copy types have been observed. The reconstructed phylogenetic relationships among copies reflected different evolutionary scenarios as follows: (1) A single and ancient origin of copies pre-dates the radiation of BOECHERA. (2) Parallel expansion and shortening of the tandem repeat within different BOECHERA lineages. (3) Conservation of the first copy, as it was the only one present in all chloroplast haplotypes.     

The structure and early evolution of recently arisen gene duplicates in the Caenorhabditis elegans genome

The significance of gene duplication in provisioning raw materials for the evolution of genomic diversity is widely recognized, but the early evolutionary dynamics of duplicate genes remain obscure. To elucidate the structural characteristics of newly arisen gene duplicates at infancy and their subsequent evolutionary properties, we analyzed gene pairs with < or =10% divergence at synonymous sites within the genome of Caenorhabditis elegans. Structural heterogeneity between duplicate copies is present very early in their evolutionary history and is maintained over longer evolutionary timescales, suggesting that duplications across gene boundaries in conjunction with shuffling events have at least as much potential to contribute to long-term evolution as do fully redundant (complete) duplicates. The median duplication span of 1.4 kb falls short of the average gene length in C. elegans (2.5 kb), suggesting that partial gene duplications are frequent. Most gene duplicates reside close to the parent copy at inception, often as tandem inverted loci, and appear to disperse in the genome as they age, as a result of reduced survivorship of duplicates located in proximity to the ancestral copy. We propose that illegitimate recombination events leading to inverted duplications play a disproportionately large role in gene duplication within this genome in comparison with other mechanisms.     

Expansion of protein domain repeats

Many proteins, especially in eukaryotes, contain tandem repeats of several domains from the same family. These repeats have a variety of binding properties and are involved in protein–protein interactions as well as binding to other ligands such as DNA and RNA. The rapid expansion of protein domain repeats is assumed to have evolved through internal tandem duplications. However, the exact mechanisms behind these tandem duplications are not well-understood. Here, we have studied the evolution, function, protein structure, gene structure, and phylogenetic distribution of domain repeats. For this purpose we have assigned Pfam-A domain families to 24 proteomes with more sensitive domain assignments in the repeat regions. These assignments confirmed previous findings that eukaryotes, and in particular vertebrates, contain a much higher fraction of proteins with repeats compared with prokaryotes. The internal sequence similarity in each protein revealed that the domain repeats are often expanded through duplications of several domains at a time, while the duplication of one domain is less common. Many of the repeats appear to have been duplicated in the middle of the repeat region. This is in strong contrast to the evolution of other proteins that mainly works through additions of single domains at either terminus. Further, we found that some domain families show distinct duplication patterns, e.g., nebulin domains have mainly been expanded with a unit of seven domains at a time, while duplications of other domain families involve varying numbers of domains. Finally, no common mechanism for the expansion of all repeats could be detected. We found that the duplication patterns show no dependence on the size of the domains. Further, repeat expansion in some families can possibly be explained by shuffling of exons. However, exon shuffling could not have created all repeats.     

Evolution of Repeated Sequence Arrays in the D-Loop Region of Bat Mitochondrial DNA

Analysis of mitochondrial DNA control region sequences from 41 species of bats representing 11 families revealed that repeated sequence arrays near the tRNA-Pro gene are present in all vespertilionine bats. Across 18 species tandem repeats varied in size from 78 to 85 bp and contained two to nine repeats. Heteroplasmy ranged from 15% to 63%. Fewer repeats among heteroplasmic than homoplasmic individuals in a species with up to nine repeats indicates selection may act against long arrays. A lower limit of two repeats and more repeats among heteroplasmic than homoplasmic individuals in two species with few repeats suggests length mutations are biased. Significant regressions of heteroplasmy, θ and π, on repeat number further suggest that repeat duplication rate increases with repeat number. Comparison of vespertilionine bat consensus repeats to mammal control region sequences revealed that tandem repeats of similar size, sequence and number also occur in shrews, cats and bighorn sheep. The presence of two conserved protein-binding sequences in all repeat units indicates that convergent evolution has occurred by duplication of functional units. We speculate that D-loop region tandem repeats may provide signal redundancy and a primitive repair mechanism in the event of somatic mutations to these binding sites.     

Intra-Genomic Variation in the Ribosomal Repeats of Nematodes

Ribosomal loci represent a major tool for investigating environmental diversity and community structure via high-throughput marker gene studies of eukaryotes (e.g. 18S rRNA). Since the estimation of species’ abundance is a major goal of environmental studies (by counting numbers of sequences), understanding the patterns of rRNA copy number across species will be critical for informing such high-throughput approaches. Such knowledge is critical, given that ribosomal RNA genes exist within multi-copy repeated arrays in a genome. Here we measured the repeat copy number for six nematode species by mapping the sequences from whole genome shotgun libraries against reference sequences for their rRNA repeat. This revealed a 6-fold variation in repeat copy number amongst taxa investigated, with levels of intragenomic variation ranging from 56 to 323 copies of the rRNA array. By applying the same approach to four C. elegans mutation accumulation lines propagated by repeated bottlenecking for an average of ~400 generations, we find on average a 2-fold increase in repeat copy number (rate of increase in rRNA estimated at 0.0285-0.3414 copies per generation), suggesting that rRNA repeat copy number is subject to selection. Within each Caenorhabditis species, the majority of intragenomic variation found across the rRNA repeat was observed within gene regions (18S, 28S, 5.8S), suggesting that such intragenomic variation is not a product of selection for rRNA coding function. We find that the dramatic variation in repeat copy number among these six nematode genomes would limit the use of rRNA in estimates of organismal abundance. In addition, the unique pattern of variation within a single genome was uncorrelated with patterns of divergence between species, reflecting a strong signature of natural selection for rRNA function. A better understanding of the factors that control or affect copy number in these arrays, as well as their rates and patterns of evolution, will be critical for informing estimates of global biodiversity.     

Heteroplasmy of Short Tandem Repeats in Mitochondrial DNA of Atlantic Cod, Gadus Morhua

The mitochondrial DNA of the Atlantic cod (Gadus morhua) contains a tandem array of 40-bp repeats in the D-loop region of the molecule. Variation among molecules in the copy number of these repeats results in mtDNA length variation and heteroplasmy (the presence of more than one form of mtDNA in an individual). In a sample of fish collected from different localities around Iceland and off George's Bank, each individual was heteroplasmic for two or more mtDNAs ranging in repeat copy number from two (common) to six (rare). An earlier report on mtDNA heteroplasmy in sturgeon (Acipenser transmontanus) presented a competitive displacement model for length mutations in mtDNAs containing tandem arrays and the cod data deviate from this model. Depending on the nature of putative secondary structures and the location of D-loop strand termination, additional mechanisms of length mutation may be needed to explain the range of mtDNA length variants maintained in these populations. The balance between genetic drift and mutation in maintaining this length polymorphism is estimated through a hierarchical analysis of diversity of mtDNA length variation in the Iceland samples. Eighty percent of the diversity lies within individuals, 8% among individuals and 12% among localities. An estimate of theta = 2N(eo) mu greater than 1 indicates that this system is characterized by a high mutation rate and is governed primarily by deterministic dynamics. The sequences of repeat arrays from fish collected in Norway, Iceland and George's Bank show no nucleotide variation suggesting that there is very little substructuring to the North Atlantic cod population.     

Evolutionary dynamics of tandem repeats in the mitochondrial DNA control region of the minnow Cyprinella spiloptera

Length variation due to tandem repeats is now recognized as a common feature of animal mitochondrial DNA; however, the evolutionary dynamics of repeated sequences are not well understood. Using phylogenetic analysis, predictions of three models of repeat evolution were tested for arrays of 260-bp repeats in the cyprinid fish Cyprinella spiloptera. Variation at different nucleotide positions in individual repeats supported different models of repeat evolution. One set of characters included several nucleotide variants found in all copies from a limited number of individuals, while the other set included an 8- bp deletion found in a limited number of copies in all individuals. The deletion and an associated nucleotide change appear to be the result of a deterministic, rather than stochastic, mutation process. Parallel origins of repeat arrays in different mitochondrial lineages, possibly coupled with a homogenization mechanism, best explain the distribution of nucleotide variation.      

Unequal crossing over is the principal pathway of homologous recombination in tandem duplications of Escherichia coli

Homologous recombination between direct DNA repeats in tandem duplications usually leads to their dissociation. An even number of crossovers between two copies of a duplication should lead to the formation of diploid segregants, i.e., to the preservation of the duplication. However, in studies of the genotype of diploid segregants in heterozygous tandem duplications of Escherichia coli, it was shown that they arise by unequal exchanges between sister chromosomes rather than by intrachromosomal exchanges. Generally, these exchanges lead to the establishment of the homozygous state of (heterozygous) duplications. Since the available data suggest that the exchange between sister chromosomes may be coupled with DNA replication, it is supposed that unequal exchanges between direct DNA repeats occur in the process of DNA replication.     

Nebulin: A Study of Protein Repeat Evolution

Protein domain repeats are common in proteins that are central to the organization of a cell, in particular in eukaryotes. They are known to evolve through internal tandem duplications. However, the understanding of the underlying mechanisms is incomplete. To shed light on repeat expansion mechanisms, we have studied the evolution of the muscle protein Nebulin, a protein that contains a large number of actin-binding nebulin domains.Nebulin proteins have evolved from an invertebrate precursor containing two nebulin domains. Repeat regions have expanded through duplications of single domains, as well as duplications of a super repeat (SR) consisting of seven nebulins. We show that the SR has evolved independently into large regions in at least three instances: twice in the invertebrate Branchiostoma floridae and once in vertebrates.In-depth analysis reveals several recent tandem duplications in the Nebulin gene. The events involve both single-domain and multidomain SR units or several SR units. There are single events, but frequently the same unit is duplicated multiple times. For instance, an ancestor of human and chimpanzee underwent two tandem duplications. The duplication junction coincides with an Alu transposon, thus suggesting duplication through Alu-mediated homologous recombination.Duplications in the SR region consistently involve multiples of seven domains. However, the exact unit that is duplicated varies both between species and within species. Thus, multiple tandem duplications of the same motif did not create the large Nebulin protein.Finally, analysis of segmental duplications in the human genome reveals that duplications are more common in genes containing domain repeats than in those coding for nonrepeated proteins. In fact, segmental duplications are found three to six times more often in long repeated genes than expected by chance.     

Digital Genotyping of Macrosatellites and Multicopy Genes Reveals Novel Biological Functions Associated with Copy Number Variation of Large Tandem Repeats

Tandem repeats are common in eukaryotic genomes, but due to difficulties in assaying them remain poorly studied. Here, we demonstrate the utility of Nanostring technology as a targeted approach to perform accurate measurement of tandem repeats even at extremely high copy number, and apply this technology to genotype 165 HapMap samples from three different populations and five species of non-human primates. We observed extreme variability in copy number of tandemly repeated genes, with many loci showing 5–10 fold variation in copy number among humans. Many of these loci show hallmarks of genome assembly errors, and the true copy number of many large tandem repeats is significantly under-represented even in the high quality ‘finished’ human reference assembly. Importantly, we demonstrate that most large tandem repeat variations are not tagged by nearby SNPs, and are therefore essentially invisible to SNP-based GWAS approaches. Using association analysis we identify many cis correlations of large tandem repeat variants with nearby gene expression and DNA methylation levels, indicating that variations of tandem repeat length are associated with functional effects on the local genomic environment. This includes an example where expansion of a macrosatellite repeat is associated with increased DNA methylation and suppression of nearby gene expression, suggesting a mechanism termed “repeat induced gene silencing”, which has previously been observed only in transgenic organisms. We also observed multiple signatures consistent with altered selective pressures at tandemly repeated loci, suggesting important biological functions. Our studies show that tandemly repeated loci represent a highly variable fraction of the genome that have been systematically ignored by most previous studies, copy number variation of which can exert functionally significant effects. We suggest that future studies of tandem repeat loci will lead to many novel insights into their role in modulating both genomic and phenotypic diversity.     

Evolution of CC chemokines in teleost fish: a case study in gene duplication and implications for immune diversity

Chemokines are a superfamily of cytokines responsible for regulating cell migration under both inflammatory and physiological conditions. CC chemokines are the largest subfamily of chemokines, with 28 members in humans. A subject of intense study in mammalian species, the known functional roles of CC chemokines ligands in both developmental and disease conditions continue to expand. They are also an important family for the study of gene copy number variation and tandem duplication in mammalian species. However, little is known regarding the evolutionary origin and status of these ligands in primitive vertebrates such as teleost fish. In this paper, we review the evolution of the teleost fish CC chemokine gene family, noting evidence of widespread tandem gene duplications and examining the implications of this phenomenon on immune diversity. Through extensive phylogenetic analysis of the CC chemokine sets of four teleost species, zebrafish, catfish, rainbow trout, and Atlantic salmon, we identified seven large groups of CC chemokines. It appeared that several major groups of CC chemokines are highly related including the CCL19/21/25 group, the CCL20 group, CCL27/28 group, and the fish-specific group. In the three remaining groups that contained the largest number of members, the CCL17/22 group, the MIP group, and the MCP group, similarities among species members were obscured by rapid, tandem duplications that may contribute to immune diversity.     

Patterns of evolution in the unique tRNA gene arrays of the genus Entamoeba

Genome sequencing of the protistan parasite Entamoeba histolytica HM-1:IMSS revealed that almost all the tRNA genes are organized into tandem arrays that make up over 10% of the genome. The 25 distinct array units contain up to 5 tRNA genes each and some also encode the 5S RNA. Between adjacent genes in array units are complex short tandem repeats (STRs) resembling microsatellites. To investigate the origins and evolution of this unique gene organization, we have undertaken a genome survey to determine the array unit organization in 4 other species of Entamoeba-Entamoeba dispar, Entamoeba moshkovskii, Entamoeba terrapinae, and Entamoeba invadens-and have explored the STR structure in other isolates of E. histolytica. The genome surveys revealed that E. dispar has the same array unit organization as E. histolytica, including the presence and numerical variation of STRs between adjacent genes. However, the individual repeat sequences are completely different to those in E. histolytica. All other species of Entamoeba studied also have tandem arrays of clustered tRNA genes, but the gene composition of the array units often differs from that in E. histolytica/E. dispar. None of the other species' arrays exhibit the complex STRs between adjacent genes although simple tandem duplications are occasionally seen. The degree of similarity in organization reflects the phylogenetic relationships among the species studied. Within individual isolates of E. histolytica most copies of the array unit are uniform in sequence with only minor variation in the number and organization of the STRs. Between isolates, however, substantial differences in STR number and organization can exist although the individual repeat sequences tend to be conserved. The origin of this unique gene organization in the genus Entamoeba clearly predates the common ancestor of the species investigated to date and their function remains unclear.     

Structure,molecular evolution and maintenance of copy number of extended repeated structures in the X-heterochromatin of Drosophila melanogaster

The 60 kb repeats located in the distal heterochromatin of the X chromosome of Drosophila melanogaster were cloned in overlapping cosmids. These regions, designated as SCLRs, comprised the following types of repeated elements Stellate genes, which are known to be involved in spermatogenesis; copia-like retrotransposons; LINE elements, including amplified Type rDNA insertions; and rDNA fragments. The following steps in SCLR formation were hypothesized: insertion of mobile elements into the rDNA and Stellate gene clusters: internal tandem duplication events; recombination between the rDNA cluster and Stellate tandem repeat; and amplification of the whole SCLR structure. There are about nine SCLR copies per haploid genome, but there is approximately a twofold variation in copy number between fly stocks. The SCLR copy number differences between closely related stocks are suggested to be the result of unequal sister chromatid exchange (USCE). The restricted variation in SCLR copy number between unrelated stocks and the absence of chromosomes free of SCLRs suggests that natural selection is active in copy number maintenance.