The utility of single nucleotide polymorphisms in inferences of population history

Taking Advantage of the Genomics Revolution for Monitoring and Conservation of Chondrichthyan Populations

Chondrichthyes (sharks, rays, skates and chimaeras) are among the oldest extant predators and are vital to top-down regulation of oceanic ecosystems. They are an ecologically diverse group occupying a wide range of habitats and are thus, exploited by coastal, pelagic and deep-water fishing industries. Chondrichthyes are among the most data deficient vertebrate species groups making design and implementation of regulatory and conservation measures challenging. High-throughput sequencing technologies have significantly propelled ecological investigations and understanding of marine and terrestrial species’ populations, but there remains a paucity of NGS based research on chondrichthyan populations. We present a brief review of current methods to access genomic and metagenomic data from Chondrichthyes and discuss applications of these datasets to increase our understanding of chondrichthyan taxonomy, evolution, ecology and population structures. Last, we consider opportunities and challenges offered by genomic studies for conservation and management of chondrichthyan populations.

Embracing heterogeneity: coalescing the Tree of Life and the future of phylogenomics

Building the Tree of Life (ToL) is a major challenge of modern biology, requiring advances in cyberinfrastructure, data collection, theory, and more. Here, we argue that phylogenomics stands to benefit by embracing the many heterogeneous genomic signals emerging from the first decade of large-scale phylogenetic analysis spawned by high-throughput sequencing (HTS). Such signals include those most commonly encountered in phylogenomic datasets, such as incomplete lineage sorting, but also those reticulate processes emerging with greater frequency, such as recombination and introgression. Here we focus specifically on how phylogenetic methods can accommodate the heterogeneity incurred by such population genetic processes; we do not discuss phylogenetic methods that ignore such processes, such as concatenation or supermatrix approaches or supertrees. We suggest that methods of data acquisition and the types of markers used in phylogenomics will remain restricted until a posteriori methods of marker choice are made possible with routine whole-genome sequencing of taxa of interest. We discuss limitations and potential extensions of a model supporting innovation in phylogenomics today, the multispecies coalescent model (MSC). Macroevolutionary models that use phylogenies, such as character mapping, often ignore the heterogeneity on which building phylogenies increasingly rely and suggest that assimilating such heterogeneity is an important goal moving forward. Finally, we argue that an integrative cyberinfrastructure linking all steps of the process of building the ToL, from specimen acquisition in the field to publication and tracking of phylogenomic data, as well as a culture that values contributors at each step, are essential for progress.

Interaction between TP53 and XRCC1 increases susceptibility to cervical cancer development: a case control study

BACKGROUND

Cervical cancer is the 4th highest cause of female reproductive tract malignancies. Multiple loci have been identified as important determinant factors for tumor susceptibility. In this report, we aimed to explore the roles of gene polymorphisms affecting x-ray repair cross complementing 1 (XRCC1), the tumor protein p53 (TP53), and fibroblast growth factor receptor 3 (FGFR3) in the context of susceptibility to cervical cancer. Additionally, we assessed the impact of single nucleotide polymorphism-single nucleotide polymorphism (SNP-SNP) interaction of these three genes in the context of cervical cancer risk in Chinese women.

METHODS

A case-control study consisted of 340 women located in Chongqing. Of these women, 121 were diagnosed with cervical cancer, 118 served as healthy controls, and 101 were specifically recruited elderly patients above the age of 80 who showed no history of cervical cancer. Three SNPs (XRCC1 rs25487, TP53 rs1042522, and FGFR3 rs121913483) were examined using mutation analysis of mismatch amplification PCR (MAMA-PCR) on samples obtained from peripheral blood.

RESULTS

Our results indicated that females from southwestern China all exhibited a wild-type phenotype at FGFR3 rs121913483. We also observed that the rs25487 mutation was significantly increased within the cervical cancer population. A 2-locus SNP-SNP interaction pattern (rs25487 and rs1042522) was significantly associated with cervical cancer risk (cases vs. negative controls: OR = 4.63, 95% CI = 1.83-11.75; cases vs. elderly group: OR = 17.61, 95% CI = 4.34-71.50).

CONCLUSIONS

This is the first study to identify a novel interaction between the XRCC1 and TP53 genes that is highly associated with susceptibility to cervical cancer risk in a female population in southwestern China.

Biotic and Abiotic Constraints in Mungbean Production-Progress in Genetic Improvement

Mungbean [Vigna radiata (L.) R. Wilczek var. radiata] is an important food and cash legume crop in Asia. Development of short duration varieties has paved the way for the expansion of mungbean into other regions such as Sub-Saharan Africa and South America. Mungbean productivity is constrained by biotic and abiotic factors. Bruchids, whitefly, thrips, stem fly, aphids, and pod borers are the major insect-pests. The major diseases of mungbean are yellow mosaic, anthracnose, powdery mildew, Cercospora leaf spot, halo blight, bacterial leaf spot, and tan spot. Key abiotic stresses affecting mungbean production are drought, waterlogging, salinity, and heat stress. Mungbean breeding has been critical in developing varieties with resistance to biotic and abiotic factors, but there are many constraints still to address that include the precise and accurate identification of resistance source(s) for some of the traits and the traits conferred by multi genes. Latest technologies in phenotyping, genomics, proteomics, and metabolomics could be of great help to understand insect/pathogen-plant, plant-environment interactions and the key components responsible for resistance to biotic and abiotic stresses. This review discusses current biotic and abiotic constraints in mungbean production and the challenges in genetic improvement.

A bioinformatic pipeline for identifying informative SNP panels for parentage assignment from RADseq data

The development of high-throughput sequencing technologies is dramatically increasing the use of single nucleotide polymorphisms (SNPs) across the field of genetics, but most parentage studies of wild populations still rely on microsatellites. We developed a bioinformatic pipeline for identifying SNP panels that are informative for parentage analysis from restriction site-associated DNA sequencing (RADseq) data. This pipeline includes options for analysis with or without a reference genome, and provides methods to maximize genotyping accuracy and select sets of unlinked loci that have high statistical power. We test this pipeline on small populations of Mexican gray wolf and bighorn sheep, for which parentage analyses are expected to be challenging due to low genetic diversity and the presence of many closely related individuals. We compare the results of parentage analysis across SNP panels generated with or without the use of a reference genome, and between SNPs and microsatellites. For Mexican gray wolf, we conducted parentage analyses for 30 pups from a single cohort where samples were available from 64% of possible mothers and 53% of possible fathers, and the accuracy of parentage assignments could be estimated because true identities of parents were known a priori based on field data. For bighorn sheep, we conducted maternity analyses for 39 lambs from five cohorts where 77% of possible mothers were sampled, but true identities of parents were unknown. Analyses with and without a reference genome produced SNP panels with ≥95% parentage assignment accuracy for Mexican gray wolf, outperforming microsatellites at 78% accuracy. Maternity assignments were completely consistent across all SNP panels for the bighorn sheep, and were 74.4% consistent with assignments from microsatellites. Accuracy and consistency of parentage analysis were not reduced when using as few as 284 SNPs for Mexican gray wolf and 142 SNPs for bighorn sheep, indicating our pipeline can be used to develop SNP genotyping assays for parentage analysis with relatively small numbers of loci.

Genome-wide SNP analyses reveal high gene flow and signatures of local adaptation among the scalloped spiny lobster (Panulirus homarus) along the Omani coastline

BACKGROUND

The scalloped spiny lobster (Panulirus homarus) is a popular seafood commodity worldwide and an important export item from Oman. Annual catches in commercial fisheries are in serious decline, which has resulted in calls for the development of an integrated stock management approach. In Oman, the scalloped spiny lobster is currently treated as a single management unit (MU) or stock and there is an absence of information on the genetic population structure of the species that can inform management decisions, particularly at a fine-scale level. This work is the first to identify genome-wide single nucleotide polymorphisms (SNPs) for P. homarus using Diversity Arrays Technology sequencing (DArT-seq) and to elucidate any stock structure in the species.

RESULTS

After stringent filtering, 7988 high utility SNPs were discovered and used to assess the genetic diversity, connectivity and structure of P. homarus populations from Al Ashkharah, Masirah Island, Duqm, Ras Madrakah, Haitam, Ashuwaymiyah, Mirbat and Dhalkut landing sites. Pairwise FST estimates revealed low differentiation among populations (pairwise FST range = - 0.0008 - 0.0021). Analysis of genetic variation using putatively directional FST outliers (504 SNPs) revealed higher and significant pairwise differentiation (p < 0.01) for all locations, with Ashuwaymiyah being the most diverged population (Ashuwaymiyah pairwise FST range = 0.0288-0.0736). Analysis of population structure using Discriminant Analysis of Principal Components (DAPC) revealed a broad admixture among P. homarus, however, Ashuwaymiyah stock appeared to be potentially under local adaptive pressures. Fine scale analysis using Netview R provided further support for the general admixture of P. homarus.

CONCLUSIONS

Findings here suggested that stocks of P. homarus along the Omani coastline are admixed. Yet, fishery managers need to treat the lobster stock from Ashuwaymiyah with caution as it might be subject to local adaptive pressures. We emphasize further study with larger number of samples to confirm the genetic status of the Ashuwaymiyah stock. The approach utilised in this study has high transferability in conservation and management of other marine stocks with similar biological and ecological attributes.

The dual role of Amazonian rivers in the generation and maintenance of avian diversity

The Amazon River and its major tributaries delimit the distributions of hundreds of terrestrial taxa. It remains unclear whether river-bounded distributions and taxon replacements reflect the historical role of rivers in generating species diversity as vicariant forces, or are the result of their role as secondary barriers, maintaining current levels of species diversity by inhibiting gene flow and population introgression. We use a community-wide comparative phylogeographic and phylogenetic approach to address the roles that the Rio Negro and the Rio Branco play in the avian speciation process in the Guiana Shield. Examining 74 pairs of ecologically similar geographic replacements that turn over across the lower Negro, we found substantial variation in the levels of genetic divergence and the inferred timing of diversification among pairs, ranging from ~0.24 to over 8 million years (Ma ago). The breadth of this variation is inconsistent with a single, shared speciation event. Coalescent simulations also rejected a simultaneous divergence scenario for pairs divided by the Rio Branco but could not reject a single diversification pulse for a subset of 12 pairs of taxa divided by the upper Negro. These results are consistent with recent geomorphological hypotheses regarding the origins of these rivers. Phylogenetically, taxon pairs represent a blend of sister (~40%) and nonsister taxa (~60%), consistent with river-associated allopatric or peripatric speciation and secondary contact, respectively. Our data provide compelling evidence that species turnover across the Rio Negro basin encompasses a mixture of histories, supporting a dual role for Amazonian rivers in the generation and maintenance of biological diversity.

Targeted resequencing of coding DNA sequences for SNP discovery in nonmodel species

Targeted capture coupled with high-throughput sequencing can be used to gain information about nuclear sequence variation at hundreds to thousands of loci. Divergent reference capture makes use of molecular data of one species to enrich target loci in other (related) species. This is particularly valuable for nonmodel organisms, for which often no a priori knowledge exists regarding these loci. Here, we have used targeted capture to obtain data for 809 nuclear coding DNA sequences (CDS) in a nonmodel organism, the Eurasian lynx Lynx lynx, using baits designed with the help of the published genome of a related model organism (the domestic cat Felis catus). Using this approach, we were able to survey intraspecific variation at hundreds of nuclear loci in L. lynx across the species' European range. A large set of biallelic candidate SNPs was then evaluated using a high-throughput SNP genotyping platform (Fluidigm), which we then reduced to a final 96 SNP-panel based on assay performance and reliability; validation was carried out with 100 additional Eurasian lynx samples not included in the SNP discovery phase. The 96 SNP-panel developed from CDS performed very successfully in the identification of individuals and in population genetic structure inference (including the assignment of individuals to their source population). In keeping with recent studies, our results show that genic SNPs can be valuable for genetic monitoring of wildlife species.

Variability in the anatoxin gene clusters of Cuspidothrix issatschenkoi from Germany, New Zealand, China and Japan

Anatoxin-a and homoanatoxin-a are neurotoxic cyanotoxins produced by benthic and planktonic cyanobacteria worldwide. These toxins are produced by the cyanobacterial genera Dolichospermum, Cuspidothrix, Phormidium, Oscillatoria, Tychonema and Cylindrospermum. In the present study the ana gene clusters (anaA-anaG; c. 21.1 kilobases) of two anatoxin producing Cuspidothrix issatschenkoi strains from Germany: (NIVA-CYA 711) and New Zealand (CAWBG02) were sequenced and compared with the ana gene clusters of two C. issatschenkoi strains from Japan (RM-6 and LBRI48) and one from China (CHABD3). All five ana gene clusters are characterized by the same gene order for anaA-anaG. Similarities were highest (99.56–99.57%) between German (NIVA-CYA 711), New Zealand (CAWBG02) and Chinese (CHABD3) strains. Similarities were lower (91.40–91.67%) when compared to the Japanese strains (RM-6 and LBRI48). Collectively, 2,037 variable sites (328 single nucleotide polymorphisms and 9 insertions/deletions, comprising 1,709 nucleotides) were found in the ana gene clusters of the German, New Zealand and Japanese strains compared to the Chinese strain (CHABD3). The ana gene clusters of the German (NIVA-CYA 711), New Zealand (CAWBG02) and Japanese (RM-6 and LBRI48) strains were characterized by 83, 84, 255 and 231 SNP’s compared to the Chinese strain (CHABD3), respectively. The anaE and anaF genes showed the highest variability in all five strains and are recommended as the best genetic markers for further phylogenetic studies of the ana gene cluster from C. issatschenkoi.

Clipperton Atoll as a model to study small marine populations: Endemism and the genomic consequences of small population size

Estimating population sizes and genetic diversity are key factors to understand and predict population dynamics. Marine species have been a difficult challenge in that respect, due to the difficulty in assessing population sizes and the open nature of such populations. Small, isolated islands with endemic species offer an opportunity to groundtruth population size estimates with empirical data and investigate the genetic consequences of such small populations. Here we focus on two endemic species of reef fish, the Clipperton damselfish, Stegastes baldwini, and the Clipperton angelfish, Holacanthus limbaughi, on Clipperton Atoll, tropical eastern Pacific. Visual surveys, performed over almost two decades and four expeditions, and genetic surveys based on genomic RAD sequences, allowed us to estimate kinship and genetic diversity, as well as to compare population size estimates based on visual surveys with effective population sizes based on genetics. We found that genetic and visual estimates of population numbers were remarkably similar. S. baldwini and H. limbaughi had population sizes of approximately 800,000 and 60,000, respectively. Relatively small population sizes resulted in low genetic diversity and the presence of apparent kinship. This study emphasizes the importance of small isolated islands as models to study population dynamics of marine organisms.

SLA-11 mutations are associated with litter size traits in Large White and Chinese DIV pigs

Litter size is an important economic traits in pigs. SLA-11 gene is a member of SLA (swine leukocyte antigen) complex. In our previous study, the SLA-11 gene was differentially expressed in PMSG-hCG stimulated preovulatory ovarian follicles of Chinese Taihu and Large White sows. Here, we identified two mutations (c.754-132 T > C and c.1421 + 38 T > C) in SLA-11 gene and analyzed the associations of two SNPs with litter size traits in Large White (n = 263) and DIV (n = 117) sows. The results showed that in Large White pigs, SLA-11 c.754-132 CC sows produced 0.74 and 0.87 more pigs per litter for TNB and NBA of all parities than did TT sows (p < .05); In DIV pigs, SLA-11 c.754-132 CC sows produced 1.17 more pigs per litter for TNB of all parities than did TC sows (p < .05). In Large White pigs, SLA-11 c.1421 + 38 CC sows produced 0.9 more pigs per litter for TNB of all parities than did TT sows (p < .05), while in DIV pigs SLA-11 c.1421 + 38 CC sows produced 0.84 and 0.7 less pigs per litter for TNB and NBA of all parities than did TT sows (p < .05). Our research indicated that SLA-11 mutations were potential molecular markers for improving the litter size traits in pigs.

Characterization of the Black Shank Pathogen, Phytophthora nicotianae, Across North Carolina Tobacco Production Areas

Black shank disease of tobacco, caused by the oomycete Phytophthora nicotianae, is a major threat to production in the United States and tobacco-producing areas worldwide. In a statewide survey of North Carolina, the rapid shift from race 0 to race 1 was documented. Collected pathogen isolates were characterized phenotypically for mating type and mefenoxam sensitivity, and genotypically by comparing sequences from three cytoplasmic and two nuclear regions. Both the A1 and A2 mating types were found throughout the state. When both mating types were recovered from the same field, pairings of isolates yielded viable oospores, indicating for the first time the potential for sexual sporulation by P. nicotianae in natural populations. Because the loss of complete resistance required a renewed use of the fungicide mefenoxam, a subset of the survey isolates was screened for sensitivity to the fungicide. All isolates were sensitive, with a mean effective concentration to inhibit 50% of hyphal growth of 0.4 μg/ml that was similar across mating types and races. Molecular characterization of 226 isolates revealed that the pathogen exists as multiple clonal types within the state. Genetic diversity among the pathogen population and the potential for sexual recombination may help explain the ability of the pathogen to rapidly adapt to host resistance genes.

Admixture, evolution, and variation in reproductive isolation in the Boechera puberula clade

BACKGROUND

Hybridization is very common in plants, and the incorporation of new alleles into existing lineages (i.e. admixture) can blur species boundaries. However, admixture also has the potential to increase standing genetic variation. With new sequencing methods, we can now study admixture and reproductive isolation at a much finer scale than in the past. The genus Boechera is an extraordinary example of admixture, with over 400 hybrid derivates of varying ploidy levels. Yet, few studies have assessed admixture in this genus on a genomic scale.

RESULTS

In this study, we used Genotyping-by-Sequencing (GBS) to clarify the evolution of the Boechera puberula clade, whose six members are scattered across the western United States. We further assessed patterns of admixture and reproductive isolation within the group, including two additional species (B. stricta and B. retrofracta) that are widespread across North America. Based on 14,815 common genetic variants, we found evidence for some cases of hybridization. We find evidence of both recent and more ancient admixture, and that levels of admixture vary across species.

CONCLUSIONS

We present evidence for a monophyletic origin of the B. puberula group, and a split of B. puberula into two subspecies. Further, when inferring reproductive isolation on the basis of presence and absence of admixture, we found that the accumulation of reproductive isolation between species does not seem to occur linearly with time since divergence in this system. We discuss our results in the context of sexuality and asexuality in Boechera.

Assessing the genetic diversity of the critically endangered Chinese sturgeon Acipenser sinensis using mitochondrial markers and genome-wide single-nucleotide polymorphisms from RAD-seq

As a living fossil, the endangered Chinese sturgeon (Acipenser sinensis) has been considered a national treasure in China. Here, the famous Gezhouba Dam and Three Gorges Dam on the Yangtze River were built in 1988 and 2006, for economic purposes. The natural population of Chinese sturgeon has declined since then, as these dams block its migratory route to the original spawning grounds in the middle reaches of the Yangtze River. In 2013 and 2014, there was an absence of spawning where it typically happened near the Gezhouba Dam. Nevertheless, from April to June in 2015, over 1,000 larvae with different body lengths (10-35 cm) were detected along the Shanghai Yangtze Estuary; but only little is currently known about the population genetic structure of the Chinese sturgeon. Herein, we inferred population genetic parameters from 462 available Chinese sturgeon specimens based on a 421-bp fragment of the mitochondrial DNA (mtDNA) D-loop region and 1,481,620 SNPs (single-nucleotide polymorphisms) generated by restriction site-associated DNA sequencing (RAD-seq). For the D-loop dataset, 15 haplotypes were determined. Randomly picked 23 individuals, representing the 15 D-loop haplotype groups, were subsequently used for further RAD-seq validation. The average nucleotide diversity calculated from the mtDNA and RAD datasets was 0.0086 and 0.000478, respectively. The overall effective female population size was calculated to be 1,255 to 2,607, and the long-term effective population size was estimated to range from 11,950 to 119,500. We observed that the genetic variability and the effective female population size of the current population in the Yangtze River are severely low, which are similar to the data reported over 10 years ago. The deduced relatively small effective population of female fish, limiting the genetic connectivity among Chinese sturgeon, should be considered a serious threat to this endangered species.

Revealing the distinct habitat ranges and hybrid zone of genetic sub-populations within Pseudo-nitzschia pungens (Bacillariophyceae) in the West Pacific area

Genetic sub-populations (clades) of cosmopolitan marine diatom Pseudo-nitzschia pungens might have distinct habitats, and their hybrid zone is suspected in higher latitude area of the West Pacific area, however, it is still unrevealed because of technical difficulties and lack of evidences in natural environments. The aim of this study is to investigate the habitat characteristics of each clade of P. pungens on geographical distribution with the habitat temperature ranges of each clade and to reveal their hybrid zone in the West Pacific area. We employed the 137 number of nucleotide sequences of P. pungens and its sampling data (spatial and temporal scale) originated from the West Pacific area, and used field application of qPCR assay for intra-specific level of P. pungens. Only two genotypes, clade I and III, were identified in the West Pacific area. Clade I was distributed from 39 to 32.3°N, and clade III were from 1.4 to 34.4°N. The estimated habitat temperature for the clade I and clade III ranges were 8.1-26.9 °C and 24.2-31.2 °C, respectively. The latitudinal distributions and temperature ranges of each clade were significantly different. The qPCR assay employed, and results suggested that the hybrid zone for clade I and III has been observed in the southern Korean coasts, and clade III might be introduced from the Southern Pacific area. The cell abundances of clade III were strongly related with the higher seawater temperature and warm current force. This study has defined distinct habitat characteristics of genetically different sub-populations of P. pungens, and revealed its hybrid zone in natural environment for the first time. We also provided strong evidences about dispersion of the population of clade III to higher latitude in the West Pacific area.

Genome-wide single nucleotide polymorphisms (SNPs) for a model invasive ascidian Botryllus schlosseri

Invasive species cause huge damages to ecology, environment and economy globally. The comprehensive understanding of invasion mechanisms, particularly genetic bases of micro-evolutionary processes responsible for invasion success, is essential for reducing potential damages caused by invasive species. The golden star tunicate, Botryllus schlosseri, has become a model species in invasion biology, mainly owing to its high invasiveness nature and small well-sequenced genome. However, the genome-wide genetic markers have not been well developed in this highly invasive species, thus limiting the comprehensive understanding of genetic mechanisms of invasion success. Using restriction site-associated DNA (RAD) tag sequencing, here we developed a high-quality resource of 14,119 out of 158,821 SNPs for B. schlosseri. These SNPs were relatively evenly distributed at each chromosome. SNP annotations showed that the majority of SNPs (63.20%) were located at intergenic regions, and 21.51% and 14.58% were located at introns and exons, respectively. In addition, the potential use of the developed SNPs for population genomics studies was primarily assessed, such as the estimate of observed heterozygosity (H _O ), expected heterozygosity (H _E ), nucleotide diversity (π), Wright's inbreeding coefficient (F _IS ) and effective population size (Ne). Our developed SNP resource would provide future studies the genome-wide genetic markers for genetic and genomic investigations, such as genetic bases of micro-evolutionary processes responsible for invasion success.

Double-digest RAD sequencing outperforms microsatellite loci at assigning paternity and estimating relatedness: A proof of concept in a highly promiscuous bird

Information on genetic relationships among individuals is essential to many studies of the behaviour and ecology of wild organisms. Parentage and relatedness assays based on large numbers of single nucleotide polymorphism (SNP) loci hold substantial advantages over the microsatellite markers traditionally used for these purposes. We present a double-digest restriction site-associated DNA sequencing (ddRAD-seq) analysis pipeline that, as such, simultaneously achieves the SNP discovery and genotyping steps and which is optimized to return a statistically powerful set of SNP markers (typically 150-600 after stringent filtering) from large numbers of individuals (up to 240 per run). We explore the trade-offs inherent in this approach through a set of experiments in a species with a complex social system, the variegated fairy-wren (Malurus lamberti) and further validate it in a phylogenetically broad set of other bird species. Through direct comparisons with a parallel data set from a robust panel of highly variable microsatellite markers, we show that this ddRAD-seq approach results in substantially improved power to discriminate among potential relatives and considerably more precise estimates of relatedness coefficients. The pipeline is designed to be universally applicable to all bird species (and with minor modifications to many other taxa), to be cost- and time-efficient, and to be replicable across independent runs such that genotype data from different study periods can be combined and analysed as field samples are accumulated.

Exploring the phylogeography of a hexaploid freshwater fish by RAD sequencing

The KwaZulu‐Natal yellowfish (Labeobarbus natalensis) is an abundant cyprinid, endemic to KwaZulu‐Natal Province, South Africa. In this study, we developed a single‐nucleotide polymorphism (SNP) dataset from double‐digest restriction site‐associated DNA (ddRAD) sequencing of samples across the distribution. We addressed several hidden challenges, primarily focusing on proper filtering of RAD data and selecting optimal parameters for data processing in polyploid lineages. We used the resulting high‐quality SNP dataset to investigate the population genetic structure of L. natalensis. A small number of mitochondrial markers present in these data had disproportionate influence on the recovered genetic structure. The presence of singleton SNPs also confounded genetic structure. We found a well‐supported division into northern and southern lineages, with further subdivision into five populations, one of which reflects north–south admixture. Approximate Bayesian Computation scenario testing supported a scenario where an ancestral population diverged into northern and southern lineages, which then diverged to yield the current five populations. All river systems showed similar levels of genetic diversity, which appears unrelated to drainage system size. Nucleotide diversity was highest in the smallest river system, the Mbokodweni, which, together with adjacent small coastal systems, should be considered as a key catchment for conservation.

Population Genetic Structure in Glyphosate-Resistant and -Susceptible Palmer Amaranth (Amaranthus palmeri) Populations Using Genotyping-by-sequencing (GBS)

Palmer amaranth (Amaranthus palmeri) is a major weed in United States cotton and soybean production systems. Originally native to the Southwest, the species has spread throughout the country. In 2004 a population of A. palmeri was identified with resistance to glyphosate, a herbicide heavily relied on in modern no-tillage and transgenic glyphosate-resistant (GR) crop systems. This project aims to determine the degree of genetic relatedness among eight different populations of GR and glyphosate-susceptible (GS) A. palmeri from various geographic regions in the United States by analyzing patterns of phylogeography and diversity to ascertain whether resistance evolved independently or spread from outside to an Arizona locality (AZ-R). Shikimic acid accumulation and EPSPS genomic copy assays confirmed resistance or susceptibility. With a set of 1,351 single nucleotide polymorphisms (SNPs), discovered by genotyping-by-sequencing (GBS), UPGMA phylogenetic analysis, principal component analysis, Bayesian model-based clustering, and pairwise comparisons of genetic distances were conducted. A GR population from Tennessee and two GS populations from Georgia and Arizona were identified as genetically distinct while the remaining GS populations from Kansas, Arizona, and Nebraska clustered together with two GR populations from Arizona and Georgia. Within the latter group, AZ-R was most closely related to the GS populations from Kansas and Arizona followed by the GR population from Georgia. GR populations from Georgia and Tennessee were genetically distinct from each other. No isolation by distance was detected and A. palmeri was revealed to be a species with high genetic diversity. The data suggest the following two possible scenarios: either glyphosate resistance was introduced to the Arizona locality from the east, or resistance evolved independently in Arizona. Glyphosate resistance in the Georgia and Tennessee localities most likely evolved separately. Thus, modern farmers need to continue to diversify weed management practices and prevent seed dispersal to mitigate herbicide resistance evolution in A. palmeri.

Efficiency of different strategies to mitigate ascertainment bias when using SNP panels in diversity studies

Background

Single nucleotide polymorphism (SNP) panels have been widely used to study genomic variations within and between populations. Methods of SNP discovery have been a matter of debate for their potential of introducing ascertainment bias, and genetic diversity results obtained from the SNP genotype data can be misleading. We used a total of 42 chicken populations where both individual genotyped array data and pool whole genome resequencing (WGS) data were available. We compared allele frequency distributions and genetic diversity measures (expected heterozygosity (H_e), fixation index (F_ST) values, genetic distances and principal components analysis (PCA)) between the two data types. With the array data, we applied different filtering options (SNPs polymorphic in samples of two Gallus gallus wild populations, linkage disequilibrium (LD) based pruning and minor allele frequency (MAF) filtering, and combinations thereof) to assess their potential to mitigate the ascertainment bias.

Results

Rare SNPs were underrepresented in the array data. Array data consistently overestimated H_e compared to WGS data, however, with a similar ranking of the breeds, as demonstrated by Spearman’s rank correlations ranging between 0.956 and 0.985. LD based pruning resulted in a reduced overestimation of H_e compared to the other filters and slightly improved the relationship with the WGS results. The raw array data and those with polymorphic SNPs in the wild samples underestimated pairwise F_ST values between breeds which had low F_ST (<0.15) in the WGS, and overestimated this parameter for high WGS F_ST (>0.15). LD based pruned data underestimated F_ST in a consistent manner. The genetic distance matrix from LD pruned data was more closely related to that of WGS than the other array versions. PCA was rather robust in all array versions, since the population structure on the PCA plot was generally well captured in comparison to the WGS data.

Conclusions

Among the tested filtering strategies, LD based pruning was found to account for the effects of ascertainment bias in the relatively best way, producing results which are most comparable to those obtained from WGS data and therefore is recommended for practical use.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12864-017-4416-9) contains supplementary material, which is available to authorized users.

Modern methods of assessing the taxonomic affiliation of honeybee colonies

At least 30 subspecies of the honeybee Apis mellifera L. were formed allopatrically during the evolution, which spreaded throughout all Africa, Europe and West Asia. The dark forest bee Apis mellifera mellifera is the only and most valuable subspecies for the Northern and Western Europe countries, adapted to productive living in the hard-continental climate of Eurasia. In the past 100 years, natural geographical isolation of subspecies has been disrupted as a result of a human activities. Mass transportations of honeybee colonies beyond the boundaries of their area have been threatened of loss the identity of gene pool of subspecies as a result of hybridization. Preservation of the gene pool of subspecies is possible only when controlling the transportation of honeybee colonies using the methods of identification of taxonomic affiliation of honeybee colonies. Now, dozens of methods have been developed to identify the taxonomic affiliation of honeybee's colony, which are based on the variability of body parts, allozyme loci, mitochondrial DNA loci, microsatellite nuclear loci, sites of single nucleotide polymorphism (SNP). The variability of microsatellite loci and the single nucleotide polymorphism sites have shown the greatest informativeness in identification of the taxonomic affiliation of honeybee's colony.

Reticulation, divergence, and the phylogeography-phylogenetics continuum

Phylogeography, and its extensions into comparative phylogeography, have their roots in the layering of gene trees across geography, a paradigm that was greatly facilitated by the nonrecombining, fast evolution provided by animal mtDNA. As phylogeography moves into the era of next-generation sequencing, the specter of reticulation at several levels-within loci and genomes in the form of recombination and across populations and species in the form of introgression-has raised its head with a prominence even greater than glimpsed during the nuclear gene PCR era. Here we explore the theme of reticulation in comparative phylogeography, speciation analysis, and phylogenomics, and ask how the centrality of gene trees has fared in the next-generation era. To frame these issues, we first provide a snapshot of multilocus phylogeographic studies across the Carpentarian Barrier, a prominent biogeographic barrier dividing faunas spanning the monsoon tropics in northern Australia. We find that divergence across this barrier is evident in most species, but is heterogeneous in time and demographic history, often reflecting the taxonomic distinctness of lineages spanning it. We then discuss a variety of forces generating reticulate patterns in phylogeography, including introgression, contact zones, and the potential selection-driven outliers on next-generation molecular markers. We emphasize the continued need for demographic models incorporating reticulation at the level of genomes and populations, and conclude that gene trees, whether explicit or implicit, should continue to play a role in the future of phylogeography.

Fast and cost-effective single nucleotide polymorphism (SNP) detection in the absence of a reference genome using semideep next-generation Random Amplicon Sequencing (RAMseq)

Biodiversity has suffered a dramatic global decline during the past decades, and monitoring tools are urgently needed providing data for the development and evaluation of conservation efforts both on a species and on a genetic level. However, in wild species, the assessment of genetic diversity is often hampered by the lack of suitable genetic markers. In this article, we present Random Amplicon Sequencing (RAMseq), a novel approach for fast and cost-effective detection of single nucleotide polymorphisms (SNPs) in nonmodel species by semideep sequencing of random amplicons. By applying RAMseq to the Eurasian otter (Lutra lutra), we identified 238 putative SNPs after quality filtering of all candidate loci and were able to validate 32 of 77 loci tested. In a second step, we evaluated the genotyping performance of these SNP loci in noninvasive samples, one of the most challenging genotyping applications, by comparing it with genotyping results of the same faecal samples at microsatellite markers. We compared (i) polymerase chain reaction (PCR) success rate, (ii) genotyping errors and (iii) Mendelian inheritance (population parameters). SNPs produced a significantly higher PCR success rate (75.5% vs. 65.1%) and lower mean allelic error rate (8.8% vs. 13.3%) than microsatellites, but showed a higher allelic dropout rate (29.7% vs. 19.8%). Genotyping results showed no deviations from Mendelian inheritance in any of the SNP loci. Hence, RAMseq appears to be a valuable tool for the detection of genetic markers in nonmodel species, which is a common challenge in conservation genetic studies.

Assessing SNP genotyping of noninvasively collected wildlife samples using microfluidic arrays

Noninvasively collected samples are a common source of DNA in wildlife genetic studies. Currently, single nucleotide polymorphism (SNP) genotyping using microfluidic arrays is emerging as an easy-to-use and cost-effective methodology. Here we assessed the performance of microfluidic SNP arrays in genotyping noninvasive samples from grey wolves, European wildcats and brown bears, and we compared results with traditional microsatellite genotyping. We successfully SNP-genotyped 87%, 80% and 97% of the wolf, cat and bear samples, respectively. Genotype recovery was higher based on SNPs, while both marker types identified the same individuals and provided almost identical estimates of pairwise differentiation. We found that samples for which all SNP loci were scored had no disagreements across the three replicates (except one locus in a wolf sample). Thus, we argue that call rate (amplification success) can be used as a proxy for genotype quality, allowing the reduction of replication effort when call rate is high. Furthermore, we used cycle threshold values of real-time PCR to guide the choice of protocols for SNP amplification. Finally, we provide general guidelines for successful SNP genotyping of degraded DNA using microfluidic technology.

Comparative assessment of SSR and SNP markers for inferring the population genetic structure of the common fungus Armillaria cepistipes

During the last years, simple sequence repeats (SSRs, also known as microsatellites) and single-nucleotide polymorphisms (SNPs) have become the most popular molecular markers for describing neutral genetic variation in populations of a wide range of organisms. However, only a limited number of studies has focused on comparing the performance of these two types of markers for describing the underlying genetic structure of wild populations. Moreover, none of these studies targeted fungi, the group of organisms with one of the most complex reproductive strategies. We evaluated the utility of SSRs and SNPs for inferring the neutral genetic structure of Armillaria cepistipes (basidiomycetes) at different spatial scales. For that, 407 samples were collected across a small (150 km²) area in the Ukrainian Carpathians and a large (41 000 km²) area in the Swiss Alps. All isolates were analyzed at 17 SSR loci distributed throughout the whole genome and at 24 SNP loci located in different single-copy conserved genes. The two markers showed different patterns of structure within the two spatial scales studied. The multi-allelic SSR markers seemed to be best suited for detecting genetic structure in indigenous fungal populations at a rather small spatial scale (radius of ~50-100 km). The pattern observed at SNP markers rather reflected ancient divergence of distant (~1000 km) populations that in addition are separated by mountain ranges. Despite these differences, both marker types were suitable for detecting the weak genetic structure of the two A. cepistipes populations investigated.

Novel efficient genome-wide SNP panels for the conservation of the highly endangered Iberian lynx

BACKGROUND

The Iberian lynx (Lynx pardinus) has been acknowledged as the most endangered felid species in the world. An intense contraction and fragmentation during the twentieth century left less than 100 individuals split in two isolated and genetically eroded populations by 2002. Genetic monitoring and management so far have been based on 36 STRs, but their limited variability and the more complex situation of current populations demand more efficient molecular markers. The recent characterization of the Iberian lynx genome identified more than 1.6 million SNPs, of which 1536 were selected and genotyped in an extended Iberian lynx sample.

METHODS

We validated 1492 SNPs and analysed their heterozygosity, Hardy-Weinberg equilibrium, and linkage disequilibrium. We then selected a panel of 343 minimally linked autosomal SNPs from which we extracted subsets optimized for four different typical tasks in conservation applications: individual identification, parentage assignment, relatedness estimation, and admixture classification, and compared their power to currently used STR panels.

RESULTS

We ascribed 21 SNPs to chromosome X based on their segregation patterns, and identified one additional marker that showed significant differentiation between sexes. For all applications considered, panels of autosomal SNPs showed higher power than the currently used STR set with only a very modest increase in the number of markers.

CONCLUSIONS

These novel panels of highly informative genome-wide SNPs provide more powerful, efficient, and flexible tools for the genetic management and non-invasive monitoring of Iberian lynx populations. This example highlights an important outcome of whole-genome studies in genetically threatened species.

Development of EPIC-PCR Markers for Lutjanus purpureus (Lutjanidae-Perciformes) and their Potential Applicability in Population Analyses

In the present study, a novel set of eight EPIC primers were developed for Lutjanus purpureus and assayed in five other marine teleosts including three lutjanids, one scianid and one anablepid. Most of the genomic regions used in this study presented genetic diversity indexes equal or greater than the intragenic regions commonly used in population genetics studies. Moreover, six out of eight markers showed cross-amplification with other taxa. Thus, the primers described here may be used to elucidate questions at the intraspecific level for a large number of taxa.

Evaluation of Bovine High-Density SNP Genotyping Array in Indigenous Dairy Cattle Breeds

In total 52 samples of Sahiwal ( 19 ), Tharparkar ( 17 ), and Gir ( 16 ) were genotyped by using BovineHD SNP chip to analyze minor allele frequency (MAF), genetic diversity, and linkage disequilibrium among these cattle. The common SNPs of BovineHD and 54K SNP Chips were also extracted and evaluated for their performance. Only 40%-50% SNPs of these arrays was found informative for genetic analysis in these cattle breeds. The overall mean of MAF for SNPs of BovineHD SNPChip was 0.248 ± 0.006, 0.241 ± 0.007, and 0.242 ± 0.009 in Sahiwal, Tharparkar and Gir, respectively, while that for 54K SNPs was on lower side. The average Reynold's genetic distance between breeds ranged from 0.042 to 0.055 based on BovineHD Beadchip, and from 0.052 to 0.084 based on 54K SNP Chip. The estimates of genetic diversity based on HD and 54K chips were almost same and, hence, low density chip seems to be good enough to decipher genetic diversity of these cattle breeds. The linkage disequilibrium started decaying (r² < 0.2) at 140 kb inter-marker distance and, hence, a 20K low density customized SNP array from HD chip could be designed for genomic selection in these cattle else the 54K Bead Chip as such will be useful.

B chromosome in Plantago lagopus Linnaeus, 1753 shows preferential transmission and accumulation through unusual processes

Plantago lagopus is a diploid (2n = 2x =12) weed belonging to family Plantaginaceae. We reported a novel B chromosome in this species composed of 5S and 45S ribosomal DNA and other repetitive elements. In the present work, presence of B chromosome(s) was confirmed through FISH on root tip and pollen mother cells. Several experiments were done to determine the transmission of B chromosome through male and female sex tracks. Progenies derived from the reciprocal crosses between plants with (1B) and without (0B) B chromosomes were studied. The frequency of B chromosome bearing plants was significantly higher than expected, in the progeny of 1B female × 0B male. Thus, the B chromosome seems to have preferential transmission through the female sex track, which may be due to meiotic drive. One of the most intriguing aspects of the present study was the recovery of plants having more chromosomes than the standard complement of 12 chromosomes. Such plants were isolated from the progenies of B chromosome carrying plants. The origin of these plants can be explained on the basis of a two step process; formation of unreduced gametes in 1B plants and fusion of unreduced gametes with the normal gametes or other unreduced gametes. Several molecular techniques were used which unequivocally confirmed similar genetic constitution of 1B (parent) and plants with higher number of chromosomes.

Karyotypic differentiation of populations of the common shrew Sorex araneus L. (Mammalia) in Belarus

The common shrews, Sorex araneus Linnaeus, 1758, inhabiting the territory of Belarus, are characterized by a significant variation in the frequency of Robertsonian (Rb) translocations. The frequency clines for translocations specific of three chromosome races: the West Dvina (gm, hk, ip, no, qr), Kiev (g/m, hi, k/o, n, p, q, r), and Białowieża (g/r, hn, ik, m/p, o, q) have already been studied in this territory. In this communication we report new data on polymorphic populations with Rb metacentrics specific of the Neroosa race (go, hi, kr, mn, p/q) in south-eastern Belarus, analyse the distribution of karyotypes in southern and central Belarus and draw particular attention to the fixation of the acrocentric variants of chromosomes in this area. The results show that certain Rb metacentrics specific of the Neroosa, West Dvina, Kiev, and Białowieża races (namely, go and pq; ip; ko; hn and ik, respectively) are absent in many polymorphic populations. Thus, the karyotypic differentiation of S. araneus in the studied area is determined by unequal spread of different Rb translocations and by fixation of acrocentric variants of specific chromosomes.

Null alleles are ubiquitous at microsatellite loci in the Wedge Clam (Donax trunculus)

Recent studies have reported an unusually high frequency of nonamplifying alleles at microsatellite loci in bivalves. Null alleles have been associated with heterozygous deficits in many studies. While several studies have tested for its presence using different analytical tools, few have empirically tested for its consequences in estimating population structure and differentiation. We characterised 16 newly developed microsatellite loci and show that null alleles are ubiquitous in the wedge clam, Donax trunculus. We carried out several tests to demonstrate that the large heterozygous deficits observed in the newly characterised loci were most likely due to null alleles. We tested the robustness of microsatellite genotyping for population assignment by showing that well-recognised biogeographic regions of the south Atlantic and south Mediterranean coast of Spain harbour genetically different populations.

Genomic Differentiation and Demographic Histories of Atlantic and Indo-Pacific Yellowfin Tuna (Thunnus albacares) Populations

Recent developments in the field of genomics have provided new and powerful insights into population structure and dynamics that are essential for the conservation of biological diversity. As a commercially highly valuable species, the yellowfin tuna (Thunnus albacares) is intensely exploited throughout its distribution in tropical oceans around the world, and is currently classified as near threatened. However, conservation efforts for this species have so far been hampered by limited knowledge of its population structure, due to incongruent results of previous investigations. Here, we use whole-genome sequencing in concert with a draft genome assembly to decipher the global population structure of the yellowfin tuna, and to investigate its demographic history. We detect significant differentiation of Atlantic and Indo-Pacific yellowfin tuna populations as well as the possibility of a third diverged yellowfin tuna group in the Arabian Sea. We further observe evidence for past population expansion as well as asymmetric gene flow from the Indo-Pacific to the Atlantic.

Genomic data detect corresponding signatures of population size change on an ecological time scale in two salamander species

Understanding the demography of species over recent history (e.g. <100 years) is critical in studies of ecology and evolution, but records of population history are rarely available. Surveying genetic variation is a potential alternative to census-based estimates of population size, and can yield insight into the demography of a population. However, to assess the performance of genetic methods, it is important to compare their estimates of population history to known demography. Here, we leveraged the exceptional resources from a wetland with 37 years of amphibian mark-recapture data to study the utility of genetically based demographic inference on salamander species with documented population declines (Ambystoma talpoideum) and expansions (A. opacum), patterns that have been shown to be correlated with changes in wetland hydroperiod. We generated ddRAD data from two temporally sampled populations of A. opacum (1993, 2013) and A. talpoideum (1984, 2011) and used coalescent-based demographic inference to compare alternate evolutionary models. For both species, demographic model inference supported population size changes that corroborated mark-recapture data. Parameter estimation in A. talpoideum was robust to our variations in analytical approach, while estimates for A. opacum were highly inconsistent, tempering our confidence in detecting a demographic trend in this species. Overall, our robust results in A. talpoideum suggest that genome-based demographic inference has utility on an ecological scale, but researchers should also be cognizant that these methods may not work in all systems and evolutionary scenarios. Demographic inference may be an important tool for population monitoring and conservation management planning.

SNPs selected by information content outperform randomly selected microsatellite loci for delineating genetic identification and introgression in the endangered dark European honeybee (Apis mellifera mellifera)

The honeybee (Apis mellifera) has been threatened by multiple factors including pests and pathogens, pesticides and loss of locally adapted gene complexes due to replacement and introgression. In western Europe, the genetic integrity of the native A. m. mellifera (M-lineage) is endangered due to trading and intensive queen breeding with commercial subspecies of eastern European ancestry (C-lineage). Effective conservation actions require reliable molecular tools to identify pure-bred A. m. mellifera colonies. Microsatellites have been preferred for identification of A. m. mellifera stocks across conservation centres. However, owing to high throughput, easy transferability between laboratories and low genotyping error, SNPs promise to become popular. Here, we compared the resolving power of a widely utilized microsatellite set to detect structure and introgression with that of different sets that combine a variable number of SNPs selected for their information content and genomic proximity to the microsatellite loci. Contrary to every SNP data set, microsatellites did not discriminate between the two lineages in the PCA space. Mean introgression proportions were identical across the two marker types, although at the individual level, microsatellites' performance was relatively poor at the upper range of Q-values, a result reflected by their lower precision. Our results suggest that SNPs are more accurate and powerful than microsatellites for identification of A. m. mellifera colonies, especially when they are selected by information content.

Variability in the sxt Gene Clusters of PSP Toxin Producing Aphanizomenon gracile Strains from Norway, Spain, Germany and North America

Paralytic shellfish poisoning (PSP) toxin production has been detected worldwide in the cyanobacterial genera Anabaena, Lyngbya, Scytonema, Cuspidothrix and Aphanizomenon. In Europe Aphanizomenon gracile and Cuspidothrix issatschenkoi are the only known producers of PSP toxins and are found in Southwest and Central European freshwater bodies. In this study the PSP toxin producing Aphanizomenon sp. strain NIVA-CYA 851 was isolated from the Norwegian Lake Hillestadvannet. In a polyphasic approach NIVA-CYA 851 was morphologically and phylogenetically classified, and investigated for toxin production. The strain NIVA-CYA 851 was identified as A. gracile using 16S rRNA gene phylogeny and was confirmed to produce neosaxitoxin, saxitoxin and gonyautoxin 5 by LC-MS. The whole sxt gene clusters (circa 27.3 kb) of four A. gracile strains: NIVA-CYA 851 (Norway); NIVA-CYA 655 & NIVA-CYA 676 (Germany); and UAM 529 (Spain), all from latitudes between 40° and 59° North were sequenced and compared with the sxt gene cluster of reference strain A. gracile NH-5 from the USA. All five sxt gene clusters are highly conserved with similarities exceeding 99.4%, but they differ slightly in the number and presence of single nucleotide polymorphisms (SNPs) and insertions/deletions (In/Dels). Altogether 178 variable sites (44 SNPs and 4 In/Dels, comprising 134 nucleotides) were found in the sxt gene clusters of the Norwegian, German and Spanish strains compared to the reference strain. Thirty-nine SNPs were located in 16 of the 27 coding regions. The sxt gene clusters of NIVA-CYA 851, NIVA-CYA 655, NIVA-CYA 676 and UAM 529, were characterized by 15, 16, 19 and 23 SNPs respectively. Only the Norwegian strain NIVA-CYA 851 possessed an insertion of 126 base pairs (bp) in the noncoding area between the sxtA and sxtE genes and a deletion of 6 nucleotides in the sxtN gene. The sxtI gene showed the highest variability and is recommended as the best genetic marker for further phylogenetic studies of the sxt gene cluster of A. gracile.

This study confirms for the first time the role of A. gracile as a PSP toxin producer in Norwegian waters, representing the northernmost occurrence of PSP toxin producing A. gracile in Europe known so far.

Global Population Structure of a Worldwide Pest and Virus Vector: Genetic Diversity and Population History of the Bemisia tabaci Sibling Species Group

The whitefly Bemisia tabaci sibling species (sibsp.) group comprises morphologically indiscernible lineages of well-known exemplars referred to as biotypes. It is distributed throughout tropical and subtropical latitudes and includes the contemporary invasive haplotypes, termed B and Q. Several well-studied B. tabaci biotypes exhibit ecological and biological diversity, however, most members are poorly studied or completely uncharacterized. Genetic studies have revealed substantial diversity within the group based on a fragment of the mitochondrial cytochrome oxidase I (mtCOI) sequence (haplotypes), with other tested markers being less useful for deep phylogenetic comparisons. The view of global relationships within the B. tabaci sibsp. group is largely derived from this single marker, making assessment of gene flow and genetic structure difficult at the population level. Here, the population structure was explored for B. tabaci in a global context using nuclear data from variable microsatellite markers. Worldwide collections were examined representing most of the available diversity, including known monophagous, polyphagous, invasive, and indigenous haplotypes. Well-characterized biotypes and other related geographic lineages discovered represented highly differentiated genetic clusters with little or no evidence of gene flow. The invasive B and Q biotypes exhibited moderate to high levels of genetic diversity, suggesting that they stemmed from large founding populations that have maintained ancestral variation, despite homogenizing effects, possibly due to human-mediated among-population gene flow. Results of the microsatellite analyses are in general agreement with published mtCOI phylogenies; however, notable conflicts exist between the nuclear and mitochondrial relationships, highlighting the need for a multifaceted approach to delineate the evolutionary history of the group. This study supports the hypothesis that the extant B. tabaci sibsp. group contains ancient genetic entities and highlights the vast cryptic diversity throughout the genome in the group.

Sex beyond species: the first genetically analyzed case of intergeneric fertile hybridization in pinnipeds

A species, according to the biological concept, is a natural group of potentially interbreeding individuals isolated by diverse mechanisms. Hybridization is considered the production of offspring resulting from the interbreeding of two genetically distinct taxa. It has been documented in over 10% of wild animals, and at least in 34 cases for Artic marine mammals. In Otariids, intergeneric hybridization has been reported though neither confirming it through genetic analyses nor presenting evidence of fertile offspring. In this study, we report the finding of a hybrid adult female between a South American fur seal (Arctocephalus australis) and a South American sea lion (Otaria byronia), and its offspring, a male pup, in Uruguay. Further based on morphological constraints and breeding seasons, sex-biased hybridization between the two species is hypothesized. Morphological and genetic (nuclear and mitochondrial) results confirm de hybrid nature of the female-pup pair. Here we discuss a genetic dilution effect, considering other hybridization events must be occurring, and how isolation mechanisms could be circumvented. Moreover, the results obtained from stable isotope analysis suggest feeding habits may be a trait transmitted maternally, leading to consider broader issues regarding hybridization as an evolutionary innovation phenomenon.