Skip to main content
Log in

Patterns of Substitution Rate Variation at Many Nuclear Loci in Two Species Trios in the Brassicaceae Partitioned with ANOVA

  • Original Article
  • Published:
Journal of Molecular Evolution Aims and scope Submit manuscript

Abstract

There are marked variations among loci and among lineages in rates of nucleotide substitution. The generation time hypothesis (GTH) is a neutral explanation for substitution rate heterogeneity that has genomewide application, predicting that species with shorter generation times accumulate DNA sequence substitutions faster than species with longer generation times do since faster genome replication provides more opportunities for mutations to occur and reach fixation by genetic drift. Relatively few studies have rigorously evaluated the GTH in plants, and there are numerous alternative hypotheses for plant substitution rate variation. One major challenge has been finding pairs of closely related plant species with contrasting generation times and appropriate outgroup taxa that all also have DNA sequence data for numerous loci. To test for causes of rate variation, we obtained sequence data for 256 genes for Arabidopsis thaliana, normally reproducing every year, and the biennial Arabidopsis lyrata with three closely related outgroup taxa (Brassica rapa, Capsella grandiflora, and Neslia paniculata) as well as the biennial Brassica oleracea and the annual B. rapa lineage with the outgroup N. paniculata. A sign test indicated that more loci than expected by chance have faster rates of substitution on the branch leading to the annual than to the perennial for one three-species trio but not another. Tajima’s 1D and 2D tests, and a likelihood ratio test that incorporated saturation correction, rejected rate homogeneity for up to 26 genes (up to 14 genes when correcting for multiple tests), consistently showing faster rates for the annual lineage in the Arabidopsis species trio. ANOVA showed significant rate heterogeneity between the Arabidopsis and Brassica species trios (about 6 % of rate variation) and among loci (about 26–32 % of rate variation). The lineage-by-locus interaction which would be caused by locus- and lineage-specific natural selection explained about 13 % of substitution rate variation in one ANOVA model using substitution rates from genes partitioned into odd and even codons but was not a significant effect without partitioned genes. Annual/perennial lineage and species trio by annual/perennial lineage each explained about 1 % of substitution rate variation.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1

Similar content being viewed by others

References

Download references

Acknowledgments

Ashley S. Timko, Nickolas Julian, Sean Burke, Marisa Wagner, Gregory Way, and Joseph R. Boland in the Braverman lab; and Jack McGuire and Christian Yon in the Hamilton lab assisted in the collection, organization, and/or analysis of the sequences. The authors acknowledge the financial support from the Howard Hughes Medical Institute Undergraduate Science education Grant. Funds were also provided from the Saint Joseph’s University Chapter of the Sigma Xi Scientific Research Society to Ashley S. Timko.

Authors’ Contributions

The JMB and MBH jointly conceived, designed, and coordinated the study. JMB and MBH assembled the sequence data, JMB performed the substitution rate analyses, and MBH performed the phylogenetic analyses. BAJ and MBH designed the ANOVAs and all authors implemented the ANOVAs. All authors contributed to writing and editing the manuscript, and approved the final manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to John M. Braverman.

Ethics declarations

Conflict of Interests

The authors declare that they have no conflicts of interests.

Electronic supplementary material

Below is the link to the electronic supplementary material.

239_2016_9752_MOESM1_ESM.xlsx

Supplementary material 1. Master results table in Microsoft Excel format with information about the sequences sampled and information calculated and summarized in this paper. (XLSX 317 kb)

Supplementary material 2. SAS script used to carry out ANOVAs. (TXT 6 kb)

Supplementary material 3. SAS output in rtf format (best viewed in Microsoft Word). (RTF 1662 kb)

Supplementary material 4. Input csv table for SAS script. (CSV 30 kb)

Supplementary material 5. Input csv table for SAS script. (CSV 63 kb)

Supplementary material 6. R script used to carry out ANOVAs (TXT 6 kb)

Supplementary material 7. Output from R script used to carry out ANOVAs. (TXT 15 kb)

Supplementary material 8 Table of accession numbers of sequences used for the Fig. 1 phylogeny. (DOCX 51 kb)

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Braverman, J.M., Hamilton, M.B. & Johnson, B.A. Patterns of Substitution Rate Variation at Many Nuclear Loci in Two Species Trios in the Brassicaceae Partitioned with ANOVA. J Mol Evol 83, 97–109 (2016). https://doi.org/10.1007/s00239-016-9752-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s00239-016-9752-x

Keywords

Navigation