Elevation,Temperature, and Aquatic Connectivity All Influence the Infection Dynamics of the Amphibian Chytrid Fungus in Adult Frogs

Infectious diseases can cause population declines and even extinctions. The amphibian chytrid fungus, Batrachochytrium dendrobatidis (Bd), has caused population declines and extinctions in amphibians on most continents. In the tropics, research on the dynamics of this disease has focused on amphibian populations in mountainous areas. In most of these areas, high and low elevation sites are connected by an assemblage of streams that may transport the infectious stage of the pathogen from high to low elevations, and, also, this pathogen, which grows well at cool temperatures, may persist better in cooler water flowing from high elevations. Thus, the dynamics of disease at low elevation sites without aquatic connections to higher elevation sites, i.e., non-contiguous low elevation sites, may differ from dynamics at contiguous low elevation sites. We sampled adult common mistfrogs (Litoria rheocola) at six sites of three types: two at high (> 400m) elevations, two at low elevations contiguous with high elevation streams, and two at low elevations non-contiguous with any high elevation site. Adults were swabbed for Bd diagnosis from June 2010 to June 2011 in each season, over a total of five sampling periods. The prevalence of Bd fluctuated seasonally and was highest in winter across all site types. Site type significantly affected seasonal patterns of prevalence of Bd. Prevalence remained well above zero throughout the year at the high elevation sites. Prevalence declined to lower levels in contiguous low sites, and reached near-zero at non-contiguous low sites. Patterns of air temperature fluctuation were very similar at both the low elevation site types, suggesting that differences in water connectivity to high sites may have affected the seasonal dynamics of Bd prevalence between contiguous and non-contiguous low elevation site types. Our results also suggest that reservoir hosts may be important in the persistence of disease at low elevations.     

Is dietary or microhabitat specialization associated with environmental heterogeneity in horned lizards (Phrynosoma)?

Niche breadth is predicted to correlate with environmental heterogeneity, such that generalists will evolve in heterogeneous environments and specialists will evolve in environments that vary less over space and time. We tested the hypothesis that lizards in a heterogeneous environment were generalists compared to lizards in a homogeneous environment. We compared niche breadths of greater short‐horned lizards by quantifying resource selection in terms of two different niche axes, diet (prey items and trophic level), and microhabitat (ground cover and shade cover) between two populations occurring at different elevations. We assessed the heterogeneity of dietary and microhabitat resources within each population's environment by quantifying the availability of prey items, ground cover, and shade cover in each environment. Overall, our results demonstrate that despite differences in resource heterogeneity between elevations, resource selection did not consistently differ between populations. Moreover, environmental heterogeneity was not associated with generalization of resource use. The low‐elevation site had a broader range of available prey items, yet lizards at the high‐elevation site demonstrated more generalization in diet. In contrast, the high‐elevation site had a broader range of available microhabitats, but the lizard populations at both sites were similarly generalized for shade cover selection and were similarly specialized for ground cover selection. Our results demonstrate that environmental heterogeneity of a particular resource does not necessarily predict the degree to which organisms specialize on that resource.     

Rapid Evolution of PARP Genes Suggests a Broad Role for ADP-Ribosylation in Host-Virus Conflicts

Post-translational protein modifications such as phosphorylation and ubiquitinylation are common molecular targets of conflict between viruses and their hosts. However, the role of other post-translational modifications, such as ADP-ribosylation, in host-virus interactions is less well characterized. ADP-ribosylation is carried out by proteins encoded by the PARP (also called ARTD) gene family. The majority of the 17 human PARP genes are poorly characterized. However, one PARP protein, PARP13/ZAP, has broad antiviral activity and has evolved under positive (diversifying) selection in primates. Such evolution is typical of domains that are locked in antagonistic ‘arms races’ with viral factors. To identify additional PARP genes that may be involved in host-virus interactions, we performed evolutionary analyses on all primate PARP genes to search for signatures of rapid evolution. Contrary to expectations that most PARP genes are involved in ‘housekeeping’ functions, we found that nearly one-third of PARP genes are evolving under strong recurrent positive selection. We identified a >300 amino acid disordered region of PARP4, a component of cytoplasmic vault structures, to be rapidly evolving in several mammalian lineages, suggesting this region serves as an important host-pathogen specificity interface. We also found positive selection of PARP9, 14 and 15, the only three human genes that contain both PARP domains and macrodomains. Macrodomains uniquely recognize, and in some cases can reverse, protein mono-ADP-ribosylation, and we observed strong signatures of recurrent positive selection throughout the macro-PARP macrodomains. Furthermore, PARP14 and PARP15 have undergone repeated rounds of gene birth and loss during vertebrate evolution, consistent with recurrent gene innovation. Together with previous studies that implicated several PARPs in immunity, as well as those that demonstrated a role for virally encoded macrodomains in host immune evasion, our evolutionary analyses suggest that addition, recognition and removal of ADP-ribosylation is a critical, underappreciated currency in host-virus conflicts.     

Degree of Functional Divergence in Duplicates Is Associated with Distinct Roles in Plant Evolution

Gene duplication is a major mechanism to create new genes. After gene duplication, some duplicated genes undergo functionalization, whereas others largely maintain redundant functions. Duplicated genes comprise various degrees of functional diversification in plants. However, the evolutionary fate of high and low diversified duplicates is unclear at genomic scale. To infer high and low diversified duplicates in Arabidopsis thaliana genome, we generated a prediction method for predicting whether a pair of duplicate genes was subjected to high or low diversification based on the phenotypes of knock-out mutants. Among 4,017 pairs of recently duplicated A. thaliana genes, 1,052 and 600 are high and low diversified duplicate pairs, respectively. The predictions were validated based on the phenotypes of generated knock-down transgenic plants. We determined that the high diversified duplicates resulting from tandem duplications tend to have lineage-specific functions, whereas the low diversified duplicates produced by whole-genome duplications are related to essential signaling pathways. To assess the evolutionary impact of high and low diversified duplicates in closely related species, we compared the retention rates and selection pressures on the orthologs of A. thaliana duplicates in two closely related species. Interestingly, high diversified duplicates resulting from tandem duplications tend to be retained in multiple lineages under positive selection. Low diversified duplicates by whole-genome duplications tend to be retained in multiple lineages under purifying selection. Taken together, the functional diversities determined by different duplication mechanisms had distinct effects on plant evolution.     

Contrasting patterns of lichen functional diversity and species richness across an elevation gradient

Major environmental gradients co‐vary with elevation and have been a longstanding natural tool allowing ecologists to study global diversity patterns at smaller scales, and to make predictions about the consequences of climate change. These analyses have traditionally studied taxonomic diversity, but new functional diversity approaches may provide a deeper understanding of the ecological mechanisms driving species assembly. We examined lichen taxonomic and functional diversity patterns on 195 plots (200 m²) together with forest structure along an elevational gradient of 1000 m in a temperate low mountain range (Bohemian Forest, Germany). Along this elevation gradient temperature decreased and precipitation increased, two macroclimatic variables critical for lichens. Elevation was more important than forest structure in driving taxonomic and functional diversity. While species richness increased with elevation, functional diversity decreased and revealed that community patterns shift with elevation from random to clustered, reflecting selection for key shared traits. Higher elevations favored species with a complex growth form (which takes advantage of high moisture) and asexual reproductive mode (facilitating establishment under low temperature conditions). Our analysis highlights the need to examine alternative forms of diversity and opens the avenue for community predictions about climate change. For a regional scenario with increasing temperature and decreasing availability of moisture, we expect a loss of specialized species with a complex growth form and those with vegetative organs at higher elevations in low mountain ranges in Europe.     

Elevational Variation in Diversity of Xanthomonas oryzae pv. oryzae in South‐West China

Virulence analysis and two polymerase chain reaction–based assays were used to evaluate the population structure of Xanthomonas oryzae pv. oryzae (Xoo) from different elevations ranging from 150 to 2600 m in south‐west China. Among the 218 isolates of Xoo, 18 pathotypes were identified using six near‐isogenic rice lines, each containing a single resistance gene. Among them, pathotype 9 predominated in low and mid‐elevations was virulent to all resistance genes, including Xa2, Xa3, xa5, xa13, Xa14 and Xa18. However, pathotype 2 was predominant at high elevation and was virulent to Xa18 only. The 18 pathotypes were grouped into four clusters. Isolates belonging to cluster 1 were mainly found at high and mid‐elevations, while those of cluster 4 were mainly found at low elevations. There were significant trends of virulence of isolates from low to high with the elevation from high to low. The ERIC and J3 primers were used to screen the genomes of 218 isolates, and 56 molecular haplotypes were found. Multiple correspondence analyses revealed that 56 haplotypes were divided into four putative genetic lineages. Lineage 2 was the most frequently detected from 150 to 2600 m; it was clearly shown that isolates from high elevation with 80% is much more than from low and mid‐elevation in the lineage. It is intriguing that genetic variation of Xoo is restricted by physical geographical barriers of elevations. This is the first report on the relationship of pathotypic and genotypic diversity of Xoo at different elevations.     

Interactive effects of disturbance and nitrogen availability on phosphorus dynamics of southern Appalachian forests

Understanding the main and interactive effects of chronically altered resource availability and disturbance on phosphorus (P) availability is increasingly important in light of the rapid pace at which human activities are altering these processes and potentially introducing P limitation. We measured P pools and fluxes in eighteen mixed forest stands at three elevations (low, mid, high) subjected to increasing atmospheric N deposition, where hemlock (Tsuga canadensis) was absent or declining due to infestation by the exotic hemlock woolly adelgid (Adelges tsugae). While total soil P was similar across the study area, phosphorus fractionation revealed distinct differences in the distribution of soil P fractions as elevation and N availability increased. Soils from high elevation plots where N availability was greatest had 139 % larger organic P pools and 55 % smaller residual and refractory P pools than soils from low elevation plots with less N availability, suggesting that increased N availability has driven the depletion of recalcitrant P pools by stimulating biotic demand and sequestration. These differences in P distribution among fractions influenced how tree mortality affected P dynamics. At high elevations, plots containing declining hemlocks had significantly greater foliar P concentrations and fluxes of P from the forest floor than reference plots at similar elevations, whereas at low and mid-elevations there were no consistent differences between plots. Across all elevation classes, hardwood foliar N:P ratios were lower in plots with declining hemlocks. Collectively, these results suggest that increased N availability enhances bioavailable P, which is sequestered in vegetation until disturbances liberate it.     

Pleistocene glacial cycles drove lineage diversification and fusion in the Yosemite toad (Anaxyrus canorus)

Species endemic to alpine environments can evolve via steep ecological selection gradients between lowland and upland environments. Additionally, many alpine environments have faced repeated glacial episodes over the past two million years, fracturing these endemics into isolated populations. In this “glacial pulse” model of alpine diversification, cycles of allopatry and ecologically divergent glacial refugia play a role in generating biodiversity, including novel admixed (“fused”) lineages. We tested for patterns of glacial pulse lineage diversification in the Yosemite toad (Anaxyrus [Bufo] canorus), an alpine endemic tied to glacially influenced meadow environments. Using double‐digest RADseq on populations densely sampled from a portion of the species range, we identified nine distinct lineages with divergence times ranging from 18 to 724 thousand years ago (ka), coinciding with multiple Sierra Nevada glacial events. Three lineages have admixed origins, and demographic models suggest these fused lineages have persisted throughout past glacial cycles. Directionality indices supported the hypothesis that some lineages recolonized Yosemite from east of the ice sheet, whereas other lineages remained in western refugia. Finally, refugial niche reconstructions suggest that low‐ and high‐elevation lineages have convergently adapted to similar climatic niches. Our results suggest glacial cycles and refugia may be important crucibles of adaptive diversity across deep evolutionary time.     

Comparative genomics reveals high rates of horizontal transfer and strong purifying selection on rhizobial symbiosis genes

Horizontal transfer (HT) alters the repertoire of symbiosis genes in rhizobial genomes and may play an important role in the on-going evolution of the rhizobia–legume symbiosis. To gain insight into the extent of HT of symbiosis genes with different functional roles (nodulation, N-fixation, host benefit and rhizobial fitness), we conducted comparative genomic and selection analyses of the full-genome sequences from 27 rhizobial genomes. We find that symbiosis genes experience high rates of HT among rhizobial lineages but also bear signatures of purifying selection (low Ka : Ks). HT and purifying selection appear to be particularly strong in genes involved in initiating the symbiosis (e.g. nodulation) and in genome-wide association candidates for mediating benefits provided to the host. These patterns are consistent with rhizobia adapting to the host environment through the loss and gain of symbiosis genes, but not with host-imposed positive selection driving divergence of symbiosis genes through recurring bouts of positive selection.     

Habitat availability for bottomland hardwood forest birds: the importance of considering elevation

ABSTRACT Estimating habitat availability at large spatial scales is critical for identifying conservation and management priorities for birds. Given the effects of patch size on habitat selection and productivity of many bird species, large‐scale approaches often focus on identifying and enumerating patches large enough to support sustainable populations. Declines in the availability of bottomland hardwood forests have made this approach imperative for species that depend on these forests. However, most remaining bottomland forests are relatively low‐elevation sites that were difficult to convert to agriculture. For species that require densely vegetated understories or well‐developed litter layers, such as Swainson's Warblers (Limnothlypis swainsonii), these relatively flood‐prone forests may not provide suitable habitat and, consequently, current prioritization methods may overestimate habitat availability. We examined the effect of elevation on estimates of habitat availability for Swainson's Warblers in a 65,000‐ha bottomland hardwood forest (White River National Wildlife Refuge) in eastern Arkansas. We detected Swainson's Warblers at 88 of 2222 sampled points, and Swainson's Warblers only used relatively high elevations in the refuge. Based on our estimates, about 25% of the refuge is at a suitable elevation for Swainson's Warblers. Without considering elevation, this refuge would support an estimated 7000 to 7600 pairs, but, when elevation is considered, the estimated number of pairs drops to 1500–2000 pairs. In reality, because of a lack of suitable seral stages at high elevations in the refuge, this area likely supports only 75–100 pairs. For Swainson's Warblers and other understory‐dependent bottomland species, elevation should be incorporated into conservation planning to obtain accurate estimates of habitat availability. In addition, management should be focused on these high‐elevation areas to maximize habitat availability for these species of concern.     

The Genomic Architecture of Adaptation to Larval Malnutrition Points to a Trade-off with Adult Starvation Resistance in Drosophila

Periods of nutrient shortage impose strong selection on animal populations. Experimental studies of genetic adaptation to nutrient shortage largely focus on resistance to acute starvation at adult stage; it is not clear how conclusions drawn from these studies extrapolate to other forms of nutritional stress. We studied the genomic signature of adaptation to chronic juvenile malnutrition in six populations of Drosophila melanogaster evolved for 150 generations on an extremely nutrient-poor larval diet. Comparison with control populations evolved on standard food revealed repeatable genomic differentiation between the two set of population, involving >3,000 candidate SNPs forming >100 independently evolving clusters. The candidate genomic regions were enriched in genes implicated in hormone, carbohydrate, and lipid metabolism, including some with known effects on fitness-related life-history traits. Rather than being close to fixation, a substantial fraction of candidate SNPs segregated at intermediate allele frequencies in all malnutrition-adapted populations. This, together with patterns of among-population variation in allele frequencies and estimates of Tajima’s D, suggests that the poor diet results in balancing selection on some genomic regions. Our candidate genes for tolerance to larval malnutrition showed a high overlap with genes previously implicated in acute starvation resistance. However, adaptation to larval malnutrition in our study was associated with reduced tolerance to acute adult starvation. Thus, rather than reflecting synergy, the shared genomic architecture appears to mediate an evolutionary trade-off between tolerances to these two forms of nutritional stress.     

Genomic Patterns of Positive Selection at the Origin of Rust Fungi

Understanding the origin and evolution of pathogenicity and biotrophic life-style of rust fungi has remained a conundrum for decades. Research on the molecular mechanisms responsible for rust fungi evolution has been hampered by their biotrophic life-style until the sequencing of some rust fungi genomes. With the availability of multiple whole genomes and EST data for this group, it is now possible to employ genome-wide surveys and investigate how natural selection shaped their evolution. In this work, we employed a phylogenomics approach to search for positive selection and genes undergoing accelerated evolution at the origin of rust fungi on an assembly of single copy genes conserved across a broad range of basidiomycetes. Up to 985 genes were screened for positive selection on the phylogenetic branch leading to rusts, revealing a pervasive signal of positive selection throughout the data set with the proportion of positively selected genes ranging between 19.6–33.3%. Additionally, 30 genes were found to be under accelerated evolution at the origin of rust fungi, probably due to a mixture of positive selection and relaxation of purifying selection. Functional annotation of the positively selected genes revealed an enrichment in genes involved in the biosynthesis of secondary metabolites and several metabolism and transporter classes.     

Gene Expression Ratios Lead to Accurate and Translatable Predictors of DR5 Agonism across Multiple Tumor Lineages

Death Receptor 5 (DR5) agonists demonstrate anti-tumor activity in preclinical models but have yet to demonstrate robust clinical responses. A key limitation may be the lack of patient selection strategies to identify those most likely to respond to treatment. To overcome this limitation, we screened a DR5 agonist Nanobody across >600 cell lines representing 21 tumor lineages and assessed molecular features associated with response. High expression of DR5 and Casp8 were significantly associated with sensitivity, but their expression thresholds were difficult to translate due to low dynamic ranges. To address the translational challenge of establishing thresholds of gene expression, we developed a classifier based on ratios of genes that predicted response across lineages. The ratio classifier outperformed the DR5+Casp8 classifier, as well as standard approaches for feature selection and classification using genes, instead of ratios. This classifier was independently validated using 11 primary patient-derived pancreatic xenograft models showing perfect predictions as well as a striking linearity between prediction probability and anti-tumor response. A network analysis of the genes in the ratio classifier captured important biological relationships mediating drug response, specifically identifying key positive and negative regulators of DR5 mediated apoptosis, including DR5, CASP8, BID, cFLIP, XIAP and PEA15. Importantly, the ratio classifier shows translatability across gene expression platforms (from Affymetrix microarrays to RNA-seq) and across model systems (in vitro to in vivo). Our approach of using gene expression ratios presents a robust and novel method for constructing translatable biomarkers of compound response, which can also probe the underlying biology of treatment response.     

Ranging of <Emphasis Type="Italic">Rhinopithecus bieti</Emphasis> in the Samage Forest,China. II. Use of Land Cover Types and Altitudes

We investigated composition and structure of a temperate montane forest called Samage at Baimaxueshan National Nature Reserve in Yunnan, one of the last refuges for the highly endangered black-and-white snub-nosed monkeys (Rhinopithecus bieti). There is a patchwork of vegetation types at Samage, and we distinguished 6 major land cover types within the home range of the focal group. We tracked the semihabituated Gehuaqing band for a full annual cycle to study their habitat utilization and altitudinal ranging. We analyzed the group’s selective use of particular habitat types via selection ratios. We calculated habitat availability from a GIS database. We found that they used mixed deciduous broadleaf/conifer forest disproportionately to its availability in all months. Subjects completely avoided meadows. Pine and evergreen broadleaf forests acted as corridors between patches of mixed forest and monkeys visited them occasionally, but at low frequencies and mostly in transit. The focal band stayed at elevations ranging from ca. 2600 m to 4000 m, and the mean elevation used is 3200 m. We found evidence for seasonal variation in use of elevational zones. The band stayed at significantly higher elevations in summer than in spring. The descent in spring was likely related to a flush of immature leaves at low-lying elevations. Availability of preferred fruits also had a highly positive influence on altitudinal ranging, i.e., during months with high fruit availability (late summer, early fall), the band stayed at medium elevations where preferred fruits were most abundant. Higher concentrations of lichens and the snub-nosed monkeys’ search for not yet depleted fruits probably caused them to remain at mid-elevations in winter. There is no significant correlation between climate parameters and elevation used. One of the main inferences of this investigation is that, contrary to previous accounts, Rhinopithecus bieti is not universally associated with high-elevation dark fir forest, but at Samage exhibits an overwhelming preference for mixed forest. Moreover, our analyses support the hypothesis that elevational migration, in this temperate-subtropical forest, is influenced by the temporal fruiting of major food trees and that climate has only a negligible effect on elevation use. Dayong Li and Cyril C. Grueter contributed equally to the paper.     

Adaptive plasticity and epigenetic variation in response to warming in an Alpine plant

Environmentally induced phenotypic plasticity may be a critical component of response to changing environments. We examined local differentiation and adaptive phenotypic plasticity in response to elevated temperature in half‐sib lines collected across an elevation gradient for the alpine herb, Wahlenbergia ceracea. Using Amplified Fragment Length Polymorphism (AFLP), we found low but significant genetic differentiation between low‐ and high‐elevation seedlings, and seedlings originating from low elevations grew faster and showed stronger temperature responses (more plasticity) than those from medium and high elevations. Furthermore, plasticity was more often adaptive for plants of low‐elevation origin and maladaptive for plants of high elevation. With methylation sensitive‐AFLP (MS‐AFLP), we revealed an increase in epigenetic variation in response to temperature in low‐elevation seedlings. Although we did not find significant direct correlations between MS‐AFLP loci and phenotypes, our results demonstrate that adaptive plasticity in temperature response to warming varies over fine spatial scales and suggest the involvement of epigenetic mechanisms in this response.     

Animal Mitochondria,Positive Selection and Cyto-Nuclear Coevolution: Insights from Pulmonates

Pulmonate snails have remarkably high levels of mtDNA polymorphism within species and divergence between species, making them an interesting group for the study of mutation and selection on mitochondrial genomes. The availability of sequence data from most major lineages – collected largely for studies of phylogeography - provides an opportunity to perform several tests of selection that may provide general insights into the evolutionary forces that have produced this unusual pattern. Several protein coding mtDNA datasets of pulmonates were analyzed towards this direction. Two different methods for the detection of positive selection were used, one based on phylogeny, and the other on the McDonald-Kreitman test. The cyto-nuclear coevolution hypothesis, often implicated to account for the high levels of mtDNA divergence of some organisms, was also addressed by assessing the divergence pattern exhibited by a nuclear gene. The McDonald-Kreitman test indicated multiple signs of positive selection in the mtDNA genes, but was significantly biased when sequence divergence was high. The phylogenetic method identified five mtDNA datasets as affected by positive selection. In the nuclear gene, the McDonald-Kreitman test provided no significant results, whereas the phylogenetic method identified positive selection as likely present. Overall, our findings indicate that: 1) slim support for the cyto-nuclear coevolution hypothesis is present, 2) the elevated rates of mtDNA polymorphims and divergence in pulmonates do not appear to be due to pervasive positive selection, 3) more stringent tests show that spurious positive selection is uncovered when distant taxa are compared and 4) there are significant examples of positive selection acting in some cases, so it appears that mtDNA evolution in pulmonates can escape from strict deleterious evolution suggested by the Muller’s ratchet effect.