-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathancistr.qmd
53 lines (43 loc) · 2.64 KB
/
ancistr.qmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
---
title: "Genetic identification of the common bristlenose"
author: "Rupert A. Collins"
theme: "cosmo"
format:
html:
self-contained: true
code-fold: true
fig-width: 11
fig-height: 12
fig-cap-location: "top"
fig-format: "svg"
fig-align: "left"
tbl-align: "left"
---
.jpg){#fig-sp3}
### INTRODUCTION
The true identity of the wild ancestor(s) of [_Ancistrus_ sp. (3)](https://www.planetcatfish.com/ancistrus_cf_cirrhosus), AKA the common bristlenose catfish (@fig-sp3), is not known. Here we use phylogenetic techniques to discover the closest potential relatives of this common aquarium species. The analysis uses maternally inherited mitochondrial DNA sequences publicly available on the [NCBI GenBank](https://www.ncbi.nlm.nih.gov/nuccore) nucleotide database. Sequence data for the cytochrome _c_ oxidase I (COI) "DNA barcode" gene were downloaded from GenBank and the data cleaned and duplicates removed. A phylogenetic tree was made using maximum likelihood, and genetic distances (sequence similarity) were calculated. An individual of the common bristlenose was obtained from an aquarium store in New Zealand in October 2009 to act as a reference for the traded species.
### RESULTS
```{r run.code}
#| eval: true
#| echo: true
#| output: false
# https://github.com/boopsboops/ancistr/blob/main/scripts/ancistr.R
source(here::here("scripts/ancistr.R"))
```
The analysis was conducted on `r today.long` using GenBank v`r gb.version`. A total of `r n.ancistrus` described _Ancistrus_ species were obtained, plus multiple undescribed or unidentified _Ancistrus_ spp. The New Zealand aquarium specimen of _Ancistrus_ sp. (3) can be seen highlighted in red in @fig-tree, clustered with its closest relatives. Sequence similarity (genetic distance) of the top ten closest species is reported in @tbl-dist.
```{r}
#| eval: true
#| echo: false
#| output: true
#| label: fig-tree
#| fig-cap: "Maximum likelihood phylogenetic tree of _Ancistrus_ COI sequence data mined from NCBI GenBank. Tree was made using the R package [Phangorn](https://klausvigo.github.io/phangorn/index.html), and was rooted on _Lasiancistrus_ (subsequently removed)."
print(p)
```
```{r}
#| eval: true
#| echo: false
#| output: true
#| label: tbl-dist
#| tbl-cap: "Uncorrected genetic distances (sequence similarity) of the top ten closest species in the analysis to the reference sequence of _Ancistrus_ sp. (3) (New Zealand trade). Accession codes can be searched for more information at [NCBI GenBank](https://www.ncbi.nlm.nih.gov/nuccore)."
knitr::kable(ancistrus.aligned.trimmed.dist.rc.filt.tidy,digits=1)
```