Elevational Goldilocks zone underlies the exceptional diversity of a large lizard radiation (Liolaemus; Liolaemidae)

Mountains are among the most biodiverse regions on the planet, and how these landforms shape diversification through the interaction of biological traits and geo-climatic dynamics is integral to understanding global biodiversity. In this study, we investigate the dual roles of climate change and mountain uplift on the evolution of a hyper-diverse radiation, Liolaemus lizards, with a spatially explicit model of diversification using a reconstruction of uplift and paleotemperature in central and southern South America. The diversification model captures a hotspot for Liolaemus around 40°S in lineages with low-dispersal ability and narrow niche breadths. Under the model, speciation rates are highest in low latitudes (<35°S) and mid elevations (~1,000 m), while extinction rates are highest at higher latitudes (>35°S) and higher elevations (>2,000 m). Temperature change through the Cenozoic explained variation in speciation and extinction rates through time and across different elevational bands. Our results point to the conditions of mid elevations being optimal for diversification (i.e., Goldilocks Zone), driven by the combination of (1) a complex topography that facilitates speciation during periods of climatic change, and (2) a relatively moderate climate that enables the persistence of ectothermic lineages and buffers species from extinction.

Evolutionary models demonstrate rapid and adaptive diversification of Australo-Papuan pythons

Lineages may diversify when they encounter available ecological niches. Adaptive divergence by ecological opportunity often appears to follow the invasion of a new environment with open ecological space. This evolutionary process is hypothesized to explain the explosive diversification of numerous Australian vertebrate groups following the collision of the Eurasian and Australian plates 25 Mya. One of these groups is the pythons, which demonstrate their greatest phenotypic and ecological diversity in Australo-Papua (Australia and New Guinea). Here, using an updated and near complete time-calibrated phylogenomic hypothesis of the group, we show that following invasion of this region, pythons experienced a sudden burst of speciation rates coupled with multiple instances of accelerated phenotypic evolution in head and body shape and body size. These results are consistent with adaptive radiation theory with an initial rapid niche-filling phase and later slow-down approaching niche saturation. We discuss these findings in the context of other Australo-Papuan adaptive radiations and the importance of incorporating adaptive diversification systems that are not extraordinarily species-rich but ecomorphologically diverse to understand how biodiversity is generated.

Ontogenetic drivers of morphological evolution in monitor lizards and allies (Squamata: Paleoanguimorpha), a clade with extreme body size disparity

Background

Heterochrony, change in the rate or timing of development, is thought to be one of the main drivers of morphological evolution, and allometry, trait scaling patterns imposed by size, is traditionally thought to represent an evolutionary constraint. However, recent studies suggest that the ontogenetic allometric trajectories describing how organisms change as they grow may be labile and adaptive. Here we investigated the role of postnatal ontogenetic development in the morphological diversification of Paleoanguimorpha, the monitor lizards and allies, a clade with extreme body size disparity. We obtained linear and geometric morphometric data for more than 1,600 specimens belonging to three families and 60 species, representing ~ 72% of extant paleoanguimorph diversity. We used these data to undertake one of the largest comparative studies of ontogenetic allometry to date.

Results

Heterochrony is likely dictating morphological divergence at shallow evolutionary scales, while changes in the magnitude and direction of ontogenetic change are found mainly between major clades. Some patterns of ontogenetic variation and morphological disparity appear to reflect ontogenetic transitions in habitat use. Generally, juveniles are more similar to each other than adults, possibly because species that differ in ecology as adults are arboreal as juveniles. The magnitude of ontogenetic change follows evolutionary models where variation is constrained around an optimal value. Conversely, the direction of ontogenetic change may follow models with different adaptive optima per habitat use category or models where interspecific interactions influence its evolution. Finally, we found that the evolutionary rates of the ontogenetic allometric trajectories are phylogenetically variable.

Conclusions

The attributes of ontogenetic allometric trajectories and their evolutionary rates are phylogenetically heterogeneous in Paleoanguimorpha. Both allometric constraints and ecological factors have shaped ontogeny in the group. Our study highlights the evolutionary lability and adaptability of postnatal ontogeny, and teases apart how different evolutionary shifts in ontogeny contribute to the generation of morphological diversity at different evolutionary scales.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01970-6.

Rapid radiation and rampant reticulation: Phylogenomics of South American Liolaemus lizards

Understanding the factors that cause heterogeneity among gene trees can increase the accuracy of species trees. Discordant signals across the genome are commonly produced by incomplete lineage sorting (ILS) and introgression, which in turn can result in reticulate evolution. Species tree inference using the multispecies coalescent is designed to deal with ILS and is robust to low levels of introgression, but extensive introgression violates the fundamental assumption that relationships are strictly bifurcating. In this study, we explore the phylogenomics of the iconic Liolaemus subgenus of South American lizards, a group of over 100 species mostly distributed in and around the Andes mountains. Using mitochondrial DNA (mtDNA) and genome-wide restriction-site associated DNA sequencing (RADseq; nDNA hereafter), we inferred a time-calibrated mtDNA gene tree, nDNA species trees, and phylogenetic networks. We found high levels of discordance between mtDNA and nDNA, which we attribute in part to extensive ILS resulting from rapid diversification. These data also reveal extensive and deep introgression, which combined with rapid diversification, explain the high level of phylogenetic discordance. We discuss these findings in the context of Andean orogeny and glacial cycles that fragmented, expanded, and contracted species distributions. Finally, we use the new phylogeny to resolve long-standing taxonomic issues in one of the most studied lizard groups in the New World.

A return-on-investment approach for prioritization of rigorous taxonomic research needed to inform responses to the biodiversity crisis

Global biodiversity loss is a profound consequence of human activity. Disturbingly, biodiversity loss is greater than realized because of the unknown number of undocumented species. Conservation fundamentally relies on taxonomic recognition of species, but only a fraction of biodiversity is described. Here, we provide a new quantitative approach for prioritizing rigorous taxonomic research for conservation. We implement this approach in a highly diverse vertebrate group-Australian lizards and snakes. Of 870 species assessed, we identified 282 (32.4%) with taxonomic uncertainty, of which 17.6% likely comprise undescribed species of conservation concern. We identify 24 species in need of immediate taxonomic attention to facilitate conservation. Using a broadly applicable return-on-investment framework, we demonstrate the importance of prioritizing the fundamental work of identifying species before they are lost.

Phylogeography, historical demography and systematics of the world's smallest pythons (Pythonidae, Antaresia)

Advances from empirical studies in phylogeography, systematics and species delimitation highlight the importance of integrative approaches for quantifying taxonomic diversity. Genomic data have greatly improved our ability to discern both systematic diversity and evolutionary history. Here we combine analyses of mitochondrial DNA sequences, thousands of genome-wide SNPs and linear and geometric morphometrics on Antaresia, a clade of four currently recognised dwarf pythons from Australia and New Guinea (Antaresia childreni, A. stimsoni, A. maculosa and A. perthensis). Our integrative analyses of phylogenetics, population structure, species delimitation, historical demography and morphometrics revealed that the true evolutionary diversity is not well reflected in the current appraisal of the diversity of the group. We find that Antaresia childreni and A. stimsoni comprise a widespread network of populations connected by gene flow and without evidence of species-level divergence among them. However, A. maculosa shows considerable genetic structuring which leads us to recognise two subspecies in northeastern Australia and a new species in Torres Strait and New Guinea. These two contrasting cases of over and under estimation of diversity, respectively, illustrate the power of thorough integrative approaches into understanding evolution of biodiversity. Furthermore, our analyses of historical demographic patterns highlight the importance of the Kimberley, Pilbara and Cape York as origins of biodiversity in Australia.

Phylogenomics, biogeography and taxonomic revision of New Guinean pythons (Pythonidae, Leiopython) harvested for international trade

The large and enigmatic New Guinean pythons in the genus Leiopython are harvested from the wild to supply the international trade in pets. Six species are currently recognized (albertisii, biakensis, fredparkeri, huonensis, meridionalis, montanus) but the taxonomy of this group has been controversial. We combined analysis of 421 nuclear loci and complete mitochondrial genomes with morphological data to construct a detailed phylogeny of this group, understand their biogeographic patterns and establish the systematic diversity of this genus. Our molecular genetic data support two major clades, corresponding to L. albertisii and L. fredparkeri, but offer no support for the other four species. Our morphological data also only support two species. We therefore recognize L. albertisii and L. fredparkeri as valid species and place L. biakensis, L. meridionalis, L. huonensis and L. montanus into synonymy. We found that L. albertisii and L. fredparkeri are sympatric in western New Guinea; an atypical pattern compared to other Papuan species complexes in which the distributions of sister taxa are partitioned to the north and south of the island's central mountain range. For the purpose of conservation management, overestimation of species diversity within Leiopython has resulted in the unnecessary allocation of resources that could have been expended elsewhere. We strongly caution against revising the taxonomy of geographically widespread species groups when little or no molecular genetic data and only small morphological samples are available.

Integrating morphological and genetic data at different spatial scales in a cosmopolitan marine turtle species: challenges for management and conservation

Patterns of genetic structure in highly mobile marine vertebrates may be accompanied by phenotypic variation. Most studies in marine turtles focused on population genetic structure have been performed at rookeries. We studied whether genetic and morphological variation of the endangered green turtle (Chelonia mydas) is consistent geographically, focusing on foraging grounds. An association between population genetic structure and body shape variation at broad (inter-lineage) and fine (foraging grounds) scales was predicted and analysed using mitochondrial DNA and geometric morphometrics. Although genetic and phenotypic differentiation patterns were congruent between lineages, no fine-scale association was found, suggesting adaptive divergence. Connectivity among Pacific foraging grounds found here suggests that temperatures of ocean surface currents may influence the genetic structure of C. mydas on a broad scale. Our results suggest that vicariance, dispersal, life-history traits and ecological conditions operating in foraging grounds have shaped the intraspecific morphology and genetic diversity of this species. Considering a range of geographic and temporal scales is useful when management strategies are required for cosmopolitan species. Integrating morphological and genetic tools at different spatial scales, conservation management is proposed based on protection of neutral and adaptive diversity. This approach opens new questions and challenges, especially regarding conservation genetics in cosmopolitan species.

No link between nymph and adult coloration in shield bugs: weak selection by predators

Many organisms use different antipredator strategies throughout their life, but little is known about the reasons or implications of such changes. For years, it has been suggested that selection by predators should favour uniformity in local warning signals. If this is the case, we would expect high resemblance in colour across life stages in aposematic animals where young and adults share similar morphology and habitat. In this study, we used shield bugs (Hemiptera: Pentatomoidea) to test whether colour and colour diversity evolve similarly at different life stages. Since many of these bugs are considered to be aposematic, we also combined multi-species analyses with predation experiments on the cotton harlequin bug to test whether there is evidence of selection for uniformity in colour across life stages. Overall, we show that the diversity of colours used by both life stages is comparable, but adults are more cryptic than nymphs. We also demonstrate that nymphs and adults of the same species do not tend to look alike. Experiments on our model system suggest that predators can generalise among life stages that look different, and exhibit strong neophobia. Altogether, our results show no evidence of selection favouring colour similarity between adults and nymphs in this speciose clade.

Phylogenomics, Biogeography, and Morphometrics Reveal Rapid Phenotypic Evolution in Pythons After Crossing Wallace’s Line

Ecological opportunities can be provided to organisms that cross stringent biogeographic barriers towards environments with new ecological niches. Wallace’s and Lyddeker’s lines are arguably the most famous biogeographic barriers, separating the Asian and Australo-Papuan biotas. One of the most ecomorphologically diverse groups of reptiles, the pythons, is distributed across these lines, and are remarkably more diverse in phenotype and ecology east of Lydekker’s line in Australo-Papua. We used an anchored hybrid enrichment approach, with near complete taxon sampling, to extract mitochondrial genomes and 376 nuclear loci to resolve and date their phylogenetic history. Biogeographic reconstruction demonstrates that they originated in Asia around 38-45 Ma and then invaded Australo-Papua around 23 Ma. Australo-Papuan pythons display a sizeable expansion in morphological space, with shifts towards numerous new adaptive optima in head and body shape, coupled with the evolution of new micro-habitat preferences. We provide an updated taxonomy of pythons and our study also demonstrates how ecological opportunity following colonization of novel environments can promote morphological diversification in a formerly ecomorphologically conservative group. [Adaptive radiation; anchored hybrid enrichment; biogeography; morphometrics; snakes.]

Chicken adaptive response to low energy diet: main role of the hypothalamic lipid metabolism revealed by a phenotypic and multi-tissue transcriptomic approach

Background

Production conditions of layer chicken can vary in terms of temperature or diet energy content compared to the controlled environment where pure-bred selection is undertaken. The aim of this study was to better understand the long-term effects of a 15%-energy depleted diet on egg-production, energy homeostasis and metabolism via a multi-tissue transcriptomic analysis. Study was designed to compare effects of the nutritional intervention in two layer chicken lines divergently selected for residual feed intake.

Results

Chicken adapted to the diet in terms of production by significantly increasing their feed intake and decreasing their body weight and body fat composition, while their egg production was unchanged. No significant interaction was observed between diet and line for the production traits. The low energy diet had no effect on adipose tissue and liver transcriptomes. By contrast, the nutritional challenge affected the blood transcriptome and, more severely, the hypothalamus transcriptome which displayed 2700 differentially expressed genes. In this tissue, the low-energy diet lead to an over-expression of genes related to endocannabinoid signaling (CN1R, NAPE-PLD) and to the complement system, a part of the immune system, both known to regulate feed intake. Both mechanisms are associated to genes related polyunsaturated fatty acids synthesis (FADS1, ELOVL5 and FADS2), like the arachidonic acid, a precursor of anandamide, a key endocannabinoid, and of prostaglandins, that mediate the regulatory effects of the complement system. A possible regulatory role of NR1H3 (alias LXRα) has been associated to these transcriptional changes. The low-energy diet further affected brain plasticity-related genes involved in the cholesterol synthesis and in the synaptic activity, revealing a link between nutrition and brain plasticity. It upregulated genes related to protein synthesis, mitochondrial oxidative phosphorylation and fatty acid oxidation in the hypothalamus, suggesting reorganization in nutrient utilization and biological synthesis in this brain area.

Conclusions

We observed a complex transcriptome modulation in the hypothalamus of chicken in response to low-energy diet suggesting numerous changes in synaptic plasticity, endocannabinoid regulation, neurotransmission, lipid metabolism, mitochondrial activity and protein synthesis. This global transcriptomic reprogramming could explain the adaptive behavioral response (i.e. increase of feed intake) of the animals to the low-energy content of the diet.

Immunome differences between porcine ileal and jejunal Peyer's patches revealed by global transcriptome sequencing of gut-associated lymphoid tissues

The epithelium of the intestinal mucosa and the gut-associated lymphoid tissues (GALT) constitute an essential physical and immunological barrier against pathogens. In order to study the specificities of the GALT transcriptome in pigs, we compared the transcriptome profiles of jejunal and ileal Peyer’s patches (PPs), mesenteric lymph nodes (MLNs) and peripheral blood (PB) of four male piglets by RNA-Seq. We identified 1,103 differentially expressed (DE) genes between ileal PPs (IPPs) and jejunal PPs (JPPs), and six times more DE genes between PPs and MLNs. The master regulator genes FOXP3, GATA3, STAT4, TBX21 and RORC were less expressed in IPPs compared to JPPs, whereas the transcription factor BCL6 was found more expressed in IPPs. In comparison between IPPs and JPPs, our analyses revealed predominant differential expression related to the differentiation of T cells into Th1, Th2, Th17 and iTreg in JPPs. Our results were consistent with previous reports regarding a higher T/B cells ratio in JPPs compared to IPPs. We found antisense transcription for respectively 24%, 22% and 14% of the transcripts detected in MLNs, PPs and PB, and significant positive correlations between PB and GALT transcriptomes. Allele-specific expression analyses revealed both shared and tissue-specific cis-genetic control of gene expression.

Deciphering the genetic regulation of peripheral blood transcriptome in pigs through expression genome-wide association study and allele-specific expression analysis

BACKGROUND

Efforts to improve sustainability in livestock production systems have focused on two objectives: investigating the genetic control of immune function as it pertains to robustness and disease resistance, and finding predictive markers for use in breeding programs. In this context, the peripheral blood transcriptome represents an important source of biological information about an individual's health and immunological status, and has been proposed for use as an intermediate phenotype to measure immune capacity. The objective of this work was to study the genetic architecture of variation in gene expression in the blood of healthy young pigs using two approaches: an expression genome-wide association study (eGWAS) and allele-specific expression (ASE) analysis.

RESULTS

The blood transcriptomes of 60-day-old Large White pigs were analyzed by expression microarrays for eGWAS (242 animals) and by RNA-Seq for ASE analysis (38 animals). Using eGWAS, the expression levels of 1901 genes were found to be associated with expression quantitative trait loci (eQTLs). We recovered 2839 local and 1752 distant associations (Single Nucleotide Polymorphism or SNP located less or more than 1 Mb from expression probe, respectively). ASE analyses confirmed the extensive cis-regulation of gene transcription in blood, and revealed allelic imbalance in 2286 SNPs, which affected 763 genes. eQTLs and ASE-genes were widely distributed on all chromosomes. By analyzing mutually overlapping eGWAS results, we were able to describe putative regulatory networks, which were further refined using ASE data. At the functional level, genes with genetically controlled expression that were detected by eGWAS and/or ASE analyses were significantly enriched in biological processes related to RNA processing and immune function. Indeed, numerous distant and local regulatory relationships were detected within the major histocompatibility complex region on chromosome 7, revealing ASE for most class I and II genes.

CONCLUSIONS

This study represents, to the best of our knowledge, the first genome-wide map of the genetic control of gene expression in porcine peripheral blood. These results represent an interesting resource for the identification of genetic markers and blood biomarkers associated with variations in immunity traits in pigs, as well as any other complex traits for which blood is an appropriate surrogate tissue.

KIRREL is differentially expressed in adipose tissue from 'fertil+' and 'fertil-' cows: in vitro role in ovary?

We have previously shown that dairy cows carrying the 'fertil-' haplotype for one quantitative trait locus affecting female fertility located on the bovine chromosome three (QTL-F-Fert-BTA3) have a significantly lower conception rate and body weight after calving than cows carrying the 'fertil+' haplotype. Here, we compared by Tiling Array the expression of genes included in the QTL-F-Fert-BTA3 in 'fertil+' and 'fertil-' adipose tissue one week after calving when plasma non-esterified fatty acid concentrations were greater in 'fertil-' animals. We observed that thirty-one genes were overexpressed whereas twelve were under-expressed in 'fertil+' as compared to 'fertil-' cows (P < 0.05). By quantitative PCR and immunoblot we confirmed that adipose tissue KIRREL mRNA and protein were significantly greater expressed in 'fertil+' than in 'fertil-'. KIRREL mRNA is abundant in bovine kidney, adipose tissue, pituitary, and ovary and detectable in hypothalamus and mammary gland. Its expression (mRNA and protein) is greater in kidney of 'fertil+' than 'fertil-' cows (P < 0.05). KIRREL (mRNA and protein) is also present in the different ovarian cells with a greater expression in granulosa cells of 'fertil+' than 'fertil-' cows. In cultured granulosa cells, recombinant KIRREL halved steroid secretion in basal state (P < 0.05). It also decreased cell proliferation (P < 0.05) and in vitro oocyte maturation (P < 0.05). These results were associated to a rapid increase in MAPK1/3 and MAPK14 phosphorylation in granulosa cells and to a decrease in MAPK1/3 phosphorylation in oocyte. Thus, KIRREL could be a potential metabolic messenger linking body composition and fertility.

Liolaemus riodamas Esquerré, Núñez &amp; Scolaro, 2013 (Reptilia: Squamata: Liolaemidae): new record for the Maule Region in Chile

Liolaemus riodamas is a recently described lizard, only known from a very restricted river canyon in the Andes of the O’Higgins Region in Chile. The species is particular for being one of the few Liolaemus lacking the pheromone secreting glands known as precloacal pores. Here we document a second locality for this species located 30 km southwest of the type locality. This is also 500 m higher than the type locality, therefore a considerable altitudinal extension. This is a very positive discovery since the type locality of L. riodamas is intervened by a hydroelectric project, and gives hope to the conservation of this unique taxon.

Two new species of the Liolaemuselongatus-kriegi complex (Iguania, Liolaemidae) from Andean highlands of southern Chile

The elongatus-kriegi complex is one of the most diverse clades of the Liolaemus (sensu stricto) subgenus of lizards. There are currently 29 species recognized in this group distributed between Chile and Argentina. Based on molecular evidence, there seem to be five main clades nested within this complex: the elongatus, leopardinus, kriegi, petrophilus and punmahuida clades. Liolaemusbuergeri and Liolaemuskriegi, both of the kriegi clade, were believed to inhabit the surroundings of the Laja Lagoon, in the Biobío Region of Chile. Moreover, this Chilean population of Liolaemuskriegi was recently recognized as an undescribed taxon called “Liolaemus sp. A” based on molecular phylogenetics. In this work, we studied these two populations of the Laja Lagoon and provided the morphological diagnosis to describe them as two new species: Liolaemusscorialissp. n. and Liolaemuszabalaisp. n., previously considered Liolaemusbuergeri and “Liolaemuskriegi/Liolaemus sp. A” respectively. Additionally, we identified another population of Liolaemusscorialis in the vicinity of La Mula Lagoon in the Araucanía Region of Chile. Liolaemusscorialis differs from almost all of the species of the elongatus-kriegi complex by its considerably smaller size. Nevertheless, without molecular data we cannot assign it to any particular subclade. Liolaemuszabalai belongs to the kriegi clade based on published molecular phylogenies. Finally, we provide some natural history data on both species and we document for the first time the presence of Liolaemusneuquensis in Chile from a museum specimen from La Mula Lagoon.

Comparison of whole-genome (13X) and capture (87X) resequencing methods for SNP and genotype callings

The number of polymorphisms identified with next-generation sequencing approaches depends directly on the sequencing depth and therefore on the experimental cost. Although higher levels of depth ensure more sensitive and more specific SNP calls, economic constraints limit the increase of depth for whole-genome resequencing (WGS). For this reason, capture resequencing is used for studies focusing on only some specific regions of the genome. However, several biases in capture resequencing are known to have a negative impact on the sensitivity of SNP detection. Within this framework, the aim of this study was to compare the accuracy of WGS and capture resequencing on SNP detection and genotype calling, which differ in terms of both sequencing depth and biases. Indeed, we have evaluated the SNP calling and genotyping accuracy in a WGS dataset (13X) and in a capture resequencing dataset (87X) performed on 11 individuals. The percentage of SNPs not identified due to a sevenfold sequencing depth decrease was estimated at 7.8% using a down-sampling procedure on the capture sequencing dataset. A comparison of the 87X capture sequencing dataset with the WGS dataset revealed that capture-related biases were leading with the loss of 5.2% of SNPs detected with WGS. Nevertheless, when considering the SNPs detected by both approaches, capture sequencing appears to achieve far better SNP genotyping, with about 4.4% of the WGS genotypes that can be considered as erroneous and even 10% focusing on heterozygous genotypes. In conclusion, WGS and capture deep sequencing can be considered equivalent strategies for SNP detection, as the rate of SNPs not identified because of a low sequencing depth in the former is quite similar to SNPs missed because of method biases of the latter. On the other hand, capture deep sequencing clearly appears more adapted for studies requiring great accuracy in genotyping.