Search

Transposable elements (TEs) are repetitive sequences of DNA that replicate and proliferate throughout genomes. Taken together, all the TEs in a genome form a diverse community of sequences, which can be studied to draw conclusions about genome evolution. TE diversity can be measured using models for ecological community diversity that consider species richness and evenness. Several models predict TE diversity decreasing as genomes expand because of selection against ectopic recombination and/or competition among TEs to garner host replicative machinery and evade host silencing mechanisms. Salamanders have some of the largest vertebrate genomes and highest TE loads. Salamanders of the genus Plethodon, in particular, have genomes that range in size from 20 to 70 Gb. Here, we use Oxford Nanopore sequencing to generate low-coverage genomic sequences for four species of Plethodon that encompass two independent genome expansion events, one in the eastern clade (P. cinereus, 29.3 Gb vs. P. glutinosus, 38.9 Gb) and one in the western clade (P. vehiculum, 46.4 Gb vs P. idahoensis, 67.0 Gb). We classified the TEs in these genomes and found ~52 TE superfamilies, accounting for 27-32% of the genomes. We calculated Simpson’s and Shannon’s diversity indices to quantify overall TE diversity. In both pairwise comparisons, the diversity index values for the smaller and larger genome were almost identical. This result indicates that, when genomes reach extremely large sizes, they maintain high levels of TE diversity at the superfamily level, in contrast to predictions made by previous studies on smaller genomes.

The Cre/loxP system is a powerful tool for generating conditional gene knockout (KO) mice and elucidate gene function in vivo. CD19-Cre and Mb1-iCre transgenic mice are commonly used for generating B cell-specific KO mice and investigate the development, as well as the physiological and pathophysiological roles of B cells. However, the CD19-Cre line low efficiency and the Mb1-iCre line occasional ectopic recombination represent challenges for their use. Thus, we developed a CD19-codon-improved Cre (CD19-iCre) knock-in mouse with the T2A-iCre sequence inserted into the Cd19 locus, just before the stop codon. The CD19-iCre mice were compared with existing models, crossed with the Rosa26-EYFP reporter mice, and their recombination activity in B cells carrying different Cre alleles was assessed. CD19-iCre mice showed more effective Cre recombination in the early B cell developmental stages compared with the CD19-Cre mice. The efficiencies of the CD19-iCre and Mb1-iCre lines were similar; however, the B lineage-specific recombination was more stringent in the CD19-iCre line. Furthermore, the utility value of the CD19-iCre model was superior than that of the CD19-Cre mice regarding deletion efficiency in IL10-floxed mice. Thus, the CD19-iCre line is a valuable tool for highly efficient gene targeting specific to the B cell compartment.

Balancers are rearranged chromosomes used in Drosophila melanogaster to maintain deleterious mutations in stable populations, preserve sets of linked genetic elements and construct complex experimental stocks. Here, we assess the phenotypes associated with breakpoint-induced mutations on commonly used third chromosome balancers and show remarkably few deleterious effects. We demonstrate that a breakpoint in p53 causes loss of radiation-induced apoptosis and a breakpoint in Fucosyltransferase A causes loss of fucosylation in nervous and intestinal tissue—the latter study providing new markers for intestinal cell identity and challenging previous conclusions about the regulation of fucosylation. We also describe thousands of potentially harmful mutations shared among X or third chromosome balancers, or unique to specific balancers, including an Ankyrin 2 mutation present on most TM3 balancers, and reiterate the risks of using balancers as experimental controls. We used long-read sequencing to confirm or refine the positions of two inversions with breakpoints lying in repetitive sequences and provide evidence that one of the inversions, In(2L)Cy, arose by ectopic recombination between foldback transposon insertions and the other, In(3R)C, cleanly separates subtelomeric and telomeric sequences and moves the subtelomeric sequences to an internal chromosome position. In addition, our characterization of In(3R)C shows that balancers may be polymorphic for terminal deletions. Finally, we present evidence that extremely distal mutations on balancers can add to the stability of stocks whose purpose is to maintain homologous chromosomes carrying mutations in distal genes. Overall, these studies add to our understanding of the structure, diversity and effectiveness of balancer chromosomes.

Regulatory T (Treg) cell-specific deletion of a gene of interest is a procedure widely used to study mechanisms controlling Treg development, homeostasis and function. Accordingly, several transgenic mouse lines have been generated that bear the Cre recombinase under control of the Foxp3 promoter either as a random transgene insertion or knocked into the endogenous Foxp3 locus, with the Foxp3YFP-Cre strain of mice being one of the most widely used. In an attempt to generate Treg cells that lacked expression of the insulin receptor (Insr), we crossed Foxp3YFP-Cre mice with Insrfl/fl mice. Using a conventional two-band PCR genotyping method we found that offspring genotypes did not correspond to the expected Mendelian ratios. We therefore developed a quantitative PCR-based genotyping method to investigate possible ectopic recombination outside the Treg lineage. With this method we found that ~50% of the F1 -generation mice showed evidence of ectopic recombination and that ~10% of the F2 -generation mice had germline Cre recombination activity leading to a high frequency of offspring with global Insr deletion. Use of the quantitative PCR genotyping method enabled accurate selection of mice without ectopic recombination and only the desired Treg cell-specific Insr deletion. Our data highlight the need to use genotyping methods that allow for assessment of possible ectopic recombination driven by the Foxp3YFP-Cre allele, particularly when studying genes that are systemically expressed.

Aspergillus flavus is a fungus that produces aflatoxin B1, one of the most carcinogenic secondary metabolites. Understanding the regulation mechanism of aflatoxin biosynthesis in this fungus requires precise methods for genomic integration of mutant alleles. To avoid the disadvantage of DNA integration into the genome by non-homologous or ectopic recombination, we developed a novel strategy for site-specific integration of foreign DNA by using a carboxin-resistant sdh2R allele (His 249 Leu). Our results demonstrated that the transformants were generated with a high efficiency (>96%) of correct integration into the sdh2-lcus of the genome of A. flavus NRRL 3357. The advantage of this method is that introduction of the eGFP expression cassette into the sdh2-locus had little effect on fungal growth and virulence while also being rapid and efficient. This system will be a valuable tool for genetic manipulation in A. flavus. To the best of our knowledge, this is the first report on the efficient site-specific integration at the sdh2-locus in the genome of Aspergillus.

The number and distribution of meiotic crossovers (COs) are highly regulated, reflecting the requirement for COs during the first round of meiotic chromosome segregation. CO control includes CO assurance and CO interference, which promote at least one CO per chromosome bivalent and evenly-spaced COs, respectively. Previous studies revealed a role for the DNA damage response (DDR) clamp and the clamp loader in CO formation by promoting interfering COs and interhomolog recombination, and also by suppressing ectopic recombination. In this study, we use classical tetrad analysis of Saccharomyces cerevisiae to show that a mutant defective in RAD24, which encodes the DDR clamp loader (RAD17 in other organisms), displayed reduced CO frequencies on two shorter chromosomes (III and V), but not on a long chromosome (chromosome VII). The residual COs in the rad24 mutant do not show interference. In contrast to rad24, mutants defective in the ATR kinase homolog Mec1, including a mec1 null and a mec1 kinase-dead mutant, show slight or few defects in CO frequency. On the other hand, mec1 COs show defects in interference, similar to the rad24 mutant. Our results support a model in which the DDR clamp and clamp-loader proteins promote interfering COs by recruiting pro-CO Zip, Mer, and Msh proteins to recombination sites, while the Mec1 kinase regulates CO distribution by a distinct mechanism. Moreover, CO formation and its control are implemented in a chromosome-specific manner, which may reflect a role for chromosome size in regulation.

In order to understand how DNA sequences of transposable elements (TEs) evolve, extensive simulations were carried out. We first used our previous model, in which the copy number of TEs is mainly controlled by selection against ectopic recombination. It was found that along a simulation run, the shape of phylogeny changes quite much, from monophyletic trees to dimorphic trees with two clusters. Our results demonstrated that the change of the phase is usually slow from a monomorphic phase to a dimorphic phase, accompanied with a growth of an internal branch by accumulation of variation between two types. Then, the phase immediately changes back to a monomorphic phase when one group gets extinct. Under this condition, monomorphic and dimorphic phases arise repeatedly, and it is very difficult to maintain two or more different types of TEs for a long time. Then, how a new subfamily can evolve? To solve this, we developed a new model, in which ectopic recombination is restricted between two types under some condition, for example, accumulation of mutations between them. Under this model, because selection works on the copy number of each types separately, two types can be maintained for a long time. As expected, our simulations demonstrated that a new type arises and persists quite stably, and that it will be recognized as a new subfamily followed by further accumulation of mutations. It is indicated that how ectopic recombination is regulated in a genome is an important factor for the evolution of a new subfamily.

Chromosomal inversions are fundamental drivers of genome evolution. In the Drosophila genus, inversions have been widely characterized cytologically, and play an important role in local adaptation. Here, we characterize chromosomal inversions in the Drosophila nasuta species group using chromosome-level, reference-quality assemblies of seven species and subspecies in this clade. Reconstruction of ancestral karyotypes allowed us to infer the order in which the 22 identified inversions occurred along the phylogeny. We found a higher rate of inversions on the X chromosome, and heterogeneity in the rate of accumulation across the phylogeny. We molecularly characterize the breakpoints of six autosomal inversions, and found that repeated sequences are associated with inversion breakpoints in four of these inversions, suggesting that ectopic recombination is an important mechanism in generating inversion. Characterization of inversions in this species group provides a foundation for future population genetic and functional studies in this recently diverged species group.

The replicative nature and generally deleterious effects of transposable elements (TEs) give rise to an outstanding question about how TE copy number is stably contained in host populations. Classic theoretical analyses predict that, when the decline in fitness due to each additional TE insertion is greater than linear, or when there is synergistic epistasis, selection against TEs can result in a stable equilibrium of TE copy number. While several mechanisms are predicted to yield synergistic deleterious effects of TEs, we lack empirical investigations of the presence of such epistatic interactions. Purifying selection with synergistic epistasis generates repulsion linkage between deleterious alleles and, accordingly, an underdispersed distribution for the number of deleterious mutations among individuals. We investigated this population genetic signal in an African Drosophila melanogaster population and found evidence for synergistic epistasis among TE insertions, especially those expected to have large fitness impacts. Curiously, even though ectopic recombination has long been predicted to generate nonlinear fitness decline with increased TE copy number, TEs predicted to suffer higher rates of ectopic recombination are not more likely to be underdispersed. On the other hand, underdispersed TE families are more likely to show signatures of deleterious epigenetic effects and stronger ping-pong signals of piRNA amplification, a hypothesized source from which synergism of TE-mediated epigenetic effects arises. Our findings set the stage for investigating the importance of epistatic interactions in the evolutionary dynamics of TEs.

Background Despite progress in genomic analysis of spiders, their chromosome evolution is not satisfactorily understood. Most information on spider chromosomes concerns the most diversified clade, entelegyne araneomorphs. Other clades are far less studied. Our study focused on haplogyne araneomorphs, which are remarkable for their unusual sex chromosome systems and for the co-evolution of sex chromosomes and nucleolus organizer regions (NORs); some haplogynes exhibit holokinetic chromosomes. To trace the karyotype evolution of haplogynes on the family level, we analysed the number and morphology of chromosomes, sex chromosomes, NORs, and meiosis in pholcids, which are among the most diverse haplogyne families. The evolution of spider NORs is largely unknown. Results Our study is based on an extensive set of species representing all major pholcid clades. Pholcids exhibit a low 2n and predominance of biarmed chromosomes, which are typical haplogyne features. Sex chromosomes and NOR patterns of pholcids are diversified. We revealed six sex chromosome systems in pholcids (X0, XY, X1X20, X1X2X30, X1X2Y, and X1X2X3X4Y). The number of NOR loci ranges from one to nine. In some clades, NORs are also found on sex chromosomes. Conclusions The evolution of cytogenetic characters was largely derived from character mapping on a recently published molecular phylogeny of the family. Based on an extensive set of species and mapping of their characters, numerous conclusions regarding the karyotype evolution of pholcids and spiders can be drawn. Our results suggest frequent autosome–autosome and autosome–sex chromosome rearrangements during pholcid evolution. Such events have previously been attributed to the reproductive isolation of species. The peculiar X1X2Y system is probably ancestral for haplogynes. Chromosomes of the X1X2Y system differ considerably in their pattern of evolution. In some pholcid clades, the X1X2Y system has transformed into the X1X20 or XY systems, and subsequently into the X0 system. The X1X2X30 system of Smeringopus pallidus probably arose from the X1X20 system by an X chromosome fission. The X1X2X3X4Y system of Kambiwa probably evolved from the X1X2Y system by integration of a chromosome pair. Nucleolus organizer regions have frequently expanded on sex chromosomes, most probably by ectopic recombination. Our data suggest the involvement of sex chromosome-linked NORs in achiasmatic pairing.

Hybrid sterility maintains reproductive isolation between species by preventing them from exchanging genetic material. The genome sequences of the yeasts Saccharomyces paradoxus and Saccharomyces paradoxus have diverged by about 12%, yet they can mate and form viable diploid hybrids, which can reproduce asexually by mitosis. But these hybrids are sexually sterile: nearly all the gametes they produce by meiosis are inviable because they are aneuploid, lacking at least one essential chromosome. Chromosomes from the different species mis-segregate during hybrid meiosis due to anti-recombination; their sequences are too diverged to efficiently recombine and form crossovers. Previous methods to increase hybrid fertility by enhancing meiotic crossing over have had counteracting detrimental effects on asexual viability, due to increased mutagenesis and ectopic recombination. We were able to increase hybrid fertility 70-fold by suppressing two genes, MSH2 (YOL090W) and SGS1 (YMR190C), specifically during meiosis, whilst maintaining their mitotic expression. This confirms that anti-recombination is the principal cause of hybrid sterility. By overcoming this barrier, we were able to generate viable euploid hybrid gametes containing recombinant genomes from these two highly diverged parent species. We used this unique resource to map intrinsic Bateson-Dobzhansky-Muller incompatibilities contributing to hybrid gamete inviability. The ability to interbreed such diverged species will allow mapping of other intrinsic traits such as hybrid depression or vigour, extrinsic ecological adaptations that differentiate species, or commercially important characteristics.

Background Chromosomal inversions have been pervasive during the evolution of the genus Drosophila, but there is significant variation between lineages in the rate of rearrangement fixation. D. mojavensis, an ecological specialist adapted to a cactophilic niche under extreme desert conditions, is a chromosomally derived species with ten fixed inversions, five of them not present in any other species. Results In order to explore the causes of the rapid chromosomal evolution in D. mojavensis, we identified and characterized all breakpoints of seven inversions fixed in chromosome 2, the most dynamic one. One of the inversions presents unequivocal evidence for its generation by ectopic recombination between transposon copies and another two harbor inverted duplications of non-repetitive DNA at the two breakpoints and were likely generated by staggered single-strand breaks and repair by non-homologous end joining. Four out of 14 breakpoints lay in the intergenic region between preexisting duplicated genes, suggesting an adaptive advantage of separating previously tightly linked duplicates. Four out of 14 breakpoints are associated with transposed genes, suggesting these breakpoints are fragile regions. Finally two inversions contain novel genes at their breakpoints and another three show alterations of genes at breakpoints with potential adaptive significance. Conclusions D. mojavensis chromosomal inversions were generated by multiple mechanisms, an observation that does not provide support for increased mutation rate as explanation for rapid chromosomal evolution. On the other hand, we have found a number of gene alterations at the breakpoints with putative adaptive consequences that directly point to natural selection as the cause of D. mojavensis rapid chromosomal evolution.

Genome stability in radioresistant bacterium Deinococcus radiodurans depends on RecA, the main bacterial recombinase. Without RecA, gross genome rearrangements occur during repair of DNA double-strand breaks. Long repeated (insertion) sequences have been identified as hot spots for ectopic recombination leading to genome rearrangements, and single-strand annealing (SSA) postulated to be the most likely mechanism involved in this process. Here, we have sequenced five isolates of D. radiodurans recA mutant carrying gross genome rearrangements to precisely characterize the rearrangements and to elucidate the underlying repair mechanism. The detected rearrangements consisted of large deletions in chromosome II in all the sequenced recA isolates. The mechanism behind these deletions clearly differs from the classical SSA; it utilized short (4–11 bp) repeats as opposed to insertion sequences or other long repeats. Moreover, it worked over larger linear DNA distances from those previously tested. Our data are most compatible with alternative end-joining, a recombination mechanism that operates in eukaryotes, but is also found in Escherichia coli. Additionally, despite the recA isolates being preselected for different rearrangement patterns, all identified deletions were found to overlap in a 35 kb genomic region. We weigh the evidence for mechanistic vs. adaptive reasons for this phenomenon.

The mitochondria of flowering plants have large and complex genomes whose structure and segregation are modulated by recombination activities. Among unresolved questions is what are the pathways responsible for the late steps of homologous recombination: while the loss of mitochondrial recombination is not viable, a deficiency in RECG1-dependent branch migration has little impact on plant development. Here we present an additional pathway required for the processing of organellar recombination intermediates, the one depending on RADA . RADA is similar in structure and activity to bacterial RadA/Sms, and in vitro it binds to ssDNA and accelerates strand-exchange reactions initiated by RecA. RADA-deficient plants are severely impacted in their development and fertility, correlating with increased mtDNA ectopic recombination and replication of recombination-generated subgenomes. The radA mutation is epistatic to recG1 , indicating that RADA drives the main branch migration pathway of plant mitochondria. In contrast, the double mutation radA recA3 is lethal, revealing the importance of an alternative RECA3-dependent pathway. Interestingly, the radA developmental phenotypes could not be correlated with obvious defects in mitochondrial gene expression. Rather, it seems that it is the activation of genes that repress cell cycle progression that is partially the cause of the stunted growth of radA mutants. SYNOPSIS Arabidopsis RADA is a main branch migration activity in plant mitochondria, whose deficiency leads to mtDNA instability by recombination, and suppression of plant growth by the activation of repressors of cell cycle progression.

During prophase I of meiosis, homologous chromosomes pair, synapse and exchange their genetic material through reciprocal homologous recombination, a phenomenon essential for faithful chromosome segregation. Partial sequence identity between non-homologous and heterologous chromosomes can also lead to recombination (ectopic recombination), a highly deleterious process that rapidly compromises genome integrity. To avoid ectopic exchange, homology recognition must be extended from the narrow position of a crossover-competent double-strand break to the entire chromosome. Here, we review advances on chromosome behaviour during meiotic prophase I in higher plants, by integrating centromere- and telomere dynamics driven by cytoskeletal motor proteins, into the processes of homologue pairing, synapsis and recombination. Centromere–centromere associations and the gathering of telomeres at the onset of meiosis at opposite nuclear poles create a spatially organised and restricted nuclear state in which homologous DNA interactions are favoured but ectopic interactions also occur. The release and dispersion of centromeres from the nuclear periphery increases the motility of chromosome arms, allowing meiosis-specific movements that disrupt ectopic interactions. Subsequent expansion of interstitial synapsis from numerous homologous interactions further corrects ectopic interactions. Movement and organisation of chromosomes, thus, evolved to facilitate the pairing process, and can be modulated by distinct stages of chromatin associations at the nuclear envelope and their collective release., Summary: We review plant meiosis, including chromosome tethering at the nuclear periphery that, we propose, defines chromosome dynamics-facilitating DNA sequence-based pairing of homologues

Mitochondrial and plastid genomes in land plants exhibit some of the slowest rates of sequence evolution observed in any eukaryotic genome, suggesting an exceptional ability to prevent or correct mutations. However, the mechanisms responsible for this extreme fidelity remain unclear. We tested seven candidate genes involved in cytoplasmic DNA replication, recombination, and repair (POLIA, POLIB, MSH1, RECA3, UNG, FPG, and OGG1) for effects on mutation rates in the model angiosperm Arabidopsis thaliana by applying a highly accurate DNA sequencing technique (duplex sequencing) that can detect newly arisen mitochondrial and plastid mutations still at low heteroplasmic frequencies. We find that disrupting MSH1 (but not the other candidate genes) leads to massive increases in the frequency of point mutations and small indels and changes to the mutation spectrum in mitochondrial and plastid DNA. We also used droplet digital PCR to show transmission of de novo heteroplasmies across generations in msh1 mutants, confirming a contribution to heritable mutation rates. This dual-targeted gene is part of an enigmatic lineage within the mutS mismatch repair family that we find is also present outside of green plants in multiple eukaryotic groups (stramenopiles, alveolates, haptophytes, and cryptomonads), as well as certain bacteria and viruses. MSH1 has previously been shown to limit ectopic recombination in plant cytoplasmic genomes. Our results point to a broader role in recognition and correction of errors in plant mitochondrial and plastid DNA sequence, leading to greatly suppressed mutation rates perhaps via initiation of double-stranded breaks and repair pathways based on faithful homologous recombination.

When DNA double-strand breaks occur, four-stranded DNA structures called Holliday junctions (HJs) form during homologous recombination. Because HJs connect homologous DNA by a covalent link, resolution of HJ is crucial to terminate homologous recombination and segregate the pair of DNA molecules faithfully. We recently identified Monokaryotic Chloroplast1 (MOC1) as a plastid DNA HJ resolvase in algae and plants. Although Cruciform cutting endonuclease1 (CCE1) was identified as a mitochondrial DNA HJ resolvase in yeasts, homologs or other mitochondrial HJ resolvases have not been identified in other eukaryotes. Here, we demonstrate that MOC1 depletion in the green alga Chlamydomonas reinhardtii and the moss Physcomitrella patens induced ectopic recombination between short dispersed repeats in ptDNA. In addition, MOC1 depletion disorganized thylakoid membranes in plastids. In some land plant lineages, such as the moss P. patens, a liverwort and a fern, MOC1 dually targeted to plastids and mitochondria. Moreover, mitochondrial targeting of MOC1 was also predicted in charophyte algae and some land plant species. Besides causing instability of plastid DNA, MOC1 depletion in P. patens induced short dispersed repeat-mediated ectopic recombination in mitochondrial DNA and disorganized cristae in mitochondria. Similar phenotypes in plastids and mitochondria were previously observed in mutants of plastid-targeted (RECA2) and mitochondrion-targeted (RECA1) recombinases, respectively. These results suggest that MOC1 functions in the double-strand break repair in which a recombinase generates HJs and MOC1 resolves HJs in mitochondria of some lineages of algae and plants as well as in plastids in algae and plants.

Oriental theileriosis is an economically important tickborne disease of bovines, caused by some members of the Theileria orientalis complex. Currently, 11 distinct operational taxonomic units (OTUs), or genotypes, are recognized based on their major piroplasm surface protein (MPSP) gene sequences. Two of these genotypes (i.e., chitose and ikeda) are recognized as pathogenic in cattle, causing significant disease in countries of the Asia-Pacific region. However, the true extent of genetic variation and associated virulence/pathogenicity within this complex is unknown. Here, we undertook a proof-of-principle study of a small panel of genomic DNAs (n = 13) from blood samples originating from individual cattle known to harbor T. orientalis, in order to assess the performance of a targeted &rdquo, next-generation&rdquo, sequencing-informatic approach to identify genotypes. Five genotypes (chitose, ikeda, buffeli, type 4, and type 5) were defined, multiple genotypes were found within individual samples, with dominant and minor sequence types representing most genotypes. This study indicates that this sequencing-informatic workflow could be useful to assess the nature and extent of genetic variation within and among populations of T. orientalis on a large scale, and to potentially employ panels of distinct gene markers for expanded molecular epidemiological investigations of socioeconomically important protistan pathogens more generally.

In various eukaryotes, supernumerary B chromosomes (Bs) are an optional genomic component that affect their integrity and functioning. In the present study, the impact of Bs on the current changes in the genome of goatgrass, Aegilops speltoides, was addressed. Individual plants from contrasting populations with and without Bs were explored using fluorescence in situ hybridization. In parallel, abundances of the Ty1-copia, Ty3-gypsy, and LINE retrotransposons (TEs), and the species-specific Spelt1 tandem repeat (TR) in vegetative and generative spike tissues were estimated by real-time quantitative PCR. The results revealed: (i) ectopic associations between Bs and the regular A chromosomes, and (ii) cell-specific rearrangements of Bs in both mitosis and microgametogenesis. Further, the copy numbers of TEs and TR varied significantly between (iii) genotypes and (iv) different spike tissues in the same plant(s). Finally, (v) in plants with and without Bs from different populations, genomic abundances and/or copy number dynamics of TEs and TR were similar. These findings indicate that fluctuations in TE and TR copy numbers are associated with DNA damage and repair processes during cell proliferation and differentiation, and ectopic recombination is one of the mechanisms by which Bs play a role in genome changes.

Malaria remains one of the most prevalent and deadly infectious diseases of the developing world, causing approximately 228 million clinical cases and nearly half a million deaths annually. The disease is caused by protozoan parasites of the genus Plasmodium, and of the five species capable of infecting humans, infections with P. falciparum are the most severe. In addition to the parasites that infect people, there are hundreds of additional species that infect birds, reptiles, and other mammals, each exquisitely evolved to meet the specific challenges inherent to survival within their respective hosts. By comparing the unique strategies that each species has evolved, key insights into host-parasite interactions can be gained, including discoveries regarding the pathogenesis of human disease. Here, we describe the surprising observation that closely related parasites with different hosts have evolved remarkably different methods for repairing their genomes. This observation has important implications for the ability of parasites to maintain chronic infections and for the development of host immunity., The protozoan parasites that cause malaria infect a wide variety of vertebrate hosts, including birds, reptiles, and mammals, and the evolutionary pressures inherent to the host-parasite relationship have profoundly shaped the genomes of both host and parasite. Here, we report that these selective pressures have resulted in unexpected alterations to one of the most basic aspects of eukaryotic biology, the maintenance of genome integrity through DNA repair. Malaria parasites that infect humans continuously generate genetic diversity within their antigen-encoding gene families through frequent ectopic recombination between gene family members, a process that is a crucial feature of the persistence of malaria globally. The continuous generation of antigen diversity ensures that different parasite isolates are antigenically distinct, thus preventing extensive cross-reactive immunity and enabling parasites to maintain stable transmission within human populations. However, the molecular basis of the recombination between gene family members is not well understood. Through computational analyses of the antigen-encoding, multicopy gene families of different Plasmodium species, we report the unexpected observation that malaria parasites that infect rodents do not display the same degree of antigen diversity as observed in Plasmodium falciparum and appear to undergo significantly less ectopic recombination. Using comparative genomics, we also identify key molecular components of the diversification process, thus shedding new light on how malaria parasites balance the maintenance of genome integrity with the requirement for continuous genetic diversification.

Retrotransposable elements are widely distributed and diverse in eukaryotes. Their copy number increases through reverse-transcription-mediated propagation, while they can be lost through recombinational processes, generating genomic rearrangements. We previously identified extensive structurally uniform retrotransposon groups in which no member contains the gag, pol, or env internal domains. Because of the lack of protein-coding capacity, these groups are non-autonomous in replication, even if transcriptionally active. The Cassandra element belongs to the non-autonomous group called terminal-repeat retrotransposons in miniature (TRIM). It carries 5S RNA sequences with conserved RNA polymerase (pol) III promoters and terminators in its long terminal repeats (LTRs). Here, we identified multiple extended tandem arrays of Cassandra retrotransposons within different plant species, including ferns. At least 12 copies of repeated LTRs (as the tandem unit) and internal domain (as a spacer), giving a pattern that resembles the cellular 5S rRNA genes, were identified. A cytogenetic analysis revealed the specific chromosomal pattern of the Cassandra retrotransposon with prominent clustering at and around 5S rDNA loci. The secondary structure of the Cassandra retroelement RNA is predicted to form super-loops, in which the two LTRs are complementary to each other and can initiate local recombination, leading to the tandem arrays of Cassandra elements. The array structures are conserved for Cassandra retroelements of different species. We speculate that recombination events similar to those of 5S rRNA genes may explain the wide variation in Cassandra copy number. Likewise, the organization of 5S rRNA gene sequences is very variable in flowering plants, part of what is taken for 5S gene copy variation may be variation in Cassandra number. The role of the Cassandra 5S sequences remains to be established.

Retrotransposon Hot Spot (RHS) is the most abundant gene family in Trypanosoma cruzi, with unknown function in this parasite. The aim of this work was to shed light on the organization and expression of RHS in T. cruzi. The diversity of the RHS protein family in T. cruzi was demonstrated by phylogenetic and recombination analyses. Transcribed sequences carrying the RHS domain were classified into ten distinct groups of monophyletic origin. We identified numerous recombination events among the RHS and traced the origins of the donors and target sequences. The transcribed RHS genes have a mosaic structure that may contain fragments of different RHS inserted in the target sequence. About 30% of RHS sequences are located in the subtelomere, a region very susceptible to recombination. The evolution of the RHS family has been marked by many events, including gene duplication by unequal mitotic crossing-over, homologous, as well as ectopic recombination, and gene conversion. The expression of RHS was analyzed by immunofluorescence and immunoblotting using anti-RHS antibodies. RHS proteins are evenly distributed in the nuclear region of T. cruzi replicative forms (amastigote and epimastigote), suggesting that they could be involved in the control of the chromatin structure and gene expression, as has been proposed for T. brucei.

In Drosophila, large variations in rearrangement rate have been reported among different lineages and among Muller’s elements. Nevertheless, the mechanisms that are involved in the generation of inversions, their increase in frequency, as well as their impact on the genome are not completely understood. This is in part due to the lack of comparative studies on species distantly related to Drosophila melanogaster. Therefore, we sequenced and assembled the genomes of two species of the virilis phylad (Drosophila novamexicana [15010-1031.00] and Drosophila americana [SF12]), which are diverging from D. melanogaster for more than 40 Myr. Based on these data, we identified the precise location of six novel inversion breakpoints. A molecular characterization provided clear evidence that DAIBAM (a miniature inverted–repeat transposable element) was involved in the generation of eight out of the nine inversions identified. In contrast to what has been previously reported for D. melanogaster and close relatives, ectopic recombination is thus the prevalent mechanism of generating inversions in species of the virilis phylad. Using pool-sequencing data for three populations of D. americana, we also show that common polymorphic inversions create a high degree of genetic differentiation between populations for chromosomes X, 4, and 5 over large physical distances. We did not find statistically significant differences in expression levels between D. americana (SF12) and D. novamexicana (15010-1031.00) strains for the three genes surveyed (CG9588, Fig 4, and fab1) flanking three inversion breakpoints. This article is a result of the project Norte-01-0145-FEDER-000008—Porto Neurosciences and Neurologic Disease Research Initiative at I3S, supported by Norte Portugal Regional Operational Programme (NORTE 2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (FEDER). N.P. and M.R. are funded by the Emmy Noether Programme of the Deutsche Forschungsgemeinschaft (Grant Number: PO 1648/3-1 to N.P.). We would like to thank the Transcriptome Analysis Lab (TAL) (University Medical Center Göttingen, UMG) in Göttingen for the Illumina sequencing.

Transposable elements (TEs) are genomic parasites that impose fitness costs on their hosts by producing deleterious mutations and disrupting gametogenesis. Host genomes avoid these costs by regulating TE activity, particularly in germline cells where new insertions are heritable and TEs are exceptionally active. However, the capacity of different TE-associated fitness costs to select for repression in the host, and the role of selection in the evolution of TE regulation more generally, remain controversial. In this study, we use forward, individual-based simulations to examine the evolution of small-RNA-mediated TE regulation, a conserved mechanism for TE repression that is employed by both prokaryotes and eukaryotes. To design and parameterize a biologically realistic model, we drew on an extensive survey of empirical studies of the transposition and regulation ofP-element DNA transposons inDrosophila melanogaster. We observed that even under conservative assumptions, where small-RNA-mediated regulation reduces transposition only, repression evolves rapidly and adaptively after the genome is invaded by a new TE. We further show that the spread of repressor alleles is greatly enhanced by two additional TE-imposed fitness costs: dysgenic sterility and ectopic recombination. Finally, we demonstrate that the mutation rate to repression (i.e., the size of the mutational target) is a critical parameter that influences both the evolutionary trajectory of host repression and the associated proliferation of TEs after invasion. Our findings suggest that adaptive evolution of TE regulation may be stronger and more prevalent than previously appreciated, and provide a framework for evaluating empirical data.

Heterochromatin mainly comprises repeated DNA sequences that are prone to ectopic recombination. In Drosophila cells, ‘safe’ repair of heterochromatic double-strand breaks by homologous recombination relies on the relocalization of repair sites to the nuclear periphery before strand invasion. The mechanisms responsible for this movement were unknown. Here we show that relocalization occurs by directed motion along nuclear actin filaments assembled at repair sites by the Arp2/3 complex. Relocalization requires nuclear myosins associated with the heterochromatin repair complex Smc5/6 and the myosin activator Unc45, which is recruited to repair sites by Smc5/6. ARP2/3, actin nucleation and myosins also relocalize heterochromatic double-strand breaks in mouse cells. Defects in this pathway result in impaired heterochromatin repair and chromosome rearrangements. These findings identify de novo nuclear actin filaments and myosins as effectors of chromatin dynamics for heterochromatin repair and stability in multicellular eukaryotes. Relocalization of heterochromatic double-strand breaks to the nuclear periphery in Drosophila cells occurs via directed motions driven by nuclear actin filaments and myosins activated by the Smc5/6 complex.

Transposable elements (TE) can constitute a large fraction of plant genomes, yet our understanding of their evolution and fitness effect is still limited. Here we tested several models of evolution that make specific predictions about differences in TE abundance between selfing and outcrossing taxa, and between small and large populations. We estimated TE abundance in multiple populations of North American Arabidopsis lyrata differing in mating system and long-term size, using transposon insertion display on several TE families. Selfing populations had higher TE copy numbers per individual and higher TE allele frequencies, supporting models which assume that selection against TEs acts predominantly against heterozygotes via the process of ectopic recombination. In outcrossing populations differing in long-term size, the data supported neither a model of density-regulated transposition nor a model of direct deleterious effect. Instead, the population structure of TEs revealed that outcrossing populations tended to split into western and eastern groups - as previously detected using microsatellite markers - whereas selfing populations from west and east were less differentiated. This, too, agrees with the model of ectopic recombination. Overall, our results suggest that TE elements are nearly neutral except for their deleterious potential to disturb meiosis in the heterozygous state.

The genome of Plasmodium falciparum, the causative agent of malaria in Africa, has been extensively studied since it was first fully sequenced in 2002. However, many open questions remain, including understanding the chromosomal context of molecular evolutionary changes (e.g., relationship between chromosome map and phylogenetic conservation, patterns of gene duplication, and patterns of selection). Here, we present PhyloChromoMap, a method that generates a phylogenomic map of chromosomes from a custom-built bioinformatics pipeline. Using P. falciparum 3D7 as a model, we analyze 2,116 genes with homologs in up to 941 diverse eukaryotic, bacterial and archaeal lineages. We estimate the level of conservation along chromosomes based on conservation across clades, and identify “young” regions (i.e., those with recent or fast evolving genes) that are enriched in subtelomeric regions as compared with internal regions. We also demonstrate that patterns of molecular evolution for paralogous genes differ significantly depending on their location as younger paralogs tend to be found in subtelomeric regions whereas older paralogs are enriched in internal regions. Combining these observations with analyses of synteny, we demonstrate that subtelomeric regions are actively shuffled among chromosome ends, which is consistent with the hypothesis that these regions are prone to ectopic recombination. We also assess patterns of selection by comparing dN/dS ratios of gene family members in subtelomeric versus internal regions, and we include the important antigenic gene family var. These analyses illustrate the highly dynamic nature of the karyotype of P. falciparum, and provide a method for exploring genome dynamics in other lineages.