Search

Background Recent evidence has linked the gut microbiome to host behavior via the gut–brain axis [1–3]; however, the underlying mechanisms remain unexplored. Here, we determined the links between host genetics, the gut microbiome and memory using the genetically defined Collaborative Cross (CC) mouse cohort, complemented with microbiome and metabolomic analyses in conventional and germ-free (GF) mice. Results A genome-wide association analysis (GWAS) identified 715 of 76,080 single-nucleotide polymorphisms (SNPs) that were significantly associated with short-term memory using the passive avoidance model. The identified SNPs were enriched in genes known to be involved in learning and memory functions. By 16S rRNA gene sequencing of the gut microbial community in the same CC cohort, we identified specific microorganisms that were significantly correlated with longer latencies in our retention test, including a positive correlation with Lactobacillus. Inoculation of GF mice with individual species of Lactobacillus (L. reuteri F275, L. plantarum BDGP2 or L. brevis BDGP6) resulted in significantly improved memory compared to uninoculated or E. coli DH10B inoculated controls. Untargeted metabolomics analysis revealed significantly higher levels of several metabolites, including lactate, in the stools of Lactobacillus-colonized mice, when compared to GF control mice. Moreover, we demonstrate that dietary lactate treatment alone boosted memory in conventional mice. Mechanistically, we show that both inoculation with Lactobacillus or lactate treatment significantly increased the levels of the neurotransmitter, gamma-aminobutyric acid (GABA), in the hippocampus of the mice. Conclusion Together, this study provides new evidence for a link between Lactobacillus and memory and our results open possible new avenues for treating memory impairment disorders using specific gut microbial inoculants and/or metabolites.

The gut microbiome produces vitamins, nutrients, and neurotransmitters, and helps to modulate the host immune system—and also plays a major role in the metabolism of many exogenous compounds, including drugs and chemical toxicants. However, the extent to which specific microbial species or communities modulate hazard upon exposure to chemicals remains largely opaque. Focusing on the effects of collateral dietary exposure to the widely used herbicide atrazine, we applied integrated omics and phenotypic screening to assess the role of the gut microbiome in modulating host resilience in Drosophila melanogaster. Transcriptional and metabolic responses to these compounds are sex-specific and depend strongly on the presence of the commensal microbiome. Sequencing the genomes of all abundant microbes in the fly gut revealed an enzymatic pathway responsible for atrazine detoxification unique to Acetobacter tropicalis. We find that Acetobacter tropicalis alone, in gnotobiotic animals, is sufficient to rescue increased atrazine toxicity to wild-type, conventionally reared levels. This work points toward the derivation of biotic strategies to improve host resilience to environmental chemical exposures, and illustrates the power of integrative omics to identify pathways responsible for adverse health outcomes., Brown et al. apply integrated omics and phenotypic screening to assess the role of the gut microbiome in modulating host resilience in Drosophila melanogaster. They find that Acetobacter tropicalis in gnotobiotic animals, is sufficient to rescue increased atrazine toxicity, which could pave the way for biotic strategies to improve host resilience to environmental chemical exposure.

Advances in engineering biology have expanded the list of renewable compounds that can be produced at scale via biological routes from plant biomass. In most cases, these chemical products have not been evaluated for effects on biological systems, defined in the present study as bioactivity, that may be relevant to their manufacture. For sustainable chemical and fuel production, the industry needs to transition from fossil to renewable carbon sources, resulting in unprecedented expansion in the production and environmental distribution of chemicals used in biomanufacturing. Further, although some chemicals have been assessed for mammalian toxicity, environmental and agricultural hazards are largely unknown. We assessed 6 compounds that are representative of the emerging biofuel and bioproduct manufacturing process for their effect on model plants (Arabidopsis thaliana, Sorghum bicolor) and show that several alter plant seedling physiology at submillimolar concentrations. However, these responses change in the presence of individual bacterial species from the A. thaliana root microbiome. We identified 2 individual microbes that change the effect of chemical treatment on root architecture and a pooled microbial community with different effects relative to its constituents individually. The present study indicates that screening industrial chemicals for bioactivity on model organisms in the presence of their microbiomes is important for biologically and ecologically relevant risk analyses. Environ Toxicol Chem 2019;38:1911–1922. © 2019 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals, Inc. on behalf of SETAC.

INTRODUCTION: To faithfully distribute genetic material to daughter cells during cell division, spindle fibers must couple to DNA by means of a structure called the kinetochore, which assembles at each chromosome’s centromere. Human centromeres are located within large arrays of tandemly repeated DNA sequences known as alpha satellite (αSat), which often span millions of base pairs on each chromosome. Arrays of αSat are frequently surrounded by other types of tandem satellite repeats, which have poorly understood functions, along with nonrepetitive sequences, including transcribed genes. Previous genome sequencing efforts have been unable to generate complete assemblies of satellite-rich regions because of their scale and repetitive nature, limiting the ability to study their organization, variation, and function. RATIONALE: Pericentromeric and centromeric (peri/centromeric) satellite DNA sequences have remained almost entirely missing from the assembled human reference genome for the past 20 years. Using a complete, telomere-to-telomere (T2T) assembly of a human genome, we developed and deployed tailored computational approaches to reveal the organization and evolutionary patterns of these satellite arrays at both large and small length scales. We also performed experiments to map precisely which αSat repeats interact with kinetochore proteins. Last, we compared peri/centromeric regions among multiple individuals to understand how these sequences vary across diverse genetic backgrounds. RESULTS: Satellite repeats constitute 6.2% of the T2T-CHM13 genome assembly, with αSat representing the single largest component (2.8% of the genome). By studying the sequence relationships of αSat repeats in detail across each centromere, we found genome-wide evidence that human centromeres evolve through “layered expansions.” Specifically, distinct repetitive variants arise within each centromeric region and expand through mechanisms that resemble successive tandem duplications, whereas older flanking sequences shrink and diverge over time. We also revealed that the most recently expanded repeats within each αSat array are more likely to interact with the inner kinetochore protein Centromere Protein A (CENP-A), which coincides with regions of reduced CpG methylation. This suggests a strong relationship between local satellite repeat expansion, kinetochore positioning, and DNA hypomethylation. Furthermore, we uncovered large and unexpected structural rearrangements that affect multiple satellite repeat types, including active centromeric αSat arrays. Last, by comparing sequence information from nearly 1600 individuals’ Xchromosomes, we observed that individuals with recent African ancestry possess the greatest genetic diversity in the region surrounding the centromere, which sometimes contains a predominantly African αSat sequence variant. CONCLUSION: The genetic and epigenetic properties of centromeres are closely interwoven through evolution. These findings raise important questions about the specific molecular mechanisms responsible for the relationship between inner kinetochore proteins, DNA hypomethylation, and layered αSat expansions. Even more questions remain about the function and evolution of non-αSat repeats. To begin answering these questions, we have produced a comprehensive encyclopedia of peri/centromeric sequences in a human genome, and we demonstrated how these regions can be studied with modern genomic tools. Our work also illuminates the rich genetic variation hidden within these formerly missing regions of the genome, which may contribute to health and disease. This unexplored variation underlines the need for more T2T human genome assemblies from genetically diverse individuals.

Exposure to ionizing radiation is marked by NASA as a major health hazard for deep space exploration missions. Ionizing radiation sensitivity is determined by both genomic and environmental factors. Understanding their contributions is crucial for designing experiments in model organisms, selecting therapeutic irradiation regimes for cancer patients and evaluating the risk of deep space radiation exposure in astronauts. We identified single nucleotide polymorphisms in 15 strains of mice associated with spontaneous and ionizing radiation-induced ex vivo DNA damage. We mapped them to pathways related to carcinogenesis, nervous system damage and immune activation, including some located within the protein coding regions of genes that were predicted to interfere with protein functions. We anticipate that the identification of genes and pathways associated with DNA damage in response to ionizing radiation will improve the selection of mouse models for ionizing radiation research, inform functional validation studies and identify targets for space travel countermeasures.

Variations in oral bacterial communities have been linked to oral cancer suggesting that the oral microbiome is an etiological factor that can influence oral cancer development. The 4-nitroquinoline 1-oxide (4-NQO)-induced murine oral and esophageal cancer model is frequently used to assess the effects of preventive and/or therapeutic agents. We used this model to assess the impact of the microbiome on tumorigenesis using axenic (germ-free) and conventionally housed mice. Increased toxicity was observed in germ-free mice, however, no difference in tumor incidence, multiplicity, and size was observed. Transcriptional profiling of liver tissue from germ-free and conventionally housed mice identified 254 differentially expressed genes including ten cytochrome p450 enzymes, the largest family of phase-1 drug metabolizing enzymes in the liver. Gene ontology revealed that differentially expressed genes were enriched for liver steatosis, inflammation, and oxidative stress in livers of germ-free mice. Our observations emphasize the importance of the microbiome in mediating chemical toxicity at least in part by altering host gene expression. Studies on the role of the microbiome in chemical-induced cancer using germ-free animal models should consider the potential difference in dose due to the microbiome-mediated changes in host metabolizing capacity, which might influence the ability to draw conclusions especially for tumor induction models that are dose dependent.

ObjectiveThe Collaborative Cross (CC) is a mouse population model with diverse and reproducible genetic backgrounds used to identify novel disease models and genes that contribute to human disease. Since spontaneous tumour susceptibility in CC mice remains unexplored, we assessed tumour incidence and spectrum.DesignWe monitored 293 mice from 18 CC strains for tumour development. Genetic association analysis and RNA sequencing were used to identify susceptibility loci and candidate genes. We analysed genomes of patients with gastric cancer to evaluate the relevance of genes identified in the CC mouse model and measured the expression levels of ISG15 by immunohistochemical staining using a gastric adenocarcinoma tissue microarray. Association of gene expression with overall survival (OS) was assessed by Kaplan-Meier analysis.ResultsCC mice displayed a wide range in the incidence and types of spontaneous tumours. More than 40% of CC036 mice developed gastric tumours within 1 year. Genetic association analysis identified Nfκb1 as a candidate susceptibility gene, while RNA sequencing analysis of non-tumour gastric tissues from CC036 mice showed significantly higher expression of inflammatory response genes. In human gastric cancers, the majority of human orthologues of the 166 mouse genes were preferentially altered by amplification or deletion and were significantly associated with OS. Higher expression of the CC036 inflammatory response gene signature is associated with poor OS. Finally, ISG15 protein is elevated in gastric adenocarcinomas and correlated with shortened patient OS.ConclusionsCC strains exhibit tremendous variation in tumour susceptibility, and we present CC036 as a spontaneous laboratory mouse model for studying human gastric tumourigenesis.

Agrobacterium sp. strain 33MFTa1.1 was isolated for functional host-microbe interaction studies from the Thlaspi arvense root-associated microbiome. The complete genome is comprised of a circular chromosome of 2,771,937 bp, a linear chromosome of 2,068,443 bp, and a plasmid of 496,948 bp, with G+C contents of 59%, 59%, and 58%, respectively.

Despite critical roles in chromosome segregation and disease, the repetitive structure and vast size of centromeres and their surrounding heterochromatic regions impede studies of genomic variation. We report here large-scale haplotypes (cenhaps) in humans that span the centromere-proximal regions of all metacentric chromosomes, including the arrays of highly repeated α-satellites on which centromeres form.Cenhapsreveal surprisingly deep diversity, including entire introgressed Neanderthal centromeres and equally ancient lineages among Africans. These centromere-spanning haplotypes contain variants, including large differences in α-satellite DNA content, which may influence the fidelity and bias of chromosome transmission. The discovery ofcenhapscreates new opportunities to investigate their contribution to phenotypic variation, especially in meiosis and mitosis, as well as to more incisively model the unexpectedly rich evolution of these challenging genomic regions.One Sentence SummaryGenomic polymorphism across centromeric regions of humans is organized into large-scale haplotypes with great diversity, including entire Neanderthal centromeres.

Exposure to thirdhand smoke (THS) is a recently described health concern that arises in many indoor environments. However, the carcinogenic potential of THS, a critical consideration in risk assessment, remains untested. Here we investigated the effects of short-term early exposure to THS on lung carcinogenesis in A/J mice. Forty weeks after THS exposure from 4 to 7 weeks of age, the mice had increased incidence of lung adenocarcinoma, tumor size and, multiplicity, compared with controls. In vitro studies using cultured human lung cancer cells showed that THS exposure induced DNA double-strand breaks and increased cell proliferation and colony formation. RNA sequencing analysis revealed that THS exposure induced endoplasmic reticulum stress and activated p53 signaling. Activation of the p53 pathway was confirmed by an increase in its targets p21 and BAX. These data indicate that early exposure to THS is associated with increased lung cancer risk.

© 2017 Elsevier Inc. Eukaryotic genomes are broadly divided between gene-rich euchromatin and the highly repetitive heterochromatin domain, which is enriched for proteins critical for genome stability and transcriptional silencing. This study shows that Drosophila KDM4A (dKDM4A), previously characterized as a euchromatic histone H3 K36 demethylase and transcriptional regulator, predominantly localizes to heterochromatin and regulates heterochromatin position-effect variegation (PEV), organization of repetitive DNAs, and DNA repair. We demonstrate that dKDM4A demethylase activity is dispensable for PEV. In contrast, dKDM4A enzymatic activity is required to relocate heterochromatic double-strand breaks outside the domain, as well as for organismal survival when DNA repair is compromised. Finally, DNA damage triggers dKDM4A-dependent changes in the levels of H3K56me3, suggesting that dKDM4A demethylates this heterochromatic mark to facilitate repair. We conclude that dKDM4A, in addition to its previously characterized role in euchromatin, utilizes both enzymatic and structural mechanisms to regulate heterochromatin organization and functions. Colmenares et al. discover that Drosophila KDM4A, previously characterized as a euchromatic histone H3K36 demethylase and transcriptional regulator, is recruited to heterochromatin to contribute non-enzymatically to position-effect variegation, a hallmark of heterochromatin integrity. Conversely, dKDM4A catalytic activity is vital to heterochromatin DNA repair and is associated with demethylation of heterochromatic H3K56me3.

The concept that a breast cancer patient's menstrual stage at the time of tumor surgery influences risk of metastases remains controversial. The scarcity of comprehensive molecular studies of menstrual stage-dependent fluctuations in the breast provides little insight in this observation. To gain a deeper understanding of the biological changes in mammary tissue and blood during the menstrual cycle and to determine the influence of environmental exposures, such as low-dose ionizing radiation (LDIR), we used the mouse to characterize estrous-cycle variations in mammary gene transcripts by RNA-sequencing, peripheral white blood cell (WBC) counts and plasma cytokine levels. We identified an estrous-variable and hormone-dependent gene cluster enriched for Type-1 interferon genes. Cox regression identified a 117-gene signature of interferon-associated genes, which correlated with lower frequencies of metastasis in breast cancer patients. LDIR (10cGy) exposure had no detectable effect on mammary transcripts. However, peripheral WBC counts varied across the estrous cycle and LDIR exposure reduced lymphocyte counts and cytokine levels in tumor-susceptible mice. Our finding of variations in mammary Type-1 interferon and immune functions across the estrous cycle provides a mechanism by which timing of breast tumor surgery during the menstrual cycle may have clinical relevance to a patient's risk for distant metastases.

Although the gut microbiome plays important roles in host physiology, health and disease1, we lack understanding of the complex interplay between host genetics and early life environment on the microbial and metabolic composition of the gut. We used the genetically diverse Collaborative Cross mouse system2 to discover that early life history impacts the microbiome composition, whereas dietary changes have only a moderate effect. By contrast, the gut metabolome was shaped mostly by diet, with specific non-dietary metabolites explained by microbial metabolism. Quantitative trait analysis identified mouse genetic trait loci (QTL) that impact the abundances of specific microbes. Human orthologues of genes in the mouse QTL are implicated in gastrointestinal cancer. Additionally, genes located in mouse QTL for Lactobacillales abundance are implicated in arthritis, rheumatic disease and diabetes. Furthermore, Lactobacillales abundance was predictive of higher host T-helper cell counts, suggesting an important link between Lactobacillales and host adaptive immunity.

This report of independent genome sequences of two natural populations of Drosophila melanogaster (37 from North America and 6 from Africa) provides unique insight into forces shaping genomic polymorphism and divergence. Evidence of interactions between natural selection and genetic linkage is abundant not only in centromere- and telomere-proximal regions, but also throughout the euchromatic arms. Linkage disequilibrium, which decays within 1 kbp, exhibits a strong bias toward coupling of the more frequent alleles and provides a high-resolution map of recombination rate. The juxtaposition of population genetics statistics in small genomic windows with gene structures and chromatin states yields a rich, high-resolution annotation, including the following: (1) 5′- and 3′-UTRs are enriched for regions of reduced polymorphism relative to lineage-specific divergence; (2) exons overlap with windows of excess relative polymorphism; (3) epigenetic marks associated with active transcription initiation sites overlap with regions of reduced relative polymorphism and relatively reduced estimates of the rate of recombination; (4) the rate of adaptive nonsynonymous fixation increases with the rate of crossing over per base pair; and (5) both duplications and deletions are enriched near origins of replication and their density correlates negatively with the rate of crossing over. Available demographic models of X and autosome descent cannot account for the increased divergence on the X and loss of diversity associated with the out-of-Africa migration. Comparison of the variation among these genomes to variation among genomes from D. simulans suggests that many targets of directional selection are shared between these species.

Evidence has emerged that suggests a link between motor deficits, obesity and many neurological disorders. However, the contributing genetic risk factors are poorly understood. Here we used the Collaborative Cross (CC), a large panel of newly inbred mice that captures 90% of the known variation among laboratory mice, to identify the genetic loci controlling rotarod performance and its relationship with body weight in a cohort of 365 mice across 16 CC strains. Body weight and rotarod performance varied widely across CC strains and were significantly negatively correlated. Genetic linkage analysis identified 14 loci that were associated with body weight. However, 45 loci affected rotarod performance, seven of which were also associated with body weight, suggesting a strong link at the genetic level. Lastly, we show that genes identified in this study overlap significantly with those related to neurological disorders and obesity found in human GWA studies. In conclusion, our results provide a genetic framework for studies of the connection between body weight, the central nervous system and behavior.

Genome function is dynamically regulated in part by chromatin, which consists of the histones, non-histone proteins and RNA molecules that package DNA. Studies in Caenorhabditis elegans and Drosophila melanogaster have contributed substantially to our understanding of molecular mechanisms of genome function in humans, and have revealed conservation of chromatin components and mechanisms. Nevertheless, the three organisms have markedly different genome sizes, chromosome architecture and gene organization. On human and fly chromosomes, for example, pericentric heterochromatin flanks single centromeres, whereas worm chromosomes have dispersed heterochromatin-like regions enriched in the distal chromosomal ‘arms’, and centromeres distributed along their lengths. To systematically investigate chromatin organization and associated gene regulation across species, we generated and analysed a large collection of genome-wide chromatin data sets from cell lines and developmental stages in worm, fly and human. Here we present over 800 new data sets from our ENCODE and modENCODE consortia, bringing the total to over 1,400. Comparison of combinatorial patterns of histone modifications, nuclear lamina-associated domains, organization of large-scale topological domains, chromatin environment at promoters and enhancers, nucleosome positioning, and DNA replication patterns reveals many conserved features of chromatin organization among the three organisms. We also find notable differences in the composition and locations of repressive chromatin. These data sets and analyses provide a rich resource for comparative and species-specific investigations of chromatin composition, organization and function., National Science Foundation (U.S.) (1122374)

An estimated 80% of genomic DNA in eukaryotes is packaged as nucleosomes, which, together with the remaining interstitial linker regions, generate higher order chromatin structures [1]. Nucleosome sequences isolated from diverse organisms exhibit ∼10 bp periodic variations in AA, TT and GC dinucleotide frequencies. These sequence elements generate intrinsically curved DNA and help establish the histone-DNA interface. We investigated an important unanswered question concerning the interplay between chromatin organization and genome evolution: do the DNA sequence preferences inherent to the highly conserved histone core exert detectable natural selection on genomic divergence and polymorphism? To address this hypothesis, we isolated nucleosomal DNA sequences from Drosophila melanogaster embryos and examined the underlying genomic variation within and between species. We found that divergence along the D. melanogaster lineage is periodic across nucleosome regions with base changes following preferred nucleotides, providing new evidence for systematic evolutionary forces in the generation and maintenance of nucleosome-associated dinucleotide periodicities. Further, Single Nucleotide Polymorphism (SNP) frequency spectra show striking periodicities across nucleosomal regions, paralleling divergence patterns. Preferred alleles occur at higher frequencies in natural populations, consistent with a central role for natural selection. These patterns are stronger for nucleosomes in introns than in intergenic regions, suggesting selection is stronger in transcribed regions where nucleosomes undergo more displacement, remodeling and functional modification. In addition, we observe a large-scale (∼180 bp) periodic enrichment of AA/TT dinucleotides associated with nucleosome occupancy, while GC dinucleotide frequency peaks in linker regions. Divergence and polymorphism data also support a role for natural selection in the generation and maintenance of these super-nucleosomal patterns. Our results demonstrate that nucleosome-associated sequence periodicities are under selective pressure, implying that structural interactions between nucleosomes and DNA sequence shape sequence evolution, particularly in introns., Author Summary In eukaryotic cells, the majority of DNA is packaged in nucleosomes comprised of ∼147 bp of DNA wound tightly around the highly conserved histone octamer. Nucleosomal DNA from diverse organisms shows an anti-correlated ∼10 bp periodicity of AT-rich and GC-rich dinucleotides. These sequence features influence DNA bending and shape, facilitating structural interactions. We asked whether natural selection mediated through the periodic sequence preferences of nucleosomes shapes the evolution of non-protein-coding regions of D. melanogaster by examining the inter- and intra-species genomic variation relative to these fundamental chromatin building blocks. The sequence changes across nucleosome-bound regions on the melanogaster lineage mirror the observed nucleosome dinucleotide periodicities. Importantly, we show that the frequencies of polymorphisms in natural populations vary across these regions, paralleling divergence, with higher frequencies of preferred alleles. These patterns are most evident for intronic regions and indicate that non-protein coding regions are evolving toward sequences that facilitate the canonical association with the histone core. This result is consistent with the hypothesis that interactions between DNA and the core have systematic impacts on function that are subject to natural selection and are not solely due to mutational bias. These ubiquitous interactions with the histone core partially account for the evolutionary constraint observed in unannotated genomic regions, and may drive broad changes in base composition.

Meiotic recombination generates crossovers, which physically link homologous chromosomes, and direct their segregation to opposite poles. Despite the benefits, crossovers are restricted both in number as well as spatially. Crossovers in the nematode C. elegans occur most frequently within the distal regions (or “arms”) of the chromosomes, which comprise ~50% of the genome, but ~90% of the genetic map length. Meiotic recombination initiates with the formation of programmed double strand breaks (DSBs) catalyzed by the SPO-11 enzyme. It has been unclear whether the crossover bias in C. elegans is a consequence of DSB distribution, or instead reflects a mechanism that biases the downstream repair outcome of DSBs. In order to address this, we generated a genome-wide map of meiotic DSBs by mapping the distribution of the sole C. elegans RecA-like recombinase, RAD-51, by ChIP-seq. We have found that the RAD-51 binding pattern is strikingly similar to the known recombination pattern in C. elegans. The majority of RAD-51 binding sites occur on chromosome arms. Comparing RAD-51 distribution relative to functional genomic elements revealed that RAD-51 is particularly enriched in intronic DNA while coding sequences are underrepresented in RAD-51 occupied regions. We find that monomethylation of histone H3 lysine 36 (H3K36me1), which has not previously been associated with meiotic recombination, is highly correlated with break site preference. H3K36me1 is particularly enriched over active genes with high intron content, suggesting that the specific chromatin environment associated with these genes favors DSB formation in C. elegans. We have also identified a sequence motif that is highly overrepresented in DSB sites. Although the DSB pattern is highly correlated with crossover frequency, there are interesting exceptions. In particular, the sub-telomeric regions are completely devoid of crossovers but are active DSB sites. We have also identified a DNA sequence motif that appears to be enriched within DSB sites suggesting that sequence specific factors may have a role in DSB specification in C. elegans. In a parallel study, we have developed an Illumina Golden Gate SNP array to rapidly measure genome-wide recombination rates in C. elegans. Using this array we measured meiotic recombination rates in met-2(n4256) mutant animals. MET-2 is one of the two histone methyltransferases responsible for H3K9 methylation, a histone mark associated with heterochromatin. We find that the genetic map is significantly altered in met-2(n4256) mutant animals. In particular, regions that are normally enriched for H3K9 methylation show elevated recombination rates in met-2(n4256) mutants. We are currently investigating whether this change in the genetic map is the result of altered DSB distribution. Taken together, our data provide new insights into the role of chromatin in meiotic recombination. Because meiotic recombination generates genetic diversity, our data implicates chromatin as a key player in the evolution of genomes.