Single nucleotide polymorphisms (SNPs) provide unprecedented resolution of species boundaries, phylogenetic relationships, and genetic diversity in the mastiff bats (Molossus)

Comparative assessment of genotyping-by-sequencing and whole-exome sequencing for estimating genetic diversity and geographic structure in small sample sizes: insights from wild jaguar populations

Biologists currently have an assortment of high-throughput sequencing techniques allowing the study of population dynamics in increasing detail. The utility of genetic estimates depends on their ability to recover meaningful approximations while filtering out noise produced by artifacts. In this study, we empirically compared the congruence of two reduced representation approaches (genotyping-by-sequencing, GBS, and whole-exome sequencing, WES) in estimating genetic diversity and population structure using SNP markers typed in a small number of wild jaguar (Panthera onca) samples from South America. Due to its targeted nature, WES allowed for a more straightforward reconstruction of loci compared to GBS, facilitating the identification of true polymorphisms across individuals. We therefore used WES-derived metrics as a benchmark against which GBS-derived indicators were compared, adjusting parameters for locus assembly and SNP filtering in the latter. We observed significant variation in SNP call rates across samples in GBS datasets, leading to a recurrent miscalling of heterozygous sites. This issue was further amplified by small sample sizes, ultimately impacting the consistency of summary statistics between genotyping methods. Recognizing that the genetic markers obtained from GBS and WES are intrinsically different due to varying evolutionary pressures, particularly selection, we consider that our empirical comparison offers valuable insights and highlights critical considerations for estimating population genetic attributes using reduced representation datasets. Our results emphasize the critical need for careful evaluation of missing data and stringent filtering to achieve reliable estimates of genetic diversity and differentiation in elusive wildlife species.

Expanding the knowledge of the bat fauna of the Brazilian Caatinga: new geographical records of molossid bats (Chiroptera, Molossidae) for the Chapada Diamantina region, with taxonomic notes

The Caatinga, an exclusive biome in Brazil, is the largest tropical dry forest area in the Americas. It is characterized by a semi-arid climate and various soils that harbor a great diversity of flora and fauna. Novel records of aerial insectivorous bat species in the family Molossidae in the Chapada Diamantina, northeastern Brazil are presented. The study is based on field sampling of 115 molossid bat specimens from six genera and 12 taxonomically confirmed species, along with four taxa requiring further evaluation for definitive species identification. All specimens were obtained using mist nets around a small freshwater lagoon surrounded by semideciduous dry forest. The verified genera were Cynomops, Eumops, Molossops, Molossus, Neoplatymops, and Nyctinomops. Our findings enhance the understanding of bat diversity in the Brazilian Caatinga, with the first records of Eumopsdelticus, E.bonariensis, and Molossuscurrentium. The most abundant species were Molossusrufus, Eumopsglaucinus, Cynomopsplanirostris, Nyctinomopslaticaudatus, and Molossusmolossus. Previously unreported morphological and morphometric variations for these Caatinga taxa were examined. Additionally, information on sexual dimorphism in craniodental characteristics of Molossopstemminckii and variations in the presence of the sagittal crest in Neoplatymopsmattogrossensis are provided. Based on the voucher specimens from this study, the recognized number of species of Molossidae known from the Caatinga has increased to 21. Our results offer new insights into the taxonomy and biogeography of Neotropical molossids, highlighting their importance as members of bat communities in dry forest ecosystems from northeastern South America.

Back to the future: A preserved specimen validates the presence of Molossuspretiosus (Molossidae, Chiroptera) in Honduras

Molossuspretiosus is a molossid bat that has been thought to exist in Honduras. While some authors have suggested its range extends all the way to Mexico, others have placed its northernmost distribution in Nicaragua. We present evidence, based on one specimen collected in 2005, confirming the presence of this species in the Caribbean of Honduras within the Islas de la Bahía department. This discovery increases the count of known species within this family to 18 in the country and raises the total bat species count for Honduras to 114. We recommend a detailed study of historical specimens to confirm the identification of species that may have been misidentified as well as a thorough examination of molossids distributed in northern Honduras.

Association analysis of production traits of Japanese quail (Coturnix japonica) using restriction-site associated DNA sequencing

This study was designed to perform an association analysis and identify SNP markers associated with production traits of Japanese quail using restriction-site-associated DNA sequencing. Weekly body weight data from 805 quail were collected from hatching to 16 weeks of age. A total number of 3990 eggs obtained from 399 female quail were used to assess egg quality traits. Egg-related traits were measured at the beginning of egg production (first stage) and at 12 weeks of age (second stage). Five eggs were analyzed at each stage. Traits, such as egg weight, egg length and short axes, eggshell strength and weight, egg equator thickness, yolk weight, diameter, and colour, albumen weight, age of first egg, total number of laid eggs, and egg production rate, were assessed. A total of 383 SNPs and 1151 associations as well as 734 SNPs and 1442 associations were identified in relation to quail production traits using general linear model (GLM) and mixed linear model (MLM) approaches, respectively. The GLM-identified SNPs were located on chromosomes 1-13, 15, 17-20, 24, 26-28, and Z, underlying phenotypic traits, except for egg and albumen weight at the first stage and yolk yellowness at the second stage. The MLM-identified SNPs were positioned on defined chromosomes associated with phenotypic traits except for the egg long axis at the second stage of egg production. Finally, 35 speculated genes were identified as candidate genes for the targeted traits based on their nearest positions. Our findings provide a deeper understanding and allow a more precise genetic improvement of production traits of Galliformes, particularly in Japanese quail.

Unmasking microsatellite deceptiveness and debunking hybridization with SNPs in four marine copepod species of Calanus

Interspecific hybridization events are on the rise in natural systems due to climate change disrupting species barriers. Across taxa, microsatellites have long been the molecular markers of choice to identify admixed individuals. However, with the advent of high-throughput sequencing easing the generation of genome-wide datasets, incorrect reports of hybridization resulting from microsatellite technical artefacts have been uncovered in a growing number of taxa. In the marine zooplankton genus Calanus (Copepoda), whose species are used as climate change indicators, microsatellite markers have suggested hybridization between C. finmarchicus and C. glacialis, while other nuclear markers (InDels) never detected any admixed individuals, leaving the scientific community divided. Here, for the first time, we investigated the potential for hybridization among C. finmarchicus, C. glacialis, C. helgolandicus and C. hyperboreus using two large and independent SNP datasets. These were derived firstly from a protocol of target-capture applied to 179 individuals collected from 17 sites across the North Atlantic and Arctic Oceans, including sympatric areas, and second from published RNA sequences. All SNP-based analyses were congruent in showing that Calanus species are distinct and do not appear to hybridize. We then thoroughly re-assessed the microsatellites showing hybrids, with the support of published transcriptomes, and identified technical issues plaguing eight out of 10 microsatellites, including size homoplasy, paralogy, potential for null alleles and even two primer pairs targeting the same locus. Our study illustrates how deceptive microsatellites can be when applied to the investigation of hybridization.

Updated list of the mammals of Costa Rica, with notes on recent taxonomic changes

Although Costa Rica occupies a mere 0.03% of the Earths land area, it nevertheless has recorded within its borders approximately 5% of the global diversity of mammals, thus making it one of the worlds megadiverse countries. Over the past ten years, 22 species have been added to the countrys inventory, bringing the total number known as here documented to 271; Chiroptera account for ten of these, having grown to 124 from 114; rodents have increased by eight species, from 47 to 55, with the caveat that we include three invasive species of Muridae that have gone feral. In contrast, the number of orders has decreased by one, by Artiodactyla incorporating the former Cetacea. Notes are provided for all taxonomic novelties since the last update. Since the first taxonomic compendium of the mammals of Costa Rica in 1869, the number of known species has grown by approximately 1.22 species year-1 (R2 = 0.96). Since 1983 however, this growth rate has been 1.64 species year-1 (R2 = 0.98). Despite this strong growth, an asymptote in the number of known species has not been reached. Conservation remains a primary need: over 60% of the countrys mammal species show population trends that are decreasing (13%), unknown (37%), or not assessed (11%), based on IUCN criteria. These analyses suggest that much remains to be known regarding the number of mammal species living in Costa Rica, but also that much more remains to be done to safeguard Costa Ricas exceptional biodiversity heritage.

The genome sequence of Molossusnigricans (Chiroptera, Molossidae; Miller, 1902)

We present a genome assembly from an individual male Molossus nigricans (Chordata; Mammalia; Chiroptera; Molossidae). The genome sequence is 2.41 gigabases in span. The majority of the assembly is scaffolded into 24 chromosomal pseudomolecules, with the X sex chromosome assembled.

Morphological and molecular analyses confirm the presence of the newly described bat species Molossus melini in the Espinal ecoregion from Argentina

Until now, Molossus melini was known only from its type locality, in the south of Santa Fe province, Argentina. Specimens of this species were collected in 2021 from a roost in a Fraxinus tree of the urban woodland of Paraná city, Entre Ríos province, Argentina. Bat identification was made by comparing external and cranial characters and measurements with those reported in the bibliography and corroborated by a phylogenetic analysis based on the cytochrome b gene. Also, multivariate morphometric analyses showed that cranial measurements, but not external ones, are informative enough to discriminate M. melini from the other Molossus species in Argentina (M. currentium, M. fluminensis, and M. molossus). This new record extends the distribution of M. melini from the south of Santa Fe province by 230 km to the northeast and represents the first record of the species in the Espinal ecoregion.

Supplementary information

The online version contains supplementary material available at 10.1007/s13364-023-00679-1.

Structural Consequences of IRS-2 nsSNPs and Implication for Insulin Receptor Substrate-2 Protein Stability

Single-Nucleotide Polymorphisms (SNPs) are common genetic variations implicated in human diseases. The non-synonymous SNPs (nsSNPs) affect the proteins' structures and their molecular interactions with other interacting proteins during the accomplishment of biochemical processes. This ultimately causes proteins functional perturbation and disease phenotypes. The Insulin receptor substrate-2 (IRS-2) protein promotes glucose absorption and participates in the biological regulation of glucose metabolism and energy production. Several IRS-2 SNPs are reported in association with type 2 diabetes and obesity in human populations. However, there are no comprehensive reports about the protein structural consequences of these nsSNPs. Keeping in view the pathophysiological consequences of the IRS-2 nsSNPs, we designed the current study to understand their possible structural impact on coding protein. The prioritized list of the deleterious IRS-2 nsSNPs was acquired from multiple bioinformatics resources, including VEP (SIFT, PolyPhen, and Condel), PROVEAN, SNPs&GO, PMut, and SNAP2. The protein structure stability assessment of these nsSNPs was performed by MuPro and I-Mutant-3.0 servers via structural modeling approaches. The atomic-level structural and molecular dynamics (MD) impact of these nsSNPs were examined using GROMACS 2019.2 software package. The analyses initially predicted 8 high-risk nsSNPs located in the highly conserved regions of IRS-2. The MD simulation analysis eventually prioritized the N232Y, R218C, and R104H nsSNPs that predicted to significantly compromise the structure stability and may affect the biological function of IRS-2. These nsSNPs are predicted as high-risk candidates for diabetes and obesity. The validation of protein structural impact of these shortlisted nsSNPs may provide biochemical insight into the IRS-2-mediated type-2 diabetes.

Out of thin air: surveying tropical bat roosts through air sampling of eDNA

Understanding roosting behaviour is essential to bat conservation and biomonitoring, often providing the most accurate methods of assessing bat population size and health. However, roosts can be challenging to survey, e.g., physically impossible to access or presenting risks for researchers. Disturbance during monitoring can also disrupt natural bat behaviour and present material risks to the population such as disrupting hibernation cycles. One solution to this is the use of non-invasive monitoring approaches. Environmental (e)DNA has proven especially effective at detecting rare and elusive species particularly in hard-to-reach locations. It has recently been demonstrated that eDNA from vertebrates is carried in air. When collected in semi-confined spaces, this airborne eDNA can provide remarkably accurate profiles of biodiversity, even in complex tropical communities. In this study, we deploy novel airborne eDNA collection for the first time in a natural setting and use this approach to survey difficult to access potential roosts in the neotropics. Using airborne eDNA, we confirmed the presence of bats in nine out of 12 roosts. The identified species matched previous records of roost use obtained from photographic and live capture methods, thus demonstrating the utility of this approach. We also detected the presence of the white-winged vampire bat (Diaemus youngi) which had never been confirmed in the area but was long suspected based on range maps. In addition to the bats, we detected several non-bat vertebrates, including the big-eared climbing rat (Ototylomys phyllotis), which has previously been observed in and around bat roosts in our study area. We also detected eDNA from other local species known to be in the vicinity. Using airborne eDNA to detect new roosts and monitor known populations, particularly when species turnover is rapid, could maximize efficiency for surveyors while minimizing disturbance to the animals. This study presents the first applied use of airborne eDNA collection for ecological analysis moving beyond proof of concept to demonstrate a clear utility for this technology in the wild.

The right tool for the right question: contrasting biogeographic patterns in the notothenioid fish Harpagifer spp. along the Magellan Province

Molecular-based analysis has become a fundamental tool to understand the role of Quaternary glacial episodes. In the Magellan Province in southern South America, ice covering during the last glacial maximum (20 ka) radically altered the landscape/seascape, speciation rates and distribution of species. For the notothenioid fishes of the genus Harpagifer, in the area are described two nominal species. Nevertheless, this genus recently colonized South America from Antarctica, providing a short time for speciation processes. Combining DNA sequences and genotyping-by-sequencing SNPs, we evaluated the role of Quaternary glaciations over the patterns of genetic structure in Harpagifer across its distribution in the Magellan Province. DNA sequences showed low phylogeographic structure, with shared and dominant haplotypes between nominal species, suggesting a single evolutionary unit. SNPs identified contrastingly two groups in Patagonia and a third well-differentiated group in the Falkland/Malvinas Islands with limited and asymmetric gene flow. Linking the information of different markers allowed us to infer the relevance of postglacial colonization mediated by the general oceanographic circulation patterns. Contrasting rough- and fine-scale genetic patterns highlights the relevance of combined methodologies for species delimitation, which, depending on the question to be addressed, allows discrimination among phylogeographic structure, discarding incipient speciation, and contemporary spatial differentiation processes.

A new species of Molossus (Chiroptera: Molossidae) from Argentina

Fourteen species of bats in the genus Molossus currently are recognized in the Neotropical region; only three are known from Argentina. Here, we describe a new species based on specimens collected in the province of Santa Fe, Argentina, in the Pampa ecoregion. The new species can be distinguished from its congeners by its general strong orange coloration, forearm length > 41 mm, dorsal hairs bicolor and long (~5 mm), infraorbital foramen laterally oriented, and long and forward-projected (pincer-like) upper incisors. The external and cranial morphology of the new species are described and comparisons made with other species of similar size and with those present in its distributional area. Wilcoxon tests and multivariate analyses (nonmetric multidimensional scaling ordination and PERMANOVA) were carried out to determine the morphometric differences between the new species and other seven species of Molossus. The species tree, estimated by *BEAST from the concatenation of mitochondrial and nuclear genes, suggests that Molossus sp. nov. is basal within the clade formed by {{M. aztecus, M. rufus}, {{M. currentium, M. pretiosus}, M. sinaloae}} with a posterior probability of 0.82.

Quantitative trait loci for growth-related traits in Japanese quail (Coturnix japonica) using restriction-site associated DNA sequencing

This study aimed to identify quantitative trait loci (QTLs) for growth-related traits by constructing a genetic linkage map based on single nucleotide polymorphism (SNP) markers in Japanese quail. A QTL mapping population of 277 F₂ birds was obtained from an intercross between a male of a large-sized strain and three females of a normal-sized strain. Body weight (BW) was measured weekly from hatching to 16 weeks of age. Non-linear regression growth models of Weibull, Logistic, Gompertz, Richards, and Brody were analyzed, and growth curve parameters of Richards was selected as the best model to describe the quail growth curve of the F₂ birds. Restriction-site associated DNA sequencing developed 125 SNP markers that were informative between their parental strains. The SNP markers were distributed on 16 linkage groups that spanned 795.9 centiMorgan (cM) with an average marker interval of 7.3 cM. QTL analysis of phenotypic traits revealed four main-effect QTLs. Detected QTLs were located on chromosomes 1 and 3 and were associated with BW from 4 to 16 weeks of age and asymptotic weight of Richards model at genome-wide significant at 1% or 5% level. No QTL was detected for BW from 0 to 3 weeks of age. This is the first report identified QTLs for asymptotic weight of the Richards parameter in Japanese quail. These results highlight that the combination of QTL studies and the RAD-seq method will aid future breeding programs identify genes underlying the QTL and the application of marker-assisted selection in the poultry industry, particularly the Japanese quail.

Phylogenomic species delimitation dramatically reduces species diversity in an Antarctic adaptive radiation

Application of genetic data to species delimitation often builds confidence in delimitations previously hypothesized using morphological, ecological, and geographic data and frequently yields recognition of previously-undescribed cryptic diversity. However, a recent critique of genomic data-based species delimitation approaches is that they have the potential to conflate population structure with species diversity, resulting in taxonomic oversplitting. The need for an integrative approach to species delimitation, in which molecular, morphological, ecological, and geographic lines of evidence are evaluated together, is becoming increasingly apparent. Here, we integrate phylogenetic, population genetic, and coalescent analyses of genome-wide sequence data with investigation of variation in multiple morphological traits to delimit species within the Antarctic barbeled plunderfishes (Artedidraconidae: Pogonophryne). Pogonophryne currently comprises 29 valid species, most of which are distinguished solely by variation in ornamentation of the mental barbel that projects from the lower jaw, a structure previously shown to vary widely within a single species. However, our genomic and phenotypic analyses result in a dramatic reduction in the number of distinct species recognized within the clade, providing evidence to support the recognition of no more than six species. We propose to synonymize 24 of the currently recognized species with five species of Pogonophryne. We find genomic and phenotypic evidence for a new species of Pogonophryne from specimens collected in the Ross Sea. Our findings represent a rare example in which application of molecular data provides evidence of taxonomic oversplitting on the basis of morphology, clearly demonstrating the utility of an integrative species delimitation framework.

Resolving species boundaries in a recent radiation with the Angiosperms353 probe set: the Lomatium packardiae/L. anomalum clade of the L. triternatum (Apiaceae) complex

PREMISE

Speciation not associated with morphological shifts is challenging to detect unless molecular data are employed. Using Sanger-sequencing approaches, the Lomatium packardiae/L. anomalum subcomplex within the larger Lomatium triternatum complex could not be resolved. Therefore, we attempt to resolve these boundaries here.

METHODS

The Angiosperms353 probe set was employed to resolve the ambiguity within Lomatium triternatum species complex using 48 accessions assigned to L. packardiae, L. anomalum, or L. triternatum. In addition to exon data, 54 nuclear introns were extracted and were complete for all samples. Three approaches were used to estimate evolutionary relationships and define species boundaries: STACEY, a Bayesian coalescent-based species tree analysis that takes incomplete lineage sorting into account; ASTRAL-III, another coalescent-based species tree analysis; and a concatenated approach using MrBayes. Climatic factors, morphological characters, and soil variables were measured and analyzed to provide additional support for recovered groups.

RESULTS

The STACEY analysis recovered three major clades and seven subclades, all of which are geographically structured, and some correspond to previously named taxa. No other analysis had full agreement between recovered clades and other parameters. Climatic niche and leaflet width and length provide some predictive ability for the major clades.

CONCLUSIONS

The results suggest that these groups are in the process of incipient speciation and incomplete lineage sorting has been a major barrier to resolving boundaries within this lineage previously. These results are hypothesized through sequencing of multiple loci and analyzing data using coalescent-based processes.

Mapping of Quantitative Trait Loci Controlling Egg-Quality and -Production Traits in Japanese Quail (Coturnix japonica) Using Restriction-Site Associated DNA Sequencing

This research was conducted to identify quantitative trait loci (QTL) associated with egg-related traits by constructing a genetic linkage map based on single nucleotide polymorphism (SNP) markers using restriction-site associated DNA sequencing (RAD-seq) in Japanese quail. A total of 138 F₂ females were produced by full-sib mating of F₁ birds derived from an intercross between a male of the large-sized strain with three females of the normal-sized strain. Eggs were investigated at two different stages: the beginning stage of egg-laying and at 12 weeks of age (second stage). Five eggs were analyzed for egg weight, lengths of the long and short axes, egg shell strength and weight, yolk weight and diameter, albumen weight, egg equator thickness, and yolk color (L*, a*, and b* values) at each stage. Moreover, the age at first egg, the cumulative number of eggs laid, and egg production rate were recorded. RAD-seq developed 118 SNP markers and mapped them to 13 linkage groups using the Map Manager QTX b20 software. Markers were spanned on 776.1 cM with an average spacing of 7.4 cM. Nine QTL were identified on chromosomes 2, 4, 6, 10, 12, and Z using the simple interval mapping method in the R/qtl package. The QTL detected affected 10 egg traits of egg weight, lengths of the long and short axes of egg, egg shell strength, yolk diameter and weight, albumen weight, and egg shell weight at the beginning stage, yellowness of the yolk color at the second stage, and age at first egg. This is the first report to perform a quail QTL analysis of egg-related traits using RAD-seq. These results highlight the effectiveness of RAD-seq associated with targeted QTL and the application of marker-assisted selection in the poultry industry, particularly in the Japanese quail.

The History of Lentil (Lens culinaris subsp. culinaris) Domestication and Spread as Revealed by Genotyping-by-Sequencing of Wild and Landrace Accessions

Protein-rich legumes accompanied carbohydrate-rich cereals since the beginning of agriculture and yet their domestication history is not as well understood. Lentil (Lens culinaris Medik. subsp. culinaris) was first cultivated in Southwest Asia (SWA) 8000-10,000 years ago but archeological evidence is unclear as to how many times it may have been independently domesticated, in which SWA region(s) this may have happened, and whether wild species within the Lens genus have contributed to the cultivated gene pool. In this study, we combined genotyping-by-sequencing (GBS) of 190 accessions from wild (67) and domesticated (123) lentils from the Old World with archeological information to explore the evolutionary history, domestication, and diffusion of lentils to different environments. GBS led to the discovery of 87,647 single-nucleotide polymorphisms (SNPs), which allowed us to infer the phylogeny of genus Lens. We confirmed previous studies proposing four groups within it. The only gene flow detected was between cultivated varieties and their progenitor (L. culinaris subsp. orientalis) albeit at very low levels. Nevertheless, a few putative hybrids or naturalized cultivars were identified. Within cultivated lentil, we found three geographic groups. Phylogenetics, population structure, and archeological data coincide in a scenario of protracted domestication of lentils, with two domesticated gene pools emerging in SWA. Admixed varieties are found throughout their range, suggesting a relaxed selection process. A small number of alleles involved in domestication and adaptation to climatic variables were identified. Both novel mutation and selection on standing variation are presumed to have played a role in adaptation of lentils to different environments. The results presented have implications for understanding the process of plant domestication (past), the distribution of genetic diversity in germplasm collections (present), and targeting genes in breeding programs (future).

Does evolution of echolocation calls and morphology in Molossus result from convergence or stasis?

Although many processes of diversification have been described to explain variation of morphological traits within clades that have obvious differentiation among taxa, not much is known about these patterns in complexes of cryptic species. Molossus is a genus of bats that is mainly Neotropical, occurring from the southeastern United States to southern Argentina, including the Caribbean islands. Molossus comprises some groups of species that are morphologically similar but phylogenetically divergent, and other groups of species that are genetically similar but morphologically distinct. This contrast allows investigation of unequal trait diversification and the evolution of morphological and behavioural characters. In this study, we assessed the role of phylogenetic history in a genus of bat with three cryptic species complexes, and evaluated if morphology and behavior are evolving concertedly. The Genotype by Sequence genomic approach was used to build a species-level phylogenetic tree for Molossus and to estimate the ancestral states of morphological and echolocation call characters. We measured the correlation of phylogenetic distances to morphological and echolocation distances, and tested the relationship between morphology and behavior when the effect of phylogeny is removed. Morphology evolved via a mosaic of convergence and stasis, whereas call design was influenced exclusively through local adaptation and convergent evolution. Furthermore, the frequency of echolocation calls is negatively correlated with the size of the bat, but other characters do not seem to be evolving in concert. We hypothesize that slight variation in both morphology and behaviour among species of the genus might result from niche specialization, and that traits evolve to avoid competition for resources in similar environments.

Leucism and updated geographic distribution of Molossus nigricans Miller, 1902 (Chiroptera: Molossidae) in Honduras

We present the first record of leucism of Molossus nigricans a chromatic disorder rarely documented in mammals of Honduras, as well as species’ biological and ecological information. In addition, we present the currently known distribution and altitudinal range of M. nigricans in Honduras, based on acoustic records, opportunistic findings and previous records. We confirm the presence of M. nigricans in sixteen departments of which four are new departmental records.

Optimization of Genotype by Sequencing data for phylogenetic purposes

• Herein we propose a framework for assembling and analyzing Genotype by Sequencing (GBS) data to better understand evolutionary relationships within a group of closely related species using the mastiff bats (Molossus) as our model system. Many species within this genus have low-levels of genetic variation within and between morphologically distinct species, and the relationships among them remain unresolved using traditional Sanger sequencing methods. Given that both de novo and reference genome pipelines can be used to assemble next generation sequences, and that several tree inference methodologies have been proposed for single nucleotide polymorphism (SNP) data, we test whether different alignments and phylogenetic approaches produce similar results. We also examined how the process of SNP identification and mapping can affect the consistency of the analyses. Different alignments and phylogenetic inferences produced consistent results, supporting the GBS approach for answering evolutionary questions on a macroevolutionary scale when the genetic distance among phenotypically identifiable clades is low. We highlight the importance of exploring the relationships among groups using different assembly assumptions and also distinct phylogenetic inference methods, particularly when addressing phylogenetic questions in genetic and morphologically conservative taxa. • The method uses the comparison of several filter settings, alignments, and tree inference approaches on Genotype by Sequencing data. • Consistent results were found among several approaches. • The methodology successfully recovered well supported species boundaries and phylogenetic relationships among species of mastiff bats not hypothesized by previous methods.

Next generation sequencing data in the phylogenetic relationships of the genus Molossus (Chiroptera, Molossidae)

The mastiff bat Molossus is a broadly distributed genus within the family Molossidae. Molossus includes groups of species that are either morphologically or genetically very similar, rendering the taxonomy of this genus confusing and unstable. In this paper, we provide inferred phylogenetic relationships of Molossus based on the genotype by sequencing approach from 189 specimens of three species of New World mastiff bats (Molossus, Promops, and Eumops). We also present data on divergent tree topologies produced by alignments using de novo and reference genome approaches and distinct phylogenetic methods (maximum likelihood and coalescent approaches). These data provide the first highly resolved phylogenetic tree for Molossus, not recovered by previous studies using Sanger sequencing. Our dataset brings new insights on relationships among species and show how different approaches might affect phylogenetic resolution and topologies.

Comparative phylogeography of mainland and insular species of Neotropical molossid bats (Molossus)

Historical events, habitat preferences, and geographic barriers might result in distinct genetic patterns in insular versus mainland populations. Comparison between these two biogeographic systems provides an opportunity to investigate the relative role of isolation in phylogeographic patterns and to elucidate the importance of evolution and demographic history in population structure. Herein, we use a genotype-by-sequencing approach (GBS) to explore population structure within three species of mastiff bats (Molossus molossus, M. coibensis, and M. milleri), which represent different ecological histories and geographical distributions in the genus. We tested the hypotheses that oceanic straits serve as barriers to dispersal in Caribbean bats and that isolated island populations are more likely to experience genetic drift and bottlenecks in comparison with highly connected ones, thus leading to different phylogeographic patterns. We show that population structures vary according to general habitat preferences, levels of population isolation, and historical fluctuations in climate. In our dataset, mainland geographic barriers played only a small role in isolation of lineages. However, oceanic straits posed a partial barrier to the dispersal for some populations within some species (M. milleri), but do not seem to disrupt gene flow in others (M. molossus). Lineages on distant islands undergo genetic bottlenecks more frequently than island lineages closer to the mainland, which have a greater exchange of haplotypes.