Aim: In theory, long-distance dispersal (LDD) outside a species’ range contributes to genetic divergence. However, previous studies have not discriminated this process from vicariant speciation in migratory bird species. We conducted an integrative phylogeographic approach to test the LDD hypothesis, which predicts that a Japanese migratory bird subspecies diverged from a population in the coastal region of the East China Sea (CRECS) via LDD over the East China Sea (ECS). Location: East Asia Taxon: Brown Shrike (Lanius cristatus) Methods: Both a haplotype network and a multi-locus gene network of its three subspecies were reconstructed to examine from which continental population the Japanese subspecies diverged. A species distribution model (SDM) for the Japanese subspecies was constructed using bioclimatic variables under the maximum entropy algorithm. It was projected to the climate of the last glacial maximum (LGM) to infer the candidate source area of colonisation. A migratory route of the Japanese subspecies, which possibly reflects a candidate past colonisation route, was tracked by light-level geolocators. Results: Molecular phylogenetic networks suggest that the Japanese subspecies diverged from a population in the CRECS. The SDM inferred that the emerged continental shelf of the ECS and the present CRECS were suitable breeding areas for the Japanese subspecies during the LGM. A major migratory route for the Japanese subspecies was inferred between the CRECS and the Japanese archipelago across the ECS. Main conclusions: Our integrative approach supported the LDD hypothesis for divergence of the Japanese subspecies of the Brown Shrike. Shrinkage and expansion of the ECS may have been responsible for successful colonisation and isolation of the new population. Vicariance was inferred for divergence of the subspecies in the northeast Asian continent from the Japanese population. Our framework provides a new phylogeographic scenario in this region, and discriminating LDD and vicariance models should improve our understanding of the phylogeographic histories of migratory species.
Comparisons of 3D shapes have recently been applied to diverse anatomical structures using landmarking techniques. However discerning evolutionary patterns can be challenging for structures lacking homologous landmarks. We used alpha shape analyses to quantify vaginal shape complexity in 40 marine mammal specimens including cetaceans, pinnipeds, and sirenians. We explored phylogenetic signal and the potential roles of natural and sexual selection on vaginal shape evolution. Complexity scores were consistent with qualitative observations. Cetaceans had a broad range of alpha complexities, while pinnipeds were comparatively simple and sirenians were complex. Intraspecific variation was found. Three-dimensional surface heat maps revealed that shape complexity was driven by invaginations and protrusions of the vaginal wall. Phylogenetic signal was weak and metrics of natural selection (relative neonate size) and sexual selection (relative testes size, sexual size dimorphism, and penis morphology) did not explain vaginal complexity patterns. Additional metrics, such as penile shape complexity, may yield interesting insights into marine mammal genital coevolution. We advocate for the use of alpha shapes to discern patterns of evolution that would otherwise not be possible in 3D anatomical structures lacking homologous landmarks.
1. Species distribution modelling, which allows users to predict the spatial distribution of species with the use of environmental covariates, has become increasingly popular, with many software platforms providing tools to fit species distribution models. However, the species observations used in species distribution models can have varying levels of quality and can have incomplete information, such as uncertain species identity. 2. In this paper, we develop two algorithms to reclassify observations with unknown species identities which simultaneously predict different species distributions using spatial point processes. We compare the performance of the different algorithms using different initializations and parameters with models fitted using only the observations with known species identity through simulations. 3. We show that performance varies with differences in correlation among species distributions, species abundance, and the proportion of observations with unknown species identities. Additionally, some of the methods developed here outperformed the models that didn't use the misspecified data. 4. These models represent an helpful and promising tool for opportunistic surveys where misidentification happens or for the distribution of species newly separated in their taxonomy.
1. Patterns in, and the underlying dynamics of, species cooccurrence is of interest in many ecological applications. Unaccounted for, imperfect detection of the species can lead to misleading inferences about the nature and magnitude of any interaction. A range of different parameterisations have been published that could be used with the same fundamental modelling framework that accounts for imperfect detection, although each parameterisation has different advantages and disadvantages. 2. We propose a parameterisation based on log-linear modelling that does not require a species hierarchy to be defined (in terms of dominance), and enables a numerically robust approach for estimating covariate effects. 3. Conceptually the parameterisation is equivalent to using the presence of species in the current, or a previous, time period as predictor variables for the current occurrence of other species. This leads to natural, ’symmetric’, interpretations of parameter estimates. 4. The parameterisation can be applied to many species, in either a maximum-likelihood or Bayesian estimation framework. We illustrate the method using camera trapping data collected on three mesocarnivore species in South Texas.
1. A time-consuming challenge faced by ecologists is the extraction of meaningful data from camera trap images to inform ecological management. Automated object detection solutions are increasingly, however, most are not sufficiently robust to be deployed on a large scale due to lack of location invariance across sites. This prevents optimal use of ecological data and results in significant resource expenditure to annotate and retrain object detectors. 2. In this study, we aimed to (a) assess the value of publicly available image datasets including FlickR and iNaturalist (FiN) when training deep learning models for camera trap object detection (b) develop a for training location invariant object detection models and (c) explore the use of small subsets of camera trap images for optimization training. 3. We collected and annotated 3 datasets of images of striped hyena, rhinoceros and pig, from FiN, and used transfer learning to train 3 object detection models in the task of animal detection. We compared the performance of these models to that of 3 models trained on the Wildlife Conservation Society and Camera CATalogue datasets, when tested on out of sample Snapshot Serengeti datasets. Furthermore, optimized the FiN models via infusion of small subsets of camera trap images to increase robustness for challenging detection cases. 4. In all experiments, the mean Average Precision (mAP) of the FiN models was significantly higher (82.33-88.59%) than that achieved by the models trained only on camera trap datasets (38.5-66.74%). The infusion of camera trap images into FiN training further improved mAP, with increases ranging from 1.78-32.08%. 5. Ecology researchers can use FiN images for training robust, location invariant, out-of-the-box, deep learning object detection solutions for camera trap image processing. This would allow AI technologies to be deployed on a large scale in ecological applications. Datasets and code related to this study are open source and available at: https://github.com/ashep29/infusion
For nearly all organisms, dispersal is a fundamental life history trait that can shape their ecology and evolution. Variation in dispersal capabilities within a species exists and can influence population genetic structure and ecological interactions. In fungus-gardening (attine) ants, co-dispersal of ants and mutualistic fungi is crucial to the success of this obligate symbiosis. Female-biased dispersal (and gene flow) may be favored in attines because virgin queens carry the responsibility of dispersing the fungi, but a paucity of research has made this conclusion difficult. Here, we investigate dispersal of the fungus-gardening ant Trachymyrmex septentrionalis using a combination of maternally- (mitochondrial DNA) and biparentally-inherited (microsatellites) markers. We found three distinct, spatially isolated mitochondrial DNA haplotypes. Two were found in the Florida panhandle and the other was found in the Florida peninsula. In contrast, biparental markers illustrated significant gene flow across this region and minimal spatial structure. The differential patterns uncovered from mitochondrial DNA and microsatellite markers suggest that most long-distance ant dispersal is male-biased and that females (and concomitantly the fungus) have more limited dispersal capabilities. Consequently, the limited female dispersal is likely an important bottleneck for the fungal symbiont. This bottleneck could slow fungal genetic diversification, which has significant implications for both ant hosts and fungal symbionts regarding population genetics, species distributions, adaptive responses to environmental change, and coevolutionary patterns.
We examined whether the presence or absence of different environmental stressors influenced the reproductive potential of a saltmarsh species - Plantago maritima. We focused on total seed output, seed quality and biomass of progeny. So far, there are no studies trying to answer the question of how different saltmarsh management affects the quality of seed in saltmarsh species. For the purposes of the study, plots subjected to light mowing, light or heavy grazing, trampling or rooting were designated in three nature reserves in Poland. On each plot, the abundance of infructescences per sq. metre was calculated. Mature infructascences were collected and their length and no of fruit capsules were measured. The seeds obtained from fruit capsules were weighted and sown in controlled conditions. The germination rate and the final germination percentage were calculated. A representative number of sprouts were grown. After a period of two months, the specimens were harvested and their total dry mass was measured. It was found that heavy grazing had the greatest effect on all of the studied characteristics. The presence of this factor resulted in shorter infructescences with a smaller number of fruit capsules. However, this phenomenon was compensated by the higher abundance of infructescences per sq. metre. At the same time seeds produced by grazed specimens were significantly lighter. Intensive trampling by people affected specimens in a similar way to heavy grazing, while mowing and rooting had less impact on the considered characteristics. Although a positive correlation between seed mass and germination success was found, the altogether lower seed mass had a negligible effect on germination parameters. Also, the differences in seed parameters did not affect dry mass of obtained progeny grown in lab conditions.
A rapid biodiversity assessment of insects and associated Laboulbeniales fungi was conducted over the course of five nights in August, 2018 at two central Florida lakes: Lake Eustis and the nearby protected and restored National Natural Landmark, Emeralda Marsh Conservation Area (EMCA), which encompasses a portion of Lake Griffin. These locations were selected because Lake Eustis was surveyed for Laboulbeniales in 1897 by mycologist Dr. Roland Thaxter, but has not since been investigated. Because Lake Eustis has been urbanized, with the lake perimeter almost entirely altered by human development, the site offers a look into Laboulbeniales diversity across a 121 year timeline, before and after human development. By surveying Lake Eustis and EMCA, a modern case study comparison of Laboulbeniales and insect diversity between a developed and unrestored system and a protected and restored system is made. A total of 4,022 insects were collected during the rapid assessment. Overall, insect abundance was greater at EMCA, with 3,001 insects collected, compared to 1,021 insects collected from Eustis. Though family level insect richness was comparable between sites, with 55 families present at EMCA and 56 at Eustis, 529 out of 3,001 (17.6%) of the insects collected at EMCA were hosts to parasitic Laboulbeniales fungi whereas only 2 out of 1,021 (0.19%) collected from Eustis were infected. There were 16 species of Laboulbeniales found at EMCA compared to only one at Eustis. The current number of Laboulbeniales species documented at Eustis was incredibly depauperate compared to the 27 species recorded by Thaxter in 1897, suggesting the possibility of utilizing Laboulbeniales as indicators of ecosystem health. A figure displaying host-parasite records and a species list of Laboulbeniales fungi is compiled and updated occurrence records for species of Ceratomyces and Hydrophilomyces are provided.
The range of the Canada lynx (Lynx canadensis) has contracted substantially from its historical range prior to the 19th century. Using harvest records, we found that the southern range of the lynx in Ontario in the late 1940s collapsed and then in a short period of time increased to its largest extent in the mid-1960s where the lynx range spread south of the boreal forest for a decade. After this expansion the southern range contracted northwards beginning in the 1970s. Most recently, there has been a slight expansion between 2010-2017. We have attributed these dynamics on the southern range periphery to the fluctuation of the boreal lynx population in the core of the species’ range. In addition, connectivity to boreal lynx populations and snow depth seemed to condition whether the lynx expanded into an area. However, we did not find any evidence that would suggest that these changes were due to anthropogenic disturbances or competition. The boreal lynx population does not reach numbers as it once did, consequently we likely will not see large expansions of the southern lynx range as in the mid-1960s. Our results suggest that southern lynx range in Ontario have been driven by the magnitude of the boreal lynx population cycle, connectivity to the boreal forest and snow conditions. Therefore, it is quite unlikely that southern lynx population in the Great Lakes will ever recover, since the warming climate and forestry practices are causing a northward contraction of the boreal forest and likely with it the core lynx populations.
Third generation sequencing technologies, such as Oxford Nanopore Technologies (ONT) and Pacific Biosciences (PacBio), have gained popularity over the last years. These platforms can generate millions of long read sequences. This is not only advantageous for genome sequencing projects, but also for amplicon-based high-throughput sequencing experiments, such as DNA barcoding. However, the relatively high error rates associated with these technologies still pose challenges for generating high quality consensus sequences. Here we present NGSpeciesID, a program which can generate highly accurate consensus sequences from long-read amplicon sequencing technologies, including ONT and PacBio. The tool includes clustering of the reads to help filter out contaminants or reads with high error rates and employs polishing strategies specific to the appropriate sequencing platform. We show that NGSpeciesID produces consensus sequences with improved usability by minimizing preprocessing and software installation and scalability by enabling rapid processing of hundreds to thousands of samples, while maintaining similar consensus accuracy as current pipelines
The Magdalena River in Colombia and its floodplain lakes are key ecosystems for the subsistence of Colombian society. Yet, hydrologic regulation, pollution, deforestation and climate change are threatening its ecological integrity. To understand how these floodplain lakes respond over decadal-centennial scales to natural and anthropogenic stressors, we selected two shallow lakes with varying degree of connectivity to the River and assessed their historical ecological and limnological change through a multi-proxy analysis of diatoms, geochemistry and lake’s morphometric variation resulting from extreme periods of high floods and droughts. The reconstruction of the more isolated San Juana Lake covered the last c.500 years. It showed riverine-flooded conditions from c.1555-1741 characterised by high detrital inputs, reductive conditions, and dominance of planktonic diatoms. From c.1758-1954, the riverine meander became disconnected, conveying into a marsh-like environment rich in aerophil diatoms and organic matter. The current lake was then formed around the mid 1960s and a diverse lake-associated diatom flora developed. Lake waters became more oxygenated, while sedimentation and nutrients increased through time since the lake formation. The reconstruction for Barbacoas Lake, a waterbody directly connected to the Magdalena River, spanned the last 60 years and showed alternating riverine-wetland-lake conditions dominated by planktonic and benthic diatoms respectively. An exception was however observed, during a prolonged period of low rainfall between 1989-1992, where the lake almost desiccated and where aerophil diatoms prevailed. Inferences of flood magnitudes and river connectivity in the lakes were supported by parallel increases in Zr/Fe (flooding) and detrital inputs (Ti/Ca) along with decreases in sedimentary OM. We proposed that lake hydrological connectivity to the Magdalena River is a main factor controlling lake long-term responses to human pressures. Highly connected lakes may respond more acutely to ENSO events while isolated lakes might be more sensitive to local land-use changes.
1. Understanding how ecosystem engineers influence other organisms has long been a goal of ecologists. Woodpeckers select nesting sites with high food availability and will excavate and then abandon multiple cavities throughout their lifetime. These cavities are crucial to secondary cavity nesting birds (SCB) that are otherwise limited by the availability of naturally occurring cavities. 2. Our study examined the influence of food resources on the nest site location and home-range size of woodpeckers, and the subsequent influence of woodpeckers on the nesting success of SCB. 3. Using five years of avian point count data to locate golden-fronted woodpeckers (GFWO), we correlated insect availability with GFWO home range size, determined differences in insect availability between GFWO occupied and unoccupied sites, and compared nesting success for the GFWO and common SCB in south Texas. We used model averaging to fit species-specific logistic regression models to predict nest success based on cavity metrics across all species. 4. Sites occupied by GFWO had a higher biomass of insects in orders Coleoptera, Hymenoptera, and Orthoptera than unoccupied sites, and there was a negative correlation between the availability of these insect orders and home-range size. GFWO nest success increased with vegetation cover and lower levels of tree decay. SCB had higher levels of nesting success in abandoned GFWO, and in trees with lower levels of nest tree decay. 5. Our results suggest that SCB may be drawn to nest in abandoned woodpecker cavities where they have higher rates of nest success compared to natural cavities. Additionally, the prevalence for GFWO to excavate cavities in trees with lower levels of decay contradicts previous literature, and may indicate a novel temperature trade-off, with live trees requiring more energy to excavate, but providing increased protection from high breeding season temperatures in arid and semi-arid areas.
Tilia cordata Mill. is a valuable tree species enriching the ecological values of the coniferous dominated boreal forests in northerly Europe. Following the historical decline, spreading of Tilia sp. is challenged by the elevated inbreeding and habitat fragmentation. We aimed to identify the main factors affecting the genetic potential of Tilia cordata for natural expansion by studying the geographical distribution of genetic diversity of Tilia cordata in semi-boreal forests of Lithuania. We used 14 genomic microsatellite markers to genotype 543 individuals from 23 wild growing populations of Tilia cordata in Lithuania. We found that Tilia cordata retained high levels of genetic diversity (population Fis = 0 to 0.15, Ho = 0.53 to 0.69, He = 0.56 to 0.75). AMOVA, Bayesian clustering and Monmonier’s barrier detection indicate weak but significant differentiation among the populations (Fst = 0.037***) into geographically interpretable clusters of (a) western Lithuania with high genetic heterogeneity but low genetic diversity, bottleneck effects, (b) peaking values of genetic diversity of Tilia cordata on rich and most soils of midland lowland, and (c) the most differentiated populations on poor soils of the coolest north-eastern highland possessing the highest rare alle frequency but elevated inbreeding and bottleneck effects, presumably, due to sub-structuring. We conclude that the genepool of Tilia cordata in Lithuania contains (a) the autochthonous populations of high genetic diversity representing the pre-historical genepools, that can be promoted, and (b) the escapes from urban sources of low diversity, that must be contained.
Exploring vegetation distribution spatial patterns facilitates understanding how biodiversity addresses the potential threat of future climate variability, especially for highly diverse and threatened tropical plant communities, but few empirical studies have been performed over various environmental scales. In this study, we used species-based and phylogeny-based methods to analyze the α- and β-diversity pattern variation in Dacrydium pectinatum communities and its key drivers along elevation and geographical gradients across three national nature reserves in Bawangling, Diaoluoshan and Jianfengling. Our study indicates that the species and phylogenetic α-diversity patterns presented consistent decreasing with elevation, with the peak occurring at low elevations. Environmental filtering caused by decreases in limiting factors, such as temperature, precipitation, soil organic matter, soil phosphorus, and light, is the main reason for the decline in diversity at high elevations, whereas low-elevation areas are affected by various factors, such as environmental filtering and similarity limitations. Species and phylogenetic β-diversity changes are closely related to environmental filtering and dispersal limitation, but the latter is key in community assembly at the heterogeneous spatial scale. In conclusion, combining species-based and phylogeny-based methods to explore the biogeographic patterns of tropical plant communities helps provide convincing evidence and confirms that the relative contributions of niche and neutral theory in the assembly process vary along environmental gradients. Though the D. pectinatum community constitute a floristically integrated unit, the genetic relationships between species are relatively far, and co-evolution to promote species coexistence is difficult when faced with habitat pressure. Hence, we believe that species coexistence in tropical plant communities requires mild environmental conditions, and low temperatures, precipitation, soil nutrients and light will aggravate environmental filtering and species competition. We also recommend strengthening the construction and management of nature reserves and the exploration of biodiversity formation mechanisms, which are crucial for biodiversity conservation in endangered tropical plant communities.
Understanding genetic variation and structure, adaptive genetic variation and its relationship with environmental factors is of great significance to understand how plants adapt to climate change and design effective conservation and management strategies. The objective of this study was to (I) investigate the genetic diversity and structure by AFLP markers in 36 populations of R. aureum from northeast China, (Ⅱ) reveal the relative contribution of geographical and environmental impacts on the distribution and genetic diﬀerentiation of R. aureum; (Ⅲ) identify outlier loci under selection and evaluate the association between outlier loci and environmental factors and (Ⅳ) exactly calculate development trend of population of R. aureum，as it is confronted with severe climate change and to provide information for designing eﬀective conservation and management strategies. We found high genetic variation (I = 0.584) and diﬀerentiation among populations (ΦST = 0.711) and moderate levels of genetic diversity within populations of R. aureum. A significant relationship between genetic distance and environmental distance was identified, which suggested that the differentiation of diﬀerent populations was the caused by environmental factors. Using BayeScan and Dfdist, 42 outlier loci identified and most of the outlier loci are associated with climate or relief factors, suggesting that these loci are linked to genes that are involved in the adaptability of R. aureum to environment. Species distribution models (SDM) showed that climate warming will cause a significant reduction of suitable area for R. aureum especially under the RCP 85 scenario. Our results help to understand the potential response of R. auruem to climatic changes, and provide new perspectives for R. auruem resource management and conservation strategies.
How and why species range size varies along spatial gradients is fundamental yet controversial topics in biogeography. To advance our understanding on these questions and to provide insight into biological conservation, we assessed the elevational variations in vascular plants range size for different life form and biogeographical affinities, and explored the main drivers underlying above variations in the longest valley in China's Himalayas---the Gyirong Valley. Elevational range sizes of vascular plants were documented by 96 sampling plots along 12 elevational bands of 300-m ranging from 1800 to 5400 m above sea level. We assessed the elevational variations in range size by averaging the range size of all species within each elevational band. We then related range size to climate, disturbance, competition factors and the mid-domain effect, and explored the relative importance of aforementioned factors in explaining the range size variations using the Random Forest model. Total 545 vascular plants were documented by our sampling plots along the elevational gradient. Out of 545 plants, 158, 387, 337 and 112 were woody, herbaceous, temperate and tropical species respectively. Range size of each groups of vascular plants shown uniform increasing trends along the elevational gradient which are in accordance with the prediction of Rapoport's rule. Climate was the main driver for the increasing trends of vascular plants range size in the Gyirong Valley. Climate variability hypothesis and mean climate condition hypothesis were both supported to jointly explain such climate-range size relationship. Our results reinforce previous notion that Rapoport's rule applies to where the influence of climate is most pronounced, and call for close attention to the impact of climate change in order to prevent range contraction and even extinction under global warming.
The mechanisms of forest seed dispersal and regeneration of various altered forest ecosystems are complemented by the action of carnivores. The objective of this study was to evaluate the role of endozoochoric and diploendozoochoric mammals in the dispersal, scarification, and germination of seeds in two different forest ecosystems: temperate forest (TF) and dry tropical forest (DTF). A direct search and scat collection were carried out to determine dispersing agents and the abundance and richness of seeds in the Protected Natural Area, Sierra Fría, Aguascalientes, Mexico. Viability was evaluated by means of X-rays and a germination test using an optical densitometry. In addition, thickness measurements and observations were made on the surfaces of the testas by a scanning electron microscopy. In the TF, four plant species were dispersed, mainly Arctostaphylos pungens (P < 0.05), by four mammals, where the gray fox dispersed the highest average (66.8 ± 68.2; P < 0.05) and diploendozoochory was detected in bobcat scats associated with rabbit hair (Sylvilagus floridanus). The DTF presented higher abundance and richness of species, where Myrtillocactus geometrizans had the highest abundance (2680 ± 4423) and the coati (P < 0.05) dispersed the highest number of seeds (8639 ± 12203). In both types of forest, endozoochory and diploendozoochory did not affect the viability, the thickness of the testas, or the germination of any species of seeds. These results suggest that dispersing carnivores adapt to the abundance and richness of seeds in the forests they inhabit, thus developing an important ecological function by dispersing, scarifying and promoting the selective germination of seeds with thick testas in TF and thin testas in DTF.
Diapause is an adaptive dormancy strategy by which arthropods endure extended periods of adverse climatic conditions. Seasonal variation in larval diapause initiation and duration in Ostrinia furnacalis influences adult mating generation number (voltinism) across different local environments. The degree of mating period overlap between sympatric voltine ecotypes influence hybridization level, but impacts on O. furnacalis population genetic structure and evolution of divergent adaptive phenotypes remains uncertain. Genetic differentiation was estimated between voltine ecotypes collected from 8 locations. Mitochondrial haplotypes were significantly different between historically allopatric univoltine and bivoltine locations. Haplotypes from sympatric locations were clustered more-closely to bivoltine locations, but influenced by local demographics. Additionally, analyses of single nucleotide polymorphism (SNP) genotypes implicate voltinism, as opposed to geographic distance, as contributing to low, but significant levels of variation among voltine ecotypes. Regardless, only 11 of 257 SNP loci were predicted to be under selection, suggesting population genetic homogenization except at loci proximal to factors putatively responsible for locally adaptive or voltinism-specific traits. These findings provide evidence that divergent voltine ecotypes may be maintained in allopatric and sympatric areas despite relatively high rates of nuclear gene flow, yet influence of voltinism on maintenance of observed haplotype divergence remains unresolved.