ClinePlotR: Visualizing genomic clines and detecting outliers in R

Patterns of multi-locus differentiation (i.e., genomic clines) often extend broadly across hybrid zones and their quantification can help diagnose how species boundaries are shaped by adaptive processes, both intrinsic and extrinsic. In this sense, the transitioning of loci across admixed individuals can be contrasted as a function of the genome-wide trend, in turn allowing an expansion of clinal theory across a much wider array of biodiversity. However, computational tools that serve to interpret and consequently visualize ‘genomic clines’ are limited. Here, we introduce the ClinePlotR R-package for visualizing genomic clines and detecting outlier loci using output generated by two popular software packages, bgc and Introgress. ClinePlotR bundles both input generation (i.e, filtering datasets and creating specialized file formats) and output processing (e.g., MCMC thinning and burn-in) with functions that directly facilitate interpretation and hypothesis testing. Tools are also provided for post-hoc analyses that interface with external packages such as ENMeval and RIdeogram. Our package increases the reproducibility and accessibility of genomic cline methods, thus allowing an expanded user base and promoting these methods as mechanisms to address diverse evolutionary questions in both model and non-model organisms.

[1]  Z. Gompert,et al.  introgress: a software package for mapping components of isolation in hybrids , 2010, Molecular ecology resources.

[2]  J. Mallet,et al.  Contrasting genomic and phenotypic outcomes of hybridization between pairs of mimetic butterfly taxa across a suture zone , 2020, Molecular ecology.

[3]  Z. Gompert,et al.  Bayesian estimation of genomic clines , 2011, Molecular ecology.

[4]  C. Buerkle,et al.  Admixture as the basis for genetic mapping. , 2008, Trends in ecology & evolution.

[5]  Robert A. Boria,et al.  ENMeval: An R package for conducting spatially independent evaluations and estimating optimal model complexity for Maxent ecological niche models , 2014 .

[6]  Isaac Overcast,et al.  ipyrad: Interactive assembly and analysis of RADseq datasets , 2020, Bioinform..

[7]  T. Papenfuss,et al.  An Annotated Chromosome-Level Reference Genome of the Red-Eared Slider Turtle (Trachemys scripta elegans) , 2020, Genome biology and evolution.

[8]  E. Mandeville,et al.  Analysis of Population Genomic Data from Hybrid Zones , 2017 .

[9]  Z. Gompert,et al.  A powerful regression‐based method for admixture mapping of isolation across the genome of hybrids , 2009, Molecular ecology.

[10]  Z. Gompert,et al.  bgc: Software for Bayesian estimation of genomic clines , 2012, Molecular ecology resources.

[11]  Robert P. Anderson,et al.  Maximum entropy modeling of species geographic distributions , 2006 .

[12]  Tyler K. Chafin,et al.  Contrasting signatures of introgression in North American box turtle (Terrapene spp.) contact zones , 2019, bioRxiv.

[13]  B. Fitzpatrick,et al.  Alternative forms for genomic clines , 2013, Ecology and evolution.

[14]  J. Fordyce,et al.  GENOMIC REGIONS WITH A HISTORY OF DIVERGENT SELECTION AFFECT FITNESS OF HYBRIDS BETWEEN TWO BUTTERFLY SPECIES , 2012, Evolution; international journal of organic evolution.

[15]  M. Braun,et al.  The genomic consequences of adaptive divergence and reproductive isolation between species of manakins , 2013, Molecular ecology.

[16]  Z. Gompert,et al.  Genomics of isolation in hybrids , 2012, Philosophical Transactions of the Royal Society B: Biological Sciences.

[17]  Jinhui Chen,et al.  RIdeogram: drawing SVG graphics to visualize and map genome-wide data on the idiograms , 2020, PeerJ Comput. Sci..

[18]  P. David,et al.  The geography of introgression in a patchy environment and the thorn in the side of ecological speciation , 2013 .

[19]  Heng Li,et al.  Minimap2: pairwise alignment for nucleotide sequences , 2017, Bioinform..

[20]  Tyler K. Chafin,et al.  Hybridization drives genetic erosion in sympatric desert fishes of western North America , 2019, Heredity.