Evolution Of Genome Architecture Research Articles

Local Adaptation and the Evolution of Genome Architecture in Threespine Stickleback.

Theory predicts that local adaptation should favor the evolution of a concentrated genetic architecture, where the alleles driving adaptive divergence are tightly clustered on chromosomes. Adaptation to marine versus freshwater environments in threespine stickleback has resulted in an architecture that seems consistent with this prediction: divergence among populations is mainly driven by a few genomic regions harboring multiple quantitative trait loci for environmentally adapted traits, as well as candidate genes with well-established phenotypic effects. One theory for the evolution of these “genomic islands” is that rearrangements remodel the genome to bring causal loci into tight proximity, but this has not been studied explicitly. We tested this theory using synteny analysis to identify micro- and macro-rearrangements in the stickleback genome and assess their potential involvement in the evolution of genomic islands. To identify rearrangements, we conducted a de novo assembly of the closely related tubesnout (Aulorhyncus flavidus) genome and compared this to the genomes of threespine stickleback and two other closely related species. We found that small rearrangements, within-chromosome duplications, and lineage-specific genes (LSGs) were enriched around genomic islands, and that all three chromosomes harboring large genomic islands have experienced macro-rearrangements. We also found that duplicates and micro-rearrangements are 9.9× and 2.9× more likely to involve genes differentially expressed between marine and freshwater genotypes. While not conclusive, these results are consistent with the explanation that strong divergent selection on candidate genes drove the recruitment of rearrangements to yield clusters of locally adaptive loci.

Inverted base composition skews and discontinuous mitochondrial genome architecture evolution in the Enoplea (Nematoda)

BackgroundWithin the class Enoplea, the earliest-branching lineages in the phylum Nematoda, the relatively highly conserved ancestral mitochondrial architecture of Trichinellida is in stark contrast to the rapidly evolving architecture of Dorylaimida and Mermithida. To better understand the evolution of mitogenomic architecture in this lineage, we sequenced the mitogenome of a fish parasite Pseudocapillaria tomentosa (Trichinellida: Capillariidae) and compared it to all available enoplean mitogenomes.ResultsP. tomentosa exhibited highly reduced noncoding regions (the largest was 98 bp), and a unique base composition among the Enoplea. We attributed the latter to the inverted GC skew (0.08) in comparison to the ancestral skew in Trichinellidae (-0.43 to -0.37). Capillariidae, Trichuridae and Longidoridae (Dorylaimida) generally exhibited low negative or low positive skews (-0.1 to 0.1), whereas Mermithidae exhibited fully inverted low skews (0 to 0.05). This is indicative of inversions in the strand replication order or otherwise disrupted replication mechanism in the lineages with reduced/inverted skews. Among the Trichinellida, Trichinellidae and Trichuridae have almost perfectly conserved architecture, whereas Capillariidae exhibit multiple rearrangements of tRNA genes. In contrast, Mermithidae (Mermithida) and Longidoridae (Dorylaimida) exhibit almost no similarity to the ancestral architecture.ConclusionsLongidoridae exhibited more rearranged mitogenomic architecture than the hypervariable Mermithidae. Similar to the Chromadorea, the evolution of mitochondrial architecture in enoplean nematodes exhibits a strong discontinuity: lineages possessing a mostly conserved architecture over tens of millions of years are interspersed with lineages exhibiting architectural hypervariability. As Longidoridae also have some of the smallest metazoan mitochondrial genomes, they contradict the prediction that compact mitogenomes should be structurally stable. Lineages exhibiting inverted skews appear to represent the intermediate phase between the Trichinellidae (ancestral) and fully derived skews in Chromadorean mitogenomes (GC skews = 0.18 to 0.64). Multiple lines of evidence (CAT-GTR analysis in our study, a majority of previous mitogenomic results, and skew disruption scenarios) support the Dorylaimia split into two sister-clades: Dorylaimida + Mermithida and Trichinellida. However, skew inversions produce strong base composition biases, which can hamper phylogenetic and other evolutionary studies, so enoplean mitogenomes have to be used with utmost care in evolutionary studies.

Comparative 3D genome architecture in vertebrates.

BackgroundThe three-dimensional (3D) architecture of the genome has a highly ordered and hierarchical nature, which influences the regulation of essential nuclear processes at the basis of gene expression, such as gene transcription. While the hierarchical organization of heterochromatin and euchromatin can underlie differences in gene expression that determine evolutionary differences among species, the way 3D genome architecture is affected by evolutionary forces within major lineages remains unclear. Here, we report a comprehensive comparison of 3D genomes, using high resolution Hi-C data in fibroblast cells of fish, chickens, and 10 mammalian species.ResultsThis analysis shows a correlation between genome size and chromosome length that affects chromosome territory (CT) organization in the upper hierarchy of genome architecture, whereas lower hierarchical features, including local transcriptional availability of DNA, are selected through the evolution of vertebrates. Furthermore, conservation of topologically associating domains (TADs) appears strongly associated with the modularity of expression profiles across species. Additionally, LINE and SINE transposable elements likely contribute to heterochromatin and euchromatin organization, respectively, during the evolution of genome architecture.ConclusionsOur analysis uncovers organizational features that appear to determine the conservation and transcriptional regulation of functional genes across species. These findings can guide ongoing investigations of genome evolution by extending our understanding of the mechanisms shaping genome architecture.Supplementary InformationThe online version contains supplementary material available at 10.1186/s12915-022-01301-7.

Genome assembly of Danaus chrysippus and comparison with the Monarch Danaus plexippus.

Milkweed butterflies in the genus Danaus are studied in a diverse range of research fields including the neurobiology of migration, biochemistry of plant detoxification, host–parasite interactions, evolution of sex chromosomes, and speciation. We have assembled a nearly chromosomal genome for Danaus chrysippus (known as the African Monarch, African Queen, and Plain Tiger) using long-read sequencing data. This species is of particular interest for the study of genome structural change and its consequences for evolution. Comparison with the genome of the North American Monarch Danaus plexippus reveals generally strong synteny but highlights 3 inversion differences. The 3 chromosomes involved were previously found to carry peaks of intraspecific differentiation in D. chrysippus in Africa, suggesting that these inversions may be polymorphic and associated with local adaptation. The D. chrysippus genome is over 40% larger than that of D. plexippus, and nearly all of the additional ∼100 Megabases of DNA comprises repeats. Future comparative genomic studies within this genus will shed light on the evolution of genome architecture.

A chromosome-level genome assembly for the eastern fence lizard (Sceloporus undulatus), a reptile model for physiological and evolutionary ecology.

BackgroundHigh-quality genomic resources facilitate investigations into behavioral ecology, morphological and physiological adaptations, and the evolution of genomic architecture. Lizards in the genus Sceloporus have a long history as important ecological, evolutionary, and physiological models, making them a valuable target for the development of genomic resources.FindingsWe present a high-quality chromosome-level reference genome assembly, SceUnd1.0 (using 10X Genomics Chromium, HiC, and Pacific Biosciences data), and tissue/developmental stage transcriptomes for the eastern fence lizard, Sceloporus undulatus. We performed synteny analysis with other snake and lizard assemblies to identify broad patterns of chromosome evolution including the fusion of micro- and macrochromosomes. We also used this new assembly to provide improved reference-based genome assemblies for 34 additional Sceloporus species. Finally, we used RNAseq and whole-genome resequencing data to compare 3 assemblies, each representing an increased level of cost and effort: Supernova Assembly with data from 10X Genomics Chromium, HiRise Assembly that added data from HiC, and PBJelly Assembly that added data from Pacific Biosciences sequencing. We found that the Supernova Assembly contained the full genome and was a suitable reference for RNAseq and single-nucleotide polymorphism calling, but the chromosome-level scaffolds provided by the addition of HiC data allowed synteny and whole-genome association mapping analyses. The subsequent addition of PacBio data doubled the contig N50 but provided negligible gains in scaffold length.ConclusionsThese new genomic resources provide valuable tools for advanced molecular analysis of an organism that has become a model in physiology and evolutionary ecology.

Chromosomal Evolution in Chiroptera.

Chiroptera is the second largest order among mammals, with over 1300 species in 21 extant families. The group is extremely diverse in several aspects of its natural history, including dietary strategies, ecology, behavior and morphology. Bat genomes show ample chromosome diversity (from 2n = 14 to 62). As with other mammalian orders, Chiroptera is characterized by clades with low, moderate and extreme chromosomal change. In this article, we will discuss trends of karyotypic evolution within distinct bat lineages (especially Phyllostomidae, Hipposideridae and Rhinolophidae), focusing on two perspectives: evolution of genome architecture, modes of chromosomal evolution, and the use of chromosome data to resolve taxonomic problems.

Adaptation despite gene flow? Low recombination helps.

About 15,000years earlier, the Northern half of Europe and North America was buried under a few kilometres of ice. Since then, many organisms have colonized and rapidly adapted to the new, vacant habitats. Some, like the threespine stickleback fish, have done so more successfully than others: from the sea, stickleback have adapted to a multitude of lake and stream habitats with a vast array of complex phenotypes and life histories. Previous studies showed that most of these "ecotypes" differ in multiple divergently selected genes throughout the genome. But how are well-adapted ecotypes of one habitat protected from maladaptive gene flow from ecotypes of another, adjacent habitat? According to a From the Cover meta-analysis in this issue of Molecular Ecology (Samuk etal., 2017), low recombination rate regions in the genome offer such protection. While inversions have often been highlighted as an efficient way to maintain linkage disequilibrium among sets of adaptive variants in the face of gene flow, Samuk etal. (2017) show that variation in recombination rate across the genome may perform a similar role in threespine stickleback. With this study, theoretical predictions for the importance of low recombination regions in adaptation are for the first time tested with a highly replicated population genomic data set. The findings from this study have implications for the adaptability of species, speciation and the evolution of genome architecture.

High virulence does not necessarily impede viral adaptation to a new host: a case study using a plant RNA virus

BackgroundTheory suggests that high virulence could hinder between-host transmission of microparasites, and that virulence therefore will evolve to lower levels. Alternatively, highly virulent microparasites could also curtail host development, thereby limiting both the host resources available to them and their own within-host effective population size. In this case, high virulence might restrain the mutation supply rate and increase the strength with which genetic drift acts on microparasite populations. Thereby, this alternative explanation limits the microparasites’ potential to adapt to the host and ultimately the ability to evolve lower virulence. As a first exploration of this hypothesis, we evolved Tobacco etch virus carrying an eGFP fluorescent marker in two semi-permissive host species, Nicotiana benthamiana and Datura stramonium, for which it has a large difference in virulence. We compared the results to those previously obtained in the natural host, Nicotiana tabacum, where we have shown that carriage of eGFP has a high fitness cost and its loss serves as a real-time indicator of adaptation.ResultsAfter over half a year of evolution, we sequenced the genomes of the evolved lineages and measured their fitness. During the evolution experiment, marker loss leading to viable virus variants was only observed in one lineage of the host for which the virus has low virulence, D. stramonium. This result was consistent with the observation that there was a fitness cost of eGFP in this host, while surprisingly no fitness cost was observed in the host for which the virus has high virulence, N. benthamiana. Furthermore, in both hosts we observed increases in viral fitness in few lineages, and host-specific convergent evolution at the genomic level was only found in N. benthamiana.ConclusionsThe results of this study do not lend support to the hypothesis that high virulence impedes microparasites’ evolution. Rather, they exemplify that jumps between host species can be game changers for evolutionary dynamics. When considering the evolution of genome architecture, host species jumps might play a very important role, by allowing evolutionary intermediates to be competitive.

Detecting Horizontal Gene Transfer between Closely Related Taxa.

Horizontal gene transfer (HGT), the transfer of genetic material between organisms, is crucial for genetic innovation and the evolution of genome architecture. Existing HGT detection algorithms rely on a strong phylogenetic signal distinguishing the transferred sequence from ancestral (vertically derived) genes in its recipient genome. Detecting HGT between closely related species or strains is challenging, as the phylogenetic signal is usually weak and the nucleotide composition is normally nearly identical. Nevertheless, there is a great importance in detecting HGT between congeneric species or strains, especially in clinical microbiology, where understanding the emergence of new virulent and drug-resistant strains is crucial, and often time-sensitive.We developed a novel, self-contained technique named Near HGT, based on the synteny index, to measure the divergence of a gene from its native genomic environment and used it to identify candidate HGT events between closely related strains. The method confirms candidate transferred genes based on the constant relative mutability (CRM). Using CRM, the algorithm assigns a confidence score based on “unusual” sequence divergence. A gene exhibiting exceptional deviations according to both synteny and mutability criteria, is considered a validated HGT product. We first employed the technique to a set of three E. coli strains and detected several highly probable horizontally acquired genes. We then compared the method to existing HGT detection tools using a larger strain data set.When combined with additional approaches our new algorithm provides richer picture and brings us closer to the goal of detecting all newly acquired genes in a particular strain.

The Genomic Challenge to Adaptationism

Since the late 1990s, the characterization of complete DNA sequences for a large and taxonomically diverse set of species has continued to gain in speed and accuracy. Sequence analyses have indicated a strikingly baroque structure for most eukaryotic genomes, with multiple repeats of DNA sequences and with very little of the DNA specifying proteins. Much of the DNA in these genomes has no known function. These results have generated strong interest in the factors that govern the evolution of genome architecture. While adaptationist ‘just so’ stories have been offered (as typically occurs in every area of biology), recent theoretical analyses based on mathematical population genetics strongly suggest that non-adaptive processes dominate genome architecture evolution. This article critically synthesizes and develops these arguments, explicating a core argument along with several variants. It provides a critical assessment of the evidence that supports these arguments’ premises. It also analyses adaptationist responses to these arguments and notes potential problems with the core argument. These theoretical analyses continue the molecular reinterpretation of evolution initiated by the neutral theory in 1968. The article ends by noting that some of these arguments can also be extended to evolution at higher levels of organization which raises questions about adaptationism in general. This remains a puzzle because there is probably little reason to doubt that many organismic features are genuine adaptations.

Patterns of rDNA and telomeric sequences diversification: contribution to repetitive DNA organization in Phyllostomidae bats

Chromosomal organization and the evolution of genome architecture can be investigated by physical mapping of the genes for 45S and 5S ribosomal DNAs (rDNAs) and by the analysis of telomeric sequences. We studied 12 species of bats belonging to four subfamilies of the family Phyllostomidae in order to correlate patterns of distribution of heterochromatin and the multigene families for rDNA. The number of clusters for 45S gene ranged from one to three pairs, with exclusively location in autosomes, except for Carollia perspicillata that had in X chromosome. The 5S gene all the species studied had only one site located on an autosomal pair. In no species the 45S and 5S genes collocated. The fluorescence in situ hybridization (FISH) probe for telomeric sequences revealed fluorescence on all telomeres in all species, except in Carollia perspicillata. Non-telomeric sites in the pericentromeric region of the chromosomes were observed in most species, ranged from one to 12 pairs. Most interstitial telomeric sequences were coincident with heterochromatic regions. The results obtained in the present work indicate that different evolutionary mechanisms are acting in Phyllostomidae genome architecture, as well as the occurrence of Robertsonian fusion during the chromosomal evolution of bats without a loss of telomeric sequences. These data contribute to understanding the organization of multigene families and telomeric sequences on bat genome as well as the chromosomal evolutionary history of Phyllostomidae bats.

Sexual and parental antagonism shape genomic architecture

Populations with two sexes are vulnerable to a pair of genetic conflicts: sexual antagonism that can arise when alleles have opposing fitness effects on females and males; and parental antagonism that arises when alleles have opposing fitness effects when maternally and paternally inherited. This paper extends previous theoretical work that found stable linkage disequilibrium (LD) between sexually antagonistic loci. We find that LD is also generated between parentally antagonistic loci, and between sexually and parentally antagonistic loci, without any requirement of epistasis. We contend that the LD in these models arises from the admixture of gene pools subject to different selective histories. We also find that polymorphism maintained by parental antagonism at one locus expands the opportunity for polymorphism at a linked locus experiencing parental or sexual antagonism. Taken together, our results predict the chromosomal clustering of loci that segregate for sexually and parentally antagonistic alleles. Thus, genetic conflict may play a role in the evolution of genomic architecture.

Bootstrapping phylogenies inferred from rearrangement data.

BackgroundLarge-scale sequencing of genomes has enabled the inference of phylogenies based on the evolution of genomic architecture, under such events as rearrangements, duplications, and losses. Many evolutionary models and associated algorithms have been designed over the last few years and have found use in comparative genomics and phylogenetic inference. However, the assessment of phylogenies built from such data has not been properly addressed to date. The standard method used in sequence-based phylogenetic inference is the bootstrap, but it relies on a large number of homologous characters that can be resampled; yet in the case of rearrangements, the entire genome is a single character. Alternatives such as the jackknife suffer from the same problem, while likelihood tests cannot be applied in the absence of well established probabilistic models.ResultsWe present a new approach to the assessment of distance-based phylogenetic inference from whole-genome data; our approach combines features of the jackknife and the bootstrap and remains nonparametric. For each feature of our method, we give an equivalent feature in the sequence-based framework; we also present the results of extensive experimental testing, in both sequence-based and genome-based frameworks. Through the feature-by-feature comparison and the experimental results, we show that our bootstrapping approach is on par with the classic phylogenetic bootstrap used in sequence-based reconstruction, and we establish the clear superiority of the classic bootstrap for sequence data and of our corresponding new approach for rearrangement data over proposed variants. Finally, we test our approach on a small dataset of mammalian genomes, verifying that the support values match current thinking about the respective branches.ConclusionsOur method is the first to provide a standard of assessment to match that of the classic phylogenetic bootstrap for aligned sequences. Its support values follow a similar scale and its receiver-operating characteristics are nearly identical, indicating that it provides similar levels of sensitivity and specificity. Thus our assessment method makes it possible to conduct phylogenetic analyses on whole genomes with the same degree of confidence as for analyses on aligned sequences. Extensions to search-based inference methods such as maximum parsimony and maximum likelihood are possible, but remain to be thoroughly tested.

Formation, regulation and evolution of Caenorhabditis elegans 3′UTRs

Posttranscriptional gene regulation frequently occurs through elements in mRNA 3′ untranslated regions (UTRs)1,2. Although crucial roles for 3′UTR-mediated gene regulation have been found in Caenorhabditis elegans3,4,5, most C. elegans genes have lacked annotated 3′UTRs6,7. Here we describe a high-throughput method to reliably identify polyadenylated RNA termini, and we apply this method, called poly(A)-position profiling by sequencing (3P-Seq), to determine C. elegans 3′UTRs. Compared to standard methods also recently applied to C. elegans UTRs8, 3P-Seq identified 8,581 additional UTRs while excluding thousands of shorter UTR isoforms that do not appear to be authentic. Analysis of this expanded and corrected dataset suggested that the high A/U content of C. elegans 3′UTRs facilitated genome compaction, since the elements specifying cleavage and polyadenylation, which are A/U-rich, can more readily emerge in A/U rich regions. Indeed, 30% of the protein-coding genes have mRNAs with alternative, partially overlapping end regions that generate another 10,498 cleavage and polyadenylation sites that had gone largely unnoticed and represent potential evolutionary intermediates of progressive UTR shortening. Moreover, a third of the convergently transcribed genes utilize palindromic arrangements of bidirectional elements to specify UTRs with convergent overlap, which also contributes to genome compaction by eliminating regions between genes. Although nematode 3′UTRs have median length only one-sixth that of mammalian 3′UTRs, they have twice the density of conserved microRNA sites, in part because additional types of seed-complementary sites are preferentially conserved. These findings reveal the influence of cleavage and polyadenylation on the evolution of genome architecture and provide resources for studying posttranscriptional gene regulation.

Origin of avian genome size and structure in non-avian dinosaurs

Avian genomes are small and streamlined compared with those of other amniotes by virtue of having fewer repetitive elements and less non-coding DNA. This condition has been suggested to represent a key adaptation for flight in birds, by reducing the metabolic costs associated with having large genome and cell sizes. However, the evolution of genome architecture in birds, or any other lineage, is difficult to study because genomic information is often absent for long-extinct relatives. Here we use a novel bayesian comparative method to show that bone-cell size correlates well with genome size in extant vertebrates, and hence use this relationship to estimate the genome sizes of 31 species of extinct dinosaur, including several species of extinct birds. Our results indicate that the small genomes typically associated with avian flight evolved in the saurischian dinosaur lineage between 230 and 250 million years ago, long before this lineage gave rise to the first birds. By comparison, ornithischian dinosaurs are inferred to have had much larger genomes, which were probably typical for ancestral Dinosauria. Using comparative genomic data, we estimate that genome-wide interspersed mobile elements, a class of repetitive DNA, comprised 5-12% of the total genome size in the saurischian dinosaur lineage, but was 7-19% of total genome size in ornithischian dinosaurs, suggesting that repetitive elements became less active in the saurischian lineage. These genomic characteristics should be added to the list of attributes previously considered avian but now thought to have arisen in non-avian dinosaurs, such as feathers, pulmonary innovations, and parental care and nesting.

Adaptive role of increased frequency of polypurine tracts in mRNA sequences of thermophilic prokaryotes.

The mechanism of an organism's adaptation to high temperatures has been investigated intensively in recent years. It was suggested that the macromolecules of thermophilic microorganisms (especially proteins) have structural features that enhance their thermostability. We compared mRNA sequences of 72 fully sequenced prokaryotic proteomes (14 thermophilic and 58 mesophilic species). Although the differences between the percentage of adenine plus guanine content of whole mRNAs of different prokaryotic species are much lower than those of guanine plus cytosine content, the thermophile purine-pyrimidine (R/Y) ratio within their mRNAs is significantly higher than that of the mesophiles. The first and third codon positions of both thermophiles and mesophiles are purine-biased, with the bias more pronounced by the thermophiles. Thermophile mRNAs that display the highest R/Y ratio (1.43-1.69) are those of the ribosomal proteins, histone-like proteins, DNA-dependent RNA polymerase subunits, and heat-shock proteins. Within mesophilic prokaryotes and five eukaryotic species, the R/Y ratio of the mRNAs of heat-shock proteins is higher than their average over coding part of the genome. Polypurine tracts (R)(n) (with n > or = 5) are much more abundant within the thermophile mRNAs compared with mesophiles. Between two sequential pure-purinic codons of thermophile mRNAs, there is a rather strong tendency for the occurrence of adenine but not guanine tracts. The data suggest that mixed adenine.guanine and polyadenine tracts in mRNAs increase the thermostability beyond the contribution of amino acids encoded by purine tracts, which highlights the importance of ecological stress in the evolution of genome architecture.

Correlational selection and the evolution of genomic architecture

We review and discuss the importance of correlational selection (selection for optimal character combinations) in natural populations. If two or more traits subject to multivariate selection are heritable, correlational selection builds favourable genetic correlations through the formation of linkage disequilibrium at underlying loci governing the traits. However, linkage disequilibria built up by correlational selection are expected to decay rapidly (ie, within a few generations), unless correlational selection is strong and chronic. We argue that frequency-dependent biotic interactions that have 'Red Queen dynamics' (eg, host-parasite interactions, predator-prey relationships or intraspecific arms races) often fuel chronic correlational selection, which is strong enough to maintain adaptive genetic correlations of the kind we describe. We illustrate these processes and phenomena using empirical examples from various plant and animal systems, including our own recent work on the evolutionary dynamics of a heritable throat colour polymorphism in the side-blotched lizard Uta stansburiana. In particular, male and female colour morphs of side-blotched lizards cycle on five- and two-generation (year) timescales under the force of strong frequency-dependent selection. Each morph refines the other morph in a Red Queen dynamic. Strong correlational selection gradients among life history, immunological and morphological traits shape the genetic correlations of the side-blotched lizard polymorphism. We discuss the broader evolutionary consequences of the buildup of co-adapted trait complexes within species, such as the implications for speciation processes.