GCalignR: An R package for aligning gas-chromatography data for ecological and evolutionary studies.

• Published in: PLOS ONE
• Main Authors: Ottensmann, M, Stoffel, MA, Nichols, HJ, Hoffman, JI
• Format: Article in Journal/Newspaper
• Language: English
• Published: Public Library of Science 2018
• Subjects: QD Chemistry; QH Natural history; QL Zoology; Antarctic; Fid; Antarc*; Antarctic Fur Seals
• Online Access: http://researchonline.ljmu.ac.uk/id/eprint/8864/; https://researchonline.ljmu.ac.uk/id/eprint/8864/1/Ottensmann%20et%20al%202018%20GCalignR.pdf; https://doi.org/10.1371/journal.pone.0198311

Chemical cues are arguably the most fundamental means of animal communication and play an important role in mate choice and kin recognition. Consequently, there is growing interest in the use of gas chromatography (GC) to investigate the chemical basis of eco-evolutionary interactions. Both GC-MS (mass spectrometry) and FID (flame ionization detection) are commonly used to characterise the chemical composition of biological samples such as skin swabs. The resulting chromatograms comprise peaks that are separated according to their retention times and which represent different substances. Across chromatograms of different samples, homologous substances are expected to elute at similar retention times. However, random and often unavoidable experimental variation introduces noise, making the alignment of homologous peaks challenging, particularly with GC-FID data where mass spectral data are lacking. Here we present GCalignR, a user-friendly R package for aligning GC-FID data based on retention times. The package was developed specifically for ecological and evolutionary studies that seek to investigate similarity patterns across multiple and often highly variable biological samples, for example representing different sexes, age classes or reproductive stages. The package also implements dynamic visualisations to facilitate inspection and fine-tuning of the resulting alignments and can be integrated within a broader workflow in R to facilitate downstream multivariate analyses. We demonstrate an example workflow using empirical data from Antarctic fur seals and explore the impact of user-defined parameter values by calculating alignment error rates for multiple datasets. The resulting alignments had low error rates for most of the explored parameter space and we could also show that GCalignR performed equally well or better than other available software. We hope that GCalignR will help to simplify the processing of chemical datasets and improve the standardization and reproducibility of chemical analyses in studies of animal chemical communication and related fields.