Structure and Conservation of Tandem Repeats in the Mitochondrial DNA Control Region of the Least Brook Lamprey (<Emphasis Type="Italic">Lampetra aepyptera</Emphasis>)

Tandemly repeated sequences are a common feature of vertebrate mitochondrial DNA control regions. However, questions still remain about their mode of evolution, function, and phylogenetic distribution. We report phylogenetic and geographic patterns of variation of control region repeat sequence and number in a nonparasitic lamprey, Lampetra aepyptera. A survey of populations from throughout the species’ range revealed remarkably low repeat sequence polymorphism but some interpopulation variation in repeat number. The high sequence similarity extended to repeats observed in other species in the genus Lampetra and other lamprey genera. The very low levels of variation suggest a high copy turnover. Our data are consistent with the illegitimate elongation model of repeat gain and loss and further suggest that repeat change occurs at internal copies. However, the limited variation across some species of lamprey suggests that functional constraints may further limit variation.     

Patterns of DNA-Sequence Divergence Between <Emphasis Type="Italic">Drosophila miranda</Emphasis> and <Emphasis Type="Italic">D. pseudoobscura</Emphasis>

Contrary to the classical view, a large amount of non-coding DNA seems to be selectively constrained in Drosophila and other species. Here, using Drosophila miranda BAC sequences and the Drosophila pseudoobscura genome sequence, we aligned coding and non-coding sequences between D. pseudoobscura and D. miranda, and investigated their patterns of evolution. We found two patterns that have previously been observed in comparisons between Drosophila melanogaster and its relatives. First, there is a negative correlation between intron divergence and intron length, suggesting that longer non-coding sequences may contain more regulatory elements than shorter sequences. Our other main finding is a negative correlation between the rate of non-synonymous substitutions (d _N) and codon usage bias (F _op), showing that fast-evolving genes have a lower codon usage bias, consistent with strong positive selection interfering with weak selection for codon usage.     

Connecting recombination,nucleotide diversity,and species divergence in Drosophila

The association between recombination rate and nucleotide diversity provides compelling evidence for the action of natural selection across much of the Drosophila melanogaster genome. This conclusion is further supported by the lack of association between recombination rate and nucleotide divergence between species. However, studies of other species, including other Drosophila, have not always yielded the same results. Our recent study measured these parameters within the D. pseudoobscura species group using next-generation sequencing and high-throughput genotyping technologies. We documented fine-scale variation in crossover rate within D. pseudoobscura, and we observed that crossover variation was strongly associated with nucleotide diversity only when measured at a fine-scale. We also observed associations between crossover rate and sequence differences between D. pseudoobscura and its close relatives. These latter associations could have been driven in part by mutagenic effects associated with double-strand break repair, but we cannot exclude the possibility that it results primarily from shared ancestral polymorphisms. Overall, this work strongly underscores the importance of scale in testing for associations of recombination rate with other parameters, and it brings us one small step closer to understanding the role of natural selection and other evolutionary forces in shaping divergence among genomes.     

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.     

Weak preservation of local neutral substitution rates across mammalian genomes

Background  

The rate at which neutral (non-functional) bases undergo substitution is highly dependent on their location within a genome. However, it is not clear how fast these location-dependent rates change, or to what extent the substitution rate patterns are conserved between lineages. To address this question, which is critical not only for understanding the substitution process but also for evaluating phylogenetic footprinting algorithms, we examine ancestral repeats: a predominantly neutral dataset with a significantly higher genomic density than other datasets commonly used to study substitution rate variation. Using this repeat data, we measure the extent to which orthologous ancestral repeat sequences exhibit similar substitution patterns in separate mammalian lineages, allowing us to ascertain how well local substitution rates have been preserved across species.     

Slip-Sliding Away: Serial Changes and Homoplasy in Repeat Number in the Drosophila yakuba Homolog of Human Cancer Susceptibility Gene BRCA2

Several recent studies have examined the function and evolution of a Drosophila homolog to the human breast cancer susceptibility gene BRCA2, named dmbrca2. We previously identified what appeared to be a recent expansion in the RAD51-binding BRC-repeat array in the ancestor of Drosophila yakuba. In this study, we examine patterns of variation and evolution of the dmbrca2 BRC-repeat array within D. yakuba and its close relatives. We develop a model of how unequal crossing over may have produced the expanded form, but we also observe short repeat forms, typical of other species in the D. melanogaster group, segregating within D. yakuba and D. santomea. These short forms do not appear to be identical-by-descent, suggesting that the history of dmbrca2 in the D. melanogaster subgroup has involved repeat unit contractions resulting in homoplasious forms. We conclude that the evolutionary history of dmbrca2 in D. yakuba and perhaps in other Drosophila species may be more complicated than can be inferred from examination of the published single genome sequences per species.     

Polymorphism of canonical and noncanonical <Emphasis Type="Italic">gypsy</Emphasis> sequences in different species of <Emphasis Type="Italic">Drosophila </Emphasis><Emphasis Type="Italic">melanogaster</Emphasis> subgroup: possible evolutionary relations

Mobile genetic elements constitute a substantial part of eukaryotic genome and play an important role in its organization and functioning. Co-evolution of retrotransposons and their hosts resulted in the establishment of control systems employing mechanisms of RNA interference that seem to be impossible to evade. However, “active” copies of endogenous retrovirus gypsy escape cellular control in some cases, while its evolutionary elder “inactive” variants do not. To clarify the evolutionary relationship between “active” and “inactive” gypsy we combined two approaches: the analysis of gypsy sequences, isolated from G32 Drosophila melanogaster strain and from different Drosophila species of the melanogaster subgroup, as well as the study of databases, available on the Internet. No signs of “intermediate” (between “active” and “inactive”) gypsy form were found in GenBank, and four full-size G32 gypsy copies demonstrated a convergence that presumably involves gene conversion. No “active” gypsy were revealed among PCR generated gypsy ORF3 sequences from the various Drosophila species indicating that “active” gypsy appeared in some population of D. melanogaster and then started to spread out. Analysis of sequences flanking gypsy variants in G32 revealed their predominantly heterochromatic location. Discrepancy between the structure of actual gypsy sites in G32 and corresponding sequences in database might indicate significant inter-strain heterochromatin diversity. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Fine-Scale Crossover Rate Heterogeneity in <Emphasis Type="Italic">Drosophila pseudoobscura</Emphasis>

Broad-scale differences in crossover rate across the genome have been characterized in most genomes studied. Fine-scale differences, however, have only been examined in a few taxa, such as Arabidopsis, yeast, humans, and mice. No prior studies have directly looked for fine-scale recombination rate heterogeneity in Drosophila. We produced 370 Drosophila pseudoobscura containing a crossover event within the 2-megabase (MB) region between the genes yellow and white. We then examined 19 intervals within this region and determined where the crossovers occurred. We found that recombination events occur nonrandomly on a small scale and that mild “hotspots“ of a few kilobases exist in Drosophila. Among the regions studied, recombination rates varied from 1.4 to 52 cM/MB. We also observed a trend toward high codon bias in regions of high recombination. Finally, we identified a significantly positive correlation between recombination rate and simple repeats, as well as the motif CACAC. These sequence features may contribute to broad-scale variation in crossover rate and, thus, shed light on features associated with crossover rate heterogeneity at a genome-wide scale. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Dmitri Petrov]     

In Vitro Repair of Gaps in Bacteriophage T7 DNA

An in vitro system based upon extracts of Escherichia coli infected with bacteriophage T7 was used to study the mechanism of double-strand break repair. Double-strand breaks were placed in T7 genomes by cutting with a restriction endonuclease which recognizes a unique site in the T7 genome. These molecules were allowed to repair under conditions where the double-strand break could be healed by (i) direct joining of the two partial genomes resulting from the break, (ii) annealing of complementary versions of 17-bp sequences repeated on either side of the break, or (iii) recombination with intact T7 DNA molecules. The data show that while direct joining and single-strand annealing contributed to repair of double-strand breaks, these mechanisms made only minor contributions. The efficiency of repair was greatly enhanced when DNA molecules that bridge the region of the double-strand break (referred to as donor DNA) were provided in the reaction mixtures. Moreover, in the presence of the donor DNA most of the repaired molecules acquired genetic markers from the donor DNA, implying that recombination between the DNA molecules was instrumental in repairing the break. Double-strand break repair in this system is highly efficient, with more than 50% of the broken molecules being repaired within 30 min under some experimental conditions. Gaps of 1,600 nucleotides were repaired nearly as well as simple double-strand breaks. Perfect homology between the DNA sequence near the break site and the donor DNA resulted in minor (twofold) improvement in the efficiency of repair. However, double-strand break repair was still highly efficient when there were inhomogeneities between the ends created by the double-strand break and the T7 genome or between the ends of the donor DNA molecules and the genome. The distance between the double-strand break and the ends of the donor DNA molecule was critical to the repair efficiency. The data argue that ends of DNA molecules formed by double-strand breaks are typically digested by between 150 and 500 nucleotides to form a gap that is subsequently repaired by recombination with other DNA molecules present in the same reaction mixture or infected cell.     

Competitive repair by naturally dispersed repetitive DNA during non-allelic homologous recombination

Genome rearrangements often result from non-allelic homologous recombination (NAHR) between repetitive DNA elements dispersed throughout the genome. Here we systematically analyze NAHR between Ty retrotransposons using a genome-wide approach that exploits unique features of Saccharomyces cerevisiae purebred and Saccharomyces cerevisiae/Saccharomyces bayanus hybrid diploids. We find that DNA double-strand breaks (DSBs) induce NAHR-dependent rearrangements using Ty elements located 12 to 48 kilobases distal to the break site. This break-distal recombination (BDR) occurs frequently, even when allelic recombination can repair the break using the homolog. Robust BDR-dependent NAHR demonstrates that sequences very distal to DSBs can effectively compete with proximal sequences for repair of the break. In addition, our analysis of NAHR partner choice between Ty repeats shows that intrachromosomal Ty partners are preferred despite the abundance of potential interchromosomal Ty partners that share higher sequence identity. This competitive advantage of intrachromosomal Tys results from the relative efficiencies of different NAHR repair pathways. Finally, NAHR generates deleterious rearrangements more frequently when DSBs occur outside rather than within a Ty repeat. These findings yield insights into mechanisms of repeat-mediated genome rearrangements associated with evolution and cancer.     

RETRACTED ARTICLE: The effect of the <Emphasis Type="Italic">mus309</Emphasis> mutation,defective in DNA double-strand break repair,on crossing over in <Emphasis Type="Italic">Drosophila melanogaster</Emphasis> suggests a mechanism for the centromere effect of crossing over

The mus309 gene in Drosophila melanogaster encodes a RecQ helicase which is involved in DNA double-strand break (DSB) repair. In a brood pattern analysis, it was observed that in mus309 mutant females, the frequency of single crossovers in the central cv–v interval of the X chromosome was reduced in young females but returned to the level of the wild type control as the females aged. In the proximal v–f interval, the frequency of single crossovers was increased during the entire experimental period. In particular, it was observed that the frequency of double crossovers, as well as the coefficient of coincidence first increased but then gradually decreased, finally reaching the level of the control flies, as the females aged. Map distances increased due to the mus309 mutation in both gene interval studies, but they did not change as the females aged, a result suggesting that the mus309 gene controls the distribution of DSBs to be repaired as crossovers instead of non-crossovers. The results suggest a mechanism for the centromere effect of crossing over in Drosophila, viz the fact the frequency of meiotic crossing over reduces with the age of the female, and that the reduction is more pronounced the closer the interval is to the proximal heterochromatin of the chromosome arm. According to the model suggested, the centromere effect is simply a matter of the balance between different pathways of the repair of the DSBs of DNA.     

Evolution of the sugar receptors in insects

Background  

Perception of sugars is an invaluable ability for insects which often derive quickly accessible energy from these molecules. A distinctive subfamily of eight proteins within the gustatory receptor (Gr) family has been identified as sugar receptors (SRs) in Drosophila melanogaster (Gr5a, Gr61a, and Gr64a-f). We examined the evolution of these SRs within the 12 available Drosophila genome sequences, as well as three mosquito, two moth, and beetle, bee, and wasp genome sequences.     

Roles of <Emphasis Type="Italic">Saccharomyces cerevisiae RAD17</Emphasis> and <Emphasis Type="Italic">CHK1</Emphasis> checkpoint genes in the repair of double-strand breaks in cycling cells

Checkpoints are components of signalling pathways involved in genome stability. We analysed the putative dual functions of Rad17 and Chk1 as checkpoints and in DNA repair using mutant strains of Saccharomyces cerevisiae. Logarithmic populations of the diploid checkpoint-deficient mutants, chk1Δ/chk1Δ and rad17Δ/rad17Δ, and an isogenic wild-type strain were exposed to the radiomimetic agent bleomycin (BLM). DNA double-strand breaks (DSBs) determined by pulsed-field electrophoresis, surviving fractions, and proliferation kinetics were measured immediately after treatments or after incubation in nutrient medium in the presence or absence of cycloheximide (CHX). The DSBs induced by BLM were reduced in the wild-type strain as a function of incubation time after treatment, with chromosomal repair inhibited by CHX. rad17Δ/rad17Δ cells exposed to low BLM concentrations showed no DSB repair, low survival, and CHX had no effect. Conversely, rad17Δ/rad17Δ cells exposed to high BLM concentrations showed DSB repair inhibited by CHX. chk1Δ/chk1Δ cells showed DSB repair, and CHX had no effect; these cells displayed the lowest survival following high BLM concentrations. Present results indicate that Rad17 is essential for inducible DSB repair after low BLM-concentrations (low levels of oxidative damage). The observations in the chk1Δ/chk1Δ mutant strain suggest that constitutive nonhomologous end-joining is involved in the repair of BLM-induced DSBs. The differential expression of DNA repair and survival in checkpoint mutants as compared to wild-type cells suggests the presence of a regulatory switch-network that controls and channels DSB repair to alternative pathways, depending on the magnitude of the DNA damage and genetic background. Nelson Bracesco and Ema C. Candreva have contributed equally to this article.     

Characterization of sequence contexts that favor alternative end joining at Cas9-induced double-strand breaks

Alternative end joining (alt-EJ) mechanisms, such as polymerase theta-mediated end joining, are increasingly recognized as important contributors to inaccurate double-strand break repair. We previously proposed an alt-EJ model whereby short DNA repeats near a double-strand break anneal to form secondary structures that prime limited DNA synthesis. The nascent DNA then pairs with microhomologous sequences on the other break end. This synthesis-dependent microhomology-mediated end joining (SD-MMEJ) explains many of the alt-EJ repair products recovered following I-SceI nuclease cutting in Drosophila. However, sequence-specific factors that influence SD-MMEJ repair remain to be fully characterized. Here, we expand the utility of the SD-MMEJ model through computational analysis of repair products at Cas9-induced double-strand breaks for 1100 different sequence contexts. We find evidence at single nucleotide resolution for sequence characteristics that drive successful SD-MMEJ repair. These include optimal primer repeat length, distance of repeats from the break, flexibility of DNA sequence between primer repeats, and positioning of microhomology templates relative to preferred primer repeats. In addition, we show that DNA polymerase theta is necessary for most SD-MMEJ repair at Cas9 breaks. The analysis described here includes a computational pipeline that can be utilized to characterize preferred mechanisms of alt-EJ repair in any sequence context.     

Differential repair of excision gaps generated by transposable elements of the ‘Ac family’

Studies on transposable elements of the Ac family have led to different models for excision gap repair in either plants or Drosophila. Excision products generated by the plant transposable elements Ac and Tam3 imply a more or less straightforward ligation of broken ends; excision products of the Drosophila P element indicate the involvement of ‘double-strand break’ (DSB) repair. Recent findings that excision products of Ac and Tam3 can also contain traces of the element ends indicate, however, that DSB repair might be an alternative repair mechanism in plants. A functional DSB repair mechanism in plants can also be deduced from the observed rapid increases of Ac copy number during plant development and from the involvement of Ac in the generation of internal Ac deletions. On the other hand, alternative repair mechanisms may also be functional in Drosophila, because some of the ‘footprints’ generated upon P excision can be explained by a mechanism that has been postulated for excision gap repair in plants. It is concluded that plants and Drosophila can use similar repair mechanisms, but that the predominance of a certain repair mechanism is determined by the host.     

Evolution and selection of Rhg1, a copy‐number variant nematode‐resistance locus

The soybean cyst nematode (SCN) resistance locus Rhg1 is a tandem repeat of a 31.2 kb unit of the soybean genome. Each 31.2‐kb unit contains four genes. One allele of Rhg1, Rhg1‐b, is responsible for protecting most US soybean production from SCN. Whole‐genome sequencing was performed, and PCR assays were developed to investigate allelic variation in sequence and copy number of the Rhg1 locus across a population of soybean germplasm accessions. Four distinct sequences of the 31.2‐kb repeat unit were identified, and some Rhg1 alleles carry up to three different types of repeat unit. The total number of copies of the repeat varies from 1 to 10 per haploid genome. Both copy number and sequence of the repeat correlate with the resistance phenotype, and the Rhg1 locus shows strong signatures of selection. Significant linkage disequilibrium in the genome outside the boundaries of the repeat allowed the Rhg1 genotype to be inferred using high‐density single nucleotide polymorphism genotyping of 15 996 accessions. Over 860 germplasm accessions were found likely to possess Rhg1 alleles. The regions surrounding the repeat show indications of non‐neutral evolution and high genetic variability in populations from different geographic locations, but without evidence of fixation of the resistant genotype. A compelling explanation of these results is that balancing selection is in operation at Rhg1.     

Molecular cloning of mei-41, a gene that influences both somatic and germline chromosome metabolism of Drosophila melanogaster

The mei-41 gene of Drosophila melanogaster plays an essential role in meiosis, in the maintenance of somatic chromosome stability, in postreplication repair and in DNA double-strand break repair. This gene has been cytogenetically localized to polytene chromosome bands 14C4-6 using available chromosomal aberrations. About 60 kb of DNA sequence has been isolated following a bidirectional chromosomal walk that extends over the cytogenetic interval 14C1-6. The breakpoints of chromosomal aberrations identified within that walk establish that the entire mei-41 gene has been cloned. Two independently derived mei-41 mutants have been shown to carry P insertions within a single 2.2 kb fragment of the walk. Since revertants of those mutants have lost the P element sequences, an essential region of the mei-41 gene is present in that fragment. A 10.5 kb genomic fragment that spans the P insertion sites has been found to restore methyl methanesulfonate resistance and female fertility of the mei-41 ^D3 mutants. The results demonstrate that all the sequences required for the proper expression of the mei-41 gene are present on this genomic fragment. This study provides the foundation for molecular analysis of a function that is essential for chromosome stability in both the germline and somatic cells.This Paper is dedicated to the memory of Professor James B. Boyd     

Extended sequence of the turkey MHC <Emphasis Type="Italic">B</Emphasis>-locus and sequence variation in the highly polymorphic B-G loci

Genetic variation in the major histocompatibility complex (MHC) is directly correlated to differences in disease resistance. Immunity is greatly dependent on highly polymorphic genes in the MHC, such as class I, class II, and class III complement genes. Preliminary studies of wild turkey populations show extreme polymorphisms in a family of genes exclusive to the avian MHC, the class IV or B-G genes. Significance of this variation is unclear as there are few and conflicting studies of the expression of these genes. Confounding understanding of B-G variation is the lack of a complete delineation of the number of loci in the turkey genome. Direct 454 sequencing of a clone from the CHORI-260 BAC library was used to extend the turkey MHC B-locus sequence, identifying five additional complete B-locus genes including two B-G loci. Sequences of the new B-G genes were compared with those of other turkey gene (BG1–3) and sequences available for other galliformes. Phylogenetic analysis shows species-specific gene evolution supporting a birth–death model of evolution for the B-G gene family. Analysis of variation within the signal peptide sequence (exon 1) found two clusters of polymorphism among the turkey B-G genes. Resequencing of exon 1 in a diverse sample including wild, heritage, and commercial turkeys confirmed multiple alleles at each B-G gene. Future studies aim to correlate B-G variation with group and individual immunological differences.     

DNA resection at chromosome breaks promotes genome stability by constraining non-allelic homologous recombination

DNA double-strand breaks impact genome stability by triggering many of the large-scale genome rearrangements associated with evolution and cancer. One of the first steps in repairing this damage is 5'→3' resection beginning at the break site. Recently, tools have become available to study the consequences of not extensively resecting double-strand breaks. Here we examine the role of Sgs1- and Exo1-dependent resection on genome stability using a non-selective assay that we previously developed using diploid yeast. We find that Saccharomyces cerevisiae lacking Sgs1 and Exo1 retains a very efficient repair process that is highly mutagenic to genome structure. Specifically, 51% of cells lacking Sgs1 and Exo1 repair a double-strand break using repetitive sequences 12-48 kb distal from the initial break site, thereby generating a genome rearrangement. These Sgs1- and Exo1-independent rearrangements depend partially upon a Rad51-mediated homologous recombination pathway. Furthermore, without resection a robust cell cycle arrest is not activated, allowing a cell with a single double-strand break to divide before repair, potentially yielding multiple progeny each with a different rearrangement. This profusion of rearranged genomes suggests that cells tolerate any dangers associated with extensive resection to inhibit mutagenic pathways such as break-distal recombination. The activation of break-distal recipient repeats and amplification of broken chromosomes when resection is limited raise the possibility that genome regions that are difficult to resect may be hotspots for rearrangements. These results may also explain why mutations in resection machinery are associated with cancer.