Do Shapes of Elevational Gradients of Species Richness Depend on the Vertical Range Studied? The Case of the Himalayas

We retrieved shapes of elevational species richness gradients (unimodal, decreasing, increasing) from 64 publications, studying Himalayan elevation patterns. We covered both plants and animals, and tested the hypothesis that unimodal gradients, explicable by the geometric mid-domain effect, prevail in the mountains, whereas decreasing or increasing gradients result from studying only short sections of entire altitudinal ranges. Multivariate canonical correspondence analysis was used to relate gradient shapes to their altitude ranges, geography positions, and taxa studied. Across taxa, most of the Himalayan altitudinal gradient display a unimodal shape, with a peak of diversity situated at ca 2500 m a.s.l. for plants, and 2200 m a.s.l. for animals. The gradient shapes were attributable to three intercorrelated predictors: vertical range, maximum elevation, and mean elevation of the gradients. Studies covering sufficiently broad altitudinal range returned unimodal gradients. Studies from the Earth’s highest mountain range reveal that surveys covering substantial parts of the elevational range of the mountains result in unimodal elevational gradients, whereas declining or increasing species richness gradients result from incomplete elevation range sampling.

Diverse Heat Tolerance of the Yeast Symbionts of Platycerus Stag Beetles in Japan

The genus Platycerus (Coleoptera: Lucanidae) is a small stag beetle group, which is adapted to cool-temperate deciduous broad-leaved forests in East Asia. Ten Platycerus species in Japan form a monophyletic clade endemic to Japan and inhabit species-specific climatic zones. They are reported to have co-evolutionary associations with their yeast symbionts of the genus Sheffersomyces based on host cytochrome oxidase subunit I (COI) and yeast intergenic spacer (IGS) phylogenies. Here we examined the heat tolerances of the yeast colonies isolated from the mycangia of 37 females belonging ten Japanese Platycerus species. The upper limits of growth and survival temperatures of each colony were decided by cultivating it at ten temperature levels between 17.5 and 40°C. Although both temperatures varied during 25.0–31.25°C, the maximum survival temperatures (MSTs) were a little higher than the maximum growth temperatures (MGTs) in 16 colonies. Pearson’s correlations between these temperatures and environmental factors (elevation and 19 bioclimatic variables from Worldclim database) of host beetle collection sites were calculated. These temperatures were significantly correlated with elevation negatively, the maximum temperature of the warmest month (Bio5) positively, and some precipitative variables, especially in the warm season (Bio12, 13, 16, 18) negatively. Sympatric Platycerus kawadai and Platycerus albisomni share the same lineage of yeast symbionts that exhibit the same heat tolerance, but the elevational lower range limit of P. kawadai is higher than that of P. albisomni. Based on the field survey in their sympatric site, the maximum temperature of host wood of P. kawadai larvae is higher about 2–3°C than that of P. albisomni larvae in the summer, which may restrict the elevational range of P. kawadai to higher area. In conclusion, it is suggested that the heat tolerance of yeast symbionts restricts the habitat range of their host Platycerus species or/and that the environmental condition that host Platycerus species prefers affect the heat tolerance of its yeast symbionts.

Seasonal distribution and habitat use preference of Barking deer (Muntiacus vaginalis) in Murree-Kotli Sattian-Kahuta National Park, Punjab Pakistan

ABSTRACT Microhabitat factors associated with the habitat of barking deer (Muntiacus vaginalis) were examined and compared between summer and winter seasons. Habitat characteristics and preferred habitat were measured by locating direct and indirect signs. To quantify the habitat utilization of barking deer, each selected study site was sampled for floral diversity from 2015 to 2017. Quadrats were deployed along transect lines to determine seasonal distribution. Barking deer were not evenly distributed across vegetation types in the study area; they occurred more often in the broad-leaved forest than in Chir pine forest, at an elevational range of 550-850 m, in thick vegetation on steep slopes. The most preferred habitat included trees and shrubs with 30% and 69% cover, respectively. Barking deer avoided thicker tree cover, possibly as it hinders movement and escape from predators. No significant difference (χ2 = 6.37, df = 3, p = 0.19) in seasonal vegetation cover was recorded.

A chromosome-scale genome assembly for the holly (Ilex polyneura) provides insights into genomic adaptations to elevation in Southwest China

Southwest China is a plant diversity hotspot. The near-cosmopolitan genus Ilex (c. 664 spp., Aquifoliaceae) reaches its maximum diversity in this region, with many narrow-range and a few widespread species. Divergent selection on widespread species leads to local adaptation, with consequences for both conservation and utilization, but is counteracted by geneflow. Many Ilex species are utilized as teas, medicines, ornamentals, honey plants, and timber, but variation below the species level is largely uninvestigated. We therefore studied the widespread Ilex polyneura, which occupies most of the elevational range available and is cultivated for its decorative leafless branches with persistent red fruits. We assembled a chromosome-scale genome using approximately 100x whole genome long-read and short-read sequencing combined with Hi-C sequencing. The genome is approximately 727.1 Mb, with a contig N50 size of 5 124 369 bp and a scaffold N50 size of 36 593 620 bp, for which the BUSCO score was 97.6%, and 98.9% of the assembly was anchored to 20 pseudochromosomes. Out of 32 838 genes predicted, 96.9% were assigned functions. Two whole genome duplication events were identified. Using this genome as a reference, we conducted a population genomics study of 112 individuals from 21 populations across the elevation range using restriction site-associated DNA sequencing (RADseq). Most populations clustered into four clades separated by distance and elevation. Selective sweep analyses identified 34 candidate genes potentially under selection at different elevations, with functions related to responses to abiotic and biotic stresses. This first high-quality genome in the Aquifoliales will facilitate the further domestication of the genus.

Rapoport’s Rule and the effect of the last glaciation upon elevational range size: An analysis using a dung beetle model (Coleoptera: Scarabaeidae: Onthophagus) in Mexican tropical mountains

Several studies have tested the Elevational Rapoport Rule (ERR) in arthropods, especially in the Neotropical mountains. Nonetheless, different approaches should be used for a more nuanced comprehension of ERR patterns and assemblage altitudinal distribution patterns, such as the biogeographical, ecological, and evolutionary contexts. This study aims to test the ERR for elevational gradients in Mexican mountains. For this study, dung beetle assemblages of the genus Onthophagus were used as a model organism, and their distribution was studied in several different mountain ranges of the Mexican tropics. Altitudinal distribution of Onthophagus species was analyzed, including ecological traits and biogeographical/phylogenetical contexts as covariables. The increase of altitude was positively correlated to the assemblage altitudinal range. Furthermore, altitudinal range, relative abundance, body size, and mountain’s topographic prominence were positively correlated to the mean altitudinal range of Onthophagus species. Nonetheless, different altitudinal relationships were observed, depending on the mountain. The results support the idea that species that inhabit higher altitudes appear to be more environmentally plastic and occur in wider altitudinal ranges than species from lower altitudes, thus supporting the ERR. The present findings stress that biogeographical, ecological, phylogenetical, and historical aspects, besides body size, are essential drivers of the altitudinal distribution of Onthophagus dung beetles.

Comparison of Ostracoda (Crustacea) species composition between lakes and streams at high elevations in Turkey

Ostracods are one of the most suitable organisms because of their ecological and morphological preferences for habitat types to study the diversity between streams and lakes. Here, 24 streams and 16 lakes were sampled in 2019 to evaluate the species compositions between streams and lakes. Environmental heterogeneity in the sense of the physico-chemical variables was found among the habitats. 24 of the 29 species with 923 individuals were found in the lakes when 13 species with 898 individuals were found in the streams. Eight species were common in both habitat types. Species contributing about 70.81% to the significant differences in the species composition between the habitats showed cosmopolitan characteristics. Species variation among the replicates of lakes were higher than streams. Results indicated that the species composition was probably associated to dissolved oxygen in lakes, while electrical conductivity in streams. Relationship amid elevational range and ostracod species richness was negative in lakes when there was no clear relationship in the streams. Overall, results suggested that determination of influential regional factors in different aquatic bodies can provide important contributions to the estimation of which habitat types can be chosen by species, and the using of species as bioindicators.

Area modulates the effect of elevation but not of land use or canopy on tropical plant species richness

One of the few general patterns in ecology is the increase of species richness with area. However, factors driving species-area relationship (SAR) are under debate, and the role of human-induced changes has been overlooked so far. Furthermore, SAR studies in tropical regions, in particular in multilayered rain forests are scarce. On the other side, studies of global change-induced impacts on biodiversity have become increasingly important, particular in the tropics, where these impacts are especially pronounced. Here, we investigated if area modulates the effect of land use, elevation and canopy on plant species richness. For the first time we studied SAR in multilayered tropical forests considering all functional groups. We selected 13 natural and disturbed habitats on Kilimanjaro in Tanzania, distributed over an elevational range of 3700 m. In each habitat type, we set up three to six modified Whittaker plots. We recorded all plant species in 64 plots and 640 subplots and described SAR using the power function. Area consistently modulated effects of elevation on plant species richness, partly effects of land use but not effects of plant canopy. Thus, area needs to be taken into account when studying elevational plant species richness patterns. In contrast to temperate regions open and forest habitats did not differ in SAR, probably due to a distinct vertical vegetation zonation in tropical forests. Therefore, it is important to consider all vegetation layers including epiphytes when studying SAR in highly structured tropical regions.

Selective sorting of ancestral introgression in maize and teosinte along an elevational cline

While often deleterious, hybridization can also be a key source of genetic variation and pre-adapted haplotypes, enabling rapid evolution and niche expansion. Here we evaluate these opposing selection forces on introgressed ancestry between maize (Zea mays ssp. mays) and its wild teosinte relative, mexicana (Zea mays ssp. mexicana). Introgression from ecologically diverse teosinte may have facilitated maize’s global range expansion, in particular to challenging high elevation regions (> 1500 m). We generated low-coverage genome sequencing data for 348 maize and mexicana individuals to evaluate patterns of introgression in 14 sympatric population pairs, spanning the elevational range of mexicana, a teosinte endemic to the mountains of Mexico. While recent hybrids are commonly observed in sympatric populations and mexicana demonstrates fine-scale local adaptation, we find that the majority of mexicana ancestry tracts introgressed into maize over 1000 generations ago. This mexicana ancestry seems to have maintained much of its diversity and likely came from a common ancestral source, rather than contemporary sympatric populations, resulting in relatively low FST between mexicana ancestry tracts sampled from geographically distant maize populations. Introgressed mexicana ancestry in maize is reduced in lower-recombination rate quintiles of the genome and around domestication genes, consistent with pervasive selection against introgression. However, we also find mexicana ancestry increases across the sampled elevational gradient and that high introgression peaks are most commonly shared among high-elevation maize populations, consistent with introgression from mexicana facilitating adaptation to the highland environment. In the other direction, we find patterns consistent with adaptive and clinal introgression of maize ancestry into sympatric mexicana at many loci across the genome, suggesting that maize also contributes to adaptation in mexicana, especially at the lower end of its elevational range. In sympatric maize, in addition to high introgression regions we find many genomic regions where selection for local adaptation maintains steep gradients in introgressed mexicana ancestry across elevation, including at least two inversions: the well-characterized 14 Mb Inv4m on chromosome 4 and a novel 3 Mb inversion Inv9f surrounding the macrohairless1 locus on chromosome 9. Most outlier loci with high mexicana introgression show no signals of sweeps or local sourcing from sympatric populations and so likely represent ancestral introgression sorted by selection, resulting in correlated but distinct outcomes of introgression in different contemporary maize populations.

Blood variation implicates respiratory limits on elevational ranges of Andean birds

The extent to which species ranges reflect intrinsic physiological tolerances is a major, unsolved question in evolutionary ecology. To date, consensus has been hindered by the limited tractability of experimental approaches across most of the tree of life. Here, we apply a macrophysiological approach to understand how hematological traits related to oxygen transport shape elevational ranges in a tropical biodiversity hotspot. Along Andean elevational gradients, we measured traits that affect blood oxygen-carrying capacity—total and cellular hemoglobin concentration and hematocrit—for 2,355 individuals of 136 bird species. We used these data to evaluate the influence of hematological traits on elevational ranges. First, we asked whether hematological plasticity is predictive of elevational range breadth. Second, we asked whether variance in hematological traits changed as a function of distance from the midpoint of the elevational range. We found that the correlation between hematological plasticity and elevational range breadth was slightly positive, consistent with a facilitative role for plasticity in elevational range expansion. We further found reduced local variation in hematological traits near elevational range limits and at high elevations, patterns consistent with intensified natural selection, reduced effective population size, or compensatory changes in other cardiohematological traits with increasing distance from species-specific optima for oxygen availability. Our findings suggest that constraints on hematological plasticity and local genetic adaptation to oxygen availability promote the evolution of the narrow elevational ranges that underpin tropical montane biodiversity.