Next Article in Journal
The Prediction of Clinical Mastitis in Dairy Cows Based on Milk Yield, Rumination Time, and Milk Electrical Conductivity Using Machine Learning Algorithms
Previous Article in Journal
Exploring the Rumen Microbiota and Serum Metabolite Profile of Hainan Black Goats with Different Body Weights before Weaning
Previous Article in Special Issue
Assessing the Seasonal and Spatial Dynamics of Zooplankton through DNA Metabarcoding in a Temperate Estuary
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Small-Mammal Genomics Highlights Viaducts as Potential Dispersal Conduits for Fragmented Populations

1
Asian School of the Environment, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798, Singapore
2
Department of Biological Sciences, National University of Singapore, 16 Science Drive 4, Singapore 117558, Singapore
3
Discipline of Biological Sciences, School of Natural Sciences, University of Tasmania, Hobart, TAS 7005, Australia
4
Faculty of Science and Marine Environment, Universiti Malaysia Terengganu, Kuala Terengganu 21030, Terengganu, Malaysia
*
Author to whom correspondence should be addressed.
Animals 2024, 14(3), 426; https://doi.org/10.3390/ani14030426
Submission received: 27 December 2023 / Revised: 15 January 2024 / Accepted: 24 January 2024 / Published: 28 January 2024
(This article belongs to the Special Issue New Tools for Monitoring Genetic Diversity in Animals)

Abstract

:

Simple Summary

Wildlife crossings are often constructed to enhance genetic connectivity among populations divided by roads (including highways). However, few studies have demonstrated the efficacy of viaducts in counteracting the barrier effects imposed by roads. We measured genetic diversity and divergence in four small mammal species commonly found in rainforests in Malaysia—Tupaia glis, Maxomys rajah, M. whiteheadi, and Niviventer cremoriventer—across three treatment types: (1) viaduct sites, at which sampling locations were separated by a highway but connected by a vegetated viaduct; (2) non-viaduct sites, at which sampling locations were separated by a highway and not connected by a viaduct; and (3) control sites, at which there was no road or highway fragmenting the forest. We found that viaducts facilitated movement in small ground-dwelling species such as M. whiteheadi and also when existing highways were relatively wide. However, despite the potential for viaducts to facilitate movement and therefore increase genetic connectivity in M. whiteheadi, the genetic distance in populations at viaduct sites was still greater than at control and/or non-viaduct sites for the other three species. Our findings highlight the importance of maintaining intact forests rather than relying solely on the construction of viaducts to connect fragmented populations.

Abstract

Wildlife crossings are implemented in many countries to facilitate the dispersal of animals among habitats fragmented by roads. However, the efficacy of different types of habitat corridors remains poorly understood. We used a comprehensive sampling regime in two lowland dipterocarp forest areas in peninsular Malaysia to sample pairs of small mammal individuals in three treatment types: (1) viaduct sites, at which sampling locations were separated by a highway but connected by a vegetated viaduct; (2) non-viaduct sites, at which sampling locations were separated by a highway and not connected by a viaduct; and (3) control sites, at which there was no highway fragmenting the forest. For four small mammal species, the common tree shrew Tupaia glis, Rajah’s spiny rat Maxomys rajah, Whitehead’s spiny rat Maxomys whiteheadi and dark-tailed tree rat Niviventer cremoriventer, we used genome-wide markers to assess genetic diversity, gene flow and genetic structure. The differences in genetic distance across sampling settings among the four species indicate that they respond differently to the presence of highways and viaducts. Viaducts connecting forests separated by highways appear to maintain higher population connectivity than forest fragments without viaducts, at least in M. whiteheadi, but apparently not in the other species.

1. Introduction

Habitat fragmentation, degradation, and loss pose the most significant threats to the structure and persistence of animal populations and communities [1,2]. Fragmentation is most rapid in developing countries where the expansion of road networks is increasing due to competing land uses such as farming, manufacturing and housing [3]. Barriers that bisect continuous habitat, particularly roads (including highways), initiate the process of habitat fragmentation and can restrict or eliminate animal movement through a landscape [4], with concomitant consequences for connectivity and gene flow [5]. The reduced connectivity may lead to a decrease in viability and persistence of isolated populations [6]. However, the extent to which habitat fragmentation has a negative effect on the genetic structure and persistence of animal populations remains debated. While some studies have demonstrated a negative impact of habitat fragmentation on species’ population genetic structure [5,7], others have failed to detect these effects [8].
Wildlife corridors are widely understood to connect habitat fragments and mediate the effects of fragmentation [9]. They include naturally occurring linear habitats such as riparian reserves and hedgerows, purpose-built structures such as wildlife overpasses and underpasses, and incidental structures such as drainage culverts. Corridors have been widely advocated as essential components of reserve design because they can connect isolated areas of suitable habitat and thus minimise the harmful effects of habitat fragmentation on animal movement [10,11,12,13]. Corridors are predicted to benefit populations in patchy habitats by promoting movement, which increases population densities, gene flow, and recolonisation of extinct patch populations [14]. However, the efficacy of such passages remains largely untested [15,16,17], and there has been much debate about their effectiveness in connecting isolated populations [18,19]. Corridors have been shown to increase connectivity, maintain biodiversity [20], increase population sizes [21,22,23], facilitate movement between fragmented patches [10,18,20,24,25,26,27,28,29,30] and promote gene flow [28,31,32]. Other studies, however, have found no significant effects of corridors or even negative impacts [33,34,35].
To demonstrate the effectiveness of wildlife corridors in increasing connectivity between populations, studies of genetic relatedness are recommended. Molecular techniques such as next-generation sequencing, which has recently been used to reveal fine-scale population structure [36,37], offer new promise in investigating the influence of fragmentation and barriers on population connectivity [38,39]. Genetic methods allow us to measure average migration rates over time, which reveals the effects of fragmentation over several generations and is not as sensitive to current population sizes as mark-recapture studies are (e.g., when populations are extremely low, mark-recapture studies may be impossible) [40]. In addition, molecular techniques measure effective dispersal, the amount of gene flow between populations [41]. As genetic techniques use a single temporal sample per population to estimate migration rather than multiple samples, these techniques require less field effort than mark-recapture [28].
The consequences of habitat fragmentation on dispersal and genetic diversity are still largely unknown for non-volant small mammals, especially in some of the world’s equatorial rainforest areas, such as peninsular Malaysia in Southeast Asia. The relative scarcity of data on the effects of fragmentation on small mammals in Southeast Asia is alarming, as these animals provide important ecosystem functions and services for their natural habitats [42]. They are important seed dispersers, pollinators, invertebrate and seed predators, as well as prey for larger predators. Small mammal communities provide a good model for studying such impacts because species in these communities generally use a wide variety of resources, have short generation times that allow for quick detection of environmental change, may be permanent residents of a site, and usually respond to disturbances in a perceptible and measurable way [43]. Smaller mammals are thought to be particularly susceptible to fragmentation due to their limited ability to travel over long distances through exposed habitats [44]. They are thus likely to suffer severe impacts of fragmentation and would benefit from any increases in connectivity brought about by the construction of wildlife corridors.
In this study, we assessed the constraints in migration and gene flow due to habitat fragmentation caused by the construction of highways, and the effectiveness of wildlife underpasses, known in Malaysia as eco-viaducts, in facilitating movement and genetic connectivity in small mammals by comparing genetic distances (a measure of genetic differences, computed by using allele frequency data from many different loci) between individuals at viaduct sites, non-viaduct sites and control sites. We assumed that genetic relationships at all sites would have been at similar levels of divergence prior to fragmentation. We hypothesized that: (1) for forest species, the genetic distance between individuals should be lowest for populations in intact forest and highest for populations fragmented by highways; (2) if eco-viaducts are effective in maintaining population connectivity, the genetic distance between individuals connected by an eco-viaduct will more closely resemble the distance between individuals in intact forest; (3) alternatively, if eco-viaducts are not effective, then pairwise genetic distances should be similar to those in pairs separated by highways. By examining the population genetic structure of small mammal species across different spatial settings in this study, we aim to determine the effectiveness of eco-viaducts across highways in maintaining genetic linkage in a fragmented landscape.

2. Materials and Methods

2.1. Study Sites

We conducted this study in Kenyir, Terengganu (Figure 1a, elevation 100–300 m) and Sungai Yu, Pahang (Figure 1b, elevation 130–210 m), Peninsular Malaysia. Terrain is hilly and consists mostly of lowland dipterocarp forest. Rainfall averages 3000 mm per year with a pronounced wet season from November to March. Flooding is common during this period. The study areas are gazetted as forest reserves and can be logged under permit. Kenyir and Sungai Yu adjoin the Taman Negara National Park, Malaysia’s first national park, to the north and west, respectively (Figure 1c). Kenyir and Sungai Yu have rich biodiversity, but are prone to illegal logging, conversion to plantations, and poaching as they do not have the same protection status as a national park. Kenyir forest is bisected by federal route 185 (henceforth highway 185) (Figure 1a) while Sungai Yu forest is bisected by federal route 34 (henceforth highway 34) (Figure 1b). Highway 185 is a two-lane single carriageway (one lane in each direction) with a width of about 8 m and shoulder width of 2 m on either side. Highway 34 is a four-lane dual carriageway (two lanes in each direction) with a width of about 20 m and shoulder width of about 1–2 m on either side. There was no fencing along these highways, although vehicle guard rails were present in certain sections. The speed limit was 90 km/h on highway 185 and 110 km/h on highway 34. Average (±s.e.) vehicular traffic (counted at two different points on each highway between 0700 and 2300 h on six separate days) was 23 ± 5 cars, 16 ± 5 motorcycles and 10 ± 4 heavy vehicles per hour on highway 185; and 52 ± 11 cars, 24 ± 8 motorcycles and 20 ± 6 heavy vehicles per hour on highway 34. Many wildlife crossings, primarily eco-viaducts, have been constructed across Malaysia to restore connectivity between highway-bisected forest fragments, including forests in Kenyir and Sungai Yu. Eco-viaducts are bridge-like elevated roads (and highways) that allow passageway beneath for wildlife to safely cross between forests on either side of the highway (Figure 1d,e). Three eco-viaducts have been constructed across highway 185 in Kenyir and another three across highway 34 in Sungai Yu. The Kenyir eco-viaducts measuring 245 m, 140 m and 245 m in length were completed in 2008 (Figure 1a), while the Sungai Yu eco-viaducts measuring 80 m, 300 m and 1000 m were completed in 2014 (Figure 1b).

2.2. Small Mammal Trapping

Eighteen sites, six of them being non-viaduct highway-side sites, another six viaduct sites and six control sites (>500 m from any roads and highways), were selected for this study (Figure 1a,b). Both Kenyir and Sungai Yu had a set of nine sites each (Figure 1a,b). Each site was characterized by a pair of grids with 20 traps, each grid within a pair on opposite sides of the highway. The 20 traps within a grid were composed of ten Elliott sheet metal traps (32 × 10 × 10 cm) and ten Tomahawk wire cage traps (48 × 15 × 15 cm) (Figure S1) (except the 80 m wide viaduct in Sungai Yu which had only ten traps as there was not enough width across the viaduct to place 20 traps 10 m apart). Alternating Elliott and cage traps were set in parallel grids, with 10 m between traps (Figure 2). Paired grids were 50 m apart from each other, separated by the highway. The distance between grids was roughly chosen to be identical to the distance required for a crossing of the highway and its verges and ditches to confirm whether the highway accounts for any inhibition of movement. To avoid pseudoreplication issues, pairs of grids were at least 500 m apart, which is more than the home range of the small mammals in our study [45]. Trapping was conducted four times at Kenyir in 2017 and four times at Sungai Yu in 2018 from March to November, during the drier inter-monsoon and southwest monsoon seasons. Traps were set, checked and rebaited with bananas and peanut butter (cage traps) and vanilla scented oats (Elliott traps) [46] during five consecutive mornings and evenings to assess the diversity and abundance of both nocturnal and diurnal mammals. Captured individuals were identified to species, sex, age class and reproductive condition, weighed, measured and ear tagged before being released at the trap site. We collected ear clips for DNA tissue sampling from all animals trapped and tagged during live trapping using a 2 mm ear punch. Tissue samples were stored in absolute ethanol.
In total, we trapped 448 individuals from 17 species. The species Whitehead’s spiny rat Maxomys whiteheadi (Figure 3a), Rajah’s spiny rat Maxomys rajah (Figure 3b), dark-tailed tree rat Niviventer cremoriventer (Figure 3c) and common tree shrew Tupaia glis (Figure 3d) were included in this study to capture a wide variety of biological and ecological characteristics. These four small mammal species are found in tall lowland forests and forest edge, and feed on insects and plant matter such as fruits and seeds. M. rajah and T. glis are similar in size and weight (100–200 g), while M. whiteheadi (35–80 g) and N. cremoriventer (50–100 g) are smaller. M. whiteheadi and M. rajah are both nocturnal ground dwelling rats; N. cremoriventer is a nocturnal tree rat which is a good climber and lives both arboreally and on the ground, T. glis is a diurnal treeshrew active on the ground and in the understory. Samples from 32 individuals (7 in Kenyir, 25 in Sungai Yu) of M. rajah, 39 (16 in Kenyir, 23 in Sungai Yu) of M. whiteheadi, 22 (20 in Kenyir, 2 in Sungai Yu) of N. cremoriventer and 45 (25 in Kenyir, 20 in Sungai Yu) of T. glis were selected based on their occurrences in pairs of sites (Table S1).

2.3. DNA Extraction and ddRAD-Seq Library Preparation

DNA extractions were performed using the DNeasy Blood & Tissue Kit (QIAGEN, Hilden, Germany) following the manufacturer’s protocol for tissue. We prepared two libraries following Ng et al.’s [47] double-digest restriction-associated DNA sequencing (ddRADseq) protocol using EcoRI and MspI. To select for 250–600 bp fragments, as well as for the clean-up steps, we used Sera-Mag magnetic beads (Thermo Scientific, Waltham, MA, USA). DNA quantifications were performed with a Qubit 2.0 High Sensitivity DNA Assay (Invitrogen, Waltham, MA, USA). Before pooling samples, we checked DNA fragment size using a Fragment Analyser (Advanced Analytical Technologies, Inc., Ames, IA, USA). The two libraries were then spiked with 5% PhiX to prevent low nucleotide diversity issues from affecting the quality of the data, and were subsequently sequenced on an Illumina HiSeq 4000 platform (150 bp paired-end run).
Reads were demultiplexed and trimmed to 145 bp with process_radtags in STACKS 1.42 [48]. Reads with one or more uncalled bases were removed. For reference-based identification of ddRAD loci, we first aligned the demultiplexed reads of T. glis to the closely related Tupaia chinensis (GCA_000334495.1) [49], using the Burrows-Wheeler Aligner (BWA) [50] to index this reference genome. We used samtools 1.3.1 to convert the sam files to bam files, sort the aligned reads according to coordinates and filter files with a minimum required mapping quality score (MAPQ) score of 20 [51]. To call and filter single nucleotide polymorphisms (SNPs), we used ref_map.pl and population in STACKS 1.42 [48] for T. glis. In population, we set stack depth to 10 and the percentage of individuals represented at each locus to 0.9 and admitted only one random SNP per locus to preclude analysis of linked SNPs. For the other three species, there was no suitable reference genome, so SNPs were called de novo. In ustacks, we set the maximum distance (in nucleotides) allowed between stacks to 2, minimum depth of coverage required to create a stack to 3 and maximum distance allowed to align secondary reads to primary stacks to 4. We checked for SNPs under selection using BayeScan 2.1 [52] and used Plink 1.90 to remove linked loci and to calculate the level of missing data [53]. We allowed a variety of filters (0% or 10% missingness, including or excluding linked loci of r2 > 0.5, and including or excluding minor allele frequency < 5%) to generate eight datasets for the preliminary testing to rule out potential sampling artefacts (e.g., non-random distribution of genotypes) (Table S2).

2.4. Statistical Analyses

To visualise genetic differentiation amongst individuals and identify potential population subdivision, we performed principal component analysis (PCA) using the R package SNPRelate 1.6.6 [54] for all eight datasets (four species at two study areas). We carried out sensitivity analysis, checking across various settings of missingness, linkage and minor allele frequencies by confirming the consistency of PCAs across eight datasets for all four species. As a consequence of this sensitivity analysis, we selected the dataset with 10% missingness, in which we excluded linked loci but included minor alleles for all subsequent analyses (Table S2). To understand the genetic differentiation between the two study areas, we calculated the Weir-Cockerham’s FST between the two study areas using VCFtools v4.1 [55] and individual-pairwise relatedness using maximum likelihood estimation as implemented in SNPRelate.
We calculated pairwise genetic distances between all individuals within a species with the R package poppr [56]. To evaluate the efficacy of viaducts in facilitating the dispersal of small mammals and test the barrier effects of highways, Kruskal–Wallis tests were conducted on the genetic distances of individuals between pairs of grids within a site for each species across the three treatment types (viaduct, non-viaduct and control) and two study areas. This means that every individual of one species on one side of the highway or pair (in the case of control sites) is compared with every other individual of the same species on the other side of the highway or pair in the same pair of grids (site).
We also compared pairwise genetic distances among all individuals within a species in relation to different spatial distances to determine how resistance to dispersal changes in different landscape types. Using the Least Cost Path function in ArcGIS 10.6.1, we modeled three types of spatial distances, namely Euclidean distance (for examining the effects of isolation by distance (IBD)), least cost distance considering roads (and highways) as agents of resistance (not considering the potential effects of viaducts), and least cost distance considering roads as agents of resistance and with viaducts facilitating movement. The resistance value for non-road areas and viaducts was set at 0 to represent low-cost distance, and for roads and water bodies set at 1 to represent the high cost of crossing roads and water bodies in the least-cost distance calculations. We calculated the correlations between genetic and spatial distances for each species at each treatment grid using Mantel tests with 99 permutations (default value) and Monte-Carlo corrected p values, using the R package ade4 [57]. One-tailed tests were used as the IBD model predicts that genetic differentiation will be positively correlated with increasing geographic distance [58]. All statistical analyses were conducted in R version 4.0.2 (R Core Team, 2020). Significance was considered at the α = 0.05 level.

3. Results

We obtained ~9000 SNPs using STACKS for the complete data set consisting of all four species and retained ~5000 SNPs, after removing SNPs under selection and disequilibrium (Table S2). For each species, PCAs across different datasets were generally consistent, therefore, we chose the SNP dataset filtering 10% missingness and absolute linkage (>0.95) for all subsequent analyses.
PCA did not reveal any clear patterns of clustering in genetic variation amongst viaduct, non-viaduct and control sites (Figure 4, Figures S2 and S3). However, PCA revealed slight genomic differentiation in M. whiteheadi between study areas along principal component 1 (Figure 5). All four species exhibited low genetic differentiation (FST) between the two study areas (M. whiteheadi FST = 0.022113378, M. rajah FST = 0.055480442, N. cremoriventer FST = 0.069197229, T. glis FST = 0.015547675). We found related individuals (relatedness > 0.1) between the two study areas in M. rajah and T. glis, suggesting occasional exchange of individuals between the two study areas.
Kruskal–Wallis tests revealed that there were differences in the ranking of treatment types by genetic distance for M. whiteheadi (χ22,79 = 28.22, p < 0.0001, Figure 6a) and M. rajah (χ22,87 = 39.01, p < 0.0001, Figure 6b). Genetic distances between individuals were lowest at viaduct sites for M. whiteheadi. For M. rajah, genetic distances were highest at viaduct sites and lowest at non-viaduct sites, with control sites in between. There were no significant differences in genetic distances between treatment types for N. cremoriventer (χ22,27 = 5.38, p = 0.068, Figure 6c) and T. glis (χ22,101 = 3.39, p = 0.18, Figure 6d).
Study area was also a significant factor influencing genetic distances between individuals within a site for M. whiteheadi (χ21,79 = 7.18, p < 0.01, Figure 7a), with individuals from Sungai Yu showing more genetic differentiation than at Kenyir. Likewise, T. glis also showed significantly more genetic differentiation at Sungai Yu (χ21,101 = 37.35, p < 0.0001, Figure 7d). M. rajah (χ21,87 = 0.72, p = 0.40, Figure 7b) did not show differences in genetic distances between study areas. Due to low sample sizes at Sungai Yu, only N. cremoriventer samples from Kenyir were included in the analysis (Figure 7c).
The dispersal of M. whiteheadi was found to be most consistent with isolation by distance (IBD) (highest correlation with Euclidean distance), indicating genetic and spatial distances between individuals were significantly positively correlated in both Kenyir and Sungai Yu, with minimal barrier effects of roads or facilitation effects of viaducts (Table 1). In Sungai Yu, the dispersal of M. rajah and T. glis correlated best with least cost distance, considering roads as resistance and viaducts as corridors facilitating movement, showing that they apparently prefer using viaducts to cross road barriers. No significant correlation with spatial distances were found in Kenyir for M. rajah, T. glis and N. cremoriventer. The paired samples of N. cremoriventer in Sungai Yu were insufficient for comparisons (Table 1).

4. Discussion

Our study showed that eco-viaducts may facilitate movement and the maintenance of gene flow in a landscape bisected by roads at least in some species such as M. whiteheadi. Presumably, the higher genetic similarity of individuals at viaduct sites in these species, relative to that at non-viaduct sites, is due to movement of individuals through the viaducts [59], thus reducing population subdivision caused by habitat fragmentation.
Highways did not seem to pose a barrier to gene flow in M. rajah (Figure 6b). Surprisingly, the greatest genetic differentiation amongst the three treatment types for M. rajah was at viaduct sites. All spatial distance model correlations were significant as well, indicating that viaducts were not essential for M. rajah to cross roads. This differential impact of roads is likely because for a less agile ground dwelling small mammal such as M. rajah, the effort taken to cross natural structures in forests such as dense vegetation and uneven ground, or a longer route to reach a viaduct and go through it, could be similar to crossing artificial structures such as roads (Figure 6b). Perhaps M. rajah only uses the viaducts for dispersal to other more favourable habitats, as viaducts may also facilitate the movements of its predators such as civets and leopard cats [60].
Movement and gene flow in M. whiteheadi were influenced by road barriers as well as the study area. The M. whiteheadi population in Kenyir was slightly differentiated from the population in Sungai Yu. Their dispersal was also most aligned with the Euclidean distance model in both study areas, pointing to IBD. A distinction between populations was not observed in the other three species. A lack of differentiation between the two study areas for the other three species is not surprising as the study areas are relatively near, there are no severe breaks in forest cover at a landscape scale, and genetic differentiation is not only the result of isolation by distance. This is why we chose to use individual-based spatial analyses to investigate the subtle differentiation.
Greater genetic connectivity was observed for both M. whiteheadi and T. glis at Kenyir than at Sungai Yu. The highway at Sungai Yu is more than double the width of the highway at Kenyir and it is probably more of a barrier to movement than the narrower Kenyir highway. Thus, the significance of viaducts in facilitating movement and increasing genetic connectivity across wider roads such as at Sungai Yu would be greater than at Kenyir. This was probably the case for M. whiteheadi as it also showed the lowest genetic differentiation in viaduct sites. Captures of M. whiteheadi individuals within the viaducts themselves suggest that M. whiteheadi does use viaducts to cross roads and may even have part of their home ranges overlapping the viaducts. Of interest was the observation that for T. glis, individuals at viaduct sites were still genetically more dissimilar than those at control sites. This shows that despite the addition of viaducts to connect fragmented populations, they will not be able to restore genetic connectivity to match that of the original intact forest.
In Kenyir, the dispersal of M. rajah, N. cremoriventer and T. glis did not correlate significantly with any of the least cost distance models (Table 1), suggesting that the dispersal ability of these three species goes beyond the spatial scale of the study in this area. The highway bisecting the forest in Kenyir is much narrower than the highway in Sungai Yu, likely presenting itself as less of a barrier to movement. M. rajah has been found as road kills in Kenyir, evidence that it does cross roads directly if narrow enough [59]. T. glis is a relatively mobile species, able to disperse >4 km [61]; one individual in this study was trapped in two viaducts, a distance separation of >5 km.
There was no clear pattern of clustering in genetic variation amongst viaduct, non-viaduct and control sites in intact forest (Figure 4, Figures S2 and S3), suggesting that there are no significant differences in genetic composition of small mammal species amongst the three habitats. However, there was slight genetic differentiation between the M. whiteheadi populations at Kenyir and Sungai Yu (Figure 5). As the smallest species in this study, it is expected to have the least dispersal ability and may be starting to show genetic divergence between populations separated by greater distances.
The upgrading of the highways and construction of the viaducts in this study were relatively recent; therefore, population genetic data may not fully reflect the evolutionary change if the generation time is long. It is possible that not enough time has passed to allow the vegetation in the viaducts to grow and establish stable corridors or habitats. Mills and Allendorf [62] suggested that only one migrant per generation is sufficient to maintain genetic diversity while allowing some divergence between populations. Rosenberg et al. [24] also argue that corridors are more effective at maintaining movement between populations for habitat specialists rather than for habitat generalists. All four species in this study are tolerant to disturbed forests and forest edge, and their diets are varied, consisting of insects, fruit and seeds. This may explain why the viaducts did not increase genetic connectivity as much as expected.
Our study shows that viaducts can increase genetic connectivity, but are not effective across all species and roads or highways. The effectiveness of the viaducts in reducing the barrier effects of roads depends on the target species’ dispersal abilities, generation time, perception of landscape characteristics, road width, viaduct age and vegetation structural maturity. Viaducts work better when they have a structure which is more conducive for movement than a direct crossing across the road. However, from this study we observed that viaducts were not able to restore genetic connectivity to match that of the original intact forest and should therefore only be used as a last resort if it is not possible to keep the forest intact.

5. Conclusions

Our study shows that wildlife crossings such as viaducts may assist in maintaining gene flow within populations of certain small mammals in areas with relatively wider roads. In a managed landscape, there must be a balance between the economic benefits of building roads and the ecological benefits to species in maintaining connectivity. We have shown that one of the benefits of constructing and maintaining viaducts is an increase in gene flow within populations in some species, which may result in an increase in population persistence and a decrease in inbreeding.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani14030426/s1, Figure S1: Photos of trap types used in this study; Figure S2: Principal component analysis of genetic differentiation among sampled individuals (PC3 and PC4); Figure S3: Principal component analysis of genetic differentiation among sampled individuals (PC5 and PC6); Table S1: Individual species samples and sites at which they were captured; Table S2: Number of SNPs at each stage of filtering.

Author Contributions

Conceptualization, T.C.Y.H., Q.T. and F.E.R.; methodology, T.C.Y.H., Q.T., E.Y.X.N. and F.E.R.; validation, Q.T., E.Y.X.N., E.M.S. and F.E.R.; formal analysis, T.C.Y.H., Q.T., E.Y.X.N. and F.E.R.; investigation, T.C.Y.H., Q.T., E.Y.X.N., J.L.C. and F.E.R.; resources, F.E.R.; data curation, T.C.Y.H. and J.L.C.; writing—original draft preparation, T.C.Y.H.; writing—review and editing, Q.T., E.Y.X.N., J.L.C., E.M.S. and F.E.R.; visualization, T.C.Y.H., Q.T., E.Y.X.N. and F.E.R.; supervision, E.M.S. and F.E.R.; project administration, J.L.C., Q.T. and F.E.R.; funding acquisition, F.E.R. All authors have read and agreed to the published version of the manuscript.

Funding

Our research was financially supported by the Rufford Foundation Grant for Nature Conservation, grant number 21698-1, an IDEA WILD Equipment Grant, Skyrail Rainforest Foundation Research Grant and Animal Behavior Society Student Research Grant. Labwork and analytical work for this project was funded by a Singapore Ministry of Education (MOE) Tier 2 grant to FER (WBS number R-154-000-C41-112).

Institutional Review Board Statement

We would like to thank the Economic Planning Unit (Permit number UPE 40/200/19/3363) of the Prime Minister’s Department and the Department of Wildlife and National Parks (HQ, Terengganu and Pahang) (Permit numbers JPHL&TN(IP):100-34/1.24 JLD.8, JPHL&TN(IP):100-34/1.24 JLD.12, JPHL&TN(IP):90-7/8 JLD.16 and NRE 600-2/2/21 JLD.6) for granting permission to undertake this project, as well as the Forestry Department of Peninsular Malaysia (Terengganu and Pahang) (Permit number JH/100 JLD.16) for granting permission to enter the forest reserves in Kenyir and Sungai Yu.

Informed Consent Statement

Not applicable.

Data Availability Statement

All available data are presented in this manuscript.

Acknowledgments

We would like to thank the staff and students of Universiti Malaysia Terengganu, Universiti Sains Malaysia, Department of Wildlife and National Parks (PERHILITAN), Rimba, Ecoteer and MYCAT staff and volunteers, Sohud bin Yaacob and Dome Nikong for providing the necessary facilities and assistance. We would also like to thank Tan Poai Ean and the late Dr Lim Boo Liat for assistance with identification of the specimens. We also acknowledge all respondents for their cooperation and knowledge sharing in this study as well as local guides that assisted us throughout this research.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Fahrig, L.; Rytwinski, T. Effects of roads on animal abundance: An empirical review and synthesis. Ecol. Soc. 2009, 14, 1. [Google Scholar] [CrossRef]
  2. Debinski, D.M.; Holt, R.D. A survey and overview of habitat fragmentation experiments. Conserv. Biol. 2000, 14, 342–355. [Google Scholar] [CrossRef]
  3. Laurance, W.F. Contemporary drivers of habitat fragmentation. In Global Forest Fragmentation; CABI: Wallingford, UK, 2014; pp. 20–27. [Google Scholar]
  4. Robert McDonald, W.; Cassady St. Clair, C. The effects of artificial and natural barriers on the movement of small mammals in Banff National Park, Canada. Oikos 2004, 105, 397–407. [Google Scholar] [CrossRef]
  5. Cros, E.; Ng, E.Y.; Oh, R.R.; Tang, Q.; Benedick, S.; Edwards, D.P.; Tomassi, S.; Irestedt, M.; Ericson, P.G.; Rheindt, F.E. Fine-scale barriers to connectivity across a fragmented South-East Asian landscape in six songbird species. Evol. Appl. 2020, 13, 1026–1036. [Google Scholar] [CrossRef] [PubMed]
  6. Rivera-Ortíz, F.; Aguilar, R.; Arizmendi, M.; Quesada, M.; Oyama, K. Habitat fragmentation and genetic variability of tetrapod populations. Anim. Conserv. 2015, 18, 249–258. [Google Scholar] [CrossRef]
  7. Chattopadhyay, B.; Garg, K.M.; Mendenhall, I.H.; Rheindt, F.E. Historic DNA reveals Anthropocene threat to a tropical urban fruit bat. Curr. Biol. 2019, 29, R1299–R1300. [Google Scholar] [CrossRef] [PubMed]
  8. Ewers, R.M.; Didham, R.K. Confounding factors in the detection of species responses to habitat fragmentation. Biol. Rev. 2006, 81, 117–142. [Google Scholar] [CrossRef]
  9. Christie, M.R.; Knowles, L.L. Habitat corridors facilitate genetic resilience irrespective of species dispersal abilities or population sizes. Evol. Appl. 2015, 8, 454–463. [Google Scholar] [CrossRef] [PubMed]
  10. Haddad, N.M.; Bowne, D.R.; Cunningham, A.; Danielson, B.J.; Levey, D.J.; Sargent, S.; Spira, T. Corridor use by diverse taxa. Ecology 2003, 84, 609–615. [Google Scholar] [CrossRef]
  11. Noss, R.F. Corridors in real landscapes: A reply to Simberloff and Cox. Conserv. Biol. 1987, 1, 159–164. [Google Scholar] [CrossRef]
  12. Bennett, A.F. Habitat corridors and the conservation of small mammals in a fragmented forest environment. Landsc. Ecol. 1990, 4, 109–122. [Google Scholar] [CrossRef]
  13. Saunders, D.A.; Hobbs, R.J.; Margules, C.R. Biological consequences of ecosystem fragmentation: A review. Conserv. Biol. 1991, 5, 18–32. [Google Scholar] [CrossRef]
  14. Mabry, K.E.; Barrett, G.W. Effects of corridors on home range sizes and interpatch movements of three small mammal species. Landsc. Ecol. 2002, 17, 629–636. [Google Scholar] [CrossRef]
  15. Clevenger, A.P.; Chruszcz, B.; Gunson, K. Drainage culverts as habitat linkages and factors affecting passage by mammals. J. Appl. Ecol. 2001, 38, 1340–1349. [Google Scholar] [CrossRef]
  16. McDonald, W.; St Clair, C.C. Elements that promote highway crossing structure use by small mammals in Banff National Park. J. Appl. Ecol. 2004, 41, 82–93. [Google Scholar] [CrossRef]
  17. Bellis, M.A.; Griffin, C.R.; Warren, P.; Jackson, S.D. Utilizing a multi-technique, multi–taxa approach to monitoring wildlife passageways in southern vermont. Oecologia Aust. 2013, 17, 111–128. [Google Scholar] [CrossRef]
  18. Beier, P.; Noss, R.F. Do habitat corridors provide connectivity? Conserv. Biol. 1998, 12, 1241–1252. [Google Scholar] [CrossRef]
  19. Haddad, N.M. Finding the corridor more traveled. Proc. Natl. Acad. Sci. USA 2008, 105, 19569–19570. [Google Scholar] [CrossRef]
  20. Gonzalez, A.; Lawton, J.; Gilbert, F.; Blackburn, T.; Evans-Freke, I. Metapopulation dynamics, abundance, and distribution in a microecosystem. Science 1998, 281, 2045–2047. [Google Scholar] [CrossRef]
  21. Dunning, J.B.; Borgella, R.; Clements, K.; Meffe, G.K. Patch isolation, corridor effects, and colonization by a resident sparrow in a managed pine woodland. Conserv. Biol. 1995, 9, 542–550. [Google Scholar] [CrossRef]
  22. Fahrig, L.; Merriam, G. Habitat patch connectivity and population survival. Ecology 1985, 66, 1762–1768. [Google Scholar] [CrossRef]
  23. Haddad, N.M.; Baum, K.A. An experimental test of corridor effects on butterfly densities. Ecol. Appl. 1999, 9, 623–633. [Google Scholar] [CrossRef]
  24. Rosenberg, D.K.; Noon, B.R.; Meslow, E.C. Biological corridors: Form, function, and efficacy. BioScience 1997, 47, 677–687. [Google Scholar] [CrossRef]
  25. Haddad, N.M. Corridor and distance effects on interpatch movements: A landscape experiment with butterflies. Ecol. Appl. 1999, 9, 612–622. [Google Scholar] [CrossRef]
  26. Tewksbury, J.J.; Levey, D.J.; Haddad, N.M.; Sargent, S.; Orrock, J.L.; Weldon, A.; Danielson, B.J.; Brinkerhoff, J.; Damschen, E.I.; Townsend, P. Corridors affect plants, animals, and their interactions in fragmented landscapes. Proc. Natl. Acad. Sci. USA 2002, 99, 12923–12926. [Google Scholar] [CrossRef] [PubMed]
  27. Haas, C.A. Dispersal and use of corridors by birds in wooded patches on an agricultural landscape. Conserv. Biol. 1995, 9, 845–854. [Google Scholar] [CrossRef]
  28. Mech, S.G.; Hallett, J.G. Evaluating the effectiveness of corridors: A genetic approach. Conserv. Biol. 2001, 15, 467–474. [Google Scholar] [CrossRef]
  29. Sutcliffe, O.L.; Thomas, C.D. Open corridors appear to facilitate dispersal by ringlet butterflies (Aphantopus hyperantus) between woodland clearings. Conserv. Biol. 1996, 10, 1359–1365. [Google Scholar] [CrossRef]
  30. Beier, P. Dispersal of juvenile cougars in fragmented habitat. J. Wildl. Manag. 1995, 59, 228–237. [Google Scholar] [CrossRef]
  31. Aars, J.; Ims, R.A. The effect of habitat corridors on rates of transfer and interbreeding between vole demes. Ecology 1999, 80, 1648–1655. [Google Scholar] [CrossRef]
  32. Hale, M.L.; Lurz, P.W.; Shirley, M.D.; Rushton, S.; Fuller, R.M.; Wolff, K. Impact of landscape management on the genetic structure of red squirrel populations. Science 2001, 293, 2246–2248. [Google Scholar] [CrossRef]
  33. Orrock, J.L. Conservation corridors affect the fixation of novel alleles. Conserv. Genet. 2005, 6, 623–630. [Google Scholar] [CrossRef]
  34. Orrock, J.L.; Damschen, E.I. Corridors cause differential seed predation. Ecol. Appl. 2005, 15, 793–798. [Google Scholar] [CrossRef]
  35. Weldon, A.J.; Haddad, N.M. The effects of patch shape on Indigo Buntings: Evidence for an ecological trap. Ecology 2005, 86, 1422–1431. [Google Scholar] [CrossRef]
  36. Kjeldsen, S.R.; Zenger, K.R.; Leigh, K.; Ellis, W.; Tobey, J.; Phalen, D.; Melzer, A.; FitzGibbon, S.; Raadsma, H.W. Genome-wide SNP loci reveal novel insights into koala (Phascolarctos cinereus) population variability across its range. Conserv. Genet. 2016, 17, 337–353. [Google Scholar] [CrossRef]
  37. Szulkin, M.; Gagnaire, P.A.; Bierne, N.; Charmantier, A. Population genomic footprints of fine-scale differentiation between habitats in Mediterranean blue tits. Mol. Ecol. 2016, 25, 542–558. [Google Scholar] [CrossRef]
  38. Angeloni, F.; Wagemaker, N.; Vergeer, P.; Ouborg, J. Genomic toolboxes for conservation biologists. Evol. Appl. 2012, 5, 130–143. [Google Scholar] [CrossRef] [PubMed]
  39. Sawaya, M.A.; Kalinowski, S.T.; Clevenger, A.P. Genetic connectivity for two bear species at wildlife crossing structures in Banff National Park. Proc. R. Soc. B Biol. Sci. 2014, 281, 20131705. [Google Scholar] [CrossRef]
  40. Braunisch, V.; Segelbacher, G.; Hirzel, A.H. Modelling functional landscape connectivity from genetic population structure: A new spatially explicit approach. Mol. Ecol. 2010, 19, 3664–3678. [Google Scholar] [CrossRef]
  41. Burkart, S.; Gugerli, F.; Senn, J.; Kuehn, R.; Bolliger, J. Evaluating the functionality of expert-assessed wildlife corridors with genetic data from roe deer. Basic Appl. Ecol. 2016, 17, 52–60. [Google Scholar] [CrossRef]
  42. Brunke, J.; Russo, I.-R.M.; Orozco-terWengel, P.; Zimmermann, E.; Bruford, M.W.; Goossens, B.; Radespiel, U. Dispersal and genetic structure in a tropical small mammal, the Bornean tree shrew (Tupaia longipes), in a fragmented landscape along the Kinabatangan River, Sabah, Malaysia. BMC Genet. 2020, 21, 1–13. [Google Scholar] [CrossRef]
  43. Steele, G.M.; Randolph, S.E. An experimental evaluation of conventional control measures against the sheep tick, Ixodes ricinus (L.)(Acari: Ixodidae). I. A unimodal seasonal activity pattern. Bull. Entomol. Res. 1985, 75, 489–500. [Google Scholar] [CrossRef]
  44. Hecht, L.; Allcock, M. Potential Effects of Habitat Fragmentation on Wild Animal Welfare. 2020. Available online: https://www.researchgate.net/publication/342292473_Potential_effects_of_habitat_fragmentation_on_wild_animal_welfare (accessed on 13 November 2023).
  45. Francis, C.M.; Barrett, P. A Guide to the Mammals of Southeast Asia; New Holland Publishers (UK) Ltd.: London, UK, 2008. [Google Scholar]
  46. Hoffmann, A.; Decher, J.; Rovero, F.; Schaer, J.; Voigt, C.; Wibbelt, G. Field methods and techniques for monitoring mammals. Man. Field Rec. Tech. Protoc. All Taxa Biodivers. Invent. 2010, 8, 482–529. [Google Scholar]
  47. Ng, E.Y.; Li, S.; Zhang, D.; Garg, K.M.; Song, G.; Martinez, J.; Hung, L.M.; Tu, V.T.; Fuchs, J.; Dong, L. Genome-wide SNPs confirm plumage polymorphism and hybridisation within a Cyornis flycatcher species complex. Zool. Scr. 2023, 52, 1–16. [Google Scholar] [CrossRef]
  48. Catchen, J.; Hohenlohe, P.A.; Bassham, S.; Amores, A.; Cresko, W.A. Stacks: An analysis tool set for population genomics. Mol. Ecol. 2013, 22, 3124–3140. [Google Scholar] [CrossRef]
  49. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef]
  50. Lerner, H.R.; Meyer, M.; James, H.F.; Hofreiter, M.; Fleischer, R.C. Multilocus resolution of phylogeny and timescale in the extant adaptive radiation of Hawaiian honeycreepers. Curr. Biol. 2011, 21, 1838–1844. [Google Scholar] [CrossRef]
  51. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; Subgroup, G.P.D.P. The sequence alignment/map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef]
  52. Foll, M.; Gaggiotti, O. A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: A Bayesian perspective. Genetics 2008, 180, 977–993. [Google Scholar] [CrossRef] [PubMed]
  53. Chang, C.C.; Chow, C.C.; Tellier, L.C.; Vattikuti, S.; Purcell, S.M.; Lee, J.J. Second-generation PLINK: Rising to the challenge of larger and richer datasets. Gigascience 2015, 4, s13742-13015-10047-13748. [Google Scholar] [CrossRef] [PubMed]
  54. Zheng, X. SNPRelate: Parrallel computing toolset for genome-wide association studies. R Package Version 2012, 95, B9. [Google Scholar]
  55. Danecek, P.; Auton, A.; Abecasis, G.; Albers, C.A.; Banks, E.; DePristo, M.A.; Handsaker, R.E.; Lunter, G.; Marth, G.T.; Sherry, S.T. The variant call format and VCFtools. Bioinformatics 2011, 27, 2156–2158. [Google Scholar] [CrossRef]
  56. Kamvar, Z.N.; Tabima, J.F.; Grünwald, N.J. Poppr: An R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ 2014, 2, e281. [Google Scholar] [CrossRef]
  57. Bougeard, S.; Dray, S. Supervised multiblock analysis in R with the ade4 package. J. Stat. Softw. 2018, 86, 1–17. [Google Scholar] [CrossRef]
  58. Wright, S. Isolation by distance. Genetics 1943, 28, 114. [Google Scholar] [CrossRef] [PubMed]
  59. Hui, T.C.; Slade, E.M.; Chong, J.L. Roadkills in Northern Peninsular Malaysia. Front. Environ. Sci. 2021, 9, 637462. [Google Scholar] [CrossRef]
  60. Hedges, L.; Clements, G.; Aziz, S.; Yap, W.; Laurance, S.; Goosem, M.; Laurance, W. Small carnivore records from a threatened habitat linkage in Terengganu, Peninsular Malaysia. Small Carniv. Conserv. 2013, 49, 9–14. [Google Scholar]
  61. Mariana, A.; Shukor, M.; Muhd, N.H.; Intan, N.B.; Ho, T. Movements and home range of a common species of tree–shrew, Tupaia glis, surrounding houses of otoacariasis cases in Kuantan, Pahang, Malaysia. Asian Pac. J. Trop. Med. 2010, 3, 427–434. [Google Scholar] [CrossRef]
  62. Mills, L.S.; Allendorf, F.W. The one-migrant-per-generation rule in conservation and management. Conserv. Biol. 1996, 10, 1509–1518. [Google Scholar] [CrossRef]
Figure 1. Sampling sites in (a) Kenyir, (b) Sungai Yu. Sites situated away from highways, along highways and adjacent to eco-viaducts are named “Control”, “Non-viaduct” and “Viaduct”, respectively. (c) Locations of study areas in Peninsular Malaysia (yellow: Kenyir; pink: Sungai Yu). Pictures of eco-viaducts (wildlife crossings) at (d) Kenyir and (e) Sungai Yu (photos by Tabitha Hui).
Figure 1. Sampling sites in (a) Kenyir, (b) Sungai Yu. Sites situated away from highways, along highways and adjacent to eco-viaducts are named “Control”, “Non-viaduct” and “Viaduct”, respectively. (c) Locations of study areas in Peninsular Malaysia (yellow: Kenyir; pink: Sungai Yu). Pictures of eco-viaducts (wildlife crossings) at (d) Kenyir and (e) Sungai Yu (photos by Tabitha Hui).
Animals 14 00426 g001
Figure 2. Schematic diagram of sampling design. The expanded grid box shows the layout of traps. Blue boxes: viaduct sites, red boxes: non-viaduct sites, green boxes: control sites.
Figure 2. Schematic diagram of sampling design. The expanded grid box shows the layout of traps. Blue boxes: viaduct sites, red boxes: non-viaduct sites, green boxes: control sites.
Animals 14 00426 g002
Figure 3. The four study species: (a) Maxomys whiteheadi, (b) Maxomys rajah, (c) Niviventer cremoriventer and (d) Tupaia glis. Photos by Tabitha Hui.
Figure 3. The four study species: (a) Maxomys whiteheadi, (b) Maxomys rajah, (c) Niviventer cremoriventer and (d) Tupaia glis. Photos by Tabitha Hui.
Animals 14 00426 g003
Figure 4. Principal component analysis (axes PC1 and PC2) of genetic differentiation among sampled individuals. (a) Maxomys whiteheadi, (b) Maxomys rajah, (c) Niviventer cremoriventer and (d) Tupaia glis. Blue: viaduct sites, red: non-viaduct sites, green: control sites. Same shape: same pair, filled/open shapes: opposite sides of the road of the same pair.
Figure 4. Principal component analysis (axes PC1 and PC2) of genetic differentiation among sampled individuals. (a) Maxomys whiteheadi, (b) Maxomys rajah, (c) Niviventer cremoriventer and (d) Tupaia glis. Blue: viaduct sites, red: non-viaduct sites, green: control sites. Same shape: same pair, filled/open shapes: opposite sides of the road of the same pair.
Animals 14 00426 g004
Figure 5. Principal component analysis (axes PC1 and PC2) of genetic differentiation among sampled individuals of Maxomys whiteheadi showing slight separation between Sungai Yu and Kenyir populations.
Figure 5. Principal component analysis (axes PC1 and PC2) of genetic differentiation among sampled individuals of Maxomys whiteheadi showing slight separation between Sungai Yu and Kenyir populations.
Animals 14 00426 g005
Figure 6. Comparisons of genetic distances between non-viaduct, viaduct and control sites. (a) Maxomys whiteheadi, (b) Maxomys rajah, (c) Niviventer cremoriventer and (d) Tupaia glis. The boxes, lines in the middle, whiskers and dots represent the interquartile range, median, minimum and maximum values (no more than 1.5 times the interquartile range) and outliers, respectively.
Figure 6. Comparisons of genetic distances between non-viaduct, viaduct and control sites. (a) Maxomys whiteheadi, (b) Maxomys rajah, (c) Niviventer cremoriventer and (d) Tupaia glis. The boxes, lines in the middle, whiskers and dots represent the interquartile range, median, minimum and maximum values (no more than 1.5 times the interquartile range) and outliers, respectively.
Animals 14 00426 g006
Figure 7. Comparisons of genetic distances between study areas. (a) Maxomys whiteheadi, (b) Maxomys rajah, (c) Niviventer cremoriventer and (d) Tupaia glis. The boxes, lines in the middle, whiskers and dots represent the interquartile range, median, minimum and maximum values (no more than 1.5 times the interquartile range) and outliers, respectively.
Figure 7. Comparisons of genetic distances between study areas. (a) Maxomys whiteheadi, (b) Maxomys rajah, (c) Niviventer cremoriventer and (d) Tupaia glis. The boxes, lines in the middle, whiskers and dots represent the interquartile range, median, minimum and maximum values (no more than 1.5 times the interquartile range) and outliers, respectively.
Animals 14 00426 g007
Table 1. Mantel’s r correlations between the genetic distances of each species and each spatial distance model (Euclidean: Euclidean distance; Road: least cost distance considering roads; Viaduct: least cost distance considering roads as resistance and viaducts facilitating movement). Results in bold are the spatial distance models with the highest correlation for the species in each study area, where multiple spatial distance models were significant.
Table 1. Mantel’s r correlations between the genetic distances of each species and each spatial distance model (Euclidean: Euclidean distance; Road: least cost distance considering roads; Viaduct: least cost distance considering roads as resistance and viaducts facilitating movement). Results in bold are the spatial distance models with the highest correlation for the species in each study area, where multiple spatial distance models were significant.
KenyirSungai Yu
Spatial Distance ModelMantel r Simulated p-ValueMantel rSimulated p-Value
Euclidean0.1180.150.3550.01
Maxomys
rajah
Road0.1170.160.3730.01
Viaduct0.1170.140.5060.01
Euclidean0.2450.010.0820.02
Maxomys whiteheadiRoad0.2370.010.1240.19
Viaduct0.2370.010.1260.22
Euclidean0.0960.1
Niviventer cremoriventerRoad0.1120.12
Viaduct0.1170.1
Euclidean0.1210.170.1130.17
Tupaia glisRoad0.1020.190.2040.09
Viaduct0.1010.230.2070.05
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hui, T.C.Y.; Tang, Q.; Ng, E.Y.X.; Chong, J.L.; Slade, E.M.; Rheindt, F.E. Small-Mammal Genomics Highlights Viaducts as Potential Dispersal Conduits for Fragmented Populations. Animals 2024, 14, 426. https://doi.org/10.3390/ani14030426

AMA Style

Hui TCY, Tang Q, Ng EYX, Chong JL, Slade EM, Rheindt FE. Small-Mammal Genomics Highlights Viaducts as Potential Dispersal Conduits for Fragmented Populations. Animals. 2024; 14(3):426. https://doi.org/10.3390/ani14030426

Chicago/Turabian Style

Hui, Tabitha C. Y., Qian Tang, Elize Y. X. Ng, Ju Lian Chong, Eleanor M. Slade, and Frank E. Rheindt. 2024. "Small-Mammal Genomics Highlights Viaducts as Potential Dispersal Conduits for Fragmented Populations" Animals 14, no. 3: 426. https://doi.org/10.3390/ani14030426

APA Style

Hui, T. C. Y., Tang, Q., Ng, E. Y. X., Chong, J. L., Slade, E. M., & Rheindt, F. E. (2024). Small-Mammal Genomics Highlights Viaducts as Potential Dispersal Conduits for Fragmented Populations. Animals, 14(3), 426. https://doi.org/10.3390/ani14030426

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop