Landscape genetics of the Alpine newt (Mesotriton alpestris) inferred from a strip-based approach

Habitat destruction and fragmentation are known to strongly affect dispersal by altering the quality of the environment between populations. As a consequence, lower landscape connectivity is expected to enhance extinction risks through a decrease in gene flow and the resulting negative effects of genetic drift, accumulation of deleterious mutations and inbreeding depression. Such phenomena are particularly harmful for amphibian species, characterized by disjunct breeding habitats. The dispersal behaviour of amphibians being poorly understood, it is crucial to develop new tools, allowing us to determine the influence of landscape connectivity on the persistence of populations. In this study, we developed a new landscape genetics approach that aims at identifying land-uses affecting genetic differentiation, without a priori assumptions about associated ecological costs. We surveyed genetic variation at seven microsatellite loci for 19 Alpine newt (Mesotriton alpestris) populations in western Switzerland. Using strips of varying widths that define a dispersal corridor between pairs of populations, we were able to identify land-uses that act as dispersal barriers (i.e. urban areas) and corridors (i.e. forests). Our results suggest that habitat destruction and landscape fragmentation might in the near future affect common species such as M. alpestris. In addition, by identifying relevant landscape variables influencing population structure without unrealistic assumptions about dispersal, our method offers a simple and flexible tool of investigation as an alternative to least-cost models and other approaches.

Combining niche modelling and landscape genetics to study local adaptation: A novel approach illustrated using alpine plants

Understanding the factors that shape adaptive genetic variation across species niches has become of paramount importance in evolutionary ecology, especially to understand how adaptation to changing climate affects the geographic range of species. The distribution of adaptive alleles in the ecological niche is determined by the emergence of novel mutations, their fitness consequences and gene flow that connects populations across species niches. Striking demographical differences and source sink dynamics of populations between the centre and the margin of the niche can play a major role in the emergence and spread of adaptive alleles. Although some theoretical predictions have long been proposed, the origin and distribution of adaptive alleles within species niches remain untested. In this paper, we propose and discuss a novel empirical approach that combines landscape genetics with species niche modelling, to test whether alleles that confer local adaptation are more likely to occur in either marginal or central populations of species niches. We illustrate this new approach by using a published data set of 21 alpine plant species genotyped with a total of 2483 amplified fragment length polymorphisms (AFLP), distributed over more than 1733 sampling sites across the Alps. Based on the assumption that alleles that were statistically associated with environmental variables were adaptive, we found that adaptive alleles in the margin of a species niche were also present in the niche centre, which suggests that adaptation originates in the niche centre. These findings corroborate models of species range evolution, in which the centre of the niche contributes to the emergence of novel adaptive alleles, which diffuse towards niche margins and facilitate niche and range expansion through subsequent local adaptation. Although these results need to be confirmed via fitness measurements in natural populations and functionally characterised genetic sequences, this study provides a first step towards understanding how adaptive genetic variation emerges and shapes species niches and geographic ranges along environmental gradients.

The influence of landscape characteristics and home-range size on the quantification of landscape-genetics relationships

A common approach used to estimate landscape resistance involves comparing correlations of ecological and genetic distances calculated among individuals of a species. However, the location of sampled individuals may contain some degree of spatial uncertainty due to the natural variation of animals moving through their home range or measurement error in plant or animal locations. In this study, we evaluate the ways that spatial uncertainty, landscape characteristics, and genetic stochasticity interact to influence the strength and variability of conclusions about landscape-genetics relationships. We used a neutral landscape model to generate 45 landscapes composed of habitat and non-habitat, varying in percent habitat, aggregation, and structural connectivity (patch cohesion). We created true and alternate locations for 500 individuals, calculated ecological distances (least-cost paths), and simulated genetic distances among individuals. We compared correlations between ecological distances for true and alternate locations. We then simulated genotypes at 15 neutral loci and investigated whether the same influences could be detected in simple Mantel tests and while controlling for the effects of isolation-by distance using the partial Mantel test. Spatial uncertainty interacted with the percentage of habitat in the landscape, but led to only small reductions in correlations. Furthermore, the strongest correlations occurred with low percent habitat, high aggregation, and low to intermediate levels of cohesion. Overall genetic stochasticity was relatively low and was influenced by landscape characteristics.

Landscape Genetics of Raccoons (Procyon lotor) Associated with Ridges and Valleys of Pennsylvania: Implications for Oral Rabies Vaccination Programs

Raccoons are the reservoir for the raccoon rabies virus variant in the United States. To combat this threat, oral rabies vaccination (ORV) programs are conducted in many eastern states. To aid in these efforts, the genetic structure of raccoons (Procyon lotor) was assessed in southwestern Pennsylvania to determine if select geographic features (i.e., ridges and valleys) serve as corridors or hindrances to raccoon gene flow (e.g., movement) and, therefore, rabies virus trafficking in this physiographic region. Raccoon DNA samples (n = 185) were collected from one ridge site and two adjacent valleys in southwestern Pennsylvania (Westmoreland, Cambria, Fayette, and Somerset counties). Raccoon genetic structure within and among these study sites was characterized at nine microsatellite loci. Results indicated that there was little population subdivision among any sites sampled. Furthermore, analyses using a model-based clustering approach indicated one essentially panmictic population was present among all the raccoons sampled over a reasonably broad geographic area (e.g., sites up to 36 km apart). However, a signature of isolation by distance was detected, suggesting that widths of ORV zones are critical for success. Combined, these data indicate that geographic features within this landscape influence raccoon gene flow only to a limited extent, suggesting that ridges of this physiographic system will not provide substantial long-term natural barriers to rabies virus trafficking. These results may be of value for future ORV efforts in Pennsylvania and other eastern states with similar landscapes.

Landscape genetics of Phaedranassa herb. (Amaryllidaceae) in Ecuador

Speciation can be understood as a continuum occurring at different levels, from population to species. The recent molecular revolution in population genetics has opened a pathway towards understanding species evolution. At the same time, speciation patterns can be better explained by incorporating a geographic context, through the use of geographic information systems (GIS). Phaedranassa (Amaryllidaceae) is a genus restricted to one of the world’s most biodiverse hotspots, the Northern Andes. I studied seven Phaedranassa species from Ecuador. Six of these species are endemic to the country. The topographic complexity of the Andes, which creates local microhabitats ranging from moist slopes to dry valleys, might explain the patterns of Phaedranassa species differentiation. With a Bayesian individual assignment approach, I assessed the genetic structure of the genus throughout Ecuador using twelve microsatellite loci. I also used bioclimatic variables and species geographic coordinates under a Maximum Entropy algorithm to generate distribution models of the species. My results show that Phaedranassa species are genetically well-differentiated. Furthermore, with the exception of two species, all Phaedranassa showed non-overlapping distributions. Phaedranassa viridiflora and P. glauciflora were the only species in which the model predicted a broad species distribution, but genetic evidence indicates that these findings are likely an artifact of species delimitation issues. Both genetic differentiation and nonoverlapping geographic distribution suggest that allopatric divergence could be the general model of genetic differentiation. Evidence of sympatric speciation was found in two geographically and genetically distinct groups of P. viridiflora. Additionally, I report the first register of natural hybridization for the genus. The findings of this research show that the genetic differentiation of species in an intricate landscape as the Andes does not necessarily show a unique trend. Although allopatric speciation is the most common form of speciation, I found evidence of sympatric speciation and hybridization. These results show that the processes of speciation in the Andes have followed several pathways. The mixture of these processes contributes to the high biodiversity of the region.

Identifying genetic signatures of selection in a non-model species, alpine gentian (Gentiana nivalis L.), using a landscape genetic approach

It is generally accepted that most plant populations are locally adapted. Yet, understanding how environmental forces give rise to adaptive genetic variation is a challenge in conservation genetics and crucial to the preservation of species under rapidly changing climatic conditions. Environmental variation, phylogeographic history, and population demographic processes all contribute to spatially structured genetic variation, however few current models attempt to separate these confounding effects. To illustrate the benefits of using a spatially-explicit model for identifying potentially adaptive loci, we compared outlier locus detection methods with a recently-developed landscape genetic approach. We analyzed 157 loci from samples of the alpine herb Gentiana nivalis collected across the European Alps. Principle coordinates of neighbor matrices (PCNM), eigenvectors that quantify multi-scale spatial variation present in a data set, were incorporated into a landscape genetic approach relating AFLP frequencies with 23 environmental variables. Four major findings emerged. 1) Fifteen loci were significantly correlated with at least one predictor variable (R (adj) (2) > 0.5). 2) Models including PCNM variables identified eight more potentially adaptive loci than models run without spatial variables. 3) When compared to outlier detection methods, the landscape genetic approach detected four of the same loci plus 11 additional loci. 4) Temperature, precipitation, and solar radiation were the three major environmental factors driving potentially adaptive genetic variation in G. nivalis. Techniques presented in this paper offer an efficient method for identifying potentially adaptive genetic variation and associated environmental forces of selection, providing an important step forward for the conservation of non-model species under global change.

Estudio de los modos de vida e historia biológica de las poblaciones humanas de la región central de argentina desde el registro bioarqueológico. Lifestyles and biological history of human populations of the central region of argentina from bioarchaeological record.

Se pretende aportar al estudio de la estructura, historia biológica y estilos de vida de las poblaciones que habitaron la región central de Argentina durante el Holoceno, desde una perspectiva que combina los aportes teóricos y metodológicos de la Genética del paisaje y la Bioarqueología. Interesa a) identificar barreras de diferenciación morfológica entre poblaciones, b) poner a prueba modelos poblacionales para explicar la variación observada e identificar las variables que contribuyan a dicha diferenciación, c) evaluar la congruencia de los resultados obtenidos, d) reconstruir los patrones de movilidad residencial de las poblaciones, e) estudiar sus patrones dietarios considerando diferencias temporales y espaciales, f) identificar indicadores de diversos tipos de estrés (nutricional, funcional), así como traumas, g) estudiar las historias tafonómicas del registro bioarqueológico regional, y h) proponer un modelo para explicar el poblamiento y la evolución local de las poblaciones que habitaron esta región, a partir de la información arqueológica y bioantropológica. Para el análisis de los patrones espaciales de variación biológica se trabajará a partir del registro de rasgos epigenéticos craneales, medidas lineales y datos obtenidos a partir de morfometría geométrica sobre fotografías en 2D sobre muestras arqueológicas procedentes de esta región y de otras regiones geográficas de la Argentina. Para el análisis de la estructura de la población se trabajará a partir del cálculo de la matriz R para datos morfológicos y sus estimaciones derivadas (distancia D², Fst, coordenadas principales) y la aplicación del modelo de Harpending y Ward. Desde la genética del paisaje, se realizarán análisis de autocorrelación espacial, barreras genéticas y análisis geoestadísticos (kriging). Para el estudio de los modos de vida a partir del registro bioarqueológico se relevarán patologías dento-alveolares y alteraciones vinculadas con la salud bucal tales como desgaste dental –a nivel micro y macroscópico- caries, abscesos, pérdidas dentales antemortem, cálculos, hipoplasias, marcadores esqueletales de salud y lesiones traumáticas. Se analizarán isótopos estables (δ13C, δ15N, 86Sr y 87Sr) en restos óseos humanos de diversos sitios arqueológicos con el objetivo de reconstruir patrones dietarios y analizar la movilidad residencial y migración de las poblaciones. Paralelamente, se establecerán procedimientos de control tafonómico de los restos óseos, y se harán análisis específicos para estudiar las historias tafonómicas y evaluar el grado de integridad de los contextos de depositación y de las colecciones en general. Estimamos que el análisis de los patrones espaciales y temporales de variabilidad morfológica craneofacial, así como el estudio de las dietas a partir de información isotópica y bioarqueológica, de las migraciones y la movilidad residencial de las poblaciones a partir de isótopos de estroncio, la reconstrucción de comportamientos y actividades cotidianas a partir de marcadores de estrés músculo-esqueletal, en un marco cronológico y espacial constituye un aporte novedoso y eficaz que permitirá incrementar de manera substancial la información sobre la evolución de las poblaciones originarias del centro del territorio argentino. The aim of this project is to study the structure, biological history and lifestyles of the people that inhabitated the central region of Argentina during the Holocene, from a perspective that combines theoretical and methodological contributions of Landscape Genetics and Bioarchaeology. To analyze the spatial patterns of biological variation we consider epigenetic cranial traits, linear measurements and data obtained from geometric morphometric on 2D photographs. Morphological variation will be focused on landscape genetics (autocorrelation, genetic barriers and geostatistical analysis –kriging-) and population structure (matrix R, D², Fst, principal coordinates, Harpendig and Ward model). For the study of lifestyles from bioarchaeological record we consider alveolar pathologies and disorders related to oral health such as tooth wear, micro and macroscopic level, caries, abscesses, antemortem tooth loss, hypoplasia, markers skeletal health and traumatic injuries, as well as taphonomic processes. Stable isotopes will be analyzed (δ13C, d15N, 86Sr and 87Sr) in human skeletal remains from various archaeological sites in order to reconstruct and analyze dietary patterns of residential mobility and migration of populations. It will be established procedures of taphonomic control on skeletal remains, analysis to study taphonomic histories and assess the degree of completeness of depositional context and collection, in general terms. We consider that analysis of spatial and temporal patterns of variability in craniofacial morphology and the study of health and diets from isotopic and bioarchaeological data, migration and residential mobility patterns from strontium isotopes, as well as activity patterns from stress markers is a novel and effective contribution that will substantially increase the information about the local evolution of populations that inhabitated the center of Argentina.

Multiple refugia and barriers explain the phylogeography of the Valais shrew, Sorex antinorii (Mammalia: Soricomorpha), Sorex antinorii (Mammalia: Soricomorpha)

The aim of the present study was to investigate the genetic structure of the Valais shrew (Sorex antinorii) by a combined phylogeographical and landscape genetic approach, and thereby to infer the locations of glacial refugia and establish the influence of geographical barriers. We sequenced part of the mitochondrial cytochrome b (cyt b) gene of 179 individuals of S. antinorii sampled across the entire species' range. Six specimens attributed to S. arunchi were included in the analysis. The phylogeographical pattern was assessed by Bayesian molecular phylogenetic reconstruction, population genetic analyses, and a species distribution modelling (SDM)-based hindcasting approach. We also used landscape genetics (including isolation-by-resistance) to infer the determinants of current intra-specific genetic structure. The phylogeographical analysis revealed shallow divergence among haplotypes and no clear substructure within S. antinorii. The starlike structure of the median-joining network is consistent with population expansion from a single refugium, probably located in the Apennines. Long branches observed on the same network also suggest that another refugium may have existed in the north-eastern part of Italy. This result is consistent with SDM, which also suggests several habitable areas for S. antinorii in the Italian peninsula during the LGM. Therefore S. antinorii appears to have occupied disconnected glacial refugia in the Italian peninsula, supporting previous data for other species showing multiple refugia within southern refugial areas. By coupling genetic analyses and SDM, we were able to infer how past climatic suitability contributed to genetic divergence of populations. The genetic differentiation shown in the present study does not support the specific status of S. arunchi.

Genetic and phenotypic population divergence on a microgeographic scale in brown trout.

Salmonid populations of many rivers are rapidly declining. One possible explanation is that habitat fragmentation increases genetic drift and reduces the populations' potential to adapt to changing environmental conditions. We measured the genetic and eco-morphological diversity of brown trout (Salmo trutta) in a Swiss stream system, using multivariate statistics and Bayesian clustering. We found large genetic and phenotypic variation within only 40 km of stream length. Eighty-eight percent of all pairwise F(ST) comparisons and 50% of the population comparisons in body shape were significant. High success rates of population assignment tests confirmed the distinctiveness of populations in both genotype and phenotype. Spatial analysis revealed that divergence increased with waterway distance, the number of weirs, and stretches of poor habitat between sampling locations, but effects of isolation-by-distance and habitat fragmentation could not be fully disentangled. Stocking intensity varied between streams but did not appear to erode genetic diversity within populations. A lack of association between phenotypic and genetic divergence points to a role of local adaptation or phenotypically plastic responses to habitat heterogeneity. Indeed, body shape could be largely explained by topographic stream slope, and variation in overall phenotype matched the flow regimes of the respective habitats.

A genetic isolation gradient of populations of the Balearic green toad (Bufo balearicus) follows rising eastward fragmentation of the rural landscapes on the island of Menorca.

We present the first approach to the genetic diversity and structure of the Balearic toad (Bufo balearicus Boettger, 1880) for the island of Menorca. Forty-one individ- uals from 21 localities were analyzed for ten microsatellite loci. We used geo-refer- enced individual multilocus genotypes and a model-based clustering method for the inference of the number of populations and of the spatial location of genetic dis- continuities between those populations.¦Only six of the microsatellites analyzed were polymorphic. We revealed a northwest- ern area inhabited by a single population with several well-connected localities and another set of populations in the southeast that includes a few unconnected small units with genetically significant differences among them as well as with the individ- uals from the northwest of the island. The observed fragmentation may be explained by shifts from agricultural to tourism practices that have been taking place on the island of Menorca since the 1960s. The abandonment of rural activities in favor of urbanization and concomitant service areas has mostly affected the southeast of the island and is currently threatening the overall geographic connectivity between the different farming areas of the island that are inhabited by the Balearic toad.

Phylogeography and dispersal in the velvet gecko (Oedura lesueurii), and potential implications for conservation of an endangered snake (Hoplocephalus bungaroides).

ABSTRACT: BACKGROUND: To conserve critically endangered predators, we also need to conserve the prey species upon which they depend. Velvet geckos (Oedura lesueurii) are a primary prey for the endangered broad-headed snake (Hoplocephalus bungaroides), which is restricted to sandstone habitats in southeastern Australia. We sequenced the ND2 gene from 179 velvet geckos, to clarify the lizards' phylogeographic history and landscape genetics. We also analysed 260 records from a longterm (3-year) capture-mark-recapture program at three sites, to evaluate dispersal rates of geckos as a function of locality, sex and body size. RESULTS: The genetic analyses revealed three ancient lineages in the north, south and centre of the species' current range. Estimates of gene flow suggest low dispersal rates, constrained by the availability of contiguous rocky habitat. Mark-recapture records confirm that these lizards are highly sedentary, with most animals moving < 30 m from their original capture site even over multi-year periods. CONCLUSION: The low vagility of these lizards suggests that they will be slow to colonise vacant habitat patches; and hence, efforts to restore degraded habitats for broad-headed snakes may need to include translocation of lizards.

The interplay of spatial and climatic landscapes in the genetic distribution of a South American parrot

AimOur aim was to understand the interplay of heterogeneous climatic and spatial landscapes in shaping the distribution of nuclear microsatellite variation in burrowing parrots, Cyanoliseus patagonus. Given the marked phenotypic differences between populations of burrowing parrots we hypothesized an important role of geographical as well climatic heterogeneity in the population structure of this species. LocationSouthern South America. MethodsWe applied a landscape genetics approach to investigate the explicit patterns of genetic spatial autocorrelation based on both geography and climate using spatial principal component analysis (sPCA). This necessitated a novel statistical estimation of the species climatic landscape, considering temperature- and precipitation-based variables separately to evaluate their weight in shaping the distribution of genetic variation in our model system. ResultsGeographical and climatic heterogeneity successfully explained molecular variance in burrowing parrots. sPCA divided the species distribution into two main areas, Patagonia and the pre-Andes, which were connected by an area of geographical and climatic transition. Moreover, sPCA revealed cryptic and conservation-relevant genetic structure: the pre-Andean populations and the transition localities were each divided into two groups, each management units for conservation. Main conclusionssPCA, a method originally developed for spatial genetics, allowed us to unravel the genetic structure related to spatial and climatic landscapes and to visualize these patterns in landscape space. These novel climatic inferences underscore the importance of our modified sPCA approach in revealing how climatic variables can drive cryptic patterns of genetic structure, making the approach potentially useful in the study of any species distributed over a climatically heterogeneous landscape.

Les impacts de la dispersion historique sur la variabilité génétique à différentes échelles spatiales : connaître l'histoire pour mieux comprendre le présent.

La variabilité génétique actuelle est influencée par une combinaison complexe de variables historiques et contemporaines. Dès lors, une interprétation juste de l’impact des processus actuels nécessite une compréhension profonde des processus historiques ayant influencé la variabilité génétique. En se basant sur la prémisse que des populations proches devraient partager une histoire commune récente, nombreuses études, effectuées à petite échelle spatiale, ne prennent pas en considération l’effet potentiel des processus historiques. Cette thèse avait pour but de vérifier la validité de cette prémisse en estimant l’effet de la dispersion historique à grande et à petite échelle spatiale. Le premier volet de cette thèse avait pour but d’évaluer l’impact de la dispersion historique sur la répartition des organismes à grande échelle spatiale. Pour ce faire, les moules d’eau douce du genre flotteurs (Pyganodon spp.) ont servies de modèle biologique. Les moules d'eau douce se dispersent principalement au stade larvaire en tant que parasites des poissons. Une série de modèles nuls ont été développés pour évaluer la co-occurrence entre des parasites et leurs hôtes potenitels. Les associations distinctes du flotteur de Terre-Neuve (P. fragilis) avec des espèces de poissons euryhalins permettent d’expliquer sa répartition. Ces associations distinctes ont également pu favoriser la différenciation entre le flotteur de Terre-Neuve et son taxon soeur : le flotteur de l’Est (P. cataracta). Cette étude a démontré les effets des associations biologiques historiques sur les répartitions à grande échelle spatiale. Le second volet de cette thèse avait pour but d’évaluer l’impact de la dispersion historique sur la variabilité génétique, à petite échelle spatiale. Cette fois, différentes populations de crapet de roche (Ambloplites rupestris) et de crapet soleil (Lepomis gibbosus), dans des drainages adjacents ont servies de modèle biologique. Les différences frappantes observées entre les deux espèces suggèrent des patrons de colonisation opposés. La faible diversité génétique observée en amont des drainages et la forte différenciation observée entre les drainages pour les populations de crapet de roche suggèrent que cette espèce aurait colonisé les drainages à partir d'une source en aval. Au contraire, la faible différenciation et la forte diversité génétique observées en amont des drainages pour les populations de crapet soleil suggèrent une colonisation depuis l’amont, induisant du même coup un faux signal de flux génique entre les drainages. La présente étude a démontré que la dispersion historique peut entraver la capacité d'estimer la connectivité actuelle, à petite échelle spatiale, invalidant ainsi la prémisse testée dans cette thèse. Les impacts des processus historiques sur la variabilité génétique ne sont pas faciles à démontrer. Le troisième volet de cette thèse avait pour but de développer une méthode permettant de les détecter. La méthode proposée est très souple et favorise la comparaison entre la variabilité génétique et plusieurs hypothèses de dispersion. La méthode pourrait donc être utilisée pour comparer des hypothèses de dispersion basées sur le paysage historique et sur le paysage actuel et ainsi permettre l’évaluation des impacts historiques et contemporains sur la variabilité génétique. Les performances de la méthode sont présentées pour plusieurs scénarios de simulations, d’une complexité croissante. Malgré un impact de la différentiation globale, du nombre d’individus ou du nombre de loci échantillonné, la méthode apparaît hautement efficace. Afin d’illustrer le potentiel de la méthode, deux jeux de données empiriques très contrastés, publiés précédemment, ont été ré analysés. Cette thèse a démontré les impacts de la dispersion historique sur la variabilité génétique à différentes échelles spatiales. Les effets historiques potentiels doivent être pris en considération avant d’évaluer les impacts des processus écologiques sur la variabilité génétique. Bref, il faut intégrer l’évolution à l’écologie.

Concordant Phylogeographies of 2 Malaria Vectors Attest to Common Spatial and Demographic Histories

The phylogeography of South American lineages is a topic of heated debate. Although a single process is unlikely to describe entire ecosystems, related species, which incur similar habitat limitations, can inform the history for a subsection of assemblages. We compared the phylogeographic patterns of the cytochrome oxidase I marker from Anopheles triannulatus (N = 72) and previous results for A. darlingi (N = 126) in a broad portion of their South American distributions. Both species share similar population subdivisions, with aggregations northeast of the Amazon River, in southern coastal Brazil and 2 regions in central Brazil. The average (ST) between these groups was 0.39 for A. triannulatus. Populations northeast of the Amazon and in southeastern Brazil are generally reciprocally monophyletic to the remaining groups. Based on these initial analyses, we constructed the a priori hypothesis that the Amazon and regions of high declivity pose geographic barriers to dispersal in these taxa. Mantel tests confirmed that these areas block gene flow for more than 1000 km for both species. The efficacy of these impediments was tested using landscape genetics, which could not reject our a priori hypothesis but did reject simpler scenarios. Results form summary statistics and phylogenetics suggest that both lineages originated in central Amazonia (south of the Amazon River) during the late Pleistocene (579 000 years ago) and that they followed the same paths of expansion into their contemporary distributions. These results may have implications for other species sharing similar ecological limitations but probably are not applicable as a general paradigm of Neotropical biogeography.