The effectiveness of selection in a species affects the direction of amino acid frequency evolution.

Nearly neutral theory predicts that species with higher effective population size (N e ) are better able to purge slightly deleterious mutations. We compare evolution in high-N e vs. low-N e vertebrates to reveal which amino acid frequencies are subject to subtle selective preferences. We take three complementary approaches, two measuring flux and one measuring outcomes. First, we fit non-stationary substitution models of amino acid flux using maximum likelihood, comparing the high-N e clade of rodents and lagomorphs to its low-N e sister clade of primates and colugos. Second, we compare evolutionary outcomes across a wider range of vertebrates, via correlations between amino acid frequencies and N e . Third, we dissect the details of flux in human, chimpanzee, mouse, and rat, as scored by parsimony - this also enables comparison to a historical paper. All three methods agree on which amino acids are preferred under more effective selection. Preferred amino acids tend to be smaller, less costly to synthesize, and to promote intrinsic structural disorder. Parsimony-induced bias in the historical study produces an apparent reduction in structural disorder, perhaps driven by slightly deleterious substitutions. Within highly exchangeable pairs of amino acids, arginine is strongly preferred over lysine, and valine over isoleucine, consistent with more effective selection preferring a marginally larger free energy of folding. These two preferences match differences between thermophiles and mesophilic relatives. These results reveal the biophysical consequences of mutation-selection-drift balance, and demonstrate the utility of nearly neutral theory for understanding protein evolution.

The protein domains of vertebrate species in which selection is more effective have greater intrinsic structural disorder.

The nearly neutral theory of molecular evolution posits variation among species in the effectiveness of selection. In an idealized model, the census population size determines both this minimum magnitude of the selection coefficient required for deleterious variants to be reliably purged, and the amount of neutral diversity. Empirically, an "effective population size" is often estimated from the amount of putatively neutral genetic diversity and is assumed to also capture a species' effectiveness of selection. A potentially more direct measure of the effectiveness of selection is the degree to which selection maintains preferred codons. However, past metrics that compare codon bias across species are confounded by among-species variation in %GC content and/or amino acid composition. Here we propose a new Codon Adaptation Index of Species (CAIS), based on Kullback-Leibler divergence, that corrects for both confounders. We demonstrate the use of CAIS correlations, as well as the Effective Number of Codons, to show that the protein domains of more highly adapted vertebrate species evolve higher intrinsic structural disorder.

nQMaker: Estimating Time Nonreversible Amino Acid Substitution Models.

Amino acid substitution models are a key component in phylogenetic analyses of protein sequences. All commonly used amino acid models available to date are time-reversible, an assumption designed for computational convenience but not for biological reality. Another significant downside to time-reversible models is that they do not allow inference of rooted trees without outgroups. In this article, we introduce a maximum likelihood approach nQMaker, an extension of the recently published QMaker method, that allows the estimation of time nonreversible amino acid substitution models and rooted phylogenetic trees from a set of protein sequence alignments. We show that the nonreversible models estimated with nQMaker are a much better fit to empirical alignments than pre-existing reversible models, across a wide range of data sets including mammals, birds, plants, fungi, and other taxa, and that the improvements in model fit scale with the size of the data set. Notably, for the recently published plant and bird trees, these nonreversible models correctly recovered the commonly estimated root placements with very high-statistical support without the need to use an outgroup. We provide nQMaker as an easy-to-use feature in the IQ-TREE software (http://www.iqtree.org), allowing users to estimate nonreversible models and rooted phylogenies from their own protein data sets. The data sets and scripts used in this article are available at https://doi.org/10.5061/dryad.3tx95x6hx. [amino acid sequence analyses; amino acid substitution models; maximum likelihood model estimation; nonreversible models; phylogenetic inference; reversible models.].

Reconstructing the evolutionary history of nitrogenases: Evidence for ancestral molybdenum-cofactor utilization.

The nitrogenase metalloenzyme family, essential for supplying fixed nitrogen to the biosphere, is one of life's key biogeochemical innovations. The three forms of nitrogenase differ in their metal dependence, each binding either a FeMo-, FeV-, or FeFe-cofactor where the reduction of dinitrogen takes place. The history of nitrogenase metal dependence has been of particular interest due to the possible implication that ancient marine metal availabilities have significantly constrained nitrogenase evolution over geologic time. Here, we reconstructed the evolutionary history of nitrogenases, and combined phylogenetic reconstruction, ancestral sequence inference, and structural homology modeling to evaluate the potential metal dependence of ancient nitrogenases. We find that active-site sequence features can reliably distinguish extant Mo-nitrogenases from V- and Fe-nitrogenases and that inferred ancestral sequences at the deepest nodes of the phylogeny suggest these ancient proteins most resemble modern Mo-nitrogenases. Taxa representing early-branching nitrogenase lineages lack one or more biosynthetic nifE and nifN genes that both contribute to the assembly of the FeMo-cofactor in studied organisms, suggesting that early Mo-nitrogenases may have utilized an alternate and/or simplified pathway for cofactor biosynthesis. Our results underscore the profound impacts that protein-level innovations likely had on shaping global biogeochemical cycles throughout the Precambrian, in contrast to organism-level innovations that characterize the Phanerozoic Eon.

Non-invasive placentation in the marsupials Macropus eugenii (Macropodidae) and Trichosurus vulpecula (Phalangeridae) involves redistribution of uterine Desmoglein-2.

In mammalian pregnancy, the uterus is remodeled to become receptive to embryonic implantation. Since non-invasive placentation in marsupials is likely derived from invasive placentation, and is underpinned by intra-uterine conflict between mother and embryo, species with non-invasive placentation may employ a variety of molecular mechanisms to maintain an intact uterine epithelium and to prevent embryonic invasion. Identifying such modifications to the uterine epithelium of marsupial species with non-invasive placentation is key to understanding how conflict is mediated during pregnancy in different mammalian groups. Desmoglein-2, involved in maintaining lateral cell-cell adhesion of the uterine epithelium, is redistributed before implantation to facilitate embryo invasion in mammals with invasive placentation. We identified localization patterns of this cell adhesion molecule throughout pregnancy in two marsupial species with non-invasive placentation, the tammar wallaby (Macropus eugenii; Macropodidae), and the brushtail possum (Trichosurus vulpecula; Phalangeridae). Interestingly, Desmoglein-2 redistribution also occurs in both M. eugenii and T. vulpecula, suggesting that cell adhesion, and thus integrity of the uterine epithelium, is reduced during implantation regardless of placental type, and may be an important component of uterine remodeling. Desmoglein-2 also localizes to the mesenchymal stromal cells of M. eugenii and to epithelial cell nuclei in T. vulpecula, suggesting its involvement in cellular processes that are independent of adhesion and may compensate for reduced lateral adhesion in the uterine epithelium. We conclude that non-invasive placentation in marsupials involves diverse and complementary strategies to maintain an intact epithelial barrier.

Uterine remodelling during pregnancy and pseudopregnancy in the brushtail possum (Trichosurus vulpecula; Phalangeridae).

The formation of a placenta is critical for successful mammalian pregnancy and requires remodelling of the uterine epithelium. In eutherian mammals, remodelling involves specific morphological changes that often correlate with the mode of embryonic attachment. Given the differences between marsupial and eutherian placentae, formation of a marsupial placenta may involve patterns of uterine remodelling that are different from those in eutherians. Here we present a detailed morphological study of the uterus of the brushtail possum (Trichosurus vulpecula; Phalangeridae) throughout pregnancy, using both scanning and transmission electron microscopy, to identify whether uterine changes in marsupials correlate with mode of embryonic attachment as they do in eutherian mammals. The uterine remodelling of T. vulpecula is similar to that of eutherian mammals with the same mode of embryonic attachment (non-invasive, epitheliochorial placentation). The morphological similarities include development of large apical projections, and a decrease in the diffusion distance for haemotrophes around the period of embryonic attachment. Importantly, remodelling of the uterus in T. vulpecula during pregnancy differs from that of a marsupial species with non-invasive attachment (Macropus eugenii; Macropodidae) but is similar to that of a marsupial with invasive attachment (Monodelphis domestica; Didelphidae). We conclude that modes of embryonic attachment may not be typified by a particular suite of uterine changes in marsupials, as is the case for eutherian mammals, and that uterine remodelling may instead reflect phylogenetic relationships between marsupial lineages.