Detection of evolutionary conserved and accelerated genomic regions related to adaptation to thermal niches in Anolis lizards

Understanding the genetic basis for adapting to thermal environments is important due to serious effects of global warming on ectothermic species. Various genes associated with thermal adaptation in lizards have been identified mainly focusing on changes in gene expression or the detection of positively selected genes using coding regions. Only a few comprehensive genome-wide analyses have included noncoding regions. This study aimed to identify evolutionarily conserved and accelerated genomic regions using whole genomes of eight Anolis lizard species that have repeatedly adapted to similar thermal environments in multiple lineages. Evolutionarily conserved genomic regions were extracted as regions with overall sequence conservation (regions with fewer base substitutions) across all lineages compared with the neutral model. Genomic regions that underwent accelerated evolution in the lineage of interest were identified as those with more base substitutions in the target branch than in the entire background branch. Conserved elements across all branches were relatively abundant in "intergenic" genomic regions among noncoding regions. Accelerated regions (ARs) of each lineage contained a significantly greater proportion of noncoding RNA genes than the entire multiple alignment. Common genes containing ARs within 5 kb of their vicinity in lineages with similar thermal habitats were identified. Many genes associated with circadian rhythms and behavior were found in hot-open and cool-shaded habitat lineages. These genes might play a role in contributing to thermal adaptation and assist future studies examining the function of genes involved in thermal adaptation via genome editing.

Ongoing declines for the world's amphibians in the face of emerging threats

Systematic assessments of species extinction risk at regular intervals are necessary for informing conservation action1,2. Ongoing developments in taxonomy, threatening processes and research further underscore the need for reassessment3,4. Here we report the findings of the second Global Amphibian Assessment, evaluating 8,011 species for the International Union for Conservation of Nature Red List of Threatened Species. We find that amphibians are the most threatened vertebrate class (40.7% of species are globally threatened). The updated Red List Index shows that the status of amphibians is deteriorating globally, particularly for salamanders and in the Neotropics. Disease and habitat loss drove 91% of status deteriorations between 1980 and 2004. Ongoing and projected climate change effects are now of increasing concern, driving 39% of status deteriorations since 2004, followed by habitat loss (37%). Although signs of species recoveries incentivize immediate conservation action, scaled-up investment is urgently needed to reverse the current trends.

Draft genome of six Cuban Anolis lizards and insights into genetic changes during their diversification

Background

Detecting genomic variants and their accumulation processes during species diversification and adaptive radiation is important for understanding the molecular and genetic basis of evolution. Anolis lizards in the West Indies are good models for studying evolutionary mechanisms because of the repeated evolution of their morphology and the ecology. We performed de novo genome assembly of six Cuban Anolis lizards with different ecomorphs and thermal habitats (Anolis isolepis, Anolis allisoni, Anolis porcatus, Anolis allogus, Anolis homolechis, and Anolis sagrei). We carried out a comparative analysis of these genome assemblies to investigate the genetic changes that occurred during their diversification.

Results

We reconstructed novel draft genomes with relatively long scaffolds and high gene completeness, with the scaffold N50 ranging from 5.56 to 39.79 Mb and vertebrate Benchmarking Universal Single-Copy Orthologs completeness ranging from 77.5% to 86.9%. Comparing the repeat element compositions and landscapes revealed differences in the accumulation process between Cuban trunk-crown and trunk-ground species and separate expansions of several families of LINE in each Cuban trunk-ground species. Duplicated gene analysis suggested that the proportional differences in duplicated gene numbers among Cuban Anolis lizards may be associated with differences in their habitat ranges. Additionally, Pairwise Sequentially Markovian Coalescent analysis suggested that the effective population sizes of each species may have been affected by Cuba’s geohistory.

Conclusions

We provide draft genomes of six Cuban Anolis lizards and detected species and lineage-specific transposon accumulation and gene copy number changes that may be involved in adaptive evolution. The change processes in the past effective population size was also estimated, and the factors involved were inferred. These results provide new insights into the genetic basis of Anolis lizard diversification and are expected to serve as a stepping stone for the further elucidation of their diversification mechanisms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02086-7.

Rootstock effects on scion gene expression in maritime pine

Pines are the dominant conifers in Mediterranean forests. As long-lived sessile organisms that seasonally have to cope with drought periods, they have developed a variety of adaptive responses. However, during last decades, highly intense and long-lasting drought events could have contributed to decay and mortality of the most susceptible trees. Among conifer species, Pinus pinaster Ait. shows remarkable ability to adapt to different environments. Previous molecular analysis of a full-sib family designed to study drought response led us to find active transcriptional activity of stress-responding genes even without water deprivation in tolerant genotypes. To improve our knowledge about communication between above- and below-ground organs of maritime pine, we have analyzed four graft-type constructions using two siblings as rootstocks and their progenitors, Gal 1056 and Oria 6, as scions. Transcriptomic profiles of needles from both scions were modified by the rootstock they were grafted on. However, the most significant differential gene expression was observed in drought-sensitive Gal 1056, while in drought-tolerant Oria 6, differential gene expression was very much lower. Furthermore, both scions grafted onto drought-tolerant rootstocks showed activation of genes involved in tolerance to abiotic stress, and is most remarkable in Oria 6 grafts where higher accumulation of transcripts involved in phytohormone action, transcriptional regulation, photosynthesis and signaling has been found. Additionally, processes, such as those related to secondary metabolism, were mainly associated with the scion genotype. This study provides pioneering information about rootstock effects on scion gene expression in conifers.

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.

The Chytrid Fungus, Batrachochytrium dendrobatidis, is Widespread Among Cuban Amphibians

The fungus Batrachochytrium dendrobatidis (Bd) is a generalist amphibian pathogen responsible for chytridiomycosis. It was documented for the first time in Cuba in 2007, the apparent cause of the decline in one species of toad. In a recent survey, Bd was reported only for the highlands of Central Cuba. In the present study, we reexamined the geographic distribution and level of impact of Bd in Cuba by conducting an island-wide sampling in 10 localities and collecting skin swabs from 18 species and 28 environmental samples. We report detection of Bd in 60% of sampled sites and in 58% of sampled taxa. We show that Bd is associated with riparian, arboreal and terrestrial species, and it was estimated to occur in approximately 30% of the aquatic habitats we sampled. In addition, we confirmed that a dying individual of the species Eleutherodactylus casparii was severely infected with Bd. We also rise concern about the endanger toad Peltophryne longinasus and about three species of endemic riparian frogs that were not detected during our surveys. This study demonstrates that this pathogen is widespread throughout Cuba and provides relevant evidence to advance our understanding of its detection in amphibians and the aquatic environment in Cuba and about the occurrence of Bd in species with different ecologies. We provide valuable baseline information for Bd risk assessment and decision-making processes to mitigate its negative impact on Cuban amphibians.

Factors affecting interspecific differences in genetic divergence among populations of Anolis lizards in Cuba

BACKGROUND

Geographical patterns and degrees of genetic divergence among populations differ between species, reflecting relative potentials for speciation or cladogenesis and differing capacities for environmental adaptation. Identification of factors that contribute to genetic divergence among populations is important to the understanding of why some species exhibit greater interpopulation genetic divergence. In this study, we calculated the mean pairwise genetic distances among populations as species' average genetic divergence by a phylogeny using nuclear and mitochondrial genes of 303 individuals from 33 Cuban Anolis species and estimated species ages by another phylogeny using nuclear and mitochondrial genes of 51 Cuban and 47 non-Cuban Anolis species. We identified factors that influence species' differences in genetic divergence among 26 species of Anolis lizards from Cuba. Species ages, environmental heterogeneity within species ranges, and ecomorph types were considered as factors affecting average genetic divergences among populations.

RESULTS

The phylogenies presented in this study provide the most comprehensive sampling of Cuban Anolis species to date. The phylogeny showed more conservative evolution of Anolis ecomorphs within Cuba and identified twig anoles as a monophyletic group. Subsequent Phylogenetic Generalized Least Squares (PGLS) analyses showed that species age was positively correlated with species' average genetic divergence among populations.

CONCLUSION

Although previous studies have focused on factors affecting genetic divergence within species, the present study showed for the first time that species differences in genetic divergence could be largely affected by species age.

Microsatellite characterization of Andean races of common bean (Phaseolus vulgaris L.)

The Andean gene pool of common bean (Phaseolus vulgaris L.) has high levels of morphological diversity in terms of seed color and size, growth habit and agro-ecological adaptation, but previously was characterized by low levels of molecular marker diversity. Three races have been described within the Andean gene pool: Chile, Nueva Granada and Peru. The objective of this study was to characterize a collection of 123 genotypes representing Andean bean diversity with 33 microsatellite markers that have been useful for characterizing race structure in common beans. The genotypes were from both the primary center of origin as well as secondary centers of diversity to which Andean beans spread and represented all three races of the gene pool. In addition we evaluated a collection of landraces from Colombia to determine if the Nueva Granada and Peru races could be distinguished in genotypes from the northern range of the primary center. Multiple correspondence analyses of the Andean race representatives identified two predominant groups corresponding to the Nueva Granada and Peru races. Some of the Chile race representatives formed a separate group but several that had been defined previously as from this race grouped with the other races. Gene flow was more notable between Nueva Granada and Peru races than between these races and the Chile race. Among the Colombian genotypes, the Nueva Granada and Peru races were identified and introgression between these two races was especially notable. The genetic diversity within the Colombian genotypes was high, reaffirming the importance of this region as an important source of germplasm. Results of this study suggest that the morphological classification of all climbing beans as Peru race genotypes and all bush beans as Nueva Granada race genotypes is erroneous and that growth habit traits have been mixed in both races, requiring a re-adjustment in the concept of morphological races in Andean beans.

Race structure within the Mesoamerican gene pool of common bean (Phaseolus vulgaris L.) as determined by microsatellite markers

Common bean (Phaseolus vulgaris L.) cultivars are distinguished morphologically, agronomically and ecologically into specific races within each of the two gene pools found for the species (Andean and Mesoamerican). The objective of this study was to describe the race structure of the Mesoamerican gene pool using microsatellite markers. A total of 60 genotypes previously described as pertaining to specific Mesoamerican races as well as two Andean control genotypes were analyzed with 52 markers. A total of 267 bands were generated with an average of 5.1 alleles per marker and 0.297 heterozygosity across all microsatellites. Correspondence analysis identified two major groups equivalent to the Mesoamerica race and a group containing both Durango and Jalisco race genotypes. Two outlying individuals were classified as potentially of the Guatemala race although this race does not have a defined structure and previously classified members of this race were classified with other races. Population structure analysis with K = 1-4 agreed with this classification. The genetic diversity based on Nei's index for the entire set of genotypes was 0.468 while this was highest for the Durango-Jalisco group (0.414), intermediate for race Mesoamerica (0.340) and low for race Guatemala (0.262). Genetic differentiation (G (ST)) between the Mesoamerican races was 0.27 while genetic distance and identity showed race Durango and Jalisco individuals to be closely related with high gene flow (N (m)) both between these two races (1.67) and between races Durango and Mesoamerica (1.58). Observed heterozygosity was low in all the races as would be expected for an inbreeding species. The analysis with microsatellite markers identified subgroups, which agreed well with commercial class divisions, and seed size was the main distinguishing factor between the two major groups identified.

Prevalence of ectopic ventricular activity after left ventricular mass regression

To assess the effects of left ventricular mass reduction on the prevalence of ventricular ectopic activity, we conducted a double-blind, placebo-controlled trial measuring ambulatory 48 h premature ventricular depolarizations in 27 patients with mild-to-moderate hypertension and an increased left ventricular mass index. Data was obtained at baseline and 6 +/- 2 months after randomization to either 25 mg captopril or placebo twice a day. Patients on captopril attained reduction in blood pressure from 167 +/- 11/103 +/- 6 to 136 +/- 10/85 +/- 5 mm Hg (P = .001), left ventricular mass index regression from 149 +/- 17 to 96 +/- 23 g/m2 (P = .001), and ventricular ectopic activity reduction from 413 +/- 172 to 77 +/- 27 ventricular extrasystoles/day (P = .001). Patients on placebo had no significant change in blood pressure (from 162 +/- 11/101 +/- 6 at baseline to 160 +/- 8/100 +/- 8 mm Hg after 6 months; P = NS). In the placebo group left ventricular mass index increased from 155 +/- 40 to 182 +/- 51 g/m2 (P = .01), and ventricular ectopic activity decreased from 634 +/- 293 to 562 +/- 260 ventricular extrasystoles/day (P = NS). Eight out of 14 patients on captopril (57%) and 1 out of 13 patients on placebo (8%) achieved reduction > 85% in ventricular ectopic activity per day (P = .004). Using multivariate logistic regression analysis, left ventricular mass index regression and reduction in systolic blood pressure were the most important correlates for this effect.(ABSTRACT TRUNCATED AT 250 WORDS)

Effects of enalapril on heart failure in hypertensive patients with diastolic dysfunction

Ten hypertensive patients with symptoms of heart failure and normal systolic function but with diastolic dysfunction were treated with 10 mg enalapril twice a day for 9 +/- 3 months to evaluate the effects of this agent alone on heart failure induced by diastolic dysfunction. After therapy, all patients improved and echocardiographic parameters of diastolic dysfunction became normalized. It is concluded that enalapril appears to be useful in the treatment of heart failure in hypertensive patients with normal systolic function and diastolic dysfunction.