Genomic evidence for the parallel regression of melatonin synthesis and signaling pathways in placental mammals

Supplementary Material for:

Emerling C.A., Springer M.S., Gatesy J., Jones Z., Hamilton D., Xia-Zhu D., Collin M.A., and Delsuc F. (2021). Genomic evidence for the parallel regression of melatonin synthesis and signaling pathways in placental mammals. Open Research Europe 1:75. doi:10.12688/openreseurope.13795.1.

Supplementary File Legends:

- Supplementary_Figure_S1.pdf: AANAT PAML ‘master model’ showing branch categories, corresponding to “Model 1: 24 ratio” in Supplementary Table S7.

- Supplementary_Figure_S2.pdf: ASMT PAML ‘master model’ showing branch categories, corresponding to “Model 2: 24 ratio” in Supplementary Table S8.

- Supplementary_Figure_S3.pdf: MTNR1A PAML ‘master model’ showing branch categories, corresponding to “Model 1: 27 ratio” in Supplementary Table S9.

- Supplementary_Figure_S4.pdf: MTNR1B PAML ‘master model’ showing branch categories, corresponding to “Model 1: 46 ratio” in Supplementary Table S10.

- Supplementary_Figure_S5.pdf: RAxML AANAT gene tree. Numbers at nodes correspond to bootstrap support values.

- Supplementary_Figure_S6.pdf: RAxML ASMT gene tree. Numbers at nodes correspond to bootstrap support values.

- Supplementary_Figure_S7.pdf: RAxML MTNR1A+MTNR1B tree. Numbers at nodes correspond to bootstrap support values.

- Supplementary_Figure_S8.pdf: Supporting data showing the inactivation of MTNR1A exon 2 in cetaceans. Read Supplementary Table S13 for further details.

- Supplementary_Figure_S9.pdf: Supporting data showing the inactivation of ASMT in spalacids and Fukomys damarensis. Read Supplementary Table S13 for further details.

- Supplementary_Figure_S10.pdf: Supporting data showing the inactivation of MTNR1A in hyracoids and Cyclopes didactylus. Read Supplementary Table S13 for further details.

- Supplementary_Figure_S11.pdf: Supporting data showing the inactivation of MTNR1A in sirenians. Read Supplementary Table S13 for further details.

- Supplementary_Figure_S12.pdf: Supporting data showing the inactivation of AANAT in sirenians and a polymorphic premature stop codon in exon 5 of ASMT in Trichechus manatus. Read Supplementary Table S13 for further details.

- Supplementary_Figure_S13.pdf: Supporting data showing the inactivation of MTNR1A in Condylura cristata. Read Supplementary Table S13 for further details.

- Supplementary_Figure_S14.pdf: Supporting data showing the inactivation of MTNR1A in Phataginus tricuspis. Read Supplementary Table S14 for further details.

- Supplementary_Figure_S15.pdf: PAML AANAT results, Model 1: 24 ratio (see Supplementary Table S7).

- Supplementary_Figure_S16.pdf: PAML ASMT results, Model 2: 24 ratio (see Supplementary Table S8).

- Supplementary_Figure_S17.pdf: PAML MTNR1A results, Model 1: 27 ratio (see Supplementary Table S9).

- Supplementary_Figure_S18.pdf: PAML MTNR1B results, Model 1: 46 ratio (see Supplementary Table S10).

- Supplementary_Table_S1.xlsx: List of species examined in this study and the sources of the genes. Source key: WGS: Sequences derived from NCBI's Whole Genome Shotgun database, with accession prefix provided; Whole Genome Sequencing of Short Reads: whole genomes were sequenced using short-read technologies. The methodologies varied for the species, and will be or have been published with other projects, so please contact the author(s) for information on the specific methodology and samples used (Xenarthrans, Proteles cristatus, Otocyon megalotis: Frédéric Delsuc, e-mail: [e-mail]; Crocodylians: John Gatesy, e-mail: [e-mail]; Dugong dugon: Mark Springer, e-mail: [e-mail]; SRA: sequences derived from NCBI's Sequence Read Archive; GenBank: sequences derived from NCBI's nucleotide collection; Bowhead Whale Genome Resource: sequences derived from http://www.bowhead-whale.org; Ensembl: sequences derived from Ensembl genome browser (www.ensembl.org)l; Discovar de novo: sequences derived genomes assembled via Discovar de novo  (https://software.broadinstitute.org/software/discovar/blog/). Coverage: indicates coverage of the whole genome (reported in NCBI or other source) or individual genes (derived from short read mapping). Scaffold and contig N50: reported in NCBI or other source.

- Supplementary_Table_S2.xlsx: Accession numbers and functionality of AANAT in species examined. If Accession # indicated as “New”, sequence generated for this study and can be found in Supplementary Dataset S1. Parentheses after accession number indicates coordinates for sequence on the contig / scaffold. Exon colors code for the following: green = putatively functional; yellow = missing (e.g., negative BLAST results, negative mapping results); pink = one or more inactivating mutations found. Abbreviations for mutations are as follows: del = deletion; ins = insertion; start = start codon mutation; stop = premature stop codon; ? = ambiguity whether the mutation is shared among all members of the clade. Abbreviations in brackets following an inactivating mutation indicate shared inactivating mutation. Key for each abbreviation follows: Bacu = Balaenoptera acutorostrata; BALA = Balaenidae; BALAEN = Balaenopteridae; Bbon = Balaenoptera bonaerensis; CAB = Cabassous; Ccap = Cebus capucinus; CETA = Cetacea; CHLAM = Chlamyphoridae; CHOL = Choloepus; Cjac = Callithrix jacchus; CING = Cingulata; DASY = Dasypodidae; DELP = Delphinidae; DERM = Dermoptera; Erob = Eschrichtius robustus; INIA = Inia; FOLI = Folivora; GALE = Galeopterus; LIPO = Lipotes; Lobl = Lagenorhynchus obliquidens; MANI = Manidae; MONO = Monodontidae; MYRM = Myrmecophagidae; MYST = Mysticeti; NPP = Not present in Platanista or Physeteroidea, but present in other Odontocetes; NPZ = Not present in Ziphiidae, but present in other Odontocetes; Oorc = Orcinus orca; PEUT = Tolypeutinae; PHOC = Phocoenidae; PHOL = Pholidota; PHOR = Chlamyphorinae; PILO = Pilosa; PHYS = Physeteroidea; PONT = Pontoporia; Schi = Sousa chinensis; SIRE = Sirenia; Tadu = Tursiops aduncus; TOLY = Tolypeutes; VERM = Vermilingua; XEN = Xenarthra.

- Supplementary_Table_S3.xlsx: Accession numbers and functionality of ASMT in species examined. See Table S2 caption for details.

- Supplementary_Table_S4.xlsx: Accession numbers and functionality of MTNR1A in species examined. See Table S2 caption for details.

- Supplementary_Table_S5.xlsx: Accession numbers and functionality of MTNR1B in species examined. See Table S2 caption for details.

- Supplementary_Table_S6.xlsx: Codon frequency model selection. These are the results from one ratio dN/dS analyses using different codon frequency models. AIC = Akaike Information Criterion.

- Supplementary_Table_S7.xlsx: Results of AANAT PAML dN/dS analyses for mammals. Model: BG = branch(es) grouped with background; fixed 1 = branch(es) fixed at 1. p’-value: p-value after Holm-Bonferroni correction for multiple testing. Model Comparison: if model comparison yields statistically significant differences (p < 0.05), model comparison bolded and given green background; if model comparison is still significant after Holm-Bonferroni correction, asterisk (*) added. For most models, w only shown for branch(es) of interest. Numbers in front of taxonomic names in first row correspond to numbers in the master model shown in Supplementary Figure S1.

- Supplementary_Table_S8.xlsx: Results of ASMT PAML dN/dS analyses for mammals. Refer to Table S7 caption for additional details. Numbers in front of taxonomic names in first row correspond to numbers in the master model shown in Supplementary Figure S2.

- Supplementary_Table_S9.xlsx: Results of MTNR1A PAML dN/dS analyses for mammals. Refer to Table S7 caption for additional details. Numbers in front of taxonomic names in first row correspond to numbers in the master model shown in Supplementary Figure S3.

- Supplementary_Table_S10.xlsx: Results of MTNR1B PAML dN/dS analyses for mammals. Refer to Table S7 caption for additional details. Numbers in front of taxonomic names in first row correspond to numbers in the master model shown in Supplementary Figure S4.

- Supplementary_Table_S11.xlsx: Results of PAML analyses for sauropsids.

- Supplementary_Table_S12.xlsx: Results of BLASTing and mapping short reads from Alligator mississippiensis RNA sequencing experiments.

- Supplementary_Table_S13.xlsx: Supporting data for validating putative inactivating mutations. Validating data came from four general sources of information: mutations shared by more than one species within a clade, mutations shared by two sources of sequencing data for the same species, mutations validated by coverage of mapped short reads and statistically elevated dN/dS ratio estimates. For additional details, see Supplementary Tables S2–S5 and S7–S10, as well as Figure 2 and Supplementary Figures S8–S18.

- Supplementary_Dataset_S1.txt: Genomic alignments in fasta format used to determine the pseudogene/functional status of all four melatonin genes in different taxonomic groups.

- Supplementary_Dataset_S2.txt: Alignment of AANAT in phylip format used in maximum likelihood phylogenetic reconstruction with RAxML. 

- Supplementary_Dataset_S3.txt: Alignment of ASMT in phylip format used in maximum likelihood phylogenetic reconstruction with RAxML. 

- Supplementary_Dataset_S4.txt: Alignment of MTNR1A and MTNR1B in phylip format used in maximum likelihood phylogenetic reconstruction with RAxML. 

- Supplementary_Dataset_S5.txt: Codon alignments of AANAT used in selection pressure analyses with PAML. 

- Supplementary_Dataset_S6.txt: Codon alignments of ASMT used in selection pressure analyses with PAML.

- Supplementary_Dataset_S7.txt: Codon alignments of MTNR1A used in selection pressure analyses with PAML.

- Supplementary_Dataset_S8.txt: Codon alignments of MTNR1B used in selection pressure analyses with PAML.

- Supplementary_Dataset_S9.txt: Tree topologies in newick format used in selection pressure analyses with PAML.