A survey of ficolin-3 activity in Systemic Lupus Erythematosus reveals a link to hematological disease manifestations and autoantibody profile.

The complement system plays a central role in the pathogenesis of Systemic Lupus Erythematosus (SLE), but most studies have focused on the classical pathway. Ficolin-3 is the main initiator of the lectin pathway of complement in humans, but its role in systemic autoimmune disease has not been conclusively determined. Here, we combined biochemical and genetic approaches to assess the contribution of ficolin-3 to SLE risk and disease manifestations. Ficolin-3 activity was measured by a functional assay in serum or plasma samples from Swedish SLE patients (n = 786) and controls matched for age and sex (n = 566). Genetic variants in an extended 300 kb genomic region spanning the FCN3 locus were analyzed for their association with ficolin-3 activity and SLE manifestations in a Swedish multicenter cohort (n = 985). Patients with ficolin-3 activity in the highest tertile showed a strong enrichment in an SLE cluster defined by anti-Sm/DNA/nucleosome antibodies (OR 3.0, p < 0.001) and had increased rates of hematological disease (OR 1.4, p = 0.078) and lymphopenia (OR = 1.6, p = 0.039). Genetic variants associated with low ficolin-3 activity mapped to an extended haplotype in high linkage disequilibrium upstream of the FCN3 gene. Patients carrying the lead genetic variant associated with low ficolin-3 activity had a lower frequency of hematological disease (OR 0.67, p = 0.018) and lymphopenia (OR 0.63, p = 0.031) and fewer autoantibodies (p = 0.0019). Loss-of-function variants in the FCN3 gene were not associated with SLE, but four (0.5 %) SLE patients developed acquired ficolin-3 deficiency where ficolin-3 activity in serum was depleted following diagnosis of SLE. Taken together, our results provide genetic and biochemical evidence that implicate the lectin pathway in hematological SLE manifestations. We also identify lectin pathway activation through ficolin-3 as a factor that contributes to the autoantibody response in SLE.

Prevalence of Wolbachia infection in field natural population of the apricot seed wasp Eurytoma samsonowi (Hymenoptera: Eurytomidae).

Obligate endosymbiont bacteria associated with insects are naturally providing their hosts with essential nutrients such as vitamins and amino acids and biological services including protection from pathogens. In this study, we aimed to investigate the presence of Wolbachia infection among males and females of the parasitic apricot seed wasp (ASW) Eurytoma samsonowi Vassiliev (Vassiliev Petrograd 11: 1-15, 1915) (Hymenoptera: Eurytomidae), a very harmful pest of apricot (Prunus armeniaca), in the oasis of Gafsa, Southern-West of Tunisia. The detection of Wolbachia infection was assessed based on the amplification of the Wolbachia surface protein (wsp) gene and a multilocus sequence typing (MLST) as a universal genotyping tool for Wolbachia involving the analyses of genes gatB, coxA, hcpA, fbpA, and ftsz. Confirming the screening results, Wolbachia was detected in the natural apricot wasp for the first time, with a significant difference between males (5%) and females (59%) based on wsp gene. All Wolbachia strains identified in E. samsonowi were clustered among supergroups B of Wolbachia.

Proteome-scale mapping of binding sites in the unstructured regions of the human proteome.

Specific protein-protein interactions are central to all processes that underlie cell physiology. Numerous studies have together identified hundreds of thousands of human protein-protein interactions. However, many interactions remain to be discovered, and low affinity, conditional, and cell type-specific interactions are likely to be disproportionately underrepresented. Here, we describe an optimized proteomic peptide-phage display library that tiles all disordered regions of the human proteome and allows the screening of ~ 1,000,000 overlapping peptides in a single binding assay. We define guidelines for processing, filtering, and ranking the results and provide PepTools, a toolkit to annotate the identified hits. We uncovered >2,000 interaction pairs for 35 known short linear motif (SLiM)-binding domains and confirmed the quality of the produced data by complementary biophysical or cell-based assays. Finally, we show how the amino acid resolution-binding site information can be used to pinpoint functionally important disease mutations and phosphorylation events in intrinsically disordered regions of the proteome. The optimized human disorderome library paired with PepTools represents a powerful pipeline for unbiased proteome-wide discovery of SLiM-based interactions.

Sexual conflict drives micro- and macroevolution of sexual dimorphism in immunity.

BACKGROUND: Sexual dimorphism in immunity is believed to reflect sex differences in reproductive strategies and trade-offs between competing life history demands. Sexual selection can have major effects on mating rates and sex-specific costs of mating and may thereby influence sex differences in immunity as well as associated host-pathogen dynamics. Yet, experimental evidence linking the mating system to evolved sexual dimorphism in immunity are scarce and the direct effects of mating rate on immunity are not well established. Here, we use transcriptomic analyses, experimental evolution and phylogenetic comparative methods to study the association between the mating system and sexual dimorphism in immunity in seed beetles, where mating causes internal injuries in females. RESULTS: We demonstrate that female phenoloxidase (PO) activity, involved in wound healing and defence against parasitic infections, is elevated relative to males. This difference is accompanied by concomitant sex differences in the expression of genes in the prophenoloxidase activating cascade. We document substantial phenotypic plasticity in female PO activity in response to mating and show that experimental evolution under enforced monogamy (resulting in low remating rates and reduced sexual conflict relative to natural polygamy) rapidly decreases female (but not male) PO activity. Moreover, monogamous females had evolved increased tolerance to bacterial infection unrelated to mating, implying that female responses to costly mating may trade off with other aspects of immune defence, an hypothesis which broadly accords with the documented sex differences in gene expression. Finally, female (but not male) PO activity shows correlated evolution with the perceived harmfulness of male genitalia across 12 species of seed beetles, suggesting that sexual conflict has a significant influence on sexual dimorphisms in immunity in this group of insects. CONCLUSIONS: Our study provides insights into the links between sexual conflict and sexual dimorphism in immunity and suggests that selection pressures moulded by mating interactions can lead to a sex-specific mosaic of immune responses with important implications for host-pathogen dynamics in sexually reproducing organisms.

Negative frequency dependent selection contributes to the maintenance of a global polymorphism in mitochondrial DNA.

BACKGROUND: Understanding the forces that maintain diversity across a range of scales is at the very heart of biology. Frequency-dependent processes are generally recognized as the most central process for the maintenance of ecological diversity. The same is, however, not generally true for genetic diversity. Negative frequency dependent selection, where rare genotypes have an advantage, is often regarded as a relatively weak force in maintaining genetic variation in life history traits because recombination disassociates alleles across many genes. Yet, many regions of the genome show low rates of recombination and genetic variation in such regions (i.e., supergenes) may in theory be upheld by frequency dependent selection. RESULTS: We studied what is essentially a ubiquitous life history supergene (i.e., mitochondrial DNA) in the fruit fly Drosophila subobscura, showing sympatric polymorphism with two main mtDNA genotypes co-occurring in populations world-wide. Using an experimental evolution approach involving manipulations of genotype starting frequencies, we show that negative frequency dependent selection indeed acts to maintain genetic variation in this region. Moreover, the strength of selection was affected by food resource conditions. CONCLUSIONS: Our work provides novel experimental support for the view that balancing selection through negative frequency dependency acts to maintain genetic variation in life history genes. We suggest that the emergence of negative frequency dependent selection on mtDNA is symptomatic of the fundamental link between ecological processes related to resource use and the maintenance of genetic variation.

Negative frequency dependent selection on sympatric mtDNA haplotypes in Drosophila subobscura.

BACKGROUND: Recent experimental evidence for selection on mitochondrial DNA (mtDNA) has prompted the question as to what processes act to maintain within-population variation in mtDNA. Balancing selection though negative frequency dependent selection (NFDS) among sympatric haplotypes is a possibility, but direct empirical evidence for this is very scarce. FINDINGS: We extend the previous findings of a multi-generation replicated cage experiment in Drosophila subobscura, where mtDNA polymorphism was maintained in a laboratory setting. First, we use a set of Monte Carlo simulations to show that the haplotype frequency dynamics observed are inconsistent with genetic drift alone and most closely match those expected under NFDS. Second, we show that haplotype frequency changes over time were significantly different from those expected under either genetic drift or positive selection but were consistent with those expected under NFSD. CONCLUSIONS: Collectively, our analyses provide novel support for NFDS on mtDNA haplotypes, suggesting that mtDNA polymorphism may at least in part be maintained by balancing selection also in natural populations. We very briefly discuss the possible mechanisms that might be involved.

Genome size correlates with reproductive fitness in seed beetles.

The ultimate cause of genome size (GS) evolution in eukaryotes remains a major and unresolved puzzle in evolutionary biology. Large-scale comparative studies have failed to find consistent correlations between GS and organismal properties, resulting in the 'C-value paradox'. Current hypotheses for the evolution of GS are based either on the balance between mutational events and drift or on natural selection acting upon standing genetic variation in GS. It is, however, currently very difficult to evaluate the role of selection because within-species studies that relate variation in life-history traits to variation in GS are very rare. Here, we report phylogenetic comparative analyses of GS evolution in seed beetles at two distinct taxonomic scales, which combines replicated estimation of GS with experimental assays of life-history traits and reproductive fitness. GS showed rapid and bidirectional evolution across species, but did not show correlated evolution with any of several indices of the relative importance of genetic drift. Within a single species, GS varied by 4-5% across populations and showed positive correlated evolution with independent estimates of male and female reproductive fitness. Collectively, the phylogenetic pattern of GS diversification across and within species in conjunction with the pattern of correlated evolution between GS and fitness provide novel support for the tenet that natural selection plays a key role in shaping GS evolution.

Comparative Genomic Analysis of the Pattern of Evolution of Male and Female Reproductive Proteins in Seed Beetles.

Male seminal fluid proteins often show signs of positive selection and divergent evolution, believed to reflect male-female coevolution. Yet, our understanding of the predicted concerted evolution of seminal fluid proteins and female reproductive proteins is limited. We sequenced, assembled, and annotated the genome of two species of seed beetles allowing a comparative analysis of four closely related species of these herbivorous insects. We compare the general pattern of evolution in genes encoding seminal fluid proteins and female reproductive proteins with those in digestive protein genes and well-conserved reference genes. We found that female reproductive proteins showed an overall ratio of nonsynonymous to synonymous substitutions (ω) similar to that of conserved genes, while seminal fluid proteins and digestive proteins exhibited higher overall ω values. Further, seminal fluid proteins and digestive proteins showed a higher proportion of sites putatively under positive selection, and explicit tests showed no difference in relaxed selection between protein types. Evolutionary rate covariation analyses showed that evolutionary rates among seminal fluid proteins were on average more closely correlated with those in female reproductive proteins than with either digestive or conserved genes. Gene expression showed the expected negative covariation with ω values, except for male-biased genes where this negative relationship was reversed. In conclusion, seminal fluid proteins showed relatively rapid evolution and signs of positive selection. In contrast, female reproductive proteins evolved at a lower rate under selective constraints, on par with genes known to be well conserved. Although our findings provide support for concerted evolution of seminal fluid proteins and female reproductive proteins, they also suggest that these two classes of proteins evolve under partly distinct selective regimes.

A chromosome-level assembly of the seed beetle Callosobruchus maculatus genome with annotation of its repetitive elements.

Callosobruchus maculatus is a major agricultural pest of legume crops worldwide and an established model system in ecology and evolution. Yet, current molecular biological resources for this species are limited. Here, we employ Hi-C sequencing to generate a greatly improved genome assembly and we annotate its repetitive elements in a dedicated in-depth effort where we manually curate and classify the most abundant unclassified repeat subfamilies. We present a scaffolded chromosome-level assembly, which is 1.01 Gb in total length with 86% being contained within the 9 autosomes and the X chromosome. Repetitive sequences accounted for 70% of the total assembly. DNA transposons covered 18% of the genome, with the most abundant superfamily being Tc1-Mariner (9.75% of the genome). This new chromosome-level genome assembly of C. maculatus will enable future genetic and evolutionary studies not only of this important species but of beetles more generally.

Zip14 expression induced by lipopolysaccharides in macrophages attenuates inflammatory response.

OBJECTIVE AND DESIGN: We investigated the role and regulation of zinc transporters in the activation of the inflammatory response in macrophages. Our exploratory computational study found that Zip14 (SLC39A14) was consistently up-regulated in activated macrophages; we therefore focused subsequently on that gene in the mechanistic study. MATERIAL: The expression and function of Zip14 was assessed in primary macrophages obtained by in-vitro differentiation of monocytes from human blood. METHODS: Primary macrophages were subjected to treatments with lipopolysaccharides, cytokines, chemicals, and pharmacological agents. SLC39A14 and inflammatory cytokine gene expressions were assessed by RT-qPCR. Zip14 siRNA knockdown was performed to explore the gene function. RESULTS: Lipopolysaccharide's inflammatory stimulus was a strong inducer of SLC39A14 mRNA expression in macrophages. This induction was dependent on calcium signaling, GC-rich DNA-binding, and NF-κB down-regulation. Impregnation of lipopolysaccharide-stimulated macrophages with the glucocorticoid dexamethasone further enhanced Zip14 expression while reducing interleukin-6 and tumor necrosis factor-α production. Zip14 knockdown in macrophages attenuated the expression and secretion of cytokines, indicating a buffering function for this zinc transporter. CONCLUSIONS: Collectively, our results identified the zinc transporter Zip14 as expressed downstream of lipopolysaccharide signals in macrophages. Zip14 induction had a regulatory function in cytokine production.

Experimental Life History Evolution Results in Sex-specific Evolution of Gene Expression in Seed Beetles.

The patterns of reproductive timing and senescence vary within and across species owing to differences in reproductive strategies, but our understanding of the molecular underpinnings of such variation is incomplete. This is perhaps particularly true for sex differences. We investigated the evolution of sex-specific gene expression associated with life history divergence in replicated populations of the seed beetle Acanthoscelides obtectus, experimentally evolving under (E)arly or (L)ate life reproduction for >200 generations which has resulted in strongly divergent life histories. We detected 1,646 genes that were differentially expressed in E and L lines, consistent with a highly polygenic basis of life history evolution. Only 30% of differentially expressed genes were similarly affected in males and females. The evolution of long life was associated with significantly reduced sex differences in expression, especially in non-reproductive tissues. The expression differences were overall more pronounced in females, in accordance with their greater phenotypic divergence in lifespan. Functional enrichment analysis revealed differences between E and L beetles in gene categories previously implicated in aging, such as mitochondrial function and defense response. The results show that divergent life history evolution can be associated with profound changes in gene expression that alter the transcriptome in a sex-specific way, highlighting the importance of understanding the mechanisms of aging in each sex.

B cell polygenic risk scores associate with anti-dsDNA antibodies and nephritis in systemic lupus erythematosus.

OBJECTIVE: B cell function and autoantibodies are important in SLE pathogenesis. In this work, we aimed to investigate the impact of cumulative SLE B cell genetics on SLE subphenotype and autoantibody profile. METHODS: Female patients with SLE (n=1248) and healthy controls (n=400) were genotyped using Illumina's Global Screening Array. Two polygenic risk scores (PRSs), one representing B cell genes and the other B cell activation genes, were calculated for each individual using risk loci for SLE in genes assigned to B cell-related pathways according to the Kyoto Encyclopedia of Genes and Genomes, Gene Ontology and Reactome Databases. RESULTS: Double-stranded DNA (dsDNA) antibodies were more prevalent among patients with a high compared with a low SLE B cell PRS (OR 1.47 (1.07 to 2.01), p=0.018), and effect sizes were augmented in patients with human leucocyte antigen (HLA) risk haplotypes HLA-DRB1*03:01 and HLA-DRB1*15:01 (DRB1*03/15 -/- (OR 0.99 (0.56 to 1.77), p=0.98; DRB1*03/15 +/- or -/+ (OR 1.64 (1.06 to 2.54), p=0.028; and DRB1*03/15 +/+ (OR 4.47 (1.21 to 16.47), p=0.024). Further, a high compared with a low B cell PRS was associated with low complement levels in DRB1*03/15 +/+ patients (OR 3.92 (1.22 to 12.64), p=0.022). The prevalence of lupus nephritis (LN) was higher in patients with a B cell activation PRS above the third quartile compared with patients below (OR 1.32 (1.00 to 1.74), p=0.048). CONCLUSIONS: High genetic burden related to B cell function is associated with dsDNA antibody development and LN. Assessing B cell PRSs may be important in order to determine immunological pathways influencing SLE and to predict clinical phenotype.

Large-scale phage-based screening reveals extensive pan-viral mimicry of host short linear motifs.

Viruses mimic host short linear motifs (SLiMs) to hijack and deregulate cellular functions. Studies of motif-mediated interactions therefore provide insight into virus-host dependencies, and reveal targets for therapeutic intervention. Here, we describe the pan-viral discovery of 1712 SLiM-based virus-host interactions using a phage peptidome tiling the intrinsically disordered protein regions of 229 RNA viruses. We find mimicry of host SLiMs to be a ubiquitous viral strategy, reveal novel host proteins hijacked by viruses, and identify cellular pathways frequently deregulated by viral motif mimicry. Using structural and biophysical analyses, we show that viral mimicry-based interactions have similar binding strength and bound conformations as endogenous interactions. Finally, we establish polyadenylate-binding protein 1 as a potential target for broad-spectrum antiviral agent development. Our platform enables rapid discovery of mechanisms of viral interference and the identification of potential therapeutic targets which can aid in combating future epidemics and pandemics.

Interaction network of the 14-3-3 protein in the ancient protozoan parasite Giardia duodenalis.

14-3-3s are phosphoserine/phosphotreonine binding proteins that play pivotal roles as regulators of multiple cellular processes in eukaryotes. The flagellated protozoan parasite Giardia duodenalis, the causing agent of giardiasis, is a valuable simplified eukaryotic model. A single 14-3-3 isoform (g14-3-3) is expressed in Giardia, and it is directly involved in the differentiation of the parasite into cyst. To define the overall functions of g14-3-3, the protein interactome has been investigated. A transgenic G. duodenalis strain was engineered to express a FLAG-tagged g14-3-3 under its own promoter. Affinity chromatography coupled with tandem mass spectrometry analysis have been used to purify and identify FLAG-g14-3-3-associated proteins from trophozoites and encysting parasites. A total of 314 putative g14-3-3 interaction partners were identified, including proteins involved in several pathways. Some interactions seemed to be peculiar of one specific stage, while others were shared among the different stages. Furthermore, the interaction of g14-3-3 with the giardial homologue of the CDC7 protein kinase (gCDC7) was characterized, leading to the identification of a multiprotein complex containing not only g14-3-3 and gCDC7 but also a newly identified and highly divergent homologue of DBF4, the putative regulatory subunit of gCDC7. The relevance of g14-3-3 interactions in G. duodenalis biology was discussed.

ELM: the status of the 2010 eukaryotic linear motif resource.

Linear motifs are short segments of multidomain proteins that provide regulatory functions independently of protein tertiary structure. Much of intracellular signalling passes through protein modifications at linear motifs. Many thousands of linear motif instances, most notably phosphorylation sites, have now been reported. Although clearly very abundant, linear motifs are difficult to predict de novo in protein sequences due to the difficulty of obtaining robust statistical assessments. The ELM resource at http://elm.eu.org/ provides an expanding knowledge base, currently covering 146 known motifs, with annotation that includes >1300 experimentally reported instances. ELM is also an exploratory tool for suggesting new candidates of known linear motifs in proteins of interest. Information about protein domains, protein structure and native disorder, cellular and taxonomic contexts is used to reduce or deprecate false positive matches. Results are graphically displayed in a 'Bar Code' format, which also displays known instances from homologous proteins through a novel 'Instance Mapper' protocol based on PHI-BLAST. ELM server output provides links to the ELM annotation as well as to a number of remote resources. Using the links, researchers can explore the motifs, proteins, complex structures and associated literature to evaluate whether candidate motifs might be worth experimental investigation.

A possible genomic footprint of polygenic adaptation on population divergence in seed beetles?

Efforts to unravel the genomic basis of incipient speciation are hampered by a mismatch between our toolkit and our understanding of the ecology and genetics of adaptation. While the former is focused on detecting selective sweeps involving few independently acting or linked speciation genes, the latter states that divergence typically occurs in polygenic traits under stabilizing selection. Here, we ask whether a role of stabilizing selection on polygenic traits in population divergence may be unveiled by using a phenotypically informed integrative approach, based on genome-wide variation segregating in divergent populations. We compare three divergent populations of seed beetles (Callosobruchus maculatus) where previous work has demonstrated a prominent role for stabilizing selection on, and population divergence in, key life history traits that reflect rate-dependent metabolic processes. We derive and assess predictions regarding the expected pattern of covariation between genetic variation segregating within populations and genetic differentiation between populations. Population differentiation was considerable (mean F ST = 0.23-0.26) and was primarily built by genes showing high selective constraints and an imbalance in inferred selection in different populations (positive Tajima's D NS in one and negative in one), and this set of genes was enriched with genes with a metabolic function. Repeatability of relative population differentiation was low at the level of individual genes but higher at the level of broad functional classes, again spotlighting metabolic genes. Absolute differentiation (d XY) showed a very different general pattern at this scale of divergence, more consistent with an important role for genetic drift. Although our exploration is consistent with stabilizing selection on polygenic metabolic phenotypes as an important engine of genome-wide relative population divergence and incipient speciation in our study system, we note that it is exceedingly difficult to firmly exclude other scenarios.

In silico promoter analysis and functional validation identify CmZFH, the co-regulator of hypoxia-responsive genes CmScylla and CmLPCAT.

Oxygen (O2) plays an essential role in aerobic organisms including terrestrial insects. Under hypoxic stress, the cowpea bruchid (Callosobruchus maculatus) ceases feeding and growth. However, larvae, particularly 4th instar larvae exhibit very high tolerance to hypoxia and can recover normal growth once brought to normoxia. To better understand the molecular mechanism that enables insects to cope with low O2 stress, we performed RNA-seq to distinguish hypoxia-responsive genes in midguts and subsequently identified potential common cis-elements in promoters of hypoxia-induced and -repressed genes, respectively. Selected elements were subjected to gel-shift and transient transfection assays to confirm their cis-regulatory function. Of these putative common cis-elements, AREB6 appeared to regulate the expression of CmLPCAT and CmScylla, two hypoxia-induced genes. CmZFH, the putative AREB6-binding protein, was hypoxia-inducible. Transient expression of CmZFH in Drosophila S2 cells activated CmLPCAT and CmScylla, and their induction was likely through interaction of CmZFH with AREB6. Binding to AREB6 was further confirmed by bacterially expressed CmZFH recombinant protein. Deletion analyses indicated that the N-terminal zinc-finger cluster of CmZFH was the key AREB6-binding domain. Through in silico and experimental exploration, we discovered novel transcriptional regulatory components associated with gene expression dynamics under hypoxia that facilitated insect survival.

Large scale discovery of coronavirus-host factor protein interaction motifs reveals SARS-CoV-2 specific mechanisms and vulnerabilities.

Viral proteins make extensive use of short peptide interaction motifs to hijack cellular host factors. However, most current large-scale methods do not identify this important class of protein-protein interactions. Uncovering peptide mediated interactions provides both a molecular understanding of viral interactions with their host and the foundation for developing novel antiviral reagents. Here we describe a viral peptide discovery approach covering 23 coronavirus strains that provides high resolution information on direct virus-host interactions. We identify 269 peptide-based interactions for 18 coronaviruses including a specific interaction between the human G3BP1/2 proteins and an ΦxFG peptide motif in the SARS-CoV-2 nucleocapsid (N) protein. This interaction supports viral replication and through its ΦxFG motif N rewires the G3BP1/2 interactome to disrupt stress granules. A peptide-based inhibitor disrupting the G3BP1/2-N interaction dampened SARS-CoV-2 infection showing that our results can be directly translated into novel specific antiviral reagents.

An experimental test of temperature-dependent selection on mitochondrial haplotypes in Callosobruchus maculatus seed beetles.

Mitochondrial DNA (mtDNA) consists of few but vital maternally inherited genes that interact closely with nuclear genes to produce cellular energy. How important mtDNA polymorphism is for adaptation is still unclear. The assumption in population genetic studies is often that segregating mtDNA variation is selectively neutral. This contrasts with empirical observations of mtDNA haplotypes affecting fitness-related traits and thermal sensitivity, and latitudinal clines in mtDNA haplotype frequencies. Here, we experimentally test whether ambient temperature affects selection on mtDNA variation, and whether such thermal effects are influenced by intergenomic epistasis due to interactions between mitochondrial and nuclear genes, using replicated experimental evolution in Callosobruchus maculatus seed beetle populations seeded with a mixture of different mtDNA haplotypes. We also test for sex-specific consequences of mtDNA evolution on reproductive success, given that mtDNA mutations can have sexually antagonistic fitness effects. Our results demonstrate natural selection on mtDNA haplotypes, with some support for thermal environment influencing mtDNA evolution through mitonuclear epistasis. The changes in male and female reproductive fitness were both aligned with changes in mtDNA haplotype frequencies, suggesting that natural selection on mtDNA is sexually concordant in stressful thermal environments. We discuss the implications of our findings for the evolution of mtDNA.

Identification of novel ejaculate proteins in a seed beetle and division of labour across male accessory reproductive glands.

The male ejaculate contains a multitude of seminal fluid proteins (SFPs), many of which are key reproductive molecules, as well as sperm. However, the identification of SFPs is notoriously difficult and a detailed understanding of this complex phenotype has only been achieved in a few model species. We employed a recently developed proteomic method involving whole-organism stable isotope labelling coupled with proteomic and transcriptomic analyses to characterize ejaculate proteins in the seed beetle Callosobruchus maculatus. We identified 317 proteins that were transferred to females at mating, and a great majority of these showed signals of secretion and were highly male-biased in expression in the abdomen. These male-derived proteins were enriched with proteins involved in general metabolic and catabolic processes but also with proteolytic enzymes and proteins involved in protection against oxidative stress. Thirty-seven proteins showed significant homology with SFPs previously identified in other insects. However, no less than 92 C. maculatus ejaculate proteins were entirely novel, receiving no significant blast hits and lacking homologs in extant data bases, consistent with a rapid and divergent evolution of SFPs. We used 3D micro-tomography in conjunction with proteomic methods to identify 5 distinct pairs of male accessory reproductive glands and to show that certain ejaculate proteins were only recovered in certain male glands. Finally, we provide a tentative list of 231 candidate female-derived reproductive proteins, some of which are likely important in ejaculate processing and/or sperm storage.