Promiscuity and the rate of molecular evolution at primate immunity genes

Gene Protein Sequence Evolution Can Predict the Rapid Divergence of Ovariole Numbers in the Drosophila melanogaster Subgroup

Ovaries play key roles in fitness and evolution: they are essential female reproductive structures that develop and house the eggs in sexually reproducing animals. In Drosophila, the mature ovary contains multiple tubular egg-producing structures known as ovarioles. Ovarioles arise from somatic cellular structures in the larval ovary called terminal filaments (TFs), formed by TF cells and subsequently enclosed by sheath (SH) cells. As in many other insects, ovariole number per female varies extensively in Drosophila. At present, however, there is a striking gap of information on genetic mechanisms and evolutionary forces that shape the well-documented rapid interspecies divergence of ovariole numbers. To address this gap, here we studied genes associated with Drosophila melanogaster ovariole number or functions based on recent experimental and transcriptional datasets from larval ovaries, including TFs and SH cells, and assessed their rates and patterns of molecular evolution in five closely related species of the melanogaster subgroup that exhibit species-specific differences in ovariole numbers. From comprehensive analyses of protein sequence evolution (dN/dS), branch-site positive selection, expression specificity (tau), and phylogenetic regressions (phylogenetic generalized least squares), we report evidence of 42 genes that showed signs of playing roles in the genetic basis of interspecies evolutionary change of Drosophila ovariole number. These included the signaling genes upd2 and Ilp5 and extracellular matrix genes vkg and Col4a1, whose dN/dS predicted ovariole numbers among species. Together, we propose a model whereby a set of ovariole-involved gene proteins have an enhanced evolvability, including adaptive evolution, facilitating rapid shifts in ovariole number among Drosophila species.

Convergent rates of protein evolution identify novel targets of sexual selection in primates

Sexual selection is the differential reproductive success of individuals, resulting from competition for mates, mate choice, or success in fertilization. In primates, this selective pressure often leads to the development of exaggerated traits which play a role in sexual competition and successful reproduction. In order to gain insight into the mechanisms driving the development of sexually selected traits, we used an unbiased genome-wide approach across 21 primate species to correlate individual rates of protein evolution to relative testes size and sexual dimorphism in body size, 2 anatomical hallmarks of sexual selection in mammals. Among species with presumed high levels of sperm competition, we detected strong conservation of testes-specific proteins responsible for spermatogenesis and ciliary form and function. In contrast, we identified accelerated evolution of female reproductive proteins expressed in the vagina, cervix, and fallopian tubes in these same species. Additionally, we found accelerated protein evolution in lymphoid tissue, indicating that adaptive immune functions may also be influenced by sexual selection. This study demonstrates the distinct complexity of sexual selection in primates revealing contrasting patterns of protein evolution between male and female reproductive tissues.

Genes underlying the evolution of tetrapod testes size

BACKGROUND

Testes vary widely in mass relative to body mass across species, but we know very little about which genes underlie and contribute to such variation. This is partly because evidence for which genes are implicated in testis size variation tends to come from investigations involving just one or a few species. Contemporary comparative phylogenetic methods provide an opportunity to test candidate genes for their role in phenotypic change at a macro-evolutionary scale-across species and over millions of years. Previous attempts to detect genotype-phenotype associations across species have been limited in that they can only detect where genes have driven directional selection (e.g. brain size increase).

RESULTS

Here, we introduce an approach that uses rates of evolutionary change to overcome this limitation to test whether any of twelve candidate genes have driven testis size evolution across tetrapod vertebrates-regardless of directionality. We do this by seeking a relationship between the rates of genetic and phenotypic evolution. Our results reveal five genes (Alkbh5, Dmrtb1, Pld6, Nlrp3, Sp4) that each have played unique and complex roles in tetrapod testis size diversity. In all five genes, we find strong significant associations between the rate of protein-coding substitutions and the rate of testis size evolution. Such an association has never, to our knowledge, been tested before for any gene or phenotype.

CONCLUSIONS

We describe a new approach to tackle one of the most fundamental questions in biology: how do individual genes give rise to biological diversity? The ability to detect genotype-phenotype associations that have acted across species has the potential to build a picture of how natural selection has sculpted phenotypic change over millions of years.

Contrasting patterns of molecular evolution in metazoan germ line genes

BACKGROUND

Germ lines are the cell lineages that give rise to the sperm and eggs in animals. The germ lines first arise from primordial germ cells (PGCs) during embryogenesis: these form from either a presumed derived mode of preformed germ plasm (inheritance) or from an ancestral mechanism of inductive cell-cell signalling (induction). Numerous genes involved in germ line specification and development have been identified and functionally studied. However, little is known about the molecular evolutionary dynamics of germ line genes in metazoan model systems.

RESULTS

Here, we studied the molecular evolution of germ line genes within three metazoan model systems. These include the genus Drosophila (N=34 genes, inheritance), the fellow insect Apis (N=30, induction), and their more distant relative Caenorhabditis (N=23, inheritance). Using multiple species and established phylogenies in each genus, we report that germ line genes exhibited marked variation in the constraint on protein sequence divergence (dN/dS) and codon usage bias (CUB) within each genus. Importantly, we found that de novo lineage-specific inheritance (LSI) genes in Drosophila (osk, pgc) and in Caenorhabditis (pie-1, pgl-1), which are essential to germ plasm functions under the derived inheritance mode, displayed rapid protein sequence divergence relative to the other germ line genes within each respective genus. We show this may reflect the evolution of specialized germ plasm functions and/or low pleiotropy of LSI genes, features not shared with other germ line genes. In addition, we observed that the relative ranking of dN/dS and of CUB between genera were each more strongly correlated between Drosophila and Caenorhabditis, from different phyla, than between Drosophila and its insect relative Apis, suggesting taxonomic differences in how germ line genes have evolved.

CONCLUSIONS

Taken together, the present results advance our understanding of the evolution of animal germ line genes within three well-known metazoan models. Further, the findings provide insights to the molecular evolution of germ line genes with respect to LSI status, pleiotropy, adaptive evolution as well as PGC-specification mode.

Infectious disease and group size: more than just a numbers game

Increased risk of infectious disease is assumed to be a major cost of group living, yet empirical evidence for this effect is mixed. We studied whether larger social groups are more subdivided structurally. If so, the social subdivisions that form in larger groups may act as barriers to the spread of infection, weakening the association between group size and infectious disease. To investigate this 'social bottleneck' hypothesis, we examined the association between group size and four network structure metrics in 43 vertebrate and invertebrate species. We focused on metrics involving modularity, clustering, distance and centralization. In a meta-analysis of intraspecific variation in social networks, modularity showed positive associations with network size, with a weaker but still positive effect in cross-species analyses. Network distance also showed a positive association with group size when using intraspecific variation. We then used a theoretical model to explore the effects of subgrouping relative to other effects that influence disease spread in socially structured populations. Outbreaks reached higher prevalence when groups were larger, but subgrouping reduced prevalence. Subgrouping also acted as a 'brake' on disease spread between groups. We suggest research directions to understand the conditions under which larger groups become more subdivided, and to devise new metrics that account for subgrouping when investigating the links between sociality and infectious disease risk.

Sexual conflict and seminal fluid proteins: a dynamic landscape of sexual interactions

Sexual reproduction requires coordinated contributions from both sexes to proceed efficiently. However, the reproductive strategies that the sexes adopt often have the potential to give rise to sexual conflict because they can result in divergent, sex-specific costs and benefits. These conflicts can occur at many levels, from molecular to behavioral. Here, we consider sexual conflict mediated through the actions of seminal fluid proteins. These proteins provide many excellent examples in which to trace the operation of sexual conflict from molecules through to behavior. Seminal fluid proteins are made by males and provided to females during mating. As agents that can modulate egg production at several steps, as well as reproductive behavior, sperm "management," and female feeding, activity, and longevity, the actions of seminal proteins are prime targets for sexual conflict. We review these actions in the context of sexual conflict. We discuss genomic signatures in seminal protein (and related) genes that are consistent with current or previous sexual conflict. Finally, we note promising areas for future study and highlight real-world practical situations that will benefit from understanding the nature of sexual conflicts mediated by seminal proteins.

Primate vaginal microbiomes exhibit species specificity without universal Lactobacillus dominance

Bacterial communities colonizing the reproductive tracts of primates (including humans) impact the health, survival and fitness of the host, and thereby the evolution of the host species. Despite their importance, we currently have a poor understanding of primate microbiomes. The composition and structure of microbial communities vary considerably depending on the host and environmental factors. We conducted comparative analyses of the primate vaginal microbiome using pyrosequencing of the 16S rRNA genes of a phylogenetically broad range of primates to test for factors affecting the diversity of primate vaginal ecosystems. The nine primate species included: humans (Homo sapiens), yellow baboons (Papio cynocephalus), olive baboons (Papio anubis), lemurs (Propithecus diadema), howler monkeys (Alouatta pigra), red colobus (Piliocolobus rufomitratus), vervets (Chlorocebus aethiops), mangabeys (Cercocebus atys) and chimpanzees (Pan troglodytes). Our results indicated that all primates exhibited host-specific vaginal microbiota and that humans were distinct from other primates in both microbiome composition and diversity. In contrast to the gut microbiome, the vaginal microbiome showed limited congruence with host phylogeny, and neither captivity nor diet elicited substantial effects on the vaginal microbiomes of primates. Permutational multivariate analysis of variance and Wilcoxon tests revealed correlations among vaginal microbiota and host species-specific socioecological factors, particularly related to sexuality, including: female promiscuity, baculum length, gestation time, mating group size and neonatal birth weight. The proportion of unclassified taxa observed in nonhuman primate samples increased with phylogenetic distance from humans, indicative of the existence of previously unrecognized microbial taxa. These findings contribute to our understanding of host-microbe variation and coevolution, microbial biogeography, and disease risk, and have important implications for the use of animal models in studies of human sexual and reproductive diseases.

Sequence diversity of Pan troglodytes subspecies and the impact of WFDC6 selective constraints in reproductive immunity

Recent efforts have attempted to describe the population structure of common chimpanzee, focusing on four subspecies: Pan troglodytes verus, P. t. ellioti, P. t. troglodytes, and P. t. schweinfurthii. However, few studies have pursued the effects of natural selection in shaping their response to pathogens and reproduction. Whey acidic protein (WAP) four-disulfide core domain (WFDC) genes and neighboring semenogelin (SEMG) genes encode proteins with combined roles in immunity and fertility. They display a strikingly high rate of amino acid replacement (d_N/d_S), indicative of adaptive pressures during primate evolution. In human populations, three signals of selection at the WFDC locus were described, possibly influencing the proteolytic profile and antimicrobial activities of the male reproductive tract. To evaluate the patterns of genomic variation and selection at the WFDC locus in chimpanzees, we sequenced 17 WFDC genes and 47 autosomal pseudogenes in 68 chimpanzees (15 P. t. troglodytes, 22 P. t. verus, and 31 P. t. ellioti). We found a clear differentiation of P. t. verus and estimated the divergence of P. t. troglodytes and P. t. ellioti subspecies in 0.173 Myr; further, at the WFDC locus we identified a signature of strong selective constraints common to the three subspecies in WFDC6—a recent paralog of the epididymal protease inhibitor EPPIN. Overall, chimpanzees and humans do not display similar footprints of selection across the WFDC locus, possibly due to different selective pressures between the two species related to immune response and reproductive biology.

Immune anticipation of mating in Drosophila: Turandot M promotes immunity against sexually transmitted fungal infections

Although it is well known that mating increases the risk of infection, we do not know how females mitigate the fitness costs of sexually transmitted infections (STIs). It has recently been shown that female fruitflies, Drosophila melanogaster, specifically upregulate two members of the Turandot family of immune and stress response genes, Turandot M and Turandot C (TotM and TotC), when they hear male courtship song. Here, we use the Gal4/UAS RNAi gene knockdown system to test whether the expression of these genes provides fitness benefits for females infected with the entomopathogenic fungus, Metarhizium robertsii under sexual transmission. As a control, we also examined the immunity conferred by Dorsal-related immunity factor (Dif), a central component of the Toll signalling pathway thought to provide immunity against fungal infections. We show that TotM, but not TotC or Dif, provides survival benefits to females following STIs, but not after direct topical infections. We also show that though the expression of TotM provides fecundity benefits for healthy females, it comes at a cost to their survival, which helps to explain why TotM is not constitutively expressed. Together, these results show that the anticipatory expression of TotM promotes specific immunity against fungal STIs and suggest that immune anticipation is more common than currently appreciated.

Evolutionary Modeling of Genotype-Phenotype Associations, and Application to Primate Coding and Non-coding mtDNA Rate Variation

Variation in substitution rates across a phylogeny can be indicative of shifts in the evolutionary dynamics of a protein or non-protein coding regions. One way to understand these signals is to seek the phenotypic correlates of rate variation. Here, we extended a previously published likelihood method designed to detect evolutionary associations between genotypic evolutionary rate and phenotype over a phylogeny. In simulation with two discrete categories of phenotype, the method has a low false-positive rate and detects greater than 80% of true-positives with a tree length of three or greater and a three-fold or greater change in substitution rate given the phenotype. In addition, we successfully extend the test from two to four phenotype categories and evaluated its performance. We then applied the method to two major hypotheses for rate variation in the mitochondrial genome of primates-longevity and generation time as well as body mass which is correlated with many aspects of life history-using three categories of phenotype through discretization of continuous values. Similar to previous results for mammals, we find that the majority of mitochondrial protein-coding genes show associations consistent with the longevity and body mass predictions and that the predominant signal of association comes from the third codon position. We also found a significant association between maximum lifespan and the evolutionary rate of the control region of the mtDNA. In contrast, 24 protein-coding genes from the nuclear genome do not show a consistent pattern of association, which is inconsistent with the generation time hypothesis. These results show the extended method can robustly identify genotype-phenotype associations up to at least four phenotypic categories, and demonstrate the successful application of the method to study factors affecting neutral evolutionary rate in protein-coding and non-coding loci.

Polyandry Has No Detectable Mortality Cost in Female Mammals

In several taxonomic groups, females mate with several males during a single reproductive cycle. Although there is evidence that polyandry provides some benefits to females, it often involves mortality costs. However, empirical evidences of mortality costs of polyandry have so far been reported only in invertebrates. Whether polyandry has mortality costs in vertebrates is currently unknown. In the present study, we aimed to fill the gap by investigating the relationships between the level of polyandry (measured either by male relative testes mass or the percentage of multiple paternities) and female patterns of mortality across mammals. While we found that the two metrics of female mortality co-varied with pace of life, we did not find any evidence that polyandry leads to either decreased median lifespan or increased aging rate in mammals. We discuss such an absence of detectable mortality costs of polyandry in female mammals in light of recent advances in the study of mammalian reproductive biology and life-history tactics.

Birth-and-death of KLK3 and KLK2 in primates: evolution driven by reproductive biology

The kallikrein (KLK) gene family comprises the largest uninterrupted locus of serine proteases in the human genome and represents a notable case for studying the evolutionary fate of duplicated genes. In primates, a recent duplication event gave rise to KLK2 and KLK3, both encoding essential proteins for the cascade of seminal plasma liquefaction. We reconstructed the evolutionary history of KLK2 and KLK3 by comparative analysis of the orthologous sequences from 22 primate species, calculated d(N)/d(S) ratios, and addressed the hypothesis of coevolution with their substrates, the semenogelins (SEMG1 and SEMG2). Our findings support the placement of the KLK2-KLK3 duplication in the Catarrhini ancestor and unveil the frequent loss of KLK2 throughout primate evolution by different genomic mechanisms, including unequal crossing-over, deletions, and pseudogenization. We provide evidences for an adaptive evolution of KLK3 toward an expanded enzymatic spectrum, with an effect on the hydrolysis of semen coagulum. Furthermore, we found associations between mating system, the number of SEMG repeat units, and the number of functional KLK2 and KLK3, suggesting complex evolutionary dynamics shaped by reproductive biology.

Is promiscuity associated with enhanced selection on MHC-DQα in mice (genus Peromyscus)?

Reproductive behavior may play an important role in shaping selection on Major Histocompatibility Complex (MHC) genes. For example, the number of sexual partners that an individual has may affect exposure to sexually transmitted pathogens, with more partners leading to greater exposure and, hence, potentially greater selection for variation at MHC loci. To explore this hypothesis, we examined the strength of selection on exon 2 of the MHC-DQα locus in two species of Peromyscus. While the California mouse (P. californicus) is characterized by lifetime social and genetic monogamy, the deer mouse (P. maniculatus) is socially and genetically promiscuous; consistent with these differences in mating behavior, the diversity of bacteria present within the reproductive tracts of females is significantly greater for P. maniculatus. To test the prediction that more reproductive partners and exposure to a greater range of sexually transmitted pathogens are associated with enhanced diversifying selection on genes responsible for immune function, we compared patterns and levels of diversity at the Class II MHC-DQα locus in sympatric populations of P. maniculatus and P. californicus. Using likelihood based analyses, we show that selection is enhanced in the promiscuous P. maniculatus. This study is the first to compare the strength of selection in wild sympatric rodents with known differences in pathogen milieu.

The molecular evolution of animal reproductive tract proteins: what have we learned from mating-system comparisons?

Postcopulatory sexual selection is thought to drive the rapid evolution of reproductive tract genes in many animals. Recently, a number of studies have sought to test this hypothesis by examining the effects of mating system variation on the evolutionary rates of reproductive tract genes. Perhaps surprisingly, there is relatively little evidence that reproductive proteins evolve more rapidly in species subject to strong postcopulatory sexual selection. This emerging trend may suggest that other processes, such as host-pathogen interactions, are the main engines of rapid reproductive gene evolution. I suggest that such a conclusion is as yet unwarranted; instead, I propose that more rigorous analytical techniques, as well as multigene and population-based approaches, are required for a full understanding of the consequences of mating system variation for the evolution of reproductive tract genes.