Search

International audience; The distribution of fitness effects (DFE) of new mutations plays a central role in molecular evolution. It is therefore crucial to be able to estimate it accurately from genomic data and to understand the factors that shape it. After a rapid overview of available methods to characterize the fitness effects of mutations, we review what is known on the factors affecting them in plants. Available data indicate that life history traits, e.g. mating system and longevity, have a major effect on the DFE. In contrast, the impact of demography within species appears to be more limited. These results remain to be confirmed and methods to estimate the joint evolution of demography, life history traits and the DFE to be developed.

Organisms use color for camouflage, sexual signaling, or as a warning sign of danger. Primates are one of the most vibrantly colored Orders of mammals. However, the genetics underlying their coat color are poorly known, limiting our ability to study molecular aspects of its evolution. The role of the melanocortin 1 receptor (MC1R) in color evolution has been implicated in studies on rocket pocket mice (Chaetodipus intermediusi), toucans (Ramphastidae), and many domesticated animals. From these studies, we know that changes in MC1R result in a yellow/red or a brown/black morphology. Here, we investigate the evolution of MC1R in Lorisidae, a monophyletic nocturnal primate family, with some genera displaying high contrast variation in color patterns and other genera being monochromatic. Even more unique, the Lorisidae family has the only venomous primate: the slow loris (Nycticebus). Research has suggested that the contrasting coat patterns of slow lorises are aposematic signals for their venom. If so, we predict the MC1R in slow lorises will be under positive selection. In our study, we found that Lorisidae MC1R is under purifying selection (ω = 0.0912). In Lorisidae MC1R, there were a total of 75 variable nucleotides, 18 of which were nonsynonymous. Six of these nonsynonymous substitutions were found on the Perodicticus branch, which our reconstructions found to be the only member of Lorisidae that has predominantly lighter coat color; no substitutions were associated with Nycticebus. Our findings generate new insight into the genetics of pelage color and evolution among a unique group of nocturnal mammals and suggest putative underpinnings of monochromatic color evolution in the Perodicticus lineage., In primates, the most colorful of mammals, few have examined the influence the melanocortin 1 receptor (MC1R) has on color variation (light or dark) in primates. Here, we examined how MC1R impacts color variation in a cryptic primate family with high contrasting color variation. We found MC1R is under purifying selection, but there appears to be a link in nucleotide changes and light coat color.

DNA barcoding based on mitochondrial (mt) nucleotide sequences is an enigma. Neutral models of mt evolution predict DNA barcoding cannot work for recently diverged taxa, and yet, mt DNA barcoding accurately delimits species for many bilaterian animals. Meanwhile, mt DNA barcoding often fails for plants and fungi. I propose that because mt gene products must cofunction with nuclear gene products, the evolution of mt genomes is best understood with full consideration of the two environments that impose selective pressure on mt genes: the external environment and the internal genomic environment. Moreover, it is critical to fully consider the potential for adaptive evolution of not just protein products of mt genes but also of mt transfer RNAs and mt ribosomal RNAs. The tight linkage of genes on mt genomes that do not engage in recombination could facilitate selective sweeps whenever there is positive selection on any element in the mt genome, leading to the purging of mt genetic diversity within a population and to the rapid fixation of novel mt DNA sequences. Accordingly, the most important factor determining whether or not mt DNA sequences diagnose species boundaries may be the extent to which the mt chromosomes engage in recombination., Neutral models of mitochondrial evolution predict mitochondrial DNA barcoding cannot work for recently diverged taxa, but mt DNA barcoding accurately delimits species across a wide spectrum of bilaterian animals. I propose that the tight linkage of genes on mt genomes that do not engage in recombination could facilitate selective sweeps whenever there is positive selection on any element in the mt genome, leading to the purging of mt genetic diversity within a population and to the rapid fixation of novel mt DNA sequences. If selective sweeps play a key role in generating mt DNA barcode gaps, then the most important factor determining whether or not mt DNA sequences diagnose species boundaries within a taxon would be the extent to which the mt chromosomes of individuals within that taxon engage in recombination.

The endosperm cell walls of mature coffee seeds accumulate large amounts of mannan storage polysaccharides, which serve as nutrient reserve for embryo and contribute to beverage quality. Our study investigated the evolutionary patterns of key galactomannan (GM) biosynthesis genes using d N/d S ratio, synteny, and phylogenetic analysis and detected heterogeneity in rate of evolution among gene copies. Selection ratio index revealed evidence of positive selection in the branch editing gene Coffea canephora alpha (α) galactosidase (Cc‐alpha Gal) at Cc11_g15950 copy (ω = 1.12), whereas strong purifying selection on deleterious mutations was observed in the Coffea canephora uridine diphosphate (UDP)‐glucose 4′‐epimerase (Cc‐UG4E) and Coffea canephora mannose‐1P guanylytransferase (Cc‐MGT) genes controlling the crucial nucleotide carbon sugar building blocks flux in the pathway. Relatively low sequence diversity and strong syntenic linkages were detected in all GM pathway genes except in Cc‐alpha Gal, which suggests a correlation between selection pressure and nucleotide diversity or synteny analysis. In addition, phylogenetic analysis revealed independent evolution or expansion of GM pathway genes in different plant species, with no obvious inferable clustering patterns according to either gene family or congruent with evolutionary plants lineages tested due to high dynamic nature and specific biochemical cell wall modification requirements. Altogether, our study shows a significant high rate of evolutionary variation among GM pathway genes in the diploid C. canephora and demonstrates the inherent variation in evolution of gene copies and their potential role in understanding selection rates in a homogenously connected metabolic pathway., This study investigated the evolutionary patterns of key galactomannan biosynthesis genes using selection ratio index, synteny, and phylogenetic analysis and detected heterogeneity in rate of evolution among gene copies.

UDP- glucose: glycoprotein glucosyltransferase (UGGT) is a protein that operates as the gatekeeper for the endoplasmic reticulum (ER) quality control mechanism of glycoprotein folding. It is known that vertebrates and Caenorhabditis genomes harbor two uggt gene copies that exhibit differences in their properties.Bayesian phylogenetic inference based on 195 UGGT and UGGT-like protein sequences of an ample spectrum of eukaryotic species showed that uggt genes went through independent duplications in Caenorhabditis and vertebrates. In both lineages, the catalytic domain of the duplicated genes was subjected to a strong purifying selective pressure, while the recognition domain was subjected to episodic positive diversifying selection. Selective relaxation in the recognition domain was more pronounced in Caenorhabditis uggt-b than in vertebrates uggt-2. Structural bioinformatics analysis revealed that Caenorhabditis UGGT-b protein lacks essential sequences proposed to be involved in the recognition of unfolded proteins. When we assayed glucosyltrasferase activity of a chimeric protein composed by Caenorhabditis uggt-b recognition domain fused to S. pombe catalytic domain expressed in yeast, no activity was detected.The present results support the conservation of the UGGT activity in the catalytic domain and a putative divergent function of the recognition domain for the UGGT2 protein in vertebrates, which would have gone through a specialization process. In Caenorhabditis, uggt-b evolved under different constraints compared to uggt-a which, by means of a putative neofunctionalization process, resulted in a non-redundant paralog. The non-canonical function of uggt-b in the worm lineage highlights the need to take precautions before generalizing gene functions in model organisms.

Recombination's omnipresence in nature is one of the most intriguing problems in evolutionary biology. The question of why recombination exhibits certain general features is no less interesting than that of why it exists at all. One such feature is recombination's fitness dependence (FD). The so far developed population genetics models have focused on the evolution of FD recombination mainly in haploids, although the empirical evidence for this phenomenon comes mostly from diploids. Using numerical analysis of modifier models for infinite panmictic populations, we show here that FD recombination can be evolutionarily advantageous in diploids subjected to purifying selection. We ascribe this advantage to the differential rate of disruption of lower‐ versus higher‐fitness genotypes, which can be manifested in selected systems with at least three loci. We also show that if the modifier is linked to such selected system, it can additionally benefit from modifying this linkage in a fitness‐dependent manner. The revealed evolutionary advantage of FD recombination appeared robust to crossover interference within the selected system, either positive or negative. Remarkably, FD recombination was often favored in situations where any constant nonzero recombination was evolutionarily disfavored, implying a relaxation of the rather strict constraints on major parameters (e.g., selection intensity and epistasis) required for the evolutionary advantage of nonzero recombination formulated by classical models., An amazing feature of recombination is fitness dependence (FD). The so far developed models have focused on the evolution of FD recombination mainly in haploids, although the empirical evidence for this phenomenon comes mostly from diploids. Using numerical analysis of modifier models for infinite panmictic populations, we show here that FD recombination can be evolutionarily advantageous in diploids subjected to purifying selection. We ascribe this advantage to the differential rate of disruption of lower‐ versus higher‐fitness genotypes, which can be manifested in selected systems with at least three loci.

Owing to a lag between a deleterious mutation’s appearance and its selective removal, gold-standard methods for mutation rate estimation assume no meaningful loss of mutations between parents and offspring. Indeed, from analysis of closely related lineages, in SARS-CoV-2, the Ka/Ks ratio was previously estimated as 1.008, suggesting no within-host selection. By contrast, we find a higher number of observed SNPs at 4-fold degenerate sites than elsewhere and, allowing for the virus’s complex mutational and compositional biases, estimate that the mutation rate is at least 49–67% higher than would be estimated based on the rate of appearance of variants in sampled genomes. Given the high Ka/Ks one might assume that the majority of such intrahost selection is the purging of nonsense mutations. However, we estimate that selection against nonsense mutations accounts for only ∼10% of all the “missing” mutations. Instead, classical protein-level selective filters (against chemically disparate amino acids and those predicted to disrupt protein functionality) account for many missing mutations. It is less obvious why for an intracellular parasite, amino acid cost parameters, notably amino acid decay rate, is also significant. Perhaps most surprisingly, we also find evidence for real-time selection against synonymous mutations that move codon usage away from that of humans. We conclude that there is common intrahost selection on SARS-CoV-2 that acts on nonsense, missense, and possibly synonymous mutations. This has implications for methods of mutation rate estimation, for determining times to common ancestry and the potential for intrahost evolution including vaccine escape.

DNA diversity varies across the genome of many species. Variation in diversity across a genome might arise from regional variation in the mutation rate, variation in the intensity and mode of natural selection, and regional variation in the recombination rate. We show that both noncoding and nonsynonymous diversity are positively correlated to a measure of the mutation rate and the recombination rate and negatively correlated to the density of conserved sequences in 50 kb windows across the genomes of humans and nonhuman homininae. Interestingly, we find that although noncoding diversity is equally affected by these three genomic variables, nonsynonymous diversity is mostly dominated by the density of conserved sequences. The positive correlation between diversity and our measure of the mutation rate seems to be largely a direct consequence of regions with higher mutation rates having more diversity. However, the positive correlation with recombination rate and the negative correlation with the density of conserved sequences suggest that selection at linked sites also affect levels of diversity. This is supported by the observation that the ratio of the number of nonsynonymous to noncoding polymorphisms is negatively correlated to a measure of the effective population size across the genome. We show these patterns persist even when we restrict our analysis to GC-conservative mutations, demonstrating that the patterns are not driven by GC biased gene conversion. In conclusion, our comparative analyses describe how recombination rate, gene density, and mutation rate interact to produce the patterns of DNA diversity that we observe along the hominine genomes. This study was supported by the Danish Council For Independent Research (grant number 4181-00358

Background The monogenean Benedenia seriolae parasitizes fishes belonging to the genus Seriola, represents a species complex, and causes substantial impact on fish welfare in aquaculture systems worldwide. This study reports, for the first time, the complete mitochondrial genome of B. humboldti n. sp., a new cryptic species from the South-East Pacific (SEP). Methods The mitogenome of B. humboldti n. sp. was assembled from short Illumina 150 bp pair-end reads. The phylogenetic position of B. humboldti n. sp. among other closely related congeneric and confamiliar capsalids was examined using mitochondrial protein-coding genes (PCGs). Morphology of B. humboldti n. sp. was examined based on fixed and stained specimens. Results The AT-rich mitochondrial genome of B. humboldti is 13,455 bp in length and comprises 12 PCGs (atp8 was absent as in other monogenean genomes), 2 ribosomal RNA genes, and 22 transfer RNA genes. All protein-coding, ribosomal RNA, and transfer RNA genes are encoded on the H-strand. The gene order observed in the mitochondrial genome of B. humboldti n. sp. was identical to that of B. seriolae from Japan but different from that of B. seriolae from Australia. The genetic distance between B. humboldti n. sp. and B. seriolae from Japan was high. Minor but reliable differences in the shape of the penis were observed between Benedenia humboldti n. sp. and congeneric species. Conclusions Phylogenetic analyses based on PCGs in association with differences in the shape of the penis permitted us to conclude that the material from the South-East Pacific represents a new species of Benedenia infecting S. lalandi off the coast of Chile. The discovery of this parasite represents the first step to improving our understanding of infestation dynamics and to develop control strategies for this pathogen infecting the farmed yellowtail kingfish, Seriola lalandi, in the South-East Pacific.

Marine microorganisms inhabiting nutrient-depleted waters play critical roles in global biogeochemical cycles due to their abundance and broad distribution. Many of these microbes share similar genomic features including small genome size, low % G + C content, short intergenic regions, and low nitrogen content in encoded amino acid residue side chains (N-ARSC), but the evolutionary drivers of these characteristics are unclear. Here, we compared the strength of purifying selection across the Marinimicrobia, a candidate phylum which encompasses a broad range of phylogenetic groups with disparate genomic features, by estimating the ratio of nonsynonymous and synonymous substitutions (dN/dS) in conserved marker genes. Our analysis reveals that epipelagic Marinimicrobia that exhibit features consistent with genome streamlining have significantly lower dN/dS values when compared with their mesopelagic counterparts. We also found a significant positive correlation between median dN/dS values and % G + C content, N-ARSC, and intergenic region length. We did not identify a significant correlation between dN/dS ratios and estimated genome size, suggesting the strength of selection is not a primary factor shaping genome size in this group. Our findings are generally consistent with genome streamlining theory, which postulates that many genomic features of abundant epipelagic bacteria are the result of adaptation to oligotrophic nutrient conditions. Our results are also in agreement with previous findings that genome streamlining is common in epipelagic waters, suggesting that microbes inhabiting this region of the ocean have been shaped by strong selection together with prevalent nutritional constraints characteristic of this environment.

In species with chromosomal sex determination, X chromosomes are predicted to evolve faster than autosomes because of positive selection on recessive alleles or weak purifying selection. We investigated X chromosome evolution in Stegodyphus spiders that differ in mating system, sex ratio, and population dynamics. We assigned scaffolds to X chromosomes and autosomes using a novel method based on flow cytometry of sperm cells and reduced representation sequencing. We estimated coding substitution patterns (dN/dS) in a subsocial outcrossing species (S. africanus) and its social inbreeding and female-biased sister species (S. mimosarum), and found evidence for faster-X evolution in both species. X chromosome-to-autosome diversity (piX/piA) ratios were estimated in multiple populations. The average piX/piA estimates of S. africanus (0.57 [95% CI: 0.55–0.60]) was lower than the neutral expectation of 0.75, consistent with more hitchhiking events on X-linked loci and/or a lower X chromosome mutation rate, and we provide evidence in support of both. The social species S. mimosarum has a significantly higher piX/piA ratio (0.72 [95% CI: 0.65–0.79]) in agreement with its female-biased sex ratio. Stegodyphus mimosarum also have different piX/piA estimates among populations, which we interpret as evidence for recurrent founder events. Simulations show that recurrent founder events are expected to decrease the piX/piA estimates in S. mimosarum, thus underestimating the true effect of female-biased sex ratios. Finally, we found lower synonymous divergence on X chromosomes in both species, and the male-to-female substitution ratio to be higher than 1, indicating a higher mutation rate in males.

Background Animal mitochondria play a central role in energy production in the cells through the oxidative phosphorylation (OXPHOS) pathway. Recent studies of selection on different mitochondrial OXPHOS genes have revealed the adaptive implications of amino acid changes in these subunits. In hares, climatic variation and/or introgression were suggested to be at the origin of such adaptation. Here we looked for evidence of positive selection in three mitochondrial OXPHOS genes, using tests of selection, protein structure modelling and effects of amino acid substitutions on the protein function and stability. We also used statistical models to test for climate and introgression effects on sites under positive selection. Results Our results revealed seven sites under positive selection in ND4 and three sites in Cytb. However, no sites under positive selection were observed in the COX1 gene. All three subunits presented a high number of codons under negative selection. Sites under positive selection were mapped on the tridimensional structure of the predicted models for the respective mitochondrial subunit. Of the ten amino acid replacements inferred to have evolved under positive selection for both subunits, six were located in the transmembrane domain. On the other hand, three codons were identified as sites lining proton translocation channels. Furthermore, four codons were identified as destabilizing with a significant variation of Δ vibrational entropy energy between wild and mutant type. Moreover, our PROVEAN analysis suggested that among all positively selected sites two fixed amino acid replacements altered the protein functioning. Our statistical models indicated significant effects of climate on the presence of ND4 and Cytb protein variants, but no effect by trans-specific mitochondrial DNA introgression, which is not uncommon in a number of hare species. Conclusions Positive selection was observed in several codons in two OXPHOS genes. We found that substitutions in the positively selected codons have structural and functional impacts on the encoded proteins. Our results are concordantly suggesting that adaptations have strongly affected the evolution of mtDNA of hare species with potential effects on the protein function. Environmental/climatic changes appear to be a major trigger of this adaptation, whereas trans-specific introgressive hybridization seems to play no major role for the occurrence of protein variants.

Climate change is considered a major threat to society and nature. UV irradiation is the most important environmental genotoxic agent. Thus, how elevated UV irradiation may influence human health and ecosystems has generated wide concern in the scientific community, as well as with policy makers and the public in general. In this study, we investigated patterns and mechanisms of UV adaptation in natural ecosystems by studying a gene-specific variation in the potato late blight pathogen, Phytophthora infestans. We compared the sequence characteristics of radiation sensitive 23 (RAD23), a gene involved in the nucleotide excision repair (NER) pathway and UV tolerance, in P. infestans isolates sampled from various altitudes. We found that lower genetic variation in the RAD23 gene was caused by natural selection. The hypothesis that UV irradiation drives this selection was supported by strong correlations between the genomic characteristics and altitudinal origin (historic UV irradiation) of the RAD23 sequences with UV tolerance of the P. infestans isolates. These results indicate that the RAD23 gene plays an important role in the adaptation of P. infestans to UV stress. We also found that different climatic factors could work synergistically to determine the evolutionary adaptation of species, making the influence of climate change on ecological functions and resilience more difficult to predict. Future attention should aim at understanding the collective impact generated by simultaneous change in several climate factors on species adaptation and ecological sustainability, using state of the art technologies such as experimental evolution, genome-wide scanning, and proteomics.

Sexual selection is a well-known biological process, yet the genomic basis and patterns of sexual selection are not fully understood. The extravagant ornamental plumage of peacock (Pavo cristatus) was instrumental in shaping Charles Darwin's theory of sexual selection and is considered to be an honest signal of its immunocompetence. Here, we used the recently generated draft genome sequence of peafowl (Pavo cristatus) and carried out a comparative analysis across 11 bird genomes that encompass a range of sexual selection and also had high-quality genomic and phenotypic data publically available to study the genomic basis of sexual selection. We found that varying degree of purifying selection was the predominant mechanism of action for sexual selection at the genome-wide scale and observed that sexual selection mostly influences genes regulating gene expression and protein processing. Specifically, the genome-wide phylogenetically corrected regression analysis supported the continuous or ongoing model of sexual selection. Genes involved in nucleic acid binding and gene expression regulation, including a specific regulator of sex-determination known as TRA2A to be under positive selection in the species with high post-copulatory sexual selection manifested as high sperm competition. We also detected specific feather-related and immune-related gene-pairs evolving under similar selection pressures across the 11 species, including peacock (Pavo cristatus), which is consistent with the Hamilton-Zuk hypothesis. The comparative genomics analysis of 11 avian taxa has provided new insights on the molecular underpinnings of sexual selection and identifies specific genomic regions for future in-depth analysis.

Most cellular functions are carried out by a dynamic network of interacting proteins. An open question is whether the network properties of protein interactomes represent phenotypes under natural selection. One proposal is that protein interactomes have evolved to be resilient, such that they tend to maintain connectivity when proteins are removed from the network. This hypothesis predicts that interactome resilience should be maintained by natural selection during long-term experimental evolution. I tested this prediction by modeling the evolution of protein-protein interaction (PPI) networks in Lenski’s long-term evolution experiment with Escherichia coli (LTEE). In this test, I removed proteins affected by nonsense, insertion, deletion, and transposon mutations in evolved LTEE strains, and measured the resilience of the resulting networks. I compared the rate of change of network resilience in each LTEE population to the rate of change of network resilience for corresponding randomized networks. The evolved PPI networks are significantly more resilient than networks in which random proteins have been deleted.Moreover, the evolved networks are generally more resilient than networks in which the random deletion of proteins was restricted to those disrupted in LTEE. These results suggest that evolution in the LTEE has favored PPI networks that are, on average, more resilient than expected from the genetic variation across the evolved strains. My findings therefore support the hypothesis that selection maintains protein interactome resilience over evolutionary time.Significance StatementUnderstanding how protein-protein interaction (PPI) networks evolve is a central goal of evolutionary systems biology. One property that has been hypothesized to be important for PPI network evolution is resilience, which means that networks tend to maintain connectivity even after many nodes (proteins in this case) have been removed. This hypothesis predicts that PPI network resilience should be maintained during long-term experimental evolution. Consistent with this prediction, I found that the PPI networks that evolved over 50,000 generations of Lenski’s long-term evolution experiment with E. coli are more resilient than expected by chance.

Due to their limited ranges and inherent isolation, island species have long been recognized as crucial systems for tackling a range of evolutionary questions, including in the early study of speciation.1,2 Such species have been less studied in the understanding of the evolutionary forces driving DNA sequence evolution. Island species usually have lower census population sizes (N) than continental species and, supposedly, lower effective population sizes (Ne). Given that both the rates of change caused by genetic drift and by selection are dependent upon Ne, island species are theoretically expected to exhibit (1) lower genetic diversity, (2) less effective natural selection against slightly deleterious mutations,3,4 and (3) a lower rate of adaptive evolution.5, 6, 7, 8 Here, we have used a large set of newly sequenced and published whole-genome sequences of Passerida species (14 insular and 11 continental) to test these predictions. We confirm that island species exhibit lower census size and Ne, supporting the hypothesis that the smaller area available on islands constrains the upper bound of Ne. In the insular species, we find lower nucleotide diversity in coding regions, higher ratios of non-synonymous to synonymous polymorphisms, and lower adaptive substitution rates. Our results provide robust evidence that the lower Ne experienced by island species has affected both the ability of natural selection to efficiently remove weakly deleterious mutations and also the adaptive potential of island species, therefore providing considerable empirical support for the nearly neutral theory. We discuss the implications for both evolutionary and conservation biology., This research was funded by the French ANR (BirdIslandGenomic project, ANR-14-CE02-0002). J.C.I. was funded by the Spanish Ministry of Science, Innovation and Universities (Ref.: PGC2018-097575-B-I00) and by a GRUPIN research grant from the Regional Government of Asturias (Ref.: IDI/2018/000151). B.M. and M.R. were partly funded by the Spanish Ministry of Science and Innovation (grant PGC2018-098897-B-I00), and M.R. was supported by a doctoral fellowship from the Spanish Ministry of Education, Culture, and Sport (FPU16/05724).

The F-box proteins function as substrate receptors to determine the specificity of Skp1-Cul1-F-box ubiquitin ligases. Genomic studies revealed large and diverse sizes of the F-box gene superfamily across plant species. Our previous studies suggested that the plant F-box gene superfamily is under genomic drift evolution promoted by epigenomic programming. However, how the size of the superfamily drifts across plant genomes is currently unknown. Through a large-scale genomic and phylogenetic comparison of the F-box gene superfamily covering 110 green plants and one red algal species, I discovered four distinct groups of plant F-box genes with diverse evolutionary processes. While the members in Clusters 1 and 2 are species/lineage-specific, those in Clusters 3 and 4 are present in over 46 plant genomes. Statistical modeling suggests that F-box genes from the former two groups are skewed toward fewer species and more paralogs compared to those of the latter two groups whose presence frequency and sizes in plant genomes follow a random statistical model. The enrichment of known Arabidopsis F-box genes in Clusters 3 and 4, along with comprehensive biochemical evidence showing that Arabidopsis members in Cluster 4 interact with the Arabidopsis Skp1-like 1 (ASK1), demonstrates over-representation of active F-box genes in these two groups. Collectively, I propose purifying and dosage balancing selection models to explain the lineage/species-specific duplications and expansions of F-box genes in plant genomes. The purifying selection model suggests that most, if not all, lineage/species-specific F-box genes are detrimental and are thus kept at low frequencies in plant genomes.

© 2021 by the authors., The current availability of complete genome sequences has allowed knowing that bacterial genomes can bear genes not present in the genome of all the strains from a specific species. So, the genes shared by all the strains comprise the core of the species, but the pangenome can be much greater and usually includes genes appearing in one only strain. Once the pangenome of a species is estimated, other studies can be undertaken to generate new knowledge, such as the study of the evolutionary selection for protein-coding genes. Most of the genes of a pangenome are expected to be subject to purifying selection that assures the conservation of function, especially those in the core group. However, some genes can be subject to selection pressure, such as genes involved in virulence that need to escape to the host immune system, which is more common in the accessory group of the pangenome. We analyzed 180 strains of Helicobacter pylori, a bacterium that colonizes the gastric mucosa of half the world population and presents a low number of genes (around 1500 in a strain and 3000 in the pangenome). After the estimation of the pangenome, the evolutionary selection for each gene has been calculated, and we found that 85% of them are subject to purifying selection and the remaining genes present some grade of selection pressure. As expected, the latter group is enriched with genes encoding for membrane proteins putatively involved in interaction to host tissues. In addition, this group also presents a high number of uncharacterized genes and genes encoding for putative spurious proteins. It suggests that they could be false positives from the gene finders used for identifying them. All these results propose that this kind of analyses can be useful to validate gene predictions and functionally characterize proteins in complete genomes.

Although it is well known that abundant proteins evolve slowly across the tree of life, there is little consensus for why this is true. Here, I report that abundant proteins evolve slowly in the hypermutator populations of Lenski’s long-term evolution experiment with Escherichia coli (LTEE). Specifically, the density of all observed mutations per gene, as measured in metagenomic time series covering 60,000 generations of the LTEE, significantly anti-correlates with mRNA abundance, protein abundance, and degree of protein-protein interaction. The same pattern holds for nonsynonymous mutation density. However, synonymous mutation density, measured across the LTEE hypermutator populations, positively correlates with protein abundance. These results show that universal constraints on protein evolution are visible in data spanning three decades of experimental evolution. Therefore, it should be possible to design experiments to answer why abundant proteins evolve slowly.Significance StatementA universal evolutionary pattern is that highly abundant and highly interacting proteins evolve slowly. This pattern was discovered in analyses that cover millions of years’ worth of sequence variation, so it is not clear how long it takes (decades, centuries, millennia) for such patterns to emerge. Here, I report that this universal evolutionary pattern emerges in metagenomic data that cover just 30 years of experimental evolution.

Atypical porcine pestivirus (APPV) is a single-stranded RNA virus from the family Flaviviridae, which is linked to congenital tremor (CT) type A-II in newborn piglets. Here, we retrospectively investigated the molecular evolution of APPV on an affected herd between 2013 and 2019. Monitoring was done at regular intervals, and the same genotype of APPV was found during the entire study period, suggesting no introductions from outside the farm. The nucleotide substitutions over time did not show substantial amino acid variation in the structural glycoproteins. Furthermore, the evolution of the virus showed mainly purifying selection, and no positive selection. The limited pressure on the virus to change at immune-dominant regions suggested that the immune pressure at the farm might be low. In conclusion, farms can have circulation of APPV for years, and massive testing and removal of infected animals are not sufficient to clear the virus from affected farms.

Pterocarpus is a genus of trees mainly distributed in tropical Asia, Africa, and South America. Some species of Pterocarpus are rosewood tree species, having important economic value for timber, and for some species, medicinal value as well. Up to now, information about this genus with regard to the genomic characteristics of the chloroplasts has been limited. Based on a combination of next-generation sequencing (Illumina Hiseq) and long-read sequencing (PacBio), the whole chloroplast genomes (cp genomes) of five species (rosewoods) in Pterocarpus (Pterocarpus macrocarpus, P. santalinus, P. indicus, P. pedatus, P. marsupium) have been assembled. The cp genomes of five species in Pterocarpus have similar structural characteristics, gene content, and sequence to other flowering plants. The cp genomes have a typical four-part structure, containing 110 unique genes (77 protein coding genes, 4 rRNAs, 29 tRNAs). Through comparative genomic analysis, abundant simple sequence repeat (SSR)loci (333&ndash, 349) were detected in Pterocarpus, among which A /T single nucleotide repeats accounted for the highest proportion (72.8&ndash, 76.4%). In the five cp genomes of Pterocarpus, eight hypervariable regions, including trnH-GUG_psbA, trnS-UGA_psbC, accD-psaI, ndhI-exon2_ndhI-exon1, ndhG_ndhi-exon2, rpoC2-exon2, ccsA, and trnfM-CAU, are proposed for use as DNA barcode regions. In the comparison of gene selection pressures (P. santalinus as the reference genome), purifying selection was inferred as the primary mode of selection in maintaining important biological functions. Phylogenetic analysis shows that Pterocarpus is a monophyletic group. The species P. tinctorius is resolved as early diverging in the genus. Pterocarpus was resolved as sister to the genus Tipuana.

The emergence of SARS-CoV-2 has resulted in more than 200,000 infections and nearly 9,000 deaths globally so far. This novel virus is thought to have originated from an animal reservoir, and acquired the ability to infect human cells using the SARS-CoV cell receptor hACE2. In the wake of a global pandemic it is essential to improve our understanding of the evolutionary dynamics surrounding the origin and spread of a novel infectious disease. One way theory predicts selection pressures should shape viral evolution is to enhance binding with host cells. We first assessed evolutionary dynamics in select betacoronavirus spike protein genes to predict where these genomic regions are under directional or purifying selection between divergent viral lineages at various scales of relatedness. With this analysis, we determine a region inside the receptor-binding domain with putative sites under positive selection interspersed among highly conserved sites, which are implicated in structural stability of the viral spike protein and its union with human receptor hACE2. Next, to gain further insights into factors associated with coronaviruses recognition of the human host receptor, we performed modeling studies of five different coronaviruses and their potential binding to hACE2. Modeling results indicate that interfering with the salt bridges at hot spot 353 could be an effective strategy for inhibiting binding, and hence for the prevention of coronavirus infections. We also propose that a glycine residue at the receptor binding domain of the spike glycoprotein can have a critical role in permitting bat variants of the coronaviruses to infect human cells.

Transposable elements (TEs) contribute to a large fraction of the expansion of many eukaryotic genomes due to the capability of TEs duplicating themselves through transposition. A first step to understanding the roles of TEs in a eukaryotic genome is to characterize the population-wide variation of TE insertions in the species. Here, we present a maximum-likelihood (ML) method for estimating allele frequencies and detecting selection on TE insertions in a diploid population, based on the genotypes at TE insertion sites detected in multiple individuals sampled from the population using paired-end (PE) sequencing reads. Tests of the method on simulated data show that it can accurately estimate the allele frequencies of TE insertions even when the PE sequencing is conducted at a relatively low coverage (=5X). The method can also detect TE insertions under strong selection, and the detection ability increases with sample size in a population, although a substantial fraction of actual TE insertions under selection may be undetected. Application of the ML method to genomic sequencing data collected from a natural Daphnia pulex population shows that, on the one hand, most (>90%) TE insertions present in the reference D. pulex genome are either fixed or nearly fixed (with allele frequencies >0.95); on the other hand, among the nonreference TE insertions (i.e., those detected in some individuals in the population but absent from the reference genome), the majority (>70%) are still at low frequencies (

Highly conserved noncoding elements (CNEs) comprise a significant proportion of the genomes of multicellular eukaryotes. The function of most CNEs remains elusive, but growing evidence indicates they are under some form of purifying selection. Noncoding regions in many species also harbor large numbers of transposable element (TE) insertions, which are typically lineage specific and depleted in exons because of their deleterious effects on gene function or expression. However, it is currently unknown whether the landscape of TE insertions in noncoding regions is random or influenced by purifying selection on CNEs. Here we combine comparative and population genomic data inDrosophila melanogasterto show that abundance of TE insertions in intronic and intergenic CNEs is reduced relative to random expectation, supporting the idea that selective constraints on CNEs eliminate a proportion of TE insertions in noncoding regions. However, we find no difference in the allele frequency spectra for polymorphic TE insertions in CNEs versus those in unconstrained spacer regions, suggesting that the distribution of fitness effects acting on observable TE insertions is similar across different functional compartments in noncoding DNA. Our results provide evidence that selective constraints on CNEs contribute to shaping the landscape of TE insertion in eukaryotic genomes, and provide further evidence supporting the conclusion that CNEs are indeed functionally constrained and not simply mutational cold spots.

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.

Background Antigen-detecting rapid diagnostic tests (RDTs) have been recommended by the World Health Organization for use in remote areas to improve malaria case management. Lactate dehydrogenase (LDH) of Plasmodium falciparum is one of the main parasite antigens employed by various commercial RDTs. It has been hypothesized that the poor detection of LDH-based RDTs is attributed in part to the sequence diversity of the gene. To test this, the present study aimed to investigate the genetic diversity of the P. falciparum ldh gene in Thailand and to construct the map of LDH sequence diversity in P. falciparum populations worldwide. Methods The ldh gene was sequenced for 50 P. falciparum isolates in Thailand and compared with hundreds of sequences from P. falciparum populations worldwide. Several indices of molecular variation were calculated, including the proportion of polymorphic sites, the average nucleotide diversity index (π), and the haplotype diversity index (H). Tests of positive selection and neutrality tests were performed to determine signatures of natural selection on the gene. Mean genetic distance within and between species of Plasmodium ldh was analysed to infer evolutionary relationships. Results Nucleotide sequences of P. falciparum ldh could be classified into 9 alleles, encoding 5 isoforms of LDH. L1a was the most common allelic type and was distributed in P. falciparum populations worldwide. Plasmodium falciparum ldh sequences were highly conserved, with haplotype and nucleotide diversity values of 0.203 and 0.0004, respectively. The extremely low genetic diversity was maintained by purifying selection, likely due to functional constraints. Phylogenetic analysis inferred the close genetic relationship of P. falciparum to malaria parasites of great apes, rather than to other human malaria parasites. Conclusions This study revealed the global genetic variation of the ldh gene in P. falciparum, providing knowledge for improving detection of LDH-based RDTs and supporting the candidacy of LDH as a therapeutic drug target. Electronic supplementary material The online version of this article (10.1186/s12936-017-2157-5) contains supplementary material, which is available to authorized users.

Phosphoglucose isomerase (PGI, EC 5.3.1.9) is an essential metabolic enzyme in all eukaryotes. An earlier study of the PgiC1 gene, which encodes cytosolic PGI in the grass Festuca ovina L., revealed a marked difference in the levels of nucleotide polymorphism between the 5’ and 3’ portions of the gene. In the present study, we characterized the sequence polymorphism in F. ovina PgiC1 in more detail and examined possible explanations for the non-uniform pattern of nucleotide polymorphism across the gene. Our study confirms that the two portions of the PgiC1 gene show substantially different levels of DNA polymorphism and also suggests that the peptide encoded by the 3’ portion of PgiC1 is functionally and structurally more important than that encoded by the 5’ portion. Although there was some evidence of purifying selection (d N/d S test) on the 5’ portion of the gene, the signature of purifying selection was considerably stronger on the 3’ portion of the gene (d N/d S and McDonald–Kreitman tests). Several tests support the action of balancing selection within the 5’ portion of the gene. Wall’s B and Q tests were significant only for the 5’ portion of the gene. There were also marked peaks of nucleotide diversity, Tajima’s D and the d N/d S ratio at or around a PgiC1 codon site (within the 5’ portion of the gene) that a previous study had suggested was subject to positive diversifying selection. Our results suggest that the two portions of the gene have been subject to different selective regimes. Purifying selection appears to have been the main force contributing to the relatively low level of polymorphism within the 3’ portion of the sequence. In contrast, it is possible that balancing selection has contributed to the maintenance of the polymorphism within the 5’ portion of the gene.

Understanding the molecular evolution of macromolecular complexes in the light of their structure, assembly, and stability is of central importance. Here, we address how the modular organization of native molecular contacts shapes the selection pressure on individual residue sites of ribosomal complexes. The bacterial ribosomal complex is represented as a residue contact network where nodes represent amino acid/nucleotide residues and edges represent their van der Waals interactions. We find statistically overrepresented native amino acid–nucleotide contacts (OaantC, one amino acid contacts one or multiple nucleotides, internucleotide contacts are disregarded). Contact number is defined as the number of nucleotides contacted. Involvement of individual amino acids in OaantCs with smaller contact numbers is more random, whereas only a few amino acids significantly contribute to OaantCs with higher contact numbers. An investigation of structure, stability, and assembly of bacterial ribosome depicts the involvement of these OaantCs in diverse biophysical interactions stabilizing the complex, including high-affinity protein–RNA contacts, interprotein cooperativity, intersubunit bridge, packing of multiple ribosomal RNA domains, etc. Amino acid–nucleotide constituents of OaantCs with higher contact numbers are generally associated with significantly slower substitution rates compared with that of OaantCs with smaller contact numbers. This evolutionary rate heterogeneity emerges from the strong purifying selection pressure that conserves the respective amino acid physicochemical properties relevant to the stabilizing interaction with OaantC nucleotides. An analysis of relative molecular orientations of OaantC residues and their interaction energetics provides the biophysical ground of purifying selection conserving OaantC amino acid physicochemical properties.

Neovascular age-related macular degeneration (nvAMD) is characterized by choroidal blood vessels growing into the subretinal space, leading to retinal pigment epithelial (RPE) cell degeneration and vision loss. Vessel growth results from an imbalance of pro-angiogenic (e.g., vascular endothelial growth factor [VEGF]) and anti-angiogenic factors (e.g., pigment epithelium-derived factor [PEDF]). Current treatment using intravitreal injections of anti-VEGF antibodies improves vision in about 30% of patients but may be accompanied by side effects and non-compliance. To avoid the difficulties posed by frequent intravitreal injections, we have proposed the transplantation of pigment epithelial cells modified to overexpress human PEDF. Stable transgene integration and expression is ensured by the hyperactive Sleeping Beauty transposon system delivered by pFAR4 miniplasmids, which have a backbone free of antibiotic resistance markers. We demonstrated efficient expression of the PEDF gene and an optimized PEDF cDNA sequence in as few as 5 × 103 primary cells. At 3 weeks post-transfection, PEDF secretion was significantly elevated and long-term follow-up indicated a more stable secretion by cells transfected with the optimized PEDF transgene. Analysis of transgene insertion sites in human RPE cells showed an almost random genomic distribution. The results represent an important contribution toward a clinical trial aiming at a non-viral gene therapy of nvAMD. Keywords: neovascular age-related macular degeneration, non-viral ex vivo gene therapy, sleeping beauty transposon system, pFAR4 miniplasmids, primary pigment epithelial cells, pigment epithelium-derived factor

Seoul virus (SEOV), which causes hemorrhagic fever with renal syndrome (HFRS) in humans, has spread all over the world, especially in mainland China. Understanding basic mechanisms of SEOV evolution is essential to better combat and prevent viral diseases. Here, we examined SEOV prevalence and evolution in the residential area of four districts in Guangzhou city, China. The carriage of SEOV was observed in 33.33% of the sampled rodents, with 35.96% of the sampled Rattus norvegicus and 13.33% of R. tanezumi. Based on the comprehensive analyses of large (L), medium (M), and small (S) segments, our study first demonstrated that the genetic characterization of urban SEOV was shaped by high nucleotide substitution rates, purifying selection, and recombination. Additionally, we detected mutational saturation in the S segment of SEOV, which may lead to the biases of genetic divergence and substitution rates in our study. Importantly, we have filled the gap of SEOV evolution in the urban area. The genetic variation of SEOV may highlight the risk of HFRS, which merits further investigation.

BackgroundSalmonid fishes are characterised by a very high level of variation in trophic, ecological, physiological, and life history adaptations. Some salmonid taxa show exceptional potential for fast, within-lake diversification into morphologically and ecologically distinct variants, often in parallel; these are the lake-resident charr and whitefish (several species in the generaSalvelinusandCoregonus). To identify selection on genes and gene categories associated with such predictable diversifications, we analysed 2702 orthogroups (4.82 Mbp total; average 4.77 genes/orthogroup; average 1783 bp/orthogroup). We did so in two charr and two whitefish species and compared to five other salmonid lineages, which do not evolve in such ecologically predictable ways, and one non-salmonid outgroup.ResultsAll selection analyses are based onCoregonusandSalvelinuscompared to non-diversifying taxa. We found more orthogroups were affected by relaxed selection than intensified selection. Of those, 122 were under significant relaxed selection, with trends of an overrepresentation of serine family amino acid metabolism and transcriptional regulation, and significant enrichment of behaviour-associated gene functions. Seventy-eight orthogroups were under significant intensified selection and were enriched for signalling process and transcriptional regulation gene ontology terms and actin filament and lipid metabolism gene sets. Ninety-two orthogroups were under diversifying/positive selection. These were enriched for signal transduction, transmembrane transport, and pyruvate metabolism gene ontology terms and often contained genes involved in transcriptional regulation and development. Several orthogroups showed signs of multiple types of selection. For example, orthogroups under relaxed and diversifying selection contained genes such asap1m2, involved in immunity and development, andslc6a8, playing an important role in muscle and brain creatine uptake. Orthogroups under intensified and diversifying selection were also found, such as genessyn3, with a role in neural processes, andctsk, involved in bone remodelling.ConclusionsOur approach pinpointed relevant genomic targets by distinguishing among different kinds of selection. We found that relaxed, intensified, and diversifying selection affect orthogroups and gene functions of ecological relevance in salmonids. Because they were found consistently and robustly across charr and whitefish and not other salmonid lineages, we propose these genes have a potential role in the replicated ecological diversifications.

Intrinsically disordered regions make up a large part of the proteome, but the sequence-to-function relationship in these regions is poorly understood, in part because the primary amino acid sequences of these regions are poorly conserved in alignments. Here we use an evolutionary approach to detect molecular features that are preserved in the amino acid sequences of orthologous intrinsically disordered regions. We find that most disordered regions contain multiple molecular features that are preserved, and we define these as “evolutionary signatures” of disordered regions. We demonstrate that intrinsically disordered regions with similar evolutionary signatures can rescue functionin vivo,and that groups of intrinsically disordered regions with similar evolutionary signatures are strongly enriched for functional annotations and phenotypes. We propose that evolutionary signatures can be used to predict function for many disordered regions from their amino acid sequences.

Background The evolution of male pregnancy is the most distinctive characteristic of syngnathids, and their specialized life history traits make syngnathid species excellent model species for many issues in biological evolution. However, the origin of syngnathids and the evolutionary divergence time of different syngnathid species remain poorly resolved. Comprehensive phylogenetic studies of the Syngnathidae will provide critical evidence to elucidate their origin, evolution, and dispersal patterns. Results We sequenced the mitochondrial genomes of eight syngnathid species in this study, and the estimated divergence times suggested that syngnathids diverged from other teleosts approximately 48.8 Mya during the Eocene period. Selection analysis showed that many mitochondrial genes of syngnathids exhibited significantly lower Ka/Ks values than those of other teleosts. The two most frequently used codons in syngnathid fishes were different from those in other teleosts, and a greater proportion of the mitochondrial simple sequence repeats (SSRs) were distributed in non-coding sequences in syngnathids compared with other teleosts. Conclusions Our study indicated that syngnathid fishes experienced an adaptive radiation process during the early explosion of species. Syngnathid mitochondrial OXPHOS genes appear to exhibit depressed Ka/Ks ratios compared with those of other teleosts, and this may suggest that their mitogenomes have experienced strong selective constraints to eliminate deleterious mutations. Electronic supplementary material The online version of this article (10.1186/s12862-019-1430-3) contains supplementary material, which is available to authorized users.

Very little is known about aquatic RNA virus populations and genome evolution. This is the first study that analyzes marine environmental RNA viral assemblages in an evolutionary and broad geographical context. This study contributes the largest marine RNA virus metagenomic data set to date, substantially increasing the sequencing space for RNA viruses and also providing a baseline for comparisons of marine RNA virus diversity. The new viruses discovered in this study are representative of the most abundant family of marine RNA viruses, the Marnaviridae, and expand our view of the diversity of this important group. Overall, our data and analyses provide a foundation for interpreting marine RNA virus diversity and evolution., RNA viruses, particularly genetically diverse members of the Picornavirales, are widespread and abundant in the ocean. Gene surveys suggest that there are spatial and temporal patterns in the composition of RNA virus assemblages, but data on their diversity and genetic variability in different oceanographic settings are limited. Here, we show that specific RNA virus genomes have widespread geographic distributions and that the dominant genotypes are under purifying selection. Genomes from three previously unknown picorna-like viruses (BC-1, -2, and -3) assembled from a coastal site in British Columbia, Canada, as well as marine RNA viruses JP-A, JP-B, and Heterosigma akashiwo RNA virus exhibited different biogeographical patterns. Thus, biotic factors such as host specificity and viral life cycle, and not just abiotic processes such as dispersal, affect marine RNA virus distribution. Sequence differences relative to reference genomes imply that virus quasispecies are under purifying selection, with synonymous single-nucleotide variations dominating in genomes from geographically distinct regions resulting in conservation of amino acid sequences. Conversely, sequences from coastal South Africa that mapped to marine RNA virus JP-A exhibited more nonsynonymous mutations, probably representing amino acid changes that accumulated over a longer separation. This biogeographical analysis of marine RNA viruses demonstrates that purifying selection is occurring across oceanographic provinces. These data add to the spectrum of known marine RNA virus genomes, show the importance of dispersal and purifying selection for these viruses, and indicate that closely related RNA viruses are pathogens of eukaryotic microbes across oceans. IMPORTANCE Very little is known about aquatic RNA virus populations and genome evolution. This is the first study that analyzes marine environmental RNA viral assemblages in an evolutionary and broad geographical context. This study contributes the largest marine RNA virus metagenomic data set to date, substantially increasing the sequencing space for RNA viruses and also providing a baseline for comparisons of marine RNA virus diversity. The new viruses discovered in this study are representative of the most abundant family of marine RNA viruses, the Marnaviridae, and expand our view of the diversity of this important group. Overall, our data and analyses provide a foundation for interpreting marine RNA virus diversity and evolution.

The mitochondrial genome is an evolutionarily persistent and cooperative component of metazoan cells that contributes to energy production and many other cellular processes. Despite sharing the same host as the nuclear genome, the multi-copy mitochondrial DNA (mtDNA) follows very different rules of replication and transmission, which translate into differences in the patterns of selection. On one hand, mtDNA is dependent on the host for its transmission, so selections would favour genomes that boost organismal fitness. On the other hand, genetic heterogeneity within an individual allows different mitochondrial genomes to compete for transmission. This intra-organismal competition could select for the best replicator, which does not necessarily give the fittest organisms, resulting in mito-nuclear conflict. In this review, we discuss the recent advances in our understanding of the mechanisms and opposing forces governing mtDNA transmission and selection in bilaterians, and what the implications of these are for mtDNA evolution and mitochondrial replacement therapy.

Background Saussurea DC. is one of the largest and most morphologically heterogeneous genera in Asteraceae. The relationships within Saussurea have been poorly resolved, probably due an early, rapid radiation. To examine plastome evolution and resolve backbone relationships within Saussurea, we sequenced the complete plastomes of 17 species representing all four subgenera. Results All Saussurea plastomes shared the gene content and structure of most Asteraceae plastomes. Molecular evolutionary analysis showed most of the plastid protein-coding genes have been under purifying selection. Phylogenomic analyses of 20 Saussurea plastomes that alternatively included nucleotide or amino acid sequences of all protein-coding genes, vs. the nucleotide sequence of the entire plastome, supported the monophyly of Saussurea and identified three clades within it. Three of the four traditional subgenera were recovered as paraphyletic. Seven plastome regions were identified as containing the highest nucleotide variability. Conclusions Our analyses reveal both the structural conservatism and power of the plastome for resolving relationships in congeneric taxa. It is very likely that differences in topology among data sets is due primarily to differences in numbers of parsimony-informative characters. Our study demonstrates that the current taxonomy of Saussurea is likely based at least partly on convergent morphological character states. Greater taxon sampling will be necessary to explore character evolution and biogeography in the genus. Our results here provide helpful insight into which loci will provide the most phylogenetic signal in Saussurea and Cardueae. Electronic supplementary material The online version of this article (10.1186/s12870-019-1896-6) contains supplementary material, which is available to authorized users.

Segmental duplications are an important class of mutations. Because a large proportion of segmental duplications may often be strongly deleterious, high frequency or fixed segmental duplications may represent only a tiny fraction of the mutational input. To understand the emergence and elimination of segmental duplications, we survey polymorphic duplications, including tandem and interspersed duplications, in natural populations of Drosophila by haploid embryo genomes. As haploid embryos are not expected to be heterozygous, the genome, sites of heterozygosity (referred to as pseudoheterozygous sites [PHS]), may likely represent recent duplications that have acquired new mutations. Among the 29 genomes of Drosophila melanogaster, we identify 2,282 polymorphic PHS duplications (linked PHS regions) in total or 154 PHS duplications per genome. Most PHS duplications are small (83.4% < 500 bp), Drosophila melanogaster lineage specific, and strain specific (72.6% singletons). The excess of the observed singleton PHS duplications deviates significantly from the neutral expectation, suggesting that most PHS duplications are strongly deleterious. In addition, these small segmental duplications are not evenly distributed in genomic regions and less common in noncoding functional element regions. The underrepresentation in RNA polymerase II binding sites and regions with active histone modifications is correlated with ages of duplications. In conclusion, small segmental duplications occur frequently in Drosophila but rapidly eliminated by natural selection.

The transcriptional regulatory structure of plant genomes is still relatively unexplored, and little is known about factors that influence expression variation in plants. We used a genetic system consisting of 10 heterozygous grape varieties with high consanguinity and high haplotypic diversity to: (i) identify regions of haplotype sharing through whole-genome resequencing and single-nucleotide polymorphism (SNP) genotyping; (ii) analyse gene expression through RNA-seq in four stages of berry development; and (iii) associate gene expression variation with genetic and epigenetic properties. We found that haplotype sharing in and around genes was positively correlated with similarity in expression and was negatively correlated with the fraction of differentially expressed genes. Genetic and epigenetic properties of the gene and the surrounding region showed significant effects on the extent of expression variation, with negative associations for the level of gene body methylation and mean expression level, and with positive associations for nucleotide diversity, structural diversity and ratio of non-synonymous to synonymous nucleotide diversity. We also observed a spatial dependency of covariation of gene expression among varieties. These results highlight relevant roles for cis-acting factors, selective constraints and epigenetic features of the gene, and the regional context in which the gene is located, in the determination of expression variation. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/bioproject/PRJNA385116; https://www.ncbi.nlm.nih.gov/bioproject/PRJNA392287; https://www.ncbi.nlm.nih.gov/bioproject/PRJNA373967 (released upon publication); https://www.ncbi.nlm.nih.gov/bioproject/PRJNA490160 (released upon publication); https://www.ncbi.nlm.nih.gov/bioproject/PRJNA265039; https://www.ncbi.nlm.nih.gov/bioproject/PRJNA265040.

[Background] Determining the factors involved in the likelihood of a gene being under adaptive selection is still a challenging goal in Evolutionary Biology. Here, we perform an evolutionary analysis of the human metabolic genes to explore the associations between network structure and the presence and strength of natural selection in the genes whose products are involved in metabolism. Purifying and positive selection are estimated at interspecific (among mammals) and intraspecific (among human populations) levels, and the connections between enzymatic reactions are differentiated between incoming (in-degree) and outgoing (out-degree) links., [Results] We confirm that purifying selection has been stronger in highly connected genes. Long-term positive selection has targeted poorly connected enzymes, whereas short-term positive selection has targeted different enzymes depending on whether the selective sweep has reached fixation in the population: genes under a complete selective sweep are poorly connected, whereas those under an incomplete selective sweep have high out-degree connectivity. The last steps of pathways are more conserved due to stronger purifying selection, with long-term positive selection targeting preferentially enzymes that catalyze the first steps. However, short-term positive selection has targeted enzymes that catalyze the last steps in the metabolic network. Strong signals of positive selection have been found for metabolic processes involved in lipid transport and membrane fluidity and permeability., [Conclusions] Our analysis highlights the importance of analyzing the same biological system at different evolutionary timescales to understand the evolution of metabolic genes and of distinguishing between incoming and outgoing links in a metabolic network. Short-term positive selection has targeted enzymes with a different connectivity profile depending on the completeness of the selective sweep, while long-term positive selection has targeted genes with fewer connections that code for enzymes that catalyze the first steps in the network., This study has been possible thanks to grant BFU2016–77961-P (AEI/FEDER, UE) awarded by the Agencia Estatal de Investigación (Ministerio de Ciencia, Innovación y Universidades, Spain) and with the support of Secretaria d’Universitats i Recerca del Departament d’Economia i Coneixement de la Generalitat de Catalunya (GRC 2017 527 SGR 702) to JB. Part of the “Unidad de Excelencia María de Maeztu” (MDM-2014-0370), funded by the Ministerio de Economía, Industria y Competividad (MINECO, Spain). JP work is supported in part by the Agencia Estatal de Investigación (Ministerio de Ciencia, Innovación y Universidades, Spain) Helios grant reference: BIO2015–66960-C3–1-R (co-financed by FEDER) and by the European Union through the BioRoboost project (H2020-NMBP-TR-IND-2018-2020/BIOTEC-01-2018 (CSA), Project ID 210491758). BD is supported by F.P.U. grant FPU13/06813 from the Ministerio de Educación, Cultura y Deporte (Spain).

BackgroundTransposable elements (TEs) can be key drivers of evolution, but the mechanisms and scope of how they impact gene and genome function are largely unknown. Previous analyses revealed that TE-mediated gene amplifications can have variable effects on fungal genomes, from inactivation of function to production of multiple active copies. For example, a DNA methyltransferase gene in the wheat pathogenZymoseptoria tritici(synonymMycosphaerella graminicola) was amplified to tens of copies, all of which were inactivated by Repeat-Induced Point mutation (RIP) including the original, resulting in loss of cytosine methylation. In another wheat pathogen,Pyrenophora tritici-repentis, a histone H3 gene was amplified to tens of copies with little evidence of RIP, leading to many potentially active copies. To further test the effects of transposon-aided gene amplifications on genome evolution and architecture, the repetitive fraction of the significantly expandedPseudocercospora fijiensisgenome was analyzed in greater detail.ResultsThese analyses identified a housekeeping gene, histone H3, which was captured and amplified to hundreds of copies by ahATDNA transposon, all of which were inactivated by RIP, except for the original. InP. fijiensisthe original H3 gene probably was not protected from RIP, but most likely was maintained intact due to strong purifying selection. Comparative analyses revealed that a similar event occurred in five additional genomes representing the fungal generaCercospora, PseudocercosporaandSphaerulina.ConclusionsThese results indicate that the interplay of TEs and RIP can result in different and unpredictable fates of amplified genes, with variable effects on gene and genome evolution.

International audience; Genetic differentiation occurs when gene flow is prevented, due to reproductive barriers or asexuality. Investigating the early barriers to gene flow is important for understanding the process of speciation. Here, we therefore investigated reproductive isolation between different genetic clusters of the fungus Penicillium roqueforti, used for maturing blue cheeses, and also occurring as food spoiler or in silage. We investigated premating and postmating fertility between and within three genetic clusters (two from cheese and one from other substrates), and we observed sexual structures under scanning electron microscopy. All intercluster types of crosses showed some fertility, suggesting that no intersterility has evolved between domesticated and wild populations despite adaptation to different environments and lack of gene flow. However, much lower fertility was found in crosses within the cheese clusters than within the noncheese cluster, suggesting reduced fertility of cheese strains, which may constitute a barrier to gene flow. Such degeneration may be due to bottlenecks during domestication and/or to the exclusive clonal replication of the strains in industry. This study shows that degeneration has occurred rapidly and independently in two lineages of a domesticated species. Altogether, these results inform on the processes and tempo of degeneration and speciation.

Although previous studies have documented a bottleneck in the transmission of mtDNA genomes from mothers to offspring, several aspects remain unclear, including the size and nature of the bottleneck. Here, we analyze the dynamics of mtDNA heteroplasmy transmission in the Genomes of the Netherlands (GoNL) data, which consists of complete mtDNA genome sequences from 228 trios, eight dizygotic (DZ) twin quartets, and 10 monozygotic (MZ) twin quartets. Using a minor allele frequency (MAF) threshold of 2%, we identified 189 heteroplasmies in the trio mothers, of which 59% were transmitted to offspring, and 159 heteroplasmies in the trio offspring, of which 70% were inherited from the mothers. MZ twin pairs exhibited greater similarity in MAF at heteroplasmic sites than DZ twin pairs, suggesting that the heteroplasmy MAF in the oocyte is the major determinant of the heteroplasmy MAF in the offspring. We used a likelihood method to estimate the effective number of mtDNA genomes transmitted to offspring under different bottleneck models; a variable bottleneck size model provided the best fit to the data, with an estimated mean of nine individual mtDNA genomes transmitted. We also found evidence for negative selection during transmission against novel heteroplasmies (in which the minor allele has never been observed in polymorphism data). These novel heteroplasmies are enhanced for tRNA and rRNA genes, and mutations associated with mtDNA diseases frequently occur in these genes. Our results thus suggest that the female germ line is able to recognize and select against deleterious heteroplasmies.

Understanding the plastid genome is extremely important for the interpretation of the genetic mechanisms associated with essential physiological and metabolic functions, the identification of possible marker regions for phylogenetic or phylogeographic analyses, and the elucidation of the modes through which natural selection operates in different regions of this genome. In the present study, we assembled the plastid genome of Artocarpus camansi, compared its repetitive structures with Artocarpus heterophyllus, and searched for evidence of synteny within the family Moraceae. We also constructed a phylogeny based on 56 chloroplast genes to assess the relationships among three families of the order Rosales, that is, the Moraceae, Rhamnaceae, and Cannabaceae. The plastid genome of A. camansi has 160,096 bp, and presents the typical circular quadripartite structure of the Angiosperms, comprising a large single copy (LSC) of 88,745 bp and a small single copy (SSC) of 19,883 bp, separated by a pair of inverted repeat (IR) regions each with a length of 25,734 bp. The total GC content was 36.0%, which is very similar to Artocarpus heterophyllus (36.1%) and other moraceous species. A total of 23,068 codons and 80 SSRs were identified in the A. camansi plastid genome, with the majority of the SSRs being mononucleotide (70.0%). A total of 50 repeat structures were observed in the A. camansi plastid genome, in contrast with 61 repeats in A. heterophyllus. A purifying selection signal was found in 70 of the 79 protein-coding genes, indicating that they have all been highly conserved throughout the evolutionary history of the genus. The comparative analysis of the structural characteristics of the chloroplast among different moraceous species found a high degree of similarity in the sequences, which indicates a highly conserved evolutionary model in these plastid genomes. The phylogenetic analysis also recovered a high degree of similarity between the chloroplast genes of A. camansi and A. heterophyllus, and reconfirmed the hypothesis of the intense conservation of the plastome in the family Moraceae.

Background: Mitochondrial DNA (mtDNA) has long been used to date historical demographic events. The idea that it is useful for molecular dating rests on the premise that its evolution is neutral. Even though this idea has long been challenged, the evidence against clock-like evolution of mtDNA is often ignored. Here, we present a particularly clear and simple example to illustrate the implications of violations of the assumption of selective neutrality. Methods: DNA sequences were generated for the mtDNA COI gene and the nuclear 28S rRNA of two closely related rocky shore snails, and species-level variation was compared. Nuclear rRNA is not usually used to study intraspecific variation in species that are not spatially structured, presumably because this marker is assumed to evolve so slowly that it is more suitable for phylogenetics. Results: Even though high inter-specific divergence reflected the faster evolutionary rate of COI, intraspecific genetic variation was similar for both markers. As a result, estimates of population expansion times based on mismatch distributions differed between the two markers by millions of years. Conclusions: Assuming that 28S evolution is more clock-like, these findings can be explained by variation-reducing purifying selection in mtDNA at the species level, and an elevated divergence rate caused by diversifying selection between the two species. Although these two selective forces together make mtDNA suitable as a marker for species identifications by means of DNA barcoding because they create a ‘barcoding gap’, estimates of demographic change based on this marker can be expected to be highly unreliable. Our study contributes to the growing evidence that the utility of mtDNA sequence data beyond DNA barcoding is limited.

Made available in DSpace on 2020-12-12T01:22:14Z (GMT). No. of bitstreams: 0 Previous issue date: 2020-04-01 Terpenoids are a diverse class of metabolites that impact plant metabolism in response to environmental cues. They are synthesized either via a predominantly cytosolic (MVA) pathway or a plastidic pathway (MEP). In Arabidopsis, several enzymes from the MVA and MEP pathways are encoded by gene families, excluding MVK and DXR, which are single-copy genes. In this study, we assess the diversity, evolution and expression of DXR and MVK genes in selected angiosperms and Coffea arabica in particular. Evolutionary analysis revealed that DXR and MVK underwent purifying selection, but the selection effect for DXR was stronger than it was for MVK. Digital gene expression (DGE) profile analysis of six species revealed that expression levels of MVK in flowers and roots were high, whereas for DXR peak values were observed in leaves. In C. arabica, both genes were highly expressed in flowers, and CaDXR was upregulated in response to methyl jasmonate. C. arabica DGE data were validated by assessing gene expression in selected organs, and by plants treated with hexanoic acid (Hx) using RT-qPCR. MVK expression was upregulated in roots treated with Hx. CaDXR was downregulated in leaves by Hx treatment in a genotype-specific manner, indicating a differential response to priming. Departamento de Biodiversidade Instituto de Biociências Universidade Estadual Paulista UNESP Laboratório de Biotecnologia Vegetal Empresa Brasileira de Pesquisa Agropecuária (Embrapa-Café) Departamento de Biodiversidade Instituto de Biociências Universidade Estadual Paulista UNESP

MicroRNA target sites are often conserved during evolution and purifying selection to maintain such sites is expected. On the other hand, comparative analyses identified a paucity of microRNA target sites in coexpressed transcripts, and novel target sites can potentially be deleterious. We proposed that selection against novel target sites pervasive. The analysis of derived allele frequencies revealed that, when the derived allele is a target site, the proportion of nontarget sites is higher than expected, particularly for highly expressed microRNAs. Thus, new alleles generating novel microRNA target sites can be deleterious and selected against. When we analyzed ancestral target sites, the derived (nontarget) allele frequency does not show statistical support for microRNA target allele conservation. We investigated the joint effects of microRNA conservation and expression and found that selection against microRNA target sites depends mostly on the expression level of the microRNA. We identified microRNA target sites with relatively high levels of population differentiation. However, when we analyze separately target sites in which the target allele is ancestral to the population, the proportion of single-nucleotide polymorphisms with high Fst significantly increases. These findings support that population differentiation is more likely in target sites that are lost than in the gain of new target sites. Our results indicate that selection against novel microRNA target sites is prevalent and, although individual sites may have a weak selective pressure, the overall effect across untranslated regions is not negligible and should be accounted when studying the evolution of genomic sequences.

De novo genes are very important for evolutionary innovation. However, how these genes originate and spread remains largely unknown. To better understand this, we rigorously searched for de novo genes in Saccharomyces cerevisiae S288C and examined their spread and fixation in the population. Here, we identified 84 de novo genes in S. cerevisiae S288C since the divergence with their sister groups. Transcriptome and ribosome profiling data revealed at least 8 (10%) and 28 (33%) de novo genes being expressed and translated only under specific conditions, respectively. DNA microarray data, based on 2-fold change, showed that 87% of the de novo genes are regulated during various biological processes, such as nutrient utilization and sporulation. Our comparative and evolutionary analyses further revealed that some factors, including single nucleotide polymorphism (SNP)/indel mutation, high GC content, and DNA shuffling, contribute to the birth of de novo genes, while domestication and natural selection drive the spread and fixation of these genes. Finally, we also provide evidence suggesting the possible parallel origin of a de novo gene between S. cerevisiae and Saccharomyces paradoxus. Together, our study provides several new insights into the origin and spread of de novo genes., IMPORTANCE Emergence of de novo genes has occurred in many lineages during evolution, but the birth, spread, and function of these genes remain unresolved. Here we have searched for de novo genes from Saccharomyces cerevisiae S288C using rigorous methods, which reduced the effects of bad annotation and genomic gaps on the identification of de novo genes. Through this analysis, we have found 84 new genes originating de novo from previously noncoding regions, 87% of which are very likely involved in various biological processes. We noticed that 10% and 33% of de novo genes were only expressed and translated under specific conditions, therefore, verification of de novo genes through transcriptome and ribosome profiling, especially from limited expression data, may underestimate the number of bona fide new genes. We further show that SNP/indel mutation, high GC content, and DNA shuffling could be involved in the birth of de novo genes, while domestication and natural selection drive the spread and fixation of these genes. Finally, we provide evidence suggesting the possible parallel origin of a new gene.

A now classical argument for the marginal thermodynamic stability of proteins explains the distribution of observed protein stabilities as a consequence of an entropic pull in protein sequence space. In particular, most sequences that are sufficiently stable to fold will have stabilities near the folding threshold. Here we extend this argument to consider its predictions for epistatic interactions for the effects of mutations on the free energy of folding. Although there is abundant evidence to indicate that the effects of mutations on the free energy of folding are nearly additive and conserved over evolutionary time, we show that these observations are compatible with the hypothesis that a non-additive contribution to the folding free energy is essential for observed proteins to maintain their native structure. In particular through both simulations and analytical results, we show that even very small departures from additivity are sufficient to drive this effect.

Significance How mutation and selection co-determine the course of cancer evolution remains an open, fundamental question. We construct a mutation-selection phase diagram, using tumor mutation load (ML) and selection strength (dN/dS) as key variables, and assess their association with clinical outcome. The results reveal a biphasic evolutionary regime whereby beyond a critical ML, tumor fitness decreases with the number of mutations, although the proteome evolves near neutrality—that is, without strong selection. Deviations from neutrality at extreme ML show how positive selection (at low ML) and purifying selection (at high ML) may act to maintain tumor fitness. These results corroborate the existence of a critical state in cancer evolution predicted by theory and have fundamental and likely clinical implications., How mutation and selection determine the fitness landscape of tumors and hence clinical outcome is an open fundamental question in cancer biology, crucial for the assessment of therapeutic strategies and resistance to treatment. Here we explore the mutation-selection phase diagram of 6,721 tumors representing 23 cancer types by quantifying the overall somatic point mutation load (ML) and selection (dN/dS) in the entire proteome of each tumor. We show that ML strongly correlates with patient survival, revealing two opposing regimes around a critical point. In low-ML cancers, a high number of mutations indicates poor prognosis, whereas high-ML cancers show the opposite trend, presumably due to mutational meltdown. Although the majority of cancers evolve near neutrality, deviations are observed at extreme MLs. Melanoma, with the highest ML, evolves under purifying selection, whereas in low-ML cancers, signatures of positive selection are observed, demonstrating how selection affects tumor fitness. Moreover, different cancers occupy specific positions on the ML–dN/dS plane, revealing a diversity of evolutionary trajectories. These results support and expand the theory of tumor evolution and its nonlinear effects on survival.

To be able to survive,Helicobacter pylorimust adhere to the gastric epithelial cells of its human host. For this purpose, the bacterium employs an array of adhesins, for example, AlpA. The adhesin AlpA has been proposed as a major adhesin because of its critical role in human stomach colonization. Therefore, understanding how AlpA evolved could be important for the development of new diagnostic strategies. However, the genetic variation and microevolutionary patterns ofalpAhave not been described in Colombia. The study aim was to describe the variation patterns and microevolutionary process ofalpAin Colombian clinical isolates ofH. pylori. The existing polymorphisms, which are deviations from the neutral model of molecular evolution, and the genetic differentiation of thealpAgene from Colombian clinical isolates ofH. pyloriwere determined. The analysis shows that gene conversion and purifying selection have shaped the evolution of three different variants ofalpAin Colombia.