Selection scan in Native Americans of Mexico identifies FADS2 rs174616: Evidence of gene-diet interactions affecting lipid levels and Delta-6-desaturase activity.

Searching for positive selection signals across genomes has identified functional genetic variants responding to environmental change. In Native Americans of Mexico, we used the fixation index (Fst) and population branch statistic (PBS) to identify SNPs suggesting positive selection. The 103 most differentiated SNPs were tested for associations with metabolic traits, the most significant association was FADS2/rs174616 with body mass index (BMI). This variant lies within a linkage disequilibrium (LD) block independent of previously reported FADS selection signals and has not been clearly associated with metabolic phenotypes. We tested this variant in two independent cohorts with cardiometabolic data. In the Genetics of Atherosclerotic Disease (GEA) cohort, the derived allele (T) was associated with increased BMI, lower LDL-C levels and a decreased risk of subclinical atherosclerosis in women. Significant gene-diet interactions affected lipid, apolipoprotein and adiponectin levels with differences according to sex, involving mainly total and complex dietary carbohydrate%. In the Genotype-related Effects of PUFA trial, the derived allele was associated with lower Δ-6 desaturase activity and erythrocyte membrane dihomo-gamma-linolenic acid (DGLA) levels, and with increased Δ-5 desaturase activity and eicosapentaenoic acid levels. This variant interacted with dietary carbohydrate% affecting Δ-6 desaturase activity. Notably, the relationship of DGLA and other erythrocyte membrane LC-PUFA indices with HOMA-IR differed according to rs174616 genotype, which has implications regarding how these indices should be interpreted. In conclusion, this observational study identified rs174616 as a signal suggesting selection in an independent linkage disequilibrium block, was associated with cardiometabolic and erythrocyte measurements of LC-PUFA in two independent Mexican cohorts and showed significant gene-diet interactions.

Genetic Signatures of Positive Selection in Human Populations Adapted to High Altitude in Papua New Guinea.

Papua New Guinea (PNG) hosts distinct environments mainly represented by the ecoregions of the Highlands and Lowlands that display increased altitude and a predominance of pathogens, respectively. Since its initial peopling approximately 50,000 years ago, inhabitants of these ecoregions might have differentially adapted to the environmental pressures exerted by each of them. However, the genetic basis of adaptation in populations from these areas remains understudied. Here, we investigated signals of positive selection in 62 highlanders and 43 lowlanders across 14 locations in the main island of PNG using whole-genome genotype data from the Oceanian Genome Variation Project (OGVP) and searched for signals of positive selection through population differentiation and haplotype-based selection scans. Additionally, we performed archaic ancestry estimation to detect selection signals in highlanders within introgressed regions of the genome. Among highland populations we identified candidate genes representing known biomarkers for mountain sickness (SAA4, SAA1, PRDX1, LDHA) as well as candidate genes of the Notch signaling pathway (PSEN1, NUMB, RBPJ, MAML3), a novel proposed pathway for high altitude adaptation in multiple organisms. We also identified candidate genes involved in oxidative stress, inflammation, and angiogenesis, processes inducible by hypoxia, as well as in components of the eye lens and the immune response. In contrast, candidate genes in the lowlands are mainly related to the immune response (HLA-DQB1, HLA-DQA2, TAAR6, TAAR9, TAAR8, RNASE4, RNASE6, ANG). Moreover, we find two candidate regions to be also enriched with archaic introgressed segments, suggesting that archaic admixture has played a role in the local adaptation of PNG populations.

Genomic decoding of Theobroma grandiflorum (cupuassu) at chromosomal scale: evolutionary insights for horticultural innovation.

BACKGROUND: Theobroma grandiflorum (Malvaceae), known as cupuassu, is a tree indigenous to the Amazon basin, valued for its large fruits and seed pulp, contributing notably to the Amazonian bioeconomy. The seed pulp is utilized in desserts and beverages, and its seed butter is used in cosmetics. Here, we present the sequenced telomere-to-telomere genome of cupuassu, disclosing its genomic structure, evolutionary features, and phylogenetic relationships within the Malvaceae family. FINDINGS: The cupuassu genome spans 423 Mb, encodes 31,381 genes distributed in 10 chromosomes, and exhibits approximately 65% gene synteny with the Theobroma cacao genome, reflecting a conserved evolutionary history, albeit punctuated with unique genomic variations. The main changes are pronounced by bursts of long-terminal repeat retrotransposons at postspecies divergence, retrocopied and singleton genes, and gene families displaying distinctive patterns of expansion and contraction. Furthermore, positively selected genes are evident, particularly among retained and dispersed tandem and proximal duplicated genes associated with general fruit and seed traits and defense mechanisms, supporting the hypothesis of potential episodes of subfunctionalization and neofunctionalization following duplication, as well as impact from distinct domestication process. These genomic variations may underpin the differences observed in fruit and seed morphology, ripening, and disease resistance between cupuassu and the other Malvaceae species. CONCLUSIONS: The cupuassu genome offers a foundational resource for both breeding improvement and conservation biology, yielding insights into the evolution and diversity within the genus Theobroma.

Comparative transcriptome analysis reveals lineage- and environment-specific adaptations in cacti from the Brazilian Atlantic Forest.

MAIN CONCLUSION: Natural selection influenced adaptive divergence between Cereus fernambucensis and Cereus insularis, revealing key genes governing abiotic stress responses and supporting neoteny in C. insularis. Uncovering the molecular mechanisms driving adaptive divergence in traits related to habitat adaptation remains a central challenge. In this study, we focused on the cactus clade, which includes Cereus sericifer F.Ritter, Cereus fernambucensis Lem., and Cereus insularis Hemsley. These allopatric species inhabit distinct relatively drier regions within the Brazilian Atlantic Forest, each facing unique abiotic conditions. We leveraged whole transcriptome data and abiotic variables datasets to explore lineage-specific and environment-specific adaptations in these species. Employing comparative phylogenetic methods, we identified genes under positive selection (PSG) and examined their association with non-synonymous genetic variants and abiotic attributes through a PhyloGWAS approach. Our analysis unveiled signatures of selection in all studied lineages, with C. fernambucensis northern populations and C. insularis showing the most PSGs. These PSGs predominantly govern abiotic stress regulation, encompassing heat tolerance, UV stress response, and soil salinity adaptation. Our exclusive observation of gene expression tied to early developmental stages in C. insularis supports the hypothesis of neoteny in this species. We also identified genes associated with abiotic variables in independent lineages, suggesting their role as environmental filters on genetic diversity. Overall, our findings suggest that natural selection played a pivotal role in the geographic range of these species in response to environmental and biogeographic transitions.

Analyses and Insights into Genetic Reassortment and Natural Selection as Key Drivers of Piscine orthoreovirus Evolution.

Piscine orthoreovirus (PRV) is a pathogen that causes heart and skeletal muscle inflammation in Salmo salar and has also been linked to circulatory disorders in other farmed salmonids, such as Oncorhynchus kisutch and Oncorhynchus mykiss. The virus has a segmented, double-stranded RNA genome, which makes it possible to undergo genetic reassortment and increase its genomic diversity through point mutations. In this study, genetic reassortment in PRV was assessed using the full genome sequences available in public databases. This study used full genome sequences that were concatenated and genome-wide reassortment events, and phylogenetic analyses were performed using the recombination/reassortment detection program version 5 (RDP5 V 5.5) software. Additionally, each segment was aligned codon by codon, and overall mean distance and selection was tested using the Molecular Evolutionary Genetics Analysis X software, version 10.2 (MEGA X version 10.2). The results showed that there were 17 significant reassortment events in 12 reassortant sequences, involving genome exchange between low and highly virulent genotypes. PRV sequences from different salmonid host species did not appear to limit the reassortment. This study found that PRV frequently undergoes reassortment events to increase the diversity of its segmented genome, leading to antigenic variation and increased virulence. This study also noted that to date, no reassortment events have been described between PRV-1 and PRV-3 genotypes. However, the number of complete genomic sequences within each genotype is uneven. This is important because PRV-3 induces cross-protection against PRV-1, making it a potential vaccine candidate.

Chloroplast Genome Provides Insights into Molecular Evolution and Species Relationship of Fleabanes (Erigeron: Tribe Astereae, Asteraceae) in the Juan Fernández Islands, Chile.

Erigeron represents the third largest genus on the Juan Fernández Islands, with six endemic species, five of which occur exclusively on the younger Alejandro Selkirk Island with one species on both islands. While its continental sister species is unknown, Erigeron on the Juan Fernández Islands appears to be monophyletic and most likely evolved from South American progenitor species. We characterized the complete chloroplast genomes of five Erigeron species, including accessions of E. fernandezia and one each from Alejandro Selkirk and Robinson Crusoe Islands, with the purposes of elucidating molecular evolution and phylogenetic relationships. We found highly conserved chloroplast genomes in size, gene order and contents, and further identified several mutation hotspot regions. In addition, we found two positively selected chloroplast genes (ccsA and ndhF) among species in the islands. The complete plastome sequences confirmed the monophyly of Erigeron in the islands and corroborated previous phylogenetic relationships among species. New findings in the current study include (1) two major lineages, E. turricola-E. luteoviridis and E. fernandezia-E. ingae-E. rupicola, (2) the non-monophyly of E. fernandezia occurring on the two islands, and (3) the non-monophyly of the alpine species E. ingae complex.

Opsin diversity and evolution in the Elateroidea superfamily: Insights from transcriptome data.

Vision plays a vital biological role in organisms, which depends on the visual pigment molecules (opsin plus chromophore). The expansion or reduction of spectral channels in the organisms is determined by distinct opsin classes and copy numbers resulting from duplication or loss. Within Coleoptera, the superfamily Elateroidea exhibits a great diversity of morphological and physiological characteristics, such as bioluminescence, making this group an important model for opsin studies. While molecular and physiological studies have been conducted in Lampyridae and Elateridae, other families remain unexplored. Here, we reused transcriptome datasets from Elateroidea species, including members of Elateridae, Lampyridae, Phengodidae, Rhagophthalmidae, Cantharidae, and Lycidae, to detect the diversity of putative opsin genes in this superfamily. In addition, we tested the signature of sites under positive selection in both ultraviolet (UV)- and long-wavelength (LW)-opsin classes. Although the visual system in Elateroidea is considered simple, we observed events of duplication in LW- and UV-opsin, as well as the absence of UV-opsin in distinct families, such as larval Phengodidae individuals. We detected different copies of LW-opsins that were highly expressed in the eyes of distinct tribes of fireflies, indicating the possible selection of each copy during the evolution of the sexual mating to avoid spectrum overlapping. In Elateridae, we found that the bioluminescent species had a distinct LW-opsin copy compared with the non-bioluminescent species, suggesting events of duplication and loss. The signature of positive selection showed only one residue associated with the chromophore binding site in the Elateroidea, which may produce a bathochromic shift in the wavelength absorption spectra in this family. Overall, this study brings important content and fills gaps regarding opsin evolution in Elateroidea.

Molecular evolution and structural analyses of the spike glycoprotein from Brazilian SARS-CoV-2 genomes: the impact of selected mutations.

The COVID-19 pandemic caused by SARS-CoV-2 has reached by February 2022 more than 380 million cases and 5.5 million deaths worldwide since its beginning in late 2019, leading to enhanced concern in the scientific community and the general population. One of the most important pieces of this host-pathogen interaction is the spike protein, which binds to the hACE2 cell receptor, mediates the membrane fusion and is the major target of neutralizing antibodies against SARS-CoV-2. The multiple amino acid substitutions observed in this region, specially in RBD have enhanced the hACE2 binding affinity and led to several modifications in the mechanisms of SARS-CoV-2 pathogenesis, improving the viral fitness and/or promoting immune evasion, with potential impact in the vaccine development. In this work, we identified 48 sites under selective pressures, 17 of them with the strongest evidence by the HyPhy tests, including VOC related mutation sites 138, 142, 222, 262, 484, 681, and 845, among others. The coevolutionary analysis identified 28 sites found not to be conditionally independent, such as E484K-N501Y. The molecular dynamics and free energy estimates showed the structural stabilizing effect and the higher impact of E484K for enhanced binding affinity between the spike RBD and hACE2 in P.1 and P.2 lineages (specially with L452V). Structural changes were also identified in the hACE molecule when interacting with B.1.1.7 RDB. Despite some destabilizing substitutions, a stabilizing effect was identified for the majority of the positively selected mutations.Communicated by Ramaswamy H. Sarma.

Comparative genomics uncovers the evolutionary history, demography, and molecular adaptations of South American canids.

The remarkable radiation of South American (SA) canids produced 10 extant species distributed across diverse habitats, including disparate forms such as the short-legged, hypercarnivorous bush dog and the long-legged, largely frugivorous maned wolf. Despite considerable research spanning nearly two centuries, many aspects of their evolutionary history remain unknown. Here, we analyzed 31 whole genomes encompassing all extant SA canid species to assess phylogenetic relationships, interspecific hybridization, historical demography, current genetic diversity, and the molecular bases of adaptations in the bush dog and maned wolf. We found that SA canids originated from a single ancestor that colonized South America 3.9 to 3.5 Mya, followed by diversification east of the Andes and then a single colonization event and radiation of Lycalopex species west of the Andes. We detected extensive historical gene flow between recently diverged lineages and observed distinct patterns of genomic diversity and demographic history in SA canids, likely induced by past climatic cycles compounded by human-induced population declines. Genome-wide scans of selection showed that disparate limb proportions in the bush dog and maned wolf may derive from mutations in genes regulating chondrocyte proliferation and enlargement. Further, frugivory in the maned wolf may have been enabled by variants in genes associated with energy intake from short-chain fatty acids. In contrast, unique genetic variants detected in the bush dog may underlie interdigital webbing and dental adaptations for hypercarnivory. Our analyses shed light on the evolution of a unique carnivoran radiation and how it was shaped by South American topography and climate change.

Genomic Signatures of Divergent Ecological Strategies in a Recent Radiation of Neotropical Wild Cats.

Ecological differentiation among diverging species is an important component of the evolutionary process and can be investigated in rapid and recent radiations. Here, we use whole genome sequences of five species from the genus Leopardus, a recently diversified Neotropical lineage with species bearing distinctive morphological, ecological, and behavioral features, to investigate genome-wide diversity, comparative demographic history and signatures of positive selection. Our results show that divergent ecological strategies are reflected in genomic features, for example a generalist species shows historically larger effective population size and higher heterozygosity than habitat specialists. The demographic history of these cats seems to have been jointly driven by climate fluctuations and habitat specialization, with different ecological adaptations leading to distinct trajectories. Finally, a gene involved in vertebrate retinal neurogenesis (POU4F2) was found to be under positive selection in the margay, a cat with notoriously large eyes that are likely associated with its nocturnal and arboreal specializations.

What We Talk About When We Talk About "Junk DNA".

"Junk DNA" is a popular yet controversial concept that states that organisms carry in their genomes DNA that has no positive impact on their fitness. Nonetheless, biochemical functions have been identified for an increasing fraction of DNA elements traditionally seen as "Junk DNA". These findings have been interpreted as fundamentally undermining the "Junk DNA" concept. Here, we reinforce previous arguments that this interpretation relies on an inadequate concept of biological function that does not consider the selected effect of a given genomic structure, which is central to the "Junk DNA" concept. Next, we suggest that another (though ignored) confounding factor is that the discussion about biological functions includes two different dimensions: a horizontal, ecological dimension that reflects how a given genomic element affects fitness in a specific time, and a vertical, temporal dimension that reflects how a given genomic element persisted along time. We suggest that "Junk DNA" should be used exclusively relative to the horizontal dimension, while for the vertical dimension, we propose a new term, "Spam DNA", that reflects the fact that a given genomic element may persist in the genome even if not selected for on their origin. Importantly, these concepts are complementary. An element can be both "Spam DNA" and "Junk DNA", and "Spam DNA" can also be recruited to perform evolved biological functions, as illustrated in processes of exaptation or constructive neutral evolution.

Evidence of episodic positive selection in Corynebacterium diphtheriae complex of species and its implementations in identification of drug and vaccine targets.

BACKGROUND: Within the pathogenic bacterial species Corynebacterium genus, six species that can produce diphtheria toxin (C. belfantii, C. diphtheriae, C. pseudotuberculosis, C. rouxii, C. silvaticum and C. ulcerans) form a clade referred to as the C. diphtheria complex. These species have been found in humans and other animals, causing diphtheria or other diseases. Here we show the results of a genome scale analysis to identify positive selection in protein-coding genes that may have resulted in the adaptations of these species to their ecological niches and suggest drug and vaccine targets. METHODS: Forty genomes were sampled to represent species, subspecies or biovars of Corynebacterium. Ten phylogenetic groups were tested for positive selection using the PosiGene pipeline, including species and biovars from the C. diphtheria complex. The detected genes were tested for recombination and had their sequences alignments and homology manually examined. The final genes were investigated for their function and a probable role as vaccine or drug targets. RESULTS: Nineteen genes were detected in the species C. diphtheriae (two), C. pseudotuberculosis (10), C. rouxii (one), and C. ulcerans (six). Those were found to be involved in defense, translation, energy production, and transport and in the metabolism of carbohydrates, amino acids, nucleotides, and coenzymes. Fourteen were identified as essential genes, and six as virulence factors. Thirteen from the 19 genes were identified as potential drug targets and four as potential vaccine candidates. These genes could be important in the prevention and treatment of the diseases caused by these bacteria.

The genome of a thorny species: comparative genomic analysis among South and North American Cactaceae.

MAIN CONCLUSION: The first South American cactus nuclear genome assembly associated with comparative genomic analyses provides insights into nuclear and plastid genomic features, such as size, transposable elements, and metabolic processes related to cactus development. Here, we assembled the partial genome, plastome, and transcriptome of Cereus fernambucensis (Cereeae, Cactaceae), a representative species of the South American core Cactoideae. We accessed other genomes and transcriptomes available for cactus species to compare the heterozygosity level, genome size, transposable elements, orthologous genes, and plastome structure. These estimates were obtained from the literature or using the same pipeline adopted for C. fermabucensis. In addition to the C. fernambucensis plastome, we also performed de novo plastome assembly of Pachycereus pringlei, Stenocereus thurberi, and Pereskia humboldtii based on the sequences available in public databases. We estimated a genome size of ~ 1.58 Gb for C. fernambucensis, the largest genome among the compared species. The genome heterozygosity was 0.88% in C. fernambucensis but ranged from 0.36 (Carnegiea gigantea) to 17.4% (Lophocereus schottii) in the other taxa. The genome lengths of the studied cacti are constituted by a high amount of transposable elements, ranging from ~ 57 to ~ 67%. Putative satellite DNAs are present in all species, excepting C. gigantea. The plastome of C. fernambucensis was ~ 104 kb, showing events of translocation, inversion, and gene loss. We observed a low number of shared unique orthologs, which may suggest gene duplication events and the simultaneous expression of paralogous genes. We recovered 37 genes that have undergone positive selection along the Cereus branch that are associated with different metabolic processes, such as improving photosynthesis during drought stress and nutrient absorption, which may be related to the adaptation to xeric areas of the Neotropics.

Adaptive evolution of Moniliophthora PR-1 proteins towards its pathogenic lifestyle.

BACKGROUND: Plant pathogenesis related-1 (PR-1) proteins belong to the CAP superfamily and have been characterized as markers of induced defense against pathogens. Moniliophthora perniciosa and Moniliophthora roreri are hemibiotrophic fungi that respectively cause the witches' broom disease and frosty pod rot in Theobroma cacao. Interestingly, a large number of plant PR-1-like genes are present in the genomes of both species and many are up-regulated during the biotrophic interaction. In this study, we investigated the evolution of PR-1 proteins from 22 genomes of Moniliophthora isolates and 16 other Agaricales species, performing genomic investigation, phylogenetic reconstruction, positive selection search and gene expression analysis. RESULTS: Phylogenetic analysis revealed conserved PR-1 genes (PR-1a, b, d, j), shared by many Agaricales saprotrophic species, that have diversified in new PR-1 genes putatively related to pathogenicity in Moniliophthora (PR-1f, g, h, i), as well as in recent specialization cases within M. perniciosa biotypes (PR-1c, k, l) and M. roreri (PR-1n). PR-1 families in Moniliophthora with higher evolutionary rates exhibit induced expression in the biotrophic interaction and positive selection clues, supporting the hypothesis that these proteins accumulated adaptive changes in response to host-pathogen arms race. Furthermore, although previous work showed that MpPR-1 can detoxify plant antifungal compounds in yeast, we found that in the presence of eugenol M. perniciosa differentially expresses only MpPR-1e, k, d, of which two are not linked to pathogenicity, suggesting that detoxification might not be the main function of most MpPR-1. CONCLUSIONS: Based on analyses of genomic and expression data, we provided evidence that the evolution of PR-1 in Moniliophthora was adaptive and potentially related to the emergence of the parasitic lifestyle in this genus. Additionally, we also discuss how fungal PR-1 proteins could have adapted from basal conserved functions to possible roles in fungal pathogenesis.

HLA-G in Mayas from Yucatan: An evolutionary approach.

HLA-G allele frequencies were studied in Yucatán (Mexico) Maya Amerindians by a direct exon DNA sequencing technique. It is described that Mayas are probably one of the first populations together with Olmecs that populated Meso America and that important HLA genetic differences between Mexican and Guatemalan Mayas support that Maya languages were imposed to several neighbouring Amerindian groups. HLA-G*01:01:02, HLA-G*01:01:01 and HLA-G*01:04:01 are the most frequent alleles in this population. It is remarkable that HLA-G*01:05N allele was not found in the population in accordance with similar results found in another Amerindians. Also, protein allele HLA-G*01:04 frequency is found not to differ to those found in another far or close living Amerindians in contrast to other World populations. It seems that while high HLA-G*01:05N frequency is found in Iran and Middle East populations, probably where this allele appeared within an ancestral HLA-A*19 group of alleles haplotype and it is maintained by unknown evolutionary forces, Amerindians do not have a high frequency because a founder effect or because required natural evolutionary forces do not exist in America. Finally, we believe useful to study HLA-G evolution for its physiopathology understanding in addition to the many papers on statistics on HLA-G and in vitro models that are yearly published.

Extremely low nucleotide diversity among thirty-six new chloroplast genome sequences from Aldama (Heliantheae, Asteraceae) and comparative chloroplast genomics analyses with closely related genera.

Aldama (Heliantheae, Asteraceae) is a diverse genus in the sunflower family. To date, nearly 200 Asteraceae chloroplast genomes have been sequenced, but the plastomes of Aldama remain undescribed. Plastomes in Asteraceae usually show little sequence divergence, consequently, our hypothesis is that species of Aldama will be overall conserved. In this study, we newly sequenced 36 plastomes of Aldama and of five species belonging to other Heliantheae genera selected as outgroups (i.e., Dimerostemma asperatum, Helianthus tuberosus, Iostephane heterophylla, Pappobolus lanatus var. lanatus, and Tithonia diversifolia). We analyzed the structure and gene content of the assembled plastomes and performed comparative analyses within Aldama and with other closely related genera. As expected, Aldama plastomes are very conserved, with the overall gene content and orientation being similar in all studied species. The length of the plastome is also consistent and the junction between regions usually contain the same genes and have similar lengths. A large ∼20 kb and a small ∼3 kb inversion were detected in the Large Single Copy (LSC) regions of all assembled plastomes, similarly to other Asteraceae species. The nucleotide diversity is very low, with only 1,509 variable sites in 127,466 bp (i.e., 1.18% of the sites in the alignment of 36 Aldama plastomes, with one of the IRs removed, is variable). Only one gene, rbcL, shows signatures of positive selection. The plastomes of the selected outgroups feature a similar gene content and structure compared to Aldama and also present the two inversions in the LSC region. Deletions of different lengths were observed in the gene ycf2. Multiple SSRs were identified for the sequenced Aldama and outgroups. The phylogenetic analysis shows that Aldama is not monophyletic due to the position of the Mexican species A. dentata. All Brazilian species form a strongly supported clade. Our results bring new understandings into the evolution and diversity of plastomes at the species level.

Detection of genes positively selected in Cuban Anolis lizards that naturally inhabit hot and open areas and currently thrive in urban areas.

Species of Anolis lizards of the West Indies that naturally inhabit hot and open areas also tend to thrive in urban areas. In this study, transcriptome was sequenced for nine species of Cuban Anolis lizards that are closely related to each other, but inhabit different thermal microhabitats. Using PAML and HyPhy software, we attempted to identify genes and amino acid sites under positive selection in the common ancestral branch of A. porcatus and A. allisoni, and the branch of A. sagrei, which inhabit hot and open areas, and thrive in urban areas. Although there were no genes where positive selection was commonly detected on both of the tested branches, positive selection was detected in genes involved in the stress response (e.g., DNA damage and oxidative stress) and cardiac function, which could be related to adaptive evolution of tolerance to heat or ultraviolet radiation, on both branches. These findings suggest that adaptive evolution of the response to stress caused by heat or ultraviolet radiation might have occurred in ancestors of Anolis species inhabiting hot and open areas and might be related to the current thriving in urban areas of them.

Examining Signatures of Natural Selection in Antifungal Resistance Genes Across Aspergillus Fungi.

Certain Aspergillus fungi cause aspergillosis, a set of diseases that typically affect immunocompromised individuals. Most cases of aspergillosis are caused by Aspergillus fumigatus, which infects millions of people annually. Some closely related so-called cryptic species, such as Aspergillus lentulus, can also cause aspergillosis, albeit at lower frequencies, and they are also clinically relevant. Few antifungal drugs are currently available for treating aspergillosis and there is increasing worldwide concern about the presence of antifungal drug resistance in Aspergillus species. Furthermore, isolates from both A. fumigatus and other Aspergillus pathogens exhibit substantial heterogeneity in their antifungal drug resistance profiles. To gain insights into the evolution of antifungal drug resistance genes in Aspergillus, we investigated signatures of positive selection in 41 genes known to be involved in drug resistance across 42 susceptible and resistant isolates from 12 Aspergillus section Fumigati species. Using codon-based site models of sequence evolution, we identified ten genes that contain 43 sites with signatures of ancient positive selection across our set of species. None of the sites that have experienced positive selection overlap with sites previously reported to be involved in drug resistance. These results identify sites that likely experienced ancient positive selection in Aspergillus genes involved in resistance to antifungal drugs and suggest that historical selective pressures on these genes likely differ from any current selective pressures imposed by antifungal drugs.

Evolution and epidemiologic dynamics of dengue virus in Nicaragua during the emergence of chikungunya and Zika viruses.

Arthropod-borne viruses (arboviruses) comprise a significant and ongoing threat to human health, infecting hundreds of millions annually. Three such arboviruses include circumtropical dengue, Zika, and chikungunya viruses, exhibiting continuous emergence primarily via Aedes mosquito vectors. Nicaragua has experienced endemic dengue virus (DENV) transmission involving multiple serotypes since 1985, with chikungunya virus (CHIKV) reported in 2014-2015, followed by Zika virus (ZIKV) first reported in 2016. In order to identify patterns of genetic variation and selection pressures shaping the evolution of co-circulating DENV serotypes in light of the arrival of CHIKV and ZIKV, we employed whole-genome sequencing on an Illumina MiSeq platform of random-amplified total RNA libraries to characterize 42 DENV low-passage isolates, derived from viremic patients in Nicaragua between 2013 and 2016. Our approach also revealed clinically undetected co-infections with CHIKV. Of the three DENV serotypes (1, 2, and 3) co-circulating during our study, we uncovered distinct patterns of evolution using comparative phylogenetic inference. DENV-1 genetic variation was structured into two distinct co-circulating lineages with no evidence of positive selection in the origins of either lineage, suggesting they are equally fit. In contrast, the evolutionary history of DENV-2 was marked by positive selection, and a unique, divergent lineage correlated with high epidemic potential emerged in 2015 to drive an outbreak in 2016. DENV-3 genetic variation remained unstructured into lineages throughout the period of study. Thus, this study reveals insights into evolutionary and epidemiologic trends exhibited during the circulation of multiple arboviruses in Nicaragua.

Chloroplast genome assembly of Handroanthus impetiginosus: comparative analysis and molecular evolution in Bignoniaceae.

MAIN CONCLUSION: Bignoniaceae species have conserved chloroplast structure, with hotspots of nucleotide diversity. Several genes are under positive selection, and can be targets for evolutionary studies. Bignoniaceae is one of the most species-rich family of woody plants in Neotropical seasonally dry forests. Here we report the assembly of Handroanthus impetiginosus chloroplast genome and evolutionary comparative analyses of ten Bignoniaceae species representing the genera for which whole-genome chloroplast sequences were available. The chloroplast genome of H. impetiginosus is 159,462 bp in size and has a similar structure compared to the other nine species. The total number of genes was slightly variable amongst the Bignoniaceae, ranging from 124 in H. impetiginosus to 144 in Anemopaegma acutifolium. The inverted repeat (IR) size was variable, ranging from 24,657 bp (Tecomaria capensis) to 40,481 bp (A. acutifolium), due to the contraction and retraction at its boundaries. However, gene boundaries were very similar among the ten species. We found 98 forward and palindromic dispersed repeats, and 85 simple sequence repeats (SSRs). In general, chloroplast sequences were highly conserved, with few nucleotide diversity hotspots in the genes accD, clpP, rpoA, ycf1, ycf2. The phylogenetic analysis based on 77 coding genes was highly consistent with Angiosperm Phylogeny Group (APG) IV. Our results also indicate that most genes are under negative selection or neutral evolution. We found no evidence of branch-site selection, implying that H. impetiginosus is not evolving faster than the other species analyzed, notwithstanding we found site positive selection signal in several genes. These genes can provide targets for evolutionary studies in Bignoniaceae and Lamiales species.