Crossing the Rift valley: using complete mitogenomes to infer the diversification and biogeographic history of ethiopian highlands Ptychadena (anura: Ptychadenidae)

The Ethiopian Highlands are considered a biodiversity hotspot, harboring a high number of endemic species. Some of the endemic species probably diversified in situ; this is, for example, the case of a monophyletic clade containing 12 known species of grass frogs of the genus Ptychadena. The different species occur at elevations ranging from 1,500 to above 3,400 m and constitute excellent models to study the process of diversification in the highlands as well as adaptations to high elevations. In this study, we sampled 294 specimens across the distribution of this clade and used complete mitogenomes and genome-wide SNP data to better understand how landscape features influenced the population structure and dispersal of these grass frogs across time and space. Using phylogenetic inference, population structure analyses, and biogeographic reconstructions, we found that the species complex probably first diversified on the south-east side of the Great Rift Valley. Later on, species dispersed to the north-west side, where more recent diversification occurred. We further demonstrate that Ptychadena species have dispersed across the Great Rift Valley at different times. Our analyses allowed for a more complete understanding of the contribution of geological events, biogeographic barriers and climatic changes as drivers of species diversification and adaptation in this important biogeographic region.

Genetic and morphometric analyses of historical type specimens clarify the taxonomy of the Ethiopian Leptopelisgramineus species complex (Anura, Arthroleptidae)

Frogs of the genus Leptopelis have diversified in the Ethiopian Highlands to occupy forests and montane grasslands both east and west of the Great Rift Valley. Genetic studies revealed that the endemic species Leptopelisgramineus (Boulenger, 1898) comprises multiple unnamed taxa. A careful examination of historical type specimens is, however, needed to fully resolve the taxonomy of the group. Here we use mitochondrial DNA and morphological analyses on a large sample of recently-collected Ethiopian Leptopelis, as well as century-old type specimens to demonstrate that the recently resurrected L.montanus Tiutenko & Zinenko, 2021 (previously Pseudocassinaocellata Ahl, 1924) is a junior synonym of L.rugosus (Ahl, 1924) and corresponds to the taxon found west of the Great Rift Valley, not east as previously thought. Our results show that populations inhabiting the mountains and plateaus east of the Rift constitute a distinct and undescribed species. We provide a re-description of L.rugosus and describe two new species inhabiting the Highlands east of the Great Rift Valley. We provide an identification key, as well as a description of the calls of the members of the Leptopelisgramineus species complex.

Genomic analyses reveal association of ASIP with a recurrently-evolving adaptive color pattern in frogs

Traits shared among distantly related lineages are indicators of common evolutionary constraints, at the ecological, physiological, or molecular level. Here, we show that the vertebral stripe, a cryptic color pattern, has evolved hundreds of times in the evolutionary history of anurans (frogs and toads) and is favored in terrestrial habitats. Using a genome-wide association study, we demonstrate that variation near the Agouti signaling protein gene (ASIP) is responsible for the different vertebral stripe phenotypes in the African grass frog Ptychadena robeensis. RNAseq and real-time quantitative PCR revealed that differential expression of the gene and an adjacent long non-coding RNA is linked to patterning in this species. Surprisingly, and although the stripe phenotypes are shared with closely related species, we found that the P. robeensis alleles are private to the species and unlikely to evolve under long-term balancing selection, thus indicating that the vertebral stripe phenotypes result from parallel evolution within the group. Our findings demonstrate that this cryptic color pattern evolved rapidly and recurrently in terrestrial anurans, and therefore constitutes an ideal system to study repeated evolution.

Linked-Read Sequencing of Eight Falcons Reveals a Unique Genomic Architecture in Flux

Falcons are diverse birds of cultural and economic importance. They have undergone major lineage-specific chromosomal rearrangements, resulting in greatly-reduced chromosome counts relative to other birds. Here, we use 10X Genomics linked reads to provide new high-contiguity genomes for two gyrfalcons, a saker falcon, a lanner falcon, three subspecies of peregrine falcons, and the common kestrel. Assisted by a transcriptome sequenced from 22 gyrfalcon tissues, we annotate these genomes for a variety of genomic features, estimate historical demography, and then investigate genomic equilibrium in the context of falcon-specific chromosomal rearrangements. We find that falcon genomes are not in AT-GC equilibrium with a bias in substitutions towards higher AT content; this bias is predominantly but not exclusively driven by hypermutability of CpG sites. Small indels and large structural variants were also biased towards insertions rather than deletions. Patterns of disequilibrium were linked to chromosomal rearrangements: falcons have lost GC content in regions that have fused to larger chromosomes from microchromosomes and gained GC content in regions of macrochromosomes that have translocated to microchromosomes. Inserted bases have accumulated on regions ancestrally belonging to microchromosomes, consistent with insertion-biased gene conversion. We also find an excess of interspersed repeats on regions of microchromosomes that have fused to macrochromosomes. Our results reveal that falcon genomes are in a state of flux. They further suggest that many of the key differences between microchromosomes and macrochromosomes are driven by differences in chromosome size, and indicate a clear role for recombination and biased-gene-conversion in determining genomic equilibrium.

Varied diversification patterns and distinct demographic trajectories in Ethiopian montane forest bird (Aves: Passeriformes) populations separated by the Great Rift Valley

Taxon-specific characteristics and extrinsic climatic and geological forces may both shape population differentiation and speciation. In geographically and taxonomically focused investigations, differentiation may occur synchronously as species respond to the same external conditions. Conversely, when evolution is investigated in taxa with largely varying traits, population differentiation and speciation is complex and shaped by interactions of Earth's template and species-specific traits. As such, it is important to characterize evolutionary histories broadly across the tree of life, especially in geographic regions that are exceptionally diverse and under pressures from human activities such as in biodiversity hotspots. Here, using whole-genome sequencing data, we characterize genomic variation in populations of six Ethiopian Highlands forest bird species separated by a lowland biogeographic barrier, the Great Rift Valley (GRV). In all six species, populations on either side of the GRV exhibited significant but varying levels of genetic differentiation. Species' dispersal ability was negatively correlated with levels of population differentiation. Isolation with migration models indicated varied patterns of population differentiation and connectivity among populations of the focal species. We found that demographic histories-estimated for each individual-varied by both species and population but were consistent between individuals of the same species and sampling region. We found that genomic diversity varied by half an order of magnitude across species, and that this variation could largely be explained by the harmonic mean of effective population size over the past 200,000 years. Overall, we found that even in highly dispersive species like birds, the GRV acts as a substantial biogeographic barrier.

Mitogenomics of historical type specimens clarifies the taxonomy of Ethiopian Ptychadena Boulenger, 1917 (Anura, Ptychadenidae)

The taxonomy of the Ptychadenaneumanni species complex, a radiation of grass frogs inhabiting the Ethiopian highlands, has puzzled scientists for decades because of the morphological resemblance among its members. Whilst molecular phylogenetic methods allowed the discovery of several species in recent years, assigning pre-existing and new names to clades was challenged by the unavailability of molecular data for century-old type specimens. We used Illumina short reads to sequence the mitochondrial DNA of type specimens in this group, as well as ddRAD-seq analyses to resolve taxonomic uncertainties surrounding the P.neumanni species complex. The phylogenetic reconstruction revealed recurrent confusion between Ptychadenaerlangeri (Ahl, 1924) and P.neumanni (Ahl, 1924) in the literature. The phylogeny also established that P.largeni Perret, 1994 represents a junior synonym of P.erlangeri (Ahl, 1924) and distinguished between two small species, P.nana Perret, 1994, restricted to the Arussi Plateau, and P.robeensis Goutte, Reyes-Velasco, Freilich, Kassie & Boissinot, 2021, which inhabits the Bale Mountains. The phylogenetic analyses of mitochondrial DNA from type specimens also corroborate the validity of seven recently described species within the group. Our study shows how modern molecular tools applied to historical type specimens can help resolve long-standing taxonomic issues in cryptic species complexes.

The Structural, Functional and Evolutionary Impact of Transposable Elements in Eukaryotes

Transposable elements (TEs) are nearly ubiquitous in eukaryotes. The increase in genomic data, as well as progress in genome annotation and molecular biology techniques, have revealed the vast number of ways mobile elements have impacted the evolution of eukaryotes. In addition to being the main cause of difference in haploid genome size, TEs have affected the overall organization of genomes by accumulating preferentially in some genomic regions, by causing structural rearrangements or by modifying the recombination rate. Although the vast majority of insertions is neutral or deleterious, TEs have been an important source of evolutionary novelties and have played a determinant role in the evolution of fundamental biological processes. TEs have been recruited in the regulation of host genes and are implicated in the evolution of regulatory networks. They have also served as a source of protein-coding sequences or even entire genes. The impact of TEs on eukaryotic evolution is only now being fully appreciated and the role they may play in a number of biological processes, such as speciation and adaptation, remains to be deciphered.

Mitogenome analyses elucidate the evolutionary relationships of a probable Eocene wet tropics relic in the xerophile lizard genus Acanthodactylus

Climate has a large impact on diversity and evolution of the world's biota. The Eocene-Oligocene transition from tropical climate to cooler, drier environments was accompanied by global species turnover. A large number of Old World lacertid lizard lineages have diversified after the Eocene-Oligocene boundary. One of the most speciose reptile genera in the arid Palearctic, Acanthodactylus, contains two sub-Saharan species with unresolved phylogenetic relationship and unknown climatic preferences. We here aim to understand how and when adaptation to arid conditions occurred in Acanthodactylus and when tropical habitats where entered. Using whole mitogenomes from fresh and archival DNA and published sequences we recovered a well-supported Acanthodactylus phylogeny and underpinned the timing of diversification with environmental niche analyses of the sub-Saharan species A. guineensis and A. boueti in comparison to all arid Acanthodactylus. We found that A. guineensis represents an old lineage that splits from a basal node in the Western clade, and A. boueti is a derived lineage and probably not its sister. Their long branches characterize them-and especially A. guineensis-as lineages that may have persisted for a long time without further diversification or have undergone multiple extinctions. Environmental niche models verified the occurrence of A. guineensis and A. boueti in hot humid environments different from the other 42 arid Acanthodactylus species. While A. guineensis probably remained in tropical habitat from periods prior to the Eocene-Oligocene boundary, A. boueti entered tropical environments independently at a later period. Our results provide an important baseline for studying adaptation and the transition from humid to arid environments in Lacertidae.

Disentangling the determinants of transposable elements dynamics in vertebrate genomes using empirical evidences and simulations

The interactions between transposable elements (TEs) and their hosts constitute one of the most profound co-evolutionary processes found in nature. The population dynamics of TEs depends on factors specific to each TE families, such as the rate of transposition and insertional preference, the demographic history of the host and the genomic landscape. How these factors interact has yet to be investigated holistically. Here we are addressing this question in the green anole (Anolis carolinensis) whose genome contains an extraordinary diversity of TEs (including non-LTR retrotransposons, SINEs, LTR-retrotransposons and DNA transposons). We observed a positive correlation between recombination rate and frequency of TEs and densities for LINEs, SINEs and DNA transposons. For these elements, there was a clear impact of demography on TE frequency and abundance, with a loss of polymorphic elements and skewed frequency spectra in recently expanded populations. On the other hand, some LTR-retrotransposons displayed patterns consistent with a very recent phase of intense amplification. To determine how demography, genomic features and intrinsic properties of TEs interact we ran simulations using SLiM3. We determined that i) short TE insertions are not strongly counter-selected, but long ones are, ii) neutral demographic processes, linked selection and preferential insertion may explain positive correlations between average TE frequency and recombination, iii) TE insertions are unlikely to have been massively recruited in recent adaptation. We demonstrate that deterministic and stochastic processes have different effects on categories of TEs and that a combination of empirical analyses and simulations can disentangle these mechanisms.

Variation in base composition underlies functional and evolutionary divergence in non-LTR retrotransposons

Background

Non-LTR retrotransposons often exhibit base composition that is markedly different from the nucleotide content of their host’s gene. For instance, the mammalian L1 element is AT-rich with a strong A bias on the positive strand, which results in a reduced transcription. It is plausible that the A-richness of mammalian L1 is a self-regulatory mechanism reflecting a trade-off between transposition efficiency and the deleterious effect of L1 on its host. We examined if the A-richness of L1 is a general feature of non-LTR retrotransposons or if different clades of elements have evolved different nucleotide content. We also investigated if elements belonging to the same clade evolved towards different base composition in different genomes or if elements from different clades evolved towards similar base composition in the same genome.

Results

We found that non-LTR retrotransposons differ in base composition among clades within the same host but also that elements belonging to the same clade differ in base composition among hosts. We showed that nucleotide content remains constant within the same host over extended period of evolutionary time, despite mutational patterns that should drive nucleotide content away from the observed base composition.

Conclusions

Our results suggest that base composition is evolving under selection and may be reflective of the long-term co-evolution between non-LTR retrotransposons and their host. Finally, the coexistence of elements with drastically different base composition suggests that these elements may be using different strategies to persist and multiply in the genome of their host.

Falcon genomics in the context of conservation, speciation, and human culture

Here, we review the diversity, evolutionary history, and genomics of falcons in the context of their conservation and interactions with humans, and provide a perspective on how new genomic approaches may be applied to expand our knowledge of these topics. For millennia, humans and falcons (genus Falco) have developed unique relationships through falconry, religious rituals, conservation efforts, and human lifestyle transitions. From an evolutionary perspective, falcons remain an enigma. Having experienced several recent radiations, they have reached an unparalleled and almost global distribution, with an intrageneric species richness that is roughly an order of magnitude higher than typical within their family (Falconidae) and across other birds (Phylum: Aves). This diversity has evolved in the context of unusual genomic architecture that includes unique chromosomal rearrangements, relatively low chromosome counts, extremely low microdeletion rates, and high levels of nuclear mitochondrial DNA segments (NUMTs). These genomic peculiarities combine with high levels of ecological and organismal diversity and a legacy of human interactions to make falcons obvious candidates for evolutionary studies, providing unique research opportunities in common topics, including chromosomal evolution, the mechanics of speciation, local adaptation, domestication, and urban adaptation.

Selection at behavioural, developmental and metabolic genes is associated with the northward expansion of a successful tropical colonizer

What makes a species able to colonize novel environments? This question is key to understand the dynamics of adaptive radiations and ecological niche shifts, but the mechanisms that underlie expansion into novel habitats remain poorly understood at a genomic scale. Lizards from the genus Anolis are typically tropical, and the green anole (Anolis carolinensis) constitutes an exception since it expanded into temperate North America from subtropical Florida. Thus, we used the green anole as a model to investigate signatures of selection associated with colonization of a new environment, namely temperate North America. To this end, we analysed 29 whole-genome sequences, covering the entire native range of the species. We used a combination of recent methods to quantify both positive and balancing selection in northern populations, including F_ST outlier methods, machine learning and ancestral recombination graphs. We naively scanned for genes of interest and assessed the overlap between multiple tests. Strikingly, we identified many genes involved in behaviour, suggesting that the recent successful colonization of northern environments may have been linked to behavioural shifts as well as physiological adaptation. Using a candidate genes strategy, we determined that genes involved in response to cold or behaviour displayed more frequently signals of selection, while controlling for local recombination rate, gene clustering and gene length. In addition, we found signatures of balancing selection at immune genes in all investigated genetic groups, but also at genes involved in neuronal and anatomical development.

Recent Secondary Contacts, Linked Selection, and Variable Recombination Rates Shape Genomic Diversity in the Model Species Anolis carolinensis

Gaining a better understanding on how selection and neutral processes affect genomic diversity is essential to gain better insights into the mechanisms driving adaptation and speciation. However, the evolutionary processes affecting variation at a genomic scale have not been investigated in most vertebrate lineages. Here, we present the first population genomics survey using whole genome resequencing in the green anole (Anolis carolinensis). Anoles have been intensively studied to understand mechanisms underlying adaptation and speciation. The green anole in particular is an important model to study genome evolution. We quantified how demography, recombination, and selection have led to the current genetic diversity of the green anole by using whole-genome resequencing of five genetic clusters covering the entire species range. The differentiation of green anole's populations is consistent with a northward expansion from South Florida followed by genetic isolation and subsequent gene flow among adjacent genetic clusters. Dispersal out-of-Florida was accompanied by a drastic population bottleneck followed by a rapid population expansion. This event was accompanied by male-biased dispersal and/or selective sweeps on the X chromosome. We show that the interaction between linked selection and recombination is the main contributor to the genomic landscape of differentiation in the anole genome.

On the Population Dynamics of Junk: A Review on the Population Genomics of Transposable Elements

Transposable elements (TEs) play an important role in shaping genomic organization and structure, and may cause dramatic changes in phenotypes. Despite the genetic load they may impose on their host and their importance in microevolutionary processes such as adaptation and speciation, the number of population genetics studies focused on TEs has been rather limited so far compared to single nucleotide polymorphisms (SNPs). Here, we review the current knowledge about the dynamics of transposable elements at recent evolutionary time scales, and discuss the mechanisms that condition their abundance and frequency. We first discuss non-adaptive mechanisms such as purifying selection and the variable rates of transposition and elimination, and then focus on positive and balancing selection, to finally conclude on the potential role of TEs in causing genomic incompatibilities and eventually speciation. We also suggest possible ways to better model TEs dynamics in a population genomics context by incorporating recent advances in TEs into the rich information provided by SNPs about the demography, selection, and intrinsic properties of genomes.

The Mobilome of Reptiles: Evolution, Structure, and Function

Transposable elements (TE) constitute one of the most variable genomic features among vertebrates, impacting genome size, structure, and composition. Despite their important role in shaping genomic diversity, they have mostly been studied in mammals, which display one of the least diverse genomes in terms of TE diversity. Recent new resources in reptilian genomics have opened a broader perspective about TE evolution in amniotes. We discuss these recent results by showing that TE diversity is high in reptiles, particularly in squamates, with strong heterogeneity in the number of TE classes retained in each lineage, even at short evolutionary scales. More research is needed to uncover the exact mechanisms that regulate TE proliferation in reptiles and to what extent these selfish elements can play a role in local adaptation or in the emergence of barriers to gene flow.

Multiple and Independent Phases of Transposable Element Amplification in the Genomes of Piciformes (Woodpeckers and Allies)

The small and conserved genomes of birds are likely a result of flight-related metabolic constraints. Recombination-driven deletions and minimal transposable element (TE) expansions have led to continually shrinking genomes during evolution of many lineages of volant birds. Despite constraints of genome size in birds, we identified multiple waves of amplification of TEs in Piciformes (woodpeckers, honeyguides, toucans, and barbets). Relative to other bird species’ genomic TE abundance (< 10% of genome), we found ∼17–30% TE content in multiple clades within Piciformes. Several families of the retrotransposon superfamily chicken repeat 1 (CR1) expanded in at least three different waves of activity. The most recent CR1 expansions (∼4–7% of genome) preceded bursts of diversification in the woodpecker clade and in the American barbets + toucans clade. Additionally, we identified several thousand polymorphic CR1 insertions (hundreds per individual) in three closely related woodpecker species. Woodpecker CR1 insertion polymorphisms are maintained at lower frequencies than single nucleotide polymorphisms indicating that purifying selection is acting against additional CR1 copies and that these elements impose a fitness cost on their host. These findings provide evidence of large scale and ongoing TE activity in avian genomes despite continual constraint on genome size.

Contrasting Rates of LINE-1 Amplification among New World Primates of the Atelidae Family

LINE-1 (L1) retrotransposons constitute the dominant category of transposons in mammalian genomes. L1 elements are active in the vast majority of mammals, and only a few cases of L1 extinction have been documented. The only possible case of extinction in primates was suggested for South American spider monkeys. However, these previous studies were based on a single species. We revisited this question with a larger phylogenetic sample, covering all 4 genera of Atelidae and 3 species of spider monkeys. We used an enrichment method to clone recently inserted L1 elements and performed an evolutionary analysis of the sequences. We were able to identify young L1 elements in all taxa, suggesting that L1 is probably still active in all Atelidae examined. However, we also detected considerable variations in the proportion of recent elements indicating that the rate of L1 amplification varies among Atelidae by a 3-fold factor. The extent of L1 amplification in Atelidae remains overall lower than in other New World monkeys. Multiple factors can affect the amplification of L1, such as the demography of the host and the control of transposition. These factors are discussed in the context of host life history.

Differential Effect of Selection against LINE Retrotransposons among Vertebrates Inferred from Whole-Genome Data and Demographic Modeling

Variation in LINE composition is one of the major determinants for the substantial size and structural differences among vertebrate genomes. In particular, the larger genomes of mammals are characterized by hundreds of thousands of copies from a single LINE clade, L1, whereas nonmammalian vertebrates possess a much greater diversity of LINEs, yet with orders of magnitude less in copy number. It has been proposed that such variation in copy number among vertebrates is due to differential effect of LINE insertions on host fitness. To investigate LINE selection, we deployed a framework of demographic modeling, coalescent simulations, and probabilistic inference against population-level whole-genome data sets for four model species: one population each of threespine stickleback, green anole, and house mouse, as well as three human populations. Specifically, we inferred a null demographic background utilizing SNP data, which was then exploited to simulate a putative null distribution of summary statistics that was compared with LINE data. Subsequently, we applied the inferred null demographic model with an additional exponential size change parameter, coupled with model selection, to test for neutrality as well as estimate the strength of either negative or positive selection. We found a robust signal for purifying selection in anole and mouse, but a lack of clear evidence for selection in stickleback and human. Overall, we demonstrated LINE insertion dynamics that are not in accordance to a mammalian versus nonmammalian dichotomy, and instead the degree of existing LINE activity together with host-specific demographic history may be the main determinants of LINE abundance.

Contrasted patterns of evolution of the LINE-1 retrotransposon in perissodactyls: the history of a LINE-1 extinction

Background

LINE-1 (L1) is the dominant autonomously replicating non-LTR retrotransposon in mammals. Although our knowledge of L1 evolution across the tree of life has considerably improved in recent years, what we know of L1 evolution in mammals is biased and comes mostly from studies in primates (mostly human) and rodents (mostly mouse). It is unclear if patterns of evolution that are shared between those two groups apply to other mammalian orders. Here we performed a detailed study on the evolution of L1 in perissodactyls by making use of the complete genome of the domestic horse and of the white rhinoceros. This mammalian order offers an excellent model to study the extinction of L1 since the rhinoceros is one of the few mammalian species to have lost active L1.

Results

We found that multiple L1 lineages, carrying different 5’UTRs, have been simultaneously active during the evolution of perissodactyls. We also found that L1 has continuously amplified and diversified in horse. In rhinoceros, L1 was very prolific early on. Two successful families were simultaneously active until ~20my ago but became extinct suddenly at exactly the same time.

Conclusions

The general pattern of L1 evolution in perissodactyls is very similar to what was previously described in mouse and human, suggesting some commonalities in the way mammalian genomes interact with L1. We confirmed the extinction of L1 in rhinoceros and we discuss several possible mechanisms.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0117-4) contains supplementary material, which is available to authorized users.

Revisiting the phylogeography, demography and taxonomy of the frog genus Ptychadena in the Ethiopian highlands with the use of genome-wide SNP data

Understanding the diversification of biological lineages is central to evolutionary studies. To properly study the process of speciation, it is necessary to link micro-evolutionary studies with macro-evolutionary mechanisms. Micro-evolutionary studies require proper sampling across a taxon's range to adequately infer genetic diversity. Here we use the grass frogs of the genus Ptychadena from the Ethiopian highlands as a model to study the process of lineage diversification in this unique biodiversity hotspot. We used thousands of genome-wide SNPs obtained from double digest restriction site associated DNA sequencing (ddRAD-seq) in populations of the Ptychadena neumanni species complex from the Ethiopian highlands in order to infer their phylogenetic relationships and genetic structure, as well as to study their demographic history. Our genome-wide phylogenetic study supports the existence of approximately 13 lineages clustered into 3 species groups. Our phylogenetic and phylogeographic reconstructions suggest that those endemic lineages diversified in allopatry, and subsequently specialized to different habitats and elevations. Demographic analyses point to a continuous decrease in the population size across the majority of lineages and populations during the Pleistocene, which is consistent with a continuous period of aridification that East Africa experienced since the Pliocene. We discuss the taxonomic implications of our analyses and, in particular, we warn against the recent practice to solely use Bayesian species delimitation methods when proposing taxonomic changes.

The Evolution of LINE-1 in Vertebrates

The abundance and diversity of the LINE-1 (L1) retrotransposon differ greatly among vertebrates. Mammalian genomes contain hundreds of thousands L1s that have accumulated since the origin of mammals. A single group of very similar elements is active at a time in mammals, thus a single lineage of active families has evolved in this group. In contrast, non-mammalian genomes (fish, amphibians, reptiles) harbor a large diversity of concurrently transposing families, which are all represented by very small number of recently inserted copies. Why the pattern of diversity and abundance of L1 is so different among vertebrates remains unknown. To address this issue, we performed a detailed analysis of the evolution of active L1 in 14 mammals and in 3 non-mammalian vertebrate model species. We examined the evolution of base composition and codon bias, the general structure, and the evolution of the different domains of L1 (5′UTR, ORF1, ORF2, 3′UTR). L1s differ substantially in length, base composition, and structure among vertebrates. The most variation is found in the 5′UTR, which is longer in amniotes, and in the ORF1, which tend to evolve faster in mammals. The highly divergent L1 families of lizard, frog, and fish share species-specific features suggesting that they are subjected to the same functional constraints imposed by their host. The relative conservation of the 5′UTR and ORF1 in non-mammalian vertebrates suggests that the repression of transposition by the host does not act in a sequence-specific manner and did not result in an arms race, as is observed in mammals.

LINE Insertion Polymorphisms are Abundant but at Low Frequencies across Populations of Anolis carolinensis

Vertebrate genomes differ considerably in size and structure. Among the features that show the most variation is the abundance of Long Interspersed Nuclear Elements (LINEs). Mammalian genomes contain 100,000s LINEs that belong to a single clade, L1, and in most species a single family is usually active at a time. In contrast, non-mammalian vertebrates (fish, amphibians and reptiles) contain multiple active families, belonging to several clades, but each of them is represented by a small number of recently inserted copies. It is unclear why vertebrate genomes harbor such drastic differences in LINE composition. To address this issue, we conducted whole genome resequencing to investigate the population genomics of LINEs across 13 genomes of the lizard Anolis carolinensis sampled from two geographically and genetically distinct populations in the Eastern Florida and the Gulf Atlantic regions of the United States. We used the Mobile Element Locator Tool to identify and genotype polymorphic insertions from five major clades of LINEs (CR1, L1, L2, RTE and R4) and the 41 subfamilies that constitute them. Across these groups we found large variation in the frequency of polymorphic insertions and the observed length distributions of these insertions, suggesting these groups vary in their activity and how frequently they successfully generate full-length, potentially active copies. Though we found an abundance of polymorphic insertions (over 45,000) most of these were observed at low frequencies and typically appeared as singletons. Site frequency spectra for most LINEs showed a significant shift toward low frequency alleles compared to the spectra observed for total genomic single nucleotide polymorphisms. Using Tajima's D, FST and the mean number of pairwise differences in LINE insertion polymorphisms, we found evidence that negative selection is acting on LINE families in a length-dependent manner, its effects being stronger in the larger Eastern Florida population. Our results suggest that a large effective population size and negative selection limit the expansion of polymorphic LINE insertions across these populations and that the probability of LINE polymorphisms reaching fixation is extremely low.

Diversification in wild populations of the model organism Anolis carolinensis: A genome-wide phylogeographic investigation

The green anole (Anolis carolinensis) is a lizard widespread throughout the southeastern United States and is a model organism for the study of reproductive behavior, physiology, neural biology, and genomics. Previous phylogeographic studies of A. carolinensis using mitochondrial DNA and small numbers of nuclear loci identified conflicting and poorly supported relationships among geographically structured clades; these inconsistencies preclude confident use of A. carolinensis evolutionary history in association with morphological, physiological, or reproductive biology studies among sampling localities and necessitate increased effort to resolve evolutionary relationships among natural populations. Here, we used anchored hybrid enrichment of hundreds of genetic markers across the genome of A. carolinensis and identified five strongly supported phylogeographic groups. Using multiple analyses, we produced a fully resolved species tree, investigated relative support for each lineage across all gene trees, and identified mito‐nuclear discordance when comparing our results to previous studies. We found fixed differences in only one clade—southern Florida restricted to the Everglades region—while most polymorphisms were shared between lineages. The southern Florida group likely diverged from other populations during the Pliocene, with all other diversification during the Pleistocene. Multiple lines of support, including phylogenetic relationships, a latitudinal gradient in genetic diversity, and relatively more stable long‐term population sizes in southern phylogeographic groups, indicate that diversification in A. carolinensis occurred northward from southern Florida.

Comparative Phylogeography of Ethiopian anurans: impact of the Great Rift Valley and Pleistocene climate change

Background

The Ethiopian highlands are a biodiversity hotspot, split by the Great Rift Valley into two distinct systems of plateaus and mountains. The Rift Valley is currently hot and dry and acts as a barrier to gene flow for highland-adapted species. It is however unlikely that the conditions in the Rift were inhospitable to highland species during the entire Pleistocene. To assess the significance of the Ethiopian Rift as a biogeographic barrier as well as the impact Pleistocene climatic changes have had on the evolution of Ethiopian organisms, we performed phylogeographic analyses and developed present and past niche models on seven anuran species with different elevational and ecological preferences.

Results

We found that highland species on the east and the west sides of the Rift are genetically differentiated and have not experienced any detectable gene flow for at least 0.4 my. In contrast, species found at elevations lower than 2500 m do not show any population structure. We also determined that highland species have lower effective population sizes than lowland species, which have experienced a large, yet gradual, demographic expansion, starting approximately half a million year ago.

Conclusions

The pattern we report here is consistent with the increasingly warmer and drier conditions of the Pleistocene in East Africa, which resulted in the expansion of savanna, the fragmentation of forests and the shrinking of highland habitats. Climatic niche models indicated that the Rift is currently non suitable for most of the studied species, but it could have been a more permeable barrier during the Last Glacial Maximum. However, considering the strong genetic structure of highland species, we hypothesize that the barrier mechanisms at the Rift are not only climatic but also topographical.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0774-1) contains supplementary material, which is available to authorized users.

Functional evolution of the OAS1 viral sensor: Insights from old world primates

Infections with viral pathogens impose considerable selective pressure on host defensive genes. Those genes at the forefront, responsible for identifying and binding exogenous molecular viral components, will carry the hallmarks of this struggle. Oligoadenylate synthetase (OAS) enzymes play a major role in the innate defense against a large number of viruses by acting as sensors of viral infections. Following their up-regulation by the interferon pathway, OASs bind viral dsRNA and then signal ribonuclease L (RNase L) to degrade RNA, shutting down viral and host protein synthesis. We have investigated the evolution of OAS1 in twenty-two Old World monkey species. We identified a total of 35 codons with the earmarks of positive selection and we performed a comprehensive analysis of their functional significance using in silico modeling of the OAS1 protein. Subdividing OAS1 into functional domains revealed intense purifying selection in the active domain but significant positive directional selection in the RNA-binding domain (RBD), the region where OAS1 binds viral dsRNA. The modeling analysis revealed a concentration of rapidly evolving residues in one region of the RBD suggestive of the sub-functionalization of different regions of the RBD. This analysis also identified several positively selected residues circumscribing the entry to the active site suggesting adaptive evasion of viral antagonism and/or selection for production of oligoadenylate of different length.

Neutral nuclear variation in Baboons (genus Papio) provides insights into their evolutionary and demographic histories

Baboons (genus Papio) are distributed over most of sub-Saharan Africa and in the southern portion of the Arabian Peninsula. Six distinct morphotypes, with clearly defined geographic distributions, are recognized (the olive, chacma, yellow, Guinea, Kinda, and hamadryas baboons). The evolutionary relationships among baboon forms have long been a controversial issue. Phylogenetic analyses based on mitochondrial DNA sequences revealed that the modern baboon morphotypes are mitochondrially paraphyletic or polyphyletic. The discordance between mitochondrial lineages and morphology is indicative of extensive introgressive hybridization between ancestral baboon populations. To gain insights into the evolutionary relationships among morphotypes and their demographic history, we performed an analysis of nuclear variation in baboons. We sequenced 13 noncoding, putatively neutral, nuclear regions, and scored the presence/absence of 18 polymorphic transposable elements in a sample of 45 baboons belonging to five of the six recognized baboon forms. We found that the chacma baboon is the sister-taxon to all other baboons and the yellow baboon is the sister-taxon to an unresolved northern clade containing the olive, Guinea, and hamadryas baboons. We estimated that the diversification of baboons occurred entirely in the Pleistocene, the earliest split dating ∼1.5 million years ago, and that baboons have experienced relatively large and constant effective population sizes for most of their evolutionary history (∼30,000 to 95,000 individuals).

Lizards and LINEs: selection and demography affect the fate of L1 retrotransposons in the genome of the green anole (Anolis carolinensis)

Autonomous retrotransposons lacking long terminal repeats (LTR) account for much of the variation in genome size and structure among vertebrates. Mammalian genomes contain hundreds of thousands of non-LTR retrotransposon copies, mostly resulting from the amplification of a single clade known as L1. The genomes of teleost fish and squamate reptiles contain a much more diverse array of non-LTR retrotransposon families, whereas copy number is relatively low. The majority of non-LTR retrotransposon insertions in nonmammalian vertebrates also appear to be very recent, suggesting strong purifying selection limits the accumulation of non-LTR retrotransposon copies. It is however unclear whether this turnover model, originally proposed in Drosophila, applies to nonmammalian vertebrates. Here, we studied the population dynamics of L1 in the green anole lizard (Anolis carolinensis). We found that although most L1 elements are recent in this genome, truncated insertions accumulate readily, and many are fixed at both the population and species level. In contrast, full-length L1 insertions are found at lower population frequencies, suggesting that the turnover model only applies to longer L1 elements in Anolis. We also found that full-length L1 inserts are more likely to be fixed in populations of small effective size, suggesting that the strength of purifying selection against deleterious alleles is highly dependent on host demographic history. Similar mechanisms seem to be controlling the fate of non-LTR retrotransposons in both Anolis and teleostean fish, which suggests that mammals have considerably diverged from the ancestral vertebrate in terms of how they interact with their intragenomic parasites.

Genetic variation in the green anole lizard (Anolis carolinensis) reveals island refugia and a fragmented Florida during the quaternary

The green anole lizard (Anolis carolinensis) is a model organism for behavior and genomics that is native to the southeastern United States. It is currently thought that the ancestors of modern green anoles dispersed to peninsular Florida from Cuba. However, the climatic changes and geological features responsible for the early diversification of A. carolinensis in North America have remained largely unexplored. This is because previous studies (1) differ in their estimates of the divergence times of populations, (2) are based on a single genetic locus or (3) did not test specific hypotheses regarding the geologic and topographic history of Florida. Here we provide a multi-locus study of green anole genetic diversity and find that the Florida peninsula contains a larger number of genetically distinct populations that are more diverse than those on the continental mainland. As a test of the island refugia hypothesis in Pleistocene Florida, we use a coalescent approach to estimate the divergence times of modern green anole lineages. We find that all demographic events occurred during or after the Upper Pliocene and suggest that green anole diversification was driven by population divergence on interglacial island refugia in Florida during the Lower Pleistocene, while the region was often separated from continental North America. When Florida reconnected to the mainland, two separate dispersal events led to the expansion of green anole populations across the Atlantic Seaboard and Gulf Coastal Plain.

Limited asymmetry dependence of correlations from single nucleon transfer

Single nucleon pickup reactions were performed with a 18.1 MeV/nucleon (14)O beam on a deuterium target. Within the coupled reaction channel framework, the measured cross sections were compared to theoretical predictions and analyzed using both phenomenological and microscopic overlap functions. The missing strength due to correlations does not show significant dependence on the nucleon separation energy asymmetry over a wide range of 37 MeV, in contrast with nucleon removal data analyzed within the sudden-eikonal formalism.

Revisiting the evolution of mouse LINE-1 in the genomic era

Background

LINE-1 (L1) is the dominant category of transposable elements in placental mammals. L1 has significantly affected the size and structure of all mammalian genomes and understanding the nature of the interactions between L1 and its mammalian host remains a question of crucial importance in comparative genomics. For this reason, much attention has been dedicated to the evolution of L1. Among the most studied elements is the mouse L1 which has been the subject of a number of studies in the 1980s and 1990s. These seminal studies, performed in the pre-genomic era when only a limited number of L1 sequences were available, have significantly improved our understanding of L1 evolution. Yet, no comprehensive study on the evolution of L1 in mouse has been performed since the completion of this genome sequence.

Results

Using the Genome Parsing Suite we performed the first evolutionary analysis of mouse L1 over the entire length of the element. This analysis indicates that the mouse L1 has recruited novel 5’UTR sequences more frequently than previously thought and that the simultaneous activity of non-homologous promoters seems to be one of the conditions for the co-existence of multiple L1 families or lineages. In addition the exchange of genetic information between L1 families is not limited to the 5’UTR as evidence of inter-family recombination was observed in ORF1, ORF2, and the 3’UTR. In contrast to the human L1, there was little evidence of rapid amino-acid replacement in the coiled-coil of ORF1, although this region is structurally unstable. We propose that the structural instability of the coiled-coil domain might be adaptive and that structural changes in this region are selectively equivalent to the rapid evolution at the amino-acid level reported in the human lineage.

Conclusions

The pattern of evolution of L1 in mouse shows some similarity with human suggesting that the nature of the interactions between L1 and its host might be similar in these two species. Yet, some notable differences, particularly in the evolution of ORF1, suggest that the molecular mechanisms involved in host-L1 interactions might be different in these two species.

The evolutionary dynamics of transposable elements in eukaryote genomes

Transposable elements (TEs) are ubiquitous components of eukaryotic genomes. They have considerably affected their size, structure and function. The sequencing of a multitude of eukaryote genomes has revealed some striking differences in the abundance and diversity of TEs among eukaryotes. Protists, plants, insects and vertebrates contain species with large numbers of TEs and species with small numbers, as well as species with diverse repertoires of TEs and species with a limited diversity of TEs. There is no apparent relationship between the complexity of organisms and their TE profile. The profile of TE diversity and abundance results from the interaction between the rate of transposition, the intensity of selection against new inserts, the demographic history of populations and the rate of DNA loss. Recent population genetics studies suggest that selection against new insertions, mostly caused by the ability of TEs to mediate ectopic recombination events, is limiting the fixation of TEs, but that reduction in effective population size, caused by population bottlenecks or inbreeding, significantly reduces the efficacy of selection. These results emphasize the importance of drift in shaping genomic architecture.

Multi-locus phylogeographic and population genetic analysis of Anolis carolinensis: historical demography of a genomic model species

The green anole (Anolis carolinensis) has been widely used as an animal model in physiology and neurobiology but has recently emerged as an important genomic model. The recent sequencing of its genome has shed new light on the evolution of vertebrate genomes and on the process that govern species diversification. Surprisingly, the patterns of genetic diversity within natural populations of this widespread and abundant North American lizard remain relatively unknown. In the present study, we use 10 novel nuclear DNA sequence loci (N = 62 to 152) and one mitochondrial locus (N = 226) to delimit green anole populations and infer their historical demography. We uncovered four evolutionarily distinct and geographically restricted lineages of green anoles using phylogenetics, bayesian clustering, and genetic distance methods. Molecular dating indicates that these lineages last shared a common ancestor ∼2 million years ago. Summary statistics and analysis of the frequency distributions of DNA polymorphisms strongly suggest range-wide expansions in population size. Using Bayesian Skyline Plots, we inferred the timing of population size expansions, which differ across lineages, and found evidence for a relatively recent and rapid westward expansion of green anoles across the Gulf Coastal Plain during the mid-Pleistocene. One surprising result is that the distribution of genetic diversity is not consistent with a latitudinal shift caused by climatic oscillations as is observed for many co-distributed taxa. This suggests that the most recent Pleistocene glacial cycles had a limited impact on the geographic distribution of the green anole at the northern limits of its range.

Accumulation and rapid decay of non-LTR retrotransposons in the genome of the three-spine stickleback

The diversity and abundance of non–long terminal repeat (LTR) retrotransposons (nLTR-RT) differ drastically among vertebrate genomes. At one extreme, the genome of placental mammals is littered with hundreds of thousands of copies resulting from the activity of a single clade of nLTR-RT, the L1 clade. In contrast, fish genomes contain a much more diverse repertoire of nLTR-RT, represented by numerous active clades and families. Yet, the number of nLTR-RT copies in teleostean fish is two orders of magnitude smaller than in mammals. The vast majority of insertions appear to be very recent, suggesting that nLTR-RT do not accumulate in fish genomes. This pattern had previously been explained by a high rate of turnover, in which the insertion of new elements is offset by the selective loss of deleterious inserts. The turnover model was proposed because of the similarity between fish and Drosophila genomes with regard to their nLTR-RT profile. However, it is unclear if this model applies to fish. In fact, a previous study performed on the puffer fish suggested that transposable element insertions behave as neutral alleles. Here we examined the dynamics of amplification of nLTR-RT in the three-spine stickleback (Gasterosteus aculeatus). In this species, the vast majority of nLTR-RT insertions are relatively young, as suggested by their low level of divergence. Contrary to expectations, a majority of these insertions are fixed in lake and oceanic populations; thus, nLTR-RT do indeed accumulate in the genome of their fish host. This is not to say that nLTR-RTs are fully neutral, as the lack of fixed long elements in this genome suggests a deleterious effect related to their length. This analysis does not support the turnover model and strongly suggests that a much higher rate of DNA loss in fish than in mammals is responsible for the relatively small number of nLTR-RT copies and for the scarcity of ancient elements in fish genomes. We further demonstrate that nLTR-RT decay in fish occurs mostly through large deletions and not by the accumulation of small deletions.

Long-term balancing selection at the antiviral gene OAS1 in Central African chimpanzees

Oligoadenylate synthetases (OAS) are interferon-induced enzymes that participate in the first line of defense against a wide range of viral infection in animals. Upon activation by viral double-stranded RNA, OAS synthesizes (2-5) oligoadenylates, which activate RNase L, leading to the nonspecific degradation of cellular and viral RNA. Some association studies in humans suggest that variation at one of the OAS genes, OAS1, could be influencing host susceptibility to viral infection. We assessed the diversity of OAS1 in hominoid primates with a focus on chimpanzees. We found that the OAS1 gene is extremely polymorphic in Central African chimpanzee and exhibits levels of silent and replacement diversity much higher than neutral regions of the chimpanzee genome. This level of variation strongly suggests that balancing selection is acting on OAS1, and indeed, this conclusion was validated by several tests of neutrality. We further demonstrated that balancing selection has been acting at this locus since the split between chimpanzees, humans, and gorillas (~8.6 Ma) and caused the persistence of two deeply divergent allelic lineages in Central African chimpanzees. These two groups of OAS1 alleles differ by a large number of amino acids (a.a.), including several a.a. putatively involved in RNA binding. It is therefore very likely that variation at the OAS1 locus affects the innate immune response of individuals to specific viral infection. Our data strongly suggest that interactions between viral RNA and OAS1 are responsible for the maintenance of ancestral polymorphisms at this locus for at least 13.2 My.

The transposable element profile of the anolis genome: How a lizard can provide insights into the evolution of vertebrate genome size and structure

The recent sequencing of the lizard genome provides a unique opportunity to examine the evolution of vertebrate genomes in a phylogenetic context. The lizard genome contains an extraordinary diversity of active transposable elements that far exceeds the diversity reported in extant mammals and birds. Retrotransposons and DNA transposons are represented by multiple active families, contributing to the very diverse repetitive landscape of the lizard. Surprisingly, ancient transposon copies are relatively rare suggesting that the transposon copy number is tightly controlled in lizard. This bias in favor of young copies results from the joint effect of purifying selection acting on novel insertions and a high rate of DNA loss. Recent analyses have revealed that the repetitive landscape of reptiles differ drastically from other extant amniotes by their diversity but also by the dynamics of amplification of their transposons. Thus, from the point of view of mobile elements, reptile genomes show more similarity to fish and amphibians than to other amniotes.

The evolution and diversity of DNA transposons in the genome of the Lizard Anolis carolinensis

DNA transposons have considerably affected the size and structure of eukaryotic genomes and have been an important source of evolutionary novelties. In vertebrates, DNA transposons are discontinuously distributed due to the frequent extinction and recolonization of these genomes by active elements. We performed a detailed analysis of the DNA transposons in the genome of the lizard Anolis carolinensis, the first non-avian reptile to have its genome sequenced. Elements belonging to six of the previously recognized superfamilies of elements (hAT, Tc1/Mariner, Helitron, PIF/Harbinger, Polinton/Maverick, and Chapaev) were identified. However, only four (hAT, Tc1/Mariner, Helitron, and Chapaev) of these superfamilies have successfully amplified in the anole genome, producing 67 distinct families. The majority (57/67) are nonautonomous and demonstrate an extraordinary diversity of structure, resulting from frequent interelement recombination and incorporation of extraneous DNA sequences. The age distribution of transposon families differs among superfamilies and reveals different dynamics of amplification. Chapaev is the only superfamily to be extinct and is represented only by old copies. The hAT, Tc1/Mariner, and Helitron superfamilies show different pattern of amplification, yet they are predominantly represented by young families, whereas divergent families are exceedingly rare. Although it is likely that some elements, in particular long ones, are subjected to purifying selection and do not reach fixation, the majority of families are neutral and accumulate in the anole genome in large numbers. We propose that the scarcity of old copies in the anole genome results from the rapid decay of elements, caused by a high rate of DNA loss.

Phloem protein partners of Cucurbit aphid borne yellows virus: possible involvement of phloem proteins in virus transmission by aphids

Poleroviruses are phytoviruses strictly transmitted by phloem-feeding aphids in a circulative and nonpropagative mode. During ingestion, aphids sample virions in sieve tubes along with sap. Therefore, any sap protein bound to virions will be acquired by the insects and could potentially be involved in the transmission process. By developing in vitro virus-overlay assays on sap proteins collected from cucumber, we observed that approximately 20 proteins were able to bind to purified particles of Cucurbit aphid borne yellows virus (CABYV). Among them, eight proteins were identified by mass spectrometry. The role of two candidates belonging to the PP2-like family (predominant lectins found in cucurbit sap) in aphid transmission was further pursued by using purified orthologous PP2 proteins from Arabidopsis. Addition of these proteins to the virus suspension in the aphid artificial diet greatly increased virus transmission rate. This shift was correlated with an increase in the number of viral genomes in insect cells and with an increase of virion stability in vitro. Surprisingly, increase of the virus transmission rate was also monitored after addition of unrelated proteins in the aphid diet, suggesting that any soluble protein at sufficiently high concentration in the diet and acquired together with virions could stimulate virus transmission.

Efficient and exact maximum likelihood quantisation of genomic features using dynamic programming

An efficient and exact dynamic programming algorithm is introduced to quantise a continuous random variable into a discrete random variable that maximises the likelihood of the quantised probability distribution for the original continuous random variable. Quantisation is often useful before statistical analysis and modelling of large discrete network models from observations of multiple continuous random variables. The quantisation algorithm is applied to genomic features including the recombination rate distribution across the chromosomes and the non-coding transposable element LINE-1 in the human genome. The association pattern is studied between the recombination rate, obtained by quantisation at genomic locations around LINE-1 elements, and the length groups of LINE-1 elements, also obtained by quantisation on LINE-1 length. The exact and density-preserving quantisation approach provides an alternative superior to the inexact and distance-based univariate iterative k-means clustering algorithm for discretisation.

Independent and parallel lateral transfer of DNA transposons in tetrapod genomes

In animals, the mode of transmission of transposable elements is generally vertical. However, recent studies have suggested that lateral transfer has occurred repeatedly in several distantly related tetrapod lineages, including mammals. Using transposons extracted from the genome of the lizard Anolis carolinensis as probes, we identified four novel families of hAT transposons that share extremely high similarity with elements in other genomes including several mammalian lineages (primates, chiropters, marsupials), one amphibian and one flatworm, the planarian Schmidtea mediterranea. The discontinuous phylogenetic distribution of these hAT families, coupled with very low synonymous divergence between species, strongly suggests that these elements were laterally transferred to these different species. This indicates that the horizontal transfer of DNA transposons in vertebrates might be more common than previously thought. Yet, it appears that the transfer of DNA transposons did not occur randomly as the same genomes have been invaded independently by different, unrelated transposon families whereas others seem to be immune to lateral transfer. This suggests that some organisms might be intrinsically more vulnerable to DNA transposon lateral transfer, possibly because of a weakened defense against transposons or because they have developed mechanisms to tolerate their impact.

Long-term balancing selection at the west nile virus resistance gene, Oas1b, maintains transspecific polymorphisms in the house mouse

Oligoadenylate synthetases (OASs) are interferon-inducible enzymes that participate in the first line of defense against a wide range of viral infection. Recent studies have determined that Oas1b, a member of the OAS gene family in the house mouse (Mus musculus), provides specific protection against flavivirus infection (e.g., West Nile virus, dengue fever virus, and yellow fever virus). We characterized the nucleotide sequence variation in coding and noncoding regions of the Oas1b gene for a large number of wild-derived strains of M. musculus and related species. Our sequence analyses determined that this gene is one of the most polymorphic genes ever described in any mammal. The level of variation in noncoding regions of Oas1b is an order of magnitude higher than the level reported for other regions of the mouse genome and is significantly different from the level of intraspecific variation expected under neutrality. Furthermore, a phylogenetic analysis of intronic sequences demonstrated that Oas1b alleles are ancient and that their divergence predates several speciation events, resulting in transspecific polymorphisms. The amino acid sequence of Oas1b is also extremely variable, with 1 out of 7 amino acid positions being polymorphic within M. musculus. Oas1b alleles are comparatively more divergent at synonymous positions than most autosomal genes and the ratio of nonsynonymous to synonymous substitution is remarkably high, suggesting that positive selection has been acting on Oas1b. The ancestry of Oas1b polymorphisms and the high level of amino acid polymorphisms strongly suggest that the allelic variation at Oas1b has been maintained in mouse populations by long-term balancing selection.

Selection against LINE-1 retrotransposons results principally from their ability to mediate ectopic recombination

LINE-1 (L1) retrotransposons constitute the most successful family of autonomous retroelements in mammals and they represent at least 17% of the size of the human genome. L1 insertions have occasionally been recruited to perform a beneficial function but the vast majority of L1 inserts are either neutral or deleterious. The basis for the deleterious effect of L1 remains a matter of debate and three possible mechanisms have been suggested: the direct effect of L1 inserts on gene activity, genetic rearrangements caused by L1-mediated ectopic recombination, or the retrotransposition process per se. We performed a genome-wide analysis of the distribution of L1 retrotransposons relative to the local recombination rate and the age and length of the elements. The proportion of L1 elements that are longer than 1.2 Kb is higher in low-recombining regions of the genome than in regions with a high recombination rate, but the genomic distributions of full-length elements (i.e. elements capable of retrotransposition) and long truncated elements were indistinguishable. We also found that the intensity of selection against long elements is proportional to the replicative success of L1 families. This suggests that the deleterious effect of L1 elements results principally from their ability to mediate ectopic recombination.

Human population genetic structure and diversity inferred from polymorphic L1(LINE-1) and Alu insertions

BACKGROUND/AIMS

The L1 retrotransposable element family is the most successful self-replicating genomic parasite of the human genome. L1 elements drive replication of Alu elements, and both have had far-reaching impacts on the human genome. We use L1 and Alu insertion polymorphisms to analyze human population structure.

METHODS

We genotyped 75 recent, polymorphic L1 insertions in 317 individuals from 21 populations in sub-Saharan Africa, East Asia, Europe and the Indian subcontinent. This is the first sample of L1 loci large enough to support detailed population genetic inference. We analyzed these data in parallel with a set of 100 polymorphic Alu insertion loci previously genotyped in the same individuals.

RESULTS AND CONCLUSION

The data sets yield congruent results that support the recent African origin model of human ancestry. A genetic clustering algorithm detects clusters of individuals corresponding to continental regions. The number of loci sampled is critical: with fewer than 50 typical loci, structure cannot be reliably discerned in these populations. The inclusion of geographically intermediate populations (from India) reduces the distinctness of clustering. Our results indicate that human genetic variation is neither perfectly correlated with geographic distance (purely clinal) nor independent of distance (purely clustered), but a combination of both: stepped clinal.

Molecular evolution and tempo of amplification of human LINE-1 retrotransposons since the origin of primates

We investigated the evolution of the families of LINE-1 (L1) retrotransposons that have amplified in the human lineage since the origin of primates. We identified two phases in the evolution of L1. From approximately 70 million years ago (Mya) until approximately 40 Mya, three distinct L1 lineages were simultaneously active in the genome of ancestral primates. In contrast, during the last 40 million years (Myr), i.e., during the evolution of anthropoid primates, a single lineage of families has evolved and amplified. We found that novel (i.e., unrelated) regulatory regions (5'UTR) have been frequently recruited during the evolution of L1, whereas the two open-reading frames (ORF1 and ORF2) have remained relatively conserved. We found that L1 families coexisted and formed independently evolving L1 lineages only when they had different 5'UTRs. We propose that L1 families with different 5'UTR can coexist because they don't rely on the same host-encoded factors for their transcription and therefore do not compete with each other. The most prolific L1 families (families L1PA8 to L1PA3) amplified between 40 and 12 Mya. This period of high activity corresponds to an episode of adaptive evolution in a segment of ORF1. The correlation between the high activity of L1 families and adaptive evolution could result from the coevolution of L1 and a host-encoded repressor of L1 activity.

Glycosylation of beet western yellows virus proteins is implicated in the aphid transmission of the virus

Beet western yellows virus relies on the aphid M. persicae for its transmission in a persistent and circulative mode. To be transmitted, the virus must cross the midgut and the accessory salivary gland epithelial barriers by a transcytosis mechanism where vector receptors interact with virions. The aphid and the peptidic viral determinants implicated in this interaction mechanism have been studied. In this paper, we report that the coat and the readthrough proteins that constitute the capsid of this virus are glycosylated. Modification of the glucidic core of these structural viral proteins by oxidation with sodium metaperiodate or deglycosylation with N-glycosidase F or alpha-D-galactosidase abrogates the aphid transmission of the virus. Aphid transmission could also be inhibited by lectins directed against alpha-D-galactose when aphids were allowed to acquire virus on artificial membranes. These results suggest that the glucidic cores of the capsid proteins of beet western yellows virus contain alpha-D-galactose residues that are implicated in virus-aphid interaction and promote aphid transmission of the virus.

The recent evolution of human L1 retrotransposons

L1 elements are the most successful retrotransposons in mammals and are responsible for at least 30% of human DNA. Far from being indolent genomic parasites, L1 elements have evolved and amplified rapidly during human evolution. Indeed during just the last 25 million years (MY) five distinct L1 families have emerged and generated tens of thousands of copies. The most recently evolved human specific L1 family is currently active and L1 copies have been accumulating in the human genome at about the same rate per generation as the currently active L1 families in Old World rats and mice. At times during the last 25 MY L1 activity constituted a significant enough genetic load to be subject to negative selection. During these same times, and in apparent response to the host, L1 underwent adaptive evolution. Understanding the molecular basis for these evolutionary changes should help illuminate one of the least understood but most important aspects of L1 biology, namely the extent and nature of the interaction between L1 and its host.

The structures of mouse and human L1 elements reflect their insertion mechanism

L1 is an abundant, interspersed repeated DNA element of mammalian genomes. It has achieved its high copy number via retrotransposition. Like other non-LTR retrotransposons, L1 insertion into chromosomal DNA apparently occurs by target-site primed reverse transcription, or TPRT. L1 retrotransposition often generates elements with 5' truncations that are flanked by a duplication of the genomic target site (TSD). It is typically assumed that the 5' truncated elements are the consequence of poor processivity of the L1 reverse transcriptase. However, we find that the majority of young L1 elements from both the human and mouse genomes are truncated at sequences that can basepair with the target site. Thus, to whatever extent truncation is a consequence of poor processivity, we suggest that truncation is likely to occur when target site sequence can basepair with L1 sequence. This finding supports a model for insertion that occurs by two sequential TPRT reactions, the second of which relies upon the homology between the target site and L1. Because perfect heteroduplex formation is not required for all insertions, a dynamic relationship between the primer, template and enzyme during reverse transcription is inferred. 5' truncation may be a successful evolutionary strategy that is exploited by L1 as a means to escape host suppression of transposition.