# Create cladogram/phylogram from newick via xelatex/tikz/forest

## A picture is worth a thousand words

![template.png](images/template.png)

## Why not FigTree/Dendroscope/MEGA?

Full control over the trees: fonts, colors, line widths, annotations, and more. And, the process can
be recorded.

Below is the resulting file opened in Adobe Illustrator.

This is a very clean vector graphics in which all text would be editable.

![clean graphics](images/clean-graphics.png)

## Start from scratch

* Let's start by manually entering some basal taxa of animals

```shell
cat <<'EOF' > newick/example.nwk
(Ctenophora,(Porifera,(Placozoa,(Bilateria,Cnidaria))))Animalia;
EOF

nwr indent --text '.   ' newick/example.nwk

```

* You can verify it with your eyes

```text
(
.   Ctenophora,
.   (
.   .   Porifera,
.   .   (
.   .   .   Placozoa,
.   .   .   (
.   .   .   .   Bilateria,
.   .   .   .   Cnidaria
.   .   .   )
.   .   )
.   )
)Animalia;
```

* Add more information to comments
    * One node at a time

```shell
cat newick/example.nwk |
    nwr comment stdin --node Ctenophora --comment 192  |
    nwr comment stdin --node Porifera --color green --comment 8579 |
    nwr comment stdin --node Placozoa --comment 1 |
    nwr comment stdin --node Cnidaria --comment 13138 |
    nwr comment stdin --lca Bilateria,Cnidaria --label Planulozoa --dot red --rec LemonChiffon |
    nwr comment stdin --lca Bilateria,Placozoa --label Parahoxozoa --color green --dot red |
    nwr comment stdin --lca Bilateria,Porifera --label NOTE --bar purple |
    nwr indent stdin |
    tee newick/example.comment.nwk
#(
#  Ctenophora[comment=192],
#  (
#    Porifera[color=green:comment=8579],
#    (
#      Placozoa[comment=1],
#      (
#        Bilateria,
#        Cnidaria[comment=13138]
#      )[label=Planulozoa:dot=red]
#    )[label=Parahoxozoa:color=green:dot=red]
#  )[label=NOTE:bar=purple]
#)Animalia;

```

* Newick file is converted to `forest` codes

```shell
nwr tex newick/example.comment.nwk --bare
#[, dot, label={Animalia}, tier=4,
#  [{Ctenophora}, comment={192}, tier=0,]
#  [, bar={purple}, label={NOTE}, tier=3,
#    [\color{green}{Porifera}, comment={8579}, tier=0,]
#    [, dot={red}, label=\color{green}{Parahoxozoa}, tier=2,
#      [{Placozoa}, comment={1}, tier=0,]
#      [, dot={red}, rec={LemonChiffon}, label={Planulozoa}, tier=1,
#        [{Bilateria}, tier=0,]
#        [{Cnidaria}, comment={13138}, tier=0,]
#      ]
#    ]
#  ]
#]

```

* Produce pdf
    * Edit the .tex file as you wish

```shell
nwr tex newick/example.comment.nwk -s -o tex/example.tex

tectonic tex/example.tex --outdir pdf

```

## From a Newick file

### hg38.30way

```shell
curl -L https://hgdownload.cse.ucsc.edu/goldenpath/hg38/multiz30way/hg38.30way.scientificNames.nh \
    > newick/hg38.30way.nwk

cat newick/hg38.30way.nwk |
    nwr comment stdin --lca Homo_sapiens,Nomascus_leucogenys --label Hominoidea --dot --rec LemonChiffon |
    nwr comment stdin --lca Macaca_mulatta,Rhinopithecus_roxellana --label Cercopithecidae  --dot --rec TeaRose |
    nwr comment stdin --lca Callithrix_jacchus,Aotus_nancymaae --label Cebidae --dot --rec Celadon |
    nwr tex stdin --bl -s -o tex/hg38.30way.tex

tectonic tex/hg38.30way.tex --outdir pdf

```

### Opisthokonta

Create `Opisthokonta.nwk` manually

```shell
cat newick/Opisthokonta.nwk |
    nwr comment stdin -n Fungi -n Metazoa --color red |
    nwr comment stdin -n Holomycota --rec TeaRose |
    nwr comment stdin -n Choanozoa --rec ElectricBlue |
    nwr tex stdin -s -o tex/Opisthokonta.tex

tectonic tex/Opisthokonta.tex --outdir pdf

```

## From taxonomy

### Algae

```shell
cd ~/Scripts/nwr/tree/

nwr common \
    "Cyanobacteria" \
    "Euglenida" \
    "Kinetoplastea" \
    "Dinophyceae" \
    "Apicomplexa" \
    "Ciliophora" \
    "Haptophyta" \
    "Cryptophyceae" \
    "Chrysophyceae" \
    "Bacillariophyta" \
    "Rhodophyta" \
    "Chlorophyta" \
    "Phaeophyceae" |
    sed 's/cellular organisms//g' |
    sed 's/\broot\b//g' |
    nwr tex stdin -s -o tex/Algae.tex

tectonic tex/Algae.tex --outdir pdf

```

### Plants

```shell
cd ~/Scripts/nwr/tree/

nwr common \
    "Anthocerotophyta" \
    "Bryophyta" \
    "Marchantiophyta" \
    \
    "Tracheophyta" \
    \
    "Lycopodiopsida" \
    "Euphyllophyta" \
    "Polypodiopsida" \
    \
    "Acrogymnospermae" \
    "Cycadopsida" \
    "Ginkgoopsida" \
    "Gnetopsida" \
    "Pinopsida" \
    \
    "Spermatophyta" \
    "Amborellales" \
    "Nymphaeales" \
    "Austrobaileyales" \
    "Mesangiospermae" \
    "Magnoliidae" \
    "Liliopsida" \
    "eudicotyledons " \
    "asterids" \
    "rosids" |
    sed 's/\broot\b//g' |
    nwr tex stdin -s -o tex/Plants.tex

tectonic tex/Plants.tex --outdir pdf

```

## From a `forest` file

### animals

```shell
# animals-simple
nwr tex forest/animals-simple.forest --forest -s -o tex/animals-simple.tex

tectonic tex/animals-simple.tex --outdir pdf

# animals-simple.trans
cat tex/animals-simple.tex |
    sed -f translation.sed \
    > tex/animals-simple.trans.tex

tectonic tex/animals-simple.trans.tex --outdir pdf

```

### chordates

```shell
nwr tex forest/chordates.forest --forest -s -o tex/chordates.tex

cat tex/chordates.tex |
    sed -f translation.sed \
    > tex/chordates.trans.tex

tectonic tex/chordates.trans.tex --outdir pdf

```

## Taxonomy

## Translation

* Create `translation.tsv` manually

* Convert .tsv to sed scripts

```shell
# sed scripts
cat translation.tsv |
    grep -v "^#" |
    grep "." |
    perl -nla -F"\t" -e '
        print q(s/{) . $F[0] . q(}/{) . $F[1] . q(}/g;);
    ' \
    > translation.sed

cat translation.tsv |
    grep -v "^#" |
    grep "." |
    perl -nla -F"\t" -e '
        print q(s/{) . $F[0] . q(}/{) . $F[0] . q( \\\\footnotesize{) . $F[1] . q(}}/g;);
    ' \
    > translation.append.sed

```

```shell
# replaced with translations
cat tex/example.tex |
    sed -f translation.sed \
    > tex/example.trans.tex

tectonic tex/example.trans.tex --outdir pdf

# append translations
cat tex/example.tex |
    sed -f translation.append.sed \
    > tex/example.trans.append.tex

tectonic tex/example.trans.append.tex --outdir pdf

```