A maximum likelihood method for reconstruction of the evolution of eukaryotic gene structure.
Spliceosomal introns are one of the principal distinctive features of eukaryotes. Nevertheless, different large-scale studies disagree about even the most basic features of their evolution. In order to come up with a more reliable reconstruction of intron evolution, we developed a model that is far more comprehensive than previous ones. This model is rich in parameters, and estimating them accurately is infeasible by straightforward likelihood maximization. Thus, we have developed an expectation-maximization algorithm that allows for efficient maximization. Here, we outline the model and describe the expectation-maximization algorithm in detail. Since the method works with intron presence-absence maps, it is expected to be instrumental for the analysis of the evolution of other binary characters as well.
Carmel
Liran
L
National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD, USA.
Rogozin
Igor B
IB
Wolf
Yuri I
YI
Koonin
Eugene V
EV
D000818:Animals; D001483:Base Sequence; D005057:Eukaryotic Cells / Q000378:metabolism* / Q000502:physiology; D019143:Evolution, Molecular; D005796:Genes; D006801:Humans; D016013:Likelihood Functions; D008962:Models, Theoretical