The influence of riverine barriers on phylogeographic patterns of Malagasy reed frogs (Heterixalus)

Genetic structure and landscape effects on gene flow in the Neotropical lizard Norops brasiliensis (Squamata: Dactyloidae)

One key research goal of evolutionary biology is to understand the origin and maintenance of genetic variation. In the Cerrado, the South American savanna located primarily in the Central Brazilian Plateau, many hypotheses have been proposed to explain how landscape features (e.g., geographic distance, river barriers, topographic compartmentalization, and historical climatic fluctuations) have promoted genetic structure by mediating gene flow. Here, we asked whether these landscape features have influenced the genetic structure and differentiation in the lizard species Norops brasiliensis (Squamata: Dactyloidae). To achieve our goal, we used a genetic clustering analysis and estimate an effective migration surface to assess genetic structure in the focal species. Optimized isolation-by-resistance models and a simulation-based approach combined with machine learning (convolutional neural network; CNN) were then used to infer current and historical effects on population genetic structure through 12 unique landscape models. We recovered five geographically distributed populations that are separated by regions of lower-than-expected gene flow. The results of the CNN showed that geographic distance is the sole predictor of genetic variation in N. brasiliensis, and that slope, rivers, and historical climate had no discernible influence on gene flow. Our novel CNN approach was accurate (89.5%) in differentiating each landscape model. CNN and other machine learning approaches are still largely unexplored in landscape genetics studies, representing promising avenues for future research with increasingly accessible genomic datasets.

Phylogenetic relationships and genetic differentiation of two Salamandrella species as revealed via COI gene from Northeastern China

Due to traditional classification methods' limitations, some cryptic species remain undiscovered. To better explore the existence of the Schrenck salamander (Salamandrella tridactyla, a cryptic species of Siberian salamander S. keyserlingii) in China, we conducted a molecular phylogenetic analysis to confirm the taxonomic relationship among Salamandrella species and investigate genetic variation. We used complete sequences of the mitochondrial COI gene from 65 specimens collected across a wide range in Northeastern China. Thirty-five haplotypes were obtained from six populations. They showed medium-high haplotype diversity (Hd) and low nucleotide polymorphism (π). The phylogenetic tree and haplotype network analysis revealed that populations from Greater Khingan Ridge (Huma: HM) and Lesser Khingan Ridge (Tieli: TL) belong to S. keyserlingii, while populations from Changbai Mountain (Shangzhi-zhuziying: SZ, Shangzhi-cuijia: SC, Hailin: HL, and Baishan: BS) belong to S. tridactyla. This indicates the monophyly of Salamandrella and each of the two species. There was a substantial level of genetic differentiation between different species and within populations of the same species. This differentiation was significantly related to geographical distance. At last, the mismatch distribution and neutrality analyses indicated that the TL populations have undergone expansion of history. The study supplements the distributional range of Schrenck salamander. And it provides a theoretical basis for species conservation of Salamandrella species.

Diversification processes in Gerp's mouse lemur demonstrate the importance of rivers and altitude as biogeographic barriers in Madagascar's humid rainforests

Madagascar exhibits exceptionally high levels of biodiversity and endemism. Models to explain the diversification and distribution of species in Madagascar stress the importance of historical variability in climate conditions which may have led to the formation of geographic barriers by changing water and habitat availability. The relative importance of these models for the diversification of the various forest-adapted taxa of Madagascar has yet to be understood. Here, we reconstructed the phylogeographic history of Gerp's mouse lemur (Microcebus gerpi) to identify relevant mechanisms and drivers of diversification in Madagascar's humid rainforests. We used restriction site associated DNA (RAD) markers and applied population genomic and coalescent-based techniques to estimate genetic diversity, population structure, gene flow and divergence times among M. gerpi populations and its two sister species M. jollyae and M. marohita. Genomic results were complemented with ecological niche models to better understand the relative barrier function of rivers and altitude. We show that M. gerpi diversified during the late Pleistocene. The inferred ecological niche, patterns of gene flow and genetic differentiation in M. gerpi suggest that the potential for rivers to act as biogeographic barriers depended on both size and elevation of headwaters. Populations on opposite sides of the largest river in the area with headwaters that extend far into the highlands show particularly high genetic differentiation, whereas rivers with lower elevation headwaters have weaker barrier functions, indicated by higher migration rates and admixture. We conclude that M. gerpi likely diversified through repeated cycles of dispersal punctuated by isolation to refugia as a result of paleoclimatic fluctuations during the Pleistocene. We argue that this diversification scenario serves as a model of diversification for other rainforest taxa that are similarly limited by geographic factors. In addition, we highlight conservation implications for this critically endangered species, which faces extreme habitat loss and fragmentation.

Genetic Diversity Relationship in Azakheli Buffalo Inferred from mtDNA and MC1R Sequences Comparison

Azakheli is relatively smaller riverine breed with a very peculiar characteristics kept under unique traditional husbandry practices in comparison with rest of the Pakistani buffalo breeds; however, milk production is comparable. The present study was conducted to evaluate the genetic diversity of Azakheli breed. A total of sixty-six blood samples were collected for the amplification of mtDNA D-loop region and MC1R gene sequencing analysis. Median-joining network analysis of 191 mtDNA D-loop sequences of Azakheli and eight Indian riverine buffalo breeds clustered into three clades. Ancient Azakheli Region 1 clade was the oldest with the highest mutation steps and was present close to the root of UPGMA phylogenetic tree. There was 5 mutated lines distance between Pakistan buffalo and Indian riverine buffaloes. The populations of neighboring countries did not share any haplotypes with Azakheli buffalo of Pakistan. Possibly, residing for so long in the cold atmosphere and high elevation regions caused the mutation in mtDNA D-loop, though these conditions did not affect the overall performance of Azakheli as milch buffalo breed of Pakistan. MC1R analyses showed high mutations in Azakheli of Albino phenotype and all the black phenotype individuals of Azakheli buffalo share haplotypes with dominant Chinese and Indian black phenotypes buffaloes in MC1R median-joining network, indicating the reason of black coat color is due to MC1R gene. The haplotype diversity and nucleotide diversity was (H. 0.923, Pi: 0.00895) in Azakheli. Current results illustrated Asian ancestry for Azakheli buffalo, and mtDNA and MC1R analyses provided further evidence. Additional genetic analyses and archeological studies may provide further insight into the domestication period and history of Azakheli buffalo breed. The further studies are required on different coat colors with different genes on Azakheli buffalo to understand the phenotype variation.

Genetic diversity in peridomiciliary populations of Triatoma mexicana (Hemiptera: Reduviidae: Triatominae) in central Mexico

Triatoma mexicana is an important vector of Trypanosoma cruzi-the etiological agent of Chagas disease. This triatomine species occurs in central Mexico, but little is known about its genetic variability. Using Cyt-b gene as a genetic marker, in this study, we determined the population genetic structure of T. mexicana collected from the States of Hidalgo, Guanajuato, and Queretaro where populations are largely peridomiciliary. A Bayesian approach was performed for the design of phylogenies, median-joining networks, and clustering among populations of T. mexicana. Our results show that the Hidalgo population was the most distinct, with the highest genetic and haplotypic variation (Hd = 0.963, π = 0.06129, and ɵ = 0.05469). Moderate gene flow (Nm) was determined among populations of Hidalgo and Queretaro. Populations from the three states showed differentiation (F_ST) values ranging from 0.22 to 0.3, suggesting an important genetic differentiation. The phylogenetic analysis showed the presence of five well-defined groups, as well as the haplotype network, where 24 haplotypes were observed forming five haplogroups with high mutational steps among them: 68 (Hgo-W2), 26 (Qto), 59 (Hgo-M), 44 (Hgo-W1), and 46 (Gto). Genetic isolation was apparently inferred in the Guanajuato population; however, the Mantel test did not show correlation between genetic (F_ST) and geographic (km) distances (p = 0.05). The STRUCTURE analyses showed seven genetic clusters and it was observed that a single cluster predominates in each sampled location. However, genetic admixture was detected in four localities. Our results show evidence that there are multiple species within the collected sampling area.

Phylogeographic inference of Sumatran ranids bearing gastromyzophorous tadpoles with regard to the Pleistocene drainage systems of Sundaland

Rivers are known to act as biogeographic barriers in several strictly terrestrial taxa, while possibly serving as conduits of dispersal for freshwater-tolerant or -dependent species. However, the influence of river systems on genetic diversity depends on taxa-specific life history traits as well as other geographic factors. In amphibians, several studies have demonstrated that river systems have only minor influence on their divergence. Here, we assess the role of the paleodrainage systems of the Sunda region (with a focus on the island of Sumatra) in shaping the evolutionary history of two genera of frogs (Sumaterana and Wijayarana) whose tadpoles are highly dependent on cascading stream habitats. Our phylogenetic results show no clear association between the genetic diversification patterns of both anurans genera and the existence of paleodrainage systems. Time-calibrated phylogenies and biogeographical models suggest that these frogs colonized Sumatra and diversified on the island before the occurrence of the Pleistocene drainage systems. Both genera demonstrate phylogenetic structuring along a north–south geographic axis, the temporal dynamics of which coincide with the geological chronology of proto Sumatran and -Javan volcanic islands. Our results also highlight the chronic underestimation of Sumatran biodiversity and call for more intense sampling efforts on the island.

Population diversification in the frog Mantidactylus bellyi on an isolated massif in northern Madagascar based on genetic, morphological, bioacoustic and ecological evidence

In the processes that give rise to new species, changes first occur at the population level. But with the continuous nature of the divergence process, change in biological properties delimiting the shift from “individuals of divergent populations” towards “individuals of distinct species”, as well as abiotic factors driving the change, remain largely ambivalent. Here we study diversification processes at the population level in a semi-aquatic frog, Mantidactylus (Brygoomantis) bellyi, across the diverse vegetation types of Montagne d’Ambre National Park (MANP), Madagascar. Genetic diversity was assessed with seven newly developed microsatellite markers as well as mitochondrial DNA sequences and concordance with patterns of ecological, morphological, and bioacoustic divergence evaluated. We found M. bellyi lacking mitochondrial differentiation within MANP, while microsatellite datasets partitioned them into three highly differentiated, geographically separated subpopulations (with indications for up to five subpopulations). The molecular grouping–primarily clustering individuals by geographic proximity–was coincident with differences in mean depth and width of waters, suggesting a possible role of fluvial characteristics in genetic exchange in this stream-breeding species. Genetic clustering not consistent with differences in call properties, except for dominant call frequencies under the two-subpopulations model. Morphological divergence was mostly consistent with the genetic clustering; subpopulations strongly differed by their snout-vent length, with individuals from high-elevation subpopulations smaller than those from populations below 1000 m above sea level. These results exemplify how mountains and environmental conditions might primarily shape genetic and morphological divergence in frog populations, without strongly affecting their calls.

Evidence for ecological processes driving speciation among endemic lizards of Madagascar

Although genetic patterns produced by population isolation during speciation are well documented, the biogeographic and ecological processes that trigger speciation remain poorly understood. Alternative hypotheses for the biogeography and ecology of speciation include geographic isolation combined with niche conservation (soft allopatry) or parapatric distribution on an environmental gradient with niche divergence (ecological speciation). Here, we use species' distributions, environmental data, and two null models (the Random Translation and Rotation and the Background Similarity Test) to test these alternative hypotheses among 28 sister pairs of microendemic lizards in Madagascar. Our results demonstrate strong bimodal peaks along a niche divergence-conservation spectrum, with at least 25 out of 28 sister pairs exhibiting either niche conservation or divergence, and the remaining pairs showing weak ecological signals. Yet despite these significant results, we do not find strong associations of niche conservation with allopatric distributions or niche divergence with parapatric distributions. Our findings thus provide strong evidence of a role for ecological processes driving speciation, rather than the classic expectation of speciation through geographic isolation, but demonstrate that the link between ecological speciation and parapatry is complex and requires further analysis of a broader taxonomic sample to fully resolve.

Isolation and reconnection: Demographic history and multiple contact zones of the green odorous frog (Odorrana margaretae) around the Sichuan Basin

The green odorous frog (Odorrana margaretae) displays a circular distribution around the Sichuan Basin of western China and possesses multiple replicate hybrid zones between lineages with high levels of divergence. To gain an understanding of the speciation process, we obtained 1540 SNPs from 29 populations and 227 individuals using ddRAD sequencing. Population structure analysis revealed three groups within the species: the West, North & South, and East groups. Demographic inference showed that they were initially isolated at ~2 million years ago, and subsequent post-glacial expansion produced the current circular distribution with four secondary contact zones. Hybridization in those zones involved lineages with various levels of divergence and produced greatly different patterns of introgression. Contact zones between the East and North & South groups (E-S and E-N) had contrast admixture levels but both showed a general lack of potential barrier loci. Meanwhile, the reconnection of the West and North & South groups produced two contact zones along the rim of the Basin. The S-W zone had extensive admixture while the N-W zone had limited admixture within a narrow geographic distance. Both showed substantial barrier effects, and a large number of potential barrier loci were shared. We also detected strong coupling among these loci. The N-W hybrid zone involved two highly-diverged lineages (F_ST = 0.704) and many loci have reached fixation around the hybrid zone. This study system offers a unique opportunity to understand the dynamics of introgression in contact zones and the architecture of reproductive isolation at different stages of speciation.

Genetic structure and demographic history of Indirana semipalmata, an endemic frog species of the Western Ghats, India

The evolutionary potential of a species mainly depends on the level of genetic variation in their populations. Maintenance of gene variation enables populations to adapt more quickly to environmental changes. The geographical gaps also influence the distribution and evolutionary history of many mountain frogs in the world. Hence, a sound knowledge in population genetic structure of a species will help understand its population dynamics and develop conservation strategies. In the context of facing threats to the amphibian fauna of Western Ghats due to habitat loss, we used both mitochondrial and nuclear DNA markers to investigate the genetic structure of an endemic frog species of the Western Ghats (Indirana semipalmata) with restricted distribution. The present study showed the importance of mountain gaps in shaping the species' structuring in the Western Ghats. Though a high genetic diversity was observed for the species when considering a single unit in the southern Western Ghats, the restricted gene flow on/between either side of the Shencottah gap with genetic clustering of the sampled populations may warrant a unique management plan for the species. The habitat fragmentation of the Western Ghats through anthropogenic activities may result in severe setbacks to the survival of the species in the future.

The genetic composition of Shina population from Gilgit-Baltistan, Pakistan based on mtDNA analyses

The present study aimed to gain insight into the genetic origin of the Shina population from Gilgit-Baltistan, Pakistan. We partially performed the mitochondrial DNA (mtDNA) control region of healthy unrelated individuals of Shina tribe residing in the remote area of Gilgit-Baltistan to investigate their maternal lineages. The present study is the first report about Shina's genetic structure, origin, and relationship with the surrounding north-western Pakistani population. The mtDNA sequences of the Shina samples were compared with the revised Cambridge Reference Sequence (rCRS) and the HVR-1 D-loop region was covered. The comparison with rCRS identified overall 38 haplotypes and 08 haplogroups for Shina samples. Among these haplotypes, 18 were shared by more than one individual of the Shina tribe. The obtained mtDNA sequences of Shina were compared with surrounding north-western Pakistani population groups, i.e. Kho, Kashmiri, and Pathan. The genetic diversity (i.e. 1.0424) and power of discrimination (i.e. 0.9266) of the Shina was found equivalent to surrounding north-western groups. The haplogroups frequencies, phylogenetic tree and network analysis identified the west Eurasian ancestral origin of Shina group with nearby maternal ancestral relationships with the Kashmiri population. However, no close genetic relationship of Shina was depicted with nearby residing Kho population group.

Caught in the act: Incipient speciation across a latitudinal gradient in a semifossorial mammal from Madagascar, the mole tenrec Oryzorictes hova (Tenrecidae)

Madagascar is one of the world's foremost biodiversity hotspots, yet a large portion of its flora and fauna remains undescribed and the driving forces of in situ diversification are not well understood. Recent studies have identified a widespread, latitudinally structured phylogeographic pattern in Madagascar's humid-forest mammals, amphibians, reptiles, and insects. Several factors may be driving this pattern, namely biogeographic barriers (i.e., rivers or valleys) or past episodes of forest contraction and expansion. In this study, we describe the phylogeographic structure of the small, semifossorial mammal Oryzorictes hova, one of Madagascar's two species of mole tenrec, found throughout Madagascar's eastern humid forest belt, from high-elevation montane forest to low-elevation forests, as well as disturbed habitat such as rice fields. Using one mitochondrial locus, four nuclear loci, and 31 craniomandibular measurements, we identified three distinct populations of O. hova associated with the northern, central, and southern regions of the island. We found little evidence of gene flow among these populations, so we treated each population as a potential species. We validated species limits using two Bayesian methods: BP&P, employing only DNA sequence data, and iBPP using both DNA and morphological data, and we assessed whether these methods are susceptible to producing false positive errors. Molecular and morphological data support the recognition of each of the three populations of O. hova as distinct species, but formal species descriptions will require additional data from type specimens. This study illustrates the importance of using integrative datasets, multiple methodological approaches, and extensive geographic sampling for species delimitation and adds evidence for a widespread phylogeographic pattern in Madagascar's humid forest taxa.

Evidence for an intrinsic factor promoting landscape genetic divergence in Madagascan leaf-litter frogs

The endemic Malagasy frog radiations are an ideal model system to study patterns and processes of speciation in amphibians. Large-scale diversity patterns of these frogs, together with other endemic animal radiations, led to the postulation of new and the application of known hypotheses of species diversification causing diversity patterns in this biodiversity hotspot. Both extrinsic and intrinsic factors have been studied in a comparative framework, with extrinsic factors usually being related to the physical environment (landscape, climate, river catchments, mountain chains), and intrinsic factors being clade-specific traits or constraints (reproduction, ecology, morphology, physiology). Despite some general patterns emerging from such large-scale comparative analyses, it became clear that the mechanism of diversification in Madagascar may vary among clades, and may be a multifactorial process. In this contribution, I test for intrinsic factors promoting population-level divergence within a clade of terrestrial, diurnal leaf-litter frogs (genus Gephyromantis) that has previously been shown to diversify according to extrinsic factors. Landscape genetic analyses of the microendemic species Gephyromantis enki and its widely distributed, larger sister species Gephyromantis boulengeri over a rugged landscape in the Ranomafana area shows that genetic variance of the smaller species cannot be explained by landscape resistance alone. Both topographic and riverine barriers are found to be important in generating this divergence. This case study yields additional evidence for the probable importance of body size in lineage diversification.

The Sepik River (Papua New Guinea) is not a dispersal barrier for lowland rain-forest frogs

Major tropical rivers have been suggested to be important dispersal barriers that increase the beta diversity of animal communities in lowland rain forests. We tested this hypothesis using assemblages of frogs in the floodplains of the Sepik River, a major river system in Papua New Guinea. We surveyed frogs at five sites within a continuous 150 × 500-km area of lowland rain forest bisected by the Sepik, using standardized visual and auditory survey techniques. We documented 769 frogs from 44 species. The similarity in species composition decreased with logarithm of geographical distance between the sites, which ranged from 82 to 465 km. The similarity decay did not depend on whether or not the compared sites were separated by the Sepik River or whether the species were aquatic or terrestrial breeders. Likewise, a DCA ordination of frog assemblages did not show separation of sites by the river as a significant factor explaining their composition. Our results suggest that even major rivers, such as the Sepik, may not act as dispersal barriers. Rivers may not limit the distribution of frogs and therefore have a limited effect on determining frog species abundance and assemblage structure in rain forests.