Divergence and introgression among the virilis group of Drosophila

Speciation with gene flow is now widely regarded as common. However, the frequency of introgression between recently diverged species and the evolutionary consequences of gene flow are still poorly understood. The virilis group of Drosophila contains around a dozen species that are geographically widespread and show varying levels of pre-zygotic and post-zygotic isolation. Here, we utilize de novo genome assemblies and whole-genome sequencing data to resolve phylogenetic relationships and describe patterns of introgression and divergence across the group. We suggest that the virilis group consists of three, rather than the traditional two, subgroups. We found evidence of pervasive phylogenetic discordance caused by ancient introgression events between distant lineages within the group, and much more recent gene flow between closely-related species. When assessing patterns of genome-wide divergence in species pairs across the group, we found no consistent genomic evidence of a disproportionate role for the X chromosome. Some genes undergoing rapid sequence divergence across the group were involved in chemical communication and may be related to the evolution of sexual isolation. We suggest that gene flow between closely-related species has potentially had an impact on lineage-specific adaptation and the evolution of reproductive barriers. Our results show how ancient and recent introgression confuse phylogenetic reconstruction, and suggest that shared variation can facilitate adaptation and speciation.

The megabase-scale crossover landscape is independent of sequence divergence

Meiotic recombination frequency varies along chromosomes and strongly correlates with sequence divergence. However, the causality underlying this correlation is unclear. To untangle the relationship between recombination landscapes and polymorphisms, we characterized the genome-wide recombination landscape in the absence of polymorphisms, using Arabidopsis thaliana homozygous inbred lines in which a few hundred genetic markers were introduced through mutagenesis. We found that megabase-scale recombination landscapes in inbred lines are strikingly similar to the recombination landscapes in hybrids, with the sole exception of heterozygous large rearrangements where recombination is prevented locally. In addition, we found that the megabase-scale recombination landscape can be accurately predicted by chromatin features. Our results show that polymorphisms are not causal for the shape of the megabase-scale recombination landscape, rather, favor alternative models in which recombination and chromatin shape sequence divergence across the genome.

Sequence Divergence and Functional Specializations of the Ancient Spliceosomal SF3b: Implications in Flexibility and Adaptations of the Multi-Protein Complex

Multi-protein assemblies are complex molecular systems that perform highly sophisticated biochemical functions in an orchestrated manner. They are subject to changes that are governed by the evolution of individual components. We performed a comparative analysis of the ancient and functionally conserved spliceosomal SF3b complex, to recognize molecular signatures that contribute to sequence divergence and functional specializations. For this, we recognized homologous sequences of individual SF3b proteins distributed across 10 supergroups of eukaryotes and identified all seven protein components of the complex in 578 eukaryotic species. Using sequence and structural analysis, we establish that proteins occurring on the surface of the SF3b complex harbor more sequence variation than the proteins that lie in the core. Further, we show through protein interface conservation patterns that the extent of conservation varies considerably between interacting partners. When we analyze phylogenetic distributions of individual components of the complex, we find that protein partners that are known to form independent subcomplexes are observed to share similar profiles, reaffirming the link between differential conservation of interface regions and their inter-dependence. When we extend our analysis to individual protein components of the complex, we find taxa-specific variability in molecular signatures of the proteins. These trends are discussed in the context of proline-rich motifs of SF3b4, functional and drug binding sites of SF3b1. Further, we report key protein-protein interactions between SF3b1 and SF3b6 whose presence is observed to be lineage-specific across eukaryotes. Together, our studies show the association of protein location within the complex and subcomplex formation patterns with the sequence conservation of SF3b proteins. In addition, our study underscores evolutionarily flexible elements that appear to confer adaptive features in individual components of the multi-protein SF3b complexes and may contribute to its functional adaptability.

A new rainfrog of the genus Pristimantis (Anura, Brachycephaloidea) from central and eastern Panama

Substantial molecular and morphological character differences lead us to the description of a new species of the genus Pristimantis from the cloud forest of Cerro Chucantí, Maje Mountains, Darien Province, as well as from several other mountain ranges in eastern and central Panama. Pristimantis gretathunbergaesp. nov. is a sister species to the allopatric P. erythropleura-penelopus group from northern Colombia with a mtDNA sequence divergence of > 4.4% at 16S and > 14.6% at COI. Its closest congener in sympatry is P. cruentus that differs by a large sequence divergence of > 9.6% in 16S mtDNA and 19.0% at COI, and from which it differs also by ventral and groin coloration, unusually prominent black eyes, a contrasting light upper lip, commonly a single conical to spine-like tubercle on the upper eyelid, and a larger head. While the habitat continuity at most sites in eastern Panama is moderate, habitats in central Panama are severely fragmented. Cerro Chucantí and the surrounding Maje Mountains are highly threatened by rapid deforestation and replaced by plantations and cattle pastures. Thus, investigations on the ecology of the new species and its population status, especially at the type locality, are highly recommended. As a flagship species, this new frog can help to preserve the Chucantí cloud forest including several recently described species known only from this isolated area in eastern Panama.

Fragmentation in mitochondrial genomes in relation to elevated sequence divergence and extreme rearrangements

Background A single circular mitochondrial (mt) genome is a common feature across most metazoans. The mt-genome includes protein-coding genes involved in oxidative phosphorylation, as well as RNAs necessary for translation of mt-RNAs, whose order and number are highly conserved across animal clades, with few known exceptions of alternative mt-gene order or mt-genome architectures. One such exception consists of the fragmented mitochondrial genome, a type of genome architecture where mt-genes are split across two or more mt-chromosomes. However, the origins of mt-genome fragmentation and its effects on mt-genome evolution are unknown. Here, we investigate these origin and potential mechanisms underlying mt-genome fragmentation, focusing on a genus of booklice, Liposcelis, which exhibits elevated sequence divergence, frequent rearrangement of mt-gene order, and fragmentation of the mt genome, and compare them to other Metazoan clades. Results We found this genus Liposcelis exhibits very low conservation of mt-gene order across species, relative to other metazoans. Levels of gene order rearrangement were, however, unrelated to whether or not mt-genomes were fragmented or intact, suggesting mitochondrial genome fragmentation is not affecting mt-gene order directly. We further investigated possible mechanisms underpinning these patterns and revealed very high conservation of non-coding sequences at the edges of multiple recombination regions across populations of one particular Liposcelis species, supportive of a hypothesis that mt-fragmentation arises from recombination errors between mt-genome copies. We propose these errors may arise as a consequence of a heightened mutation rate in clades exhibiting mt-fragmentation. Consistent with this, we observed a striking pattern across three Metazoan phyla (Arthropoda, Nematoda, Cnidaria) characterised by members exhibiting high levels of mt-gene order rearrangement and cases of mt-fragmentation, whereby the mt-genomes of species more closely related to species with fragmented mt-genomes diverge more rapidly despite experiencing strong purifying selection. Conclusions We showed that contrary to expectations, mt-genome fragmentation is not correlated with the increase in mt-genome rearrangements. Furthermore, we present evidence that fragmentation of the mt-genome may be part of a general relaxation of a natural selection on the mt-genome, thus providing new insights into the origins of mt-genome fragmentation and evolution.

The complete chloroplast genome sequence of Quercus ningangensis and its phylogenetic implication

Quercus ningangensis is an economically and ecologically important tree species belonging to the family Fagaceae. In this study, the complete chloroplast (cp) genome of Q. ningangensis was sequenced and assembled, and 18 published cp genomes of Quercus were retrieved for genomic analyses (including sequence divergence, repeat elements, and structure) and phylogenetic inference. With this study, we found that complete cp genomes in Quercus are conserved, and we discovered a codon composition bias, which may be related to genomic content and genetic characteristics. In addition, we detected considerable structural variations in the expansion and contraction of inverted repeat regions. Six regions with relatively high variable (matK-rps16, psbC, ycf3 intron, rbcL, petA-psbJ, and ycf1) were detected by conducting a sliding window analysis, which has a high potential for developing effective genetic markers. Phylogenetic analysis based on Bayesian inference and maximum likelihood methods resulted in a robust phylogenetic tree of Quercus with high resolution for nearly all identified nodes. The phylogenetic relationships showed that the phylogenetic position of Q. ningangensis was located between Q. sichourensis and Q. acuta. The results of this study contribute to future research into the phylogenetic evolution of Quercus section Cyclobalanopsis (Fagaceae).

Evolutionary pathways to SARS-CoV-2 resistance are opened and closed by epistasis acting on ACE2

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infects a broader range of mammalian species than previously predicted, binding a diversity of angiotensin converting enzyme 2 (ACE2) orthologs despite extensive sequence divergence. Within this sequence degeneracy, we identify a rare sequence combination capable of conferring SARS-CoV-2 resistance. We demonstrate that this sequence was likely unattainable during human evolution due to deleterious effects on ACE2 carboxypeptidase activity, which has vasodilatory and cardioprotective functions in vivo. Across the 25 ACE2 sites implicated in viral binding, we identify 6 amino acid substitutions unique to mouse—one of the only known mammalian species immune to SARS-CoV-2. Substituting human variants at these positions is sufficient to confer binding of the SARS-CoV-2 S protein to mouse ACE2, facilitating cellular infection. Conversely, substituting mouse variants into either human or dog ACE2 abolishes viral binding, diminishing cellular infection. However, these same substitutions decrease human ACE2 activity by 50% and are predicted as pathogenic, consistent with the extreme rarity of human polymorphisms at these sites. This trade-off can be avoided, however, depending on genetic background; if substituted simultaneously, these same mutations have no deleterious effect on dog ACE2 nor that of the rodent ancestor estimated to exist 70 million years ago. This genetic contingency (epistasis) may have therefore opened the road to resistance for some species, while making humans susceptible to viruses that use these ACE2 surfaces for binding, as does SARS-CoV-2.

Transposable Prophages in Leptospira: An Ancient, Now Diverse, Group Predominant in Causative Agents of Weil’s Disease

The virome associated with the corkscrew shaped bacterium Leptospira, responsible for Weil’s disease, is scarcely known, and genetic tools available for these bacteria remain limited. To reduce these two issues, potential transposable prophages were searched in Leptospiraceae genomes. The 236 predicted transposable prophages were particularly abundant in the most pathogenic leptospiral clade, being potentially involved in the acquisition of virulent traits. According to genomic similarities and phylogenies, these prophages are distantly related to known transposable phages and are organized into six groups, one of them encompassing prophages with unusual TA-TA ends. Interestingly, structural and transposition proteins reconstruct different relationships between groups, suggesting ancestral recombinations. Based on the baseplate phylogeny, two large clades emerge, with specific gene-contents and high sequence divergence reflecting their ancient origin. Despite their high divergence, the size and overall genomic organization of all prophages are very conserved, a testimony to the highly constrained nature of their genomes. Finally, similarities between these prophages and the three known non-transposable phages infecting L. biflexa, suggest gene transfer between different Caudovirales inside their leptospiral host, and the possibility to use some of the transposable prophages in that model strain.

Identification of Novel Phage Resistance Mechanisms in Campylobacter jejuni by Comparative Genomics

Phages infecting Campylobacter jejuni are considered a promising intervention strategy at broiler farms, yet phage sensitivity of naturally occurring poultry isolates is not well studied. Here, we investigated phage sensitivity and identified resistance mechanisms of C. jejuni strains originating from Danish broilers belonging to the most prevalent MLST (ST) types. Determining plaque formation of 51 phages belonging to Fletchervirus or Firehammervirus showed that 21 out of 31 C. jejuni strains were susceptible to at least one phage. While C. jejuni ST-21 strains encoded the common phase variable O-methyl phosphoramidate (MeOPN) receptor of the Fletchervirus and were only infected by these phages, ST-45 strains did not encode this receptor and were exclusively infected by Firehammervirus phages. To identify internal phage resistance mechanism in ST-21 strains, we performed comparative genomics of two strains, CAMSA2002 sensitive to almost all Fletchervirus phages and CAMSA2038, resistant to all 51 phages. The strains encoded diverse clustered regularly interspaced short palindromic repeats (CRISPR) spacers but none matched the tested phages. Sequence divergence was also observed in a predicted SspE homolog and putative restriction modification systems including a methyl-specific McrBC endonuclease. Furthermore, when mcrB was deleted, CAMSA2038 became sensitive to 17 out of 43 phages, three being Firehammervirus phages that otherwise did not infect any ST-21 strains. Yet, 16 phages demonstrated significantly lower efficiencies of plating on the mcrB mutant suggesting additional resistance mechanism still restricting phage propagation in CAMSA2038. Thus, our work demonstrates that C. jejuni isolates originating from broilers may have acquired several resistance mechanisms to successfully prevent phage infection in their natural habitat.